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

Changed stabilisation from rootVSmall to small for to avoid FPEs with extended precision

Browse Source
This commit is contained in:
Henry Weller
2023-06-16 21:36:05 +01:00
parent 025b96e217
commit ef85d538e9
4 changed files with 11 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/modules/multiphaseEuler/interfacialModels/liftModels/SaffmanMei/SaffmanMei.C
View File

@ -63,18 +63,18 @@ Foam::liftModels::SaffmanMei::~SaffmanMei()
Foam::tmp<Foam::volScalarField> Foam::liftModels::SaffmanMei::Cl() const
{
volScalarField Re(max(interface_.Re(), residualRe_));
volScalarField Rew
const volScalarField Re(max(interface_.Re(), residualRe_));
const volScalarField Rew
(
mag(fvc::curl(interface_.continuous().U()))
*sqr(interface_.dispersed().d())
/(
interface_.continuous().thermo().nu()
+ dimensionedScalar(dimViscosity, rootVSmall)
+ dimensionedScalar(dimViscosity, small)
)
);
volScalarField Cld
const volScalarField Cld
(
neg0(Re - 40)*6.46
*(
@ -84,7 +84,7 @@ Foam::tmp<Foam::volScalarField> Foam::liftModels::SaffmanMei::Cl() const
+ pos(Re - 40)*6.46*0.0524*sqrt(0.5*(Rew/Re)*Re)
);
return 3/(twoPi*sqrt(Rew + rootVSmall))*Cld;
return 3/(twoPi*sqrt(Rew + small))*Cld;
}