[Alczem]

Hello everyone!


I am experimenting with icoReactingMultiphaseInterFoam to try and simulate the melting of a snowman made of ice surrounded by air, and with an incoming radiant flux from the top boundary to provide heat.


The issue is that as soon as the simulation starts, my snowman seems to collapse on itself, like a fluid would, even without any mass transfer. I think I defined the properties of the different phases properly and I can't figure out the issue. I tried to play with gravity and the simulation works fine without any gravity, so I must have defined the boundary conditions or the solid phase poorly.


Does someone have an idea? Thanks!




Here are some files that may be relevant


fvSolution.solvers

    "alpha.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 2;
        cAlphas         ((liquid and solid)  0
                    (liquid and gas)  1
                    (gas and solid)  0);

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          0;
    }

phaseProperties (without melting):

type    massTransferMultiphaseSystem;

phases    (liquid gas solid);

liquid
{
    type            pureMovingPhaseModel;
}

solid
{
    type            pureStaticSolidPhaseModel;
}

gas
{
    type            pureMovingPhaseModel;
}

thermophysicalProperties.solid

thermoType
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        nMoles          1;
        molWeight       18;
    }
    equationOfState
    {
        rho             1000;
    }
    thermodynamics
    {
        Cp              2093;
        Hf              0; //[J/Kg];
    }
    transport
    {
        kappa           2.22;
    }
}

[Krapf]

Could you please upload the whole case?