applications/solvers/solidMechanics/elasticAcpSolidFoam/waveCourantNo.H

   1 {
   2     scalar E0 = rheology.law().E()().internalField()[0];
   3     scalar nu0 = rheology.law().nu()().internalField()[0];
   4     scalar rho0 = rho.internalField()[0];
   5
   6     scalar waveVelocity =
   7         ::sqrt(E0*(1 - nu0)/(rho0*(1 + nu0)*(1 - 2*nu0)));
   8
   9     // Courant number
  10     scalarField Co =
  11         waveVelocity*runTime.deltaT().value()
  12        *mesh.surfaceInterpolation::deltaCoeffs().internalField();
  13
  14     Info<< "Courant Number mean: " << gAverage(Co)
  15         << " max: " << gMax(Co)
  16         << " wave velocity magnitude: " << waveVelocity << endl;
  17 }