[sharonyue]

Hi guys,

In twoPhaseEulerFoam , we see:

while (pimple.loop())
        {
            fluid.solve();
            fluid.correct();

            volScalarField contErr1

I looked into the code, it returns like this:

void Foam::diameterModel::correct()
{}

So I think it calculates nothing....so why add this here? Im confused. Any ideas? Thanks,

[alexeym]

Hi,

If you take a look at description of the class

Class
    Foam::diameterModel
        
Description
    A2stract base-class for dispersed-phase particle diameter models.

you'll see that you should look at correct method in children of diameterModel class, such as IATE. In this class correct method is not so short.

[anuragm]

Hi Alexey,

I have a similar question about fluid.correct() and I do not follow your answer. Could you please explain it a little more?

Thanks

[ron_OFuser]

Anurag, You need to dig into that. The function is a member of the referred base class.

[alexeym]

Hi,

Originally Posted by anuragm

I have a similar question about fluid.correct() and I do not follow your answer. Could you please explain it a little more?

Well, you can start with ron_OFuser's answer.

More detailed version is: what is fluid? It is variable of type twoPhaseSystem defined in createFields.H:

twoPhaseSystem fluid(mesh, g);

correct method of twoPhaseSystem does this:

void Foam::twoPhaseSystem::correct()
{
    phase1_.correct();
    phase2_.correct();
}

what are phase1_ and phase2_? They are properties of twoPhaseSystem of type phaseModel:

class twoPhaseSystem
:
    public IOdictionary
{
...
        //- Phase model 1
        phaseModel phase1_;
...
}

What correct method of phaseModel does? It calls correct method of dPtr_

void Foam::phaseModel::correct()
{
    return dPtr_->correct();
}

and dPtr_ is a pointer to diameter model

class phaseModel
:
    public volScalarField,
    public transportModel
{
...
        //- Diameter model
        autoPtr<diameterModel> dPtr_;
...
}

Currently there are 3 diameter models in twoPhaseSystem:

constantDiameter
isothermalDiameter
IATE

First two models really have empty correct method. The last however has non-empty correct method.

Is it a little bit more clear this time?

[anuragm]

Many thanks for the detailed explanation. I really appreciate it.