applications/utilities/preProcessing/molConfig/velocityDistributions/uniform.H

   1 scalar initVelMag =
   2     sqrt
   3     (
   4         3.0*(1.0 - 1.0 / totalZoneMols)
   5        *moleculeCloud::kb*temperature
   6        /mass
   7     );
   8
   9 for (molN = totalMols; molN < totalMols + totalZoneMols; molN++)
  10 {
  11     // Assign velocity: random direction, magnitude determined by desired
  12     // temperature
  13
  14     // Temperature gradients could be created by specifying a gradient in the
  15     // zone subDict, or by reading a field from a mesh.
  16
  17     // The velocities are treated on a zone-by-zone basis for the purposes of
  18     // removal of bulk momentum - hence nMols becomes totalZoneMols
  19
  20     velocity = (2.0*rand.vector01() - vector::one);
  21
  22     velocity *= initVelMag/mag(velocity);
  23
  24     momentumSum += mass*velocity;
  25
  26     initialVelocities.append(velocity);
  27 }