[gaza]

Hi All,
I coded function

tmp<volVectorField> fun()
{
    return tmp<volVectorField> 
             (
                  fvc::div(rho1/rho2*B*B*G*G)
             );
};

and when I run the solver I got this message like

keyword div((((rho1|rho2)*B*B*G*G) is undefined in fvSchemes

how to avoid this?

[guin]

Try following in Solver:

fvc::div(rho1/rho2*B*B*G*G,"weirdDivergenceTerm")

and in fvSchemes:

weirdDivergenceTerm Gauss upwind; //as an example

[gaza]

Hi Rodrigo,
Thank you so much, it works