Pressure statistics from the path integral for Darcy flow through random porous media

Abstract

The path integral for classical statistical dynamics is used to determine the properties of one-dimensional Darcy flow through a porous medium with a correlated stochastic permeability for several spatial correlation lengths. Pressure statistics are obtained from the numerical evaluation of the path integral by using the Markov chain Monte Carlo method. Comparisons between these pressure distributions and those calculated from the classic finite-volume method for the corresponding stochastic differential equation show excellent agreement for Dirichlet and Neumann boundary conditions. The evaluation of the variance of the pressure based on a continuum description of the medium provides an estimate of the effects of discretization. Log-normal and Gaussian fits to the pressure distributions as a function of position within the porous medium are discussed in relation to the spatial extent of the correlations of the permeability fluctuations.

References in corpus (7)

• Optimal non-reversible linear drift for the convergence to equilibrium of a diffusion
• Variance Reduction using Nonreversible Langevin Samplers
• Using Perturbed Underdamped Langevin Dynamics to Efficiently Sample from Probability Distributions
• The Monte Carlo method in quantum field theory
• Renormalization Group theory outperforms other approaches in statistical comparison between upscaling techniques for porous media
• User's guide to Monte Carlo methods for evaluating path integrals
• Evaluation of the path integral for flow through random porous media