Published on

Ogata and Banks Analytical Solution

Authors

Author

  • Name
    Albert Nardi
    mail

Background

Numerical modelling is nowadays the reference tool to simulate and predict the behaviour of hydrogeological systems. This fact, however, has not stopped the research and the continuous development of analytical solutions for the different porous media equations. Analytical solutions are still useful tools today for system understanding and fast prediction scenarios.

At Amphos21, we also use analytical solutions to verify code developments. They are valuable resources to validate numerical implementations and explore the numerical method limits for handling certain physics (for example advection predominant problems when using FEM).

In the last years, we have been involved in the development of different reactive transport codes. For their verification, we have been using intensively the analyitical solution Differential Equation of Longitudinal Dispersion in Porous Media provided by Akio Ogata and R. B. Banks (1961). This solution does not involve chemical reactions, we use it mainly to verify that the conservative tracers are handled correctly. This way we have fast feedback to determine if possible differences in reactive species are caused or not by errors or wrong settings in the solute transport equation resolution.

See this use on one of the iCP verification problems (Figure 1). The chloride, which is a conservative species in this case, fits to the analytical solution.

Figure 1. Comparison between iCP, Phreeqc and a numerical solution (for Cl) for the cation exchange (example 11) problem extracted from the Phreeqc Manual.

Ogata and Banks 1961

Ogata and Banks analytical solution dates from 1961 and can be found in the following USGS report.

It solves the ADE equation for a Continuous Source of Infinite Duration and a 1D domain:

R∂c∂t=−u∂c∂x+D∂2c∂xR \frac{\partial c}{\partial t} = - u \frac{ \partial c}{\partial x} + D\frac{ \partial ^2 c}{\partial x}

with the following boundary and initial conditions:

c(0,t)=co;t≥0c(∞,t)=0;t≥0c(x,0)=0;x≥0\begin{aligned} c(0, t) = c_o; t \geq 0 \\ c(\infty, t) = 0; t \geq 0 \\ c(x, 0) = 0; x \geq 0 \\ \end{aligned}

where cc is the concentration [ML−3][ML^{- 3}], xx is the distance [L][L], RR is the retardation factor [−][-], DD the effective dispersion/difusion [L2T−1][L^2 T^{- 1}], uu the flow velocity [LT−1][LT^{-1}] and coc_o the concentration at the upstream boundary [ML−3][ML^{- 3}].

The analytical solution is:

c(x,t)=co2[erfc(x−ut2Dt/R)+exp(uxD/R)erfc(x+ut2Dt/R)]\begin{aligned} c(x, t) = \frac{c_o}{2} \left[ \text{erfc} \left( \frac{x-ut}{2 \sqrt{Dt/R} } \right) +\text{exp} \left( \frac{ux}{D/R } \right) \text{erfc} \left( \frac{x+ut}{2 \sqrt{Dt/R} } \right) \right] \end{aligned}

Experiment

We have implemented the Ogata and Banks solution in a web application to predict quickly the arrival concentrations at a distance xx. We use it to have first estimations of tracer arrival times and the role that play different parameters on modifying the arrival curve shape.

Just set your simulation parameters and play Add Curve to get your results. We hope you enjoyed it!

References

Ogata, A., & Banks, R. B. (1961). A solution of the differential equation of longitudinal dispersion in porous media: fluid movement in earth materials. US Government Printing Office.