[JFM]

Good day to All
I would like to confirm with the community that the standard utilities in OpenFOAM 2.3

• R
• yPlusRAS
• wallShearStress

do not work with interFoam as they are not designed to interpret multiphase model results. When I have tried to initiate these utilities I receive the following error message.

Reading/calculating face flux field phi



--> FOAM FATAL IO ERROR: 
keyword transportModel is undefined in dictionary "/media/TERRA/MEng/MEng2014/OpenFOAM/working/5-k-e_60-90s/1448kMesh_ke_60-90s.org/constant/transportProperties"

file: /media/TERRA/MEng/MEng2014/OpenFOAM/working/5-k-e_60-90s/1448kMesh_ke_60-90s.org/constant/transportProperties from line 21 to line 73.

    From function dictionary::lookupEntry(const word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 437.

FOAM exiting

After reading a number of posts on the forum I have modified the transportProperties dictionary by including an additional single phase transportModel (i.e. water in this instance), as below, and the utilities will then output the respective output for the 'dummy' phase.

    transportModel  Newtonian;
    // Kinematic viscosity
    nu              nu [ 0 2 -1 0 0 0 0 ] 1e-06;
    // Density
    rho             rho [ 1 -3 0 0 0 0 0 ] 1000;

phases (water air);

water
{
    transportModel  Newtonian;
    nu              nu [ 0 2 -1 0 0 0 0 ] 1e-06;
    rho             rho [ 1 -3 0 0 0 0 0 ] 1000;
    CrossPowerLawCoeffs
    {
        nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
        nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        m               m [ 0 0 1 0 0 0 0 ] 1;
        n               n [ 0 0 0 0 0 0 0 ] 0;
    }

    BirdCarreauCoeffs
    {
        nu0             nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
        nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        k               k [ 0 0 1 0 0 0 0 ] 99.6;
        n               n [ 0 0 0 0 0 0 0 ] 0.1003;
    }
}

air
{
    transportModel  Newtonian;
    nu              nu [ 0 2 -1 0 0 0 0 ] 1.48e-05;
    rho             rho [ 1 -3 0 0 0 0 0 ] 1;
    CrossPowerLawCoeffs
    {
        nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
        nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        m               m [ 0 0 1 0 0 0 0 ] 1;
        n               n [ 0 0 0 0 0 0 0 ] 0;
    }

    BirdCarreauCoeffs
    {
        nu0             nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;
        nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
        k               k [ 0 0 1 0 0 0 0 ] 99.6;
        n               n [ 0 0 0 0 0 0 0 ] 0.1003;
    }
}

At this point in time I cannot code (but will be learning over the next few months) however this leads me to ask several questions:

• How accurate are the results
• What happens at the interface
• Is there a different utility that can be used in OF2.3
• Has someone created their own utilities for compilation into OpenFOAM

Any feedback will be greatly appreciated.