LK81

Hi all,

I am applying some Dirichlet conditions (with max constraint) to the locations of some river nodes in a steady state Richard´s model. The fixed head conditions are equal to the respective node elevation. Although the calculated hydraulic head in the steady state solution for most of these nodes is higher than the corresponding elevation, the BCs are deactivated. Has anyone had a similar experience or knows what can be the reason for this? This causes my model to be saturated and "pressurized" at the fixed BC nodes.

I have already checked that there is a gradient from the neighboring nodes towards the boundary condition, so this could not be the problem. Also, I do not want to use fluid-transfer in this case which could be an alternative.

Many thanks,
L