View Issue Details

IDProjectCategoryView StatusLast Update
0000565OpenFOAMBugpublic2012-06-28 16:30
Reporteruser411Assigned Touser2 
PriorityhighSeveritymajorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version10.04
Summary0000565: tMom calculation in calcBreakup function in SprayParcel.C
Descriptionline 248 in SprayParcel.C:
scalar tMom = 1.0/(Fcp.Sp() + Fncp.Sp());
tMom is the reciprocal of implicit force coefficient.

line 125 in ReitzKHRT.C:
vector acceleration = Urel/tMom;

In fact Urel/tMom is not acceleration. It is force. It has to be divided
by drolet mass to get acceleration. So the tMom should be

tMom = 1.0/(Fcp.Sp() + Fncp.Sp) * mass;

same in line 305, child->tMom() should be
child->tMom() = 1.0/(Fcp.Sp() + Fncp.Sp()) * massChild;

Thank you.
TagsNo tags attached.

Activities

user2

2012-06-28 16:30

  ~0001430

Thanks for the report - bug corrected by commit 3ee36d

Issue History

Date Modified Username Field Change
2012-06-27 06:16 user411 New Issue
2012-06-28 16:30 user2 Note Added: 0001430
2012-06-28 16:30 user2 Status new => resolved
2012-06-28 16:30 user2 Fixed in Version => 2.1.x
2012-06-28 16:30 user2 Resolution open => fixed
2012-06-28 16:30 user2 Assigned To => user2