Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

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.

Switch the boundary condition from flux to pressure head (Richards’Equation)

Please login with a confirmed email address before reporting spam

I read some paper, in the old version of COMSOL, we can choose a Mixed boundary condition (mix of Dirichlet, Newman and Cauchy conditions). in order to automatically transfer the boundary condition form pressure head to Flux, or Flux to pressure head.

But in the new version of COMSOL 4.2 or 4.3, I can not found that kind of setting,we would like to know how to switch the boundary condition from flux to pressure head during the calculation. We started from a flux boundary for rainfall infiltration (with negative pressure head), and when the pressure head at the surface boundary becomes positive, we would like to swith it to the pressure head boundary. We have no idea how to implement it in COMSOL4.3(Richards’Equation).

Could any expert have some suggestion? I appreciate for your help.

3 Replies Last Post Jun 18, 2014, 6:13 a.m. EDT
COMSOL Moderator

Hello Andy Huang

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.


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Jun 6, 2014, 10:39 a.m. EDT
Hi,
I know this discussion is quite old but I was wondering if you found a solution (or if someone else has any suggestion) as I'm having exactly the same problem. I am using the same module (Richards' equation) and the same version of COMSOL (4.3); I am analysing water flow through an embankment during a flood event and I would like to switch boundary condition from "no flow" to "pressure head" when the pressure head at the toe of the embankment becomes positive.

I would also like to automatically change the boundary condition from "hydraulic head" to "no flow" according to the elevation of the points on a boundary. I'll try to explain this point more clearly: to represent the rising water level on the river side of the embankment I am applying a boundary condition in terms of hydraulic head varying over time (i.e. the hydrograph of my flood event); anyway on the portion of the embankment boundary which is not wet by the river, i.e. when the elevation y is more than the hydraulic head that I am applying, I need a "no flow" boundary condition. I have managed to change the value of the hydraulic head according to the elevation with a simple "if" statement (if(y<H(t),H(t),0)), but in the "else" scenario I need to change the type of boundary condition rather than its value.

Is there any way to do this? I was thinkink about a "conditional enabling", where I could enable the boundary condition "no flow" on the boundary only if y>H(t) and override the hydraulic head boundary condition only in this case, but I've found nothing about it. Could anyone have suggestion about this?
I would really appreciate any help. Thank you.
Hi, I know this discussion is quite old but I was wondering if you found a solution (or if someone else has any suggestion) as I'm having exactly the same problem. I am using the same module (Richards' equation) and the same version of COMSOL (4.3); I am analysing water flow through an embankment during a flood event and I would like to switch boundary condition from "no flow" to "pressure head" when the pressure head at the toe of the embankment becomes positive. I would also like to automatically change the boundary condition from "hydraulic head" to "no flow" according to the elevation of the points on a boundary. I'll try to explain this point more clearly: to represent the rising water level on the river side of the embankment I am applying a boundary condition in terms of hydraulic head varying over time (i.e. the hydrograph of my flood event); anyway on the portion of the embankment boundary which is not wet by the river, i.e. when the elevation y is more than the hydraulic head that I am applying, I need a "no flow" boundary condition. I have managed to change the value of the hydraulic head according to the elevation with a simple "if" statement (if(yH(t) and override the hydraulic head boundary condition only in this case, but I've found nothing about it. Could anyone have suggestion about this? I would really appreciate any help. Thank you.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Jun 10, 2014, 5:15 a.m. EDT

Hi,
I know this discussion is quite old but I was wondering if you found a solution (or if someone else has any suggestion) as I'm having exactly the same problem. I am using the same module (Richards' equation) and the same version of COMSOL (4.3); I am analysing water flow through an embankment during a flood event and I would like to switch boundary condition from "no flow" to "pressure head" when the pressure head at the toe of the embankment becomes positive.

I would also like to automatically change the boundary condition from "hydraulic head" to "no flow" according to the elevation of the points on a boundary. I'll try to explain this point more clearly: to represent the rising water level on the river side of the embankment I am applying a boundary condition in terms of hydraulic head varying over time (i.e. the hydrograph of my flood event); anyway on the portion of the embankment boundary which is not wet by the river, i.e. when the elevation y is more than the hydraulic head that I am applying, I need a "no flow" boundary condition. I have managed to change the value of the hydraulic head according to the elevation with a simple "if" statement (if(y<H(t),H(t),0)), but in the "else" scenario I need to change the type of boundary condition rather than its value.

Is there any way to do this? I was thinkink about a "conditional enabling", where I could enable the boundary condition "no flow" on the boundary only if y>H(t) and override the hydraulic head boundary condition only in this case, but I've found nothing about it. Could anyone have suggestion about this?
I would really appreciate any help. Thank you.


It is better to use mixed boundary condition or you can define step function for the embankment.
[QUOTE] Hi, I know this discussion is quite old but I was wondering if you found a solution (or if someone else has any suggestion) as I'm having exactly the same problem. I am using the same module (Richards' equation) and the same version of COMSOL (4.3); I am analysing water flow through an embankment during a flood event and I would like to switch boundary condition from "no flow" to "pressure head" when the pressure head at the toe of the embankment becomes positive. I would also like to automatically change the boundary condition from "hydraulic head" to "no flow" according to the elevation of the points on a boundary. I'll try to explain this point more clearly: to represent the rising water level on the river side of the embankment I am applying a boundary condition in terms of hydraulic head varying over time (i.e. the hydrograph of my flood event); anyway on the portion of the embankment boundary which is not wet by the river, i.e. when the elevation y is more than the hydraulic head that I am applying, I need a "no flow" boundary condition. I have managed to change the value of the hydraulic head according to the elevation with a simple "if" statement (if(yH(t) and override the hydraulic head boundary condition only in this case, but I've found nothing about it. Could anyone have suggestion about this? I would really appreciate any help. Thank you. [/QUOTE] It is better to use mixed boundary condition or you can define step function for the embankment.

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Jun 18, 2014, 6:13 a.m. EDT
Thanks for your reply, but what do you mean with mixed boundary condition? How can I apply this in version 4.3?
Thanks for your reply, but what do you mean with mixed boundary condition? How can I apply this in version 4.3?

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.