[cfdfsi]

Hi,

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

density = 100;
sumFlow = 0;
face_t f;
   begin_f_loop(f, t)
      {
       sumFlow = sumFlow + F_FLUX(f,t);
       }
   end_f_loop(f, t)
sumFlow = sumFlow/density;

Now, what is the best way to find P(t-1) on all the faces of the outlet ?
Is there a macro like F_P_M1(f,t) ?

[nacaairfoil]

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

[PranjalNewton]

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?

[behest]

Dear all,
did you manage the problem? I have aslo the same issue, I need previous time step preesure on the wall, please let me know if you found any solutions