[bijan darbari]

Hello


I want to calculate new incompressible face flux "phi" based on a different velocity ( named "Ve" ) and use this flux in solving a special equation.


So I make a new file "CreatePhie.H" in the solver folder and called it in "createfields.H":

surfaceScalarField phie
(
    IOobject
    (
        "phie",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    linearInterpolate(Ve) & mesh.Sf()
 );

#include "createPhi.H"
#include "createPhie.H"

but the results are incorrect!!


for testing, I defined "phie" based on velociy "U":

linearInterpolate(U) & mesh.Sf()

but the results are much different from that solved by "phi".


What's the problem?? and what's the way for definition new flux "phie"??

[openfoam_aero]

Originally Posted by bijan darbari

Hello


I want to calculate new incompressible face flux "phi" based on a different velocity ( named "Ve" ) and use this flux in solving a special equation.


So I make a new file "CreatePhie.H" in the solver folder and called it in "createfields.H":

surfaceScalarField phie
(
    IOobject
    (
        "phie",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    linearInterpolate(Ve) & mesh.Sf()
 );

#include "createPhi.H"
#include "createPhie.H"

but the results are incorrect!!


for testing, I defined "phie" based on velociy "U":

linearInterpolate(U) & mesh.Sf()

but the results are much different from that solved by "phi".


What's the problem?? and what's the way for definition new flux "phie"??

Where have you created the createPhie.H file? It should be in the new solver directory which is different from the directory of the library (~/OpenFOAM/src/finiteVolume/lnInclude/). If you create your createPhie.H in this directory, then it will change the existing directory file createPhi.H which contains the flux calculation based on the velocity U.
So copy the file to the directory of your new solver, change the name and then compile.

Another thing is that you should not change the object identifier "phi" to something of your choice. The reason is that phi is used later in continuityErrors.H to calculate the local and global continuity errors. Keep phi as it is (and do not change it to "phie") in createPhie.H and then try recompiling. This way, both files createPhi.H and createPhie.H will be run by the pre-processor directive and the fluxes will be processed separately by the continuityErrors.H header file.