CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent UDF and Scheme Programming

Outlet Pressure-flow coupled boundary condition (Windkessel)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 16, 2010, 18:45
Default Outlet Pressure-flow coupled boundary condition (Windkessel)
  #1
New Member
 
CFDFSI
Join Date: Dec 2009
Posts: 4
Rep Power: 16
cfdfsi is on a distinguished road
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

Code:
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) ?
cfdfsi is offline   Reply With Quote

Old   October 12, 2011, 11:01
Default
  #2
New Member
 
Join Date: Sep 2009
Posts: 9
Rep Power: 16
nacaairfoil is on a distinguished road
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
nacaairfoil is offline   Reply With Quote

Old   December 11, 2015, 23:48
Default
  #3
Member
 
PranjalNewton's Avatar
 
Pranjal Singh
Join Date: Sep 2015
Posts: 34
Rep Power: 10
PranjalNewton is on a distinguished road
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?
PranjalNewton is offline   Reply With Quote

Old   July 21, 2017, 13:30
Unhappy
  #4
Member
 
Ali.E
Join Date: Sep 2010
Location: Lisboa
Posts: 83
Rep Power: 15
behest is on a distinguished road
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
behest is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
pressure outlet boundary condition Sastry FLUENT 4 February 19, 2011 02:33
Turbulent intensity for pressure Outlet Boundary condition Mohsin FLUENT 1 April 30, 2010 11:36
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Pressure Boundary Condition Matt Umbel Main CFD Forum 0 January 11, 2002 11:06
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 14:19


All times are GMT -4. The time now is 19:58.