CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Solid collapses in icoReactingMultiphaseInterFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   October 26, 2022, 04:42
Default Solid collapses in icoReactingMultiphaseInterFoam
Join Date: Dec 2021
Posts: 89
Rep Power: 3
Alczem is on a distinguished road
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


        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);

    type            pureMovingPhaseModel;

    type            pureStaticSolidPhaseModel;

    type            pureMovingPhaseModel;
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;

        nMoles          1;
        molWeight       18;
        rho             1000;
        Cp              2093;
        Hf              0; //[J/Kg];
        kappa           2.22;
Attached Images
File Type: jpg snowman_0s.JPG (21.6 KB, 5 views)
File Type: jpg snowman_1s.JPG (26.7 KB, 5 views)
File Type: jpg snowman_2s.JPG (29.6 KB, 6 views)
Alczem is offline   Reply With Quote

Old   October 26, 2022, 13:05
Senior Member
Join Date: Oct 2017
Posts: 111
Rep Power: 7
Krapf is on a distinguished road
Could you please upload the whole case?
Krapf is offline   Reply With Quote


ice, icoreacting, melting, multiphase

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Solid combustion problem.... godspeedsam FLUENT 3 May 24, 2022 04:29
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 05:15
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
Moving solid problem Mudblood FLUENT 0 June 10, 2010 12:11
How to calculate density of solid phase zhou Main CFD Forum 0 December 17, 1999 19:06

All times are GMT -4. The time now is 12:30.