From 48c79c0bf50c710a8677c26b3d1c483b872512b5 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Fri, 25 May 2018 15:46:19 +0100 Subject: [PATCH] reactingTwoPhaseEulerFoam: Added evaporation test cases Two single-cell test cases have been added for reactingTwoPhaseEulerFoam with an interface composition phase system. These are droplet evaporation cases; one single- and one multi-component. The cases run for every possible inert specie, and check that the results between the runs are broadly similar. The multi-component case shows some unphysical changes at the start due to non-convergence of the pimple iteration during the initial transient. This can be mitigated by reducing the time-step. --- .../0/C3H8O.liquid | 24 +++ .../0/H2O.liquid | 24 +++ .../waterAndIsopropanolEvaporation/0/T.gas | 24 +++ .../waterAndIsopropanolEvaporation/0/T.liquid | 24 +++ .../waterAndIsopropanolEvaporation/0/U.gas | 24 +++ .../waterAndIsopropanolEvaporation/0/U.liquid | 24 +++ .../0/Ydefault.gas | 24 +++ .../0/Ydefault.liquid | 24 +++ .../waterAndIsopropanolEvaporation/0/air.gas | 24 +++ .../0/alpha.gas | 24 +++ .../0/alpha.liquid | 24 +++ .../waterAndIsopropanolEvaporation/0/p | 24 +++ .../waterAndIsopropanolEvaporation/0/p_rgh | 24 +++ .../waterAndIsopropanolEvaporation/Allclean | 7 + .../waterAndIsopropanolEvaporation/Allrun | 100 +++++++++ .../waterAndIsopropanolEvaporation/constant/g | 21 ++ .../constant/phaseProperties | 189 ++++++++++++++++++ .../thermophysicalProperties.gas.orig | 92 +++++++++ .../thermophysicalProperties.liquid.orig | 83 ++++++++ .../constant/turbulenceProperties.gas | 20 ++ .../constant/turbulenceProperties.liquid | 20 ++ .../system/blockMeshDict | 55 +++++ .../system/controlDict | 80 ++++++++ .../system/fvSchemes | 62 ++++++ .../system/fvSolution | 69 +++++++ .../waterEvaporation/0/T.gas | 24 +++ .../waterEvaporation/0/T.liquid | 24 +++ .../waterEvaporation/0/U.gas | 24 +++ .../waterEvaporation/0/U.liquid | 24 +++ .../waterEvaporation/0/Ydefault.gas | 24 +++ .../waterEvaporation/0/air.gas | 24 +++ .../waterEvaporation/0/alpha.gas | 24 +++ .../waterEvaporation/0/alpha.liquid | 24 +++ .../interfaceComposition/waterEvaporation/0/p | 24 +++ .../waterEvaporation/0/p_rgh | 24 +++ .../waterEvaporation/Allclean | 7 + .../waterEvaporation/Allrun | 78 ++++++++ .../waterEvaporation/constant/g | 21 ++ .../waterEvaporation/constant/phaseProperties | 151 ++++++++++++++ .../thermophysicalProperties.gas.orig | 80 ++++++++ .../constant/thermophysicalProperties.liquid | 52 +++++ .../constant/turbulenceProperties.gas | 20 ++ .../constant/turbulenceProperties.liquid | 20 ++ .../waterEvaporation/system/blockMeshDict | 55 +++++ .../waterEvaporation/system/controlDict | 77 +++++++ .../waterEvaporation/system/fvSchemes | 62 ++++++ .../waterEvaporation/system/fvSolution | 69 +++++++ 47 files changed, 2042 insertions(+) create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/C3H8O.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/H2O.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/air.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p_rgh create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allclean create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allrun create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/g create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/phaseProperties create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.gas.orig create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.liquid.orig create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/blockMeshDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/controlDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSchemes create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSolution create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/Ydefault.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/air.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p_rgh create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allclean create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allrun create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/g create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/phaseProperties create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.gas.orig create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.gas create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.liquid create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/blockMeshDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/controlDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSchemes create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSolution diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/C3H8O.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/C3H8O.liquid new file mode 100644 index 000000000..d8b0b16c9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/C3H8O.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object C3H8O.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.5; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/H2O.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/H2O.liquid new file mode 100644 index 000000000..e2f848236 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/H2O.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object H2O.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.5; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.gas new file mode 100644 index 000000000..c48c4f412 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object T.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 350; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.liquid new file mode 100644 index 000000000..1818cff50 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/T.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object T.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 350; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.gas new file mode 100644 index 000000000..c9ba322ef --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volVectorField; + location "0"; + object U.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.liquid new file mode 100644 index 000000000..288a3cd85 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/U.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volVectorField; + location "0"; + object U.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.gas new file mode 100644 index 000000000..89051d2cf --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object Ydefault.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.liquid new file mode 100644 index 000000000..9b0718a56 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/Ydefault.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object Ydefault.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/air.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/air.gas new file mode 100644 index 000000000..6e05699b9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/air.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object air.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.gas new file mode 100644 index 000000000..7642cd558 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object alpha.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.99; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.liquid new file mode 100644 index 000000000..b85446d4a --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/alpha.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object alpha.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.01; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p new file mode 100644 index 000000000..e9f87dd28 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p_rgh b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p_rgh new file mode 100644 index 000000000..4ed773893 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/0/p_rgh @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allclean b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allclean new file mode 100755 index 000000000..4087e9b46 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allclean @@ -0,0 +1,7 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase && rm -rf postProcessing_* postProcessing.pdf diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allrun b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allrun new file mode 100755 index 000000000..412588ec4 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/Allrun @@ -0,0 +1,100 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +gasSpecies="air H2O C3H8O none" +liquidSpecies="H2O C3H8O none" + +runApplication blockMesh + +function setArg() +{ + [ $1 != "none" ] && echo -set $1 || echo -remove +} + +for g in $gasSpecies +do + foamDictionary -entry inertSpecie $(setArg $g) constant/thermophysicalProperties.gas + for l in $liquidSpecies + do + foamDictionary -entry inertSpecie $(setArg $l) constant/thermophysicalProperties.liquid + runApplication -s $g.gas_$l.liquid $(getApplication) + mv postProcessing postProcessing_$g.gas_$l.liquid + done +done + +gnuplot << EOF + +set terminal eps enhanced size 5.83,8.27 +set output 'postProcessing.eps' + +set lmargin at screen 0.15 +set rmargin at screen 0.84 + +gasSpecies = '$gasSpecies' +liquidSpecies = '$liquidSpecies' + +set multiplot layout 4,1 + +set xlabel "Time (s)" + +set ytics nomirror +set y2tics +set ylabel 'Gas volume fraction' +set y2label 'Liquid volume fraction' + +plot \ + 'postProcessing_none.gas_none.liquid/plot/0/alpha.gas' w l lc 1 t 'Gas', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/alpha.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas_none.liquid/plot/0/alpha.liquid' w l axes x1y2 lc 2 t 'Liquid', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/alpha.liquid' w l axes x1y2 lc 2 notitle + +set ytics mirror +unset y2tics +set ylabel 'Temperature (K)' +unset y2label + +plot \ + 'postProcessing_none.gas_none.liquid/plot/0/T.gas' w l lc 1 t 'Gas', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/T.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas_none.liquid/plot/0/T.liquid' w l lc 2 t 'Liquid', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/T.liquid' w l lc 2 notitle + +set ytics nomirror +set y2tics +set ylabel "Vapour mass fraction" +set y2label "Air mass fraction" + +plot \ + 'postProcessing_none.gas_none.liquid/plot/0/H2O.gas' w l lc 1 t 'H2O', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/H2O.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas_none.liquid/plot/0/C3H8O.gas' w l lc 2 t 'C3H8O', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/C3H8O.gas' w l lc 2 notitle, \ + 'postProcessing_none.gas_none.liquid/plot/0/air.gas' w l axes x1y2 lc 3 t 'Air', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/air.gas' w l axes x1y2 lc 3 notitle + +set ytics mirror +unset y2tics +set ylabel "Liquid mass fraction" +unset y2label + +plot \ + 'postProcessing_none.gas_none.liquid/plot/0/H2O.liquid' w l lc 1 t 'H2O', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/H2O.liquid' w l lc 1 notitle, \ + 'postProcessing_none.gas_none.liquid/plot/0/C3H8O.liquid' w l lc 2 t 'C3H8O', \ + for [g in gasSpecies] for [l in liquidSpecies] \ + 'postProcessing_'.g.'.gas_'.l.'.liquid/plot/0/C3H8O.liquid' w l lc 2 notitle + +unset multiplot + +EOF diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/g b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/g new file mode 100644 index 000000000..d182180d9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/g @@ -0,0 +1,21 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/phaseProperties b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/phaseProperties new file mode 100644 index 000000000..e10341367 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/phaseProperties @@ -0,0 +1,189 @@ +/*--------------------------------*- 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 phaseProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +type interfaceCompositionPhaseChangeTwoPhaseSystem; + +phases (gas liquid); + +nInterfaceCorrectors 1; + +gas +{ + type multiComponentPhaseModel; + diameterModel constant; + d 1e-3; + residualAlpha 1e-6; + Sc 0.7; +} + +liquid +{ + type multiComponentPhaseModel; + diameterModel constant; + d 5e-4; + residualAlpha 1e-6; + Sc 1; +} + +blending +{ + default + { + type none; + continuousPhase gas; + } +} + +surfaceTension +( + (gas and liquid) + { + type constant; + sigma 0.07; + } +); + +aspectRatio +( + (gas in liquid) + { + type constant; + E0 1.0; + } + + (liquid in gas) + { + type constant; + E0 1.0; + } +); + +drag +( +); + +virtualMass +( +); + +interfaceComposition +( + (gas in liquid) + { + type nonRandomTwoLiquid; + species (H2O C3H8O); + Le 1.0; + H2O + { + alpha 0.3; + beta 0.0; + interaction + { + type Antoine; + A 6.8284; + B -1283.46; + C 0; + } + type saturated; + species (H2O); + Le 1.0; + saturationPressure + { + type ArdenBuck; + } + } + C3H8O + { + alpha 0.3; + beta 0.0; + interaction + { + type Antoine; + A -1.3115; + B 426.398; + C 0; + } + type saturated; + species (C3H8O); + Le 1.0; + saturationPressure + { + type Antoine; + A 22.7057916; + B -3125.59117503; + C -75.814; + } + } + } +); + +heatTransfer.gas +( + (liquid in gas) + { + type RanzMarshall; + residualAlpha 1e-4; + } +); + +heatTransfer.liquid +( + (liquid in gas) + { + type spherical; + residualAlpha 1e-4; + } +); + +massTransfer.gas +( + (liquid in gas) + { + type Frossling; + Le 1.0; + } +); + +massTransfer.liquid +( + (liquid in gas) + { + type spherical; + Le 1.0; + } +); + +phaseTransfer +( +); + +lift +( +); + +wallLubrication +( +); + +turbulentDispersion +( +); + +// Minimum allowable pressure +pMin 10000; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.gas.orig b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.gas.orig new file mode 100644 index 000000000..05fe3eaa6 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.gas.orig @@ -0,0 +1,92 @@ +/*--------------------------------*- 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 thermophysicalProperties.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture multiComponentMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +species +( + air + H2O + C3H8O +); + +inertSpecie air; + +air +{ + specie + { + molWeight 28.9596; + } + thermodynamics + { + Hf 0; + Cp 1012.5; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + +H2O +{ + specie + { + molWeight 18.0153; + } + thermodynamics + { + Hf -1.3435e+07; + Cp 1857.8; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + +C3H8O +{ + specie + { + molWeight 60; + } + thermodynamics + { + Hf -4.538e+06; + Cp 1490; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.liquid.orig b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.liquid.orig new file mode 100644 index 000000000..a11ee89b2 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/thermophysicalProperties.liquid.orig @@ -0,0 +1,83 @@ +/*--------------------------------*- 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 thermophysicalProperties.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture multiComponentMixture; + transport const; + thermo hConst; + equationOfState perfectFluid; + specie specie; + energy sensibleInternalEnergy; +} + +species +( + H2O + C3H8O +); + +inertSpecie H2O; + +H2O +{ + specie + { + molWeight 18.0153; + } + equationOfState + { + R 3000; + rho0 1027; + } + thermodynamics + { + Hf -1.5879e+07; + Cp 4195; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + +C3H8O +{ + specie + { + molWeight 60; + } + equationOfState + { + R 3000; + rho0 786; + } + thermodynamics + { + Hf -5.303e+06; + Cp 2680; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.gas new file mode 100644 index 000000000..1db1a2670 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.gas @@ -0,0 +1,20 @@ +/*--------------------------------*- 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 turbulenceProperties.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.liquid new file mode 100644 index 000000000..8a2a00b53 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/constant/turbulenceProperties.liquid @@ -0,0 +1,20 @@ +/*--------------------------------*- 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 turbulenceProperties.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/blockMeshDict b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/blockMeshDict new file mode 100644 index 000000000..d801695b3 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/blockMeshDict @@ -0,0 +1,55 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +vertices +( + (-1 -1 -1) + (1 -1 -1) + (1 1 -1) + (-1 1 -1) + + (-1 -1 1) + (1 -1 1) + (1 1 1) + (-1 1 1) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (1 1 1) simpleGrading (1 1 1) +); + +edges +( +); + +defaultPatch +{ + name defaultPatch; + type empty; +} + +patches +( +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/controlDict b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/controlDict new file mode 100644 index 000000000..3c7704970 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/controlDict @@ -0,0 +1,80 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application reactingTwoPhaseEulerFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 1; + +deltaT 0.005; + +writeControl runTime; + +writeInterval 1; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 10; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 0.5; + +maxDeltaT 1; + +functions +{ + plot + { + lib ("libsampling.so"); + type probes; + writeControl timeStep; + writeInterval 1; + fields + ( + air.gas + alpha.gas + alpha.liquid + H2O.gas + H2O.liquid + C3H8O.gas + C3H8O.liquid + p + T.gas + T.liquid + ); + probeLocations ((0 0 0)); + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSchemes b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSchemes new file mode 100644 index 000000000..4e5905487 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSchemes @@ -0,0 +1,62 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + "div\(phi,alpha.*\)" Gauss vanLeer; + "div\(phir,alpha.*\)" Gauss vanLeer; + + "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; + "div\(phi.*,U.*\)" Gauss limitedLinearV 1; + + "div\(alphaRhoPhi.*,Yi\)" Gauss limitedLinear 1; + "div\(alphaRhoPhi.*,(h|e).*\)" Gauss limitedLinear 1; + "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; + "div\(alphaPhi.*,p\)" Gauss limitedLinear 1; + + "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear uncorrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default uncorrected; +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSolution b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSolution new file mode 100644 index 000000000..3f9f83487 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterAndIsopropanolEvaporation/system/fvSolution @@ -0,0 +1,69 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "alpha.*" + { + nAlphaCorr 1; + nAlphaSubCycles 3; + } + + p_rgh + { + solver PCG; + preconditioner none; + tolerance 1e-11; + relTol 0.001; + } + + p_rghFinal + { + $p_rgh; + relTol 0; + } + + "(U|e|h|Yi).*" + { + solver smoothSolver; + smoother symGaussSeidel; + tolerance 1e-12; + relTol 0; + minIter 1; + residualAlpha 0; + } +} + +PIMPLE +{ + nOuterCorrectors 5; + nCorrectors 1; + nNonOrthogonalCorrectors 0; + pRefCell 0; + pRefValue 1e5; +} + +relaxationFactors +{ + equations + { + ".*" 1; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.gas new file mode 100644 index 000000000..6f10c2c58 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object T.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 350; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.liquid new file mode 100644 index 000000000..9dd39afbb --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/T.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object T.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 350; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.gas new file mode 100644 index 000000000..c9ba322ef --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volVectorField; + location "0"; + object U.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.liquid new file mode 100644 index 000000000..288a3cd85 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/U.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class volVectorField; + location "0"; + object U.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/Ydefault.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/Ydefault.gas new file mode 100644 index 000000000..89051d2cf --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/Ydefault.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object Ydefault.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/air.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/air.gas new file mode 100644 index 000000000..6e05699b9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/air.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object air.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.gas new file mode 100644 index 000000000..7642cd558 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.gas @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object alpha.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.99; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.liquid new file mode 100644 index 000000000..b85446d4a --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/alpha.liquid @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object alpha.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.01; + +boundaryField {} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p new file mode 100644 index 000000000..e9f87dd28 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p_rgh b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p_rgh new file mode 100644 index 000000000..4ed773893 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/0/p_rgh @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + location "0"; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField {} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allclean b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allclean new file mode 100755 index 000000000..4087e9b46 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allclean @@ -0,0 +1,7 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase && rm -rf postProcessing_* postProcessing.pdf diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allrun b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allrun new file mode 100755 index 000000000..9b2c10df4 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/Allrun @@ -0,0 +1,78 @@ +#!/bin/sh + +cd ${0%/*} || exit 1 + +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +gasSpecies="air H2O none" + +runApplication blockMesh + +function setArg() +{ + [ $1 != "none" ] && echo -set $1 || echo -remove +} + +for g in $gasSpecies +do + foamDictionary -entry inertSpecie $(setArg $g) constant/thermophysicalProperties.gas + runApplication -s $g.gas $(getApplication) + mv postProcessing postProcessing_$g.gas +done + +gnuplot << EOF + +set terminal eps enhanced size 5.83,8.27 +set output 'postProcessing.eps' + +set lmargin at screen 0.15 +set rmargin at screen 0.84 + +gasSpecies = '$gasSpecies' + +set multiplot layout 3,1 + +set xlabel "Time (s)" + +set ytics nomirror +set y2tics +set ylabel 'Gas volume fraction' +set y2label 'Liquid volume fraction' + +plot \ + 'postProcessing_none.gas/plot/0/alpha.gas' w l lc 1 t 'Gas', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/alpha.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas/plot/0/alpha.liquid' w l axes x1y2 lc 2 t 'Liquid', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/alpha.liquid' w l axes x1y2 lc 2 notitle + +set ytics mirror +unset y2tics +set ylabel 'Temperature (K)' +unset y2label + +plot \ + 'postProcessing_none.gas/plot/0/T.gas' w l lc 1 t 'Gas', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/T.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas/plot/0/T.liquid' w l lc 2 t 'Liquid', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/T.liquid' w l lc 2 notitle + +set ytics nomirror +set y2tics +set ylabel "Vapour mass fraction" +set y2label "Air mass fraction" + +plot \ + 'postProcessing_none.gas/plot/0/H2O.gas' w l lc 1 t 'H2O', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/H2O.gas' w l lc 1 notitle, \ + 'postProcessing_none.gas/plot/0/air.gas' w l axes x1y2 lc 2 t 'Air', \ + for [g in gasSpecies] \ + 'postProcessing_'.g.'.gas/plot/0/air.gas' w l axes x1y2 lc 2 notitle + +unset multiplot + +EOF diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/g b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/g new file mode 100644 index 000000000..d182180d9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/g @@ -0,0 +1,21 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/phaseProperties b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/phaseProperties new file mode 100644 index 000000000..ce6bfa882 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/phaseProperties @@ -0,0 +1,151 @@ +/*--------------------------------*- 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 phaseProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +type interfaceCompositionPhaseChangeTwoPhaseSystem; + +phases (gas liquid); + +nInterfaceCorrectors 1; + +gas +{ + type multiComponentPhaseModel; + diameterModel constant; + d 1e-3; + residualAlpha 1e-6; + Sc 0.7; +} + +liquid +{ + type purePhaseModel; + diameterModel constant; + d 5e-4; + residualAlpha 1e-6; +} + +blending +{ + default + { + type none; + continuousPhase gas; + } +} + +surfaceTension +( + (gas and liquid) + { + type constant; + sigma 0.07; + } +); + +aspectRatio +( + (gas in liquid) + { + type constant; + E0 1.0; + } + + (liquid in gas) + { + type constant; + E0 1.0; + } +); + +drag +( +); + +virtualMass +( +); + +interfaceComposition +( + (gas in liquid) + { + type saturated; + species ( H2O ); + Le 1.0; + saturationPressure + { + type ArdenBuck; + } + } +); + +heatTransfer.gas +( + (liquid in gas) + { + type RanzMarshall; + residualAlpha 1e-4; + } +); + +heatTransfer.liquid +( + (liquid in gas) + { + type spherical; + residualAlpha 1e-4; + } +); + +massTransfer.gas +( + (liquid in gas) + { + type Frossling; + Le 1.0; + } +); + +massTransfer.liquid +( + (liquid in gas) + { + type spherical; + Le 1.0; + } +); + +phaseTransfer +( +); + +lift +( +); + +wallLubrication +( +); + +turbulentDispersion +( +); + +// Minimum allowable pressure +pMin 10000; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.gas.orig b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.gas.orig new file mode 100644 index 000000000..68649a6af --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.gas.orig @@ -0,0 +1,80 @@ +/*--------------------------------*- 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 thermophysicalProperties.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture multiComponentMixture; + transport sutherland; + thermo janaf; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +inertSpecie air; + +species +( + air + H2O +); + +H2O +{ + specie + { + molWeight 18.0153; + } + thermodynamics + { + Tlow 200; + Thigh 3500; + Tcommon 1000; + highCpCoeffs ( 3.03399 0.00217692 -1.64073e-07 -9.7042e-11 1.68201e-14 -30004.3 4.96677 ); + lowCpCoeffs ( 4.19864 -0.00203643 6.5204e-06 -5.48797e-09 1.77198e-12 -30293.7 -0.849032 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + +air +{ + specie + { + molWeight 28.9596; + } + thermodynamics + { + Tlow 200; + Thigh 3500; + Tcommon 1000; + highCpCoeffs ( 3.57304 -7.24383e-04 1.67022e-06 -1.26501e-10 -4.20580e-13 -1047.41 3.12431 ); + lowCpCoeffs ( 3.09589 1.22835e-03 -4.14267e-07 6.56910e-11 -3.87021e-15 -983.191 5.34161 ); + } + transport + { + As 1.67212e-06; + Ts 170.672; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.liquid new file mode 100644 index 000000000..3772231dd --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/thermophysicalProperties.liquid @@ -0,0 +1,52 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectFluid; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + molWeight 18.0153; + } + equationOfState + { + R 3000; + rho0 1027; + } + thermodynamics + { + Hf -1.5879e+07; + Cp 4195; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.gas b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.gas new file mode 100644 index 000000000..1db1a2670 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.gas @@ -0,0 +1,20 @@ +/*--------------------------------*- 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 turbulenceProperties.gas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.liquid b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.liquid new file mode 100644 index 000000000..8a2a00b53 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/constant/turbulenceProperties.liquid @@ -0,0 +1,20 @@ +/*--------------------------------*- 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 turbulenceProperties.liquid; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/blockMeshDict b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/blockMeshDict new file mode 100644 index 000000000..d801695b3 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/blockMeshDict @@ -0,0 +1,55 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +vertices +( + (-1 -1 -1) + (1 -1 -1) + (1 1 -1) + (-1 1 -1) + + (-1 -1 1) + (1 -1 1) + (1 1 1) + (-1 1 1) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (1 1 1) simpleGrading (1 1 1) +); + +edges +( +); + +defaultPatch +{ + name defaultPatch; + type empty; +} + +patches +( +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/controlDict b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/controlDict new file mode 100644 index 000000000..e8b4b9ab6 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/controlDict @@ -0,0 +1,77 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application reactingTwoPhaseEulerFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 1; + +deltaT 0.005; + +writeControl runTime; + +writeInterval 1; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 10; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 0.5; + +maxDeltaT 1; + +functions +{ + plot + { + lib ("libsampling.so"); + type probes; + writeControl timeStep; + writeInterval 1; + fields + ( + air.gas + alpha.gas + alpha.liquid + H2O.gas + p + T.gas + T.liquid + ); + probeLocations ((0 0 0)); + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSchemes b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSchemes new file mode 100644 index 000000000..4e5905487 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSchemes @@ -0,0 +1,62 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + "div\(phi,alpha.*\)" Gauss vanLeer; + "div\(phir,alpha.*\)" Gauss vanLeer; + + "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; + "div\(phi.*,U.*\)" Gauss limitedLinearV 1; + + "div\(alphaRhoPhi.*,Yi\)" Gauss limitedLinear 1; + "div\(alphaRhoPhi.*,(h|e).*\)" Gauss limitedLinear 1; + "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; + "div\(alphaPhi.*,p\)" Gauss limitedLinear 1; + + "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear uncorrected; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default uncorrected; +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSolution b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSolution new file mode 100644 index 000000000..3f9f83487 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/interfaceComposition/waterEvaporation/system/fvSolution @@ -0,0 +1,69 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "alpha.*" + { + nAlphaCorr 1; + nAlphaSubCycles 3; + } + + p_rgh + { + solver PCG; + preconditioner none; + tolerance 1e-11; + relTol 0.001; + } + + p_rghFinal + { + $p_rgh; + relTol 0; + } + + "(U|e|h|Yi).*" + { + solver smoothSolver; + smoother symGaussSeidel; + tolerance 1e-12; + relTol 0; + minIter 1; + residualAlpha 0; + } +} + +PIMPLE +{ + nOuterCorrectors 5; + nCorrectors 1; + nNonOrthogonalCorrectors 0; + pRefCell 0; + pRefValue 1e5; +} + +relaxationFactors +{ + equations + { + ".*" 1; + } +} + + +// ************************************************************************* //