Abstract
A boundary integral equation method for steady unsaturated flow in non-homogeneous porous media is presented. Steady unsaturated flow in porous media is described by the steady form of the so-called Richards equation, a highly nonlinear Fokker-Planck equation. By applying a Kirchhoff transformation and employing an exponential model for the relation between capillary pressure and hydraulic conductivity, the flow equation is rendered linear in each subdomain of a piece-wise homogeneous material. Unfortunately, the transformation results in nonlinear conditions along material interfaces, giving rise to a jump in the potential along these boundaries. An algorithm developed to solve the nonhomogeneous flow problem is described and verified by comparison to analytical and numerical solutions. The code is applied to examine the moisture distribution in a layered porous medium due to infiltration from a strip source, a model for infiltration from shallow ponds and washes in arid regions.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call