Unformatted text preview: NUMERICAL ACCURACY AND ROBUSTNESS
Continuation method for highly heterogeneous media Some special considerations were given to the non-linear aspects of the SWIM problem.
Indeed, Seawater Intrusion Modeling is “doubly non-linear” because of the non-linearity of
the Dupuit Boussinesq equation (head-dependent transmissivity) and because of the
nonlinearity due to the interception of the salt interface by the aquifer’s substratum. This
non-linearity is even stronger when the domain is highly heterogeneous (K(x,y). For these
cases, we used a special iterative “re-start method”, also known as “continuation” or
“homotopy” method. The continuation parameter is chosen to be σlnK, the standard
deviation of lnK. Modélisation stochastique de l'intrusion saline en 2D plan 73 The method is implemented iteratively as follows. Step 1: we generate an aquifer with
mildly heterogeneous log-permeability (lnK(x,y)) characterized by its mean, standard
deviation, and covariance function, u...
View Full Document
This document was uploaded on 01/19/2014.
- Winter '14
- The Land, In French