From d27dbbf096df91926ff3aa5afbd2dd3933177144 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Tue, 9 May 2017 08:48:38 +0100 Subject: [PATCH] Euler-Euler: turbulentDispersionModels: Burns: Minor re-formulation to improve stability. Resolves bug report --- .../turbulentDispersionModels/Burns/Burns.C | 8 ++++++-- .../turbulentDispersionModels/Burns/Burns.C | 8 ++++++-- 2 files changed, 12 insertions(+), 4 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C index 1f0fc837cb..d06a7929f2 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2014-2015 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2014-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -101,7 +101,11 @@ Foam::turbulentDispersionModels::Burns::D() const *sqr(pair_.dispersed().d()) ) *pair_.continuous().rho() - *(1.0 + pair_.dispersed()/max(pair_.continuous(), residualAlpha_)); + *pair_.dispersed() + *( + 1.0/max(pair_.dispersed(), residualAlpha_) + + 1.0/max(pair_.continuous(), residualAlpha_) + ); } diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C index 093013cca3..33db7656cb 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Burns/Burns.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2014-2015 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2014-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -101,7 +101,11 @@ Foam::turbulentDispersionModels::Burns::D() const *sqr(pair_.dispersed().d()) ) *pair_.continuous().rho() - *(1.0 + pair_.dispersed()/max(pair_.continuous(), residualAlpha_)); + *pair_.dispersed() + *( + 1.0/max(pair_.dispersed(), residualAlpha_) + + 1.0/max(pair_.continuous(), residualAlpha_) + ); }