[hisham007]

Hi foamers,

Im a new OP user.

Now I create a new source term (body force) named "f" and add to the N-S equation in icoFoam solver. This f represented by using scalar potential. I want to know how can I add this f? Any expert can help me?

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;

[hisham007]

Hi guyz,

I've successful add body force term but the result don't appear any effect from adding this body force? Anyone can help me?