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

BUG: FriedrichModel: correct relative-speed calculation (fixes #3436)

Browse Source
This commit is contained in:
Kutalmis Bercin
2025-09-19 10:55:05 +01:00
parent f7fd9f8186
commit 2d522e921f
1 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
src/regionFaModels/liquidFilm/subModels/kinematic/injectionModel/filmSeparation/filmSeparationModels/FriedrichModel/FriedrichModel.C
View File

@ -5,7 +5,7 @@
\\ / A nd | www.openfoam.com \\ / A nd | www.openfoam.com
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2024 OpenCFD Ltd. Copyright (C) 2024-2025 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -401,8 +401,8 @@ tmp<scalarField> FriedrichModel::Fratio() const
const scalar Re = h[i]*mag(Uf[i])*rho[i]/mu[i]; const scalar Re = h[i]*mag(Uf[i])*rho[i]/mu[i];
// Weber number (FLW:Eq. 17) // Weber number (FLW:Eq. 17)
const scalar We = const vector Urel(Up[i] - Uf[i]);
h[i]*rhop_*sqr(mag(Up[i]) - mag(Uf[i]))/(2.0*sigma[i]); const scalar We = h[i]*rhop_*sqr(mag(Urel))/(2.0*sigma[i]);
// Characteristic breakup length (FLW:Eq. 15) // Characteristic breakup length (FLW:Eq. 15)
const scalar Lb = const scalar Lb =
@ -455,8 +455,8 @@ tmp<scalarField> FriedrichModel::Fratio() const
const scalar Re = hp[i]*mag(Ufp[i])*rhop[i]/mup[i]; const scalar Re = hp[i]*mag(Ufp[i])*rhop[i]/mup[i];
// Weber number (FLW:Eq. 17) // Weber number (FLW:Eq. 17)
const scalar We = const vector Urelp(Upp[i] - Ufp[i]);
hp[i]*rhop_*sqr(mag(Upp[i]) - mag(Ufp[i]))/(2.0*sigmap[i]); const scalar We = hp[i]*rhop_*sqr(mag(Urelp))/(2.0*sigmap[i]);
// Characteristic breakup length (FLW:Eq. 15) // Characteristic breakup length (FLW:Eq. 15)
const scalar Lb = const scalar Lb =