CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

turbulentTemperatureRadCoupledMixedFvPatch

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

Like Tree2Likes
  • 1 Post By JoeFriend
  • 1 Post By chriss85

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 12, 2017, 14:41
Question turbulentTemperatureRadCoupledMixedFvPatch
  #1
New Member
 
Join Date: Oct 2016
Posts: 20
Rep Power: 10
JoeFriend is on a distinguished road
Hi, I've been reading the script code of the file turbulentTemperatureRadCoupledMixedFvPatch.C and I can't understand how the radiation term is transmitted from a solid to a fluid region...I believe it's in this part of the code:

scalarField Qr(Tp.size(), 0.0);
if (QrName_ != "none")
{
Qr = patch().lookupPatchField<volScalarField, scalar>(QrName_);
}

scalarField QrNbr(Tp.size(), 0.0);
if (QrNbrName_ != "none")
{
QrNbr = nbrPatch.lookupPatchField<volScalarField, scalar>(QrNbrName_);
mpp.distribute(QrNbr);
}

scalarField alpha(KDeltaNbr - (Qr + QrNbr)/Tp);

I also can't understand the meaning of "alpha".

I would appreciate any help.
Zhiheng Wang likes this.
JoeFriend is offline   Reply With Quote

Old   January 17, 2017, 03:43
Default
  #2
Senior Member
 
Join Date: Oct 2013
Posts: 397
Rep Power: 19
chriss85 will become famous soon enough
I suggest that you look at the analytical solution of the 1D heat conduction equation for three nodes (two boundaries with fixed value, one center node, two different, constant heat conductivities). You can easily solve it analytically and will find the value for the temperature at the patch using the temperatures of the cells.
You can then expand it with the condition that all heat fluxes (thermal and radiative) must be equal on both sides of the patch.
You also need to keep in mind that the patch values in mixedValueFvPatchFields are determined like this:
Value=valueFraction()*refValue()+(1-valueFraction())(Value_cell+refGrad()/deltaCoeffs())
deltaCoeffs equals 1/dx in the one dimensional case, i.e. it's more or less the inverse patch-cell centre distance.
I recently checked the code of this boundary condition and found that it matches with this derivation.
atulkjoy likes this.
chriss85 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 06:50.