zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

kOmega: Added support for k and omega sources in derived models

Browse Source 
Updated sources to use internal fields.
This commit is contained in:
Henry Weller
2019-01-10 10:39:52 +00:00
parent 66c26af922
commit f5d684b034
2 changed files with 18 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/TurbulenceModels/turbulenceModels/RAS/kOmega/kOmega.C
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -195,13 +195,16 @@ void kOmega<BasicTurbulenceModel>::correct()
eddyViscosity<RASModel<BasicTurbulenceModel>>::correct(); eddyViscosity<RASModel<BasicTurbulenceModel>>::correct();
volScalarField divU(fvc::div(fvc::absolute(this->phi(), U))); volScalarField::Internal divU
(
fvc::div(fvc::absolute(this->phi(), U))().v()
);
tmp<volTensorField> tgradU = fvc::grad(U); tmp<volTensorField> tgradU = fvc::grad(U);
volScalarField G volScalarField::Internal G
( (
this->GName(), this->GName(),
nut*(tgradU() && dev(twoSymm(tgradU()))) nut.v()*(dev(twoSymm(tgradU().v())) && tgradU().v())
); );
tgradU.clear(); tgradU.clear();
@ -215,9 +218,10 @@ void kOmega<BasicTurbulenceModel>::correct()
+ fvm::div(alphaRhoPhi, omega_) + fvm::div(alphaRhoPhi, omega_)
- fvm::laplacian(alpha*rho*DomegaEff(), omega_) - fvm::laplacian(alpha*rho*DomegaEff(), omega_)
== ==
gamma_*alpha*rho*G*omega_/k_ gamma_*alpha()*rho()*G*omega_()/k_()
- fvm::SuSp(((2.0/3.0)*gamma_)*alpha*rho*divU, omega_) - fvm::SuSp(((2.0/3.0)*gamma_)*alpha()*rho()*divU, omega_)
- fvm::Sp(beta_*alpha*rho*omega_, omega_) - fvm::Sp(beta_*alpha()*rho()*omega_(), omega_)
+ omegaSource()
+ fvOptions(alpha, rho, omega_) + fvOptions(alpha, rho, omega_)
); );
@ -236,9 +240,10 @@ void kOmega<BasicTurbulenceModel>::correct()
+ fvm::div(alphaRhoPhi, k_) + fvm::div(alphaRhoPhi, k_)
- fvm::laplacian(alpha*rho*DkEff(), k_) - fvm::laplacian(alpha*rho*DkEff(), k_)
== ==
alpha*rho*G alpha()*rho()*G
- fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_) - fvm::SuSp((2.0/3.0)*alpha()*rho()*divU, k_)
- fvm::Sp(Cmu_*alpha*rho*omega_, k_) - fvm::Sp(Cmu_*alpha()*rho()*omega_(), k_)
+ kSource()
+ fvOptions(alpha, rho, k_) + fvOptions(alpha, rho, k_)
); );

4
src/TurbulenceModels/turbulenceModels/RAS/kOmega/kOmega.H
View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org \\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -97,6 +97,8 @@ protected:
// Protected Member Functions // Protected Member Functions
virtual void correctNut(); virtual void correctNut();
virtual tmp<fvScalarMatrix> kSource() const;
virtual tmp<fvScalarMatrix> omegaSource() const;
public: public: