[MamboJambo]

Hello,

I would like to add a mass source term to simpleFOAM but do not know how to do it.

Can I just create a field called mass term:

volScalarField massSource
(
    IOobject
    (
        "massSource",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

And add it to pEqn.H line 28 as:

        fvScalarMatrix pEqn
        (
            fvm::laplacian(rAtU(), p) == fvc::div(phiHbyA) - massSource 
        );

Is this correct?

I was going to set the values of massSource from `setExprFields`.


Is this the correct procedure? Does the momentum equation need to changed?

[DuanYB]

Hello, I encountered the same problem as you, and I corrected it in the same way as you, but the program reported an error. Have you solved it now?