2D Unconfined Aquifer Modelling

Hi all,

I'm trying to model a shallow unconfined aquifer with FEFlow, in Planar 2D. The mean-sea level is at 0, and there's approximately 30 meters of unsaturated aquifer above.

Using the Standard (saturated) flow equation doesn't make sense to me, given the unconfined (and unsaturated) portion of the aquifer.

Do I need to use Richards Equations to accurately model this situation, or are there settings I can use to better represent this with standard flow equation?

Cheers,

For me it seems solving the Richards equation is the way to go.

Thanks Bjorn - Could you elaborate why, from a professionals point-of-view? I need to convince my colleagues why it's a better idea, and I'm not certain of the correct theory to apply..

In a 3D model with uncofined conditions, you have the option of using a simplified approach to represent the unsaturated zone (e.g. phreatic). This option is not provided in 2D vertical models. In 2D vertical models you may either solve the saturated Darcy equation or Richards. If you have an unsaturated zone, then Richard seems to me the better option.