A Boundary Integral Equation Method for Steady Two-Dimensional Flow in Partially Saturated Media
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
The governing equation for steady flow in a partially saturated, porous medium can be written in a linear form if one adopts a hydraulic conductivity function that is exponential in the capillary-pressure head. The resulting linear field equation is well suited to numerical solution by the boundary integral equation method (BIEM). The exponential conductivity function is compared to a more complex form often assumed for tuffs, and is found to be a reasonable approximation over limited ranges of pressure head. A computer code based on the BIEM is described and tested. The BIEM is found to exhibit quadratic convergence with element size reduction on smooth solutions and on singular problems, if mesh grading is used. Agreement between results from the BIEM code an a finite-element code that solves the fully nonlinear problem is excellent, and is achieved at a substantial advantage in computer processing time.
- Research Organization:
- Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE
- DOE Contract Number:
- AC04-76DP00789
- OSTI ID:
- 138166
- Report Number(s):
- SAND--90-0253; ON: DE91016056
- Country of Publication:
- United States
- Language:
- English
Similar Records
A boundary integral method for steady unsaturated flow in nonhomogeneous media
LLUVIA: A program for one-dimensional, steady-state flow through partially saturated porous media