Dear all,
I have run Hydrus 1D for infiltration of water into the soil profile. I use pressure head boundary condition at the Upper Boundary Condition. The input time Variable Boundary Conditions is shown in the picture. After I run the calculations I see the wetting boundary in the Profile information of the Water Content are all becomes equal to the saturated water content. but if we see on the input data when the times equal 1 day to 4 days, 5 days the pressure head is equal to Zero opposite to the Water content that all the print times the wetting front becomes saturated all the print times? So what wrong with this problems? I have try so many way but it remain the same!
In my opinion, this water content should be become unsaturated when the input pressure head become to zero. Could you help me to solve this problems?
Thank in advance
Best Regards
K. Makara
Water Content Postprocessing
Water Content Postprocessing
 Attachments

 Capture.PNG (37.7 KiB) Viewed 286 times

 Capture1.PNG (15.96 KiB) Viewed 286 times
Re: Water Content Postprocessing
Quote “So what wrong with this problems?”
Problem seems to be your basic understanding of soil physics. A zero pressure head BC does not mean a zero flux BC, it means a zero pressure head at the boundary. Note that the zero pressure head corresponds to full saturation (see the retention curve).
If you want to alternate between pressure head and flux BCs, you need to use the “Variably Pressure Head/Flux” BC and tell the code when the head or flux (which can be specified to be zero) is to be considered. A Dirichlet (head) BC never (or rarely) leads to zero flux.
J.
Problem seems to be your basic understanding of soil physics. A zero pressure head BC does not mean a zero flux BC, it means a zero pressure head at the boundary. Note that the zero pressure head corresponds to full saturation (see the retention curve).
If you want to alternate between pressure head and flux BCs, you need to use the “Variably Pressure Head/Flux” BC and tell the code when the head or flux (which can be specified to be zero) is to be considered. A Dirichlet (head) BC never (or rarely) leads to zero flux.
J.