[amirhoshangm]

I'm using ReactingMultiphaseParcel in my simulation. I found that in a particle either volume or density can change .

    // Update particle density or diameter
    if (td.cloud().constProps().constantVolume())
    {
        this->rho_ = mass1/this->volume();
    }
    else
    {
        this->d_ = cbrt(mass1/this->rho_*6.0/pi);
    }

This is very limiting. for example, in the case of reaction/drying in a particle, both density and size might change.
Is'n it better to calculate effective density based on mass fraction of different components (similar to CpEff) and then calculate particle diameter? In this way both density and size can change at the same time. Does anyone have any idea about this?

[amirhoshangm]

Some thing like this

    // Update particle density or diameter
    if (td.cloud().constProps().constantVolume())
    {
        this->rho_ = mass1/this->volume();
    }
    else
    {
        this->rho_ = rhoEff(td, pc, T0, idG, idL, idS);

        this->d_ = cbrt(mass1/this->rho_*6.0/pi);
    }