Page 126 - Fundamentals of Computational Geoscience Numerical Methods and Algorithms

114 5 Simulating Chemical Dissolution Front Instabilities in Fluid-Saturated Porous Rocks

Dimensionless concentration 1 Dimensionless concentration 1

0
0 10 0 0 10
x x
(τ = 0.25) (τ = 0.625)

Dimensionless concentration 1

0
0 10
x
(τ = 0.8)
Fig. 5.3 Comparison of numerical solutions with analytical ones at different time instants (Dimen-
sionless concentration): the thick line shows the numerical results, while the thin line shows the
corresponding analytical solutions

Dimensionless pressure 180 Dimensionless pressure 180

150

100 100
0 10
0 10
x x
(τ = 0.25) (τ = 0.625)
Dimensionless pressure 180

140

100
0 10
x
(τ = 0.8)
Fig. 5.4 Comparison of numerical solutions with analytical ones at different time instants (Dimen-
sionless pressure): the thick line shows the numerical results, while the thin line shows the corre-
sponding analytical solutions