diff --git a/bin/tools/CleanFunctions b/bin/tools/CleanFunctions index 3ac00c2865..16b45561c5 100644 --- a/bin/tools/CleanFunctions +++ b/bin/tools/CleanFunctions @@ -39,7 +39,9 @@ cleanTimeDirectories() rm -rf ./${timeDir} ./-${timeDir} > /dev/null 2>&1 zeros="0$zeros" done - rm -rf ./[1-9]* ./-[1-9]* ./log ./log.* ./log-* ./logSummary.* ./*.xml ./ParaView* ./paraFoam* ./*.OpenFOAM ./*.blockMesh ./.setSet > /dev/null 2>&1 + rm -rf ./[1-9]* ./-[1-9]* ./log ./log.* ./log-* ./logSummary.* ./*.xml \ + ./*.dat ./*.eps ./ParaView* ./paraFoam* ./*.OpenFOAM ./*.blockMesh \ + ./.setSet > /dev/null 2>&1 } cleanDynamicCode() diff --git a/test/Allclean b/test/Allclean new file mode 120000 index 0000000000..8775d6c59f --- /dev/null +++ b/test/Allclean @@ -0,0 +1 @@ +../tutorials/Allclean \ No newline at end of file diff --git a/test/Allrun b/test/Allrun new file mode 120000 index 0000000000..4920f4195c --- /dev/null +++ b/test/Allrun @@ -0,0 +1 @@ +../tutorials/Allrun \ No newline at end of file diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/README b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/README new file mode 100644 index 0000000000..abdaf191b9 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/README @@ -0,0 +1,11 @@ +The cases contained in this directory can be used for testing the +populationBalanceModel class which is a part of the reactingEulerFoam framework. +They represent single-cell setups, i.e. they solve the population balance +equation without spatial transport. The setups already represent the +inhomogeneous case, i.e. the populationBalance is split over multiple +velocityGroups. This can help to identify wether the populationBalance-induced +mass transfer terms are handled correctly. The cases can also be used to study +the influence of certain parameters on the form of the size distribution, before +running a simulation with spatial transport. Further test cases are contained in + +test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air1 new file mode 100644 index 0000000000..9e24e22827 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air1 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air2 new file mode 100644 index 0000000000..2deff50b59 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air2 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air3 new file mode 100644 index 0000000000..2ab28ff567 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air3 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.water new file mode 100644 index 0000000000..c288479186 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air1 new file mode 100644 index 0000000000..dc7beaee58 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air1 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air2 new file mode 100644 index 0000000000..c298549e6c --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air2 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air3 new file mode 100644 index 0000000000..3386be15c4 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air3 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.water new file mode 100644 index 0000000000..65c08b1231 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air1 new file mode 100644 index 0000000000..4986d6f250 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air1 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air2 new file mode 100644 index 0000000000..883a29cd81 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air2 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air3 new file mode 100644 index 0000000000..7d63260a5f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air3 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.05; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.water new file mode 100644 index 0000000000..7ca1240054 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object alpha.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.95; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air1.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air1.bubbles new file mode 100644 index 0000000000..a0d485d3db --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air1.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air1.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air2.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air2.bubbles new file mode 100644 index 0000000000..62c6b255f5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air2.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air2.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air3.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air3.bubbles new file mode 100644 index 0000000000..40817be164 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air3.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air3.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p new file mode 100644 index 0000000000..dbbc5c2dad --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p_rgh b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p_rgh new file mode 100644 index 0000000000..c76174de22 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p_rgh @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/Allrun b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/Allrun new file mode 100755 index 0000000000..076ba12d80 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/Allrun @@ -0,0 +1,39 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # run from this directory + +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +# Set application name +application=`getApplication` + +runApplication blockMesh +runApplication $application + +# Get and transpose results +cp postProcessing/box.all.numberDensity.volume.bubbles/0/box.all.numberDensity.volume.bubbles.dat . + +sed -i -e 1d box.all.numberDensity.volume.bubbles.dat +sed -i -e 's/# /#/g' box.all.numberDensity.volume.bubbles.dat + +awk ' +{ + for (i=1; i<=NF; i++) { + a[NR,i] = $i + } +} +NF>p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/g b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/phaseProperties b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/phaseProperties new file mode 100644 index 0000000000..8ecbc1b6cf --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/phaseProperties @@ -0,0 +1,211 @@ +/*--------------------------------*- 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 populationBalanceMultiphaseSystem; + +phases (air1 air2 air3 water); + +populationBalances (bubbles); + +air1 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 2.2723832645E-02; value 0.0;} + f1 {d 2.6012294874E-02; value 0.0;} + f2 {d 2.9776644423E-02; value 0.0;} + f3 {d 3.4085748967E-02; value 0.0;} + f4 {d 3.9018442311E-02; value 0.0;} + f5 {d 4.4664966635E-02; value 0.0;} + f6 {d 5.1128623451E-02; value 0.0;} + f7 {d 5.8527663466E-02; value 0.0;} + f8 {d 6.6997449953E-02; value 0.0;} + f9 {d 7.6692935176E-02; value 0.0;} + f10{d 8.7791495199E-02; value 0.0;} + f11{d 1.0049617493E-01; value 0.0;} + f12{d 1.1503940276E-01; value 0.0;} + f13{d 1.3168724280E-01; value 0.0;} + f14{d 1.5074426239E-01; value 0.0;} + f15{d 1.7255910415E-01; value 0.0;} + f16{d 1.9753086420E-01; value 0.0;} + f17{d 2.2611639359E-01; value 0.0;} + f18{d 2.5883865622E-01; value 0.0;} + f19{d 2.9629629630E-01; value 1.0;} + ); + } + + residualAlpha 1e-6; +} + +air2 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f20{d 3.3917459039E-01; value 1.0;} + f21{d 3.8825798433E-01; value 0.0;} + f22{d 4.4444444444E-01; value 0.0;} + f23{d 5.0876188558E-01; value 0.0;} + f24{d 5.8238697649E-01; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +air3 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f25{d 6.6666666667E-01; value 0.0;} + f26{d 7.6314282837E-01; value 0.0;} + f27{d 8.7358046474E-01; value 0.0;} + f28{d 1.0; value 1.0;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + (); + + binaryBreakupModels + ( + powerLawUniformBinary + { + power 2.0; + } + ); + + breakupModels + (); + + driftModels + (); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +(); + +virtualMass +(); + +heatTransfer.air1 +(); + +heatTransfer.air2 +(); + +heatTransfer.air3 +(); + +heatTransfer.water +(); + +massTransfer.air1 +(); + +massTransfer.air2 +(); + +massTransfer.air3 +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air1 new file mode 100644 index 0000000000..5254f4787a --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air1 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air2 new file mode 100644 index 0000000000..5dacc2e880 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air2 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air3 new file mode 100644 index 0000000000..0219c745d5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air3 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..bd7ad0b326 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.water @@ -0,0 +1,53 @@ +/*--------------------------------*- 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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 18; + } + equationOfState + { + rho 998.2; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airI b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air1 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airI rename to test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air1 index 6bc136989b..fda3c5b86f 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airI +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air1 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object turbulenceProperties.airI; + object turbulenceProperties.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airI b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air2 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airI rename to test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air2 index 6bc136989b..f6bb961594 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airI +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air2 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object turbulenceProperties.airI; + object turbulenceProperties.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airII b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air3 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airII rename to test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air3 index 8c3b8d3691..24af3b4c66 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airII +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air3 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object turbulenceProperties.airII; + object turbulenceProperties.air3; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airII b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.water similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airII rename to test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.water index 8c3b8d3691..7f0d75a82f 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airII +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.water @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object turbulenceProperties.airII; + object turbulenceProperties.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/plot b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/plot new file mode 100755 index 0000000000..9a081b8735 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/plot @@ -0,0 +1,27 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "numberDensity.eps" +set decimalsign "." + +set format xy '%g' +set xtics 1e1 mirror +set xlabel 'v(m^{3})' +set ytics 1e1 mirror +set ylabel 'n(m^{-3}m^{-3})' + +set logscale xy +set xrange [1e-5:1.5] +set yrange [1e-5:10] +set key center + +N0 = 0.05 +v0 = 1.0 + +n(x,t) = N0*exp(-t*x**2.0)*2.0*t*v0 + +plot "p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/g b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/phaseProperties b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/phaseProperties new file mode 100644 index 0000000000..5ec6f774c7 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/phaseProperties @@ -0,0 +1,212 @@ +/*--------------------------------*- 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 populationBalanceMultiphaseSystem; + +phases (air1 air2 air3 water); + +populationBalances (bubbles); + +air1 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 2.2723832645E-02; value 0.0;} + f1 {d 2.6012294874E-02; value 0.0;} + f2 {d 2.9776644423E-02; value 0.0;} + f3 {d 3.4085748967E-02; value 0.0;} + f4 {d 3.9018442311E-02; value 0.0;} + f5 {d 4.4664966635E-02; value 0.0;} + f6 {d 5.1128623451E-02; value 0.0;} + f7 {d 5.8527663466E-02; value 0.0;} + f8 {d 6.6997449953E-02; value 0.0;} + f9 {d 7.6692935176E-02; value 0.0;} + f10{d 8.7791495199E-02; value 0.0;} + f11{d 1.0049617493E-01; value 0.0;} + f12{d 1.1503940276E-01; value 0.0;} + f13{d 1.3168724280E-01; value 0.0;} + f14{d 1.5074426239E-01; value 0.0;} + f15{d 1.7255910415E-01; value 0.0;} + f16{d 1.9753086420E-01; value 0.0;} + f17{d 2.2611639359E-01; value 0.0;} + f18{d 2.5883865622E-01; value 0.0;} + f19{d 2.9629629630E-01; value 1.0;} + ); + } + + residualAlpha 1e-6; +} + +air2 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f20{d 3.3917459039E-01; value 1.0;} + f21{d 3.8825798433E-01; value 0.0;} + f22{d 4.4444444444E-01; value 0.0;} + f23{d 5.0876188558E-01; value 0.0;} + f24{d 5.8238697649E-01; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +air3 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f25{d 6.6666666667E-01; value 0.0;} + f26{d 7.6314282837E-01; value 0.0;} + f27{d 8.7358046474E-01; value 0.0;} + f28{d 1.0; value 1.0;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + (); + + binaryBreakupModels + (); + + breakupModels + ( + powerLaw + { + power 2.0; + daughterSizeDistributionModel uniformBinary; + } + ); + + driftModels + (); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +(); + +virtualMass +(); + +heatTransfer.air1 +(); + +heatTransfer.air2 +(); + +heatTransfer.air3 +(); + +heatTransfer.water +(); + +massTransfer.air1 +(); + +massTransfer.air2 +(); + +massTransfer.air3 +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air1 new file mode 100644 index 0000000000..5254f4787a --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air1 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air2 new file mode 100644 index 0000000000..5dacc2e880 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air2 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air3 new file mode 100644 index 0000000000..0219c745d5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air3 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..bd7ad0b326 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.water @@ -0,0 +1,53 @@ +/*--------------------------------*- 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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 18; + } + equationOfState + { + rho 998.2; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air1 new file mode 100644 index 0000000000..fda3c5b86f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air1 @@ -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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air2 new file mode 100644 index 0000000000..f6bb961594 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air2 @@ -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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air3 new file mode 100644 index 0000000000..24af3b4c66 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air3 @@ -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.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.water new file mode 100644 index 0000000000..7f0d75a82f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/plot b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/plot new file mode 100755 index 0000000000..9a081b8735 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/plot @@ -0,0 +1,27 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "numberDensity.eps" +set decimalsign "." + +set format xy '%g' +set xtics 1e1 mirror +set xlabel 'v(m^{3})' +set ytics 1e1 mirror +set ylabel 'n(m^{-3}m^{-3})' + +set logscale xy +set xrange [1e-5:1.5] +set yrange [1e-5:10] +set key center + +N0 = 0.05 +v0 = 1.0 + +n(x,t) = N0*exp(-t*x**2.0)*2.0*t*v0 + +plot "p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/g b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/phaseProperties b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/phaseProperties new file mode 100644 index 0000000000..e6b4288c8b --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/phaseProperties @@ -0,0 +1,253 @@ +/*--------------------------------*- 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 populationBalanceMultiphaseSystem; + +phases (air1 air2 air3 water); + +populationBalances (bubbles); + +air1 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 0.021336488 ; value 0.0000000005 ;} + f1 {d 0.022984035 ; value 0.0000000002 ;} + f2 {d 0.024758801 ; value 0.0000000004 ;} + f3 {d 0.02667061 ; value 0.0000000008 ;} + f4 {d 0.028730043 ; value 0.0000000015 ;} + f5 {d 0.030948501 ; value 0.0000000029 ;} + f6 {d 0.033338262 ; value 0.0000000057 ;} + f7 {d 0.035912554 ; value 0.0000000110 ;} + f8 {d 0.038685626 ; value 0.0000000215 ;} + f9 {d 0.041672828 ; value 0.0000000420 ;} + f10 {d 0.044890693 ; value 0.0000000819 ;} + f11 {d 0.048357033 ; value 0.0000001595 ;} + f12 {d 0.052091034 ; value 0.0000003107 ;} + f13 {d 0.056113366 ; value 0.0000006046 ;} + f14 {d 0.060446291 ; value 0.0000011757 ;} + f15 {d 0.065113793 ; value 0.0000022836 ;} + f16 {d 0.070141707 ; value 0.0000044295 ;} + f17 {d 0.075557864 ; value 0.0000085770 ;} + f18 {d 0.081392241 ; value 0.0000165724 ;} + f19 {d 0.087677134 ; value 0.0000319345 ;} + f20 {d 0.09444733 ; value 0.0000613299 ;} + f21 {d 0.101740302 ; value 0.0001172884 ;} + f22 {d 0.109596418 ; value 0.0002231264 ;} + f23 {d 0.118059162 ; value 0.0004216851 ;} + f24 {d 0.127175377 ; value 0.0007904101 ;} + f25 {d 0.136995522 ; value 0.0014663924 ;} + f26 {d 0.147573953 ; value 0.0026857456 ;} + f27 {d 0.158969221 ; value 0.0048406166 ;} + f28 {d 0.171244403 ; value 0.0085509222 ;} + f29 {d 0.184467441 ; value 0.0147306117 ;} + f30 {d 0.198711527 ; value 0.0245922381 ;} + f31 {d 0.214055503 ; value 0.0394763633 ;} + f32 {d 0.230584301 ; value 0.0603362728 ;} + f33 {d 0.248389408 ; value 0.0867356518 ;} + f34 {d 0.267569379 ; value 0.1154885945 ;} + f35 {d 0.288230376 ; value 0.1397279858 ;} + f36 {d 0.310486761 ; value 0.1499792809 ;} + f37 {d 0.334461724 ; value 0.1386064986 ;} + f38 {d 0.36028797 ; value 0.1062408386 ;} + f39 {d 0.388108451 ; value 0.0644531062 ;} + f40 {d 0.418077155 ; value 0.0291911897 ;} + f41 {d 0.450359963 ; value 0.0091744326 ;} + f42 {d 0.485135563 ; value 0.0018262321 ;} + f43 {d 0.522596444 ; value 0.0002053972 ;} + f44 {d 0.562949953 ; value 0.0000113164 ;} + f45 {d 0.606419454 ; value 0.0000002555 ;} + f46 {d 0.653245554 ; value 0.0000000019 ;} + ); + } + + residualAlpha 1e-6; +} + +air2 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f47 {d 0.703687442 ; value 1.0000000000 ;} + f48 {d 0.758024318 ; value 0.0000000000 ;} + f49 {d 0.816556943 ; value 0.0000000000 ;} + f50 {d 0.879609302 ; value 0.0000000000 ;} + f51 {d 0.947530397 ; value 0.0000000000 ;} + f52 {d 1.020696179 ; value 0.0000000000 ;} + f53 {d 1.099511628 ; value 0.0000000000 ;} + f54 {d 1.184412996 ; value 0.0000000000 ;} + f55 {d 1.275870223 ; value 0.0000000000 ;} + ); + } + + residualAlpha 1e-6; +} + +air3 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f56 {d 1.374389535 ; value 1.0000000000 ;} + f57 {d 1.480516246 ; value 0.0000000000 ;} + f58 {d 1.594837779 ; value 0.0000000000 ;} + f59 {d 1.717986918 ; value 0.0000000000 ;} + f60 {d 1.850645307 ; value 0.0000000000 ;} + f61 {d 1.993547224 ; value 0.0000000000 ;} + f62 {d 2.147483648 ; value 0.0000000000 ;} + f63 {d 2.313306634 ; value 0.0000000000 ;} + f64 {d 2.49193403 ; value 0.0000000000 ;} + f65 {d 2.68435456 ; value 0.0000000000 ;} + f66 {d 2.891633292 ; value 0.0000000000 ;} + f67 {d 3.114917538 ; value 0.0000000000 ;} + f68 {d 3.3554432 ; value 0.0000000000 ;} + f69 {d 3.614541615 ; value 0.0000000000 ;} + f70 {d 3.893646922 ; value 0.0000000000 ;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + ( + constant + { + rate 1.0; + } + ); + + binaryBreakupModels + (); + + breakupModels + (); + + driftModels + (); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +(); + +virtualMass +(); + +heatTransfer.air1 +(); + +heatTransfer.air2 +(); + +heatTransfer.air3 +(); + +heatTransfer.water +(); + +massTransfer.air1 +(); + +massTransfer.air2 +(); + +massTransfer.air3 +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air1 new file mode 100644 index 0000000000..5254f4787a --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air1 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air2 new file mode 100644 index 0000000000..5dacc2e880 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air2 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air3 new file mode 100644 index 0000000000..0219c745d5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air3 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..bd7ad0b326 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.water @@ -0,0 +1,53 @@ +/*--------------------------------*- 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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 18; + } + equationOfState + { + rho 998.2; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air1 new file mode 100644 index 0000000000..fda3c5b86f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air1 @@ -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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air2 new file mode 100644 index 0000000000..f6bb961594 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air2 @@ -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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air3 new file mode 100644 index 0000000000..24af3b4c66 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air3 @@ -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.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.water new file mode 100644 index 0000000000..7f0d75a82f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/plot b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/plot new file mode 100755 index 0000000000..e43f9fc684 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/plot @@ -0,0 +1,47 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "numberDensity.eps" +set decimalsign "." + +set format xy '%g' +set xtics 1e2 mirror +set xlabel 'v(m^{3})' +set ytics 1e5 mirror +set ylabel 'n(m^{-3}m^{-3})' + +set logscale xy +set xrange [1e-5:1e2] +set yrange [1e-15:100] +set key center left + +C = 1 +N0 = 2.5 +v0 = 0.01 + +# Dimensionless volume +X(x) = x/v0 + +# Initial condition +n0(x) = (N0/v0)*X(x)*exp(-X(x)) + +T(t) = C*N0*t + +# For solution of quadratic saddle point equation +p(x) = -1/X(x) +q(t) = -(T(t)/(T(t) + 2)) + +# Saddle point calculation +y_s(t,x) = -p(x)/2 + sqrt((p(x)/2)**2 - q(T(t))) + +# Dimensionless spectrum function +phi(x,t) = (8*exp((y_s(t,x) - 1)*X(x)*2))/(((T(t) + 2)**2)*y_s(t,x)*(4*pi*(2 - 1/(y_s(t,x)*X(x))))) + +# Number density at time t +n(x,t) = (N0/v0)*phi(x,t) + +plot n0(x) ls -1 t 'Initial Condition',\ +n(x,10.0) ls 2 lc rgb 'black' t 'Scott (1968)',\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::0::46 u 1:2 w p pt 1 t 'air1',\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::47::55 u 1:2 w p pt 5 t 'air2',\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::56::70 u 1:2 w p pt 9 t 'air3' diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/blockMeshDict b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/blockMeshDict new file mode 100644 index 0000000000..c925e3542d --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/blockMeshDict @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +vertices +( + (0 0 0) + (1 0 0) + (1 1 0) + (0 1 0) + (0 0 1) + (1 0 1) + (1 1 1) + (0 1 1) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (1 2 1) simpleGrading (1 1 1) +); + +edges +( +); + +patches +( +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/controlDict b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/controlDict new file mode 100644 index 0000000000..c430fd85df --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/controlDict @@ -0,0 +1,72 @@ +/*--------------------------------*- 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 reactingMultiphaseEulerFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 10; + +deltaT 0.2; + +writeControl runTime; + +writeInterval 10; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 0.5; + +maxDeltaT 1; + +functions +{ + box.all.numberDensity.volume.bubbles + { + type sizeDistribution; + functionObjectLibs ("libreactingEulerFoamFunctionObjects.so"); + + writeControl outputTime; + writeInterval 1; + log yes; + + functionType numberDensity; + abszissaType volume; + regionType all; + populationBalance bubbles; + } +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSchemes b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSchemes new file mode 100644 index 0000000000..809d86bb65 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSchemes @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + "div\(phi,alpha.*)" Gauss vanLeer; + "div\(phir,alpha.*,alpha.*)" Gauss vanLeer; + "div\(phir,alpha.*)" Gauss vanLeer; + + "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; + "div\(phi.*,U.*\)" Gauss limitedLinearV 1; + + "div\(alphaRhoPhi.*,(k|epsilon|omega|f).*\)" Gauss limitedLinear 1; + + "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear uncorrected; +} + +interpolationSchemes +{ + default linear; +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSolution b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSolution new file mode 100644 index 0000000000..db2be3b84b --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSolution @@ -0,0 +1,74 @@ +/*--------------------------------*- 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 1; + } + + bubbles + { + nCorr 1; + tolerance 1e-4; + renormalize false; + } + + p_rgh + { + solver PCG; + preconditioner none; + tolerance 1e-8; + relTol 0; + } + + p_rghFinal + { + $p_rgh; + relTol 0; + } + + "f.*" + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-6; + relTol 0; + } +} + +PIMPLE +{ + nOuterCorrectors 1; + nCorrectors 1; + + pRefCell 0; + pRefValue 1e5; +} + +relaxationFactors +{ + equations + { + f 1.0; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air1 new file mode 100644 index 0000000000..9e24e22827 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air1 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air2 new file mode 100644 index 0000000000..2deff50b59 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air2 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air3 new file mode 100644 index 0000000000..2ab28ff567 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air3 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.water new file mode 100644 index 0000000000..c288479186 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air1 new file mode 100644 index 0000000000..dc7beaee58 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air1 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air2 new file mode 100644 index 0000000000..c298549e6c --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air2 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air3 new file mode 100644 index 0000000000..3386be15c4 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air3 @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.water new file mode 100644 index 0000000000..65c08b1231 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air1 new file mode 100644 index 0000000000..25f2298937 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air1 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air2 new file mode 100644 index 0000000000..883a29cd81 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air2 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air3 new file mode 100644 index 0000000000..e12aa013f5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air3 @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.water new file mode 100644 index 0000000000..9b1c583191 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object alpha.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air1.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air1.bubbles new file mode 100644 index 0000000000..a0d485d3db --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air1.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air1.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air2.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air2.bubbles new file mode 100644 index 0000000000..62c6b255f5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air2.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air2.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air3.bubbles b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air3.bubbles new file mode 100644 index 0000000000..40817be164 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air3.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air3.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p new file mode 100644 index 0000000000..dbbc5c2dad --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p_rgh b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p_rgh new file mode 100644 index 0000000000..c76174de22 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p_rgh @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/Allrun b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/Allrun new file mode 100755 index 0000000000..6d4e4d38af --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/Allrun @@ -0,0 +1,59 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # run from this directory + +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +# Set application name +application=`getApplication` + +runApplication blockMesh +runApplication $application + +# Get and transpose results +cp postProcessing/box.all.numberDensity.volume.bubbles/0/box.all.numberDensity.volume.bubbles.dat . +cp postProcessing/box.all.volumeDensity.volume.bubbles/0/box.all.volumeDensity.volume.bubbles.dat . + +sed -i -e 1d box.all.numberDensity.volume.bubbles.dat +sed -i -e 's/# /#/g' box.all.numberDensity.volume.bubbles.dat +sed -i -e 1d box.all.volumeDensity.volume.bubbles.dat +sed -i -e 's/# /#/g' box.all.volumeDensity.volume.bubbles.dat + +awk ' +{ + for (i=1; i<=NF; i++) { + a[NR,i] = $i + } +} +NF>p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.volumeDensity.volume.bubbles.dat > box.all.volumeDensity.volume.bubbles.transposed.dat + +awk ' +{ + for (i=1; i<=NF; i++) { + a[NR,i] = $i + } +} +NF>p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/g b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/phaseProperties b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/phaseProperties new file mode 100644 index 0000000000..a43d0afd93 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/phaseProperties @@ -0,0 +1,217 @@ +/*--------------------------------*- 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 populationBalanceMultiphaseSystem; + +phases (air1 air2 air3 water); + +populationBalances (bubbles); + +air1 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 1; value 1.0;} + f1 {d 1.080059739; value 0.0;} + f2 {d 1.16652904; value 0.0;} + f3 {d 1.25992105; value 0.0;} + f4 {d 1.36079; value 0.0;} + f5 {d 1.469734492; value 0.0;} + f6 {d 1.587401052; value 0.0;} + f7 {d 1.714487966; value 0.0;} + f8 {d 1.851749425; value 0.0;} + f9 {d 2; value 0.0;} + f10{d 2.160119478; value 0.0;} + f11{d 2.333058079; value 0.0;} + f12{d 2.5198421; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +air2 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f13{d 2.72158; value 1.0;} + f14{d 2.939468985; value 0.0;} + f15{d 3.174802104; value 0.0;} + f16{d 3.428975931; value 0.0;} + f17{d 3.703498849; value 0.0;} + f18{d 4; value 0.0;} + f19{d 4.320238956; value 0.0;} + f20{d 4.666116158; value 0.0;} + f21{d 5.0396842; value 0.0;} + f22{d 5.443160001; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +air3 +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f23{d 5.878937969; value 1.0;} + f24{d 6.349604208; value 0.0;} + f25{d 6.857951863; value 0.0;} + f26{d 7.406997698; value 0.0;} + f27{d 8; value 0.0;} + f28{d 8.640477911; value 0.0;} + f29{d 9.332232317; value 0.0;} + f30{d 10.0793684 ; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + ( + hydrodynamic{} + ); + + binaryBreakupModels + (); + + breakupModels + ( + exponential + { + C 0.1; + exponent 0.01; + daughterSizeDistributionModel uniformBinary; + } + ); + + driftModels + (); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +(); + +virtualMass +(); + +heatTransfer.air1 +(); + +heatTransfer.air2 +(); + +heatTransfer.air3 +(); + +heatTransfer.water +(); + +massTransfer.air1 +(); + +massTransfer.air2 +(); + +massTransfer.air3 +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air1 new file mode 100644 index 0000000000..5254f4787a --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air1 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air2 new file mode 100644 index 0000000000..5dacc2e880 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air2 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air3 new file mode 100644 index 0000000000..0219c745d5 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air3 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState perfectGas; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 28.9; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..bd7ad0b326 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.water @@ -0,0 +1,53 @@ +/*--------------------------------*- 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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + nMoles 1; + molWeight 18; + } + equationOfState + { + rho 998.2; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air1 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air1 new file mode 100644 index 0000000000..fda3c5b86f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air1 @@ -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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air2 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air2 new file mode 100644 index 0000000000..f6bb961594 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air2 @@ -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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air3 b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air3 new file mode 100644 index 0000000000..24af3b4c66 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air3 @@ -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.air3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.water b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.water new file mode 100644 index 0000000000..7f0d75a82f --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/plot b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/plot new file mode 100755 index 0000000000..a888f94708 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/plot @@ -0,0 +1,25 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "probabilityDensity.eps" +set decimalsign "." + +set format x '%g' +set format y "%.1f" +set xtics 1 mirror +set xlabel 'd(m)' +set ytics 0.1 mirror +set ylabel 'Probability density function (m^{-3})' + +set xrange [1:10] +set yrange [0:0.4] +set key left + +plot 'rigorous.numberDensity.dat' w l t "Number density, Vanni (2000)",\ +'rigorous.volumeDensity.dat' w l t "Volume density, Vanni (2000)",\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::0::12 u 1:2 w p pt 1 t 'Number density, air1',\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::13::22 u 1:2 w p pt 5 t 'Number density, air2',\ +'box.all.numberDensity.volume.bubbles.transposed.dat' every ::23::30 u 1:2 w p pt 9 t 'Number density, air3',\ +'box.all.volumeDensity.volume.bubbles.transposed.dat' every ::0::12 u 1:2 w p pt 2 t 'Volume density, air1',\ +'box.all.volumeDensity.volume.bubbles.transposed.dat' every ::13::22 u 1:2 w p pt 6 t 'Volume density, air2',\ +'box.all.volumeDensity.volume.bubbles.transposed.dat' every ::23::30 u 1:2 w p pt 10 t 'Volume density, air3' diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/blockMeshDict b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/blockMeshDict new file mode 100644 index 0000000000..c925e3542d --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/blockMeshDict @@ -0,0 +1,48 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +vertices +( + (0 0 0) + (1 0 0) + (1 1 0) + (0 1 0) + (0 0 1) + (1 0 1) + (1 1 1) + (0 1 1) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (1 2 1) simpleGrading (1 1 1) +); + +edges +( +); + +patches +( +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/controlDict b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/controlDict new file mode 100644 index 0000000000..b1ac53709b --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/controlDict @@ -0,0 +1,89 @@ +/*--------------------------------*- 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 reactingMultiphaseEulerFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 10; + +deltaT 0.1; + +writeControl runTime; + +writeInterval 10; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 0.5; + +maxDeltaT 1; + +functions +{ + box.all.numberDensity.volume.bubbles + { + type sizeDistribution; + functionObjectLibs ("libreactingEulerFoamFunctionObjects.so"); + + writeControl outputTime; + writeInterval 1; + log yes; + + functionType numberDensity; + abszissaType diameter; + regionType all; + populationBalance bubbles; + normalize true; + } + + box.all.volumeDensity.volume.bubbles + { + type sizeDistribution; + functionObjectLibs ("libreactingEulerFoamFunctionObjects.so"); + + writeControl outputTime; + writeInterval 1; + log yes; + + functionType volumeDensity; + abszissaType diameter; + regionType all; + populationBalance bubbles; + normalize true; + } +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSchemes b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSchemes new file mode 100644 index 0000000000..809d86bb65 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSchemes @@ -0,0 +1,54 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + "div\(phi,alpha.*)" Gauss vanLeer; + "div\(phir,alpha.*,alpha.*)" Gauss vanLeer; + "div\(phir,alpha.*)" Gauss vanLeer; + + "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; + "div\(phi.*,U.*\)" Gauss limitedLinearV 1; + + "div\(alphaRhoPhi.*,(k|epsilon|omega|f).*\)" Gauss limitedLinear 1; + + "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear uncorrected; +} + +interpolationSchemes +{ + default linear; +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSolution b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSolution new file mode 100644 index 0000000000..d0808e00b9 --- /dev/null +++ b/test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSolution @@ -0,0 +1,74 @@ +/*--------------------------------*- 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 1; + } + + bubbles + { + nCorr 1; + tolerance 1e-4; + renormalize true; + } + + p_rgh + { + solver PCG; + preconditioner none; + tolerance 1e-8; + relTol 0; + } + + p_rghFinal + { + $p_rgh; + relTol 0; + } + + "f.*" + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-6; + relTol 0; + } +} + +PIMPLE +{ + nOuterCorrectors 1; + nCorrectors 1; + + pRefCell 0; + pRefValue 1e5; +} + +relaxationFactors +{ + equations + { + f 1.0; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/README b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/README new file mode 100644 index 0000000000..52d47def23 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/README @@ -0,0 +1,8 @@ +The cases contained in this directory can be used for testing the +populationBalanceModel class which is a part of the reactingEulerFoam framework. +They represent single-cell setups, i.e. they solve the population balance +equation without spatial transport. The cases can also be used to study the +influence of certain parameters on the form of the size distribution, before +running a simulation with spatial transport. Further test cases are contained in + +test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.air new file mode 100644 index 0000000000..854dfc6185 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.air @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.water new file mode 100644 index 0000000000..c288479186 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object T.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 1 0 0 0]; + +internalField uniform 300; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.air new file mode 100644 index 0000000000..726323aa3e --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.air @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.water new file mode 100644 index 0000000000..65c08b1231 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.water @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object U.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/alpha.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/alpha.air new file mode 100644 index 0000000000..f9d89cc088 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/alpha.air @@ -0,0 +1,26 @@ +/*--------------------------------*- 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.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 0.0751237; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/f.air.bubbles b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/f.air.bubbles new file mode 100644 index 0000000000..61a9ba40e9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/f.air.bubbles @@ -0,0 +1,26 @@ +/*--------------------------------*- 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 f.air.bubbles; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 0 0 0 0 0 0]; + +internalField uniform 1.0; + +boundaryField +{ +} + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p new file mode 100644 index 0000000000..dbbc5c2dad --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p_rgh b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p_rgh new file mode 100644 index 0000000000..c76174de22 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p_rgh @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + object p_rgh; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [1 -1 -2 0 0 0 0]; + +internalField uniform 1e5; + +boundaryField +{ +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/Allrun b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/Allrun new file mode 100755 index 0000000000..076ba12d80 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/Allrun @@ -0,0 +1,39 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # run from this directory + +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +# Set application name +application=`getApplication` + +runApplication blockMesh +runApplication $application + +# Get and transpose results +cp postProcessing/box.all.numberDensity.volume.bubbles/0/box.all.numberDensity.volume.bubbles.dat . + +sed -i -e 1d box.all.numberDensity.volume.bubbles.dat +sed -i -e 's/# /#/g' box.all.numberDensity.volume.bubbles.dat + +awk ' +{ + for (i=1; i<=NF; i++) { + a[NR,i] = $i + } +} +NF>p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/fvOptions b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/fvOptions new file mode 100644 index 0000000000..0455a1a643 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/fvOptions @@ -0,0 +1,43 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object fvOptions; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +injector1 +{ + timeStart 0; + duration 6; + selectionMode all; +} + +options +{ + massSource1 + { + type scalarSemiImplicitSource; + + $injector1; + + volumeMode specific; + injectionRateSuSp + { + thermo:rho.air (0.024535 0); // kg/s/m³ + thermo:rho.water (-24.535 0); // kg/s/m³ + } + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/g b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/phaseProperties b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/phaseProperties new file mode 100644 index 0000000000..7dfd21d74a --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/phaseProperties @@ -0,0 +1,172 @@ +/*--------------------------------*- 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 populationBalanceTwoPhaseSystem; + +phases (air water); + +populationBalances (bubbles); + +air +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 1.0000000 ; value 0.0000000 ;} + f1 {d 1.0500000 ; value 0.0000000 ;} + f2 {d 1.1000000 ; value 0.0000000 ;} + f3 {d 1.1500000 ; value 0.0000000 ;} + f4 {d 1.2000000 ; value 0.0000000 ;} + f5 {d 1.2500000 ; value 0.0609671 ;} + f6 {d 1.3000000 ; value 0.0741729 ;} + f7 {d 1.3500000 ; value 0.0895739 ;} + f8 {d 1.4000000 ; value 0.1074333 ;} + f9 {d 1.4500000 ; value 0.1280343 ;} + f10 {d 1.5000000 ; value 0.1516810 ;} + f11 {d 1.5500000 ; value 0.1786995 ;} + f12 {d 1.6000000 ; value 0.2094380 ;} + f13 {d 1.6500000 ; value 0.0000000 ;} + f14 {d 1.7000000 ; value 0.0000000 ;} + f15 {d 1.7500000 ; value 0.0000000 ;} + f16 {d 1.8000000 ; value 0.0000000 ;} + f17 {d 1.8500000 ; value 0.0000000 ;} + f18 {d 1.9000000 ; value 0.0000000 ;} + f19 {d 1.9500000 ; value 0.0000000 ;} + f20 {d 2.0000000 ; value 0.0000000 ;} + f21 {d 2.0500000 ; value 0.0000000 ;} + f22 {d 2.1000000 ; value 0.0000000 ;} + f23 {d 2.1500000 ; value 0.0000000 ;} + f24 {d 2.2000000 ; value 0.0000000 ;} + f25 {d 2.2500000 ; value 0.0000000 ;} + f26 {d 2.3000000 ; value 0.0000000 ;} + f27 {d 2.3500000 ; value 0.0000000 ;} + f28 {d 2.4000000 ; value 0.0000000 ;} + f29 {d 2.4500000 ; value 0.0000000 ;} + f30 {d 2.5000000 ; value 0.0000000 ;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + (); + + binaryBreakupModels + (); + + breakupModels + (); + + driftModels + ( + constant{} + ); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +( + (air in water) + { + type SchillerNaumann; + residualRe 1e-3; + swarmCorrection + { + type none; + } + } +); + +virtualMass +(); + +heatTransfer.air +(); + +heatTransfer.water +(); + +massTransfer.air +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.air new file mode 100644 index 0000000000..4c8bb3629e --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.air @@ -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.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + molWeight 28.9; + } + equationOfState + { + rho 1; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..4a566672f6 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + molWeight 18; + } + equationOfState + { + rho 1000; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.air new file mode 100644 index 0000000000..1296429b72 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.air @@ -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.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.water new file mode 100644 index 0000000000..7f0d75a82f --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/plot b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/plot new file mode 100755 index 0000000000..517570ccd9 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/plot @@ -0,0 +1,35 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "numberDensity.eps" +set decimalsign "." + +set format xy '%g' +set xtics 1 mirror +set xlabel 'v(m^3)' +set ytics 0.005 mirror +set ylabel 'n(m^{-3}m^{-3})' + +set xrange [1:15] +set yrange [0:1.5e-2] +set key top right + +plot "p { p = NF } +END { + for(j=1; j<=p; j++) { + str=a[1,j] + for(i=2; i<=NR; i++){ + str=str" "a[i,j]; + } + print str + } +}' box.all.numberDensity.volume.bubbles.dat > box.all.numberDensity.volume.bubbles.transposed.dat + +# Create plot +./plot + +#------------------------------------------------------------------------------ diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/fvOptions b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/fvOptions new file mode 100644 index 0000000000..c203d8c1a2 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/fvOptions @@ -0,0 +1,43 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object fvOptions; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +injector1 +{ + timeStart 0; + duration 4; + selectionMode all; +} + +options +{ + massSource1 + { + type scalarSemiImplicitSource; + + $injector1; + + volumeMode specific; + injectionRateSuSp + { + thermo:rho.air (-0.0258575 0); // kg/s/m³ + thermo:rho.water (25.8575 0); // kg/s/m³ + } + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/g b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/g new file mode 100644 index 0000000000..508d658494 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: dev | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 0); + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/phaseProperties b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/phaseProperties new file mode 100644 index 0000000000..14e2b7db72 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/phaseProperties @@ -0,0 +1,166 @@ +/*--------------------------------*- 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 populationBalanceTwoPhaseSystem; + +phases (air water); + +populationBalances (bubbles); + +air +{ + type pureIsothermalPhaseModel; + diameterModel velocityGroup; + velocityGroupCoeffs + { + populationBalance bubbles; + + formFactor 1; + + sizeGroups + ( + f0 {d 1.0000000 ; value 0.0000000 ;} + f1 {d 1.0500000 ; value 0.0000000 ;} + f2 {d 1.1000000 ; value 0.0000000 ;} + f3 {d 1.1500000 ; value 0.0000000 ;} + f4 {d 1.2000000 ; value 0.0000000 ;} + f5 {d 1.2500000 ; value 0.0000000 ;} + f6 {d 1.3000000 ; value 0.0000000 ;} + f7 {d 1.3500000 ; value 0.0000000 ;} + f8 {d 1.4000000 ; value 0.0000000 ;} + f9 {d 1.4500000 ; value 0.0000000 ;} + f10 {d 1.5000000 ; value 0.0000000 ;} + f11 {d 1.5500000 ; value 0.0000000 ;} + f12 {d 1.6000000 ; value 0.0000000 ;} + f13 {d 1.6500000 ; value 0.0000000 ;} + f14 {d 1.7000000 ; value 0.0000000 ;} + f15 {d 1.7500000 ; value 0.0000000 ;} + f16 {d 1.8000000 ; value 0.0000000 ;} + f17 {d 1.8500000 ; value 0.0000000 ;} + f18 {d 1.9000000 ; value 0.0000000 ;} + f19 {d 1.9500000 ; value 0.0000000 ;} + f20 {d 2.0000000 ; value 0.2064732 ;} + f21 {d 2.0500000 ; value 0.2336032 ;} + f22 {d 2.1000000 ; value 0.2635129 ;} + f23 {d 2.1500000 ; value 0.2964107 ;} + f24 {d 2.2000000 ; value 0.0000000 ;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + diameterModel constant; + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + (); + + binaryBreakupModels + (); + + breakupModels + (); + + driftModels + ( + constant{} + ); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +(); + +interfaceCompression +(); + +aspectRatio +(); + +drag +( + (air in water) + { + type SchillerNaumann; + residualRe 1e-3; + swarmCorrection + { + type none; + } + } +); + +virtualMass +(); + +heatTransfer.air +(); + +heatTransfer.water +(); + +massTransfer.air +(); + +massTransfer.water +(); + +lift +(); + +wallLubrication +(); + +turbulentDispersion +(); + +bubbleInducedTurbulence +(); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.air new file mode 100644 index 0000000000..4c8bb3629e --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.air @@ -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.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + molWeight 28.9; + } + equationOfState + { + rho 1; + } + thermodynamics + { + Cp 1007; + Hf 0; + } + transport + { + mu 1.84e-05; + Pr 0.7; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.water new file mode 100644 index 0000000000..4a566672f6 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type heRhoThermo; + mixture pureMixture; + transport const; + thermo hConst; + equationOfState rhoConst; + specie specie; + energy sensibleInternalEnergy; +} + +mixture +{ + specie + { + molWeight 18; + } + equationOfState + { + rho 1000; + } + thermodynamics + { + Cp 4195; + Hf 0; + } + transport + { + mu 3.645e-4; + Pr 2.289; + } +} + + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.air b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.air new file mode 100644 index 0000000000..1296429b72 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.air @@ -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.air; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.water b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.water new file mode 100644 index 0000000000..7f0d75a82f --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.water @@ -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.water; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/plot b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/plot new file mode 100755 index 0000000000..ec7e486289 --- /dev/null +++ b/test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/plot @@ -0,0 +1,35 @@ +#!/usr/bin/gnuplot + +set terminal postscript eps monochrome +set output "numberDensity.eps" +set decimalsign "." + +set format xy '%g' +set xtics 1 mirror +set xlabel 'v(m^3)' +set ytics 0.005 mirror +set ylabel 'n(m^{-3}m^{-3})' + +set xrange [1:15] +set yrange [0:1.5e-2] +set key top right + +plot " # Extracts useful info from log file. logReport() @@ -67,7 +69,7 @@ logReport() } -# Recursively run all tutorials +# Recursively run all cases foamRunTutorials -test -skipFirst diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air1 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air1 index 3d377d0950..ec92afac21 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air1 @@ -10,7 +10,7 @@ FoamFile version 2.0; format ascii; class volScalarField; - object T.airI; + object T.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -27,7 +27,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airI; + phi phi.air1; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air2 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air2 index 3d377d0950..d1cafcb7e0 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.air2 @@ -10,7 +10,7 @@ FoamFile version 2.0; format ascii; class volScalarField; - object T.airI; + object T.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -27,7 +27,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airI; + phi phi.air2; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air1 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air1 index da969ff410..ced8937c0a 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air1 @@ -10,7 +10,7 @@ FoamFile version 2.0; format binary; class volVectorField; - object U.airI; + object U.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -28,7 +28,7 @@ boundaryField outlet { type pressureInletOutletVelocity; - phi phi.airI; + phi phi.air1; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air2 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air2 index da969ff410..a83d1299ea 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.air2 @@ -10,7 +10,7 @@ FoamFile version 2.0; format binary; class volVectorField; - object U.airI; + object U.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -28,7 +28,7 @@ boundaryField outlet { type pressureInletOutletVelocity; - phi phi.airI; + phi phi.air2; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airI.orig b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airI.orig rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig index 203cc22fc5..4a30789394 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airI.orig +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object alpha.airI; + object alpha.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airI; + phi phi.air1; inletValue uniform 0; value uniform 0; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airII.orig b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air2.orig similarity index 94% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airII.orig rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air2.orig index 019c483a0e..5c8b4c69cc 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.airII.orig +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air2.orig @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object alpha.airII; + object alpha.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airII; + phi phi.air2; inletValue uniform 1; value uniform 1; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/p_rgh b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/p_rgh index 8d1b65707c..5dcbfa90da 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/p_rgh +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/p_rgh @@ -29,9 +29,9 @@ boundaryField { type prghTotalPressure; p0 $internalField; - U U.airII; - phi phi.airII; - rho thermo:rho.airII; + U U.air2; + phi phi.air2; + rho thermo:rho.air2; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties index 62d61e05e5..6c09d3376f 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties @@ -17,9 +17,9 @@ FoamFile type heatAndMomentumTransferMultiphaseSystem; -phases (airI airII water); +phases (air1 air2 water); -airI +air1 { type purePhaseModel; diameterModel isothermal; @@ -32,7 +32,7 @@ airI residualAlpha 1e-6; } -airII +air2 { type purePhaseModel; diameterModel isothermal; @@ -62,10 +62,10 @@ blending default { type linear; - minFullyContinuousAlpha.airI 0.7; - minPartlyContinuousAlpha.airI 0.5; - minFullyContinuousAlpha.airII 0.7; - minPartlyContinuousAlpha.airII 0.5; + minFullyContinuousAlpha.air1 0.7; + minPartlyContinuousAlpha.air1 0.5; + minFullyContinuousAlpha.air2 0.7; + minPartlyContinuousAlpha.air2 0.5; minFullyContinuousAlpha.water 0.7; minPartlyContinuousAlpha.water 0.5; } @@ -73,13 +73,13 @@ blending surfaceTension ( - (airI and water) + (air1 and water) { type constant; sigma 0.07; } - (airII and water) + (air2 and water) { type constant; sigma 0.07; @@ -91,25 +91,25 @@ interfaceCompression aspectRatio ( - (airI in water) + (air1 in water) { type constant; E0 1.0; } - (water in airI) + (water in air1) { type constant; E0 1.0; } - (airII in water) + (air2 in water) { type constant; E0 1.0; } - (water in airII) + (water in air2) { type constant; E0 1.0; @@ -118,7 +118,7 @@ aspectRatio drag ( - (airI in water) + (air1 in water) { type SchillerNaumann; residualRe 1e-3; @@ -128,7 +128,7 @@ drag } } - (water in airI) + (water in air1) { type SchillerNaumann; residualRe 1e-3; @@ -138,7 +138,7 @@ drag } } - (airI and water) + (air1 and water) { type segregated; m 0.5; @@ -149,7 +149,7 @@ drag } } - (airII in water) + (air2 in water) { type SchillerNaumann; residualRe 1e-3; @@ -159,7 +159,7 @@ drag } } - (water in airII) + (water in air2) { type SchillerNaumann; residualRe 1e-3; @@ -169,7 +169,7 @@ drag } } - (airII and water) + (air2 and water) { type segregated; m 0.5; @@ -180,7 +180,7 @@ drag } } - (airI in airII) + (air1 in air2) { type SchillerNaumann; residualRe 1e-3; @@ -190,7 +190,7 @@ drag } } - (airII in airI) + (air2 in air1) { type SchillerNaumann; @@ -201,7 +201,7 @@ drag } } - (airI and airII) + (air1 and air2) { type segregated; m 0.5; @@ -215,37 +215,37 @@ drag virtualMass ( - (airI in water) + (air1 in water) { type constantCoefficient; Cvm 0.5; } - (water in airI) + (water in air1) { type constantCoefficient; Cvm 0.5; } - (airII in water) + (air2 in water) { type constantCoefficient; Cvm 0.5; } - (water in airII) + (water in air2) { type constantCoefficient; Cvm 0.5; } - (airI in airII) + (air1 in air2) { type constantCoefficient; Cvm 0.5; } - (airI in airII) + (air1 in air2) { type constantCoefficient; Cvm 0.5; @@ -254,37 +254,37 @@ virtualMass heatTransfer ( - (airI in water) + (air1 in water) { type RanzMarshall; residualAlpha 1e-4; } - (water in airI) + (water in air1) { type RanzMarshall; residualAlpha 1e-4; } - (airII in water) + (air2 in water) { type RanzMarshall; residualAlpha 1e-4; } - (water in airII) + (water in air2) { type RanzMarshall; residualAlpha 1e-4; } - (airI in airII) + (air1 in air2) { type RanzMarshall; residualAlpha 1e-4; } - (airII in airI) + (air2 in air1) { type RanzMarshall; residualAlpha 1e-4; diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 similarity index 96% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 index e11ccef1bd..7426132d89 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object thermophysicalProperties.airI; + object thermophysicalProperties.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airI b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air2 similarity index 96% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airI rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air2 index e11ccef1bd..62bff79a0a 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airI +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air2 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object thermophysicalProperties.airI; + object thermophysicalProperties.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air1 b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air1 new file mode 100644 index 0000000000..fda3c5b86f --- /dev/null +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air1 @@ -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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air2 b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air2 new file mode 100644 index 0000000000..f6bb961594 --- /dev/null +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air2 @@ -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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/controlDict b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/controlDict index 4aa6dd1f74..a364cd0244 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/controlDict +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/controlDict @@ -60,14 +60,14 @@ functions0 writeControl writeTime; fields ( - U.airI + U.air1 { mean on; prime2Mean off; base time; } - U.airII + U.air2 { mean on; prime2Mean off; @@ -81,14 +81,14 @@ functions0 base time; } - alpha.airI + alpha.air1 { mean on; prime2Mean off; base time; } - alpha.airII + alpha.air2 { mean on; prime2Mean off; diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/setFieldsDict b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/setFieldsDict index 6bd134c535..6c8e902e31 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/setFieldsDict +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/system/setFieldsDict @@ -17,8 +17,8 @@ FoamFile defaultFieldValues ( - volScalarFieldValue alpha.airI 0 - volScalarFieldValue alpha.airII 1 + volScalarFieldValue alpha.air1 0 + volScalarFieldValue alpha.air2 1 volScalarFieldValue alpha.water 0 ); @@ -29,7 +29,7 @@ regions box (0 0 -0.1) (0.15 0.701 0.1); fieldValues ( - volScalarFieldValue alpha.airII 0 + volScalarFieldValue alpha.air2 0 volScalarFieldValue alpha.water 1 ); } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air1 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air1 index 03b1f9a28e..ec92afac21 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/T.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air1 @@ -10,7 +10,7 @@ FoamFile version 2.0; format ascii; class volScalarField; - object T.airII; + object T.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -27,7 +27,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airII; + phi phi.air1; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air2 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air2 index 03b1f9a28e..d1cafcb7e0 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/T.air2 @@ -10,7 +10,7 @@ FoamFile version 2.0; format ascii; class volScalarField; - object T.airII; + object T.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -27,7 +27,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airII; + phi phi.air2; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air1 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air1 index 8c2d21224c..ced8937c0a 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air1 @@ -10,7 +10,7 @@ FoamFile version 2.0; format binary; class volVectorField; - object U.airII; + object U.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -28,7 +28,7 @@ boundaryField outlet { type pressureInletOutletVelocity; - phi phi.airII; + phi phi.air1; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air2 similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air2 index 8c2d21224c..a83d1299ea 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/U.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/U.air2 @@ -10,7 +10,7 @@ FoamFile version 2.0; format binary; class volVectorField; - object U.airII; + object U.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -28,7 +28,7 @@ boundaryField outlet { type pressureInletOutletVelocity; - phi phi.airII; + phi phi.air2; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airI.orig b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig similarity index 95% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airI.orig rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig index 203cc22fc5..4a30789394 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airI.orig +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object alpha.airI; + object alpha.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airI; + phi phi.air1; inletValue uniform 0; value uniform 0; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airII.orig b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air2.orig similarity index 94% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airII.orig rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air2.orig index 019c483a0e..5c8b4c69cc 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.airII.orig +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air2.orig @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object alpha.airII; + object alpha.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airII; + phi phi.air2; inletValue uniform 1; value uniform 1; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airI.bubbles b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles similarity index 94% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airI.bubbles rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles index 9311ec38fc..c3a6d0ec90 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airI.bubbles +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object f.airI.bubbles; + object f.air1.bubbles; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airI; + phi phi.air1; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airII.bubbles b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles similarity index 94% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airII.bubbles rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles index 643e75c2b1..5494f7f5f3 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.airII.bubbles +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles @@ -11,7 +11,7 @@ FoamFile format ascii; class volScalarField; location "0"; - object f.airII.bubbles; + object f.air2.bubbles; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -29,7 +29,7 @@ boundaryField outlet { type inletOutlet; - phi phi.airII; + phi phi.air2; inletValue $internalField; value $internalField; } diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/p_rgh b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/p_rgh index 8d1b65707c..5dcbfa90da 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/p_rgh +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/p_rgh @@ -29,9 +29,9 @@ boundaryField { type prghTotalPressure; p0 $internalField; - U U.airII; - phi phi.airII; - rho thermo:rho.airII; + U U.air2; + phi phi.air2; + rho thermo:rho.air2; value $internalField; } walls diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties index 96469c3f75..35f0c3544c 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties @@ -17,11 +17,11 @@ FoamFile type populationBalanceMultiphaseSystem; -phases (airI airII water); +phases (air1 air2 water); populationBalances (bubbles); -airI +air1 { type purePhaseModel; diameterModel velocityGroup; @@ -51,7 +51,7 @@ airI residualAlpha 1e-6; } -airII +air2 { type purePhaseModel; diameterModel velocityGroup; @@ -133,10 +133,10 @@ blending default { type linear; - minFullyContinuousAlpha.airI 0.7; - minPartlyContinuousAlpha.airI 0.5; - minFullyContinuousAlpha.airII 0.7; - minPartlyContinuousAlpha.airII 0.5; + minFullyContinuousAlpha.air1 0.7; + minPartlyContinuousAlpha.air1 0.5; + minFullyContinuousAlpha.air2 0.7; + minPartlyContinuousAlpha.air2 0.5; minFullyContinuousAlpha.water 0.7; minPartlyContinuousAlpha.water 0.5; } @@ -144,13 +144,13 @@ blending surfaceTension ( - (airI and water) + (air1 and water) { type constant; sigma 0.07; } - (airII and water) + (air2 and water) { type constant; sigma 0.07; @@ -162,25 +162,25 @@ interfaceCompression aspectRatio ( - (airI in water) + (air1 in water) { type constant; E0 1.0; } - (water in airI) + (water in air1) { type constant; E0 1.0; } - (airII in water) + (air2 in water) { type constant; E0 1.0; } - (water in airII) + (water in air2) { type constant; E0 1.0; @@ -189,7 +189,7 @@ aspectRatio drag ( - (airI in water) + (air1 in water) { type SchillerNaumann; residualRe 1e-3; @@ -199,7 +199,7 @@ drag } } - (water in airI) + (water in air1) { type SchillerNaumann; residualRe 1e-3; @@ -209,7 +209,7 @@ drag } } - (airI and water) + (air1 and water) { type segregated; m 0.5; @@ -220,7 +220,7 @@ drag } } - (airII in water) + (air2 in water) { type SchillerNaumann; residualRe 1e-3; @@ -230,7 +230,7 @@ drag } } - (water in airII) + (water in air2) { type SchillerNaumann; residualRe 1e-3; @@ -240,7 +240,7 @@ drag } } - (airII and water) + (air2 and water) { type segregated; m 0.5; @@ -251,7 +251,7 @@ drag } } - (airI in airII) + (air1 in air2) { type SchillerNaumann; residualRe 1e-3; @@ -261,7 +261,7 @@ drag } } - (airII in airI) + (air2 in air1) { type SchillerNaumann; @@ -272,7 +272,7 @@ drag } } - (airI and airII) + (air1 and air2) { type segregated; m 0.5; @@ -286,91 +286,91 @@ drag virtualMass ( - (airI in water) + (air1 in water) { type constantCoefficient; Cvm 0.5; } - (water in airI) + (water in air1) { type constantCoefficient; Cvm 0.5; } - (airII in water) + (air2 in water) { type constantCoefficient; Cvm 0.5; } - (water in airII) + (water in air2) { type constantCoefficient; Cvm 0.5; } - (airI in airII) + (air1 in air2) { type constantCoefficient; Cvm 0.5; } - (airI in airII) + (air1 in air2) { type constantCoefficient; Cvm 0.5; } ); -heatTransfer.airI +heatTransfer.air1 ( - (airI in water) + (air1 in water) { type spherical; residualAlpha 1e-3; } - (water in airI) + (water in air1) { type RanzMarshall; residualAlpha 1e-3; } - (airI in airII) + (air1 in air2) { type spherical; residualAlpha 1e-3; } - (airII in airI) + (air2 in air1) { type spherical; residualAlpha 1e-3; } ); -heatTransfer.airII +heatTransfer.air2 ( - (airII in water) + (air2 in water) { type spherical; residualAlpha 1e-3; } - (water in airII) + (water in air2) { type RanzMarshall; residualAlpha 1e-3; } - (airI in airII) + (air1 in air2) { type spherical; residualAlpha 1e-3; } - (airII in airI) + (air2 in air1) { type spherical; residualAlpha 1e-3; @@ -379,35 +379,35 @@ heatTransfer.airII heatTransfer.water ( - (airI in water) + (air1 in water) { type RanzMarshall; residualAlpha 1e-3; } - (water in airI) + (water in air1) { type spherical; residualAlpha 1e-3; } - (airII in water) + (air2 in water) { type RanzMarshall; residualAlpha 1e-3; } - (water in airII) + (water in air2) { type spherical; residualAlpha 1e-3; } ); -massTransfer.airI +massTransfer.air1 (); -massTransfer.airII +massTransfer.air2 (); massTransfer.water diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air1 similarity index 96% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air1 index a190c7d7a9..7426132d89 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air1 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object thermophysicalProperties.airII; + object thermophysicalProperties.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airII b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air2 similarity index 96% rename from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airII rename to tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air2 index a190c7d7a9..62bff79a0a 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airII +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/thermophysicalProperties.air2 @@ -11,7 +11,7 @@ FoamFile format ascii; class dictionary; location "constant"; - object thermophysicalProperties.airII; + object thermophysicalProperties.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air1 b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air1 new file mode 100644 index 0000000000..fda3c5b86f --- /dev/null +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air1 @@ -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.air1; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air2 b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air2 new file mode 100644 index 0000000000..f6bb961594 --- /dev/null +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air2 @@ -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.air2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType laminar; + +// ************************************************************************* // diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/controlDict b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/controlDict index 1dcc62d121..3bda239ba6 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/controlDict +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/controlDict @@ -60,14 +60,14 @@ functions0 writeControl writeTime; fields ( - U.airI + U.air1 { mean on; prime2Mean off; base time; } - U.airII + U.air2 { mean on; prime2Mean off; @@ -81,14 +81,14 @@ functions0 base time; } - alpha.airI + alpha.air1 { mean on; prime2Mean off; base time; } - alpha.airII + alpha.air2 { mean on; prime2Mean off; diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/setFieldsDict b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/setFieldsDict index 6bd134c535..6c8e902e31 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/setFieldsDict +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/setFieldsDict @@ -17,8 +17,8 @@ FoamFile defaultFieldValues ( - volScalarFieldValue alpha.airI 0 - volScalarFieldValue alpha.airII 1 + volScalarFieldValue alpha.air1 0 + volScalarFieldValue alpha.air2 1 volScalarFieldValue alpha.water 0 ); @@ -29,7 +29,7 @@ regions box (0 0 -0.1) (0.15 0.701 0.1); fieldValues ( - volScalarFieldValue alpha.airII 0 + volScalarFieldValue alpha.air2 0 volScalarFieldValue alpha.water 1 ); }