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

Mixed boundary condition

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

LinkBack Thread Tools Display Modes
Old   December 19, 2010, 06:05
Default Mixed boundary condition
New Member
Join Date: Aug 2010
Posts: 14
Rep Power: 7
NewFoamer is on a distinguished road

For a scalar field T, I want to implement a boundary condition which ensures following:

A*T + B*(dT/dn) = C

I'm having difficulties understanding how to translate this equation to expressions for valueInternalCoeffs, valueBoundaryCoeffs, gradientInternalCoeffs and gradientBoundaryCoeffs.

Any help would be appreciated.
NewFoamer is offline   Reply With Quote

Old   January 2, 2011, 14:35
New Member
Join Date: Aug 2010
Posts: 14
Rep Power: 7
NewFoamer is on a distinguished road
I can easily understand the fixedValue where:

VIC = 0;
VBC = value;
GIC = -1/delta;
GBC = value/delta;

with delta being the distance from centrum to face and value being the boundary value. The fixedGradient I can also understand:

VIC = 1;
VBC = gradient * delta;
GIC = 0;
GBC = gradient;

but I can not see how to derive expressions for an equation where both the gradient and value occur, such as: A*T + B*(dT/dn) = C

Can anyone help?
NewFoamer is offline   Reply With Quote

Old   January 11, 2011, 10:36
Senior Member
Ben K
Join Date: Feb 2010
Location: Ottawa, Canada
Posts: 140
Rep Power: 10
benk is on a distinguished road
Are you trying to understand the implementation or are you just trying to implement it?

If you just need to implement it, I think the easiest way is to install groovyBC and then:

type groovyBC;
variables "Tinf=298;";
valueExpression "Tinf";
fractionExpression "1.0/(1.0 + K/(mag(delta())*h))";

or you can take a look at the wallHeatTransfer boundary condition which does the same thing.

This is for a convective heat transfer boundary condition which looks like:

K*dT/dx = h(T-Tinf) or
h/K*T - dT/dx = h/K*Tinf
I really don't understand this implementation though and I'm hoping somebody can take the time to explain this to me too!
benk is offline   Reply With Quote

Old   February 17, 2011, 11:35
New Member
Raimonds Vilums
Join Date: Oct 2010
Posts: 17
Rep Power: 7
argonaut is on a distinguished road
I tried to explain this implementation here.

It took me several days of thinking and browsing the forum to understand this implementation. Hope this will help for others to find answers quicker.
argonaut 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
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 05:05
asking for Boundary condition in FLUENT Destry FLUENT 0 July 27, 2010 00:55
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 02:53
problem with boundary condition??? smn CFX 5 November 24, 2009 07:37
Fluent Treatment of mixed boundary condition HELP Amr FLUENT 0 May 26, 2006 05:46

All times are GMT -4. The time now is 11:38.