def createFields( runTime, mesh ): ext_Info() << "Reading field U\n" << nl from Foam.finiteVolume import volVectorField from Foam.OpenFOAM import IOobject, fileName, word U = volVectorField( IOobject( word( "U" ), fileName( runTime.timeName() ), mesh, IOobject.MUST_READ, IOobject.AUTO_WRITE ), mesh ) from Foam.finiteVolume import volSymmTensorField from Foam.OpenFOAM import dimensionedSymmTensor, symmTensor from Foam.OpenFOAM import dimForce, dimArea sigma = volSymmTensorField( IOobject( word( "sigma" ), fileName( runTime.timeName() ), mesh, IOobject.READ_IF_PRESENT, IOobject.AUTO_WRITE ), mesh, dimensionedSymmTensor( word( "zero" ), dimForce/dimArea, symmTensor.zero) ) from materialModels.rheologyModel import rheologyModel rheology = rheologyModel( sigma ) return U, sigma, rheology
def createFields(runTime, mesh): ext_Info() << "Reading field U\n" << nl from Foam.finiteVolume import volVectorField from Foam.OpenFOAM import IOobject, fileName, word U = volVectorField( IOobject(word("U"), fileName(runTime.timeName()), mesh, IOobject.MUST_READ, IOobject.AUTO_WRITE), mesh) from Foam.finiteVolume import volSymmTensorField from Foam.OpenFOAM import dimensionedSymmTensor, symmTensor from Foam.OpenFOAM import dimForce, dimArea sigma = volSymmTensorField( IOobject(word("sigma"), fileName(runTime.timeName()), mesh, IOobject.READ_IF_PRESENT, IOobject.AUTO_WRITE), mesh, dimensionedSymmTensor(word("zero"), dimForce / dimArea, symmTensor.zero)) from materialModels.rheologyModel import rheologyModel rheology = rheologyModel(sigma) return U, sigma, rheology
def main_standalone(argc, argv): from Foam.OpenFOAM.include import setRootCase args = setRootCase(argc, argv) from Foam.OpenFOAM.include import createTime runTime = createTime(args) from Foam.OpenFOAM.include import createMesh mesh = createMesh(runTime) U, phi, laminarTransport, turbulence, Ubar, wallNormal, flowDirection, flowMask, y, gradP = _createFields( runTime, mesh) from Foam.OpenFOAM import ext_Info, nl ext_Info() << "\nStarting time loop\n" << nl while runTime.loop(): ext_Info() << "\nTime = " << runTime.timeName() << nl << nl divR = turbulence.divDevReff(U) from Foam.OpenFOAM import vectorField tmp = flowMask & divR.source() divR.source().ext_assign(tmp) UEqn = divR == gradP UEqn.relax() UEqn.solve() # Correct driving force for a constant mass flow rate UbarStar = flowMask & U.weightedAverage(mesh.V()) U.ext_assign(U + (Ubar - UbarStar)) gradP += (Ubar - UbarStar) / (1.0 / UEqn.A()).weightedAverage(mesh.V()) id_ = y.size() - 1 wallShearStress = flowDirection & turbulence.R()()[id_] & wallNormal from Foam.OpenFOAM import mag from math import sqrt yplusWall = sqrt( mag(wallShearStress)) * y()[id_] / turbulence.nuEff()()[id_] ext_Info() << "Uncorrected Ubar = " << ( flowDirection & UbarStar.value() )<< " " \ << "pressure gradient = " << ( flowDirection & gradP.value() ) << " " \ << "min y+ = " << yplusWall << nl turbulence.correct() if runTime.outputTime(): from Foam.finiteVolume import volSymmTensorField from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName R = volSymmTensorField( IOobject(word("R"), fileName(runTime.timeName()), mesh, IOobject.NO_READ, IOobject.AUTO_WRITE), turbulence.R()) runTime.write() gFormat = runTime.graphFormat() from Foam.sampling import makeGraph makeGraph(y, flowDirection & U, word("Uf"), gFormat) makeGraph(y, laminarTransport.ext_nu(), gFormat) makeGraph(y, turbulence.ext_k(), gFormat) makeGraph(y, turbulence.ext_epsilon(), gFormat) from Foam.OpenFOAM import tensor makeGraph(y, R.component(tensor.XY), word("uv"), gFormat) from Foam import fvc makeGraph(y, fvc.grad(U).mag(), word("gammaDot"), gFormat) runTime.write() pass ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \ " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << nl pass ext_Info() << "End\n" << nl import os return os.EX_OK
def main_standalone( argc, argv ): from Foam.OpenFOAM.include import setRootCase args = setRootCase( argc, argv ) from Foam.OpenFOAM.include import createTime runTime = createTime( args ) from Foam.OpenFOAM.include import createMesh mesh = createMesh( runTime ) U, phi, laminarTransport, turbulence, Ubar, wallNormal, flowDirection, flowMask, y, gradP = _createFields( runTime, mesh ) from Foam.OpenFOAM import ext_Info, nl ext_Info() << "\nStarting time loop\n" << nl while runTime.loop() : ext_Info() << "\nTime = " << runTime.timeName() << nl << nl divR = turbulence.divDevReff( U ) divR.source().ext_assign( flowMask & divR.source() ) UEqn = divR == gradP UEqn.relax() UEqn.solve() # Correct driving force for a constant mass flow rate UbarStar = flowMask & U.weightedAverage(mesh.V()) U.ext_assign( U + ( Ubar - UbarStar ) ) gradP += ( Ubar - UbarStar ) / ( 1.0 / UEqn.A() ).weightedAverage( mesh.V() ) id_ = y.size() - 1 wallShearStress = flowDirection & turbulence.R()()[id_] & wallNormal from Foam.OpenFOAM import mag from math import sqrt yplusWall = sqrt( mag( wallShearStress ) ) * y()[ id_ ] / turbulence.nuEff()()[ id_ ] ext_Info() << "Uncorrected Ubar = " << ( flowDirection & UbarStar.value() )<< " " \ << "pressure gradient = " << ( flowDirection & gradP.value() ) << " " \ << "min y+ = " << yplusWall << nl turbulence.correct() if runTime.outputTime(): from Foam.finiteVolume import volSymmTensorField from Foam.OpenFOAM import IOdictionary, IOobject, word, fileName R = volSymmTensorField( IOobject( word( "R" ), fileName( runTime.timeName() ), mesh, IOobject.NO_READ, IOobject.AUTO_WRITE ), turbulence.R() ) runTime.write() gFormat = runTime.graphFormat() from Foam.sampling import makeGraph makeGraph( y, flowDirection & U, word( "Uf" ), gFormat ) makeGraph( y, laminarTransport.ext_nu(), gFormat ) makeGraph( y, turbulence.ext_k(), gFormat ) makeGraph( y, turbulence.ext_epsilon(), gFormat ) from Foam.OpenFOAM import tensor makeGraph( y, R.component( tensor.XY ), word( "uv" ), gFormat ) from Foam import fvc makeGraph( y, fvc.grad(U).mag(), word( "gammaDot" ), gFormat ) runTime.write() pass ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \ " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << nl pass ext_Info() << "End\n" << nl import os return os.EX_OK