Understanding of two-phase flow through porous medium with intense phase change is of interest in many situations, including nuclear, chemical or geophysical applications. Intense boiling occurs when the liquid is injected into a highly superheated medium. Under such conditions, the heat flux extracted by the fluid from the porous medium is mainly governed by the nucleation of bubbles and by the evaporation of thin liquid films. Both configurations are possible, depending on local flow conditions and on the ratio of bubble size to pore size. The present study is motivated by the safety evaluation of light water nuclear reactors in case of a severe accident scenario, such as the one that happened in Fukushima Dai-ichi plant in March, 2011. If water sources are not available for a long period of time, the reactor core heats up due to the residual power and eventually becomes significantly damaged due to intense oxidation of metals and fragmentation of fuel rods resulting in the formation of a porous medium where the particles have a characteristic length-scale of 1 to 5 mm The coolability of the porous medium will depend on the water flow rate which can enter the medium under the available driving head and on the geometrical features of the porous matrix (average pore size, porosity). Therefore, it is of high interest to evaluate the conditions for which the injection of water in such porous medium is likely to stop the progression of the accident. The present paper addresses the issue of modelling two-phase flow and heat transfers in a porous medium initially dry, where water is injected. The medium is initially at a temperature well above the saturation temperature of water. In a first part, a summary of existing knowledge is provided, showing the scarcity of models and experimental data. In a second part, new experimental results obtained in an IRSN facility are analysed. The experiment consists in a bed of steel particles that are heated up to 700 degrees Cbefore injecting water. The facility is briefly described. The velocity of the "quench front" (location where particles are quickly cooled down) and the total pressure drop across the medium are estimated. The dependencies of those quantities with respect to the inlet water flow rate, the initial temperature of the medium and the diameter of particles are obtained. In a third part, a model is proposed, based on a previously developed model which is improved in order to take into account intense boiling regimes (in particular nucleate boiling). The model includes a function that takes into account the contact area between water and the particles which depends on the temperature of particles and on the void fraction. That function affects the local intensity of phase change. The model involves a few parameters which cannot be evaluated analytically. Those parameters are bounded, following the analysis of experimental data. Finally, the model is assessed by comparison of calculations with those new experimental data. The satisfactory agreement shows that the model is almost predictive in the range of parameters studied. The experimental results also show that the quench front becomes unstable under certain conditions. This is also analysed and compared with the predictions of the model.
