[AHRabiee]

Hi
I am very new in OpenFOAM.
I want to change icoFoam solver in order to solve a special equation.
my equation is very similar to Navier Stokes equation with small differences. one of these differences is in convection term.
in my equation, convection term has an additional parameter which is a volScalarField , therefore I must implement div(U*U*eps) instead of div(U*U) -------> eps = volScalarField

but I can't write it in OpenFOAM language with no syntax errors

please help me

thanks!

[Mehdi3031]

Originally Posted by AHRabiee

Hi
I am very new in OpenFOAM.
I want to change icoFoam solver in order to solve a special equation.
my equation is very similar to Navier Stokes equation with small differences. one of these differences is in convection term.
in my equation, convection term has an additional parameter which is a volScalarField , therefore I must implement div(U*U*eps) instead of div(U*U) -------> eps = volScalarField

but I can't write it in OpenFOAM language with no syntax errors

please help me

thanks!

what you want is exactly like fvm::div(rhoPhi, U) where rho is a volScalar ,
have a look at the UEqn.H and creatField.H files in interFoam solver in this directory:
applications/solvers/multiphase/interFoam

and follow the same path.