From 963ef7e6887700c6df243d3e05b6e528bbf3875b Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 12 Apr 2011 12:00:39 +0100
Subject: [PATCH] ENH: Updated cloud splash model
---
.../ThermoSurfaceFilm/ThermoSurfaceFilm.C | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
index 7249846b138..541e9831556 100644
--- a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
+++ b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
@@ -391,8 +391,8 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
const scalar dBarSplash = 1/cbrt(6.0)*cbrt(mRatio/Ns)*d + ROOTVSMALL;
// cumulative diameter splash distribution
- const scalar dMax = cbrt(mRatio)*d;
- const scalar dMin = 0.001*dMax;
+ const scalar dMax = 0.9*cbrt(mRatio)*d;
+ const scalar dMin = 0.1*dMax;
const scalar K = exp(-dMin/dBarSplash) - exp(-dMax/dBarSplash);
// surface energy of secondary parcels [J]
@@ -437,7 +437,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
// magnitude of the normal velocity of the first splashed parcel
const scalar magUns0 =
- sqrt(2.0*parcelsPerSplash_*EKs/mSplash/(1 + coeff1/sqr(coeff2)));
+ sqrt(2.0*parcelsPerSplash_*EKs/mSplash/(1.0 + coeff1/sqr(coeff2)));
// Set splashed parcel properties
forAll(dNew, i)
@@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
// Apply correction to velocity for 2-D cases
meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U());
-
+Info<< "NEW PARTICLE: " << *pPtr << endl;
// Add the new parcel
this->owner().addParticle(pPtr);
--
GitLab