Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Unsaturated flow (Richard's Equation) convergence issue with very dry initial condition

Please login with a confirmed email address before reporting spam

Hi,

I am trying to simulate a 2D infiltration experiment on a dry homogeneous soil column. The initial volumetric water content (wc) before applying any infiltration was 0.05. To achieve this initial wc, I had to simulate a time-dependent study (study 1) with initial pressure values assigned as -1500000 Pa. All the details (i.e., model geometry, variables, and boundary conditions) can be found in the attached model file. Now, when I am applying inflow at the top boundary in the second time-dependent study (in which the initial values of variables solved for and not solved for are assigned as the last value obtained from the solution of study 1), the model isn't converging. I suppose, as water starts to come into the soil the pressure becomes >=0 and COMSOL is having issues handling this large pressure gradient. I would very much appreciate if anyone could take a look at my model and suggest me anything to work it out and run.

Thank you for your time and effort. I really appreciate it.

Abdullah



0 Replies Last Post 23 mars 2021, 05:38 UTC−4
COMSOL Moderator

Hello Abdullah Al Mamun

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.