|
[Sponsors] |
October 3, 2020, 13:29 |
Melting using icoReactingMultiPhaseInterFoam
|
#1 |
Member
Join Date: Mar 2019
Posts: 81
Rep Power: 7 |
Dear Foamers,
I am trying to simulate the melting process of a solid body (phase change material) in 3D using icoReactingMultiPhaseInterFoam (first question: is this the best solver for melting? ). I included the buoyant melted movement by Boussinesq as the equation of state for the liquid. The problem is that when I run the simulation, it goes fine until towards the end of melting when suddenly the Courant number goes higher and higher until the alpha value becomes negative... I tried this in 2D and it never happened. But in 3D despite using a relatively fine mesh this happens... I also tried with adjustableRunTime which went down to 10^-6... I also checked the mesh quality: Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Time = 0 Mesh stats points: 206283 faces: 2290188 internal faces: 2234676 cells: 1131216 faces per cell: 4 boundary patches: 3 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 0 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 1131216 polyhedra: 0 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology Symmetry 19091 9950 ok (non-closed singly connected) Tube 18819 9679 ok (non-closed singly connected) Wall 17602 9007 ok (non-closed singly connected) Checking faceZone topology for multiply connected surfaces... No faceZones found. Checking basic cellZone addressing... No cellZones found. Checking geometry... Overall domain bounding box (0 0 0) (0.175 0.175 0.5) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (1.30135e-17 -3.18322e-16 2.85069e-16) OK. Max cell openness = 2.33158e-16 OK. Max aspect ratio = 4.71177 OK. Minimum face area = 6.46877e-07. Maximum face area = 4.8044e-05. Face area magnitudes OK. Min volume = 2.68775e-10. Max volume = 8.59991e-08. Total volume = 0.01183. Cell volumes OK. Mesh non-orthogonality Max: 52.2945 average: 14.4077 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.623961 OK. Coupled point location match (average 0) OK. Mesh OK. End Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1812 | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object phaseProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // type massTransferMultiphaseSystem; phases (solid liquid); liquid { type pureMovingPhaseModel; } solid { type pureStaticSolidPhaseModel; } interfacePorous ( (solid and liquid) { type VollerPrakash; solidPhase alpha.solid; Cu 1e5; } ); massTransferModel ( (solid to liquid) { type Lee; C 40; Tactivate 358.15; } ); // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1812 | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application icoReactingMultiphaseInterFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 4000; deltaT 0.01; writeControl runTime; writeInterval 5; purgeWrite 0; writeFormat binary; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable no; adjustTimeStep no; //maxDeltaT 0.01; //maxCo 1.; //maxAlphaCo 1.; //maxAlphaDdt 2.3; //maxDi 25; // ************************************************************************* // Thanks, MJ PS: Now I am trying to further refine the mesh, hoping it is possible to do so... Last edited by mm66; October 4, 2020 at 10:19. |
|
Tags |
alpha, melting, negative, openfoam 1812, phase change material |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Stefan Problem : Tin Melting : Moving Boundary | swparth | OpenFOAM | 0 | September 21, 2017 14:53 |
Melting and Solidification | Mansoor_shad | FLUENT | 0 | April 23, 2017 13:37 |
melting in horizontal shell in tube, need help friends | thermal energy | FLUENT | 0 | January 9, 2014 17:33 |
melting of a pcm | flex00 | FLUENT | 2 | April 8, 2013 13:33 |
Simple Melting Problem | flex00 | FLUENT | 0 | March 15, 2013 05:10 |