CFD Online Discussion Forums - InterFoam UEqn understanding

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)

- - InterFoam UEqn understanding (https://www.cfd-online.com/Forums/openfoam-programming-development/196031-interfoam-ueqn-understanding.html)

InterFoam UEqn understanding

Hi Dear Foamers,

I've some questions about UEqn.H interFoam.

Code:



    if (pimple.momentumPredictor())

    {

        solve

        (

            UEqn

         ==

            fvc::reconstruct

            (

                (

                    fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)

                  - ghf*fvc::snGrad(rho)

                  - fvc::snGrad(p_rgh)

                ) * mesh.magSf()

            )

        );

    }

1- Why right hand side is multiplied by area (mesh.magSf())?
2- Why every thing is in surface form (surface normal gradient, fvc::interpolate) and then reconstructed to cell values? couldn't every thing be in cell centered values from first place?

I'll be appreciated if someone make clear these to me.

Thanks in advance,

Behzad

Quote:

Originally Posted by gaza (Post 672637)

Hi Bezhad,
You can read this post

https://www.cfd-online.com/Forums/op...-pressure.html

Thanks Przemek. It was exactly what i wanted.

One resource I found helpful when looking into interFoam was the following:

Title: Description and utilization of interFoam multiphase solver

https://www.google.nl/url?sa=t&sourc...SuBIwv7yJMAI2P