zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Updated cloud splash model

Browse Source
This commit is contained in:
andy
2011-04-12 12:00:39 +01:00
parent 1291caa0dc
commit 963ef7e688
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C
View File

@ -391,8 +391,8 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
const scalar dBarSplash = 1/cbrt(6.0)*cbrt(mRatio/Ns)*d + ROOTVSMALL; const scalar dBarSplash = 1/cbrt(6.0)*cbrt(mRatio/Ns)*d + ROOTVSMALL;
// cumulative diameter splash distribution // cumulative diameter splash distribution
const scalar dMax = cbrt(mRatio)*d; const scalar dMax = 0.9*cbrt(mRatio)*d;
const scalar dMin = 0.001*dMax; const scalar dMin = 0.1*dMax;
const scalar K = exp(-dMin/dBarSplash) - exp(-dMax/dBarSplash); const scalar K = exp(-dMin/dBarSplash) - exp(-dMax/dBarSplash);
// surface energy of secondary parcels [J] // 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 // magnitude of the normal velocity of the first splashed parcel
const scalar magUns0 = 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 // Set splashed parcel properties
forAll(dNew, i) forAll(dNew, i)
@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
// Apply correction to velocity for 2-D cases // Apply correction to velocity for 2-D cases
meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U()); meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U());
Info<< "NEW PARTICLE: " << *pPtr << endl;
// Add the new parcel // Add the new parcel
this->owner().addParticle(pPtr); this->owner().addParticle(pPtr);