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 regionProperties import regionProperties rp = regionProperties( runTime ) from fluid import createFluidMeshes fluidRegions = createFluidMeshes( rp, runTime ) from solid import createSolidMeshes solidRegions = createSolidMeshes( rp, runTime ) from fluid import createFluidFields pdf, thermof, rhof, Kf, Uf, phif, turb, DpDtf, ghf, initialMassf, pRef = createFluidFields( fluidRegions, runTime, rp ) from solid import createSolidField rhos, cps, rhosCps, Ks, Ts = createSolidField( solidRegions, runTime ) from Foam.finiteVolume.cfdTools.general.include import initContinuityErrs cumulativeContErr = initContinuityErrs() from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls( runTime ) if fluidRegions.size() : from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo( fluidRegions, runTime, rhof, phif ) from fluid import setInitialDeltaT runTime = setInitialDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum ) pass from Foam.OpenFOAM import ext_Info, nl while runTime.run(): from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls( runTime ) if fluidRegions.size() : from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo( fluidRegions, runTime, rhof, phif ) from Foam.finiteVolume.cfdTools.general.include import setDeltaT runTime = setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum ) pass runTime.increment() ext_Info()<< "Time = " << runTime.timeName() << nl << nl for i in range( fluidRegions.size() ): ext_Info() << "\nSolving for fluid region " << fluidRegions[ i ].name() << nl from fluid import readFluidMultiRegionPISOControls piso, nCorr, nNonOrthCorr, momentumPredictor, transonic, nOuterCorr = readFluidMultiRegionPISOControls( fluidRegions[ i ] ) from fluid import solveFluid cumulativeContErr = solveFluid( i, fluidRegions, pdf, thermof, rhof, Kf, Uf, phif, turb, DpDtf, ghf, initialMassf, pRef,\ nCorr, nNonOrthCorr, momentumPredictor, transonic, nOuterCorr, cumulativeContErr ) pass for i in range( solidRegions.size() ): ext_Info() << "\nSolving for solid region " << solidRegions[ i ].name() << nl from solid import readSolidMultiRegionPISOControls piso, nNonOrthCorr = readSolidMultiRegionPISOControls( solidRegions[ i ] ) from solid import solveSolid solveSolid( i, rhosCps, Ks, Ts, nNonOrthCorr ) pass runTime.write(); ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \ << " ClockTime = " << runTime.elapsedClockTime() << " s" \ << nl << nl ext_Info() << "End\n" pass 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 regionProperties import regionProperties rp = regionProperties(runTime) from fluid import createFluidMeshes fluidRegions = createFluidMeshes(rp, runTime) from solid import createSolidMeshes solidRegions = createSolidMeshes(rp, runTime) from fluid import createFluidFields pdf, thermof, rhof, Kf, Uf, phif, turb, DpDtf, ghf, initialMassf, pRef = createFluidFields( fluidRegions, runTime, rp) from solid import createSolidField rhos, cps, rhosCps, Ks, Ts = createSolidField(solidRegions, runTime) from Foam.finiteVolume.cfdTools.general.include import initContinuityErrs cumulativeContErr = initContinuityErrs() from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls(runTime) if fluidRegions.size(): from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo(fluidRegions, runTime, rhof, phif) from fluid import setInitialDeltaT runTime = setInitialDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum) pass from Foam.OpenFOAM import ext_Info, nl while runTime.run(): from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls(runTime) if fluidRegions.size(): from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo(fluidRegions, runTime, rhof, phif) from Foam.finiteVolume.cfdTools.general.include import setDeltaT runTime = setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum) pass runTime.increment() ext_Info() << "Time = " << runTime.timeName() << nl << nl for i in range(fluidRegions.size()): ext_Info( ) << "\nSolving for fluid region " << fluidRegions[i].name() << nl from fluid import readFluidMultiRegionPISOControls piso, nCorr, nNonOrthCorr, momentumPredictor, transonic, nOuterCorr = readFluidMultiRegionPISOControls( fluidRegions[i]) from fluid import solveFluid cumulativeContErr = solveFluid( i, fluidRegions, pdf, thermof, rhof, Kf, Uf, phif, turb, DpDtf, ghf, initialMassf, pRef,\ nCorr, nNonOrthCorr, momentumPredictor, transonic, nOuterCorr, cumulativeContErr ) pass for i in range(solidRegions.size()): ext_Info( ) << "\nSolving for solid region " << solidRegions[i].name() << nl from solid import readSolidMultiRegionPISOControls piso, nNonOrthCorr = readSolidMultiRegionPISOControls( solidRegions[i]) from solid import solveSolid solveSolid(i, rhosCps, Ks, Ts, nNonOrthCorr) pass runTime.write() ext_Info() << "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \ << " ClockTime = " << runTime.elapsedClockTime() << " s" \ << nl << nl ext_Info() << "End\n" pass 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 regionProperties import regionProperties rp = regionProperties( runTime ) from fluid import createFluidMeshes fluidRegions = createFluidMeshes( rp, runTime ) from solid import createSolidMeshes,createSolidField solidRegions=createSolidMeshes( rp,runTime ) from fluid import createFluidFields thermoFluid, rhoFluid, KFluid, UFluid, phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid = createFluidFields( fluidRegions, runTime ) from solid import createSolidField rhos, cps, rhosCps, Ks, Ts = createSolidField( solidRegions, runTime ) from fluid import initContinuityErrs cumulativeContErr = initContinuityErrs( fluidRegions.size() ) from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls( runTime ) from Foam.OpenFOAM import ext_Info, nl if fluidRegions.size() : from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo( fluidRegions, runTime, rhoFluid, phiFluid ) from Foam.finiteVolume.cfdTools.general.include import setInitialDeltaT runTime = setInitialDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum ) pass while runTime.run() : adjustTimeStep, maxCo, maxDeltaT = readTimeControls(runTime) nOuterCorr = readPIMPLEControls( runTime ) if fluidRegions.size() : from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo( fluidRegions, runTime, rhoFluid, phiFluid ) from Foam.finiteVolume.cfdTools.general.include import setDeltaT runTime = setDeltaT( runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum ) pass runTime.increment() ext_Info()<< "Time = " << runTime.timeName() << nl << nl if nOuterCorr != 1 : for i in range( fluidRegions.size() ): from fluid import setRegionFluidFields mesh, thermo, rho, K, U, phi, g, turb, DpDt, p, psi, h, massIni = \ setRegionFluidFields( i, fluidRegions, thermoFluid, rhoFluid, KFluid, UFluid, \ phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid ) from fluid import storeOldFluidFields storeOldFluidFields( p, rho ) pass pass # --- PIMPLE loop for oCorr in range( nOuterCorr ): for i in range( fluidRegions.size() ): ext_Info() << "\nSolving for fluid region " << fluidRegions[ i ].name() << nl from fluid import setRegionFluidFields mesh, thermo, rho, K, U, phi, g, turb, DpDt, p, psi, h, massIni = \ setRegionFluidFields( i, fluidRegions, thermoFluid, rhoFluid, KFluid, UFluid, \ phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid ) from fluid import readFluidMultiRegionPIMPLEControls pimple, nCorr, nNonOrthCorr, momentumPredictor = readFluidMultiRegionPIMPLEControls( mesh ) from fluid import solveFluid cumulativeContErr = solveFluid( i, mesh, thermo, thermoFluid, rho, K, U, phi, g, h, turb, DpDt, p, psi, \ massIni, oCorr, nCorr, nOuterCorr, nNonOrthCorr, momentumPredictor, cumulativeContErr ) pass for i in range( solidRegions.size() ): ext_Info() << "\nSolving for solid region " << solidRegions[ i ].name() << nl from solid import setRegionSolidFields mesh, rho, cp, K, T = setRegionSolidFields( i, solidRegions, rhos, cps, Ks, Ts ) from solid import readSolidMultiRegionPIMPLEControls pimple, nNonOrthCorr = readSolidMultiRegionPIMPLEControls( mesh ) from solid import solveSolid solveSolid( mesh, rho, cp, K, T, nNonOrthCorr ) pass pass pass runTime.write() ext_Info()<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \ << " ClockTime = " << runTime.elapsedClockTime() << " s" \ << nl << nl ext_Info() << "End\n" 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 regionProperties import regionProperties rp = regionProperties(runTime) from fluid import createFluidMeshes fluidRegions = createFluidMeshes(rp, runTime) from solid import createSolidMeshes, createSolidField solidRegions = createSolidMeshes(rp, runTime) from fluid import createFluidFields thermoFluid, rhoFluid, KFluid, UFluid, phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid = createFluidFields( fluidRegions, runTime) from solid import createSolidField rhos, cps, rhosCps, Ks, Ts = createSolidField(solidRegions, runTime) from fluid import initContinuityErrs cumulativeContErr = initContinuityErrs(fluidRegions.size()) from Foam.finiteVolume.cfdTools.general.include import readTimeControls adjustTimeStep, maxCo, maxDeltaT = readTimeControls(runTime) from Foam.OpenFOAM import ext_Info, nl if fluidRegions.size(): from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo(fluidRegions, runTime, rhoFluid, phiFluid) from Foam.finiteVolume.cfdTools.general.include import setInitialDeltaT runTime = setInitialDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum) pass while runTime.run(): adjustTimeStep, maxCo, maxDeltaT = readTimeControls(runTime) nOuterCorr = readPIMPLEControls(runTime) if fluidRegions.size(): from fluid import compressubibleMultiRegionCourantNo CoNum = compressubibleMultiRegionCourantNo(fluidRegions, runTime, rhoFluid, phiFluid) from Foam.finiteVolume.cfdTools.general.include import setDeltaT runTime = setDeltaT(runTime, adjustTimeStep, maxCo, maxDeltaT, CoNum) pass runTime.increment() ext_Info() << "Time = " << runTime.timeName() << nl << nl if nOuterCorr != 1: for i in range(fluidRegions.size()): from fluid import setRegionFluidFields mesh, thermo, rho, K, U, phi, g, turb, DpDt, p, psi, h, massIni = \ setRegionFluidFields( i, fluidRegions, thermoFluid, rhoFluid, KFluid, UFluid, \ phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid ) from fluid import storeOldFluidFields storeOldFluidFields(p, rho) pass pass # --- PIMPLE loop for oCorr in range(nOuterCorr): for i in range(fluidRegions.size()): ext_Info() << "\nSolving for fluid region " << fluidRegions[ i].name() << nl from fluid import setRegionFluidFields mesh, thermo, rho, K, U, phi, g, turb, DpDt, p, psi, h, massIni = \ setRegionFluidFields( i, fluidRegions, thermoFluid, rhoFluid, KFluid, UFluid, \ phiFluid, gFluid, turbulence, DpDtFluid, initialMassFluid ) from fluid import readFluidMultiRegionPIMPLEControls pimple, nCorr, nNonOrthCorr, momentumPredictor = readFluidMultiRegionPIMPLEControls( mesh) from fluid import solveFluid cumulativeContErr = solveFluid( i, mesh, thermo, thermoFluid, rho, K, U, phi, g, h, turb, DpDt, p, psi, \ massIni, oCorr, nCorr, nOuterCorr, nNonOrthCorr, momentumPredictor, cumulativeContErr ) pass for i in range(solidRegions.size()): ext_Info() << "\nSolving for solid region " << solidRegions[ i].name() << nl from solid import setRegionSolidFields mesh, rho, cp, K, T = setRegionSolidFields( i, solidRegions, rhos, cps, Ks, Ts) from solid import readSolidMultiRegionPIMPLEControls pimple, nNonOrthCorr = readSolidMultiRegionPIMPLEControls(mesh) from solid import solveSolid solveSolid(mesh, rho, cp, K, T, nNonOrthCorr) pass pass pass runTime.write() ext_Info()<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" \ << " ClockTime = " << runTime.elapsedClockTime() << " s" \ << nl << nl ext_Info() << "End\n" import os return os.EX_OK