[shipman]

Dear All,

Maybe it is a bit boring question but I just wondered that how the code read the alpha1 phaseName as alpha.water in the compressibleInterFoam solver?

When I looked at the code, I could not see anywhere alpha1 phaseName assigned as alpha.water.

First I looked at createFields where alpha1 is seen:

Info<< "Constructing twoPhaseMixtureThermo\n" << endl;
twoPhaseMixtureThermo mixture(mesh); 

volScalarField& alpha1(mixture.alpha1());
volScalarField& alpha2(mixture.alpha2());

then looked at twoPhaseMixtureThermo.H which in inherited twoPhaseMixute:

class twoPhaseMixtureThermo
:
    public psiThermo,
    public twoPhaseMixture

When I checked the twoPhaseMixture.H file:

class twoPhaseMixture
{
protected:

    // Protected data

        word phase1Name_;
        word phase2Name_;

        volScalarField alpha1_;
        volScalarField alpha2_;

Finally Phase names are defined in twoPhaseMixture.C file as follows:

 phase1Name_(wordList(dict.lookup("phases"))[0]),
    phase2Name_(wordList(dict.lookup("phases"))[1]),

    alpha1_
    (
        IOobject
        (
            IOobject::groupName("alpha", phase1Name_),
            mesh.time().timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    ),

    alpha2_
    (
        IOobject
        (
            IOobject::groupName("alpha", phase2Name_),
            mesh.time().timeName(),
            mesh
        ),
        1.0 - alpha1_
    )

So, my question is how the alpha.water is assigned to name of alpha1?

I may overlook something.

Thank you

[shipman]

Ok sorry for boring question. I found it finally. It reads from const/thermophyPro dic.

[cyw]

Originally Posted by shipman

Ok sorry for boring question. I found it finally. It reads from const/thermophyPro dic.

Hello, where did you find it? I cannot find the dict in the code.