def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
pimple = man.pimpleControl(mesh)
pThermo, p, h, psi, rho, U, phi, rhoMax, rhoMin, turbulence, dpdt, K = createFields( runTime, mesh, pimple )
mrfZones, pZones, pressureImplicitPorosity = createZones( mesh, U )
cumulativeContErr = ref.initContinuityErrs()
ref.ext_Info()<< "\nStarting time loop\n" << ref.nl;
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.compressibleCourantNo( mesh, phi, rho, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info()<< "Time = " << runTime.timeName() << ref.nl << ref.nl;
ref.rhoEqn( rho, phi )
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = fun_Ueqn( pimple, rho, p, U, phi, turbulence, mrfZones, pZones )
fun_hEqn(pThermo, rho, p, h, phi, turbulence, dpdt, K )
# --- PISO loop
while (pimple.correct()):
cumulativeContErr = fun_pEqn( mesh, runTime, pimple, pThermo, rho, p, h, psi, U, phi, mrfZones,
turbulence, UEqn, dpdt, K, cumulativeContErr, rhoMax, rhoMin )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info()<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \
<< " ClockTime = " << runTime.elapsedClockTime() << " s" \
<< ref.nl << ref.nl
pass
ref.ext_Info()<< "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
p, U, phi, turbulence, pRefCell, pRefValue, laminarTransport = _createFields( runTime, mesh )
cumulativeContErr = ref.initContinuityErrs()
pimple = man.pimpleControl( mesh )
ref.ext_Info() << "\nStarting time loop\n" <<ref.nl
while runTime.run() :
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
# --- Pressure-velocity PIMPLE corrector loop
pimple.start()
while pimple.loop():
if pimple.nOuterCorr() != 1 :
p.storePrevIter()
pass
UEqn, rAU = Ueqn( mesh, pimple, phi, U, p, turbulence )
# --- PISO loop
for corr in range( pimple.nCorr() ):
cumulativeContErr = pEqn( runTime, mesh, pimple, U, rAU, UEqn, phi, p, corr, pRefCell, pRefValue, cumulativeContErr )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pimple.increment()
pass
runTime.write();
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
p, U, phi, turbulence, pRefCell, pRefValue, laminarTransport, sources = _createFields(
runTime, mesh)
cumulativeContErr = ref.initContinuityErrs()
pimple = man.pimpleControl(mesh)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
runTime = ref.setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT,
CoNum)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn, rAU = Ueqn(mesh, pimple, phi, U, p, turbulence, sources)
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = pEqn(runTime, mesh, pimple, U, rAU, UEqn,
phi, p, pRefCell, pRefValue,
cumulativeContErr, sources)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
p, U, Urel, SRF, phi, turbulence, pRefCell, pRefValue, laminarTransport, sources = _createFields( runTime, mesh )
cumulativeContErr = ref.initContinuityErrs()
pimple = man.pimpleControl( mesh )
ref.ext_Info() << "\nStarting time loop\n" <<ref.nl
while runTime.run() :
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UrelEqn = _UrelEqn( mesh, pimple, phi, Urel, p, turbulence, SRF, sources )
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = pEqn( runTime, mesh, pimple, Urel, UrelEqn, phi, p, pRefCell, pRefValue, cumulativeContErr, sources )
pass
# Update the absolute velocity
U << Urel() + SRF.U() # mixed calculations
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write();
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createDynamicFvMesh(runTime)
cumulativeContErr = ref.initContinuityErrs()
p_rgh, p, alpha1, U, phi, rho1, rho2, rho, rhoPhi, twoPhaseProperties, pRefCell, pRefValue, interface, turbulence, g, gh, ghf = _createFields(
runTime, mesh
)
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
pimple = man.pimpleControl(mesh)
phiAbs = ref.surfaceScalarField(ref.word("phiAbs"), phi)
ref.fvc.makeAbsolute(phiAbs, U)
cumulativeContErr = fun_correctPhi(
runTime, mesh, phi, phiAbs, p, p_rgh, rho, U, cumulativeContErr, pimple, pRefCell, pRefValue
)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
runTime = ref.setInitialDeltaT(runTime, adjustTimeStep, maxCo, CoNum)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT, correctPhi, checkMeshCourantNo = readControls(runTime, mesh, pimple)
maxAlphaCo, alphaCoNum, meanAlphaCoNum = alphaCourantNo(runTime, mesh, alpha1, phi)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
runTime = setDeltaT(runTime, adjustTimeStep, maxCo, CoNum, maxAlphaCo, alphaCoNum, maxDeltaT)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
timeBeforeMeshUpdate = runTime.elapsedCpuTime()
fun_Urel(mesh, U)
if mesh.changing():
ref.ext_Info() << "Execution time for mesh.update() = " << runTime.elapsedCpuTime() - timeBeforeMeshUpdate << " s" << ref.nl
gh << (g & mesh.C())
ghf << (g & mesh.Cf())
pass
if mesh.changing() and correctPhi:
cumulativeContErr = fun_correctPhi(
runTime, mesh(), phi, phiAbs, p, p_rgh, rho, U, cumulativeContErr, pimple, pRefCell, pRefValue
)
pass
if mesh.changing() and checkMeshCourantNo:
meshCoNum, meanMeshCoNum = ref.meshCourantNo(runTime, mesh(), phi)
pass
twoPhaseProperties.correct()
alphaEqnSubCycle(runTime, pimple, mesh, phi, alpha1, rho, rhoPhi, rho1, rho2, interface)
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = _UEqn(
mesh(), alpha1, U, p, p_rgh, ghf, rho, rhoPhi, turbulence, g, twoPhaseProperties, interface, pimple
)
# --- PISO loop
while pimple.correct():
cumulativeContErr = _pEqn(
runTime,
mesh,
UEqn,
U,
p,
p_rgh,
gh,
ghf,
phi,
phiAbs,
alpha1,
rho,
g,
interface,
pimple,
pRefCell,
pRefValue,
cumulativeContErr,
)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
pimple = man.pimpleControl(mesh)
maxDeltaT, rDeltaT = setInitialrDeltaT(runTime, mesh, pimple)
pThermo, p, h, psi, rho, U, phi, rhoMax, rhoMin, turbulence, dpdt, K = createFields(
runTime, mesh, pimple)
mrfZones, pZones, pressureImplicitPorosity = createZones(mesh, U)
cumulativeContErr = ref.initContinuityErrs()
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.compressibleCourantNo(mesh, phi, rho, runTime)
runTime = ref.setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT,
CoNum)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
setrDeltaT(runTime, mesh, pimple, phi, psi, U, rho, rDeltaT, maxDeltaT)
ref.rhoEqn(rho, phi)
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
turbulence.correct()
UEqn = fun_Ueqn(pimple, rho, p, U, phi, turbulence, mrfZones,
pZones)
fun_hEqn(pThermo, rho, p, h, phi, turbulence, dpdt, K)
# --- PISO loop
while (pimple.correct()):
cumulativeContErr = fun_pEqn(mesh, runTime, pimple, pThermo,
rho, p, h, psi, U, phi, mrfZones,
turbulence, UEqn, dpdt, K,
cumulativeContErr, rhoMax, rhoMin)
pass
pass
runTime.write()
ref.ext_Info()<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \
<< " ClockTime = " << runTime.elapsedClockTime() << " s" \
<< ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createDynamicFvMesh(runTime)
cumulativeContErr = ref.initContinuityErrs()
p, U, phi, laminarTransport, turbulence, rAU, pRefCell, pRefValue, sources = _createFields(
runTime, mesh)
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
pimple = man.pimpleControl(mesh)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT, pimpleDic, correctPhi, checkMeshCourantNo, ddtPhiCorr = readControls(
runTime, mesh, pimple)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
# Make the fluxes absolute
ref.fvc.makeAbsolute(phi, U)
runTime = ref.setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT,
CoNum)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
mesh.update()
if mesh.changing() and correctPhi:
cumulativeContErr = _correctPhi(runTime, mesh, pimple, p, U, rAU,
phi, pRefCell, pRefValue,
cumulativeContErr)
pass
# Make the fluxes relative to the mesh motion
ref.fvc.makeRelative(phi, U)
if mesh.changing() and checkMeshCourantNo:
meshCoNum, meanMeshCoNum = ref.meshCourantNo(runTime, mesh, phi)
pass
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = fun_UEqn(mesh, phi, U, p, rAU, turbulence, pimple)
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = fun_pEqn(mesh, runTime, pimple, U, phi,
turbulence, p, rAU, UEqn,
pRefCell, pRefValue,
cumulativeContErr, ddtPhiCorr)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
pimple = man.pimpleControl( mesh )
cumulativeContErr = ref.initContinuityErrs()
p_rgh, p, alpha1, U, phi, rho1, rho2, rho, rhoPhi, twoPhaseProperties, pRefCell, \
pRefValue, interface, turbulence, g, gh, ghf = _createFields( runTime, mesh )
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
cumulativeContErr = correctPhi( runTime, mesh, phi, p, p_rgh, rho, U, cumulativeContErr, pimple, pRefCell, pRefValue)
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setInitialDeltaT( runTime, adjustTimeStep, maxCo, CoNum )
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run() :
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
maxAlphaCo, alphaCoNum, meanAlphaCoNum = alphaCourantNo( runTime, mesh, alpha1, phi )
runTime = setDeltaT( runTime, adjustTimeStep, maxCo, CoNum, maxAlphaCo, alphaCoNum, maxDeltaT )
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
twoPhaseProperties.correct()
alphaEqnSubCycle( runTime, pimple, mesh, phi, alpha1, rho, rhoPhi, rho1, rho2, interface )
while pimple.loop():
UEqn = _UEqn( mesh, alpha1, U, p, p_rgh, ghf, rho, rhoPhi, turbulence, g, twoPhaseProperties, interface, pimple )
# --- PISO loop
while pimple.correct():
cumulativeContErr = _pEqn( runTime, mesh, UEqn, U, p, p_rgh, gh, ghf, phi, alpha1, rho, g,
interface, pimple, pRefCell, pRefValue, cumulativeContErr )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createDynamicFvMesh( runTime )
cumulativeContErr = ref.initContinuityErrs()
p, U, phi, laminarTransport, turbulence, rAU, pRefCell, pRefValue, sources = _createFields( runTime, mesh )
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
pimple = man.pimpleControl( mesh )
ref.ext_Info() << "\nStarting time loop\n" <<ref.nl
while runTime.run() :
adjustTimeStep, maxCo, maxDeltaT, pimpleDic, correctPhi, checkMeshCourantNo, ddtPhiCorr = readControls( runTime, mesh, pimple )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
# Make the fluxes absolute
ref.fvc.makeAbsolute(phi, U)
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
mesh.update()
if mesh.changing() and correctPhi :
cumulativeContErr = _correctPhi( runTime, mesh, pimple, p, U, rAU, phi, pRefCell, pRefValue, cumulativeContErr )
pass
# Make the fluxes relative to the mesh motion
ref.fvc.makeRelative( phi, U )
if mesh.changing() and checkMeshCourantNo :
meshCoNum, meanMeshCoNum = ref.meshCourantNo( runTime, mesh, phi )
pass
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = fun_UEqn( mesh, phi, U, p, rAU, turbulence, pimple, sources )
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = fun_pEqn( mesh, runTime, pimple, U, phi, turbulence, p, rAU, UEqn, pRefCell, pRefValue, cumulativeContErr, ddtPhiCorr, sources )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
pThermo, p, h, psi, rho, U, phi, rhoMax, rhoMin, turbulence, DpDt = createFields( runTime, mesh )
cumulativeContErr = ref.initContinuityErrs()
pimple = man.pimpleControl(mesh)
ref.ext_Info()<< "\nStarting time loop\n" << ref.nl;
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.compressibleCourantNo( mesh, phi, rho, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info()<< "Time = " << runTime.timeName() << ref.nl << ref.nl;
ref.rhoEqn( rho, phi )
# --- Pressure-velocity PIMPLE corrector loop
pimple.start()
while pimple.loop():
if pimple.nOuterCorr() != 1:
p.storePrevIter()
rho.storePrevIter()
pass
UEqn, rAU = fun_Ueqn( pimple, rho, p, U, phi, turbulence )
fun_hEqn( pThermo, rho, p, h, phi, turbulence, DpDt )
# --- PISO loop
for corr in range( pimple.nCorr() ):
cumulativeContErr = fun_pEqn( mesh, runTime, pimple, pThermo, rho, p, h, psi, U, phi,
turbulence, UEqn, rAU, DpDt, cumulativeContErr, corr, rhoMax, rhoMin )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pimple.increment()
pass
runTime.write()
ref.ext_Info()<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \
<< " ClockTime = " << runTime.elapsedClockTime() << " s" \
<< ref.nl << ref.nl
pass
ref.ext_Info()<< "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
g = ref.readGravitationalAcceleration(runTime, mesh)
thermo, p, rho, h, psi, U, phi, turbulence, gh, ghf, p_rgh, dpdt, K = createFields(
runTime, mesh, g)
cumulativeContErr = ref.initContinuityErrs()
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.compressibleCourantNo(mesh, phi, rho, runTime)
runTime = ref.setInitialDeltaT(runTime, adjustTimeStep, maxCo, CoNum)
pimple = man.pimpleControl(mesh)
# * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.compressibleCourantNo(mesh, phi, rho, runTime)
runTime = ref.setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT,
CoNum)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
ref.rhoEqn(rho, phi)
# --- Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = fun_Ueqn(pimple, mesh, rho, U, phi, turbulence, ghf, p_rgh,
K)
fun_hEqn(thermo, rho, p, h, phi, turbulence, dpdt, K)
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = fun_pEqn(mesh, runTime, pimple, thermo,
rho, p, h, psi, U, phi,
turbulence, gh, ghf, p_rgh, UEqn,
dpdt, K, cumulativeContErr)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
rho << thermo.rho()
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \
<< " ClockTime = " << runTime.elapsedClockTime() << " s" \
<< ref.nl << ref.nl
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
g = ref.readGravitationalAcceleration(runTime, mesh)
thermo, p, rho, h, psi, U, phi, turbulence, gh, ghf, p_rgh, DpDt = createFields(runTime, mesh, g)
cumulativeContErr = ref.initContinuityErrs()
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.compressibleCourantNo(mesh, phi, rho, runTime)
runTime = ref.setInitialDeltaT(runTime, adjustTimeStep, maxCo, CoNum)
pimple = man.pimpleControl(mesh)
# * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.compressibleCourantNo(mesh, phi, rho, runTime)
runTime = ref.setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
ref.rhoEqn(rho, phi)
# --- Pressure-velocity PIMPLE corrector loop
pimple.start()
while pimple.loop():
if pimple.nOuterCorr() != 1:
p_rgh.storePrevIter()
pass
UEqn = fun_Ueqn(pimple, mesh, rho, U, phi, turbulence, ghf, p_rgh)
fun_hEqn(thermo, rho, p, h, phi, turbulence, DpDt)
# --- PISO loop
for corr in range(pimple.nCorr()):
cumulativeContErr = fun_pEqn(
mesh,
runTime,
pimple,
thermo,
rho,
p,
h,
psi,
U,
phi,
turbulence,
gh,
ghf,
p_rgh,
UEqn,
DpDt,
cumulativeContErr,
corr,
)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pimple.increment()
pass
rho << thermo.rho()
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime = man.createTime( args )
mesh = man.createMesh( runTime )
g = man.readGravitationalAcceleration( runTime, mesh )
cumulativeContErr = ref.initContinuityErrs()
p_rgh, alpha1, U, phi, twoPhaseProperties, rho1, rho2, Dab, \
alphatab, rho, rhoPhi, turbulence, gh, ghf, p, pRefCell, pRefValue = _createFields( runTime, mesh, g )
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setInitialDeltaT( runTime, adjustTimeStep, maxCo, CoNum )
pimple = man.pimpleControl( mesh )
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
# --- Pressure-velocity PIMPLE corrector loop
pimple.start()
while pimple.loop():
twoPhaseProperties.correct()
alphaEqn( mesh, phi, alpha1, alphatab, Dab, rhoPhi, rho, rho1, rho2, turbulence )
UEqn = fun_UEqn( mesh, U, p_rgh, ghf, rho, rhoPhi, turbulence, twoPhaseProperties, pimple )
# --- PISO loop
for corr in range( pimple.nCorr() ):
cumulativeContErr = fun_pEqn( runTime, mesh, UEqn, U, p, p_rgh, gh, ghf, phi, rho, pimple, \
corr, pRefCell, pRefValue, cumulativeContErr )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pimple.increment()
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
gravitationalProperties, g, rotating, Omega, magg, gHat = readGravitationalAcceleration(
runTime, mesh)
h, h0, U, hU, hTotal, phi, F = _createFields(runTime, mesh, Omega, gHat)
pimple = man.pimpleControl(mesh)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.loop():
ref.ext_Info() << "\n Time = " << runTime.timeName(
) << ref.nl << ref.nl
CourantNo(runTime, mesh, h, phi, magg)
pimple.start()
while pimple.loop():
phiv = ref.surfaceScalarField(
ref.word("phiv"),
phi() / ref.fvc.interpolate(h)) # mixed calculations
hUEqn = ref.fvm.ddt(hU) + ref.fvm.div(phiv, hU)
hUEqn.relax()
if pimple.momentumPredictor():
if rotating:
ref.solve(hUEqn + (F ^ hU) == -magg * h *
ref.fvc.grad(h + h0))
pass
else:
ref.solve(hUEqn == -magg * h * ref.fvc.grad(h + h0))
pass
# Constrain the momentum to be in the geometry if 3D geometry
if mesh.nGeometricD() == 3:
hU -= (gHat & hU) * gHat
hU.correctBoundaryConditions()
pass
for corr in range(pimple.nCorr()):
hf = ref.fvc.interpolate(h)
rUA = 1.0 / hUEqn.A()
ghrUAf = magg * ref.fvc.interpolate(h * rUA)
phih0 = ghrUAf * mesh.magSf() * ref.fvc.snGrad(h0)
if rotating:
hU << rUA * (hUEqn.H() - (F ^ hU))
pass
else:
hU << rUA * hUEqn.H()
pass
phi << (ref.fvc.interpolate(hU) & mesh.Sf()
) + ref.fvc.ddtPhiCorr(rUA, h, hU, phi) - phih0
for nonOrth in range(pimple.nNonOrthCorr() + 1):
hEqn = ref.fvm.ddt(h) + ref.fvc.div(
phi) - ref.fvm.laplacian(ghrUAf, h)
hEqn.solve(
mesh.solver(
h.select(pimple.finalInnerIter(corr, nonOrth))))
if nonOrth == pimple.nNonOrthCorr():
phi += hEqn.flux()
pass
hU -= rUA * h * magg * ref.fvc.grad(h + h0)
#Constrain the momentum to be in the geometry if 3D geometry
if mesh.nGeometricD() == 3:
hU -= (gHat & hU) * gHat
pass
hU.correctBoundaryConditions()
pass
pimple.increment()
pass
U == hU / h
hTotal == h + h0
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
gravitationalProperties, g, rotating, Omega, magg, gHat = readGravitationalAcceleration(runTime, mesh)
h, h0, U, hU, hTotal, phi, F = _createFields(runTime, mesh, Omega, gHat)
pimple = man.pimpleControl(mesh)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.loop():
ref.ext_Info() << "\n Time = " << runTime.timeName() << ref.nl << ref.nl
CourantNo(runTime, mesh, h, phi, magg)
pimple.start()
while pimple.loop():
phiv = ref.surfaceScalarField(ref.word("phiv"), phi() / ref.fvc.interpolate(h)) # mixed calculations
hUEqn = ref.fvm.ddt(hU) + ref.fvm.div(phiv, hU)
hUEqn.relax()
if pimple.momentumPredictor():
if rotating:
ref.solve(hUEqn + (F ^ hU) == -magg * h * ref.fvc.grad(h + h0))
pass
else:
ref.solve(hUEqn == -magg * h * ref.fvc.grad(h + h0))
pass
# Constrain the momentum to be in the geometry if 3D geometry
if mesh.nGeometricD() == 3:
hU -= (gHat & hU) * gHat
hU.correctBoundaryConditions()
pass
for corr in range(pimple.nCorr()):
hf = ref.fvc.interpolate(h)
rUA = 1.0 / hUEqn.A()
ghrUAf = magg * ref.fvc.interpolate(h * rUA)
phih0 = ghrUAf * mesh.magSf() * ref.fvc.snGrad(h0)
if rotating:
hU << rUA * (hUEqn.H() - (F ^ hU))
pass
else:
hU << rUA * hUEqn.H()
pass
phi << (ref.fvc.interpolate(hU) & mesh.Sf()) + ref.fvc.ddtPhiCorr(rUA, h, hU, phi) - phih0
for nonOrth in range(pimple.nNonOrthCorr() + 1):
hEqn = ref.fvm.ddt(h) + ref.fvc.div(phi) - ref.fvm.laplacian(ghrUAf, h)
hEqn.solve(mesh.solver(h.select(pimple.finalInnerIter(corr, nonOrth))))
if nonOrth == pimple.nNonOrthCorr():
phi += hEqn.flux()
pass
hU -= rUA * h * magg * ref.fvc.grad(h + h0)
# Constrain the momentum to be in the geometry if 3D geometry
if mesh.nGeometricD() == 3:
hU -= (gHat & hU) * gHat
pass
hU.correctBoundaryConditions()
pass
pimple.increment()
pass
U == hU / h
hTotal == h + h0
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone(argc, argv):
args = ref.setRootCase(argc, argv)
runTime = man.createTime(args)
mesh = man.createMesh(runTime)
pimple = man.pimpleControl(mesh)
cumulativeContErr = ref.initContinuityErrs()
p_rgh, p, alpha1, U, phi, rho1, rho2, rho, rhoPhi, twoPhaseProperties, pRefCell, \
pRefValue, interface, turbulence, g, gh, ghf = _createFields( runTime, mesh )
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
cumulativeContErr = correctPhi(runTime, mesh, phi, p, p_rgh, rho, U,
cumulativeContErr, pimple, pRefCell,
pRefValue)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
runTime = ref.setInitialDeltaT(runTime, adjustTimeStep, maxCo, CoNum)
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.run():
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls(runTime)
CoNum, meanCoNum = ref.CourantNo(mesh, phi, runTime)
maxAlphaCo, alphaCoNum, meanAlphaCoNum = alphaCourantNo(
runTime, mesh, alpha1, phi)
runTime = setDeltaT(runTime, adjustTimeStep, maxCo, CoNum, maxAlphaCo,
alphaCoNum, maxDeltaT)
runTime.increment()
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
twoPhaseProperties.correct()
alphaEqnSubCycle(runTime, pimple, mesh, phi, alpha1, rho, rhoPhi, rho1,
rho2, interface)
while pimple.loop():
UEqn = _UEqn(mesh, alpha1, U, p, p_rgh, ghf, rho, rhoPhi,
turbulence, g, twoPhaseProperties, interface, pimple)
# --- PISO loop
while pimple.correct():
cumulativeContErr = _pEqn(runTime, mesh, UEqn, U, p, p_rgh, gh,
ghf, phi, alpha1, rho, g, interface,
pimple, pRefCell, pRefValue,
cumulativeContErr)
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << \
" ClockTime = " << runTime.elapsedClockTime() << " s" << ref.nl << ref.nl
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
def main_standalone( argc, argv ):
args = ref.setRootCase( argc, argv )
runTime=man.createTime( args )
mesh = man.createMesh( runTime )
g = readGravitationalAcceleration( runTime, mesh );
T, p_rgh, U, phi, laminarTransport, turbulence, rhok, \
kappat, gh, ghf, p, pRefCell, pRefValue, beta, TRef, Pr, Prt = createFields( runTime, mesh, g )
cumulativeContErr = ref.initContinuityErrs()
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setInitialDeltaT( runTime, adjustTimeStep, maxCo, CoNum )
pimple = man.pimpleControl( mesh )
# * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * #
ref.ext_Info() << "\nStarting time loop\n" << ref.nl
while runTime.loop():
ref.ext_Info() << "Time = " << runTime.timeName() << ref.nl << ref.nl
adjustTimeStep, maxCo, maxDeltaT = ref.readTimeControls( runTime )
CoNum, meanCoNum = ref.CourantNo( mesh, phi, runTime )
runTime = ref.setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum )
# Pressure-velocity PIMPLE corrector loop
while pimple.loop():
UEqn = fun_UEqn( mesh, pimple, U, phi, turbulence, p, rhok, p_rgh, ghf )
fun_TEqn( phi, turbulence, kappat, T, rhok, beta, TRef, Prt, Pr )
# --- Pressure corrector loop
while pimple.correct():
cumulativeContErr = fun_pEqn( mesh, runTime, pimple, p, rhok, U, phi, turbulence, gh, ghf, p_rgh, UEqn, pRefCell, pRefValue, cumulativeContErr )
pass
if pimple.turbCorr():
turbulence.correct()
pass
pass
runTime.write()
ref.ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \
<< " ClockTime = " << runTime.elapsedClockTime() << " s" \
<< ref.nl << ref.nl
pass
pass
ref.ext_Info() << "End\n" << ref.nl
import os
return os.EX_OK
