Please wait...
×

Error

Posted Wed, 30 Apr 2008 10:45:35 GMT by Juliette
Hello,
I am working on a 20 slices model under unsaturated flow conditions (and transient regim). A problem appears around the wells : the head drasticaly decreases (actually it goes out of the physical domain). I tried a horizontal mesh refinement around the well but it didn't work. Note that the same model defined under saturated flow conditions well-runs !!

Any ideas about where the problem comes from?

Thanks in advances,
Juliette
Posted Fri, 23 May 2008 09:15:04 GMT by Alexander Renz
if the head falls below the model domain, the cells around the well will surely be dry (negative pressure). Dry cells become less permeable, together with a constant pumping rate of the well this will further decrease the pressure/head.

If you trust the numerical result, you might want to switch of the wells if they fall dry; you can use boundary complemtary constraints (BCC) for this.
Posted Fri, 23 May 2008 11:37:49 GMT by Christopherus Braun
Maybe your numerical problems came from the parametric model you use for unsaturated flow. if you define some very sharp curves FEFFLOW may not produce meaningful results.

You have to keep in mind, that for the stability of the numerical simulation not the total value of capillar pressure or relative  conductivity are important, but there deviation over saturation. So avoid any steep functions.

I made some good experience with the modified Van Genuchten approach. Here you have two possibilities to smooth your simulation:
1. Reduce the parameter A: this doews introduce some capillara difussivity into the simulation and helps a lot.
2. Use a value for the parameter delta that does not make the relative permeability curye to steep (e.g. 1.1 - 1.3)

Of course the parameters you use should have some physical meanings, but there maybe a large range in data.


The other thing you could try is a mesh refinement. Not so much in hirizontal direction, but in vertical direction. This shoud help much more the unsaturated processes take place in the vertical direction.

The vertical element length you should use has also something to do with the parametric model you use. You have to assure that an unsaturated profile could be computed an not just " wet or dry ". There is a rule of the thumb: if you use a (modified) Van Genuchten model the vertical spacing of your slices should be in the range of 1/A

Zebra
Posted Mon, 26 May 2008 15:55:48 GMT by Juliette
thanks for the advices!

You must be signed in to post in this forum.