[dsh5400]

I'm working with the v2f RAS model. Ultimately, I'd like to be able to modify the model's coefficients without having to recompile.

The two files, v2f.C and v2f.H, found in the directory, /src/turbulenceModels/incompressible/RAS/v2f/, seem to contain everything I need to make the changes.

v2f.C contains a section called "Constructors" and looks like this

// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

v2f::v2f
(
    const volVectorField& U,
    const surfaceScalarField& phi,
    transportModel& transport,
    const word& turbulenceModelName,
    const word& modelName
)
:
    RASModel(modelName, U, phi, transport, turbulenceModelName),

    Cmu_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "Cmu",
            coeffDict_,
            0.22
        )
    ),
    CmuKEps_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "CmuKEps",
            coeffDict_,
            0.09
        )
    ),
...

Which leads me to believe I should just modify these values as needed. However, the OpenFOAM User Guide states:

7.2.1 Model coefficients
The coefficients for the RAS turbulence models are given default values in their respective source code. If the user wishes to override these default values, then they can do so by adding a sub-dictionary entry to the RAS sub-dictionary file, whose keyword name is that of the model with Coeffs appended, e.g. kEpsilonCoeffs for the kEpsilon model. If the printCoeffs switch is on in the RAS sub-dictionary, an example of the relevant …Coeffs dictionary is printed to standard output when the model is created at the beginning of a run. The user can simply copy this into the RAS sub-dictionary file and edit the entries as required.

And the v2f.H file contains the following in the header:

    \verbatim
        v2fCoeffs
        {
            Cmu         0.22;
            CmuKEps     0.09;
            C1          1.4;
            C2          0.3;
            CL          0.23;
            Ceta        70;
            Ceps2       1.9;
            sigmaEps    1.3;
            sigmaK      1;
        }
    \endverbatim

Which has me wondering if this is the "sub-directory entry" mentioned by the User Guide. Is something that I should move into the body of the code, and modify the values?

So, I guess my real question is: what is the easiest way to modify these coefficients? And, does this require that I compile the RAS directory (using the Allwmake found in the /src/turbulenceModels/incompressible/ directory, parallel to the /RAS/ directory)?

Thanks!

[akidess]

You don't need to change the code at all to modify these values. Just add the v2fCoeffs into the file where you select the turbulence model. Note that a dictionary never refers to a code file, but a configuration file.

[dsh5400]

Originally Posted by akidess

You don't need to change the code at all to modify these values. Just add the v2fCoeffs into the file where you select the turbulence model. Note that a dictionary never refers to a code file, but a configuration file.

Thank you! I added the code to the RASProperties file found in the /constant/ directory as follows:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RASModel        v2f;

v2fCoeffs
{
    Cmu         0.221;
    CmuKEps     0.091;
    C1          1.41;
    C2          0.31;
    CL          0.231;
    Ceta        701;
    Ceps2       1.91;
    sigmaEps    1.31;
    sigmaK      1.1;
}

turbulence      on;

printCoeffs     on;


// ************************************************************************* //

[yeya]

Just wondering but in which cases would you need to modify them?

[dsh5400]

Originally Posted by yeya

Just wondering but in which cases would you need to modify them?

Just looking to see if there is some room to improve the values of the constants in various RANS models. Working with flows that are transitioning in/out of turbulent states, where RANS models can have some trouble predicting the physics of the flow.

[yeya]

sounds good thanks!