compressibleTwoPhaseEulerFoam tutorials: Updated fvSchemes and improved naming of kineticTheory stress correction

This commit is contained in:
Henry
2013-07-28 22:59:22 +01:00
parent 0df1fa0642
commit 7a3181c593
6 changed files with 37 additions and 46 deletions

View File

@ -506,7 +506,8 @@ Foam::tmp<Foam::fvVectorMatrix> Foam::kineticTheoryModel::divDevRhoReff
- fvc::div - fvc::div
( (
(this->rho_*this->nut_)*dev2(T(fvc::grad(U))) (this->rho_*this->nut_)*dev2(T(fvc::grad(U)))
+ ((this->rho_*lambda_)*fvc::div(this->phi_))*symmTensor::I + ((this->rho_*lambda_)*fvc::div(this->phi_))
*dimensioned<symmTensor>("I", dimless, symmTensor::I)
) )
); );
} }

View File

@ -27,18 +27,17 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phi,alpha.air) Gauss vanLeer; div(phi,alpha.air) Gauss vanLeer;
div(phir,alpha.air) Gauss vanLeer; div(phir,alpha.air) Gauss vanLeer;
"div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1;
"div\(\(alpha.*Rc\)\)" Gauss linear; "div\(phi.*,.*rho.*\)" Gauss linear;
"div\(phi.*,.*rho.*\)" Gauss linear;
"div\(alphaPhi.*,h.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,K.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
"div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; "div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
} }

View File

@ -27,25 +27,19 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
"div\(phi,alpha.*\)" Gauss vanLeer; "div\(phi,alpha.*\)" Gauss vanLeer;
"div\(phir,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss vanLeer;
"div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1;
"div\(\(alpha.*Rc\)\)" Gauss linear;
"div\(phi.*,.*rho.*\)" Gauss linear; "div\(phi.*,.*rho.*\)" Gauss linear;
div(alphaPhi.particles,Theta.particles) Gauss limitedLinear 1;
div(alphaPhi.air,k.air) Gauss limitedLinear 1;
div(alphaPhi.air,epsilon.air) Gauss limitedLinear 1;
"div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
"div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; "div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
div(((nut.particles*dev2(T(grad(U.particles))))+((lambda.particles*div(phi.particles))*(1,0,0,1,0,1)))) Gauss linear;
} }
laplacianSchemes laplacianSchemes

View File

@ -27,18 +27,17 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phi,alpha.air) Gauss vanLeer; div(phi,alpha.air) Gauss vanLeer;
div(phir,alpha.air) Gauss vanLeer; div(phir,alpha.air) Gauss vanLeer;
"div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1;
"div\(\(alpha.*Rc\)\)" Gauss linear; "div\(phi.*,.*rho.*\)" Gauss linear;
"div\(phi.*,.*rho.*\)" Gauss linear;
"div\(alphaPhi.*,h.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,K.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
"div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; "div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
} }
@ -62,7 +61,6 @@ fluxRequired
{ {
default no; default no;
p ; p ;
alpha1;
} }

View File

@ -27,18 +27,17 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phi,alpha.air) Gauss vanLeer; div(phi,alpha.air) Gauss vanLeer;
div(phir,alpha.air) Gauss vanLeer; div(phir,alpha.air) Gauss vanLeer;
"div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1;
"div\(\(alpha.*Rc\)\)" Gauss linear; "div\(phi.*,.*rho.*\)" Gauss linear;
"div\(phi.*,.*rho.*\)" Gauss linear;
"div\(alphaPhi.*,h.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,K.*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1;

View File

@ -27,25 +27,25 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
"div\(phi,alpha.*\)" Gauss vanLeer; "div\(phi,alpha.*\)" Gauss vanLeer;
"div\(phir,alpha.*\)" Gauss vanLeer; "div\(phir,alpha.*\)" Gauss vanLeer;
"div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1; "div\(alphaPhi.*,U.*\)" Gauss limitedLinearV 1;
"div\(phi.*,U.*\)" Gauss limitedLinearV 1; "div\(phi.*,U.*\)" Gauss limitedLinearV 1;
"div\(\(alpha.*Rc\)\)" Gauss linear;
"div\(phi.*,.*rho.*\)" Gauss linear; "div\(phi.*,.*rho.*\)" Gauss linear;
div(alphaPhi.particles,Theta.particles) Gauss limitedLinear 1;
div(alphaPhi.air,k.air) Gauss limitedLinear 1;
div(alphaPhi.air,epsilon.air) Gauss limitedLinear 1;
"div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(h|e).*\)" Gauss limitedLinear 1;
"div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1; "div\(alphaPhi.*,(K.*|p)\)" Gauss limitedLinear 1;
div(alphaPhi.particles,Theta.particles) Gauss limitedLinear 1;
"div\(alphaPhi.*,(k|epsilon).*\)" Gauss limitedLinear 1;
"div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear; "div\(\(\(alpha.*nuEff.*\)*dev2\(T\(grad\(U.*\)\)\)\)\)" Gauss linear;
div(((nut.particles*dev2(T(grad(U.particles))))+((lambda.particles*div(phi.particles))*(1,0,0,1,0,1)))) Gauss linear; div(((nut.particles*dev2(T(grad(U.particles))))+((lambda.particles*div(phi.particles))*I))) Gauss linear;
} }
laplacianSchemes laplacianSchemes