[praveen]

In order to do steady state calculations with rhoCentralFoam, I am trying to use localEuler ddt scheme. From what I understood, this requires specifying rDeltaT as a volScalarField. I have made following additions to rhoCentralFoam in version 2.0.0

createField.H:

volScalarField rDeltaT
(
    IOobject
    (
        "rDeltaT",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    1/dimensionedScalar("maxDeltaT", dimTime, GREAT),
    zeroGradientFvPatchScalarField::typeName
);

setrDeltaT.H:

{       
        
    // Set the reciprocal time-step from the local Courant number
    rDeltaT.dimensionedInternalField() = max
    (   
        1/dimensionedScalar("maxDeltaT", dimTime, maxDeltaT),
        fvc::surfaceSum(amaxSf)().dimensionedInternalField()
       /(maxCo*mesh.V())
    );  

    // Update tho boundary values of the reciprocal time-step
    rDeltaT.correctBoundaryConditions();
        
}

rhoCentralFoam.C:

        surfaceScalarField aphiv_pos(phiv_pos - aSf);
        surfaceScalarField aphiv_neg(phiv_neg + aSf);

        // Reuse amaxSf for the maximum positive and negative fluxes
        // estimated by the central scheme
        amaxSf = max(mag(aphiv_pos), mag(aphiv_neg));

        #include "compressibleCourantNo.H"
        #include "readTimeControls.H"
        //#include "setDeltaT.H"
        #include "setrDeltaT.H"

        runTime++;

        Info<< "Time = " << runTime.timeName() << nl << endl;

Then in system/fvSchemes, I specify

ddtSchemes
{
    default         localEuler  rDeltaT;
}

The code compiles and runs but blows up after some iterations. With the usual rhoCentralFoam it runs fine, but it takes forever to converge. I use a maxCo of 0.1 and it is a subsonic flow. Limiters are also used (Gamma scheme). I could not find anything on localEuler in the forums so I want to know if anybody has succeeded in using it, or if you can find any mistake in my modifications.