diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/H2O b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/H2O new file mode 100644 index 0000000000..4b2f7e370e --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/H2O @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object H2O; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedValue; + value $internalField; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type zeroGradient; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/T b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/T new file mode 100644 index 0000000000..3aabadd19a --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/T @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object T; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedValue; + value $internalField; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type zeroGradient; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/U b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/U new file mode 100644 index 0000000000..c41e3689f0 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/U @@ -0,0 +1,55 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + object U; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + inlet1 + { + type fixedValue; + value uniform (0.2 0 0); + } + + inlet2 + { + type fixedValue; + value uniform (0 -0.1 0); + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type fixedValue; + value uniform (0 0 0); + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/air b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/air new file mode 100644 index 0000000000..a57f142526 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/air @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedValue; + value $internalField; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type zeroGradient; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/alphat b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/alphat new file mode 100644 index 0000000000..dca4298d3b --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/alphat @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object alphat; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + "(inlet1|inlet2)" + { + type calculated; + value $internalField; + } + + outlet + { + type calculated; + value $internalField; + } + + walls + { + type compressible::alphatWallFunction; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/epsilon b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/epsilon new file mode 100644 index 0000000000..20146f0f38 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/epsilon @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object epsilon; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -3 0 0 0 0]; + +internalField uniform 3e-3; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedValue; + value $internalField; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type epsilonWallFunction; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/k b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/k new file mode 100644 index 0000000000..7186f8cdce --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/k @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object k; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0.015; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedValue; + value $internalField; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + } + + walls + { + type kqRWallFunction; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/nut b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/nut new file mode 100644 index 0000000000..5966631b88 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/nut @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object alphat; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + "(inlet1|inlet2)" + { + type calculated; + value $internalField; + } + + outlet + { + type calculated; + value $internalField; + } + + walls + { + type nutkWallFunction; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p new file mode 100644 index 0000000000..26ccba8a6b --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 100000; + +boundaryField +{ + "(inlet1|inlet2)" + { + type calculated; + value $internalField; + } + + outlet + { + type calculated; + value $internalField; + } + + walls + { + type calculated; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p_rgh b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p_rgh new file mode 100644 index 0000000000..61bfd2b24c --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/0.orig/p_rgh @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 100000; + +boundaryField +{ + "(inlet1|inlet2)" + { + type fixedFluxPressure; + value $internalField; + } + + outlet + { + type fixedValue; + value $internalField; + } + + walls + { + type fixedFluxPressure; + value $internalField; + } + + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allclean b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allclean new file mode 100755 index 0000000000..fb1f384730 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allclean @@ -0,0 +1,8 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/CleanFunctions # Tutorial clean functions +#------------------------------------------------------------------------------ + +cleanCase0 + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun new file mode 100755 index 0000000000..5f30e6c3bb --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun @@ -0,0 +1,12 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +restore0Dir + +runApplication blockMesh + +runApplication $(getApplication) + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun-parallel b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun-parallel new file mode 100755 index 0000000000..9726aa17d4 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/Allrun-parallel @@ -0,0 +1,16 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +restore0Dir + +runApplication blockMesh + +runApplication decomposePar + +runParallel $(getApplication) + +runApplication reconstructPar + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/chemistryProperties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/chemistryProperties new file mode 100644 index 0000000000..704ca0b751 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/chemistryProperties @@ -0,0 +1,25 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object chemistryProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +chemistryType +{ + solver noChemistrySolver; +} + +chemistry off; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/combustionProperties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/combustionProperties new file mode 100644 index 0000000000..6ebdf3d7fa --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/combustionProperties @@ -0,0 +1,20 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object combustionProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +combustionModel none; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/g b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/g new file mode 100644 index 0000000000..2cc65c945e --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; + +value (0 0 0); + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/radiationProperties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/radiationProperties new file mode 100644 index 0000000000..c83a13ecf5 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/radiationProperties @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object radiationProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +radiation off; + +radiationModel none; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactingCloud1Properties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactingCloud1Properties new file mode 100644 index 0000000000..4d3bd24396 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactingCloud1Properties @@ -0,0 +1,176 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object reactingCloud1Properties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active true; + coupled true; + transient yes; + cellValueSourceCorrection on; + maxCo 0.3; + + sourceTerms + { + schemes + { + rho explicit 1; + U explicit 1; + Yi explicit 1; + h explicit 1; + radiation explicit 1; + } + } + + interpolationSchemes + { + rho cell; + U cellPoint; + thermo:mu cell; + T cell; + Cp cell; + kappa cell; + p cell; + } + + integrationSchemes + { + U Euler; + T analytical; + } +} + +constantProperties +{ + rho0 1000; + T0 300; + Cp0 4100; + + constantVolume false; +} + +subModels +{ + particleForces + { + sphereDrag; + } + + injectionModels + { + patchInjector + { + type patchInjection; + parcelBasisType fixed; + patch inlet1; + U0 (1 0.0 0.0); + parcelsPerSecond 10; + nParticle 1; + sizeDistribution + { + type uniform; + uniformDistribution + { + minValue 0.0001; + maxValue 0.001; + } + } + flowRateProfile constant 1; + massTotal 1; + SOI 0.0; + duration 100.0; + } + } + + dispersionModel none; + + patchInteractionModel multiInteraction; + + heatTransferModel RanzMarshall; + + compositionModel singleMixtureFraction; + + phaseChangeModel liquidEvaporation; + + devolatilisationModel none; + + surfaceReactionModel none; + + stochasticCollisionModel none; + + surfaceFilmModel none; + + radiation off; + + multiInteractionCoeffs + { + oneInteractionOnly no; + + model2 + { + patchInteractionModel recycleInteraction; + recycleInteractionCoeffs + { + recyclePatches ((outlet inlet2)); + recycleFraction 0.8; + } + } + model1 + { + patchInteractionModel standardWallInteraction; + standardWallInteractionCoeffs + { + type rebound; + } + + writeToFile yes; + } + } + + RanzMarshallCoeffs + { + BirdCorrection true; + } + + singleMixtureFractionCoeffs + { + phases + ( + gas + { + } + liquid + { + H2O 1; + } + solid + { + } + ); + YGasTot0 0; + YLiquidTot0 1; + YSolidTot0 0; + } + + liquidEvaporationCoeffs + { + enthalpyTransfer enthalpyDifference; + + activeLiquids ( H2O ); + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactions b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactions new file mode 100644 index 0000000000..228f5f836b --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/reactions @@ -0,0 +1,8 @@ +species +( + air + H2O +); + +reactions +{} diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermo.incompressiblePoly b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermo.incompressiblePoly new file mode 100644 index 0000000000..e6c93a7113 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermo.incompressiblePoly @@ -0,0 +1,163 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object thermo.incompressiblePoly; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +/* +in order to use the Boussinesq approximation for the density is sufficient to use the icoPolynominial model by setting its coefficients in the following way +[>> BOUSSINESQ TO icoPolynomial converter] + +Rearranging Boussinesq equation in the form of rho = A + B * T +rho = rho_0 - rho_0 * Beta * (T - T_0) +rho = rho_0 + rho_0 * Beta * T_0 - rho_0 * Beta * T + +By using the following values : + rho_0 = 1; + T_0 = 298; + Beta = 0.0034; // thermal expansion coefficient of air : At normal standard conditions of 25 degree Celsius or 298 Kelvin, Thermal expansion coefficient of air is said to be around 0.0034/K. + +A = rho_0 + rho_0 * Beta * T_0 = 1+1*0.0034*298 = 2.0132 +B = -rho_0 * Beta = -0.0034 + +Therefore in the icoPolynomial subdictionaries in constant/thermo.incompressiblePoly we will set : + + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); // rho = A + B * T + } +*/ + +N2 +{ + specie + { + molWeight 28.0134; + } + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); + } + thermodynamics + { + Hf 0; + Sf 0; + CpCoeffs<8> ( 979.08 0.41787 -0.0011761 1.6742e-06 -7.2559e-10 0 0 0 ); + } + transport + { + muCoeffs<8> ( 1.5068e-06 6.1598e-08 -1.8188e-11 0 0 0 0 0 ); + kappaCoeffs<8> ( 0.0031494 8.4997e-05 -1.2621e-08 0 0 0 0 0 ); + } +} + +air +{ + $N2; +} + +O2 +{ + specie + { + molWeight 31.9988; + } + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); + } + thermodynamics + { + Hf 0; + Sf 0; + CpCoeffs<8> ( 834.84 0.29297 -0.00014959 3.4143e-07 -2.2786e-10 0 0 0 ); + } + transport + { + muCoeffs<8> ( 1.5068e-06 6.1598e-08 -1.8188e-11 0 0 0 0 0 ); + kappaCoeffs<8> ( 0.00016082 8.5301e-05 -1.4998e-08 0 0 0 0 0 ); + } +} + +H2O +{ + specie + { + molWeight 18.0153; + } + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); + } + thermodynamics + { + Hf -13423000; + Sf 10482; + CpCoeffs<8> ( 1563.1 1.604 -0.0029334 3.2168e-06 -1.1571e-09 0 0 0 ); + } + transport + { + muCoeffs<8> ( 1.5068e-06 6.1598e-08 -1.8188e-11 0 0 0 0 0 ); + kappaCoeffs<8> ( 0.0037972 0.00015336 -1.1859e-08 0 0 0 0 0 ); + } +} + +CO2 +{ + specie + { + molWeight 44.01; + } + // to be updated: following CO2 coefficients taken from 02 subdictionary + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); + } + thermodynamics + { + Hf 0; + Sf 0; + CpCoeffs<8> ( 834.84 0.29297 -0.00014959 3.4143e-07 -2.2786e-10 0 0 0 ); + } + transport + { + muCoeffs<8> ( 1.5068e-06 6.1598e-08 -1.8188e-11 0 0 0 0 0 ); + kappaCoeffs<8> ( 0.00016082 8.5301e-05 -1.4998e-08 0 0 0 0 0 ); + } +} + +air +{ + specie + { + molWeight 28.85; + } + equationOfState + { + rhoCoeffs<8> ( 2.0132 -0.0034 0 0 0 0 0 0 ); + } + thermodynamics + { + Hf 0; + Sf 0; + CpCoeffs<8> ( 948.76 0.39171 -0.00095999 1.393e-06 -6.2029e-10 0 0 0 ); + } + transport + { + muCoeffs<8> ( 1.5061e-06 6.16e-08 -1.819e-11 0 0 0 0 0 ); + kappaCoeffs<8> ( 0.0025219 8.506e-05 -1.312e-08 0 0 0 0 0 ); + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermophysicalProperties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermophysicalProperties new file mode 100644 index 0000000000..4b2bf263e6 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/thermophysicalProperties @@ -0,0 +1,47 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture reactingMixture; + transport polynomial; + thermo hPolynomial; + energy sensibleEnthalpy; + equationOfState icoPolynomial; + specie specie; +} + +dpdt no; + +chemistryReader foamChemistryReader; + +foamChemistryFile "/reactions"; + +foamChemistryThermoFile "/thermo.incompressiblePoly"; + +liquids +{ + H2O; +} + +solids +{} + +inertSpecie air; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/turbulenceProperties b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/turbulenceProperties new file mode 100644 index 0000000000..3f370460ea --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/constant/turbulenceProperties @@ -0,0 +1,28 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType RAS; + +RAS +{ + RASModel kEpsilon; + + printCoeffs no; + turbulence on; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/blockMeshDict b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/blockMeshDict new file mode 100644 index 0000000000..9c153e2614 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/blockMeshDict @@ -0,0 +1,105 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +scale 1; + +vertices +( + ( 0 0 0 ) + ( 1 0 0 ) + ( 2 0 0 ) + ( 4 0 0 ) + ( 4 1 0 ) + ( 2 1 0 ) + ( 2 2 0 ) + ( 1 2 0 ) + ( 1 1 0 ) + ( 0 1 0 ) + + ( 0 0 0.1 ) + ( 1 0 0.1 ) + ( 2 0 0.1 ) + ( 4 0 0.1 ) + ( 4 1 0.1 ) + ( 2 1 0.1 ) + ( 2 2 0.1 ) + ( 1 2 0.1 ) + ( 1 1 0.1 ) + ( 0 1 0.1 ) +); + +blocks +( + hex (0 1 8 9 10 11 18 19) ( 10 10 1 ) simpleGrading ( 1 1 1 ) + hex (1 2 5 8 11 12 15 18) ( 10 10 1 ) simpleGrading ( 1 1 1 ) + hex (2 3 4 5 12 13 14 15) ( 20 10 1 ) simpleGrading ( 1 1 1 ) + hex (8 5 6 7 18 15 16 17) ( 10 10 1 ) simpleGrading ( 1 1 1 ) +); + +edges +( +); + +defaultPatch +{ + name frontAndBack; + type empty; +} + +boundary +( + inlet1 + { + type patch; + faces + ( + (0 10 19 9) + ); + } + inlet2 + { + type patch; + faces + ( + (7 17 16 6) + ); + } + outlet + { + type patch; + faces + ( + (3 4 14 13) + ); + } + walls + { + type wall; + faces + ( + (0 1 11 10) + (1 2 12 11) + (2 3 13 12) + (9 19 18 8) + (8 18 17 7) + (5 6 16 15) + (5 15 14 4) + ); + } +); + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/controlDict b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/controlDict new file mode 100644 index 0000000000..f3245d3018 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/controlDict @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application reactingParcelFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 60; + +deltaT 1e-3; + +writeControl adjustable; + +writeInterval 0.5; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 8;; + +writeCompression off; + +timeFormat general; + +timePrecision 8; + +runTimeModifiable true; + +adjustTimeStep yes; + +maxCo 1.0; + +maxDeltaT 1; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/decomposeParDict b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/decomposeParDict new file mode 100644 index 0000000000..6509bf67c8 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/decomposeParDict @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 2; + +method scotch; + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSchemes b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSchemes new file mode 100644 index 0000000000..c8fd195733 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSchemes @@ -0,0 +1,60 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ 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; + + turbulence Gauss upwind; + div(phi,k) $turbulence; + div(phi,epsilon) $turbulence; + + energy Gauss linear; + div(phi,K) $energy; + div(phi,h) $energy; + + div(phi,Yi_h) Gauss upwind; + div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear orthogonal; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default orthogonal; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSolution b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSolution new file mode 100644 index 0000000000..c2c8727594 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/recycleParticles/system/fvSolution @@ -0,0 +1,97 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2011 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "rho.*" + { + solver diagonal; + } + + "(U|k|epsilon)" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-06; + relTol 0.1; + } + + "(U|k|epsilon)Final" + { + $U; + tolerance 1e-06; + relTol 0; + } + + p_rgh + { + solver GAMG; + smoother DICGaussSeidel; + tolerance 0; + relTol 0.1; + } + + p_rghFinal + { + $p_rgh; + tolerance 1e-06; + relTol 0; + } + + Phi + { + $p_rgh; + } + + Yi + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0; + } + + h + { + $Yi; + relTol 0.1; + } + + hFinal + { + $Yi; + } +} + +PIMPLE +{ + transonic no; + nOuterCorrectors 1; + nCorrectors 2; + nNonOrthogonalCorrectors 0; + momentumPredictor yes; +} + +relaxationFactors +{ + equations + { + ".*Final" 1; + } +} + + +// ************************************************************************* //