View Issue Details

IDProjectCategoryView StatusLast Update
0003394OpenFOAMBugpublic2019-11-26 08:53
ReporteraatxutegiAssigned Towill 
PrioritynormalSeverityminorReproducibilityalways
Status closedResolutionno change required 
PlatformGNU/LinuxOSUbuntuOS Version15.04
Product Version7 
Fixed in Version 
Summary0003394: Reitz & Diwakar coefficients
DescriptionThe original paper (attached) defined the stripping breakage cutoff value as We/ sqrt(Re) while the code is implemented as We * sqrt(Re). In addition it seems that the stripping breakup time should be multiplied by 0.5 (line 118) in order to be cosnistent with the authors as they used the radiuos instead of the droplet diamater. I might be missing
something but any insight would be appreciated.
Steps To Reproducesrc/lagrangian/spray/submodels/BreakupModel/ReitzDiwakar/ReitzDiwakar.C file should be checked in lines 106 and 118
TagsNo tags attached.

Activities

aatxutegi

2019-11-21 12:58

reporter  

DiwakarCoeffs.pdf (1,916,220 bytes)

will

2019-11-21 16:28

manager   ~0010932

Last edited: 2019-11-21 16:30

View 3 revisions

The code has the condition as "We > 0.5*sqrt(Re)". The paper has "We/sqrt(Re) > 0.5". They are the same.

I think line 118 should have "(d^3)/8" rather than "d". It's "r^3" in the paper, and it's not dimensionally consistent as it is.

Do you agree?

aatxutegi

2019-11-21 20:26

reporter   ~0010934

About the We condition you are right , my bad.

The other line d^(3/2) is obtained by d * sqrt(d). The original papers seems to set Cb to pi and Cs to 20 but these parameters assimilate the constant values by Cb_ = (Cb/4) = 0.785 and Cs_ = (Cs /2) = 10 which are the constants already under use. Sorry, you can close the issue.

will

2019-11-26 08:53

manager   ~0010939

Yes, you are right. I didn't spot the "d" outside the square root. It's all OK then.

Issue History

Date Modified Username Field Change
2019-11-21 12:58 aatxutegi New Issue
2019-11-21 12:58 aatxutegi File Added: DiwakarCoeffs.pdf
2019-11-21 16:28 will Note Added: 0010932
2019-11-21 16:30 will Note Edited: 0010932 View Revisions
2019-11-21 16:30 will Note Edited: 0010932 View Revisions
2019-11-21 20:26 aatxutegi Note Added: 0010934
2019-11-26 08:53 will Note Added: 0010939
2019-11-26 08:53 will Assigned To => will
2019-11-26 08:53 will Status new => closed
2019-11-26 08:53 will Resolution open => no change required