[Ahyar]

I tried to create a polydisperse simulation in openfoam using multiphaseEulerFoam solver. As far as I know, the difference between the solid and gas phase in multiphase simulation on openfoam lies in the thermoPhysicalProperties file in constant directory. so I tried to resemble the thermoPhysicalProperties of my solid phase into the thermoPhysicalProperties of titaniaSynthesis case located at $FOAM_TUTORIALS/multiphase/multiphaseEulerFoam/laminar/titaniaSynthesis.
But my simulation went divergent when I apply the EquationOfState to rhoConst and the equationOfState.rho value into a number that bigger than water density.



here's the code from my thermoPhysicalProperties

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      thermophysicalProperties.particle;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture; // awalnya puremixture
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        molWeight  58.6934;
    }
    equationOfState
    {
        rho         1009;
    }
    thermodynamics
    {
        Cp          684;
        Hf          -1.18307e+07;
    }
    transport
    {
        mu          0;
        Pr          8.0e-6;
    }
}

// ************************************************************************* //

and here's the one from titaniaSynthesis tutorial case

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties.particles;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            heRhoThermo;
    mixture         multiComponentMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

species (TiO2);

defaultSpecie TiO2;

TiO2
{
    specie
    {
        molWeight   79.87880;
    }
    equationOfState
    {
        rho         4230;
    }
    thermodynamics
    {
        Hf          -1.18307e+07;
        Cp          684.246;
    }
    transport
    {
        mu          1e-7;
        Pr          8.0e-6;
    }
}


// ************************************************************************* //

and this is the error





what should I change to make the simulation runs well?