During a core meltdown severe accident (SA) of light water reactors (LWRs), melt penetration in porous media (so-called melt infiltration) may occur due to melt relocation on top of a debris bed, or after dryout and re-melting of metallic particles containing steel and Zr, and then remelting of more refractory oxidic particles containing fuel and oxidized Zr in a hot zone of a debris bed. Both in-vessel and ex-vessel phases of SA progression may involve melt infiltration. For instance, if a debris bed in the lower head of reactor pressure vessel (RPV) is uncoolable, melt infiltration occurs due to corium debris remelting and can affect RPV failure mode and place. After the lower head failure, molten corium discharges from the RPV, leading to melt spreading in the reactor cavity if the cavity floor is dry or covered by a shallow water pool. In the case of deep water pool below the RPV, remelting of uncoolable debris bed can also result in melt spreading under the water. Clearly, modeling of melt infiltration and spreading are paramount not only to the prediction of severe accident progression, but also for the analyses of corium coolability and retention. In spite of their importance, insufficient work has been done in mechanistic modeling of the phenomena.

The objective of this doctoral thesis is to provide high-fidelity predictive capabilities for melt infiltration and spreading, which can be employed to substantiate the understanding of melt infiltration in the RPV and melt spreading in the reactor cavity. The focus of the work is to develop a computational code with multi-physics models for modeling specific phenomena important to melt infiltration and melt spreading. The models and code are based on the Moving Particle Semi-implicit (MPS) method, starting from two-dimensional representation and progressing with three-dimensional simulations.

The thesis first describes the MPS method which is a mesh free method suitable for free-surface flow. The governing equations including mass, momentum and energy conservation equations are discretized with particle interaction models, such as gradient model and Laplacian model. Furthermore, models for viscosity and phase change, surface tension and wettability are implemented in the code. Additional models of multi-phase flow, crust formation, and film boiling heat transfer are added for relevant applications.

The developed MPS code is then validated against various experiments for analyses of in-vessel melt infiltration, and ex-vessel melt spreading in dry and underwater conditions. In particular, three numerical studies have been conducted to investigate the capabilities of the MPS code for prediction of melt infiltration, melt spreading under dry condition and melt spreading under water. The key points from the numerical studies are as follows.  

• The MPS code is applied to predict melt infiltration phenomena in various particulate beds employed in the REMCOD experiments carried out at KTH-NPS with corium simulant materials. The wettability model is implemented, which is characterized by the contact angle between the melt and the debris bed. The REMCOD-E09-C4 and E09-C2 tests are calculated in which melt penetrates through hot debris beds with spherical particles and cylindrical particles, respectively.  Then, the MPS code is applied to simulate the REMCOD-E08-C4 test in which the debris bed temperature is below the melting point of melt, thereby the solidification occurs. The results suggest that the melt infiltration process in the experiment with cylindrical particles can be predicted by using Sauter mean diameter. Additionally, the simulation shows a good agreement with the tests.

• The MPS code is further extended and employed for the dry spreading test. A modified crust formation model has been proposed based on rigid body assumption to keep the dynamic shape of crust, which is appropriate for the crust formed at the top surface. Two tests conducted at KTH are simulated as one-dimensional and two-dimensional spreading schemes, respectively. These simulations concentrate on hydrodynamic motion and heat transfer while ignoring the associated phenomena such as MCCI and generation of gaseous concrete decomposition products. The results illustrate the capability of the MPS code predicting melt leading-edge progression and spread thickness in dry cavity.

• The MPS code capabilities are further extended to model the specific phenomena appearing during melt spreading over substrates under a water layer, such as multiphase flow and film boiling heat transfer. A smoothing scheme is proposed to replace the real properties of interfacial particles for multiphase flow. The film boiling heat transfer is calculated along the film boiling regime of the boiling curve. For validation, the updated MPS code is then employed to simulate the S3E-2MWS-Ox-1 and PULiMS-E9 tests, which are conducted at KTH. These simulations focus on investigating the thermal-hydraulic characteristics of 1D and 2D melt underwater spreading, and they predict the leading edge progression and terminal spread thickness well.

In general, the developed MPS code offers a novel approach to predict melt infiltration in the RPV and melt spreading in the reactor cavity with high fidelity. The insight from the simulations helps the understanding of corium coolability and retention.

Vid svåra haverier med härdsmälta i lättvattenreaktorer kan smältan penetrera porösa material (infiltration av smälta) då smältan förflyttas ner ovanpå fragmenterade härdrester eller efter torrkokning och återsmältning av metalliska partiklar med stål och zirkonium eller mer eldfasta oxidpartiklar med bränsle och oxiderad zirkonium i en het region av en bädd av fragmenterade härdrester. Infiltration av smälta kan uppstå både i och utanför reaktortryckkärlet under förloppet av ett svårt haveri. Om en bädd av härdrester i reaktortryckkärlets nedre plenum exempelvis är okylbar, uppstår infiltration av smälta på grund av återsmältning av härdrester som kan inverka på reaktortryckkärlets brottmekanism och brottställe. Efter brott rinner smältan ur reaktortryckkärlet och sprids ut på golvet i utrymmet under reaktortanken om golvet är torrt eller täckt av grunt vatten. Med en djup vattennivå under reaktortanken kan återsmältning av okylbara härdrester också resultera i att smältan sprids ut under vattnet. Det är uppenbart att modellering av smältinfiltration och smältspridning är avgörande inte bara för att förutsäga svåra haveriers utveckling, utan också för analyser av härdresternas kylbarhet och kvarhållning i reaktortryckkärlet. Trots betydelsen av infiltrering och spridning är det arbete som hittills genomförts inom mekanistisk modellering av infiltrationsfenomenen otillräckligt.

Syftet med denna doktorsavhandling är att utveckla verklighetsnära prognosverktyg för smältinfiltration och -spridning, vilka kan användas för att underbygga förståelsen av smältinfiltration i reaktortanken och i utrymmet under den. Fokus för arbetet är att utveckla en beräkningskod med multifysik-modeller för modellering av specifika fenomen som är viktiga för smältinfiltration och smältspridning. Modellerna och koden är baserade på metoden eng. Moving Particle Semi-implicit (MPS), med utgångspunkt från tvådimensionell representation och vidare utveckling till tredimensionella simuleringar.

Avhandlingen beskriver först MPS-metoden som är en nätfri (eng. meshfree) metod lämplig för flöde med fri yta (eng. free-surface flow). De tillämpade ekvationerna, inklusive ekvationer för massans, rörelsemängdens och energins bevarande, diskretiseras med partikelinteraktionsmodeller, såsom gradientmodell och Laplace-operator-modell (eng. Laplacian model). Vidare är modeller för viskositet och fasförändring, ytspänning och vätbarhet implementerade i koden. Ytterligare modeller läggs till för flerfasflöde, formation av skorpa och värmeöverföring genom filmkokning för relevanta tillämpningar.

Den utvecklade MPS-koden valideras sedan mot olika experiment för analyser av smältinfiltration i reaktortanken, och i utrymmet under reaktortanken under torra och våta förhållanden. I synnerhet har tre separata numeriska studier genomförts för att undersöka MPS-kodens möjligheter för prognoser av smältinfiltration, smältspridning under torrt tillstånd och smältspridning under vatten. Huvudpunkterna från de tre numeriska studierna är följande.

• MPS-koden tillämpas för att förutsäga smältinfiltrationsfenomen i olika partikelbäddar som använts i REMCOD-experimenten utförda vid KTH-NPS med simulerade härdmaterial. Vätbarhetsmodellen, som kännetecknas av kontaktvinkeln mellan smältan och bädden av fragmenterade härdrester, är implementerad. REMCOD-E09-C4- och E09-C2-testerna beräknas där smältan tränger igenom heta fragmentbäddar med sfäriska respektive cylindriska partiklar. Sedan tillämpas MPS-koden för att simulera REMCOD-E08-C4-testet där bäddens temperatur understiger smältans smältpunkt, varigenom smältan stelnar. Resultaten tyder på att smältinfiltrationsprocessen i experimentet med cylindriska partiklar kan förutsägas genom att använda Sauterdiametern. Dessutom visar simuleringen en god överensstämmelse med testerna.

• MPS-koden utökas ytterligare och tillämpas för testet med spridning på torr yta. En modifierad modell för bildning av skorpa har föreslagits baserat på ett stelkroppsantagande för att behålla skorpans dynamiska form, vilket är lämpligt för den skorpa som bildas på den övre ytan. Två tester utförda vid KTH simuleras som endimensionella respektive tvådimensionella spridningssystem. Dessa simuleringar koncentrerar sig på hydrodynamisk rörelse och värmeöverföring samtidigt som man bortser från tillhörande fenomen som reaktioner mellan smälta och betong (MCCI) och produktion av gasformiga betongnedbrytningsprodukter. Resultaten påvisar förmågan hos MPS-koden att förutsäga smältfrontens utbredning och spridningstjocklek i ett torrt utrymme.

• MPS-kodfunktionerna utökas ytterligare för att modellera de specifika fenomen som uppstår under smältspridning över substrat under ett vattenskikt, såsom flerfasflöde och värmeöverföring vid filmkokning. En utjämningsmetod (eng. smoothing scheme) föreslås för att ersätta de verkliga egenskaperna hos partiklar i gränsytan för flerfasflöde. Värmeöverföring genom filmkokning beräknas längs kokkurvans filmkokningsregion. För validering används sedan den uppdaterade MPS-koden för att simulera testerna S3E-2MWS-Ox-1 och PULiMS-E9, som genomförts på KTH. Dessa simuleringar fokuserar på att undersöka de termohydrauliska karaktäristika hos smältspridning under vatten i 1D och 2D, och de förutsäger smältfrontens utbredning och den slutliga tjockleken hos smältan väl.

Generellt sett erbjuder den utvecklade MPS-koden ett nytt sätt att förutsäga smältinfiltration och smältspridning i reaktortryckkärlet respektive utrymmet under reaktortanken med hög tillförlitlighet. Insikten från simuleringarna bidrar till förståelsen av härdsmältans kylbarhet och kvarhållning i reaktortryckkärlet.