CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Mixed boundary condition

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

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

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
Default
  #2
New Member
 
Join Date: Aug 2010
Posts: 14
Rep Power: 6
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
Default
  #3
Senior Member
 
Ben K
Join Date: Feb 2010
Location: Ottawa, Canada
Posts: 140
Rep Power: 9
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:

Quote:
rightWall
{
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:

Code:
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
Default
  #4
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

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
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 06:05
asking for Boundary condition in FLUENT Destry FLUENT 0 July 27, 2010 01:55
External Radiation Boundary Condition for Grid Interface CFD XUE FLUENT 0 July 9, 2010 03: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 06:46


All times are GMT -4. The time now is 16:56.