View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0002682 | OpenFOAM | Bug | public | 2017-08-31 12:53 | 2017-09-05 14:03 |
Reporter | shildenbrand | Assigned To | will | ||
Priority | high | Severity | minor | Reproducibility | always |
Status | resolved | Resolution | fixed | ||
Platform | HP Z640 | OS | Debian Linux | OS Version | 9 |
Product Version | dev | ||||
Summary | 0002682: Overprediction of energy of secondary droplets in BaiGosman spray-wall impingement model | ||||
Description | In src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C there is a mistake in the computation of the energy of the secondary droplets after splash: Right now, the incident kinetic energy is computed as const scalar EKIn = 0.5*m*magSqr(Urel); which should be const scalar EKIn = 0.5*m*magSqr(Un); Please see original article by Bai and Gosman (eqn. 8) where it is clearly stated that the wall-normal velocity has to be used. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.31.6230&rep=rep1&type=pdf | ||||
Tags | Lagrangian, wallInteraction | ||||
|
BaiGosman.patch (453 bytes)
--- src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C 2017-08-31 13:20:55.297970475 +0200 +++ /tmp/ThermoSurfaceFilm.C 2017-08-31 13:33:47.229357887 +0200 @@ -400,7 +400,7 @@ } // Incident kinetic energy [J] - const scalar EKIn = 0.5*m*magSqr(Urel); + const scalar EKIn = 0.5*m*magSqr(Un); // Incident surface energy [J] const scalar ESigmaIn = np*sigma*p.areaS(d); |
|
I cannot access the article |
|
Thanks for the report, and the patch. Resolved in dev by commit d31bd913. |
Date Modified | Username | Field | Change |
---|---|---|---|
2017-08-31 12:53 | shildenbrand | New Issue | |
2017-08-31 12:53 | shildenbrand | Tag Attached: Lagrangian | |
2017-08-31 12:53 | shildenbrand | Tag Attached: wallInteraction | |
2017-08-31 13:03 | shildenbrand | File Added: BaiGosman.patch | |
2017-09-05 09:27 | will | Note Added: 0008696 | |
2017-09-05 14:03 | will | Assigned To | => will |
2017-09-05 14:03 | will | Status | new => resolved |
2017-09-05 14:03 | will | Resolution | open => fixed |
2017-09-05 14:03 | will | Note Added: 0008697 |