capillaryCourantNo.H
Go to the documentation of this file.
1 {
2  scalar CoNumSigma = max
3  (
4  sqrt
5  (
6  2*M_PI*sigma*sqr(aMesh.edgeInterpolation::deltaCoeffs())
7  *mag(aMesh.edgeInterpolation::deltaCoeffs())
8  /rhol
9  )
10  ).value()*runTime.deltaT().value();
11 
12  Info<< "Max Capillary Courant Number = " << CoNumSigma << '\n' << endl;
13 }
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
dimensionedSymmTensor sqr(const dimensionedVector &dv)
engineTime & runTime
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:487
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
dimensionedScalar rhol("rhol", dimDensity, transportProperties)
messageStream Info
Information stream (stdout output on master, null elsewhere)
faMesh aMesh(mesh)