Posted Tue, 10 Sep 2019 18:44:41 GMT by Francesconi, Matteo

I am here to ask a question about the free surface setting called "constrained head".

I have a model with that setting applied to the top of the model domain (including the prevent inflow option). As far as I understood, that setting should enable the activation of hydraulic head BC with a max flow rate set to 0 at any node where groundwater level reaches the topographic surface, in order to remove exceeding water.

However, some nodes show a hydraulic head value higher than the topographic elevation. Why does it happen? Shouldn't this setting allow the hydraulic head value be at most equal to the topopgrahic elevation?

Thanks in advance,
Posted Wed, 11 Sep 2019 06:40:46 GMT by Peter Schätzl

in general you are right, the setting enables a head bc with max Q = 0. This BC is applied iteratively, as setting the bc at one node may influence the neighboring ones. To avoid endless iteration, though, there's a max. number of loops - if I remember correctly this is 30. Whatever nodes remain not matching the conditions after this number of iterations can show a higher head.

Hope this helps to understand what's going on.

Posted Wed, 11 Sep 2019 07:37:43 GMT by Francesconi, Matteo
Hi Peter,

Thank you for the quick reply.

Yes, that helps me understanding more what FEFLOW does.

According to this, the hydraulic head values higher than elevation that come out from calculation should be taken just as a qualitative indicator (in those nodes flooding has occured) and not as a quantitative indicator (the higher the hydraulic head is compared to the elevation, the higher is the amount of flooding water in that node), is that right?


You must be signed in to post in this forum.