View Issue Details

IDProjectCategoryView StatusLast Update
0002682OpenFOAMBugpublic2017-09-05 14:03
Reportershildenbrand Assigned Towill  
PriorityhighSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
PlatformHP Z640OSDebian LinuxOS Version9
Product Versiondev 
Summary0002682: Overprediction of energy of secondary droplets in BaiGosman spray-wall impingement model
DescriptionIn 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

TagsLagrangian, wallInteraction

Activities

shildenbrand

2017-08-31 13:03

reporter  

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);
BaiGosman.patch (453 bytes)   

will

2017-09-05 09:27

manager   ~0008696

I cannot access the article

will

2017-09-05 14:03

manager   ~0008697

Thanks for the report, and the patch. Resolved in dev by commit d31bd913.

Issue History

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