diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.air deleted file mode 100644 index 0a372e6175..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.air +++ /dev/null @@ -1,45 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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 -{ - walls - { - type zeroGradient; - } - outlet - { - type inletOutlet; - phi phi.air; - inletValue $internalField; - value $internalField; - } - inlet - { - type fixedValue; - value $internalField; - } - frontAndBackPlanes - { - type empty; - } -} - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.water deleted file mode 100644 index 26f3e79683..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/T.water +++ /dev/null @@ -1,45 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - object T.water; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 0 0 1 0 0 0]; - -internalField uniform 350; - -boundaryField -{ - walls - { - type zeroGradient; - } - outlet - { - type inletOutlet; - phi phi.water; - inletValue $internalField; - value $internalField; - } - inlet - { - type fixedValue; - value $internalField; - } - frontAndBackPlanes - { - type empty; - } -} - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/Theta b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/Theta deleted file mode 100644 index 377979b988..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/Theta +++ /dev/null @@ -1,47 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - object Theta; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 2 -2 0 0 0 0]; - -internalField uniform 0.0; - -boundaryField -{ - inlet - { - type fixedValue; - value uniform 1.0e-7; - } - - outlet - { - type inletOutlet; - inletValue uniform 1.0e-7; - value uniform 1.0e-7; - } - - walls - { - type zeroGradient; - } - - defaultFaces - { - type empty; - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.air deleted file mode 100644 index 6738552ad8..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.air +++ /dev/null @@ -1,41 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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.1 0); - -boundaryField -{ - inlet - { - type fixedValue; - value $internalField; - } - outlet - { - type pressureInletOutletVelocity; - phi phi.air; - value $internalField; - } - walls - { - type fixedValue; - value uniform (0 0 0); - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.water deleted file mode 100644 index edafc657f1..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/U.water +++ /dev/null @@ -1,41 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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 -{ - inlet - { - type fixedValue; - value $internalField; - } - outlet - { - type pressureInletOutletVelocity; - phi phi.water; - value $internalField; - } - walls - { - type fixedValue; - value uniform (0 0 0); - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.air.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.air.orig deleted file mode 100644 index ddbf5788cc..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.air.orig +++ /dev/null @@ -1,42 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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; - -boundaryField -{ - inlet - { - type fixedValue; - value uniform 0.5; - } - outlet - { - type inletOutlet; - phi phi.air; - inletValue uniform 1; - value uniform 1; - } - walls - { - type zeroGradient; - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.air deleted file mode 100644 index 6dcf851a8e..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.air +++ /dev/null @@ -1,49 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - object k.air; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 2 -2 0 0 0 0]; - -internalField uniform 3.75e-5; - -boundaryField -{ - inlet - { - type fixedValue; - value $internalField; - } - - outlet - { - type inletOutlet; - phi phi.air; - inletValue $internalField; - value $internalField; - } - - walls - { - type kqRWallFunction; - value $internalField; - } - - defaultFaces - { - type empty; - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p deleted file mode 100644 index 334c4205ed..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p +++ /dev/null @@ -1,40 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class volScalarField; - object p; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [1 -1 -2 0 0 0 0]; - -internalField uniform 1e5; - -boundaryField -{ - inlet - { - type calculated; - value $internalField; - } - outlet - { - type calculated; - value $internalField; - } - walls - { - type calculated; - value $internalField; - } -} - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p_rgh b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p_rgh deleted file mode 100644 index f3becb1f99..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/p_rgh +++ /dev/null @@ -1,41 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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 -{ - inlet - { - type fixedFluxPressure; - value $internalField; - } - outlet - { - type prghPressure; - p $internalField; - value $internalField; - } - walls - { - type fixedFluxPressure; - value $internalField; - } -} - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/Allrun b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/Allrun deleted file mode 100755 index c9bcf24f33..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/Allrun +++ /dev/null @@ -1,12 +0,0 @@ -#!/bin/sh -cd ${0%/*} || exit 1 # Run from this directory - -# Source tutorial run functions -. $WM_PROJECT_DIR/bin/tools/RunFunctions - -runApplication blockMesh -runApplication topoSet -runApplication setFields -runApplication $(getApplication) - -#------------------------------------------------------------------------------ diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/g b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/g deleted file mode 100644 index e75aeb5349..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/g +++ /dev/null @@ -1,22 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class uniformDimensionedVectorField; - location "constant"; - object g; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 1 -2 0 0 0 0]; -value (0 -9.81 0); - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.air deleted file mode 100644 index bfc76d624d..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.air +++ /dev/null @@ -1,34 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "constant"; - object momentumTransport.air; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -simulationType RAS; - -RAS -{ - model mixtureKEpsilon; // continuousGasKEpsilon; - - turbulence on; - printCoeffs on; - - mixtureKEpsilonCoeffs - { - Cp 1; - C3 1; - } -} - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties deleted file mode 100644 index c715eeaebc..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties +++ /dev/null @@ -1,233 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "constant"; - object phaseProperties; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -type populationBalanceMultiphaseSystem; - -phases (air water); - -populationBalances (bubbles); - -air -{ - type purePhaseModel; - diameterModel velocityGroup; - velocityGroupCoeffs - { - populationBalance bubbles; - - shapeModel spherical; - - sizeGroups - ( - f0{dSph 1.00e-3; value 0;} - f1{dSph 1.08e-3; value 0;} - f2{dSph 1.16e-3; value 0;} - f3{dSph 1.25e-3; value 0;} - f4{dSph 1.36e-3; value 0;} - f5{dSph 1.46e-3; value 0;} - f6{dSph 1.58e-3; value 0;} - f7{dSph 1.71e-3; value 0;} - f8{dSph 1.85e-3; value 0;} - f9{dSph 2.00e-3; value 0;} - f10{dSph 2.16e-3; value 0;} - f11{dSph 2.33e-3; value 0;} - f12{dSph 2.51e-3; value 0;} - f13{dSph 2.72e-3; value 0.25;} - f14{dSph 2.93e-3; value 0.5;} - f15{dSph 3.17e-3; value 0.25 ;} - f16{dSph 3.42e-3; value 0;} - f17{dSph 3.70e-3; value 0;} - f18{dSph 4.00e-3; value 0;} - f19{dSph 4.32e-3; value 0;} - f20{dSph 4.66e-3; value 0;} - f21{dSph 5.03e-3; value 0;} - ); - } - - residualAlpha 1e-6; -} - -water -{ - type purePhaseModel; - diameterModel constant; - constantCoeffs - { - d 1e-4; - } - - residualAlpha 1e-6; -} - -populationBalanceCoeffs -{ - bubbles - { - continuousPhase water; - - coalescenceModels - ( - LehrMilliesMewes{} - ); - - binaryBreakupModels - ( - LehrMilliesMewes{} - ); - - breakupModels - (); - - driftModels - ( - densityChange{} - ); - - nucleationModels - (); - } -} - -blending -{ - default - { - type linear; - minFullyContinuousAlpha.air 0.7; - minPartlyContinuousAlpha.air 0.3; - minFullyContinuousAlpha.water 0.7; - minPartlyContinuousAlpha.water 0.3; - } - - drag - { - type linear; - minFullyContinuousAlpha.air 0.7; - minPartlyContinuousAlpha.air 0.5; - minFullyContinuousAlpha.water 0.7; - minPartlyContinuousAlpha.water 0.5; - } -} - -surfaceTension -( - (air and water) - { - type constant; - sigma 0.07; - } -); - -aspectRatio -( - (air in water) - { - type constant; - E0 1.0; - } - - (water in air) - { - type constant; - E0 1.0; - } -); - -drag -( - (air in water) - { - type SchillerNaumann; - residualRe 1e-3; - swarmCorrection - { - type none; - } - } - - (water in air) - { - type SchillerNaumann; - residualRe 1e-3; - swarmCorrection - { - type none; - } - } - - (air and water) - { - type segregated; - m 0.5; - n 8; - swarmCorrection - { - type none; - } - } -); - -virtualMass -( - (air in water) - { - type constantCoefficient; - Cvm 0.5; - } - - (water in air) - { - type constantCoefficient; - Cvm 0.5; - } -); - -heatTransfer -( - (air in water) - { - type RanzMarshall; - residualAlpha 1e-4; - } - - (water in air) - { - type RanzMarshall; - residualAlpha 1e-4; - } -); - -phaseTransfer -(); - -lift -(); - -wallLubrication -(); - -turbulentDispersion -(); - -interfaceCompression -(); - -// Minimum allowable pressure -pMin 10000; - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.air deleted file mode 100644 index aca85e96f4..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.air +++ /dev/null @@ -1,48 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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 perfectGas; - specie specie; - energy sensibleInternalEnergy; -} - -mixture -{ - specie - { - molWeight 28.9; - } - thermodynamics - { - Cp 1007; - Hf 0; - } - transport - { - mu 1.84e-05; - Pr 0.7; - } -} - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.water deleted file mode 100644 index f745ebe9c5..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/thermophysicalProperties.water +++ /dev/null @@ -1,52 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "constant"; - object thermophysicalProperties.water; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -thermoType -{ - type heRhoThermo; - mixture pureMixture; - transport const; - thermo eConst; - equationOfState rPolynomial; - specie specie; - energy sensibleInternalEnergy; -} - -mixture -{ - specie - { - molWeight 18; - } - equationOfState - { - C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16); - } - thermodynamics - { - Cv 4195; - Hf 0; - } - transport - { - mu 3.645e-4; - Pr 2.289; - } -} - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/blockMeshDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/blockMeshDict deleted file mode 100644 index 2aaa8f88ec..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/blockMeshDict +++ /dev/null @@ -1,61 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - object blockMeshDict; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -convertToMeters 1; - -vertices -( - (0 0 0) - (0.15 0 0) - (0.15 1 0) - (0 1 0) - (0 0 0.1) - (0.15 0 0.1) - (0.15 1 0.1) - (0 1 0.1) -); - -blocks -( - hex (0 1 2 3 4 5 6 7) (25 75 1) simpleGrading (1 1 1) -); - -edges -( -); - -patches -( - patch inlet - ( - (1 5 4 0) - ) - patch outlet - ( - (3 7 6 2) - ) - wall walls - ( - (0 4 7 3) - (2 6 5 1) - ) -); - -mergePatchPairs -( -); - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/controlDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/controlDict deleted file mode 100644 index ff6d2d0219..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/controlDict +++ /dev/null @@ -1,76 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "system"; - object controlDict; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -application multiphaseEulerFoam; - -startFrom startTime; - -startTime 0; - -stopAt endTime; - -endTime 100; - -deltaT 0.005; - -writeControl runTime; - -writeInterval 1; - -purgeWrite 0; - -writeFormat ascii; - -writePrecision 6; - -writeCompression off; - -timeFormat general; - -timePrecision 6; - -runTimeModifiable yes; - -adjustTimeStep no; - -maxCo 0.5; - -maxDeltaT 1; - -functions -{ - #includeFunc fieldAverage(U.air, U.water, alpha.air, p) - - numberDensity - { - type sizeDistribution; - functionObjectLibs ("libmultiphaseEulerFoamFunctionObjects.so"); - - writeControl outputTime; - writeInterval 1; - - populationBalance bubbles; - regionType cellZone; - name zone; - functionType number; - coordinateType diameter; - densityFunction yes; - } -} - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSchemes b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSchemes deleted file mode 100644 index c16f288a39..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSchemes +++ /dev/null @@ -1,67 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ 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.air) Gauss vanLeer; - div(phi,alpha.water) Gauss vanLeer; - div(phir,alpha.water,alpha.air) Gauss vanLeer; - div(phir,alpha.air,alpha.water) Gauss vanLeer; - - "div\(alphaRhoPhi.*,U.*\)" Gauss limitedLinearV 1; - "div\(phi.*,U.*\)" Gauss limitedLinearV 1; - - "div\(alphaRhoPhi.*,(h|e).*\)" Gauss limitedLinear 1; - "div\(alphaPhi.*,f.*\)" Gauss limitedLinear 1; - "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; - "div\(alphaPhi.*,p\)" Gauss limitedLinear 1; - - "div\(alphaRhoPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1; - "div\(phim,(k|epsilon)m\)" Gauss limitedLinear 1; - - "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; -} - -laplacianSchemes -{ - default Gauss linear uncorrected; -} - -interpolationSchemes -{ - default linear; -} - -snGradSchemes -{ - default uncorrected; -} - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSolution b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSolution deleted file mode 100644 index d91095e200..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/fvSolution +++ /dev/null @@ -1,101 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "system"; - object fvSolution; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -solvers -{ - "alpha.*" - { - nAlphaCorr 1; - nAlphaSubCycles 2; - } - - bubbles - { - nCorr 1; - tolerance 1e-4; - scale true; - sourceUpdateInterval 1; - solveOnFinalIterOnly true; - } - - p_rgh - { - solver GAMG; - smoother DIC; - tolerance 1e-8; - relTol 0; - } - - p_rghFinal - { - $p_rgh; - relTol 0; - } - - "U.*" - { - solver smoothSolver; - smoother symGaussSeidel; - tolerance 1e-7; - relTol 0; - minIter 1; - } - - "e.*" - { - solver smoothSolver; - smoother symGaussSeidel; - tolerance 1e-7; - relTol 0; - minIter 1; - } - - "(k|epsilon|Theta).*" - { - solver smoothSolver; - smoother symGaussSeidel; - tolerance 1e-7; - relTol 0; - minIter 1; - } - - "f.*" - { - solver PBiCGStab; - preconditioner DILU; - tolerance 1e-6; - relTol 0; - } -} - -PIMPLE -{ - nOuterCorrectors 3; - nCorrectors 1; - nNonOrthogonalCorrectors 0; -} - -relaxationFactors -{ - equations - { - ".*" 1; - } -} - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/setFieldsDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/setFieldsDict deleted file mode 100644 index 572fdd0f28..0000000000 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/setFieldsDict +++ /dev/null @@ -1,38 +0,0 @@ -/*--------------------------------*- C++ -*----------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Version: dev - \\/ M anipulation | -\*---------------------------------------------------------------------------*/ -FoamFile -{ - version 2.0; - format ascii; - class dictionary; - location "system"; - object setFieldsDict; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -defaultFieldValues -( - volScalarFieldValue alpha.air 1 - volScalarFieldValue alpha.water 0 -); - -regions -( - boxToCell - { - box (0 0 -0.1) (0.15 0.701 0.1); - fieldValues - ( - volScalarFieldValue alpha.air 0 - volScalarFieldValue alpha.water 1 - ); - } -); - - -// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/f.air.bubbles b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air1 similarity index 76% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/f.air.bubbles rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air1 index eabe6c247d..6413c104ff 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/f.air.bubbles +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air1 @@ -10,29 +10,29 @@ FoamFile version 2.0; format ascii; class volScalarField; - location "0"; - object f.air.bubbles; + object T.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 0 0 0 0 0 0]; +dimensions [0 0 0 1 0 0 0]; -internalField uniform 1; +internalField uniform 300; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; value $internalField; } + outlet { - type inletOutlet; - phi phi.air; - inletValue $internalField; - value $internalField; + type zeroGradient; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.water.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air2 similarity index 75% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.water.orig rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air2 index c7c9938b2f..e9cb49ffda 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alpha.water.orig +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.air2 @@ -10,28 +10,29 @@ FoamFile version 2.0; format ascii; class volScalarField; - object alpha.water; + object T.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 0 0 0 0 0 0]; +dimensions [0 0 0 1 0 0 0]; -internalField uniform 1; +internalField uniform 300; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; - value uniform 0.5; + value $internalField; } + outlet { - type inletOutlet; - phi phi.water; - inletValue uniform 0; - value uniform 0; + type zeroGradient; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.water.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.water similarity index 75% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.water.orig rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.water index c7c9938b2f..0d31cca9f1 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.water.orig +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/T.water @@ -10,28 +10,29 @@ FoamFile version 2.0; format ascii; class volScalarField; - object alpha.water; + object T.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 0 0 0 0 0 0]; +dimensions [0 0 0 1 0 0 0]; -internalField uniform 1; +internalField uniform 300; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; - value uniform 0.5; + value $internalField; } + outlet { - type inletOutlet; - phi phi.water; - inletValue uniform 0; - value uniform 0; + type zeroGradient; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilonm b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air1 similarity index 63% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilonm rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air1 index fc18e67dbc..5d2855187e 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilonm +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air1 @@ -9,40 +9,39 @@ FoamFile { version 2.0; format ascii; - class volScalarField; - object epsilonm; + class volVectorField; + object U.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 2 -3 0 0 0 0]; +dimensions [0 1 -1 0 0 0 0]; -internalField uniform 1.5e-4; +internalField uniform (0 0 0); boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type fixedValue; + type mapped; + field U.air1; + setAverage no; + average (0 0 0); + interpolationScheme cell; value $internalField; } outlet { - type inletOutlet; - phi phim; - inletValue $internalField; + type pressureInletOutletVelocity; + phi phi.air1; value $internalField; } walls { - type zeroGradient; - value $internalField; - } - - defaultFaces - { - type empty; + type slip; } } diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air2 similarity index 63% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.air rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air2 index 9eb8e5eb39..8a2fea9f55 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.air +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.air2 @@ -9,40 +9,39 @@ FoamFile { version 2.0; format ascii; - class volScalarField; - object epsilon.air; + class volVectorField; + object U.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 2 -3 0 0 0 0]; +dimensions [0 1 -1 0 0 0 0]; -internalField uniform 1.5e-4; +internalField uniform (0 0 0); boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type fixedValue; + type mapped; + field U.air2; + setAverage no; + average (0 0 0); + interpolationScheme cell; value $internalField; } outlet { - type inletOutlet; - phi phi.air; - inletValue $internalField; + type pressureInletOutletVelocity; + phi phi.air2; value $internalField; } walls { - type epsilonWallFunction; - value $internalField; - } - - defaultFaces - { - type empty; + type slip; } } diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.water similarity index 61% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.water index b2c0cd7d21..be155f4703 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/epsilon.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/U.water @@ -7,43 +7,41 @@ \*---------------------------------------------------------------------------*/ FoamFile { - version 2.0; - format ascii; - class volScalarField; - object epsilon.water; + version 2; + format ascii; + class volVectorField; + object U.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 2 -3 0 0 0 0]; +dimensions [0 1 -1 0 0 0 0]; -internalField uniform 1.5e-4; +internalField uniform (0 0 1.6); boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type fixedValue; + type mapped; + field U.water; + setAverage yes; + average (0 0 1.6); + interpolationScheme cell; value $internalField; } - outlet { - type inletOutlet; + type pressureInletOutletVelocity; phi phi.water; - inletValue $internalField; value $internalField; } - walls { - type epsilonWallFunction; - value $internalField; - } - - defaultFaces - { - type empty; + type noSlip; } } + // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air1 similarity index 92% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air1 index 76ec8c113a..d9ea52e2e3 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/alpha.air1.orig +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air1 @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class volScalarField; - location "0"; object alpha.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -21,11 +20,14 @@ internalField uniform 0; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; - value uniform 0.25; + value uniform 0; } + outlet { type inletOutlet; @@ -33,6 +35,7 @@ boundaryField inletValue uniform 0; value uniform 0; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air2 similarity index 86% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air2 index 76ec8c113a..c50c4a6ec4 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.air1.orig +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.air2 @@ -10,8 +10,7 @@ FoamFile version 2.0; format ascii; class volScalarField; - location "0"; - object alpha.air1; + object alpha.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -21,18 +20,22 @@ internalField uniform 0; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; - value uniform 0.25; + value uniform 0; } + outlet { type inletOutlet; - phi phi.air1; + phi phi.air2; inletValue uniform 0; value uniform 0; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.water.orig b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.water similarity index 86% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.water.orig rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.water index c7c9938b2f..eff1a2d720 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/alpha.water.orig +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alpha.water @@ -20,18 +20,22 @@ internalField uniform 1; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; - value uniform 0.5; + value uniform 1; } + outlet { type inletOutlet; phi phi.water; - inletValue uniform 0; - value uniform 0; + inletValue uniform 1; + value uniform 1; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alphat.water similarity index 94% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alphat.water index d71c6f334c..b067b6a740 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/alphat.water @@ -20,6 +20,8 @@ internalField uniform 0; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type calculated; @@ -38,11 +40,6 @@ boundaryField Prt 0.85; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air1.bubbles similarity index 82% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air1.bubbles index 2164bcd380..456fb5900a 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air1.bubbles +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air1.bubbles @@ -10,29 +10,29 @@ 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; +internalField uniform 1.0; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; value $internalField; } + outlet { - type inletOutlet; - phi phi.air1; - inletValue $internalField; - value $internalField; + type zeroGradient; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air2.bubbles similarity index 82% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air2.bubbles index 7b81d4e221..43b2db8362 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/f.air2.bubbles +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/f.air2.bubbles @@ -10,29 +10,29 @@ 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; +internalField uniform 1.0; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type fixedValue; value $internalField; } + outlet { - type inletOutlet; - phi phi.air2; - inletValue $internalField; - value $internalField; + type zeroGradient; } + walls { type zeroGradient; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/k.water similarity index 76% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/k.water index 0c7cc3d01b..e6e6820354 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/k.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/k.water @@ -16,22 +16,25 @@ FoamFile dimensions [0 2 -2 0 0 0 0]; -internalField uniform 3.75e-5; +internalField uniform 0.01; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type fixedValue; + type mapped; + field k.water; + setAverage no; + average 0; + interpolationScheme cell; value $internalField; } outlet { - type inletOutlet; - phi phi.water; - inletValue $internalField; - value $internalField; + type zeroGradient; } walls @@ -39,11 +42,6 @@ boundaryField type kqRWallFunction; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/nut.water similarity index 91% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/nut.water index 4f22baf76e..4dc1afc5a2 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/nut.water @@ -16,10 +16,12 @@ FoamFile dimensions [0 2 -1 0 0 0 0]; -internalField uniform 1e-8; +internalField uniform 1e-4; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type calculated; @@ -37,11 +39,6 @@ boundaryField type nutkWallFunction; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/km b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/omega.water similarity index 71% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/km rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/omega.water index 1db01cbce9..bd3398476d 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/km +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/omega.water @@ -10,40 +10,38 @@ FoamFile version 2.0; format ascii; class volScalarField; - object km; + object omega.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 2 -2 0 0 0 0]; +dimensions [0 0 -1 0 0 0 0]; -internalField uniform 3.75e-5; +internalField uniform 100; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type fixedValue; + type mapped; + field omega.water; + setAverage no; + average 0; + interpolationScheme cell; value $internalField; } outlet { - type inletOutlet; - phi phim; - inletValue $internalField; - value $internalField; + type zeroGradient; } walls { - type zeroGradient; + type omegaWallFunction; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p similarity index 82% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.air rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p index b2bbd25a14..73cced3063 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/nut.air +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p @@ -10,16 +10,18 @@ FoamFile version 2.0; format ascii; class volScalarField; - object nut.air; + object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [0 2 -1 0 0 0 0]; +dimensions [1 -1 -2 0 0 0 0]; -internalField uniform 1e-8; +internalField uniform 1e5; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { type calculated; @@ -34,14 +36,9 @@ boundaryField walls { - type nutkWallFunction; + type calculated; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.air b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p_rgh similarity index 73% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.air rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p_rgh index 91661525ab..6143b21dca 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/0/alphat.air +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/0/p_rgh @@ -10,39 +10,36 @@ FoamFile version 2.0; format ascii; class volScalarField; - object alphat.air; + object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -dimensions [1 -1 -1 0 0 0 0]; +dimensions [1 -1 -2 0 0 0 0]; -internalField uniform 0; +internalField uniform 1e5; boundaryField { + #includeEtc "caseDicts/setConstraintTypes" + inlet { - type calculated; + type fixedFluxPressure; value $internalField; } outlet { - type calculated; + type prghPressure; + p $internalField; value $internalField; } walls { - type compressible::alphatWallFunction; - Prt 0.85; + type fixedFluxPressure; value $internalField; } - - defaultFaces - { - type empty; - } } // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allclean b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allclean new file mode 100755 index 0000000000..fa2b9901cf --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allclean @@ -0,0 +1,13 @@ +#!/bin/sh + +# Run from this directory +cd "${0%/*}" || exit 1 + +# Source clean functions +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase && rm -f *.eps + +(cd 0 && rm -rf uniform d.* *Force.water) + +#------------------------------------------------------------------------------ diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allrun b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allrun new file mode 100755 index 0000000000..83a37ba4d5 --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/Allrun @@ -0,0 +1,18 @@ +#!/bin/sh + +# Run from this directory +cd "${0%/*}" || exit 1 + +# Source run functions +. "$WM_PROJECT_DIR/bin/tools/RunFunctions" + +# Run case +runApplication blockMesh +runApplication topoSet +runApplication decomposePar +runParallel $(getApplication) +runApplication reconstructPar + +( cd validation && ./createGraphs ) + +#------------------------------------------------------------------------------ diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/topoSetDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/fvOptions similarity index 52% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/topoSetDict rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/fvOptions index 9cb2253ee4..a32133c7e3 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/system/topoSetDict +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/fvOptions @@ -10,36 +10,35 @@ FoamFile version 2.0; format ascii; class dictionary; - object topoSetDict; + object fvOptions; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -actions -( +options +{ + massSource { - name zone; - type cellSet; - action new; - source boxToCell; - sourceInfo + type semiImplicitSource; + + selectionMode cellZone; + cellZone injection; + + volumeMode absolute; + + sources { - boxes - ( - (0 0.4 0) (0.15 0.5 0.1) - ); + thermo:rho.air1 + { + explicit 6e-07; + implicit 0; + } + + f3.air1.bubbles + { + value 1.0; + explicit #calc "$value*$../thermo:rho.air1/explicit"; + implicit 0; + } } } - - { - name zone; - type cellZoneSet; - action new; - source setToCellZone; - sourceInfo - { - set zone; - } - } -); - -// ************************************************************************* // +} diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/g b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/g similarity index 92% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/g rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/g index e75aeb5349..88fccf33b0 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/g +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/g @@ -10,13 +10,11 @@ FoamFile version 2.0; format ascii; class uniformDimensionedVectorField; - location "constant"; object g; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -2 0 0 0 0]; -value (0 -9.81 0); - +value (0 0 -9.81); // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air1 b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air1 similarity index 96% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air1 rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air1 index 3c3b2564e9..3bafec10ec 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air1 +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air1 @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; object momentumTransport.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.air2 b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air2 similarity index 92% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.air2 rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air2 index 217f7e40ef..3172370612 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.air2 +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.air2 @@ -10,11 +10,10 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; object momentumTransport.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -simulationType laminar; +#include "momentumTransport.air1" // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.water similarity index 90% rename from tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.water index 0dcd319014..bfc97b6399 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/momentumTransport.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/momentumTransport.water @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; object momentumTransport.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -19,7 +18,7 @@ simulationType RAS; RAS { - model mixtureKEpsilon; // LaheyKEpsilon; + model kOmegaSSTSato; turbulence on; printCoeffs on; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/phaseProperties b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/phaseProperties new file mode 100644 index 0000000000..d72d179c96 --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/phaseProperties @@ -0,0 +1,285 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: dev + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object phaseProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +type populationBalanceMultiphaseSystem; + +phases (air1 air2 water); + +populationBalances (bubbles); + +air1 +{ + type pureIsothermalPhaseModel; + + diameterModel velocityGroup; + + velocityGroupCoeffs + { + populationBalance bubbles; + + shapeModel spherical; + + sizeGroups + ( + f1 {dSph 1e-3; value 1.0;} + f2 {dSph 2e-3; value 0.0;} + f3 {dSph 3e-3; value 0.0;} + f4 {dSph 4e-3; value 0.0;} + f5 {dSph 5e-3; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +air2 +{ + type pureIsothermalPhaseModel; + + diameterModel velocityGroup; + + velocityGroupCoeffs + { + populationBalance bubbles; + + shapeModel spherical; + + sizeGroups + ( + f6 {dSph 6e-3; value 1.0;} + f7 {dSph 7e-3; value 0.0;} + f8 {dSph 8e-3; value 0.0;} + f9 {dSph 9e-3; value 0.0;} + f10{dSph 10e-3; value 0.0;} + f11{dSph 11e-3; value 0.0;} + f12{dSph 12e-3; value 0.0;} + ); + } + + residualAlpha 1e-6; +} + +water +{ + type pureIsothermalPhaseModel; + + diameterModel constant; + + constantCoeffs + { + d 1e-4; + } + + residualAlpha 1e-6; +} + +populationBalanceCoeffs +{ + bubbles + { + continuousPhase water; + + coalescenceModels + ( + LehrMilliesMewes{} + ); + + binaryBreakupModels + ( + LehrMilliesMewes{} + ); + + breakupModels + (); + + driftModels + ( + densityChange{} + ); + + nucleationModels + (); + } +} + +blending +{ + default + { + type none; + continuousPhase water; + } +} + +surfaceTension +( + (air1 and water) + { + type constant; + sigma 0.083; + } + + (air2 and water) + { + type constant; + sigma 0.083; + } +); + +interfaceCompression +(); + +aspectRatio +( + (air1 in water) + { + type Wellek; + } + + (air2 in water) + { + type Wellek; + } +); + +drag +( + (air1 in water) + { + type IshiiZuber; + + residualRe 1e-3; + swarmCorrection + { + type none; + } + } + + (air2 in water) + { + type IshiiZuber; + + residualRe 1e-3; + swarmCorrection + { + type none; + } + } +); + +virtualMass +( + (air1 in water) + { + type constantCoefficient; + Cvm 0.5; + } + + (air2 in water) + { + type constantCoefficient; + Cvm 0.5; + } +); + +heatTransfer +(); + +phaseTransfer +(); + +lift +( + (air1 in water) + { + type wallDamped; + + wallDamping + { + type cosine; + Cd 3.0; + } + + lift + { + type Tomiyama; + + swarmCorrection + { + type none; + } + } + } + + (air2 in water) + { + type wallDamped; + + wallDamping + { + type cosine; + Cd 3.0; + } + + lift + { + type Tomiyama; + + swarmCorrection + { + type none; + } + } + } +); + +wallLubrication +( + (air1 in water) + { + type Antal; + Cw1 -0.01; + Cw2 0.05; + } + + (air2 in water) + { + type Antal; + Cw1 -0.01; + Cw2 0.05; + } +); + +turbulentDispersion +( + (air1 in water) + { + type Burns; + sigma 0.9; + } + + (air2 in water) + { + type Burns; + sigma 0.9; + } +); + +// Minimum allowable pressure +pMin 10000; + + +// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air1 similarity index 97% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air1 index 8402edfc63..9d0a92e152 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1 +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air1 @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; object thermophysicalProperties.air1; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -44,5 +43,4 @@ mixture } } - // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air2 b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air2 similarity index 87% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air2 rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air2 index 217f7e40ef..e12a79c6e9 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/momentumTransport.air2 +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.air2 @@ -10,11 +10,10 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; - object momentumTransport.air2; + object thermophysicalProperties.air2; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -simulationType laminar; +#include "thermophysicalProperties.air1" // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.water similarity index 85% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.water index f745ebe9c5..d9162f3385 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/constant/thermophysicalProperties.water @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; object thermophysicalProperties.water; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -20,8 +19,8 @@ thermoType type heRhoThermo; mixture pureMixture; transport const; - thermo eConst; - equationOfState rPolynomial; + thermo hConst; + equationOfState rhoConst; specie specie; energy sensibleInternalEnergy; } @@ -34,11 +33,11 @@ mixture } equationOfState { - C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16); + rho 997; } thermodynamics { - Cv 4195; + Cp 4195; Hf 0; } transport @@ -48,5 +47,4 @@ mixture } } - // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/blockMeshDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/blockMeshDict new file mode 100644 index 0000000000..c6f0a023d9 --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/blockMeshDict @@ -0,0 +1,99 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: dev + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 1; + +vertices +( + (0.0 0.0 -0.1) + (0.025 -2.2e-4 -0.1) + (0.025 2.2e-4 -0.1) + (0.0 0.0 -0.1) + (0.0 0.0 0.9) + (0.025 -2.2e-4 0.9) + (0.025 2.2e-4 0.9) + (0.0 0.0 0.9) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (30 1 200) simpleGrading (0.375 1 1) +); + +edges +(); + +boundary +( + inlet + { + type mappedPatch; + offset (0 0 0.05); + sampleRegion region0; + sampleMode nearestCell; + samplePatch none; + faces + ( + (0 3 2 1) + ); + } + outlet + { + type patch; + faces + ( + (4 5 6 7) + ); + } + walls + { + type wall; + faces + ( + (1 2 6 5) + ); + } + front + { + type wedge; + faces + ( + (0 1 5 4) + ); + } + back + { + type wedge; + faces + ( + (2 3 7 6) + ); + } + axis + { + type empty; + faces + ( + (3 0 4 7) + ); + } +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/controlDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/controlDict new file mode 100644 index 0000000000..2118231bd2 --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/controlDict @@ -0,0 +1,97 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: dev + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application multiphaseEulerFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 4; + +deltaT 1e-4; + +writeControl adjustableRunTime; + +writeInterval 0.5; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep yes; + +maxCo 0.8; + +maxDeltaT 1; + +functions +{ + #includeFunc residuals(p_rgh, k.water, omega.water) + #includeFunc timeStep + #includeFunc writeObjects(d.air1) + #includeFunc writeObjects(d.air2) + #includeFunc phaseForces(phase=water) + #includeFunc singleGraph + + probabilityDensity.injection + { + type sizeDistribution; + functionObjectLibs ("libmultiphaseEulerFoamFunctionObjects.so"); + + writeControl outputTime; + writeInterval 1; + + populationBalance bubbles; + regionType cellZone; + name injection; + functionType volume; + coordinateType diameter; + densityFunction yes; + normalize yes; + } + + probabilityDensity.outlet + { + type sizeDistribution; + functionObjectLibs ("libmultiphaseEulerFoamFunctionObjects.so"); + + writeControl outputTime; + writeInterval 1; + + populationBalance bubbles; + regionType cellZone; + name outlet; + functionType volume; + coordinateType diameter; + densityFunction yes; + normalize yes; + } +} + +// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.water b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/decomposeParDict similarity index 80% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.water rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/decomposeParDict index 8f42ef3c8a..9e0c338807 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/momentumTransport.water +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/decomposeParDict @@ -10,11 +10,18 @@ FoamFile version 2.0; format ascii; class dictionary; - location "constant"; - object momentumTransport.water; + object decomposeParDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -simulationType laminar; +numberOfSubdomains 4; + +method simple; + +simpleCoeffs +{ + n (1 1 4); + delta 0.001; +} // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSchemes b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSchemes similarity index 62% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSchemes rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSchemes index da72928048..65cec97471 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSchemes +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSchemes @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -23,25 +22,22 @@ ddtSchemes gradSchemes { default Gauss linear; + limited cellLimited Gauss linear 1; } divSchemes { - default none; + default none; - "div\(phi,alpha.*)" Gauss vanLeer; + "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.*,(h|e).*\)" Gauss limitedLinear 1; + "div\(alphaRhoPhi.*,U.*\)" Gauss linearUpwind limited; + "div\(phi.*,U.*\)" Gauss linearUpwind limited; + "div\(alphaRhoPhi.*,(k|epsilon|omega).*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,f.*\)" Gauss limitedLinear 1; - "div\(alphaRhoPhi.*,K.*\)" Gauss limitedLinear 1; - "div\(alphaPhi.*,p\)" Gauss limitedLinear 1; - "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; + "div\(\(\(\(alpha.*\*thermo:rho.*\)*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; } laplacianSchemes @@ -54,10 +50,10 @@ interpolationSchemes default linear; } -snGradSchemes +wallDist { - default uncorrected; + method Poisson; + nRequired true; } - // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSolution b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSolution similarity index 77% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSolution rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSolution index 3958fede5c..3dfb8a953b 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/fvSolution +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/fvSolution @@ -10,7 +10,6 @@ FoamFile version 2.0; format ascii; class dictionary; - location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -28,11 +27,11 @@ solvers nCorr 1; tolerance 1e-4; scale true; - sourceUpdateInterval 1; solveOnFinalIterOnly true; + sourceUpdateInterval 1; } - p_rgh + "p_rgh.*" { solver GAMG; smoother DIC; @@ -40,26 +39,11 @@ solvers relTol 0; } - p_rghFinal - { - $p_rgh; - relTol 0; - } - - "U.*" + "(k|omega).*" { solver smoothSolver; smoother symGaussSeidel; - tolerance 1e-5; - relTol 0; - minIter 1; - } - - "e.*" - { - solver smoothSolver; - smoother symGaussSeidel; - tolerance 1e-8; + tolerance 1e-7; relTol 0; minIter 1; } @@ -71,23 +55,26 @@ solvers tolerance 1e-6; relTol 0; } + + yPsi + { + solver PCG; + preconditioner none; + tolerance 1e-10; + relTol 0; + } } PIMPLE { - nOuterCorrectors 3; - nCorrectors 1; + nOuterCorrectors 2; + nCorrectors 2; nNonOrthogonalCorrectors 0; - faceMomentum yes; + faceMomentum true; } relaxationFactors { - equations - { - ".*" 1; - } } - // ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/singleGraph b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/singleGraph new file mode 100644 index 0000000000..e6c1b6545b --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/singleGraph @@ -0,0 +1,41 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: dev + \\/ M anipulation | +------------------------------------------------------------------------------- +Description + Writes graph data for specified fields along a line, specified by start + and end points. + +\*---------------------------------------------------------------------------*/ + +start (0 0 0.89); +end (0.025 0 0.89); +fields +( + alpha.air1 + alpha.air2 + alpha.bubbles + liftForce.water + wallLubricationForce.water + turbulentDispersionForce.water +); + +// Sampling and I/O settings +#includeEtc "caseDicts/postProcessing/graphs/sampleDict.cfg" + +// Override settings here, e.g. +/* +setConfig +{ + type lineCell; + axis x; // y, z, xyz +} +*/ + +// Must be last entry +#includeEtc "caseDicts/postProcessing/graphs/graph.cfg" + +// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/topoSetDict b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/topoSetDict similarity index 64% rename from tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/topoSetDict rename to tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/topoSetDict index 9cb2253ee4..15047166c5 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/bubbleColumnPolydisperse/system/topoSetDict +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/system/topoSetDict @@ -17,7 +17,7 @@ FoamFile actions ( { - name zone; + name injection; type cellSet; action new; source boxToCell; @@ -25,19 +25,44 @@ actions { boxes ( - (0 0.4 0) (0.15 0.5 0.1) + (0 -1 0.0) (0.005 1 0.01) ); } } { - name zone; + name injection; type cellZoneSet; action new; source setToCellZone; sourceInfo { - set zone; + set injection; + } + } + + { + name outlet; + type cellSet; + action new; + source boxToCell; + sourceInfo + { + boxes + ( + (0 -1 0.89) (0.025 1 0.9) + ); + } + } + + { + name outlet; + type cellZoneSet; + action new; + source setToCellZone; + sourceInfo + { + set outlet; } } ); diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/validation/createGraphs b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/validation/createGraphs new file mode 100755 index 0000000000..9841f06f3f --- /dev/null +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/bubblePipe/validation/createGraphs @@ -0,0 +1,47 @@ +#!/bin/sh + +if ! which gnuplot > /dev/null 2>&1 +then + echo 'gnuplot not found - skipping graph creation' >&2 + exit 1 +fi + +gnuplot< 5.5mm' + + set output '../nonDragForces.eps' + + set ylabel 'F_{r} (N/m^3)' + + plot forces u (\$1*1000):(-\$2) w l lw 2 t 'Lift',\ + forces u (\$1*1000):(-\$5) w l lw 2 t 'Wall lubrication',\ + forces u (\$1*1000):(-\$8) w l lw 2 t 'Turbulent dispersion' +EOF + +gnuplot<