From b6b146058df39b6f139820889c72f70ded2f6ff3 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 22 Mar 2018 16:42:11 +0000 Subject: [PATCH] test: New directory for model tests and validation The initial set of cases in the test directory are aimed at testing the reactingEulerFoam populationBalance functionality. Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum Dresden - Rossendorf (HZDR) and VTT Technical Research Centre of Finland Ltd. Integrated with the "tutorials" functionality by CFD Direct Ltd. --- bin/tools/CleanFunctions | 4 +- test/Allclean | 1 + test/Allrun | 1 + .../populationBalanceModeling/README | 11 + .../binaryBreakup/0/T.air1 | 25 ++ .../binaryBreakup/0/T.air2 | 25 ++ .../binaryBreakup/0/T.air3 | 25 ++ .../binaryBreakup/0/T.water | 25 ++ .../binaryBreakup/0/U.air1 | 25 ++ .../binaryBreakup/0/U.air2 | 25 ++ .../binaryBreakup/0/U.air3 | 25 ++ .../binaryBreakup/0/U.water | 25 ++ .../binaryBreakup/0/alpha.air1 | 26 ++ .../binaryBreakup/0/alpha.air2 | 26 ++ .../binaryBreakup/0/alpha.air3 | 26 ++ .../binaryBreakup/0/alpha.water | 25 ++ .../binaryBreakup/0/f.air1.bubbles | 26 ++ .../binaryBreakup/0/f.air2.bubbles | 26 ++ .../binaryBreakup/0/f.air3.bubbles | 26 ++ .../binaryBreakup/0/p | 25 ++ .../binaryBreakup/0/p_rgh | 25 ++ .../binaryBreakup/Allrun | 39 +++ .../binaryBreakup/constant/g | 22 ++ .../binaryBreakup/constant/phaseProperties | 211 +++++++++++++++ .../constant/thermophysicalProperties.air1 | 49 ++++ .../constant/thermophysicalProperties.air2 | 49 ++++ .../constant/thermophysicalProperties.air3 | 49 ++++ .../constant/thermophysicalProperties.water | 53 ++++ .../constant/turbulenceProperties.air1 | 2 +- .../constant/turbulenceProperties.air2 | 2 +- .../constant/turbulenceProperties.air3 | 2 +- .../constant/turbulenceProperties.water | 2 +- .../binaryBreakup/plot | 27 ++ .../binaryBreakup/system/blockMeshDict | 48 ++++ .../binaryBreakup/system/controlDict | 72 +++++ .../binaryBreakup/system/fvSchemes | 54 ++++ .../binaryBreakup/system/fvSolution | 74 +++++ .../breakup/0/T.air1 | 25 ++ .../breakup/0/T.air2 | 25 ++ .../breakup/0/T.air3 | 25 ++ .../breakup/0/T.water | 25 ++ .../breakup/0/U.air1 | 25 ++ .../breakup/0/U.air2 | 25 ++ .../breakup/0/U.air3 | 25 ++ .../breakup/0/U.water | 25 ++ .../breakup/0/alpha.air1 | 26 ++ .../breakup/0/alpha.air2 | 26 ++ .../breakup/0/alpha.air3 | 26 ++ .../breakup/0/alpha.water | 25 ++ .../breakup/0/f.air1.bubbles | 26 ++ .../breakup/0/f.air2.bubbles | 26 ++ .../breakup/0/f.air3.bubbles | 26 ++ .../populationBalanceModeling/breakup/0/p | 25 ++ .../populationBalanceModeling/breakup/0/p_rgh | 25 ++ .../populationBalanceModeling/breakup/Allrun | 39 +++ .../breakup/constant/g | 22 ++ .../breakup/constant/phaseProperties | 212 +++++++++++++++ .../constant/thermophysicalProperties.air1 | 49 ++++ .../constant/thermophysicalProperties.air2 | 49 ++++ .../constant/thermophysicalProperties.air3 | 49 ++++ .../constant/thermophysicalProperties.water | 53 ++++ .../constant/turbulenceProperties.air1 | 20 ++ .../constant/turbulenceProperties.air2 | 20 ++ .../constant/turbulenceProperties.air3 | 20 ++ .../constant/turbulenceProperties.water | 20 ++ .../populationBalanceModeling/breakup/plot | 27 ++ .../breakup/system/blockMeshDict | 48 ++++ .../breakup/system/controlDict | 72 +++++ .../breakup/system/fvSchemes | 54 ++++ .../breakup/system/fvSolution | 74 +++++ .../coalescence/0/T.air1 | 25 ++ .../coalescence/0/T.air2 | 25 ++ .../coalescence/0/T.air3 | 25 ++ .../coalescence/0/T.water | 25 ++ .../coalescence/0/U.air1 | 25 ++ .../coalescence/0/U.air2 | 25 ++ .../coalescence/0/U.air3 | 25 ++ .../coalescence/0/U.water | 25 ++ .../coalescence/0/alpha.air1 | 26 ++ .../coalescence/0/alpha.air2 | 26 ++ .../coalescence/0/alpha.air3 | 26 ++ .../coalescence/0/alpha.water | 25 ++ .../coalescence/0/f.air1.bubbles | 26 ++ .../coalescence/0/f.air2.bubbles | 26 ++ .../coalescence/0/f.air3.bubbles | 26 ++ .../populationBalanceModeling/coalescence/0/p | 25 ++ .../coalescence/0/p_rgh | 25 ++ .../coalescence/Allrun | 39 +++ .../coalescence/constant/g | 22 ++ .../coalescence/constant/phaseProperties | 253 ++++++++++++++++++ .../constant/thermophysicalProperties.air1 | 49 ++++ .../constant/thermophysicalProperties.air2 | 49 ++++ .../constant/thermophysicalProperties.air3 | 49 ++++ .../constant/thermophysicalProperties.water | 53 ++++ .../constant/turbulenceProperties.air1 | 20 ++ .../constant/turbulenceProperties.air2 | 20 ++ .../constant/turbulenceProperties.air3 | 20 ++ .../constant/turbulenceProperties.water | 20 ++ .../coalescence/plot | 47 ++++ .../coalescence/system/blockMeshDict | 48 ++++ .../coalescence/system/controlDict | 72 +++++ .../coalescence/system/fvSchemes | 54 ++++ .../coalescence/system/fvSolution | 74 +++++ .../0/T.air1 | 25 ++ .../0/T.air2 | 25 ++ .../0/T.air3 | 25 ++ .../0/T.water | 25 ++ .../0/U.air1 | 25 ++ .../0/U.air2 | 25 ++ .../0/U.air3 | 25 ++ .../0/U.water | 25 ++ .../0/alpha.air1 | 26 ++ .../0/alpha.air2 | 26 ++ .../0/alpha.air3 | 26 ++ .../0/alpha.water | 25 ++ .../0/f.air1.bubbles | 26 ++ .../0/f.air2.bubbles | 26 ++ .../0/f.air3.bubbles | 26 ++ .../simultaneousCoalescenceAndBreakup/0/p | 25 ++ .../simultaneousCoalescenceAndBreakup/0/p_rgh | 25 ++ .../simultaneousCoalescenceAndBreakup/Allrun | 59 ++++ .../constant/g | 22 ++ .../constant/phaseProperties | 217 +++++++++++++++ .../constant/thermophysicalProperties.air1 | 49 ++++ .../constant/thermophysicalProperties.air2 | 49 ++++ .../constant/thermophysicalProperties.air3 | 49 ++++ .../constant/thermophysicalProperties.water | 53 ++++ .../constant/turbulenceProperties.air1 | 20 ++ .../constant/turbulenceProperties.air2 | 20 ++ .../constant/turbulenceProperties.air3 | 20 ++ .../constant/turbulenceProperties.water | 20 ++ .../simultaneousCoalescenceAndBreakup/plot | 25 ++ .../system/blockMeshDict | 48 ++++ .../system/controlDict | 89 ++++++ .../system/fvSchemes | 54 ++++ .../system/fvSolution | 74 +++++ .../populationBalanceModeling/README | 8 + .../populationBalanceModeling/drift/0/T.air | 25 ++ .../populationBalanceModeling/drift/0/T.water | 25 ++ .../populationBalanceModeling/drift/0/U.air | 25 ++ .../populationBalanceModeling/drift/0/U.water | 25 ++ .../drift/0/alpha.air | 26 ++ .../drift/0/f.air.bubbles | 26 ++ .../populationBalanceModeling/drift/0/p | 25 ++ .../populationBalanceModeling/drift/0/p_rgh | 25 ++ .../populationBalanceModeling/drift/Allrun | 39 +++ .../drift/constant/fvOptions | 43 +++ .../drift/constant/g | 22 ++ .../drift/constant/phaseProperties | 172 ++++++++++++ .../constant/thermophysicalProperties.air | 52 ++++ .../constant/thermophysicalProperties.water | 52 ++++ .../drift/constant/turbulenceProperties.air | 20 ++ .../drift/constant/turbulenceProperties.water | 20 ++ .../populationBalanceModeling/drift/plot | 35 +++ .../drift/system/blockMeshDict | 48 ++++ .../drift/system/controlDict | 88 ++++++ .../drift/system/fvSchemes | 54 ++++ .../drift/system/fvSolution | 74 +++++ .../negativeDrift/0/T.air | 25 ++ .../negativeDrift/0/T.water | 25 ++ .../negativeDrift/0/U.air | 25 ++ .../negativeDrift/0/U.water | 25 ++ .../negativeDrift/0/alpha.air | 26 ++ .../negativeDrift/0/f.air.bubbles | 26 ++ .../negativeDrift/0/p | 25 ++ .../negativeDrift/0/p_rgh | 25 ++ .../negativeDrift/Allrun | 39 +++ .../negativeDrift/constant/fvOptions | 43 +++ .../negativeDrift/constant/g | 22 ++ .../negativeDrift/constant/phaseProperties | 166 ++++++++++++ .../constant/thermophysicalProperties.air | 52 ++++ .../constant/thermophysicalProperties.water | 52 ++++ .../constant/turbulenceProperties.air | 20 ++ .../constant/turbulenceProperties.water | 20 ++ .../negativeDrift/plot | 35 +++ .../negativeDrift/system/blockMeshDict | 48 ++++ .../negativeDrift/system/controlDict | 88 ++++++ .../negativeDrift/system/fvSchemes | 54 ++++ .../negativeDrift/system/fvSolution | 74 +++++ tutorials/Allclean | 4 +- tutorials/Allrun | 6 +- .../0/{T.airI => T.air1} | 4 +- .../0/T.air2} | 4 +- .../0/U.air1} | 4 +- .../0/{U.airI => U.air2} | 4 +- .../0/{alpha.airI.orig => alpha.air1.orig} | 4 +- .../0/{alpha.airII.orig => alpha.air2.orig} | 4 +- .../bubbleColumnFixedPolydisperse/0/p_rgh | 6 +- .../constant/phaseProperties | 68 ++--- .../constant/thermophysicalProperties.air1} | 2 +- ...ies.airI => thermophysicalProperties.air2} | 2 +- .../constant/turbulenceProperties.air1 | 20 ++ .../constant/turbulenceProperties.air2 | 20 ++ .../system/controlDict | 8 +- .../system/setFieldsDict | 6 +- .../0/T.air1} | 4 +- .../0/{T.airII => T.air2} | 4 +- .../0/{U.airII => U.air1} | 4 +- .../0/U.air2} | 4 +- .../0/{alpha.airI.orig => alpha.air1.orig} | 4 +- .../0/{alpha.airII.orig => alpha.air2.orig} | 4 +- .../0/{f.airI.bubbles => f.air1.bubbles} | 4 +- .../0/{f.airII.bubbles => f.air2.bubbles} | 4 +- .../laminar/bubbleColumnPolydisperse/0/p_rgh | 6 +- .../constant/phaseProperties | 88 +++--- ...es.airII => thermophysicalProperties.air1} | 2 +- .../constant/thermophysicalProperties.air2} | 2 +- .../constant/turbulenceProperties.air1 | 20 ++ .../constant/turbulenceProperties.air2 | 20 ++ .../system/controlDict | 8 +- .../system/setFieldsDict | 6 +- 211 files changed, 7137 insertions(+), 138 deletions(-) create mode 120000 test/Allclean create mode 120000 test/Allrun create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/README create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/T.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/U.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/alpha.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air1.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air2.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/f.air3.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/0/p_rgh create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/Allrun create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/g create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/phaseProperties create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/thermophysicalProperties.water rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airI => test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air1 (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airI => test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air2 (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.airII => test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.air3 (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.airII => test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/constant/turbulenceProperties.water (95%) create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/plot create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/system/blockMeshDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/system/controlDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/system/fvSchemes create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/binaryBreakup/system/fvSolution create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/T.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/T.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/T.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/T.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/U.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/U.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/U.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/U.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/alpha.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/alpha.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/alpha.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/alpha.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/f.air1.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/f.air2.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/f.air3.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/p create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/0/p_rgh create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/Allrun create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/g create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/phaseProperties create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/thermophysicalProperties.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/constant/turbulenceProperties.water create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/plot create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/system/blockMeshDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/system/controlDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/system/fvSchemes create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/breakup/system/fvSolution create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/T.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/T.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/T.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/T.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/U.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/U.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/U.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/U.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/alpha.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/alpha.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/alpha.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/alpha.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/f.air1.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/f.air2.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/f.air3.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/p create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/0/p_rgh create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/Allrun create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/g create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/phaseProperties create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/thermophysicalProperties.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/constant/turbulenceProperties.water create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/plot create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/blockMeshDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/controlDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSchemes create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/coalescence/system/fvSolution create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/T.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/U.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/alpha.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air1.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air2.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/f.air3.bubbles create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/0/p_rgh create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/Allrun create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/g create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/phaseProperties create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/thermophysicalProperties.water create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air1 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air2 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.air3 create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/constant/turbulenceProperties.water create mode 100755 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/plot create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/blockMeshDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/controlDict create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSchemes create mode 100644 test/multiphase/reactingMultiphaseEulerFoam/populationBalanceModeling/simultaneousCoalescenceAndBreakup/system/fvSolution create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/README create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/T.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/U.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/alpha.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/f.air.bubbles create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/0/p_rgh create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/Allrun create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/fvOptions create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/g create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/phaseProperties create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/thermophysicalProperties.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/constant/turbulenceProperties.water create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/plot create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/system/blockMeshDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/system/controlDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/system/fvSchemes create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/drift/system/fvSolution create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/T.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/T.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/U.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/U.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/alpha.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/f.air.bubbles create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/p create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/0/p_rgh create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/Allrun create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/fvOptions create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/g create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/phaseProperties create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/thermophysicalProperties.water create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.air create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/constant/turbulenceProperties.water create mode 100755 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/plot create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/system/blockMeshDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/system/controlDict create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/system/fvSchemes create mode 100644 test/multiphase/reactingTwoPhaseEulerFoam/populationBalanceModeling/negativeDrift/system/fvSolution rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/{T.airI => T.air1} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnPolydisperse/0/T.airI => bubbleColumnFixedPolydisperse/0/T.air2} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnPolydisperse/0/U.airI => bubbleColumnFixedPolydisperse/0/U.air1} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/{U.airI => U.air2} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/{alpha.airI.orig => alpha.air1.orig} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/0/{alpha.airII.orig => alpha.air2.orig} (94%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnPolydisperse/constant/thermophysicalProperties.airI => bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.air1} (96%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/{thermophysicalProperties.airI => thermophysicalProperties.air2} (96%) create mode 100644 tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air1 create mode 100644 tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/turbulenceProperties.air2 rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnFixedPolydisperse/0/T.airII => bubbleColumnPolydisperse/0/T.air1} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{T.airII => T.air2} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{U.airII => U.air1} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnFixedPolydisperse/0/U.airII => bubbleColumnPolydisperse/0/U.air2} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{alpha.airI.orig => alpha.air1.orig} (95%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{alpha.airII.orig => alpha.air2.orig} (94%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{f.airI.bubbles => f.air1.bubbles} (94%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/0/{f.airII.bubbles => f.air2.bubbles} (94%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/{thermophysicalProperties.airII => thermophysicalProperties.air1} (96%) rename tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/{bubbleColumnFixedPolydisperse/constant/thermophysicalProperties.airII => bubbleColumnPolydisperse/constant/thermophysicalProperties.air2} (96%) create mode 100644 tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air1 create mode 100644 tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/turbulenceProperties.air2 diff --git a/bin/tools/CleanFunctions b/bin/tools/CleanFunctions index 3ac00c286..16b45561c 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 000000000..8775d6c59 --- /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 000000000..4920f4195 --- /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 000000000..abdaf191b --- /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 000000000..9e24e2282 --- /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 000000000..2deff50b5 --- /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 000000000..2ab28ff56 --- /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 000000000..c28847918 --- /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 000000000..dc7beaee5 --- /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 000000000..c298549e6 --- /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 000000000..3386be15c --- /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 000000000..65c08b123 --- /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 000000000..4986d6f25 --- /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 000000000..883a29cd8 --- /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 000000000..7d63260a5 --- /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 000000000..7ca124005 --- /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 000000000..a0d485d3d --- /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 000000000..62c6b255f --- /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 000000000..40817be16 --- /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 000000000..dbbc5c2da --- /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 000000000..c76174de2 --- /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 000000000..076ba12d8 --- /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 000000000..508d65849 --- /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 000000000..8ecbc1b6c --- /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 000000000..5254f4787 --- /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 000000000..5dacc2e88 --- /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 000000000..0219c745d --- /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 000000000..bd7ad0b32 --- /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 6bc136989..fda3c5b86 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 6bc136989..f6bb96159 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 8c3b8d369..24af3b4c6 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 8c3b8d369..7f0d75a82 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 000000000..9a081b873 --- /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 000000000..508d65849 --- /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 000000000..5ec6f774c --- /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 000000000..5254f4787 --- /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 000000000..5dacc2e88 --- /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 000000000..0219c745d --- /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 000000000..bd7ad0b32 --- /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 000000000..fda3c5b86 --- /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 000000000..f6bb96159 --- /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 000000000..24af3b4c6 --- /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 000000000..7f0d75a82 --- /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 000000000..9a081b873 --- /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 000000000..508d65849 --- /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 000000000..e6b4288c8 --- /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 000000000..5254f4787 --- /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 000000000..5dacc2e88 --- /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 000000000..0219c745d --- /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 000000000..bd7ad0b32 --- /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 000000000..fda3c5b86 --- /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 000000000..f6bb96159 --- /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 000000000..24af3b4c6 --- /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 000000000..7f0d75a82 --- /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 000000000..e43f9fc68 --- /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 000000000..c925e3542 --- /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 000000000..c430fd85d --- /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 000000000..809d86bb6 --- /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 000000000..db2be3b84 --- /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 000000000..9e24e2282 --- /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 000000000..2deff50b5 --- /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 000000000..2ab28ff56 --- /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 000000000..c28847918 --- /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 000000000..dc7beaee5 --- /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 000000000..c298549e6 --- /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 000000000..3386be15c --- /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 000000000..65c08b123 --- /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 000000000..25f229893 --- /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 000000000..883a29cd8 --- /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 000000000..e12aa013f --- /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 000000000..9b1c58319 --- /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 000000000..a0d485d3d --- /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 000000000..62c6b255f --- /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 000000000..40817be16 --- /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 000000000..dbbc5c2da --- /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 000000000..c76174de2 --- /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 000000000..6d4e4d38a --- /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 000000000..508d65849 --- /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 000000000..a43d0afd9 --- /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 000000000..5254f4787 --- /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 000000000..5dacc2e88 --- /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 000000000..0219c745d --- /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 000000000..bd7ad0b32 --- /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 000000000..fda3c5b86 --- /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 000000000..f6bb96159 --- /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 000000000..24af3b4c6 --- /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 000000000..7f0d75a82 --- /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 000000000..a888f9470 --- /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 000000000..c925e3542 --- /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 000000000..b1ac53709 --- /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 000000000..809d86bb6 --- /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 000000000..d0808e00b --- /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 000000000..52d47def2 --- /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 000000000..854dfc618 --- /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 000000000..c28847918 --- /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 000000000..726323aa3 --- /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 000000000..65c08b123 --- /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 000000000..f9d89cc08 --- /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 000000000..61a9ba40e --- /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 000000000..dbbc5c2da --- /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 000000000..c76174de2 --- /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 000000000..076ba12d8 --- /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 000000000..0455a1a64 --- /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 000000000..508d65849 --- /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 000000000..7dfd21d74 --- /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 000000000..4c8bb3629 --- /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 000000000..4a566672f --- /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 000000000..1296429b7 --- /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 000000000..7f0d75a82 --- /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 000000000..517570ccd --- /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 000000000..c203d8c1a --- /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 000000000..508d65849 --- /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 000000000..14e2b7db7 --- /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 000000000..4c8bb3629 --- /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 000000000..4a566672f --- /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 000000000..1296429b7 --- /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 000000000..7f0d75a82 --- /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 000000000..ec7e48628 --- /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 3d377d095..ec92afac2 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 3d377d095..d1cafcb7e 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 da969ff41..ced8937c0 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 da969ff41..a83d1299e 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 203cc22fc..4a3078939 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 019c483a0..5c8b4c69c 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 8d1b65707..5dcbfa90d 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 62d61e05e..6c09d3376 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 e11ccef1b..7426132d8 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 e11ccef1b..62bff79a0 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 000000000..fda3c5b86 --- /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 000000000..f6bb96159 --- /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 4aa6dd1f7..a364cd024 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 6bd134c53..6c8e902e3 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 03b1f9a28..ec92afac2 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 03b1f9a28..d1cafcb7e 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 8c2d21224..ced8937c0 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 8c2d21224..a83d1299e 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 203cc22fc..4a3078939 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 019c483a0..5c8b4c69c 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 9311ec38f..c3a6d0ec9 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 643e75c2b..5494f7f5f 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 8d1b65707..5dcbfa90d 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 96469c3f7..35f0c3544 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 a190c7d7a..7426132d8 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 a190c7d7a..62bff79a0 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 000000000..fda3c5b86 --- /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 000000000..f6bb96159 --- /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 1dcc62d12..3bda239ba 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 6bd134c53..6c8e902e3 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 ); }