Experimental and Numerical Investigation of Solidification of Gallium in an Initially Emptied Horizontal Pipe Flow

Author(s):  
Shawn Somers-Neal ◽  
Alex Pegarkov ◽  
Edgar Matida ◽  
Vinh Tang ◽  
Tarik Kaya

Abstract In a reactor core meltdown under postulated severe accidents, the molten material (corium) could be ejected or relocated through existing vessel penetrations (cooling pipe connections), thus potentially contaminating other locations in the power plant. There exists, however, a potential for plugging of melt flow due to its complete solidification, providing the availability of an adequate heat sink. Therefore, a numerical model was created to simulate the flow of molten metal through an initially empty horizontal pipe. The numerical model was verified using a previously developed analytical model and validated against experimental tests with gallium (low melting temperature) as a substitute for corium. The numerical model was able to predict the penetration length (length of distance travelled by the molten metal) after a complete blockage occurred with an average percent error range of 9%. Since the numerical model has been verified and validated, the model can be updated to predict the penetration length in the cooling pipe in case of a severe accident.

Author(s):  
Rida S. N. Mahmudah ◽  
Masahiro Kumabe ◽  
Takahito Suzuki ◽  
LianCheng Guo ◽  
Koji Morita ◽  
...  

Understanding the freezing behavior of molten metal in flow channels is of importance for severe accident analysis of liquid metal reactors. In order to simulate its fundamental behavior, a 3D fluid dynamics code was developed using Finite Volume Particle (FVP) method, which is one of the moving particle methods. This method, which is fully Lagrangian particle method, assumes that each moving particle occupies certain volume. The governing equations that determine the phase change process are solved by discretizing its gradient and Laplacian terms with the moving particles. The motions of each particle and heat transfer between particles are calculated through interaction with its neighboring particles. A series of experiments for fundamental freezing behavior of molten metal during penetration on to a metal structure was also performed to provide data for the validation of the developed code. The comparison between simulation and experimental results indicates that the present 3D code using the FVP method can successfully reproduce the observed freezing process such as molten metal temperature profile, frozen molten metal shape and its penetration length on the metal structure.


Author(s):  
Osamu Kawabata ◽  
Masao Ogino

When the primary reactor system remain pressurized during core meltdown for a typical PWR plant, loop seals formed in the primary reactor system would lead to natural circulations in hot leg and steam generator. In this case, the hot gas released from the reactor core moves to a steam generator, and a steam generator tube would be failed with cumulative creep damage. From such phenomena, a high-pressure scenario during core meltdown may lead to large release of fission products to the environment. In the present study, natural circulation and creep damage in the primary reactor system accompanying the hot gas generation in the reactor core were discussed and the combining analysis with MELCOR and FLUENT codes were performed to examine the natural circulation behavior. For a typical 4 loop PWR plant, MELCOR code which can analyze for the severe accident progression was applied to the accident analyses from accident initiation to reactor vessel failure for the accident sequence of the main steam pipe break which is maintained at high pressure during core meltdown. In addition, using the CFD code FLUENT, fluid dynamics in the reactor vessel plenum, hot leg and steam generator of one loop were simulated with three-dimensional coordinates. And the hot gas natural circulation flow and the heat transfer to adjoining structures were analyzed using results provided by the MELCOR code as boundary conditions. The both ratios of the natural circulation flow calculated in the hot leg and the steam generator using MELCOR code and FLUENT code were obtained to be about 2 (two). And using analytical results of thermal hydraulic analysis with both codes, creep damage analysis at hottest temperature points of steam generator tube and hot leg were carried out. The results in both cases showed that a steam generator tube would be failed with creep rupture earlier than that of hot leg rupture.


Author(s):  
Juanhua Zhang ◽  
Jiming Lin ◽  
Shishun Zhang

Reactor Pit Flooding System (RPF) is adopted under the severe accidents situation in CPR1000+ units. It can move the heat generated from the reactor core via external reactor vessel cooling (ERVC) to keep the integrity of RPV and achieve the in-vessel corium retention (IVR). But if IVR function of RPF is failed, there is Ex-Vessel Steam Explosion (EX-SE) risk. The Ex-Vessel Steam Explosion is analyzed by MC3D software which is for fuel and cooling interaction (FCI). The physical model of CPR1000+ for Steam Explosion is built firstly and then the phenomenon of Ex-Vessel Steam Explosion under typical severe accident is analyzed. The conclusion of this study is that the impulse load of pressure on the cavity wall induced by steam explosion is about 310KPas ∼ 440KPas. Referencing the structure capacity of AP600 containment, if the structural capacity of CPR1000+ containment is equal to AP600, the impulse load of pressure is lower than it. So it could be preliminarily estimated that steam explosion will not threaten the integrality of CPR1000+ containment.


Author(s):  
Peiqi Liu ◽  
Tao Yu ◽  
Hongyan Yang

A typical 1000MW pressurized-water reactor (PWR) unit model of China’s living nuclear power plant (NPP) units is built based on MAAP4[1] in this paper. Different severe accidents cases caused by different LOCA area on hot leg of primary loop are studied. And different mitigation measures are focused to evaluated their effectiveness. The study indicates that during the accident, the larger broken area LOCA case caused the more severe rector core damaged. However, it is important to inject water into the reactor core in good time. And that can mitigate the severe accident progress effectively.


Author(s):  
Sunil Nijhawan ◽  
YongMann Song

Abstract As analysts still grapple with understanding core damage accident progression at Three Mile Island and Fukushima that caught the nuclear industry off-guard once too many times, one notices the very limited detail with which the large reactor cores of these subject reactors have been modelled in their severe accident simulation code packages. At the same time, modelling of CANDU severe accidents have largely borrowed from and suffered from the limitations of the same LWR codes (see IAEA TECDOC 1727) whose applications to PHWRs have poorly caught critical PHWR design specifics and vulnerabilities. As a result, accident management measures that have been instituted at CANDU PHWRs, while meeting the important industry objective of publically seeming to be doing something about lessons learnt from say Fukushima and showing that the reactor designs are oh so close to perfect and the off-site consequences of severe accidents happily benign. Integrated PHWR severe accident progression and consequence assessment code ROSHNI can make a significant contribution to actual, practical understanding of severe accident progression in CANDU PHWRs, improving significantly on the other PHWR specific computer codes developed three decades ago when modeling decisions were constrained by limited computing power and poor understanding of and interest in severe core damage accidents. These codes force gross simplifications in reactor core modelling and do not adequately represent all the right CANDU core details, materials, fluids, vessels or phenomena. But they produce results that are familiar and palatable. They do, however to their credit, also excel in their computational speed, largely because they model and compute so little and with such un-necessary simplifications. ROSHNI sheds most previous modelling simplifications and represents each of the 380 channels, 4560 bundle, 37 elements in four concentric ring, Zircaloy clad fuel geometry, materials and fluids more faithfully in a 2000 MW(Th) CANDU6 reactor. It can be used easily for other PHWRs with different number of fuel channels and bundles per each channel. Each of horizontal PHWR reactor channels with all their bundles, fuel rings, sheaths, appendages, end fittings and feeders are modelled and in detail that reflects large across core differences. While other codes model at best a few hundred core fuel entities, thermo-chemical transient behaviour of about 73,000 different fuel channel entities within the core is considered by ROSHNI simultaneously along with other 15,000 or so other flow path segments. At each location all known thermo-chemical and hydraulic phenomena are computed. With such detail, ROSHNI is able to provide information on their progressive and parallel thermo-chemical contribution to accident progression and a more realistic fission product release source term that would belie the miniscule one (100 TBq of Cs-137 or 0.15% of core inventory) used by EMOs now in Canada on recommendation of our national regulator CNSC. ROSHNI has an advanced, more CANDU specific consideration of each bundle transitioning to a solid debris behaviour in the Calandria vessel without reverting to a simplified molten corium formulation that happily ignores interaction of debris with vessel welds, further vessel failures and energetic interactions. The code is able to follow behaviour of each fuel bundle following its disassembly from the fuel channel and thus demonstrate that the gross assumption of a core collapse made in some analyses is wrong and misleading. It is able to thus demonstrate that PHWR core disassembly is not only gradual, it will be also be incomplete with a large number of low power, peripheral fuel channels never disassembling under most credible scenarios. The code is designed to grow into and use its voluminous results in a severe accident simulator for operator training. It’s phenomenological models are able to examine design inadequacies / issues that affect accident progression and several simple to implement design improvements that have a profound effect on results. For example, an early pressure boundary failure due to inadequacy of heat sinks in a station blackout scenario can be examined along with the effect of improved and adequate over pressure protection. A best effort code such as ROSHNI can be instrumental in identifying the risk reduction benefits of undertaking certain design, operational and accidental management improvements for PHWRs, with some of the multi-unit ones handicapped by poor pressurizer placement and leaky containments with vulnerable materials, poor overpressure protection, ad-hoc mitigation measures and limited instrumentation common to all CANDUs. Case in point is the PSA supported design and installed number of Hydrogen recombiners that are neither for the right gas (designed mysteriously for H2 instead of D2) or its potential release quantity (they are sparse and will cause explosions). The paper presents ROSHNI results of simulations of a postulated station blackout scenario and sheds a light on the challenges ahead in minimizing risk from operation of these otherwise unique power reactors.


1996 ◽  
Vol 118 (1) ◽  
pp. 164-172 ◽  
Author(s):  
C. H. Amon ◽  
K. S. Schmaltz ◽  
R. Merz ◽  
F. B. Prinz

A molten metal droplet landing and bonding to a solid substrate is investigated with combined analytical, numerical, and experimental techniques. This research supports a novel, thermal spray shape deposition process, referred to as microcasting, capable of rapidly manufacturing near netshape, steel objects. Metallurgical bonding between the impacting droplet and the previous deposition layer improves the strength and material property continuity between the layers, producing high-quality metal objects. A thorough understanding of the interface heat transfer process is needed to optimize the microcast object properties by minimizing the impacting droplet temperature necessary for superficial substrate remelting, while controlling substrate and deposit material cooling rates, remelt depths, and residual thermal stresses. A mixed Lagrangian–Eulerian numerical model is developed to calculate substrate remelting and temperature histories for investigating the required deposition temperatures and the effect of operating conditions on remelting. Experimental and analytical approaches are used to determine initial conditions for the numerical simulations, to verify the numerical accuracy, and to identify the resultant microstructures. Numerical results indicate that droplet to substrate conduction is the dominant heat transfer mode during remelting and solidification. Furthermore, a highly time-dependent heat transfer coefficient at the droplet/substrate interface necessitates a combined numerical model of the droplet and substrate for accurate predictions of the substrate remelting. The remelting depth and cooling rate numerical results are also verified by optical metallography, and compare well with both the analytical solution for the initial deposition period and the temperature measurements during droplet solidification.


Sign in / Sign up

Export Citation Format

Share Document