setRDeltaT.H

Go to the documentation of this file.

{
    volScalarField& rDeltaT = trDeltaT.ref();
 
    scalar rDeltaTSmoothingCoeff
    (
        runTime.controlDict().getOrDefault<scalar>
        (
            "rDeltaTSmoothingCoeff",
            0.02
        )
    );
 
    // Set the reciprocal time-step from the local Courant number
    rDeltaT.ref() = max
    (
        1/dimensionedScalar("maxDeltaT", dimTime, maxDeltaT),
        fvc::surfaceSum(amaxSf)()()
       /((2*maxCo)*mesh.V())
    );
 
    // Update tho boundary values of the reciprocal time-step
    rDeltaT.correctBoundaryConditions();
 
    fvc::smooth(rDeltaT, rDeltaTSmoothingCoeff);
 
    {
        auto limits = gMinMax(rDeltaT.primitiveField());
        limits.reset(1/(limits.max()+VSMALL), 1/(limits.min()+VSMALL));
 
        Info<< "Flow time scale min/max = "
            << limits.min() << ", " << limits.max() << endl;
    }
}