From ef885b407c870eb330ebc9139519da386ef0acc2 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Thu, 22 Mar 2018 12:17:57 +0000 Subject: [PATCH] reactingEulerFoam: Corrected blending in tutorials Sub-model blending should be set such that the sum of all the blending coefficients equals one. If there are three models specified for a phase pair (e.g., (air in water), (water in air) and (air and water)), then the sum-to-one constraint is guaranteed by the blending functions. Frequently, however, the symmetric model ((air and water) in this example) is omitted. In that case, the blending coefficients should be selected so that the sum of just the two non-symmetric coefficients equal one. In the case of linear blending, this means setting the minimum partially continuous alpha to one-minus the fully continuous value of the opposite phase. For example: blending { default { type linear; minFullyContinuousAlpha.air 0.7; minPartlyContinuousAlpha.air 0.3; minFullyContinuousAlpha.water 0.7; minPartlyContinuousAlpha.water 0.3; } } The reactingTwoPhaseEulerFoam and reactingMultiPhaseEulerFoam tutorials have been modified to adhere to this principle. --- .../wallBoiling1D_2phase/constant/phaseProperties | 12 ------------ .../wallBoiling1D_3phase/constant/phaseProperties | 12 ------------ .../laminar/bubbleColumn/constant/phaseProperties | 9 +++++++++ .../constant/phaseProperties | 11 +++++++++++ .../constant/phaseProperties | 11 +++++++++++ .../laminar/mixerVessel2D/constant/phaseProperties | 13 +++++++++++++ .../LES/bubbleColumn/constant/phaseProperties | 9 +++++++++ .../RAS/bubbleColumn/constant/phaseProperties | 9 +++++++++ .../constant/phaseProperties | 10 +++------- .../constant/phaseProperties | 9 +++++++++ .../laminar/bubbleColumn/constant/phaseProperties | 9 +++++++++ .../constant/phaseProperties | 14 +++++--------- .../constant/phaseProperties | 14 +++++--------- .../bubbleColumnIATE/constant/phaseProperties | 9 +++++++++ .../laminar/injection/constant/phaseProperties | 9 +++++++++ .../laminar/mixerVessel2D/constant/phaseProperties | 9 +++++++++ 16 files changed, 120 insertions(+), 49 deletions(-) diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_2phase/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_2phase/constant/phaseProperties index 083971043..698ff5793 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_2phase/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_2phase/constant/phaseProperties @@ -56,18 +56,6 @@ blending type none; continuousPhase liquid; } - - heatTransfer - { - type none; - continuousPhase liquid; - } - - massTransfer - { - type none; - continuousPhase liquid; - } } surfaceTension diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_3phase/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_3phase/constant/phaseProperties index f3e0f3388..b9f16eaf5 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_3phase/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/RAS/wallBoiling1D_3phase/constant/phaseProperties @@ -71,18 +71,6 @@ blending type none; continuousPhase liquid; } - - heatTransfer - { - type none; - continuousPhase liquid; - } - - massTransfer - { - type none; - continuousPhase liquid; - } } surfaceTension diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties index 98192e477..a7e648285 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties index 6c09d3376..4653dff1a 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnFixedPolydisperse/constant/phaseProperties @@ -60,6 +60,17 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air1 0.7; + minPartlyContinuousAlpha.air1 0.3; + minFullyContinuousAlpha.air2 0.7; + minPartlyContinuousAlpha.air2 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air1 0.7; diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties index 35f0c3544..f376397ec 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/bubbleColumnPolydisperse/constant/phaseProperties @@ -140,6 +140,17 @@ blending minFullyContinuousAlpha.water 0.7; minPartlyContinuousAlpha.water 0.5; } + + drag + { + type linear; + minFullyContinuousAlpha.air1 0.7; + minPartlyContinuousAlpha.air1 0.3; + minFullyContinuousAlpha.air2 0.7; + minPartlyContinuousAlpha.air2 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } } surfaceTension diff --git a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties index 371c656ea..1a63c74cb 100644 --- a/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties +++ b/tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties @@ -71,6 +71,19 @@ air blending { default + { + type linear; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + minFullyContinuousAlpha.oil 0.7; + minPartlyContinuousAlpha.oil 0.3; + minFullyContinuousAlpha.mercury 0.7; + minPartlyContinuousAlpha.mercury 0.3; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + } + + drag { type linear; minFullyContinuousAlpha.water 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties index db138464e..2e2cdb3ce 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties index db138464e..2e2cdb3ce 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnEvaporatingReacting/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnEvaporatingReacting/constant/phaseProperties index 23358cfec..0dcc5e83d 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnEvaporatingReacting/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnEvaporatingReacting/constant/phaseProperties @@ -51,9 +51,9 @@ blending { type linear; minFullyContinuousAlpha.gas 0.7; - minPartlyContinuousAlpha.gas 0.5; + minPartlyContinuousAlpha.gas 0.3; minFullyContinuousAlpha.liquid 0.7; - minPartlyContinuousAlpha.liquid 0.5; + minPartlyContinuousAlpha.liquid 0.3; } heatTransfer @@ -67,11 +67,7 @@ blending massTransfer { - type linear; - minFullyContinuousAlpha.gas 1; - minPartlyContinuousAlpha.gas 0; - minFullyContinuousAlpha.liquid 1; - minPartlyContinuousAlpha.liquid 0; + $heatTransfer; } } diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties index f2017e44d..80e345cbe 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/bubbleColumnPolydisperse/constant/phaseProperties @@ -113,6 +113,15 @@ populationBalanceCoeffs blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties index 3f7a06fbc..92508e67a 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating/constant/phaseProperties index eb0bbd970..91f7bf3a3 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporating/constant/phaseProperties @@ -50,13 +50,13 @@ blending default { type linear; - minPartlyContinuousAlpha.gas 0.5; + minPartlyContinuousAlpha.gas 0.3; minFullyContinuousAlpha.gas 0.7; - minPartlyContinuousAlpha.liquid 0.5; + minPartlyContinuousAlpha.liquid 0.3; minFullyContinuousAlpha.liquid 0.7; } - heatTransferModel + heatTransfer { type linear; minPartlyContinuousAlpha.gas 0; @@ -65,13 +65,9 @@ blending minFullyContinuousAlpha.liquid 1; } - massTransferModel + massTransfer { - type linear; - minPartlyContinuousAlpha.gas 0; - minFullyContinuousAlpha.gas 1; - minPartlyContinuousAlpha.liquid 0; - minFullyContinuousAlpha.liquid 1; + $heatTransfer; } } diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporatingDissolving/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporatingDissolving/constant/phaseProperties index 6b628569c..17a75c8bb 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporatingDissolving/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnEvaporatingDissolving/constant/phaseProperties @@ -52,12 +52,12 @@ blending { type linear; minFullyContinuousAlpha.gas 0.7; - minPartlyContinuousAlpha.gas 0.5; + minPartlyContinuousAlpha.gas 0.3; minFullyContinuousAlpha.liquid 0.7; - minPartlyContinuousAlpha.liquid 0.5; + minPartlyContinuousAlpha.liquid 0.3; } - heatTransferModel + heatTransfer { type linear; minFullyContinuousAlpha.gas 1; @@ -66,13 +66,9 @@ blending minPartlyContinuousAlpha.liquid 0; } - massTransferModel + massTransfer { - type linear; - minFullyContinuousAlpha.gas 1; - minPartlyContinuousAlpha.gas 0; - minFullyContinuousAlpha.liquid 1; - minPartlyContinuousAlpha.liquid 0; + $heatTransfer; } } diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties index e46eefe31..2e3ef5368 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties @@ -69,6 +69,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/injection/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/injection/constant/phaseProperties index efcc99353..1412341ad 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/injection/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/injection/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7; diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties index e374b0285..02325b51f 100644 --- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties +++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties @@ -47,6 +47,15 @@ water blending { default + { + type linear; + minFullyContinuousAlpha.air 0.7; + minPartlyContinuousAlpha.air 0.3; + minFullyContinuousAlpha.water 0.7; + minPartlyContinuousAlpha.water 0.3; + } + + drag { type linear; minFullyContinuousAlpha.air 0.7;