Modeling Dewatered Domains in Multilayer Analytic Element Models

This proposed technique allows sensible and numerically stable behavior in multilayer analytic element models when layers dewater. When saturated thickness approaches zero in an unconfined or fresh/salt interface domain, the domain transitions to a very thin confined domain with a minimum saturated thickness M. M is an adjustable input parameter, so you can make the horizontal flow in dewatered domains negligibly small by making the minimum saturated thickness very small. Vertical flows can pass through a dewatered domain, whether it is near the surface or at depth. For example, recharge may pass through a shallow dewatered layer to a deeper layer that is not dewatered. This approach is examined in detail in an example multilayer model of mine dewatering.