diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/U b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/U new file mode 100644 index 0000000000..a40b782ffe --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/U @@ -0,0 +1,95 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + object U; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +Uinlet (1 0 0); + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + inlet + { + type fixedValue; + value uniform $Uinlet; + } + + "(outlet|side|top)" + { + type zeroGradient; + value uniform (0 0 0); + } + + film + { + type velocityFilmShell; + active true; + U U; + pRef 1e5; + T0 300; + + thermo + { + H2O; + } + + turbulence laminar; + + laminarCoeffs + { + friction ManningStrickler; + n 0.1; // Manning number + Cf 0.9; // Gas friction + } + + injectionModels (); + + forces (perturbedTemperatureDependentContactAngle); + + perturbedTemperatureDependentContactAngleCoeffs + { + Ccf 0.085; + theta constant 45; + distribution + { + type normal; + normalDistribution + { + minValue 50; + maxValue 100; + expectation 75; + variance 100; + } + } + } + + region film; + liquidFilmModel kinematicThinFilm; + value uniform (0 0 0); + } + + cylinder + { + type fixedValue; + value uniform (0 0 0); + } + + #includeEtc "caseDicts/setConstraintTypes" +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/Uf_film b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/Uf_film new file mode 100644 index 0000000000..f8dfa4a37e --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/Uf_film @@ -0,0 +1,53 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class areaVectorField; + object Uf_film; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + inlet + { + type zeroGradient; + } + + outlet + { + type zeroGradient; + } + + side + { + type fixedValue; + value $internalField; + } + + symmetry + { + type fixedValue; + value $internalField; + } + + cylinder + { + type fixedValue; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/hf_film b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/hf_film new file mode 100644 index 0000000000..becb44a6d9 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/hf_film @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class areaScalarField; + object hf_film; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 0 0 0 0 0]; + +internalField uniform 1e-5; + +boundaryField +{ + inlet + { + type fixedValue; + value uniform 0.000141; + } + + outlet + { + type inletOutlet; + inletValue $internalField; + value $internalField; + phi phi2s_film; + } + + side + { + type zeroGradient; + } + + symmetry + { + type zeroGradient; + } + + cylinder + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/p b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/p new file mode 100644 index 0000000000..92c2d6259d --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/0/p @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + inlet + { + type zeroGradient; + } + + "(outlet|side|top)" + { + type fixedValue; + value uniform 0; + } + + film + { + type zeroGradient; + } + + cylinder + { + type zeroGradient; + } + + #includeEtc "caseDicts/setConstraintTypes" +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allclean b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allclean new file mode 100755 index 0000000000..c77ab8663b --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allclean @@ -0,0 +1,10 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/CleanFunctions # Tutorial clean functions +#------------------------------------------------------------------------------ + +cleanCase + +cleanFaMesh + +#------------------------------------------------------------------------------ diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun new file mode 100755 index 0000000000..85eceb74ad --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun @@ -0,0 +1,12 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +runApplication blockMesh + +runApplication makeFaMesh + +runApplication $(getApplication) + +#------------------------------------------------------------------------------ diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun-parallel b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun-parallel new file mode 100755 index 0000000000..3c324be8ba --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/Allrun-parallel @@ -0,0 +1,16 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +runApplication blockMesh + +runApplication decomposePar + +runParallel makeFaMesh + +runParallel $(getApplication) + +runApplication reconstructPar + +#------------------------------------------------------------------------------ diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/g b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/g new file mode 100644 index 0000000000..6987bdc182 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/g @@ -0,0 +1,21 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 (4.905 0 -8.4957); + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/transportProperties b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/transportProperties new file mode 100644 index 0000000000..cfc1f72c62 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/transportProperties @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object transportProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +transportModel Newtonian; + +nu 1e-5; + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/turbulenceProperties b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/turbulenceProperties new file mode 100644 index 0000000000..60227a5725 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/constant/turbulenceProperties @@ -0,0 +1,20 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/blockMeshDict b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/blockMeshDict new file mode 100644 index 0000000000..232fdc2f60 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/blockMeshDict @@ -0,0 +1,162 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +scale 0.016; + +vertices +( + (0.5 0 0) + (1 0 0) + (2 0 0) + (2 0.707107 0) + (0.707107 0.707107 0) + (0.353553 0.353553 0) + (2 2 0) + (0.707107 2 0) + (0 2 0) + (0 1 0) + (0 0.5 0) + (-0.5 0 0) + (-1 0 0) + (-2 0 0) + (-2 0.707107 0) + (-0.707107 0.707107 0) + (-0.353553 0.353553 0) + (-2 2 0) + (-0.707107 2 0) + + (0.5 0 0.5) + (1 0 0.5) + (2 0 0.5) + (2 0.707107 0.5) + (0.707107 0.707107 0.5) + (0.353553 0.353553 0.5) + (2 2 0.5) + (0.707107 2 0.5) + (0 2 0.5) + (0 1 0.5) + (0 0.5 0.5) + (-0.5 0 0.5) + (-1 0 0.5) + (-2 0 0.5) + (-2 0.707107 0.5) + (-0.707107 0.707107 0.5) + (-0.353553 0.353553 0.5) + (-2 2 0.5) + (-0.707107 2 0.5) +); + +blocks +( + hex (5 4 9 10 24 23 28 29) (10 10 1) simpleGrading (1 1 1) + hex (0 1 4 5 19 20 23 24) (10 10 1) simpleGrading (1 1 1) + hex (1 2 3 4 20 21 22 23) (20 10 1) simpleGrading (1 1 1) + hex (4 3 6 7 23 22 25 26) (20 20 1) simpleGrading (1 1 1) + hex (9 4 7 8 28 23 26 27) (10 20 1) simpleGrading (1 1 1) + hex (15 16 10 9 34 35 29 28) (10 10 1) simpleGrading (1 1 1) + hex (12 11 16 15 31 30 35 34) (10 10 1) simpleGrading (1 1 1) + hex (13 12 15 14 32 31 34 33) (20 10 1) simpleGrading (1 1 1) + hex (14 15 18 17 33 34 37 36) (20 20 1) simpleGrading (1 1 1) + hex (15 9 8 18 34 28 27 37) (10 20 1) simpleGrading (1 1 1) +); + +edges +( + arc 0 5 (0.469846 0.17101 0) + arc 5 10 (0.17101 0.469846 0) + arc 1 4 (0.939693 0.34202 0) + arc 4 9 (0.34202 0.939693 0) + arc 19 24 (0.469846 0.17101 0.5) + arc 24 29 (0.17101 0.469846 0.5) + arc 20 23 (0.939693 0.34202 0.5) + arc 23 28 (0.34202 0.939693 0.5) + arc 11 16 (-0.469846 0.17101 0) + arc 16 10 (-0.17101 0.469846 0) + arc 12 15 (-0.939693 0.34202 0) + arc 15 9 (-0.34202 0.939693 0) + arc 30 35 (-0.469846 0.17101 0.5) + arc 35 29 (-0.17101 0.469846 0.5) + arc 31 34 (-0.939693 0.34202 0.5) + arc 34 28 (-0.34202 0.939693 0.5) +); + +patches +( + symmetry bottom + ( + (0 1 20 19) + (1 2 21 20) + (12 11 30 31) + (13 12 31 32) + ) + patch outlet + ( + (2 3 22 21) + (3 6 25 22) + ) + patch side + ( + (7 8 27 26) + (6 7 26 25) + (8 18 37 27) + (18 17 36 37) + ) + patch inlet + ( + (14 13 32 33) + (17 14 33 36) + ) + wall cylinder + ( + (10 5 24 29) + (5 0 19 24) + (16 10 29 35) + (11 16 35 30) + ) + patch film + ( + (5 4 1 0) + (4 3 2 1) + (7 6 3 4) + (8 7 4 9) + (9 4 5 10) + (18 8 9 15) + (15 9 10 16) + (17 18 15 14) + (14 15 12 13) + (15 16 11 12) + ) + patch top + ( + (19 20 23 24) + (20 21 22 23) + (23 22 25 26) + (29 24 23 28) + (28 23 26 27) + (35 29 28 34) + (34 28 27 37) + (32 31 34 33) + (33 34 37 36) + (31 30 35 34) + ) +); + +mergePatchPairs +( +); + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/controlDict b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/controlDict new file mode 100644 index 0000000000..e25436ac45 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/controlDict @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application pimpleFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 0.5; // please increase it to 2.5 or more, if need be. + +deltaT 0.0002; + +writeControl timeStep; + +writeInterval 300; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression off; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 1; + +maxDeltaT 0.1; + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/decomposeParDict b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/decomposeParDict new file mode 100644 index 0000000000..0e5089d2b5 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/decomposeParDict @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 4; + +method simple; + +coeffs +{ + n ( 2 2 1 ); +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faMeshDefinition b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faMeshDefinition new file mode 100644 index 0000000000..2138897bd0 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faMeshDefinition @@ -0,0 +1,58 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object faMeshDefinition; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +polyMeshPatches ( film ); + +boundary +{ + inlet + { + type patch; + ownerPolyPatch film; + neighbourPolyPatch inlet; + } + + outlet + { + type patch; + ownerPolyPatch film; + neighbourPolyPatch outlet; + } + + side + { + type patch; + ownerPolyPatch film; + neighbourPolyPatch side; + } + + symmetry + { + type patch; + ownerPolyPatch film; + neighbourPolyPatch bottom; + } + + cylinder + { + type patch; + ownerPolyPatch film; + neighbourPolyPatch cylinder; + } +} + + +// ************************************************************************** // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSchemes b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSchemes new file mode 100644 index 0000000000..5a134f2062 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSchemes @@ -0,0 +1,56 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object faSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + div(phif_film,hf_film) Gauss Gamma 0.5; + div(phi2s_film,Uf_film) Gauss linearUpwind; +} + +laplacianSchemes +{ + default none; + laplacian(hf_film) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default corrected; +} + +fluxRequired +{ + hf_film true; +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSolution b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSolution new file mode 100644 index 0000000000..d08d1911b0 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/faSolution @@ -0,0 +1,52 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object faSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + Uf_film + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-08; + relTol 0.0; + } + + hf_film + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-08; + relTol 0.0; + } +} + +PIMPLE +{ + momentumPredictor true; + nOuterCorr 4; + nCorr 1; + nNonOrthCorr 0; + nFilmCorr 1; +} + +relaxationFactors +{ + hf_Film 0.9; + Uf_Film 0.9; +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSchemes b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSchemes new file mode 100644 index 0000000000..f2a58e1724 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSchemes @@ -0,0 +1,53 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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; + grad(p) Gauss linear; + grad(U) Gauss linear; +} + +divSchemes +{ + default none; + div(phi,U) Gauss linearUpwind grad(U); + div((nuEff*dev2(T(grad(U))))) Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; + interpolate(HbyA) linear; +} + +snGradSchemes +{ + default corrected; +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSolution b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSolution new file mode 100644 index 0000000000..bbc9377e10 --- /dev/null +++ b/tutorials/incompressible/pimpleFoam/laminar/inclinedPlaneFilm/system/fvSolution @@ -0,0 +1,86 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "pcorr.*" + { + solver GAMG; + smoother GaussSeidel; + cacheAgglomeration no; + tolerance 0.02; + relTol 0; + } + + p + { + $pcorr; + tolerance 1e-06; + relTol 0.01; + } + + pFinal + { + $p; + tolerance 1e-06; + relTol 0; + } + + U + { + solver smoothSolver; + smoother GaussSeidel; + tolerance 1e-05; + relTol 0.01; + } + + UFinal + { + $U; + tolerance 1e-06; + relTol 0; + } + + cellMotionUx + { + solver PCG; + preconditioner DIC; + tolerance 1e-08; + relTol 0; + } +} + +PIMPLE +{ + correctPhi no; + nOuterCorrectors 1; + nCorrectors 1; + nNonOrthogonalCorrectors 0; + + pRefCell 0; + pRefValue 0; +} + +relaxationFactors +{ + equations + { + "U.*" 1; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/T b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/T new file mode 100644 index 0000000000..8f0f5eb983 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/T @@ -0,0 +1,31 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 +{ + ".*" + { + type fixedValue; + value uniform 300; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/U b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/U new file mode 100644 index 0000000000..e646f5924a --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/U @@ -0,0 +1,87 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 (3 0 0); + +boundaryField +{ + inlet + { + type fixedValue; + value uniform (3 0 0); + inletValue uniform (0 0 0); + } + + outlet + { + type pressureInletOutletVelocity; + inletValue uniform (0 0 0); + value uniform (0 0 0); + } + + "(top|sides)" + { + type noSlip; + } + + base + { + type velocityFilmShell; + active true; + infoOutput false; + deltaWet 1e-6; + U U; + pRef 1e5; // Reference pressure for thermo + T0 300; // Isothermal film temperature + + zeroWallVelocity true; + + thermo + { + H2O; + } + + turbulence laminar; + + laminarCoeffs + { + friction ManningStrickler; // Wall friction model + n 0.005; // Manning number + Cf 0; // Gas friction + } + + injectionModels + ( + curvatureSeparation + ); + + forces (); + + curvatureSeparationCoeffs + { + definedPatchRadii 0; + } + + region film; + liquidFilmModel kinematicThinFilm; + value uniform (0 0 0); + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Uf_film b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Uf_film new file mode 100644 index 0000000000..dbd1b2b06c --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Uf_film @@ -0,0 +1,41 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class areaVectorField; + object Uf_film; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + inlet + { + type fixedValue; + value uniform (5 0 0); + } + + outlet + { + type zeroGradient; + } + + sides + { + type slip; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Ydefault b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Ydefault new file mode 100644 index 0000000000..82d34bdd1c --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/Ydefault @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object Ydefault; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1; + +boundaryField +{ + inlet + { + type fixedValue; + value uniform 1; + inletValue uniform 1; + } + + ".*" + { + type zeroGradient; + } + + outlet + { + type inletOutlet; + value uniform 1; + inletValue uniform 1; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/alphat b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/alphat new file mode 100644 index 0000000000..35e61fe774 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/alphat @@ -0,0 +1,43 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 1e-3; + +boundaryField +{ + inlet + { + type calculated; + value uniform 1e-3; + } + + outlet + { + type calculated; + value uniform 1e-3; + } + + ".*" + { + type compressible::alphatWallFunction;; + value uniform 1e-3; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/epsilon b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/epsilon new file mode 100644 index 0000000000..36db91fd77 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/epsilon @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 200; + +boundaryField +{ + inlet + { + type turbulentIntensityKineticEnergyInlet; + intensity 0.05; + value uniform 1; + } + + outlet + { + type inletOutlet; + inletValue uniform 200; + value uniform 200; + } + + ".*" + { + type epsilonWallFunction; + value uniform 200; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/hf_film b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/hf_film new file mode 100644 index 0000000000..1118d28bf8 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/hf_film @@ -0,0 +1,41 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class areaScalarField; + object hf_film; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 0 0 0 0 0]; + +internalField uniform 1e-5; + +boundaryField +{ + inlet + { + type fixedValue; + value uniform 1e-4; + } + + outlet + { + type zeroGradient; + } + + sides + { + type zeroGradient; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/k b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/k new file mode 100644 index 0000000000..366fa08842 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/k @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 1; + +boundaryField +{ + inlet + { + type turbulentIntensityKineticEnergyInlet; + intensity 0.05; + value uniform 1; + } + + outlet + { + type inletOutlet; + inletValue uniform 1; + value uniform 1; + } + + ".*" + { + type kqRWallFunction; + value uniform 1; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/nut b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/nut new file mode 100644 index 0000000000..d5b8a2c36e --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/nut @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object nut; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -1 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + inlet + { + type calculated; + value uniform 0; + } + + outlet + { + type inletOutlet; + inletValue uniform 0; + value uniform 0; + } + + ".*" + { + type nutkWallFunction; + value uniform 0; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p new file mode 100644 index 0000000000..39ec1e0fa5 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p @@ -0,0 +1,31 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 1e5; + +boundaryField +{ + ".*" + { + type calculated; + value $internalField; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p_rgh b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p_rgh new file mode 100644 index 0000000000..2600309f00 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/0/p_rgh @@ -0,0 +1,46 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 +{ + inlet + { + type fixedFluxPressure; + } + + base + { + type fixedFluxPressure; + } + + outlet + { + type totalPressure; + p0 uniform 100000; + } + + "(sides|top)" + { + type fixedFluxPressure; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allclean b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allclean new file mode 100755 index 0000000000..c77ab8663b --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allclean @@ -0,0 +1,10 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/CleanFunctions # Tutorial clean functions +#------------------------------------------------------------------------------ + +cleanCase + +cleanFaMesh + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun new file mode 100755 index 0000000000..85eceb74ad --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun @@ -0,0 +1,12 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +runApplication blockMesh + +runApplication makeFaMesh + +runApplication $(getApplication) + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun-parallel b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun-parallel new file mode 100755 index 0000000000..b5e535e97e --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/Allrun-parallel @@ -0,0 +1,16 @@ +#!/bin/sh +cd "${0%/*}" || exit # Run from this directory +. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions +#------------------------------------------------------------------------------ + +runApplication blockMesh + +runParallel -s decompose redistributePar -decompose + +runParallel makeFaMesh + +runParallel $(getApplication) + +runParallel -s reconstruct redistributePar -reconstruct + +#------------------------------------------------------------------------------ diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.dat b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.dat new file mode 100644 index 0000000000..482bba7c53 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.dat @@ -0,0 +1,36 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object foam.dat; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +air +{ + specie + { + molWeight 28.9; + } + thermodynamics + { + Cp 1005; + Hf 0; + } + transport + { + mu 1.82e-05; + Pr 0.71; + } +} + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.inp b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.inp new file mode 100644 index 0000000000..76841b250c --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/foam.inp @@ -0,0 +1,7 @@ +species +( + air +); + +reactions +{} diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/g b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/g new file mode 100644 index 0000000000..e340b6cb92 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/g @@ -0,0 +1,21 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 -9.8); + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/reactingCloud1Properties b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/reactingCloud1Properties new file mode 100644 index 0000000000..00cb119a3d --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/reactingCloud1Properties @@ -0,0 +1,172 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object reactingCloud1Properties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active true; + coupled no; + transient yes; + cellValueSourceCorrection no; + 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 4187; + + constantVolume true; +} + +subModels +{ + particleForces + { + sphereDrag; + gravity; + } + + injectionModels + { + model1 + { + type coneInjection; + SOI 5.00; + duration 10.00; + positionAxis + ( + ((-0.5 0.1 0.7) (1 0 -1)) + ((-0.5 0.2 0.7) (1 0 -1)) + ((-0.5 0.3 0.7) (1 0 -1)) + ((-0.5 0.4 0.7) (1 0 -1)) + ((-0.5 0.5 0.7) (1 0 -1)) + ((-0.5 0.6 0.7) (1 0 -1)) + ((-0.5 0.7 0.7) (1 0 -1)) + ((-0.5 0.8 0.7) (1 0 -1)) + ((-0.5 0.9 0.7) (1 0 -1)) + ); + + massTotal 10; + parcelsPerInjector 20000; + parcelBasisType mass; + flowRateProfile constant 1; + Umag constant 0.2; + thetaInner constant 0; + thetaOuter constant 55; + + sizeDistribution + { + type RosinRammler; + RosinRammlerDistribution + { + minValue 5e-04; + maxValue 0.0012; + d 7.5e-05; + n 0.5; + } + } + } + } + + dispersionModel none; + + patchInteractionModel standardWallInteraction; + + heatTransferModel none; + + compositionModel singleMixtureFraction; + + phaseChangeModel none; + + devolatilisationModel none; + + surfaceReactionModel none; + + stochasticCollisionModel none; + + surfaceFilmModel thermoSurfaceFilm; + + radiation off; + + standardWallInteractionCoeffs + { + type rebound; + } + + singleMixtureFractionCoeffs + { + phases + ( + gas + { + } + liquid + { + H2O 1; + } + solid + { + } + ); + YGasTot0 0; + YLiquidTot0 1; + YSolidTot0 0; + } + + thermoSurfaceFilmCoeffs + { + interactionType absorb; + injectionOffset 5; // Times film thickness + minDiameter 1e-4; + } +} + + +cloudFunctions +{} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/thermophysicalProperties b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/thermophysicalProperties new file mode 100644 index 0000000000..3d70875bcb --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/thermophysicalProperties @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture reactingMixture; + transport const; + thermo hConst; + energy sensibleEnthalpy; + equationOfState perfectGas; + specie specie; +} + +chemistryReader foamChemistryReader; + +foamChemistryFile "/foam.inp"; + +foamChemistryThermoFile "/foam.dat"; + +inertSpecie air; + +liquids +{ + H2O; +} + +solids +{} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/turbulenceProperties b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/turbulenceProperties new file mode 100644 index 0000000000..0b743e5074 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/constant/turbulenceProperties @@ -0,0 +1,25 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType RAS; + +RAS +{ + RASModel kEpsilon; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/blockMeshDict b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/blockMeshDict new file mode 100644 index 0000000000..2fcb06b7ec --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/blockMeshDict @@ -0,0 +1,104 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +scale 1; + +vertices +( + (0 0 0) //0 + (5 0 0) //1 + (5 1 0) //2 + (0 1 0) //3 + (-1 0 0.5) //4 + (-1 1 0.5) //5 + (-1 1 1) //6 + (-1 0 1) //7 + (5 1 1) //8 + (0 1 0.5) //9 + (5 0 0.5) //10 + (5 1 0.5) //11 + (0 0 1) //12 + (0 1 1) //13 + (0 0 0.5) //14 + (5 0 1) //15 +); + +blocks +( + hex (14 10 11 9 12 15 8 13 ) (100 20 10) simpleGrading (1 1 1) + hex (0 1 2 3 14 10 11 9) (100 20 10) simpleGrading (1 1 1) + hex (4 14 9 5 7 12 13 6) (20 20 10) simpleGrading (1 1 1) +); + +edges +( +); + +boundary +( + base + { + type wall; + faces + ( + (0 1 2 3) + (0 3 9 14) + (4 14 9 5) + ); + } + outlet + { + type patch; + faces + ( + (1 2 11 10) + (15 10 11 8) + ); + } + inlet + { + type patch; + faces + ( + (4 5 6 7) + ); + } + sides + { + type wall; + faces + ( + (0 1 10 14) + (14 10 15 12) + (4 14 12 7) + (3 2 11 9) + (9 11 8 13) + (5 9 13 6) + ); + } + top + { + type wall; + faces + ( + (12 15 8 13) + (7 12 13 6) + ); + } +); + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/controlDict b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/controlDict new file mode 100644 index 0000000000..fa5759a96b --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/controlDict @@ -0,0 +1,64 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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 1; // please increase it to 10 or more, if need be. + +deltaT 0.005; + +writeControl adjustable; + +writeInterval 0.05; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression off; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep yes; + +maxCo 0.9; + +regionFaMaxCo 4; // max finiteArea Courant number + +maxDeltaT 0.1; + +functions +{ + faRegionCo + { + type setTimeStepFaRegion; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/decomposeParDict b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/decomposeParDict new file mode 100644 index 0000000000..68dbf1e8ad --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/decomposeParDict @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 4; + +method hierarchical; + +coeffs +{ + n ( 2 2 1 ); +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faMeshDefinition b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faMeshDefinition new file mode 100644 index 0000000000..7133740322 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faMeshDefinition @@ -0,0 +1,45 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant/faMesh"; + object faMeshDefinition; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +polyMeshPatches 1( base ); + +boundary +{ + inlet + { + type patch; + ownerPolyPatch base; + neighbourPolyPatch inlet; + } + + outlet + { + type patch; + ownerPolyPatch base; + neighbourPolyPatch outlet; + } + + sides + { + type patch; + ownerPolyPatch base; + neighbourPolyPatch sides; + } +} + + +// ************************************************************************** // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSchemes b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSchemes new file mode 100644 index 0000000000..e97836da2c --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSchemes @@ -0,0 +1,57 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object faSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + div(phif_film,hf_film) Gauss upwind; + div(phi2s_film,Uf_film) Gauss upwind; + div(faceAreaNormals) Gauss linear; +} + +laplacianSchemes +{ + default none; + laplacian(hf_film) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default corrected; +} + +fluxRequired +{ + hf_film true; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSolution b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSolution new file mode 100644 index 0000000000..b44bbd9d87 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/faSolution @@ -0,0 +1,52 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object faSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + Uf_film + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-08; + relTol 0.0; + } + + hf_film + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-08; + relTol 0.0; + } +} + +PIMPLE +{ + momentumPredictor true; + nOuterCorr 4; + nCorr 1; + nNonOrthCorr 0; + nFilmCorr 1; +} + +relaxationFactors +{ + hf_Film 1; + Uf_Film 1; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSchemes b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSchemes new file mode 100644 index 0000000000..aa871a54a1 --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSchemes @@ -0,0 +1,59 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / 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; + grad(p) Gauss linear; + grad(U) Gauss linear; +} + +divSchemes +{ + default none; + div(phi,U) Gauss upwind; + div(phi,K) Gauss upwind; + div(phi,h) Gauss upwind; + div(phi,k) Gauss upwind; + div(phi,epsilon) Gauss upwind; + div(U) 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; + interpolate(HbyA) linear; +} + +snGradSchemes +{ + default corrected; +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSolution b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSolution new file mode 100644 index 0000000000..071788055a --- /dev/null +++ b/tutorials/lagrangian/reactingParcelFoam/liquidFilmStepWithSprinkles/system/fvSolution @@ -0,0 +1,102 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: v2106 | +| \\ / A nd | Website: www.openfoam.com | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "(rho|G)" + { + solver PCG; + preconditioner DIC; + tolerance 0; + relTol 0.1; + } + + "(rho|G)Final" + { + $rho; + tolerance 1e-05; + relTol 0; + } + + p_rgh + { + solver PBiCGStab; + tolerance 0; + relTol 0.1; + smoother GaussSeidel; + preconditioner DIC; + maxIter 100; + } + + p_rghFinal + { + $p_rgh; + tolerance 1e-06; + relTol 0; + maxIter 1000; + } + + "(h|Yi|O2|N2|H2O|k|epsilon)" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0.1; + } + + "(h|Yi|O2|N2|H2O|k|epsilon)Final" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0; + } + + U + { + solver smoothSolver; + smoother GaussSeidel; + tolerance 1e-05; + relTol 0.01; + } + + UFinal + { + $U; + tolerance 1e-06; + relTol 0; + } +} + +PIMPLE +{ + momentumPredictor yes; + correctPhi no; + nOuterCorrectors 1; + nCorrectors 3; + nNonOrthogonalCorrectors 0; +} + +relaxationFactors +{ + equations + { + ".*" 1; + } +} + + +// ************************************************************************* //