[purnp2]

Hello,
I am using one RAS model (Hellsten-EARSM from Furst's GIT repository) to make some changes into it.
In the transport equation maker function, I want to change a few coefficients (zeta_, beta_, and alphaD_) from its type [dimensionedScalar] to another type [volScalarField] to make it compatible with the rest of the code in its original code.
Original Code:

        volScalarField gamma( this->gamma(fMix) );
        volScalarField beta( this->beta(fMix) );
        volScalarField alphaD( this->alphaD(fMix) );

Here, this->gamma(fMix) is of the type volScalarField.
Changes I made which are producing error:

        volScalarField zeta( this->zeta_ );
        volScalarField beta( this->beta_ );
        volScalarField alphaD= this->alphaD_ ;

In the changed code, zeta_ [dimensionedScalar] replaces this->gamma(fMix) [volScalarField] and like wise in two other cases.
I am attaching changed files also:
So, how can I change its type?

Thank you for your time for reading this.
Purn

[purnp2]

One possible solution (which, I guess, is not correct even though it compiles and run) is as follows:
I used a formula:

    tmp<volScalarField> blend
    (
        const volScalarField& fMix,
        const dimensionedScalar& psi
    ) const
    {
        return fMix*(psi)/fMix;
    }

were fMix is of type [volScalarFeild].
Is it an acceptable solution?

I am hoping for a solution where I can use any standard function of foam:: like:

typedef GeometricField<scalar, fvPatchField, volMesh> volScalarField

[jherb]

Have a look e.g. at https://github.com/OpenFOAM/OpenFOAM...ewtonian.C#L54

     nu_
     (
         IOobject
         (
             name,
             U_.time().timeName(),
             U_.db(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
         U_.mesh(),
         nu0_
     )

nu_ is a volScalarField and nu0_ a dimensioned scalar. So if you use this in a function (and not in a constructor), it should look like this (not tested):

     volScalarField nu_
     (
         IOobject
         (
             "name",
             U_.time().timeName(),
             U_.db(),
             IOobject::NO_READ,
             IOobject::NO_WRITE
         ),
         U_.mesh(),
         nu0_
     )

you need the variable containing a reference to the mesh or the internal database in for the call to IOobject. (see e.g. the solver source codes or boundary conditions)

[purnp2]

Thank you very much, sir. I'll check this code and update this reply if it works for me or not. I appreciate your help.