CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Coupled Pressure-Flow Outlet Boundary Condition (Windkessel)

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

LinkBack Thread Tools Display Modes
Old   January 23, 2012, 10:34
Default Coupled Pressure-Flow Outlet Boundary Condition (Windkessel)
New Member
Join Date: Jan 2012
Posts: 1
Rep Power: 0
torben_w is on a distinguished road
Hi there,

is it possible to use a Windkessel type outlet boundary condition with OpenFoam?

Here, P(t) depends on P(t-1) and Phi(t). How can this be done as a UDF?

The question is similar to an unsolved thread

torben_w is offline   Reply With Quote

Old   February 29, 2012, 06:51
New Member
Michael Bruckner
Join Date: Apr 2009
Location: France
Posts: 27
Rep Power: 10
michaelb is on a distinguished road

have you find a way of implementing the Windkessel boundary conditions in OpenFoam ?
I'd be very interested in knowing how you proceeded.


michaelb is offline   Reply With Quote

Old   March 2, 2012, 12:12
Join Date: Jan 2012
Posts: 32
Rep Power: 8
Vitaly is on a distinguished road
Hi, the previous thread read:


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) ?

I am also starting to work on this issue. I think that you can find P(t-1) using a memory macro, and I think you can find Q(t) using the code supplied. If you only have one outlet for an incompressible fluid, then Q(t) should not be a problem.
In my case, I have multiple outlets. But, I think that if each outlet has the same R and C values, then the flow split should stay the same.
Vitaly is offline   Reply With Quote


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
total pressure boundary problem ==> flow from outlet to inlet!! mrshb4 OpenFOAM 0 November 20, 2010 13:41
Pressure outlet in two-phase flow in horizontal 2D channel AlmostSurelyRob Main CFD Forum 0 November 17, 2010 08:32
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 05:31.