[Dark Helmet]

Hi,

I'm new to OpenFOAM, so I'm guessing this is simple but could use really some help. I'd like to add a momentum source via fvOptions to an interFoam simulation. I think that SemiImplicitSource is what I should be using, along the lines of this:

bodyForce
{
    type vectorSemiImplicitSource;
    active true;
    selectionMode cellZone;
    cellZone heatStrip;

    vectorSemiImplicitSourceCoeffs
    {
        volumeMode specific;
        injectionRateSuSp
        {
            U ((fx fy fz) 0);
        }
    }
}

Hopefully that is the correct syntax. Anyway, the "tricky" part is I want the source to be a function of position, fx = F0*(z-z0)+y etc. Does anyone know how I could do this?

[pbachant]

You may need to write your own source fvOption, or use a codedSource, putting your code in the codeAddSup entry. You'll have access to the mesh object, and therefore cell centers (mesh_.C()), from which you can calculate your source term field to add to eqn.source().

See https://github.com/OpenFOAM/OpenFOAM.../CodedSource.H for the source code and dictionary syntax.