A novel multiphase MPS algorithm for modeling crust formation by highly viscous fluid for simulating corium spreading

• Published in: Nuclear engineering and design 2019-03, Vol.343, p.218-231
• Main Authors: Duan, Guangtao, Yamaji, Akifumi, Koshizuka, Seiichi
• Format: Article
• Language: English
• Subjects: Accidents; Algorithms; Computer simulation; Corium spreading; Creep (materials); Crust formation; Dam failure; Fluid-solid interaction; Free surfaces; Lava; Mathematical models; MPS method; Nuclear accidents & safety; Nuclear energy; Nuclear power plants; Numerical creeping; Numerical methods; Particle method; Particle methods (mathematics); Solidification; Spreading; Velocity; Viscosity; Viscous fluids
• ISSN: 0029-5493; 1872-759X
• DOI: 10.1016/j.nucengdes.2019.01.005

•A novel algorithm calculating viscosity after pressure is proposed for MPS.•Numerical creeping of highly viscous fluid is eliminated by the new algorithm.•Viscosity escalation can reasonably represent crust formation with the new algorithm.•VULCANO VE-U7 experiment demonstrates improvement of the new algorithm. Corium (lava-like mixture of fissile material) spreading prediction is of great significance in the severe accidents of nuclear power plants. Crust formation due to solidification distinguishes corium spreading from common isothermal spreading. The Lagrangian moving particle semi-implicit (MPS) method is potential for such spreading flow with both free surface and crust-melt interface. Crust formation is usually represented by viscosity escalation, but crust creeping is an associated problem. In the original MPS algorithm, creeping velocity cannot be reduced steadily by the continuous increase of viscosity, owing to the numerical creeping. A new solution algorithm is proposed for particle methods to eliminate such numerical creeping, so that creeping velocity decreases proportionally with viscosity rise. In this situation, high enough viscosity can effectively represent crust behaviors. Three numerical examples, leakage flow with high viscosity, dam break flow with low viscosity and the VULCANO VE-U7 corium spreading experiment with both high and low viscosities simultaneously, are investigated to contrast the performance difference between the original and new algorithms. It is demonstrated that the current algorithm is suitable for crust formation in corium spreading.