def updateCoeffs( self ):
try:
if self.updated():
return
from Foam.finiteVolume import volScalarField
psip = volScalarField.ext_lookupPatchField(self.patch(), word( "psi" ) )
pp = volScalarField.ext_lookupPatchField( self.patch(), word( "p" ) )
self == ( psip * pp )
fixedValueFvPatchScalarField.updateCoeffs( self )
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def updateCoeffs( self ):
try:
if self.updated():
return
from Foam.finiteVolume import volScalarField
psip = volScalarField.ext_lookupPatchField(self.patch(), word( "psi" ) )
pp = volScalarField.ext_lookupPatchField( self.patch(), word( "p" ) )
self.gradient().ext_assign( psip * pp.ext_snGrad() + psip.ext_snGrad() * pp )
fixedGradientFvPatchScalarField.updateCoeffs( self )
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def updateCoeffs( self ) :
try:
if self.updated() :
return
from Foam.finiteVolume import volScalarField
from Foam.OpenFOAM import word
rhop = volScalarField.ext_lookupPatchField( self.patch(), word( "rho" ) )
from Foam.finiteVolume import volVectorField
rhoUp =volVectorField.ext_lookupPatchField( self.patch(), word( "rhoU" ) )
T = volScalarField.ext_lookupObject( self.db(), word( "T" ) )
patchi = self.patch().index()
Tp = T.ext_boundaryField()[patchi]
Tp.evaluate()
from Foam.OpenFOAM import IOdictionary
thermodynamicProperties = IOdictionary.ext_lookupObject( self.db(), word( "thermodynamicProperties" ) )
from Foam.OpenFOAM import dimensionedScalar
Cv = dimensionedScalar( thermodynamicProperties.lookup( word( "Cv" ) ) )
self.valueFraction().ext_assign( rhop.ext_snGrad() / ( rhop.ext_snGrad() - rhop * self.patch().deltaCoeffs() ) )
self.refValue().ext_assign( 0.5 * rhop * ( rhoUp / rhop ).magSqr() )
self.refGrad().ext_assign( rhop * Cv.value() * Tp.ext_snGrad() +\
( self.refValue() - ( 0.5 * rhop.patchInternalField()\
* ( rhoUp.patchInternalField() /rhop.patchInternalField() ).magSqr() ) ) * self.patch().deltaCoeffs() )
mixedFvPatchScalarField.updateCoeffs( self )
pass
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def flux(self):
try:
from Foam.finiteVolume import volScalarField
Kw = volScalarField.ext_lookupPatchField(self.patch(), self.KName_)
Tw = self
return Tw.ext_snGrad() * self.patch().magSf() * Kw
except Exception, exc:
import sys, traceback
traceback.print_exc(file=sys.stdout)
raise exc
def flux( self ):
try:
from Foam.finiteVolume import volScalarField
Kw = volScalarField.ext_lookupPatchField( self.patch(), self.KName_ )
Tw = self
return Tw.ext_snGrad()*self.patch().magSf()*Kw
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def updateCoeffs( self ):
try:
if self.updated():
return
from Foam.OpenFOAM import scalar
neighbourField = self.coupleManager_.neighbourPatchField( scalar )
from Foam.finiteVolume import volScalarField
K = volScalarField.ext_lookupPatchField( self.patch(), self.KName_ )
from chtMultiRegionFlux.r1_5.derivedFvPatchFields import solidWallTemperatureCoupledFvPatchScalarField
self.gradient().ext_assign( solidWallTemperatureCoupledFvPatchScalarField.ext_refCast( neighbourField ).flux() / K )
fixedGradientFvPatchScalarField.updateCoeffs( self )
pass
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def updateCoeffs(self):
try:
if self.updated():
return
from Foam.OpenFOAM import scalar
neighbourField = self.coupleManager_.neighbourPatchField(scalar)
from Foam.finiteVolume import volScalarField
K = volScalarField.ext_lookupPatchField(self.patch(), self.KName_)
from chtMultiRegionFlux.r1_5.derivedFvPatchFields import solidWallTemperatureCoupledFvPatchScalarField
self.gradient().ext_assign(
solidWallTemperatureCoupledFvPatchScalarField.ext_refCast(
neighbourField).flux() / K)
fixedGradientFvPatchScalarField.updateCoeffs(self)
pass
except Exception, exc:
import sys, traceback
traceback.print_exc(file=sys.stdout)
raise exc
def updateCoeffs(self):
try:
if self.updated():
return
from Foam.finiteVolume import volScalarField
from Foam.OpenFOAM import word
rhop = volScalarField.ext_lookupPatchField(self.patch(),
word("rho"))
from Foam.finiteVolume import volVectorField
rhoUp = volVectorField.ext_lookupPatchField(
self.patch(), word("rhoU"))
T = volScalarField.ext_lookupObject(self.db(), word("T"))
patchi = self.patch().index()
Tp = T.ext_boundaryField()[patchi]
Tp.evaluate()
from Foam.OpenFOAM import IOdictionary
thermodynamicProperties = IOdictionary.ext_lookupObject(
self.db(), word("thermodynamicProperties"))
from Foam.OpenFOAM import dimensionedScalar
Cv = dimensionedScalar(thermodynamicProperties.lookup(word("Cv")))
self.valueFraction().ext_assign(
rhop.ext_snGrad() /
(rhop.ext_snGrad() - rhop * self.patch().deltaCoeffs()))
self.refValue().ext_assign(0.5 * rhop * (rhoUp / rhop).magSqr())
self.refGrad().ext_assign( rhop * Cv.value() * Tp.ext_snGrad() +\
( self.refValue() - ( 0.5 * rhop.patchInternalField()\
* ( rhoUp.patchInternalField() /rhop.patchInternalField() ).magSqr() ) ) * self.patch().deltaCoeffs() )
mixedFvPatchScalarField.updateCoeffs(self)
pass
except Exception, exc:
import sys, traceback
traceback.print_exc(file=sys.stdout)
raise exc
def updateCoeffs( self ) :
try:
if self.updated() :
return
from Foam.meshTools import directMappedPatchBase
mpp = directMappedPatchBase.ext_refCast( self.patch().patch() )
nbrMesh = mpp.sampleMesh()
intFld = self.patchInternalField()
if self.interfaceOwner( nbrMesh ):
# Note: other side information could be cached - it only needs
# to be updated the first time round the iteration (i.e. when
# switching regions) but unfortunately we don't have this information.
distMap = mpp.map()
from Foam.finiteVolume import fvMesh
nbrPatch = fvMesh.ext_refCast( nbrMesh ).boundary()[ mpp.samplePolyPatch().index() ]
# Calculate the temperature by harmonic averaging
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
from Foam.finiteVolume import volScalarField
nbrField = solidWallMixedTemperatureCoupledFvPatchScalarField.ext_refCast( volScalarField.ext_lookupPatchField( nbrPatch, self.neighbourFieldName_ ) )
#Swap to obtain full local values of neighbour internal field
nbrIntFld = nbrField.patchInternalField()
from Foam.OpenFOAM import mapDistribute, Pstream
mapDistribute.distribute( Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrIntFld() )
# Swap to obtain full local values of neighbour K*delta
nbrKDelta = nbrField.K()*nbrPatch.deltaCoeffs()
mapDistribute.distribute( Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrKDelta() )
myKDelta = self.K()*self.patch().deltaCoeffs()
# Calculate common wall temperature. Reuse *this to store common value.
Twall = ( myKDelta*intFld + nbrKDelta*nbrIntFld ) / ( myKDelta + nbrKDelta )
# Assign to me
from Foam.finiteVolume import fvPatchScalarField
fvPatchScalarField.ext_assign( self, Twall )
mapDistribute.distribute( Pstream.defaultCommsType.fget(),
distMap.schedule(),
nbrField.size(),
distMap.constructMap(), # reverse : what to send
distMap.subMap(),
Twall() )
fvPatchScalarField.ext_assign( nbrField, Twall )
pass
# Switch between fixed value (of harmonic avg) or gradient
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
nFixed = 0
#Like snGrad but bypass switching on refValue/refGrad.
normalGradient = ( self-intFld )*self.patch().deltaCoeffs()
if self.debug:
Q = ( self.K() * self.patch().magSf() * normalGradient() ).gSum()
ext_Info ()<< "solidWallMixedTemperatureCoupledFvPatchScalarField::" << "updateCoeffs() :"\
<< " patch:" << self.patch().name()<< " heatFlux:" << Q << " walltemperature "\
<< " min:" << self.gMin() << " max:" << self.gMax() << " avg:" << self.gAverage() << nl
pass
for i in range( self.size() ):
# if outgoing flux use fixed value.
if normalGradient()[i] < 0.0:
self.refValue()[i] = self[i]
self.refGrad()[i] = 0.0 # not used
self.valueFraction()[i] = 1.0
nFixed+=1
pass
else:
self.refValue()[i]= 0.0 # not used
self.refGrad()[i] = normalGradient()[i]
self.valueFraction()[i] = 0.0
pass
from Foam.OpenFOAM import ext_reduce, sumOp_label
nFixed = ext_reduce( nFixed, sumOp_label() )
self.fixesValue_ = ( nFixed > 0 )
if (self.debug):
from Foam.OpenFOAM import returnReduce
nTotSize = returnReduce(self.size(), sumOp_label())
ext_Info() << "solidWallMixedTemperatureCoupledFvPatchScalarField::" << "updateCoeffs() :" \
<< " patch:" << self.patch().name() << " out of:" << nTotSize << " fixedBC:" << nFixed \
<< " gradient:" << nTotSize-nFixed << nl
pass
mixedFvPatchScalarField.updateCoeffs( self )
pass
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc
def K( self ) : from Foam.finiteVolume import volScalarField return volScalarField.ext_lookupPatchField( self.patch(), self.KName_ )
def updateCoeffs(self):
try:
if self.updated():
return
from Foam.meshTools import directMappedPatchBase
mpp = directMappedPatchBase.ext_refCast(self.patch().patch())
nbrMesh = mpp.sampleMesh()
from Foam.finiteVolume import fvMesh
nbrPatch = fvMesh.ext_refCast(nbrMesh).boundary()[
mpp.samplePolyPatch().index()]
# Force recalculation of mapping and schedule
distMap = mpp.map()
intFld = self.patchInternalField()
from Foam.finiteVolume import volScalarField
nbrField = solidWallMixedTemperatureCoupledFvPatchScalarField.ext_refCast(
volScalarField.ext_lookupPatchField(nbrPatch,
self.neighbourFieldName_))
#Swap to obtain full local values of neighbour internal field
nbrIntFld = nbrField.patchInternalField()
from Foam.OpenFOAM import mapDistribute, Pstream
mapDistribute.distribute(
Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrIntFld())
# Swap to obtain full local values of neighbour K*delta
nbrKDelta = nbrField.K() * nbrPatch.deltaCoeffs()
mapDistribute.distribute(
Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrKDelta())
myKDelta = self.K() * self.patch().deltaCoeffs()
# 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())
self.refValue().ext_assign(nbrIntFld)
self.refGrad().ext_assign(0.0)
self.valueFraction().ext_assign(nbrKDelta /
(nbrKDelta + myKDelta()))
mixedFvPatchScalarField.updateCoeffs(self)
if self.debug:
Q = (self.K() * self.patch().magSf() * self.snGrad()).gSum()
ext_Info ()<< patch().boundaryMesh().mesh().name() << ":" \
<< patch().name() << ':' << self.dimensionedInternalField().name() << " <- " \
<< nbrMesh.name() << ':' << nbrPatch.name() << ':' << self.dimensionedInternalField().name() << " :" \
<< " heat[W]:" << Q << " walltemperature "\
<< " min:" << self.gMin() << " max:" << self.gMax() << " avg:" << self.gAverage() << nl
pass
pass
except Exception, exc:
import sys, traceback
traceback.print_exc(file=sys.stdout)
raise exc
def updateCoeffs( self ) :
try:
if self.updated() :
return
from Foam.meshTools import directMappedPatchBase
mpp = directMappedPatchBase.ext_refCast( self.patch().patch() )
nbrMesh = mpp.sampleMesh()
from Foam.finiteVolume import fvMesh
nbrPatch = fvMesh.ext_refCast( nbrMesh ).boundary()[ mpp.samplePolyPatch().index() ]
# Force recalculation of mapping and schedule
distMap = mpp.map()
intFld = self.patchInternalField()
from Foam.finiteVolume import volScalarField
nbrField = solidWallMixedTemperatureCoupledFvPatchScalarField.ext_refCast( volScalarField.ext_lookupPatchField( nbrPatch, self.neighbourFieldName_ ) )
#Swap to obtain full local values of neighbour internal field
nbrIntFld = nbrField.patchInternalField()
from Foam.OpenFOAM import mapDistribute, Pstream
mapDistribute.distribute( Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrIntFld() )
# Swap to obtain full local values of neighbour K*delta
nbrKDelta = nbrField.K()*nbrPatch.deltaCoeffs()
mapDistribute.distribute( Pstream.defaultCommsType.fget(),
distMap.schedule(),
distMap.constructSize(),
distMap.subMap(), #what to send
distMap.constructMap(), #what to receive
nbrKDelta() )
myKDelta = self.K()*self.patch().deltaCoeffs()
# 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())
self.refValue().ext_assign( nbrIntFld )
self.refGrad().ext_assign( 0.0 )
self.valueFraction().ext_assign( nbrKDelta / ( nbrKDelta + myKDelta() ) )
mixedFvPatchScalarField.updateCoeffs( self )
if self.debug:
Q = ( self.K() * self.patch().magSf() * self.snGrad() ).gSum()
ext_Info ()<< patch().boundaryMesh().mesh().name() << ":" \
<< patch().name() << ':' << self.dimensionedInternalField().name() << " <- " \
<< nbrMesh.name() << ':' << nbrPatch.name() << ':' << self.dimensionedInternalField().name() << " :" \
<< " heat[W]:" << Q << " walltemperature "\
<< " min:" << self.gMin() << " max:" << self.gMax() << " avg:" << self.gAverage() << nl
pass
pass
except Exception, exc:
import sys, traceback
traceback.print_exc( file = sys.stdout )
raise exc