[Andreas M.]

Hi,

i ran into an issue using fvOptions and since I haven't found any solution yet:

I want to use fvOptions in order to add a precalculated sourceterm to my UEqn.

Despite the examples i found yet, fvOptions should contain a vector/vectorfield, which are calculated in advance for the whole interior region.

Examplecode:

tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
      + turbulence->divDevReff(U)
      ==
        fvOptions(U)
    );

where "fvOptions(U)" is a (most likely) different vector for every cell.
Just adding the vector/vectorfield into the equation results into a Foam Fatal Error.

Im glad for every response.

Thanks in advance.


Im sorry for grammatical errors.

[Zeppo]

fvVectorMatrix UEqnLHS
(
    fvm::div(phi, U)
 + turbulence->divDevReff(U)
);

IOobject io
(
    "USource",
    U.instance(),
    U.mesh(),
    NO_READ,
    NO_WRITE,
    false
);

DimensionedVectorField USource
(
    io,
    U.mesh(),
    UEqnLHS.dimensions()/dimVolume;
);

forAll(USource, i)
{
    USource[i] = vector(i, 2*i, 3*i);
}

fvVectorMatrix UEqn
(
    UEqnLHS 
 == USource
);