Welcome to our discussion forums. The purpose of these forums is to provide help and technical
support to Hydrus2D/Meshgen2D users. This service is provided for free and therefore we can't
guarantee that all questions will be answered and/or how long time it will take. We do offer
additional customer support and other services - for more information visit the
Hydrus Services page.
If you want to contribute to these discussion forums you must be registered. Registration is free
and very simple - just click Register and fill the registration form.
We recommend you to provide your full name, affiliations, address, and, especially in the
HYDRUS (2D/3D) forum, the name of the company/institution that purchased the Hydrus/Meshgen program.
Questions from fully registered members will be answered with higher priority than other questions.
Questions from anonymous users will have lower priority.
There is an option that allows you to get automated e-mail notifications
about new posting to a specific topic or about all new contributions.
Please login and then click on the "Subscribe to this forum" or "Subscribe to this topic" icon.
(These icons are not visible if you are not logged on)
I wonder if water storage is considered in hydrus2D in a completely saturated domain (horizontal and vertical plane). The vG Model yields C=0 below the phreatic surface however one would expect some small value depending on the materials compressibility.
If there is no consideration, is there a possibility to assign some value to the storage (or capacity in the terminology of Richard's equation) via the soil parameters?
Neither the old HYDRUS-2D, nor the new HYDRUS (2D/3D) consider the storage coefficient in the saturated zone. However, I do have it implemented in the new code and only have it disabled. If you have a license for the new code, I could send you an executable module with this option. Sorry, not for the old one. While we still support that, we are not developing it further for some last 5 years.
thank you for your hints. I decided to perform the modifications in the code (I guess the old one) since I am interested in the incorporation of a pressure dependent storage coefficient (or capacity term) to represent air entrapment below the phreatic surface. Even a very small volume of entrapped gas may influence water pressure propagation in the "saturated"-zone.
My Hydrus code now reads a table of hydraulic properties not only for suction but also for pressure (below the phreatic surface). I assigned C a constant value for pressure and checked that against a quasi 1-dimensional analytical solution (was a bit surprised that Excel has a complementary error function already built in ;-). The results showed quite a strong dependence on temporal discretization but yielded for small time steps almost perfect agreement.
Apparently a non-zero C-term at saturation yields much "smoother" runs (saturated/unsaturated) than the conventional approach with C(h=0)=0 - at least for certain cases and soil properties.
Storage coefficient [1/L] in saturated zone
Posted - 07/27/2010 : 16:40:43
