CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > 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 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: 7
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, 10:01
Default
  #2
New Member
 
Join Date: Sep 2009
Posts: 9
Rep Power: 7
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

Reply

Thread Tools
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 On
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 10: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 21:11.