Outlet Pressure-flow coupled boundary condition (Windkessel)
I'm trying to apply a windkessel type boundary condition at the outlet, which is of the form:
P(t) + k*P'(t) = Q(t), where k is a known constant
To do this, I model the outlet as a 'pressure-outlet' and use a UDF to impose P(t) on the outlet faces.
To calculate P(t), I approximate P'(t) by [P(t) - P(t-1)]/delta_t, and arrive at the closed form expression
P(t) = Function[P(t-1), Q(t), k, delta_t]
My question is:
how do I specify this equation in the UDF ?
At present, I'm using the DEFINE_PROFILE(pressure, t, i) macro,
calculating Q(t) by
Is there a macro like F_P_M1(f,t) ?
I am also facing the same problem. Does anyone has answer for this. Also since Q is evaluated at current time step. I need to iterate to get p. How should I do the iteration using UDF?
Thanks a lot
I'm also stuck on this. I made UDM for storing P(t-1) and initialised it. So, it should work but in the results, I see that the compliance phenomenon isn't observed. I've modelled wall material as elastic. Does it need to be hyperelastic or something to get FSI coupling behavior?
|All times are GMT -4. The time now is 20:40.|