CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

chtMultiRegionFoam boundary condition

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

LinkBack Thread Tools Display Modes
Old   March 14, 2014, 05:28
Default chtMultiRegionFoam boundary condition
Join Date: May 2012
Location: Dresden, Germany
Posts: 32
Rep Power: 7
dl6tud is on a distinguished road
I wonder, why the boundary condition for chtMultiRegionFoam is implemented as mixed and not as fixedValue.

How the boundary "turbulentTemperatureCoupledBaffleMixed" works, is explained here

For explanation see picture. Temperatur (T) has to be equal at the boundary as well as the normal heat flow (q). Here, they combine the two boundary conditions into one. Tw is wall Temperature, T1 internal field of phase one T2 of phase two. Delta is the distance cell center to face center, n normal vector.

So basically we have a wall temperature given, fullfilling both boundary conditions. I would base this on fixedValue, but it does not work (does not converge).
In cht they base it on mixed with
Tw = refValue*valueFraction + (1-valueFraction)*(Tinternal + refGradient*delta).
If you use the values provided in turbulentTemperatureCoupledBaffleMixed, you get exactly the formula, but why so complicated? Any hints?

    // Both sides agree on
    // - temperature : (myKDelta*fld + nbrKDelta*nbrFld)/(myKDelta+nbrKDelta)
    // - gradient    : (temperature-fld)*delta
    // We've got a degree of freedom in how to implement this in a mixed bc.
    // (what gradient, what fixedValue and mixing coefficient)
    // Two reasonable choices:
    // 1. specify above temperature on one side (preferentially the high side)
    //    and above gradient on the other. So this will switch between pure
    //    fixedvalue and pure fixedgradient
    // 2. specify gradient and temperature such that the equations are the
    //    same on both sides. This leads to the choice of
    //    - refGradient = zero gradient
    //    - refValue = neighbour value
    //    - mixFraction = nbrKDelta / (nbrKDelta + myKDelta())
Attached Images
File Type: jpg bild.jpg (17.7 KB, 22 views)
dl6tud 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
Radiation interface hinca CFX 15 January 26, 2014 18:11
conjugate boundary condition Daniel_Khazaei OpenFOAM Programming & Development 0 December 31, 2013 14:11
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
Opening Boundary Condition andreachan Main CFD Forum 11 March 19, 2013 17:46

All times are GMT -4. The time now is 15:22.