From 42612f2db58b0aaec1b876c75eed194083cd2576 Mon Sep 17 00:00:00 2001 From: Henry Date: Wed, 4 Apr 2012 17:44:18 +0100 Subject: [PATCH] turbulenceModels/compressible/RAS: ensure transport equations produce bounded solutions even if continuity is not obeyed --- src/turbulenceModels/compressible/RAS/LRR/LRR.C | 4 +++- .../compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C | 4 +++- .../compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C | 4 +++- .../compressible/RAS/RNGkEpsilon/RNGkEpsilon.C | 4 +++- .../compressible/RAS/SpalartAllmaras/SpalartAllmaras.C | 1 + src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C | 4 +++- src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C | 4 +++- .../compressible/RAS/realizableKE/realizableKE.C | 4 +++- 8 files changed, 22 insertions(+), 7 deletions(-) diff --git a/src/turbulenceModels/compressible/RAS/LRR/LRR.C b/src/turbulenceModels/compressible/RAS/LRR/LRR.C index e0958e376e..be15b76474 100644 --- a/src/turbulenceModels/compressible/RAS/LRR/LRR.C +++ b/src/turbulenceModels/compressible/RAS/LRR/LRR.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -351,6 +351,7 @@ void LRR::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) //- fvm::laplacian(Ceps*rho_*(k_/epsilon_)*R_, epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == @@ -393,6 +394,7 @@ void LRR::correct() ( fvm::ddt(rho_, R_) + fvm::div(phi_, R_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), R_) //- fvm::laplacian(Cs*rho_*(k_/epsilon_)*R_, R_) - fvm::laplacian(DREff(), R_) + fvm::Sp(Clrr1_*rho_*epsilon_/k_, R_) diff --git a/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C b/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C index 6aef97660a..fb9bee6570 100644 --- a/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C +++ b/src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -389,6 +389,7 @@ void LaunderGibsonRSTM::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) //- fvm::laplacian(Ceps*rho_*(k_/epsilon_)*R_, epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == @@ -432,6 +433,7 @@ void LaunderGibsonRSTM::correct() ( fvm::ddt(rho_, R_) + fvm::div(phi_, R_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), R_) //- fvm::laplacian(Cs*rho_*(k_/epsilon_)*R_, R_) - fvm::laplacian(DREff(), R_) + fvm::Sp(Clg1_*rho_*epsilon_/k_, R_) diff --git a/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C b/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C index 996f0e5d2c..d03f6b5de9 100644 --- a/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C +++ b/src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -314,6 +314,7 @@ void LaunderSharmaKE::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == C1_*G*epsilon_/k_ + fvm::SuSp((C3_ - 2.0/3.0*C1_)*rho_*divU, epsilon_) @@ -333,6 +334,7 @@ void LaunderSharmaKE::correct() ( fvm::ddt(rho_, k_) + fvm::div(phi_, k_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), k_) - fvm::laplacian(DkEff(), k_) == G - fvm::SuSp(2.0/3.0*rho_*divU, k_) diff --git a/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C b/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C index d1d751e54d..96155d40a4 100644 --- a/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C +++ b/src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -321,6 +321,7 @@ void RNGkEpsilon::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == (C1_ - R)*G*epsilon_/k_ @@ -342,6 +343,7 @@ void RNGkEpsilon::correct() ( fvm::ddt(rho_, k_) + fvm::div(phi_, k_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), k_) - fvm::laplacian(DkEff(), k_) == G - fvm::SuSp(2.0/3.0*rho_*divU, k_) diff --git a/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C b/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C index 358e6810ec..b76872cfa1 100644 --- a/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C +++ b/src/turbulenceModels/compressible/RAS/SpalartAllmaras/SpalartAllmaras.C @@ -416,6 +416,7 @@ void SpalartAllmaras::correct() ( fvm::ddt(rho_, nuTilda_) + fvm::div(phi_, nuTilda_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), nuTilda_) - fvm::laplacian(DnuTildaEff(), nuTilda_) - Cb2_/sigmaNut_*rho_*magSqr(fvc::grad(nuTilda_)) == diff --git a/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C b/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C index 2db4878a4c..aa63cf9f23 100644 --- a/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C +++ b/src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -292,6 +292,7 @@ void kEpsilon::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == C1_*G*epsilon_/k_ @@ -313,6 +314,7 @@ void kEpsilon::correct() ( fvm::ddt(rho_, k_) + fvm::div(phi_, k_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), k_) - fvm::laplacian(DkEff(), k_) == G diff --git a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C index e4032d854c..e1ffc23d8d 100644 --- a/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C +++ b/src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -411,6 +411,7 @@ void kOmegaSST::correct() ( fvm::ddt(rho_, omega_) + fvm::div(phi_, omega_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), omega_) - fvm::laplacian(DomegaEff(F1), omega_) == rhoGammaF1*GbyMu @@ -435,6 +436,7 @@ void kOmegaSST::correct() ( fvm::ddt(rho_, k_) + fvm::div(phi_, k_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), k_) - fvm::laplacian(DkEff(F1), k_) == min(G, (c1_*betaStar_)*rho_*k_*omega_) diff --git a/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C b/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C index c169831c51..bb0f3041fa 100644 --- a/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C +++ b/src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -331,6 +331,7 @@ void realizableKE::correct() ( fvm::ddt(rho_, epsilon_) + fvm::div(phi_, epsilon_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), epsilon_) - fvm::laplacian(DepsilonEff(), epsilon_) == C1*rho_*magS*epsilon_ @@ -355,6 +356,7 @@ void realizableKE::correct() ( fvm::ddt(rho_, k_) + fvm::div(phi_, k_) + - fvm::Sp(fvc::ddt(rho_) + fvc::div(phi_), k_) - fvm::laplacian(DkEff(), k_) == G - fvm::SuSp(2.0/3.0*rho_*divU, k_)