HR: 0830h
AN: H21D-0884    [PDF]
TI: An Operator Splitting Algorithm for Numerical Simulation of Transport of Reactive Compounds in Porous Media
AU: * Hernandez -Rendon, M
EM: carmenhr@geofisica.unam.mx
AF: Instituto de Geofisica, UNAM, Circuito de la Investigacion S/N, Cd. Universitaria, Mexico, D.F., 04510 Mexico
AB: Chemical and/or biological processes are of great importance in modeling transport of pollutants in many hydrogeological settings. When considering these phenomena, a system of partial differential equations (PDE's) nonlinear and coupled in the reaction terms is obtained. Besides these numerical difficulties and the non-physical oscillations in advection-dominated flows, stiffness may be introduced into the system since the range of characteristic reaction time is huge. In this work an operator splitting method is applied to simulate transport of species undergoing biodegradation in porous media. It is the $\Theta$ scheme which is second order. The mathematical model takes into account the mass balance of all the chemical species involved in the successive transformations. A Monod type kinetics is assumed for biodegradation. The advection-diffusion operator is solved implicitly in the first and third subintervals and the biochemical terms are solved implicitly in the intermediate one. In this way, stiff matrices are avoided. By a special treatment of the source terms, the system of PDE's is totally decoupled. Numerical results are obtained in one spatial dimension.
DE: 1832 Groundwater transport
DE: 3220 Nonlinear dynamics
DE: 3230 Numerical solutions
SC: Hydrology [H]
MN: 2003 Fall Meeting