[anuragm]

Hello everyone,

I need to implement a custom diameter model in the multiphaseEulerFoam solver. For this purpose, I modified the isothermalDiameter model in the following fashion.

The following variables were declared in the header file:

    // Private data

        //- Reference diameter of the phase
        volScalarField d_;

        //- Initial number of particles
        dimensionedScalar No_;

        //- Initial global dispersed void fraction
        dimensionedScalar Vo_;

        //- Number assigned to the phase
        dimensionedScalar pNum_;

These variables were initialized in the constructor as follows:

    diameterModel(dict, phase),
    d_
    (
        IOobject
        (
            IOobject::groupName("d", phase.name()),
            phase.time().timeName(),
            phase.db(),
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        phase.mesh()
    ),
    No_("No", dimless/dimVol, dict_.lookup("No")),
    Vo_("Vo", dimless, dict_.lookup("Vo")),
    pNum_("pNum", dimless, dict_.lookup("pNum"))

And finally, the function d() was modified in the following fashion:

Foam::tmp<Foam::volScalarField> Foam::diameterModels::DQMOM::d() const
{
    scalar oneDivThree = 1.0/3.0;
    const dimensionedScalar& CC 
      = Foam::pow( 6.0 * Vo_ / constant::mathematical::pi / No_, oneDivThree );
            
    label n = pNum_.value();
    word fieldName = "ABSC" + Foam::name(n);
    const volScalarField& absc 
      = phase_.db().lookupObject<volScalarField>(fieldName);

    return CC * Foam::pow(absc, oneDivThree);
}

The compilation completed without errors but when I select this model for the dispersed phase and execute my simulation, I get no change in the diameter for any of the dispersed phases.

[anuragm]

I forgot to add that the volScalarfield variables ABSC1/ABSC2/ABSC3 are non-uniform and change with time (verified by my simulation results). However, the diameter field does not get changed at all despite the fact that the value returned by the function d() depends on ABSC1/ABSC2/ABSC3