[aliyah.]

Hello

I want to modify the thermoPhysicalModels library in a way that by changing alpha, viscosity change.
The problem is alpha should be read from mesh and mesh should be defined in defining mu.

The original mu definition that I want to change is located in

openfoam2006/src/thermophysicalModels/specie/transport/polynomial/polynomialTransportI.H

as follow:

template<class Thermo, int PolySize>
inline Foam::scalar Foam::polynomialTransport<Thermo, PolySize>::mu
(
    const scalar p,
    const scalar T
) const
{
    return muCoeffs_.value(T);
}

So I wan to change it to something like

template<class Thermo, int PolySize>
inline Foam::scalar Foam::polynomialTransport<Thermo, PolySize>::mu
(
    const scalar p,
    const scalar T
) const
{

const fvMesh& mesh;
const volScalarField& alphaL_ = mesh.lookupObject<volScalarField>("alpha.liquid");
scalar liquidVolume = fvc::domainIntegrate(alphaL_).value(); 

    return muCoeffs_.value(T)/liquidVolume;
}

It will be compiled, but the mesh must be initialized; without initialization, the solver will not work when we attempt to use it.

I looked at some other codes in openfoam and initializing mesh can be done by a code like

const fvMesh& mesh = alpha1_.mesh();

I tried

const fvMesh& mesh = this->mesh();

it is not working because this does not have a member named mesh, or .

const fvMesh& mesh = T.mesh();

which obviously wrong because T is scalar.

So any suggestion how I can initialize mesh here?