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

GradientDispersionRAS: Updated to use the new GaussNormal function

Browse Source 
Provided by Timo Niemi
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1650
This commit is contained in:
Henry
2015-04-20 17:48:16 +01:00
parent ded8cf365d
commit 971e5a67d5
1 changed files with 6 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
src/lagrangian/turbulence/submodels/Kinematic/DispersionModel/GradientDispersionRAS/GradientDispersionRAS.C
View File

@ -124,33 +124,22 @@ Foam::vector Foam::GradientDispersionRAS<CloudType>::update
{ {
tTurb = 0.0; tTurb = 0.0;
scalar sigma = sqrt(2.0*k/3.0); const scalar sigma = sqrt(2.0*k/3.0);
vector dir = -gradk/(mag(gradk) + SMALL); const vector dir = -gradk/(mag(gradk) + SMALL);
// Numerical Recipes... Ch. 7. Random Numbers... scalar fac = 0.0;
scalar x1 = 0.0;
scalar x2 = 0.0;
scalar rsq = 10.0;
while ((rsq > 1.0) || (rsq == 0.0))
{
x1 = 2.0*rnd.sample01<scalar>() - 1.0;
x2 = 2.0*rnd.sample01<scalar>() - 1.0;
rsq = x1*x1 + x2*x2;
}
scalar fac = sqrt(-2.0*log(rsq)/rsq);
// In 2D calculations the -grad(k) is always // In 2D calculations the -grad(k) is always
// away from the axis of symmetry // away from the axis of symmetry
// This creates a 'hole' in the spray and to // This creates a 'hole' in the spray and to
// prevent this we let x1 be both negative/positive // prevent this we let fac be both negative/positive
if (this->owner().mesh().nSolutionD() == 2) if (this->owner().mesh().nSolutionD() == 2)
{ {
fac *= x1; fac = rnd.GaussNormal<scalar>();
} }
else else
{ {
fac *= mag(x1); fac = mag(rnd.GaussNormal<scalar>());
} }
UTurb = sigma*fac*dir; UTurb = sigma*fac*dir;