zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

INT: Updated code related to change in pos(x) function

Browse Source
This commit is contained in:
Andrew Heather
2017-09-07 14:57:10 +01:00
parent 56bfc75949
commit f549e28fb8
4 changed files with 16 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
src/OpenFOAM/orientedType/orientedType.C
View File

@ -362,12 +362,24 @@ Foam::orientedType Foam::pos(const orientedType& ot)
} }
Foam::orientedType Foam::pos0(const orientedType& ot)
{
return ot;
}
Foam::orientedType Foam::neg(const orientedType& ot) Foam::orientedType Foam::neg(const orientedType& ot)
{ {
return ot; return ot;
} }
Foam::orientedType Foam::neg0(const orientedType& ot)
{
return ot;
}
Foam::orientedType Foam::posPart(const orientedType& ot) Foam::orientedType Foam::posPart(const orientedType& ot)
{ {
return ot; return ot;

2
src/OpenFOAM/orientedType/orientedType.H
View File

@ -169,7 +169,9 @@ orientedType magSqr(const orientedType& ot);
orientedType mag(const orientedType& ot); orientedType mag(const orientedType& ot);
orientedType sign(const orientedType& ot); orientedType sign(const orientedType& ot);
orientedType pos(const orientedType& ot); orientedType pos(const orientedType& ot);
orientedType pos0(const orientedType& ot);
orientedType neg(const orientedType& ot); orientedType neg(const orientedType& ot);
orientedType neg0(const orientedType& ot);
orientedType posPart(const orientedType& ot); orientedType posPart(const orientedType& ot);
orientedType negPart(const orientedType& ot); orientedType negPart(const orientedType& ot);
orientedType inv(const orientedType& ot); orientedType inv(const orientedType& ot);

2
src/TurbulenceModels/turbulenceModels/DES/SpalartAllmarasIDDES/SpalartAllmarasIDDES.C
View File

@ -135,7 +135,7 @@ tmp<volScalarField> SpalartAllmarasIDDES<BasicTurbulenceModel>::dTilda
tmp<volScalarField> fB = min(2*pow(expTerm, -9.0), scalar(1)); tmp<volScalarField> fB = min(2*pow(expTerm, -9.0), scalar(1));
tmp<volScalarField> fe1 = tmp<volScalarField> fe1 =
2*(pos(alpha)*pow(expTerm, -11.09) + neg(alpha)*pow(expTerm, -9.0)); 2*(pos0(alpha)*pow(expTerm, -11.09) + neg(alpha)*pow(expTerm, -9.0));
tmp<volScalarField> fe2 = 1 - max(ft(magGradU), fl(magGradU)); tmp<volScalarField> fe2 = 1 - max(ft(magGradU), fl(magGradU));
tmp<volScalarField> fe = max(fe1 - 1, scalar(0))*psi*fe2; tmp<volScalarField> fe = max(fe1 - 1, scalar(0))*psi*fe2;

2
src/TurbulenceModels/turbulenceModels/DES/kOmegaSSTIDDES/kOmegaSSTIDDES.C
View File

@ -134,7 +134,7 @@ tmp<volScalarField> kOmegaSSTIDDES<BasicTurbulenceModel>::dTilda
tmp<volScalarField> fB = min(2*pow(expTerm, -9.0), scalar(1)); tmp<volScalarField> fB = min(2*pow(expTerm, -9.0), scalar(1));
tmp<volScalarField> fe1 = tmp<volScalarField> fe1 =
2*(pos(alpha)*pow(expTerm, -11.09) + neg(alpha)*pow(expTerm, -9.0)); 2*(pos0(alpha)*pow(expTerm, -11.09) + neg(alpha)*pow(expTerm, -9.0));
tmp<volScalarField> fe2 = 1 - max(ft(magGradU), fl(magGradU)); tmp<volScalarField> fe2 = 1 - max(ft(magGradU), fl(magGradU));
tmp<volScalarField> fe = max(fe1 - 1, scalar(0))*fe2; tmp<volScalarField> fe = max(fe1 - 1, scalar(0))*fe2;