Paper detail

On reconstruction of dynamic permeability and tortuosity from data at distinct frequencies

This article focuses on the mathematical problem of reconstructing dynamic permeability $K(ω)$ of two-phase composites from data at different frequencies, utilizing the analytic structure of the Stieltjes function representation of $K(ω)$ derived by Avellaneda and Tortquato in \cite{avellaneda1991rigorous-link-b}, which is valid for all pore space geometry. The integral representation formula (IRF) for dynamic tortuosity is derived and its analytic structure exploited for reconstructing the function from a finite data set. All information of pore-space microstructure is contained in the measure of the IRF. The theory of multipoint Padé approximates for Stieltjes functions guarantees the existence of relaxation kernels that can approximate the dynamic permeability function and the dynamic tortuosity function with high accuracy. In this paper, a numerical algorithm is proposed for computing the relaxation time and the corresponding strength for each element in the relaxation kernels; in the frequency domain, the approximation can be regarded as approximating the Stieltjes function by rational functions with simple poles and positive residues. With the idea from dehomogenization, we also established the exact relations between the moments of the positive measures in the IRFs of permeability and tortuosity with two important parameters in the theory of poroelasticity: the infinite-frequency tortuosity $α_\infty$ for general case and the weighted volume-to-surface ratio $Λ$ for the JKD model. From these relations, we suggest a new way for evaluating these two microstructure-dependent parameters from a finite data set of permeability and given values of porosity $ϕ$ and pore fluid kinetic viscosity $ν$. Numerical results for JKD permeability and tortuosity are presented.