Merged reactingParcelFilmFoam into reactingParcelFoam

The combined solver includes the most advanced and general functionality from
each solver including:

    Continuous phase
    Lagrangian multiphase parcels
    Optional film
    Continuous and Lagrangian phase reactions
    Radiation
    Strong buoyancy force support by solving for p_rgh

The reactingParcelFoam and reactingParcelFilmFoam tutorials have been combined
and updated.
This commit is contained in:
Henry Weller
2017-08-29 09:33:45 +01:00
parent 940e28f636
commit 82675f0976
340 changed files with 1020 additions and 3025 deletions

View File

@ -1,20 +0,0 @@
/*--------------------------------*- 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 additionalControls;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvePrimaryRegion true; // false; // true;
// ************************************************************************* //

View File

@ -1,20 +0,0 @@
/*--------------------------------*- 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 additionalControls;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvePrimaryRegion true;
// ************************************************************************* //

View File

@ -1,20 +0,0 @@
/*--------------------------------*- 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 additionalControls;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvePrimaryRegion false; // true;
// ************************************************************************* //

View File

@ -1,20 +0,0 @@
/*--------------------------------*- 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 additionalControls;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvePrimaryRegion true; // false; // true;
// ************************************************************************* //

View File

@ -23,17 +23,22 @@ boundaryField
{
fuel
{
type zeroGradient;
type calculated;
value $internalField;
}
air
{
type zeroGradient;
type calculated;
value $internalField;
}
outlet
{
type totalPressure;
p0 $internalField;
type calculated;
value $internalField;
}
frontAndBack
{
type empty;

View File

@ -10,50 +10,37 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
location "1";
object H2O;
location "0";
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
internalField uniform 1e5;
boundaryField
{
walls
fuel
{
type zeroGradient;
type fixedFluxPressure;
}
inlet
air
{
type fixedValue;
value uniform 0;
type fixedFluxPressure;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
type prghTotalPressure;
p0 $internalField;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}

View File

@ -22,7 +22,7 @@ solvers
solver diagonal;
}
p
p_rgh
{
solver GAMG;
tolerance 1e-6;
@ -30,9 +30,9 @@ solvers
smoother GaussSeidel;
}
pFinal
p_rghFinal
{
$p;
$p_rgh;
relTol 0;
}

View File

@ -14,11 +14,11 @@ FoamFile
object scalarListList;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// (x y z) (u v w) d rho mDot T Cp (Y0..Y2)
// (x y z) (u v w) d rho mDot T Cp (Y0..Y2) (Yg0..YgN) (Yl0..YlN) (Ys0..YsN)
(
(0 1.95 -0.2) (0 -5 0) 0.001 1000 0.002 300 4200 (1)
(0 1.95 0) (0 -5 0) 0.001 1000 0.002 300 4200 (1)
(0 1.95 0.2) (0 -5 0) 0.001 1000 0.002 300 4200 (1)
(0 1.95 -0.2) (0 -5 0) 0.001 1000 0.002 300 4200 (0 1 0) () (1) ()
(0 1.95 0) (0 -5 0) 0.001 1000 0.002 300 4200 (0 1 0) () (1) ()
(0 1.95 0.2) (0 -5 0) 0.001 1000 0.002 300 4200 (0 1 0) () (1) ()
);

View File

@ -74,7 +74,7 @@ subModels
{
model1
{
type reactingLookupTableInjection;
type reactingMultiphaseLookupTableInjection;
massTotal 0.5;
parcelBasisType mass;
SOI 0;
@ -91,10 +91,14 @@ subModels
heatTransferModel none;
compositionModel singlePhaseMixture;
compositionModel singleMixtureFraction;
phaseChangeModel none;
devolatilisationModel none;
surfaceReactionModel none;
stochasticCollisionModel none;
surfaceFilmModel thermoSurfaceFilm;
@ -106,15 +110,24 @@ subModels
type rebound;
}
singlePhaseMixtureCoeffs
singleMixtureFractionCoeffs
{
phases
(
gas
{
}
liquid
{
H2O 1;
H2O 1;
}
solid
{
}
);
YGasTot0 0;
YLiquidTot0 1;
YSolidTot0 0;
}
thermoSurfaceFilmCoeffs

View File

@ -15,9 +15,9 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application reactingParcelFilmFoam;
application reactingParcelFoam;
startFrom latestTime; // startTime;
startFrom latestTime;
startTime 0;

View File

@ -97,10 +97,6 @@ PIMPLE
relaxationFactors
{
fields
{
".*Final" 1;
}
equations
{
".*Final" 1;

View File

@ -10,7 +10,6 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
location "1";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -23,33 +22,42 @@ boundaryField
{
walls
{
type zeroGradient;
type calculated;
value $internalField;
}
inlet
{
type zeroGradient;
type calculated;
value $internalField;
}
outlet
{
type fixedValue;
value uniform 100000;
type calculated;
value $internalField;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;

View File

@ -10,8 +10,7 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
location "1";
object p;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -23,33 +22,40 @@ boundaryField
{
walls
{
type zeroGradient;
type fixedFluxPressure;
}
inlet
{
type zeroGradient;
type fixedFluxPressure;
}
outlet
{
type fixedValue;
value uniform 100000;
type prghPressure;
p uniform 100000;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;

View File

@ -1,71 +0,0 @@
/*--------------------------------*- 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 "1";
object G;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 0 -3 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type MarshakRadiation;
emissivityMode lookup;
emissivity uniform 1;
value uniform 0;
refValue uniform 0;
refGradient uniform 0;
valueFraction uniform 0;
}
inlet
{
type MarshakRadiation;
emissivityMode lookup;
emissivity uniform 1;
value uniform 0;
refValue uniform 0;
refGradient uniform 0;
valueFraction uniform 0;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,60 +0,0 @@
/*--------------------------------*- 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 "1";
object N2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0.79;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform 0.79;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,60 +0,0 @@
/*--------------------------------*- 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 "1";
object O2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0.21;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform 0.21;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,61 +0,0 @@
/*--------------------------------*- 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 "1";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 350;
boundaryField
{
walls
{
type fixedValue;
value uniform 400;
}
inlet
{
type fixedValue;
value uniform 350;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,60 +0,0 @@
/*--------------------------------*- 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 "1";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type noSlip;
}
inlet
{
type fixedValue;
value uniform (5 0 0);
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,62 +0,0 @@
/*--------------------------------*- 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 "1";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type compressible::alphatWallFunction;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,61 +0,0 @@
/*--------------------------------*- 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 "1";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 0.0449;
boundaryField
{
walls
{
type epsilonWallFunction;
value uniform 0.0449;
}
inlet
{
type fixedValue;
value uniform 0.0449;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,61 +0,0 @@
/*--------------------------------*- 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 "1";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.0938;
boundaryField
{
walls
{
type kqRWallFunction;
value uniform 0.0938;
}
inlet
{
type fixedValue;
value uniform 0.0938;
}
outlet
{
type zeroGradient;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -1,62 +0,0 @@
/*--------------------------------*- 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 "1";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type nutkWallFunction;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
cycLeft_half0
{
type cyclic;
}
cycRight_half0
{
type cyclic;
}
frontAndBack
{
type empty;
}
cycLeft_half1
{
type cyclic;
}
cycRight_half1
{
type cyclic;
}
}
// ************************************************************************* //

View File

@ -12,7 +12,4 @@ rm -rf 0
# remove post-processing directory
rm -rf postProcessing
# copy 0.orig to 0
cp -r 0.orig 0
#------------------------------------------------------------------------------

View File

@ -6,6 +6,9 @@ cd ${0%/*} || exit 1 # Run from this directory
application=`getApplication`
# copy 0.orig to 0
cp -r 0.orig 0
# create mesh
runApplication blockMesh

View File

@ -47,7 +47,7 @@ solvers
relTol 0;
}
p
p_rgh
{
solver GAMG;
tolerance 0;
@ -55,9 +55,9 @@ solvers
smoother GaussSeidel;
}
pFinal
p_rghFinal
{
$p;
$p_rgh;
tolerance 1e-06;
relTol 0;
}

View File

@ -15,4 +15,3 @@ runParallel $application
runApplication -s wallFilmRegion reconstructPar -region wallFilmRegion
runApplication -s primaryRegion reconstructPar

View File

@ -112,10 +112,14 @@ subModels
heatTransferModel none;
compositionModel singlePhaseMixture;
compositionModel singleMixtureFraction;
phaseChangeModel none;
devolatilisationModel none;
surfaceReactionModel none;
stochasticCollisionModel none;
surfaceFilmModel thermoSurfaceFilm;
@ -127,15 +131,24 @@ subModels
type rebound;
}
singlePhaseMixtureCoeffs
singleMixtureFractionCoeffs
{
phases
(
gas
{
}
liquid
{
H2O 1;
}
solid
{
}
);
YGasTot0 0;
YLiquidTot0 1;
YSolidTot0 0;
}
thermoSurfaceFilmCoeffs

View File

@ -25,6 +25,6 @@ faceSet cubeFaceSet add boxToFace (0.399 0.599 0.399) (0.601 0.601 0.601)
faceSet cubeFaceSet add boxToFace (0.399 0.399 0.399) (0.601 0.401 0.601)
cellSet cubeFacesCells new faceToCell cubeFaceSet owner
faceZoneSet cubeFaces new setsToFaceZone cubeFaceSet cubeFacesCells
faceSet floorFaces new boxToFace (-0.001 -0.001 -0.001) (1.001 1.001 0.001)
cellSet floorCells new faceToCell floorFaces owner
faceZoneSet floorFaces new setsToFaceZone floorFaces floorCells
faceSet FloorFaces new boxToFace (-0.001 -0.001 -0.001) (1.001 1.001 0.001)
cellSet floorCells new faceToCell FloorFaces owner
faceZoneSet floorFaces new setsToFaceZone FloorFaces floorCells

View File

@ -15,7 +15,7 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application reactingParcelFilmFoam;
application reactingParcelFoam;
startForm latestTime;

View File

@ -78,10 +78,6 @@ PIMPLE
relaxationFactors
{
fields
{
".*Final" 1;
}
equations
{
".*Final" 1;

Some files were not shown because too many files have changed in this diff Show More