[linyanx]

Dear OpenFOAMers,

I have a quick question regarding the volVectorField. How can we implement the volVectorField term into the equation?

Like, for example, UEqn.H in interMixingFoam. After defining a volVectorField in createField.H named 'F'. I want to add the force term(which is a vector term) into the UEqn.H, written as 'interpolate(F)'. But the system will always complain about the format of this vector term.

My UEqn.H is partially quoted below:

UEqn
==
fvc::reconstruct
(
(
mixture.surfaceTensionForce()
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
+ interpolate(F) ----The term I want to add.
) * mesh.magSf()
)

Hence, I look over all the source code in openfoam's solver package, but still have no idea of the solution. Really appreciate your help.

Regards,
Linyan

[linyanx]

Thanks to Fanny's kindly answer here. The volVectorField F is supposed to be ok when it was written at the right side of the UEqn. Hence, the code is corrected as follow:

PHP Code:

solve
        (
            UEqn
         ==
            fvc::reconstruct
            (
                (
                    mixture.surfaceTensionForce()
                  - ghf*fvc::snGrad(rho)
                  - fvc::snGrad(p_rgh)
                ) * mesh.magSf()
            )
//*               Modified on Mar/16            *//
            + F
//*               Modified on Mar/16            *//
        );