mean¶

poreana.density.mean(density, is_print=True, int_limit=2.5)¶

This function uses the desnity slope obtained from function poreana.density.bins() to calculate a weighted density inside the pore

\[\langle\rho\rangle =\frac{\int\rho(r)dA(r)}{\int dA(r)}.\]

In a discrete form, following formula is evaluated

\[\langle\rho\rangle=\frac{\sum_{i=1}^n\rho(r_i)A(r_i)}{\sum_{i=1}^nA(r_i)}\]

with the partial area

\[A(r_i)=\pi(r_i^2-r_{i-1}^2)\]

of radial bin \(i\).

Parameters:

Returns:

meandictionary: Weighted mean density in \(\frac{\text{#}}{\text{nm}^3}\) num_dens and \(\frac{\text{kg}}{\text{m}^3}\) dens