[abhiv]

Hello OFoamers, I'm relatively new to OpenFOAM community. I'm facing an issue while defining boundary condition to one of the variables in a custom solver using openfoam.

In particular, I'm calculating Flux of a quantity, which is function of two other independent variables. But I need to provide boundary conditions for the flux but not the independent variables. Independent variables may have "calculated" BC.

Here is a part of createFields.H, showing one possible situation. "Flux" requires a boundary condition to solve the physics correctly, but the solver is not reading boundaryField while running. I think that it occurs because I'm not providing 'mesh' condition instead of equation for "Flux" in IOobject. The "Flux" is also being used in another PDE solution.

    volVectorField Flux
    (
        IOobject
        (
            "Flux",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        N*U
    );

Is there a solution already exists in OpenFOAM for this problem? Please advice me if I am interpreting the openfoam guidelines incorrectly for createFields.H. How to approach if I have to come up with an utility to perform this job?

Thanks in advance,
ABHI

[Sugajen]

Hi Abhi,

I have a similar situation. Did you figure out a way?

Thanks