[3nhu1]

Hello,

Please could I ask for some assistance, and thank you especially for all the help previously in the last few months

I have a previous boundary condition implemented that had a 2 element windkessel (see picture).
It required me to create an additional variable for pressure at the previous timestep i.e. P(t-1)
What we had done for this was create an additional variable called pprev. Set it to scalar. In the default domain Modified field, pprev was activated as an Algebraic Equation and set to Pressure.
This additional variable in the expression could then be used in the expression for the pout boundary condition using the final equation for P(t) in the picture by calling for P(t-1) as areaAve(pprev)@out1 (the model has multiple outlets that does similar)

I now need to implement a 3 element windkessel (see 2nd picture). The methodology is fairly similar, but I'm having some difficulty as it requires creating an additional variable for flow rate Q at the previous timestep i.e. Q(t-1).
Mainly as there is no global variable for flow. (Although density is constant in the model).
A colleague suggested using velocity. First I have checked and tried the 2 inbuilt functions, and by using a normal velocity at outlet give a reasonable answer.

Additional variable
NormalVelocity = VelDotn
and
CEL
VelDotn = (Velocity u*Normal X + Velocity v*Normal Y + Velocity w*Normal Z)

Then option 1 for q -
areaInt(NormalVelocity)@OutAorta
and option 2 for q -
Qout_option5=areaAve(NormalVelocity)@OutAorta * area()@OutAorta

Option and Option 2 are equal and very similar to the massflow()@OutAorta / Density when post processed.

Knowing this - What would be a reasonable way to implement Q(t-1)?


Ref 2017 Acedo- Optimisation of a windkessel system for arterial simulation (thesis)