[immortality]

whats the mathematical equivalence of this formula?
which equation is used in rhoCentralFoam?

// --- Solve momentum
        solve(fvm::ddt(rhoU) + fvc::div(phiUp));

        U.dimensionedInternalField() =
            rhoU.dimensionedInternalField()
           /rho.dimensionedInternalField();
        U.correctBoundaryConditions();
        rhoU.boundaryField() = rho.boundaryField()*U.boundaryField();

        volScalarField rhoBydt(rho/runTime.deltaT());

        if (!inviscid)
        {
            solve
            (
                fvm::ddt(rho, U) - fvc::ddt(rho, U)
              - fvm::laplacian(muEff, U)
              - fvc::div(tauMC)
            );
            rhoU = rho*U;
        }

[wyldckat]

Hi Ehsan,

I don't know the answer, but I've picked on the following breadcrumbs on this topic:

• http://www.cfd-online.com/Forums/ope...algorithm.html
• http://www.cfd-online.com/Forums/ope...am-detail.html
• http://www.cfd-online.com/Forums/ope...ow-solver.html
• http://openfoamwiki.net/index.php/TestLucaG

Check the PDF files on the last link.

Best regards,
Bruno