Abstract

Nonlinear systems of convection-dominated diffusion equations are used as the mathematical model of contamination transport problem which is an important topic in environ mental protection science. An elliptic equation defines the pressure, a convection-diffusion equation expresses the concentration of contamination, and an ordinary differential equation interprets the surface absorption concentration. The transport pressure appears in the equation of the concentration which determines the Darcy velocity and also controls the physical process. The method of conservative mixed volume element is used to solve the flow equation which improves the computational accuracy of Darcy velocity by one order. We use the mixed volume element with the characteristic to approximate the concentration. This method of characteristic not only preserves the strong computational stability at sharp front, but also eliminates numerical dispersion and nonphysical oscillation. In the present scheme, we could adopt a large step without losing accuracy. The diffusion is approximated by the mixed volume element. The concentration and its adjoint vector function are obtained simultaneously, and the locally conservative law is preserved. An optimal second order estimates in l2-norm is derived.

Highlights

  • Numerical simulation of the contamination transport problem is an efficient way to find how to protect the environment and treat the pollution

  • Nonlinear systems of convection-dominated diffusion equations are used as the mathematical model of contamination transport problem which is an important topic in environ mental protection science

  • An elliptic equation defines the pressure, a convection-diffusion equation expresses the concentration of contamination, and an ordinary differential equation interprets the surface absorption concentration

Read more

Summary

Introduction

Numerical simulation of the contamination transport problem is an efficient way to find how to protect the environment and treat the pollution. The problem (1)-(7) is assumed to be Ω-periodic (Dawson, Van Duiji & Wheeler, 1994; Dawson, 1998; Ewing, 1983; Hornung, 1988; Vogt, 1982), that is, all the functions are Ω-periodic This assumption seems physically reasonable, since no-flow boundaries are generally treated by reflection, and because interior flow factors are much more important than boundary effects.

The Preparations
The Procedures
The Conservation of Mass
Convergence Analysis
Conclusions and Discussions
Full Text
Published version (Free)

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