From f37ec93404a3b92b19a4bb775aaa0d33b2a776c9 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Wed, 13 Dec 2017 08:37:25 +0000 Subject: [PATCH] ENH: chtMultiRegionFoam: Added support for reactions chtMultiRegionFoam now supports reaction/combustion modelling in fluid regions in the same way as reactingFoam. TUT: chtMultiRegionFoam: Added reverseBurner tutorial This tutorial demonstrates chtMultiRegionFoam's combustion capability --- .../chtMultiRegionFoam/Make/options | 11 +- .../chtMultiRegionFoam/chtMultiRegionFoam.C | 3 +- .../chtMultiRegionFoam/fluid/EEqn.H | 3 +- .../chtMultiRegionFoam/fluid/UEqn.H | 2 +- .../chtMultiRegionFoam/fluid/YEqn.H | 61 ++++++++ .../fluid/createFluidFields.H | 64 ++++++-- .../fluid/setRegionFluidFields.H | 29 +++- .../chtMultiRegionFoam/fluid/solveFluid.H | 3 +- .../bottomWater/thermophysicalProperties | 2 +- .../constant/topAir/thermophysicalProperties | 2 +- .../reverseBurner/0/gas/CH4 | 50 +++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/N2 | 50 +++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/O2 | 50 +++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/T | 47 ++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/U | 48 ++++++ .../reverseBurner/0/gas/Ydefault | 50 +++++++ .../reverseBurner/0/gas/alphat | 44 ++++++ .../reverseBurner/0/gas/epsilon | 45 ++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/k | 45 ++++++ .../reverseBurner/0/gas/nut | 44 ++++++ .../chtMultiRegionFoam/reverseBurner/0/gas/p | 34 +++++ .../reverseBurner/0/gas/p_rgh | 44 ++++++ .../reverseBurner/0/solid/T | 45 ++++++ .../reverseBurner/0/solid/p | 34 +++++ .../chtMultiRegionFoam/reverseBurner/Allclean | 9 ++ .../chtMultiRegionFoam/reverseBurner/Allrun | 20 +++ .../constant/gas/chemistryProperties | 40 +++++ .../constant/gas/combustionProperties | 28 ++++ .../reverseBurner/constant/gas/fvOptions | 38 +++++ .../reverseBurner/constant/gas/g | 21 +++ .../constant/gas/radiationProperties | 23 +++ .../reverseBurner/constant/gas/reactions | 40 +++++ .../constant/gas/thermo.compressibleGas | 124 ++++++++++++++++ .../constant/gas/thermophysicalProperties | 38 +++++ .../constant/gas/turbulenceProperties | 28 ++++ .../reverseBurner/constant/regionProperties | 24 +++ .../constant/solid/radiationProperties | 23 +++ .../constant/solid/thermophysicalProperties | 52 +++++++ .../reverseBurner/system/blockMeshDict | 138 ++++++++++++++++++ .../reverseBurner/system/controlDict | 54 +++++++ .../reverseBurner/system/decomposeParDict | 32 ++++ .../reverseBurner/system/fvSchemes | 42 ++++++ .../reverseBurner/system/fvSolution | 22 +++ .../reverseBurner/system/gas/decomposeParDict | 1 + .../reverseBurner/system/gas/fvSchemes | 57 ++++++++ .../reverseBurner/system/gas/fvSolution | 81 ++++++++++ .../reverseBurner/system/gas/topoSetDict | 31 ++++ .../system/solid/decomposeParDict | 1 + .../reverseBurner/system/solid/fvSchemes | 49 +++++++ .../reverseBurner/system/solid/fvSolution | 48 ++++++ .../bottomAir/thermophysicalProperties | 2 +- 51 files changed, 1854 insertions(+), 22 deletions(-) create mode 100644 applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/CH4 create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/N2 create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/O2 create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/T create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/U create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/Ydefault create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/alphat create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/epsilon create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/k create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/nut create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/T create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p create mode 100755 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allclean create mode 100755 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allrun create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/chemistryProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/combustionProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/g create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/radiationProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/reactions create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/turbulenceProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/regionProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/radiationProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/thermophysicalProperties create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/blockMeshDict create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/controlDict create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/decomposeParDict create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSchemes create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSolution create mode 120000 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/decomposeParDict create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/topoSetDict create mode 120000 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/decomposeParDict create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSchemes create mode 100644 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSolution diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/Make/options b/applications/solvers/heatTransfer/chtMultiRegionFoam/Make/options index fd53d162e6..deef74acce 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/Make/options +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/Make/options @@ -1,4 +1,3 @@ - EXE_INC = \ -I. \ -I./fluid \ @@ -12,7 +11,11 @@ EXE_INC = \ -I$(LIB_SRC)/transportModels/compressible/lnInclude \ -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \ -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \ + -I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \ -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \ + -I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \ + -I$(LIB_SRC)/ODE/lnInclude \ + -I$(LIB_SRC)/combustionModels/lnInclude \ -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \ -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \ -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \ @@ -22,8 +25,12 @@ EXE_INC = \ EXE_LIBS = \ -lcompressibleTransportModels \ -lfluidThermophysicalModels \ - -lsolidThermo \ -lspecie \ + -lreactionThermophysicalModels \ + -lsolidThermo \ + -lchemistryModel \ + -lODE \ + -lcombustionModels \ -lturbulenceModels \ -lcompressibleTurbulenceModels \ -lmeshTools \ diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C index ab34e7873e..1bababea09 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C @@ -37,8 +37,9 @@ Description \*---------------------------------------------------------------------------*/ #include "fvCFD.H" -#include "rhoThermo.H" #include "turbulentFluidThermoModel.H" +#include "rhoReactionThermo.H" +#include "CombustionModel.H" #include "fixedGradientFvPatchFields.H" #include "regionProperties.H" #include "compressibleCourantNo.H" diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H index d8f4516f5e..8f30905c40 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H @@ -15,10 +15,11 @@ ) : -dpdt ) - - fvm::laplacian(turb.alphaEff(), he) + - fvm::laplacian(turbulence.alphaEff(), he) == rho*(U&g) + rad.Sh(thermo, he) + + Qdot + fvOptions(rho, he) ); diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H index a12d0e8146..9b5ce0e396 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/UEqn.H @@ -6,7 +6,7 @@ ( fvm::ddt(rho, U) + fvm::div(phi, U) + MRF.DDt(rho, U) - + turb.divDevRhoReff(U) + + turbulence.divDevRhoReff(U) == fvOptions(rho, U) ); diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H new file mode 100644 index 0000000000..c32c1fdab2 --- /dev/null +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H @@ -0,0 +1,61 @@ +tmp> mvConvection(nullptr); + +if (Y.size()) +{ + mvConvection = tmp> + ( + fv::convectionScheme::New + ( + mesh, + fields, + phi, + mesh.divScheme("div(phi,Yi_h)") + ) + ); +} + +{ + reaction.correct(); + Qdot = reaction.Qdot(); + volScalarField Yt + ( + IOobject("Yt", runTime.timeName(), mesh), + mesh, + dimensionedScalar("Yt", dimless, 0) + ); + + forAll(Y, i) + { + if (i != inertIndex && composition.active(i)) + { + volScalarField& Yi = Y[i]; + + fvScalarMatrix YiEqn + ( + fvm::ddt(rho, Yi) + + mvConvection->fvmDiv(phi, Yi) + - fvm::laplacian(turbulence.muEff(), Yi) + == + reaction.R(Yi) + + fvOptions(rho, Yi) + ); + + YiEqn.relax(); + + fvOptions.constrain(YiEqn); + + YiEqn.solve(mesh.solver("Yi")); + + fvOptions.correct(Yi); + + Yi.max(0.0); + Yt += Yi; + } + } + + if (Y.size()) + { + Y[inertIndex] = scalar(1) - Yt; + Y[inertIndex].max(0.0); + } +} diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H index 5bcd4b4999..be1e3e57c0 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H @@ -1,5 +1,5 @@ // Initialise fluid field pointer lists -PtrList thermoFluid(fluidRegions.size()); +PtrList thermoFluid(fluidRegions.size()); PtrList rhoFluid(fluidRegions.size()); PtrList UFluid(fluidRegions.size()); PtrList phiFluid(fluidRegions.size()); @@ -7,11 +7,15 @@ PtrList gFluid(fluidRegions.size()); PtrList hRefFluid(fluidRegions.size()); PtrList ghFluid(fluidRegions.size()); PtrList ghfFluid(fluidRegions.size()); -PtrList turbulence(fluidRegions.size()); +PtrList turbulenceFluid(fluidRegions.size()); +PtrList> reactionFluid(fluidRegions.size()); PtrList p_rghFluid(fluidRegions.size()); PtrList radiation(fluidRegions.size()); PtrList KFluid(fluidRegions.size()); PtrList dpdtFluid(fluidRegions.size()); +PtrList::fieldTable> + fieldsFluid(fluidRegions.size()); +PtrList QdotFluid(fluidRegions.size()); List initialMassFluid(fluidRegions.size()); List frozenFlowFluid(fluidRegions.size(), false); @@ -29,12 +33,7 @@ forAll(fluidRegions, i) << fluidRegions[i].name() << nl << endl; Info<< " Adding to thermoFluid\n" << endl; - - thermoFluid.set - ( - i, - rhoThermo::New(fluidRegions[i]).ptr() - ); + thermoFluid.set(i, rhoReactionThermo::New(fluidRegions[i]).ptr()); Info<< " Adding to rhoFluid\n" << endl; rhoFluid.set @@ -156,8 +155,8 @@ forAll(fluidRegions, i) ) ); - Info<< " Adding to turbulence\n" << endl; - turbulence.set + Info<< " Adding to turbulenceFluid\n" << endl; + turbulenceFluid.set ( i, compressible::turbulenceModel::New @@ -169,6 +168,17 @@ forAll(fluidRegions, i) ).ptr() ); + Info<< " Adding to reactionFluid\n" << endl; + reactionFluid.set + ( + i, + CombustionModel::New + ( + thermoFluid[i], + turbulenceFluid[i] + ) + ); + p_rghFluid.set ( i, @@ -191,6 +201,7 @@ forAll(fluidRegions, i) fluidRegions[i].setFluxRequired(p_rghFluid[i].name()); + Info<< " Adding to radiationFluid\n" << endl; radiation.set ( i, @@ -227,6 +238,37 @@ forAll(fluidRegions, i) ) ); + Info<< " Adding to fieldsFluid\n" << endl; + fieldsFluid.set + ( + i, + new multivariateSurfaceInterpolationScheme::fieldTable + ); + forAll(thermoFluid[i].composition().Y(), j) + { + fieldsFluid[i].add(thermoFluid[i].composition().Y()[j]); + } + fieldsFluid[i].add(thermoFluid[i].he()); + + Info<< " Adding to QdotFluid\n" << endl; + QdotFluid.set + ( + i, + new volScalarField + ( + IOobject + ( + "Qdot", + runTime.timeName(), + fluidRegions[i], + IOobject::READ_IF_PRESENT, + IOobject::AUTO_WRITE + ), + fluidRegions[i], + dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0) + ) + ); + const dictionary& pimpleDict = fluidRegions[i].solutionDict().subDict("PIMPLE"); pimpleDict.readIfPresent("frozenFlow", frozenFlowFluid[i]); @@ -245,7 +287,7 @@ forAll(fluidRegions, i) new fv::options(fluidRegions[i]) ); - turbulence[i].validate(); + turbulenceFluid[i].validate(); pRefCellFluid[i] = 0; pRefValueFluid[i] = 0.0; diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H index e725968af3..6c67f6b9c0 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H @@ -1,13 +1,33 @@ fvMesh& mesh = fluidRegions[i]; - rhoThermo& thermo = thermoFluid[i]; + CombustionModel& reaction = reactionFluid[i]; + + rhoReactionThermo& thermo = reaction.thermo(); thermo.validate(args.executable(), "h", "e"); + basicMultiComponentMixture& composition = thermo.composition(); + PtrList& Y = composition.Y(); + + label inertIndex = -1; + if (Y.size()) + { + const word inertSpecie(thermo.lookup("inertSpecie")); + if (!composition.species().found(inertSpecie)) + { + FatalIOErrorIn(args.executable().c_str(), thermo) + << "Inert specie " << inertSpecie + << " not found in available species " + << composition.species() + << exit(FatalIOError); + } + inertIndex = composition.species()[inertSpecie]; + } + volScalarField& rho = rhoFluid[i]; volVectorField& U = UFluid[i]; surfaceScalarField& phi = phiFluid[i]; - compressible::turbulenceModel& turb = turbulence[i]; + compressible::turbulenceModel& turbulence = turbulenceFluid[i]; volScalarField& K = KFluid[i]; volScalarField& dpdt = dpdtFluid[i]; @@ -20,6 +40,11 @@ const volScalarField& gh = ghFluid[i]; const surfaceScalarField& ghf = ghfFluid[i]; + multivariateSurfaceInterpolationScheme::fieldTable& fields = + fieldsFluid[i]; + + volScalarField& Qdot = QdotFluid[i]; + radiation::radiationModel& rad = radiation[i]; IOMRFZoneList& MRF = MRFfluid[i]; diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H index 6f511e0292..d4d18e9fdd 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H @@ -15,6 +15,7 @@ else } #include "UEqn.H" + #include "YEqn.H" #include "EEqn.H" // --- PISO loop @@ -23,7 +24,7 @@ else #include "pEqn.H" } - turb.correct(); + turbulence.correct(); rho = thermo.rho(); } diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/bottomWater/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/bottomWater/thermophysicalProperties index c1602d4d09..03a2776252 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/bottomWater/thermophysicalProperties +++ b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/bottomWater/thermophysicalProperties @@ -17,7 +17,7 @@ FoamFile thermoType { type heRhoThermo; - mixture pureMixture; + mixture singleComponentMixture; transport const; thermo hConst; equationOfState rhoConst; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/topAir/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/topAir/thermophysicalProperties index 3f94dc8ad7..bbd0bb688e 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/topAir/thermophysicalProperties +++ b/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/constant/topAir/thermophysicalProperties @@ -17,7 +17,7 @@ FoamFile thermoType { type heRhoThermo; - mixture pureMixture; + mixture singleComponentMixture; transport const; thermo hConst; equationOfState perfectGas; diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/CH4 b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/CH4 new file mode 100644 index 0000000000..6395bf6ebf --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/CH4 @@ -0,0 +1,50 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object CH4; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + inletFuel + { + type fixedValue; + value uniform 1; + } + inletAir + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + + } + gas_to_solid + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/N2 b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/N2 new file mode 100644 index 0000000000..7d570a6484 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/N2 @@ -0,0 +1,50 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object N2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.79; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + inletFuel + { + type fixedValue; + value uniform 0; + } + inletAir + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + + } + gas_to_solid + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/O2 b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/O2 new file mode 100644 index 0000000000..a1316b69ca --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/O2 @@ -0,0 +1,50 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object O2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.21; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + inletFuel + { + type fixedValue; + value uniform 0; + } + inletAir + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + + } + gas_to_solid + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/T b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/T new file mode 100644 index 0000000000..0bb20553ae --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/T @@ -0,0 +1,47 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object T; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 0 0 0 1 0 0 0 ]; + +internalField uniform 300; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + value $internalField; + inletValue $internalField; + } + gas_to_solid + { + type compressible::turbulentTemperatureCoupledBaffleMixed; + value $internalField; + Tnbr T; + kappaMethod fluidThermo; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/U b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/U new file mode 100644 index 0000000000..d2deb67078 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/U @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + location "0/gas"; + object U; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 0 1 -1 0 0 0 0 ]; + +internalField uniform (0 0 0); + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + inletFuel + { + type fixedValue; + value uniform (0 0.5 0); + } + inletAir + { + type fixedValue; + value uniform (0 1.5 0); + } + outlet + { + type pressureInletOutletVelocity; + value $internalField; + } + gas_to_solid + { + type noSlip; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/Ydefault b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/Ydefault new file mode 100644 index 0000000000..9d658a3b59 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/Ydefault @@ -0,0 +1,50 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object Ydefault; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + inletFuel + { + type fixedValue; + value $internalField; + } + inletAir + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + + } + gas_to_solid + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/alphat b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/alphat new file mode 100644 index 0000000000..b022543181 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/alphat @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volScalarField; + location "0/gas"; + object alphat; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type calculated; + value $internalField; + } + outlet + { + type calculated; + value $internalField; + } + gas_to_solid + { + type compressible::alphatWallFunction; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/epsilon b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/epsilon new file mode 100644 index 0000000000..ca483a1e1b --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/epsilon @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object epsilon; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 0 2 -3 0 0 0 0 ]; + +internalField uniform 0.01; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + value $internalField; + inletValue $internalField; + } + gas_to_solid + { + type epsilonWallFunction; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/k b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/k new file mode 100644 index 0000000000..e12d41f0a4 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/k @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object k; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 0 2 -2 0 0 0 0 ]; + +internalField uniform 0.1; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type fixedValue; + value $internalField; + } + outlet + { + type inletOutlet; + value $internalField; + inletValue $internalField; + } + gas_to_solid + { + type kqRWallFunction; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/nut b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/nut new file mode 100644 index 0000000000..a726e214c7 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/nut @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volScalarField; + location "0/gas"; + object nutt; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type calculated; + value $internalField; + } + outlet + { + type calculated; + value $internalField; + } + gas_to_solid + { + type nutkWallFunction; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p new file mode 100644 index 0000000000..1fa18d55c1 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p @@ -0,0 +1,34 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 1 -1 -2 0 0 0 0 ]; + +internalField uniform 1e5; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + ".*" + { + type calculated; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh new file mode 100644 index 0000000000..cf2d9d0fba --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/gas"; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 1 -1 -2 0 0 0 0 ]; + +internalField uniform 1e5; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + "inlet.*" + { + type fixedFluxPressure; + value $internalField; + } + outlet + { + type fixedValue; + value $internalField; + } + gas_to_solid + { + type fixedFluxPressure; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/T b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/T new file mode 100644 index 0000000000..1f6bbe5e61 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/T @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/solid"; + object T; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [ 0 0 0 1 0 0 0 ]; + +internalField uniform 300; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + external + { + type externalWallHeatFluxTemperature; + kappaMethod solidThermo; + mode coefficient; + h uniform 1e3; + Ta $internalField; + value $internalField; + } + solid_to_gas + { + type compressible::turbulentTemperatureCoupledBaffleMixed; + value $internalField; + Tnbr T; + kappaMethod solidThermo; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p new file mode 100644 index 0000000000..4968950d50 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p @@ -0,0 +1,34 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0/solid"; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ + #includeEtc "caseDicts/setConstraintTypes" + + ".*" + { + type calculated; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allclean b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allclean new file mode 100755 index 0000000000..b95f77b46c --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allclean @@ -0,0 +1,9 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase + +rm -rf constant/*/polyMesh 0/cellToRegion 0/*/cellToRegion diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allrun b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allrun new file mode 100755 index 0000000000..054d3cff03 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allrun @@ -0,0 +1,20 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +runApplication blockMesh + +runApplication splitMeshRegions -cellZonesOnly -overwrite + +runApplication topoSet -region gas + +paraFoam -region gas -touch +paraFoam -region solid -touch + +runApplication decomposePar -allRegions + +runParallel $(getApplication) + +runApplication reconstructPar -allRegions -newTimes diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/chemistryProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/chemistryProperties new file mode 100644 index 0000000000..809a143c23 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/chemistryProperties @@ -0,0 +1,40 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object chemistryProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +chemistryType +{ + solver EulerImplicit; +} + +chemistry on; + +initialChemicalTimeStep 1e-07; + +odeCoeffs +{ + solver seulex; + eps 0.05; +} + +EulerImplicitCoeffs +{ + cTauChem 1; + equilibriumRateLimiter off; +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/combustionProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/combustionProperties new file mode 100644 index 0000000000..36dc402e0f --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/combustionProperties @@ -0,0 +1,28 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object combustionProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +combustionModel PaSR; + +active yes; + +PaSRCoeffs +{ + Cmix 1.0; +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions new file mode 100644 index 0000000000..64385fa9c6 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions @@ -0,0 +1,38 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object fvOptions; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +options +{ + energySource + { + type scalarSemiImplicitSource; + + timeStart 0.2; + duration 0.3; + selectionMode cellSet; + cellSet ignition; + volumeMode specific; + + injectionRateSuSp + { + h (2e7 0); // kg/m/s^3 + } + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/g b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/g new file mode 100644 index 0000000000..3708f63c33 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/g @@ -0,0 +1,21 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; + +value (0 -9.81 0); + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/radiationProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/radiationProperties new file mode 100644 index 0000000000..94cc3749b1 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/radiationProperties @@ -0,0 +1,23 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object radiationProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +radiation off; + +radiationModel none; + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/reactions b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/reactions new file mode 100644 index 0000000000..562407030f --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/reactions @@ -0,0 +1,40 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant/gas"; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +species +( + O2 + H2O + CH4 + CO2 + N2 +); + +reactions +{ + methaneReaction + { + type irreversibleArrheniusReaction; + reaction "CH4 + 2O2 = CO2 + 2H2O"; + A 5.2e16; + beta 0; + Ta 14906; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas new file mode 100644 index 0000000000..1bacad52c3 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas @@ -0,0 +1,124 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermo.compressibleGas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +O2 +{ + specie + { + molWeight 31.9988; + } + thermodynamics + { + Tlow 200; + Thigh 5000; + Tcommon 1000; + highCpCoeffs ( 3.69758 0.00061352 -1.25884e-07 1.77528e-11 -1.13644e-15 -1233.93 3.18917 ); + lowCpCoeffs ( 3.21294 0.00112749 -5.75615e-07 1.31388e-09 -8.76855e-13 -1005.25 6.03474 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + +H2O +{ + specie + { + molWeight 18.0153; + } + thermodynamics + { + Tlow 200; + Thigh 5000; + Tcommon 1000; + highCpCoeffs ( 2.67215 0.00305629 -8.73026e-07 1.201e-10 -6.39162e-15 -29899.2 6.86282 ); + lowCpCoeffs ( 3.38684 0.00347498 -6.3547e-06 6.96858e-09 -2.50659e-12 -30208.1 2.59023 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + +CH4 +{ + specie + { + molWeight 16.0428; + } + thermodynamics + { + Tlow 200; + Thigh 6000; + Tcommon 1000; + highCpCoeffs ( 1.63543 0.0100844 -3.36924e-06 5.34973e-10 -3.15528e-14 -10005.6 9.9937 ); + lowCpCoeffs ( 5.14988 -0.013671 4.91801e-05 -4.84744e-08 1.66694e-11 -10246.6 -4.64132 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + +CO2 +{ + specie + { + molWeight 44.01; + } + thermodynamics + { + Tlow 200; + Thigh 5000; + Tcommon 1000; + highCpCoeffs ( 4.45362 0.00314017 -1.27841e-06 2.394e-10 -1.66903e-14 -48967 -0.955396 ); + lowCpCoeffs ( 2.27572 0.00992207 -1.04091e-05 6.86669e-09 -2.11728e-12 -48373.1 10.1885 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + +N2 +{ + specie + { + molWeight 28.0134; + } + thermodynamics + { + Tlow 200; + Thigh 5000; + Tcommon 1000; + highCpCoeffs ( 2.92664 0.00148798 -5.68476e-07 1.0097e-10 -6.75335e-15 -922.798 5.98053 ); + lowCpCoeffs ( 3.29868 0.00140824 -3.96322e-06 5.64152e-09 -2.44486e-12 -1020.9 3.95037 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties new file mode 100644 index 0000000000..b71964fd54 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties @@ -0,0 +1,38 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant/gas"; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture reactingMixture; + transport sutherland; + thermo janaf; + equationOfState perfectGas; + specie specie; + energy sensibleEnthalpy; +} + +inertSpecie N2; + +chemistryReader foamChemistryReader; + +foamChemistryFile "$FOAM_CASE/constant/gas/reactions"; + +foamChemistryThermoFile "$FOAM_CASE/constant/gas/thermo.compressibleGas"; + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/turbulenceProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/turbulenceProperties new file mode 100644 index 0000000000..0e69ddc759 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/turbulenceProperties @@ -0,0 +1,28 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType RAS; + +RAS +{ + RASModel kEpsilon; + + turbulence on; + + printCoeffs on; +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/regionProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/regionProperties new file mode 100644 index 0000000000..b31bac6ee5 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/regionProperties @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object regionProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +regions +( + fluid (gas) + solid (solid) +); + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/radiationProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/radiationProperties new file mode 100644 index 0000000000..94cc3749b1 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/radiationProperties @@ -0,0 +1,23 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object radiationProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +radiation off; + +radiationModel none; + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/thermophysicalProperties new file mode 100644 index 0000000000..97bdf46290 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/thermophysicalProperties @@ -0,0 +1,52 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heSolidThermo; + mixture pureMixture; + transport constIso; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleEnthalpy; +} + +mixture +{ + specie + { + molWeight 50; + } + + transport + { + kappa 80; + } + + thermodynamics + { + Hf 0; + Cp 450; + } + + equationOfState + { + rho 8000; + } +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/blockMeshDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/blockMeshDict new file mode 100644 index 0000000000..4b37b4b926 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/blockMeshDict @@ -0,0 +1,138 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 0.01; + +vertices +( + (0 0 -1) (1 0 -1) (2 0 -1) (10 0 -1) (11 0 -1) (21 0 -1) (22 0 -1) + (0 10 -1) (1 10 -1) (2 10 -1) (10 10 -1) (11 10 -1) (21 10 -1) (22 10 -1) + (0 40 -1) (1 40 -1) (2 40 -1) (10 40 -1) (11 40 -1) (21 40 -1) (22 40 -1) + (0 50 -1) (1 50 -1) (2 50 -1) (10 50 -1) (11 50 -1) (21 50 -1) (22 50 -1) + (0 51 -1) (1 51 -1) (2 51 -1) (10 51 -1) (11 51 -1) (21 51 -1) (22 51 -1) + + (0 0 1) (1 0 1) (2 0 1) (10 0 1) (11 0 1) (21 0 1) (22 0 1) + (0 10 1) (1 10 1) (2 10 1) (10 10 1) (11 10 1) (21 10 1) (22 10 1) + (0 40 1) (1 40 1) (2 40 1) (10 40 1) (11 40 1) (21 40 1) (22 40 1) + (0 50 1) (1 50 1) (2 50 1) (10 50 1) (11 50 1) (21 50 1) (22 50 1) + (0 51 1) (1 51 1) (2 51 1) (10 51 1) (11 51 1) (21 51 1) (22 51 1) +); + +blocks +( + hex ( 0 1 8 7 35 36 43 42) gas ( 4 40 1) simpleGrading (1 1 1) + hex ( 1 2 9 8 36 37 44 43) solid ( 4 40 1) simpleGrading (1 1 1) + hex ( 2 3 10 9 37 38 45 44) gas (32 40 1) simpleGrading (1 1 1) + hex ( 3 4 11 10 38 39 46 45) solid ( 4 40 1) simpleGrading (1 1 1) + hex ( 4 5 12 11 39 40 47 46) gas (40 40 1) simpleGrading (1 1 1) + hex ( 5 6 13 12 40 41 48 47) solid ( 4 40 1) simpleGrading (1 1 1) + + hex ( 7 8 15 14 42 43 50 49) gas ( 4 120 1) simpleGrading (1 1 1) + hex ( 8 9 16 15 43 44 51 50) gas ( 4 120 1) simpleGrading (1 1 1) + hex ( 9 10 17 16 44 45 52 51) gas (32 120 1) simpleGrading (1 1 1) + hex (10 11 18 17 45 46 53 52) solid ( 4 120 1) simpleGrading (1 1 1) + hex (11 12 19 18 46 47 54 53) gas (40 120 1) simpleGrading (1 1 1) + hex (12 13 20 19 47 48 55 54) solid ( 4 120 1) simpleGrading (1 1 1) + + hex (14 15 22 21 49 50 57 56) gas ( 4 40 1) simpleGrading (1 1 1) + hex (15 16 23 22 50 51 58 57) gas ( 4 40 1) simpleGrading (1 1 1) + hex (16 17 24 23 51 52 59 58) gas (32 40 1) simpleGrading (1 1 1) + hex (17 18 25 24 52 53 60 59) gas ( 4 40 1) simpleGrading (1 1 1) + hex (18 19 26 25 53 54 61 60) gas (40 40 1) simpleGrading (1 1 1) + hex (19 20 27 26 54 55 62 61) solid ( 4 40 1) simpleGrading (1 1 1) + + hex (21 22 29 28 56 57 64 63) solid ( 4 4 1) simpleGrading (1 1 1) + hex (22 23 30 29 57 58 65 64) solid ( 4 4 1) simpleGrading (1 1 1) + hex (23 24 31 30 58 59 66 65) solid (32 4 1) simpleGrading (1 1 1) + hex (24 25 32 31 59 60 67 66) solid ( 4 4 1) simpleGrading (1 1 1) + hex (25 26 33 32 60 61 68 67) solid (40 4 1) simpleGrading (1 1 1) + hex (26 27 34 33 61 62 69 68) solid ( 4 4 1) simpleGrading (1 1 1) +); + +edges +( +); + +defaultPatch +{ + type empty; + name frontAndBack; +} + +boundary +( + inletFuel + { + type patch; + faces + ( + (0 1 36 35) + ); + } + inletAir + { + type patch; + faces + ( + (2 3 38 37) + ); + } + outlet + { + type patch; + faces + ( + (4 5 40 39) + ); + } + axis + { + type symmetryPlane; + faces + ( + (0 7 42 35) + (7 14 49 42) + (14 21 56 49) + (21 28 63 56) + ); + } + external + { + type patch; + faces + ( + (1 2 37 36) + (3 4 39 38) + (5 6 41 40) + (6 13 48 41) + (13 20 55 48) + (20 27 62 55) + (27 34 69 62) + (28 29 64 63) + (29 30 65 64) + (30 31 66 65) + (31 32 67 66) + (32 33 68 67) + (33 34 69 68) + ); + } +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/controlDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/controlDict new file mode 100644 index 0000000000..e0ce92c8a7 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/controlDict @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application chtMultiRegionFoam; + +startFrom latestTime; + +startTime 0; + +stopAt endTime; + +endTime 1; + +deltaT 1e-5; + +writeControl adjustableRunTime; + +writeInterval 1e-2; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 8; + +writeCompression off; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +maxCo 0.2; + +maxDi 5.0; + +adjustTimeStep yes; + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/decomposeParDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/decomposeParDict new file mode 100644 index 0000000000..2429fa4fad --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/decomposeParDict @@ -0,0 +1,32 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + note "mesh decomposition control dictionary"; + location "system"; + object decomposeParDict; +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 4; + +method hierarchical; + +hierarchicalCoeffs +{ + n (1 4 1); + delta 0.001; + order xyz; +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSchemes b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSchemes new file mode 100644 index 0000000000..a61bdb780e --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSchemes @@ -0,0 +1,42 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ +} + +gradSchemes +{ +} + +divSchemes +{ +} + +laplacianSchemes +{ +} + +interpolationSchemes +{ +} + +snGradSchemes +{ +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSolution new file mode 100644 index 0000000000..cde9e49189 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/fvSolution @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +PIMPLE +{ + nOuterCorrectors 1; +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/decomposeParDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/decomposeParDict new file mode 120000 index 0000000000..46638e60f6 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/decomposeParDict @@ -0,0 +1 @@ +../decomposeParDict \ No newline at end of file diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes new file mode 100644 index 0000000000..74b4dd4a89 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes @@ -0,0 +1,57 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + div(phi,U) Gauss upwind; + div(phi,K) Gauss linear; + div(phi,h) Gauss upwind; + div(phi,k) Gauss upwind; + div(phi,epsilon) Gauss upwind; + div(phi,R) Gauss upwind; + div(R) Gauss linear; + div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; + div(phi,Yi_h) Gauss upwind; +} + +laplacianSchemes +{ + default Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default corrected; +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution new file mode 100644 index 0000000000..1befa3157a --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution @@ -0,0 +1,81 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + rho + { + solver PCG; + preconditioner DIC; + tolerance 1e-7; + relTol 0.1; + } + + rhoFinal + { + $rho; + tolerance 1e-7; + relTol 0; + } + + p_rgh + { + solver GAMG; + tolerance 1e-7; + relTol 0.01; + smoother GaussSeidel; + } + + p_rghFinal + { + $p_rgh; + tolerance 1e-7; + relTol 0; + } + + "(U|h|k|epsilon|R|Yi)" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-7; + relTol 0.1; + } + + "(U|h|k|epsilon|R|Yi)Final" + { + $U; + tolerance 1e-7; + relTol 0; + } +} + +PIMPLE +{ + momentumPredictor yes; + nCorrectors 2; + nNonOrthogonalCorrectors 0; +} + +relaxationFactors +{ + equations + { + "h.*" 1; + "U.*" 1; + } +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/topoSetDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/topoSetDict new file mode 100644 index 0000000000..5d236e5202 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/topoSetDict @@ -0,0 +1,31 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object topoSetDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +actions +( + { + name ignition; + type cellSet; + action new; + source boxToCell; + sourceInfo + { + box (0.01 0.1 -0.01) (0.02 0.11 0.01); + } + } +); + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/decomposeParDict b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/decomposeParDict new file mode 120000 index 0000000000..46638e60f6 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/decomposeParDict @@ -0,0 +1 @@ +../decomposeParDict \ No newline at end of file diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSchemes b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSchemes new file mode 100644 index 0000000000..062f0ff436 --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSchemes @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default steadyState; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; +} + +laplacianSchemes +{ + default none; + laplacian(alpha,h) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default corrected; +} + + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSolution b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSolution new file mode 100644 index 0000000000..28faba64aa --- /dev/null +++ b/tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/solid/fvSolution @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + h + { + solver PCG; + preconditioner DIC; + tolerance 1e-06; + relTol 0.1; + } + + hFinal + { + $h; + tolerance 1e-06; + relTol 0; + } +} + +PIMPLE +{ + nNonOrthogonalCorrectors 0; +} + +relaxationFactors +{ + equations + { + h 0.7; + } +} + +// ************************************************************************* // diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/bottomAir/thermophysicalProperties b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/bottomAir/thermophysicalProperties index 3f94dc8ad7..bbd0bb688e 100644 --- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/bottomAir/thermophysicalProperties +++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/bottomAir/thermophysicalProperties @@ -17,7 +17,7 @@ FoamFile thermoType { type heRhoThermo; - mixture pureMixture; + mixture singleComponentMixture; transport const; thermo hConst; equationOfState perfectGas;