[chris Liu]

Hi, foamers:

I am a new OpenFOAM user, and i want to add a source term (body force) to icoFoam.

I creat a new field named "F" and add F to the N-S equation. After complie, i test it in the lid-driven-cavity example, and use the force on the lid to drive the fluid instead of the lid velocity. However, the fluid stay still and the force fill the cavity.(I want the force on the lid only and the fluid is driven.)

Any suggestions will be appreciated!
Best regards!
chris Liu

createFields.H
// add source term;
Info<< "Reading field F\n" << endl;
volVectorField F
(
IOobject
(
"F",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

icoFoam-force.C
fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)

);

solve(UEqn == F -fvc::grad(p)); //add source term;


Detailed source file in the attachment.