I will create a new CFD model based on Aneurysm example and have 2 questions about the constraint of outvar ‘integral_continuity’.

*Environment
Simnet version : 21.06

*Questions

Is the value of ‘integral_continuity’ calculated in IntegralContinuity.make_node() in navier_stokes.py?
If it was wrong, please let me know where it is calculated.

Why is the constraint value of ic_2 integral_continuity negative(-2.540) as line 119 in aneurysm.py?
I think the value should be positive(2.540) so that it is equivalent as ic_1.

That is correct. The integral continuity value is calculated in IntegralContinuity.make_node().

This is because for the second integral continuity plane, the plane normal is in the opposite direction of the flow, but the plane normal for the first integral continuity plane is in the flow direction.

That is correct. The integral continuity value is calculated in IntegralContinuity.make_node() .

understood.

This is because for the second integral continuity plane, the plane normal is in the opposite direction of the flow, but the plane normal for the first integral continuity plane is in the flow direction.

I have a new question.
Where is the normal direction of the second plane written?
Could you please let me know how to understand that the normal is in the opposite direction of the flow?

@mnabian Hi, If we plan to reconstruct the field and we didn’t know the value of flow rate, can we define a loss term based on equality of volume flow rate of outlet and inlet based on integral constraint? How?

What’s interesting is that the code from GitHub did seem to replicate the figure in the documentation accurately, but I wonder if this is the result of a mistake (or maybe it doesn’t affect the output too much, but it could be a typo nonetheless).