def compressibleCreatePhi(runTime, mesh, rhoU):
from Foam.OpenFOAM import IOobject, word, fileName
phiHeader = IOobject(word("phi"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ)
from Foam.OpenFOAM import ext_Info, nl
if phiHeader.headerOk():
ext_Info() << "Reading face flux field phi\n" << nl
from Foam.finiteVolume import surfaceScalarField
phi = surfaceScalarField(
IOobject(word("phi"), fileName(runTime.timeName()), mesh,
IOobject.MUST_READ, IOobject.AUTO_WRITE), mesh)
pass
else:
ext_Info() << "Calculating face flux field phi\n" << nl
from Foam.OpenFOAM import wordList
from Foam.finiteVolume import calculatedFvPatchScalarField
phiTypes = wordList(2, calculatedFvPatchScalarField.typeName)
from Foam.finiteVolume import surfaceScalarField, linearInterpolate
phi = surfaceScalarField(
IOobject(word("phi"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ, IOobject.AUTO_WRITE),
linearInterpolate(rhoU) & mesh.Sf(), phiTypes)
pass
return phiHeader, phi, phiTypes
def createPhi(runTime, hU, mesh):
from Foam.OpenFOAM import ext_Info, nl
from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName
from Foam.finiteVolume import surfaceScalarField
ext_Info() << "Reading/calculating face flux field phi\n" << nl
from Foam.finiteVolume import linearInterpolate
phi = surfaceScalarField(
IOobject(word("phi"), fileName(runTime.timeName()), mesh,
IOobject.READ_IF_PRESENT, IOobject.AUTO_WRITE),
linearInterpolate(hU) & mesh.Sf())
return phi
def createPhi( runTime, hU, mesh ):
from Foam.OpenFOAM import ext_Info, nl
from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName
from Foam.finiteVolume import surfaceScalarField
ext_Info() << "Reading/calculating face flux field phi\n" << nl
from Foam.finiteVolume import linearInterpolate
phi = surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.READ_IF_PRESENT,
IOobject.AUTO_WRITE ),
linearInterpolate( hU ) & mesh.Sf() )
return phi
def createPhi(runTime, mesh, U):
print "Reading/calculating face flux field phi\n"
from Foam.OpenFOAM import Time
from Foam.OpenFOAM import word
from Foam.OpenFOAM import fileName
from Foam.OpenFOAM import IOobject
from Foam.finiteVolume import fvMesh
from Foam.finiteVolume import volScalarField
from Foam.finiteVolume import surfaceScalarField
from Foam.finiteVolume import linearInterpolate
phi = surfaceScalarField(
IOobject(word("phi"), fileName(runTime.timeName()), mesh,
IOobject.READ_IF_PRESENT, IOobject.AUTO_WRITE),
linearInterpolate(U) & mesh.Sf())
return phi
def createPhi( runTime, mesh, U ):
print "Reading/calculating face flux field phi\n"
from Foam.OpenFOAM import Time
from Foam.OpenFOAM import word
from Foam.OpenFOAM import fileName
from Foam.OpenFOAM import IOobject
from Foam.finiteVolume import fvMesh
from Foam.finiteVolume import volScalarField
from Foam.finiteVolume import surfaceScalarField
from Foam.finiteVolume import linearInterpolate
phi = surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.READ_IF_PRESENT,
IOobject.AUTO_WRITE ),
linearInterpolate( U ) & mesh.Sf() )
return phi
def compressibleCreatePhi( runTime, mesh, rhoU ):
from Foam.OpenFOAM import IOobject, word, fileName
phiHeader = IOobject( word( "phi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ )
from Foam.OpenFOAM import ext_Info, nl
if phiHeader.headerOk():
ext_Info() << "Reading face flux field phi\n" << nl
from Foam.finiteVolume import surfaceScalarField
phi = surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
mesh )
pass
else:
ext_Info() << "Calculating face flux field phi\n" << nl
from Foam.OpenFOAM import wordList
from Foam.finiteVolume import calculatedFvPatchScalarField
phiTypes =wordList( 2, calculatedFvPatchScalarField.typeName )
from Foam.finiteVolume import surfaceScalarField, linearInterpolate
phi = surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
linearInterpolate( rhoU ) & mesh.Sf(),
phiTypes )
pass
return phiHeader, phi, phiTypes
def createFluidFields( fluidRegions, runTime ) :
# Load boundary conditions
from .. import derivedFvPatchFields
# Initialise fluid field pointer lists
from Foam.finiteVolume import PtrList_volScalarField
rhoFluid = PtrList_volScalarField( fluidRegions.size() )
KFluid = PtrList_volScalarField( fluidRegions.size() )
from Foam.finiteVolume import PtrList_volVectorField
UFluid = PtrList_volVectorField( fluidRegions.size() )
from Foam.finiteVolume import PtrList_surfaceScalarField
phiFluid = PtrList_surfaceScalarField( fluidRegions.size() )
DpDtFluid = PtrList_volScalarField( fluidRegions.size() )
from Foam.OpenFOAM import PtrList_uniformDimensionedVectorField
gFluid = PtrList_uniformDimensionedVectorField( fluidRegions.size() )
from Foam.compressible import PtrList_compressible_turbulenceModel
turbulence = PtrList_compressible_turbulenceModel( fluidRegions.size() )
from Foam.thermophysicalModels import PtrList_basicRhoThermo
thermoFluid = PtrList_basicRhoThermo( fluidRegions.size() )
p_rghFluid = PtrList_volScalarField( fluidRegions.size() )
ghFluid = PtrList_volScalarField( fluidRegions.size() )
ghfFluid = PtrList_surfaceScalarField(fluidRegions.size())
from Foam.OpenFOAM import scalarList
initialMassFluid = scalarList( fluidRegions.size() )
#Populate fluid field pointer lists
for index in range( fluidRegions.size() ) :
from Foam.OpenFOAM import ext_Info, nl
ext_Info() << "*** Reading fluid mesh thermophysical properties for region " \
<< fluidRegions[ index ].name() << nl << nl
ext_Info()<< " Adding to thermoFluid\n" << nl
from Foam.thermophysicalModels import autoPtr_basicRhoThermo, basicRhoThermo
thermo= basicRhoThermo.New( fluidRegions[ index ] )
thermoFluid.ext_set( index, thermo )
ext_Info()<< " Adding to rhoFluid\n" << nl
from Foam.OpenFOAM import word, fileName, IOobject
from Foam.finiteVolume import volScalarField
rhoFluid.ext_set( index, volScalarField( IOobject( word( "rho" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
thermoFluid[ index ].rho() ) )
ext_Info()<< " Adding to KFluid\n" << nl
KFluid.ext_set( index, volScalarField( IOobject( word( "K" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.NO_READ,
IOobject.NO_WRITE ),
thermoFluid[ index ].Cp().ptr() * thermoFluid[ index ].alpha() ) )
ext_Info()<< " Adding to UFluid\n" << nl
from Foam.finiteVolume import volVectorField
UFluid.ext_set( index, volVectorField( IOobject( word( "U" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
fluidRegions[ index ] ) )
ext_Info()<< " Adding to phiFluid\n" << nl
from Foam.finiteVolume import surfaceScalarField
from Foam.finiteVolume import linearInterpolate
phiFluid.ext_set( index, surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.READ_IF_PRESENT,
IOobject.AUTO_WRITE),
linearInterpolate( rhoFluid[ index ] * UFluid[ index ] ) & fluidRegions[ index ].Sf() ) )
ext_Info()<< " Adding to gFluid\n" << nl
from Foam.OpenFOAM import uniformDimensionedVectorField
gFluid.ext_set( index, uniformDimensionedVectorField( IOobject( word( "g" ),
fileName( runTime.constant() ),
fluidRegions[ index ],
IOobject.MUST_READ,
IOobject.NO_WRITE ) ) )
ext_Info()<< " Adding to turbulence\n" << nl
from Foam import compressible
turbulence.ext_set( index, compressible.turbulenceModel.New( rhoFluid[ index ],
UFluid[ index ],
phiFluid[ index ],
thermoFluid[ index ] ) )
ext_Info() << " Adding to ghFluid\n" << nl
ghFluid.ext_set( index, volScalarField( word( "gh" ) , gFluid[ index ] & fluidRegions[ index ].C() ) )
ext_Info() << " Adding to ghfFluid\n" << nl
ghfFluid.ext_set( index, surfaceScalarField( word( "ghf" ), gFluid[ index ] & fluidRegions[ index ].Cf() ) )
p_rghFluid.ext_set( index, volScalarField( IOobject( word( "p_rgh" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
fluidRegions[ index ] ) )
# Force p_rgh to be consistent with p
p_rghFluid[ index ].ext_assign( thermoFluid[ index ].p() - rhoFluid[ index ] * ghFluid[ index ] )
from Foam import fvc
initialMassFluid[ index ] = fvc.domainIntegrate( rhoFluid[ index ] ).value()
ext_Info()<< " Adding to DpDtFluid\n" << nl
DpDtFluid.ext_set( index, volScalarField( word( "DpDt" ), fvc.DDt( surfaceScalarField( word( "phiU" ),
phiFluid[ index ] / fvc.interpolate( rhoFluid[ index ] ) ),
thermoFluid[ index ].p() ) ) )
return thermoFluid, rhoFluid, KFluid, UFluid, phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid, ghFluid, ghfFluid, p_rghFluid
def createFluidFields(fluidRegions, runTime):
# Load boundary conditions
from .. import derivedFvPatchFields
# Initialise fluid field pointer lists
from Foam.finiteVolume import PtrList_volScalarField
rhoFluid = PtrList_volScalarField(fluidRegions.size())
KFluid = PtrList_volScalarField(fluidRegions.size())
from Foam.finiteVolume import PtrList_volVectorField
UFluid = PtrList_volVectorField(fluidRegions.size())
from Foam.finiteVolume import PtrList_surfaceScalarField
phiFluid = PtrList_surfaceScalarField(fluidRegions.size())
DpDtFluid = PtrList_volScalarField(fluidRegions.size())
from Foam.OpenFOAM import PtrList_uniformDimensionedVectorField
gFluid = PtrList_uniformDimensionedVectorField(fluidRegions.size())
from Foam.compressible import PtrList_compressible_turbulenceModel
turbulence = PtrList_compressible_turbulenceModel(fluidRegions.size())
from Foam.thermophysicalModels import PtrList_basicPsiThermo
thermoFluid = PtrList_basicPsiThermo(fluidRegions.size())
p_rghFluid = PtrList_volScalarField(fluidRegions.size())
ghFluid = PtrList_volScalarField(fluidRegions.size())
ghfFluid = PtrList_surfaceScalarField(fluidRegions.size())
from Foam.OpenFOAM import scalarList
initialMassFluid = scalarList(fluidRegions.size())
#Populate fluid field pointer lists
for index in range(fluidRegions.size()):
from Foam.OpenFOAM import ext_Info, nl
ext_Info() << "*** Reading fluid mesh thermophysical properties for region " \
<< fluidRegions[ index ].name() << nl << nl
ext_Info() << " Adding to thermoFluid\n" << nl
from Foam.thermophysicalModels import autoPtr_basicPsiThermo, basicPsiThermo
thermo = basicPsiThermo.New(fluidRegions[index])
thermoFluid.ext_set(index, thermo)
ext_Info() << " Adding to rhoFluid\n" << nl
from Foam.OpenFOAM import word, fileName, IOobject
from Foam.finiteVolume import volScalarField
rhoFluid.ext_set(
index,
volScalarField(
IOobject(word("rho"), fileName(runTime.timeName()),
fluidRegions[index], IOobject.NO_READ,
IOobject.AUTO_WRITE), thermoFluid[index].rho()))
ext_Info() << " Adding to KFluid\n" << nl
KFluid.ext_set(
index,
volScalarField(
IOobject(word("K"), fileName(runTime.timeName()),
fluidRegions[index], IOobject.NO_READ,
IOobject.NO_WRITE),
thermoFluid[index].Cp().ptr() * thermoFluid[index].alpha()))
ext_Info() << " Adding to UFluid\n" << nl
from Foam.finiteVolume import volVectorField
UFluid.ext_set(
index,
volVectorField(
IOobject(word("U"), fileName(runTime.timeName()),
fluidRegions[index], IOobject.MUST_READ,
IOobject.AUTO_WRITE), fluidRegions[index]))
ext_Info() << " Adding to phiFluid\n" << nl
from Foam.finiteVolume import surfaceScalarField
from Foam.finiteVolume import linearInterpolate
phiFluid.ext_set(
index,
surfaceScalarField(
IOobject(word("phi"), fileName(runTime.timeName()),
fluidRegions[index], IOobject.READ_IF_PRESENT,
IOobject.AUTO_WRITE),
linearInterpolate(rhoFluid[index] * UFluid[index])
& fluidRegions[index].Sf()))
ext_Info() << " Adding to gFluid\n" << nl
from Foam.OpenFOAM import uniformDimensionedVectorField
gFluid.ext_set(
index,
uniformDimensionedVectorField(
IOobject(word("g"), fileName(runTime.constant()),
fluidRegions[index], IOobject.MUST_READ,
IOobject.NO_WRITE)))
ext_Info() << " Adding to turbulence\n" << nl
from Foam import compressible
turbulence.ext_set(
index,
compressible.turbulenceModel.New(rhoFluid[index], UFluid[index],
phiFluid[index],
thermoFluid[index]))
ext_Info() << " Adding to ghFluid\n" << nl
ghFluid.ext_set(
index,
volScalarField(word("gh"),
gFluid[index] & fluidRegions[index].C()))
ext_Info() << " Adding to ghfFluid\n" << nl
ghfFluid.ext_set(
index,
surfaceScalarField(word("ghf"),
gFluid[index] & fluidRegions[index].Cf()))
p_rghFluid.ext_set(
index,
volScalarField(
IOobject(word("p_rgh"), fileName(runTime.timeName()),
fluidRegions[index], IOobject.MUST_READ,
IOobject.AUTO_WRITE), fluidRegions[index]))
# Force p_rgh to be consistent with p
p_rghFluid[index].ext_assign(thermoFluid[index].p() -
rhoFluid[index] * ghFluid[index])
from Foam import fvc
initialMassFluid[index] = fvc.domainIntegrate(rhoFluid[index]).value()
ext_Info() << " Adding to DpDtFluid\n" << nl
DpDtFluid.ext_set(
index,
volScalarField(
word("DpDt"),
fvc.DDt(
surfaceScalarField(
word("phiU"),
phiFluid[index] / fvc.interpolate(rhoFluid[index])),
thermoFluid[index].p())))
return thermoFluid, rhoFluid, KFluid, UFluid, phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid, ghFluid, ghfFluid, p_rghFluid
def _createFields(runTime, mesh, R, Cv):
from Foam.OpenFOAM import ext_Info, nl
ext_Info() << "Reading field p\n" << nl
from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName
from Foam.finiteVolume import volScalarField
p = volScalarField(
IOobject(word("p"), fileName(runTime.timeName()), mesh,
IOobject.MUST_READ, IOobject.AUTO_WRITE), mesh)
p.oldTime()
ext_Info() << "Reading field T\n" << nl
T = volScalarField(
IOobject(word("T"), fileName(runTime.timeName()), mesh,
IOobject.MUST_READ, IOobject.AUTO_WRITE), mesh)
T.correctBoundaryConditions()
psi = volScalarField(
IOobject(word("psi"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ, IOobject.AUTO_WRITE), 1.0 / (R * T))
psi.oldTime()
pbf, rhoBoundaryTypes = _rhoBoundaryTypes(p)
rho = volScalarField(
IOobject(word("rho"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ, IOobject.AUTO_WRITE), p * psi,
rhoBoundaryTypes)
ext_Info() << "Reading field U\n" << nl
from Foam.finiteVolume import volVectorField
U = volVectorField(
IOobject(word("U"), fileName(runTime.timeName()), mesh,
IOobject.MUST_READ, IOobject.AUTO_WRITE), mesh)
Ubf, rhoUboundaryTypes = _rhoUboundaryTypes(U)
rhoU = volVectorField(
IOobject(word("rhoU"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ, IOobject.AUTO_WRITE), rho * U,
rhoUboundaryTypes)
Tbf, rhoEboundaryTypes = _rhoEboundaryTypes(T)
rhoE = volScalarField(
IOobject(word("rhoE"), fileName(runTime.timeName()), mesh,
IOobject.NO_READ, IOobject.AUTO_WRITE),
rho * Cv * T + 0.5 * rho * (rhoU / rho).magSqr(), rhoEboundaryTypes)
phiHeader, phi, phiTypes = compressibleCreatePhi(runTime, mesh, rhoU)
phi.oldTime()
from Foam.finiteVolume import surfaceScalarField, linearInterpolate
phiv = surfaceScalarField(
IOobject(word("phiv"), fileName(runTime.timeName()), mesh),
phi / linearInterpolate(rho),
phi.ext_boundaryField().types())
rhoU.correctBoundaryConditions()
from Foam.finiteVolume import TfieldTable_scalar
fields = TfieldTable_scalar()
magRhoU = rhoU.mag()
H = volScalarField(word("H"), (rhoE + p) / rho)
fields.add(rho)
fields.add(magRhoU)
fields.add(H)
return p, T, psi, pbf, rhoBoundaryTypes, rho, U, Ubf, rhoUboundaryTypes, rhoU, Tbf, rhoEboundaryTypes, rhoE, phi, phiv, rhoU, fields, magRhoU, H
def _createFields( runTime, mesh, R, Cv ):
from Foam.OpenFOAM import ext_Info, nl
ext_Info() << "Reading field p\n" << nl
from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName
from Foam.finiteVolume import volScalarField
p = volScalarField( IOobject( word( "p" ),
fileName( runTime.timeName() ),
mesh,
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
mesh )
p.oldTime()
ext_Info() << "Reading field T\n" << nl
T = volScalarField( IOobject( word( "T" ),
fileName( runTime.timeName() ),
mesh,
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
mesh )
T.correctBoundaryConditions()
psi = volScalarField( IOobject( word( "psi" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
1.0 / ( R * T ) )
psi.oldTime()
pbf, rhoBoundaryTypes = _rhoBoundaryTypes( p )
rho = volScalarField( IOobject( word( "rho" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
p * psi,
rhoBoundaryTypes )
ext_Info() << "Reading field U\n" << nl
from Foam.finiteVolume import volVectorField
U = volVectorField( IOobject( word( "U" ),
fileName( runTime.timeName() ),
mesh,
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
mesh )
Ubf, rhoUboundaryTypes = _rhoUboundaryTypes( U )
rhoU = volVectorField( IOobject( word( "rhoU" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
rho * U,
rhoUboundaryTypes )
Tbf, rhoEboundaryTypes = _rhoEboundaryTypes( T )
rhoE = volScalarField( IOobject( word( "rhoE" ),
fileName( runTime.timeName() ),
mesh,
IOobject.NO_READ,
IOobject.AUTO_WRITE ),
rho * Cv * T + 0.5 * rho * ( rhoU / rho ).magSqr(),
rhoEboundaryTypes )
phiHeader, phi, phiTypes = compressibleCreatePhi( runTime, mesh, rhoU )
phi.oldTime()
from Foam.finiteVolume import surfaceScalarField, linearInterpolate
phiv = surfaceScalarField( IOobject( word( "phiv" ),
fileName( runTime.timeName() ),
mesh ),
phi / linearInterpolate( rho ),
phi.ext_boundaryField().types() )
rhoU.correctBoundaryConditions()
from Foam.finiteVolume import TfieldTable_scalar
fields = TfieldTable_scalar()
magRhoU = rhoU.mag()
H = volScalarField( word( "H" ) , ( rhoE + p ) / rho )
fields.add( rho )
fields.add( magRhoU )
fields.add( H )
return p, T, psi, pbf, rhoBoundaryTypes, rho, U, Ubf, rhoUboundaryTypes, rhoU, Tbf, rhoEboundaryTypes, rhoE, phi, phiv, rhoU, fields, magRhoU, H
ext_Info() << " Adding to Uf\n" << nl
Uf.ext_set( index, volVectorField( IOobject( word( "U" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.MUST_READ,
IOobject.AUTO_WRITE ),
fluidRegions[ index ] ) )
ext_Info() << " Adding to phif\n" << nl
from Foam.finiteVolume import linearInterpolate
phif.ext_set( index, surfaceScalarField( IOobject( word( "phi" ),
fileName( runTime.timeName() ),
fluidRegions[ index ],
IOobject.READ_IF_PRESENT,
IOobject.AUTO_WRITE ),
linearInterpolate( rhof[ index ]*Uf[ index ] ) & fluidRegions[ index ].Sf() ) )
ext_Info() << " Adding to turb\n" << nl
from Foam import compressible
turb.ext_set( index, compressible.RASModel.New( rhof[ index ], Uf[ index ], phif[ index ], thermof[ index ] ) )
ext_Info() << " Adding to DpDtf\n" << nl
from Foam import fvc
DpDtf.ext_set( index, volScalarField( fvc.DDt( surfaceScalarField( word( "phiU" ),
phif[ index ] / fvc.interpolate( rhof[ index ] ) ),
thermof[ index ].p() ) ) )
from Foam.OpenFOAM import IOdictionary
environmentalProperties = IOdictionary.ext_lookupObject( fluidRegions[ index ], word( "environmentalProperties" ) )
