Numerical performance and sensitivity of the thermo-hygro-corrosive model of the French underground radioactive waste storage tunnels
In recent years, a multiphysics model that represents the complex physical phenomena that influence the accumulation of rust in the storage tunnels (alveoli) has been numerically implemented in the finite element method (FEM) software Cast3M. Seeking to estimate the thermo-hygro-corrosive properties of the alveoli at long time scales, recent (and ongoing) works on improving the execution time of the resolution algorithm were undertaken. However, for a FEM numerical model to be considered as a rigorous engineering/scientific tool, error bars must be associated to all computational results; therefore, the careful quantification of the plethora of modeling uncertainties is primordial. To undertake such an endeavor, multiple issues must first be tackled, begininning with an improvement of the physical representativity of the multiphysics model, following with an improvement of the computational performance of the numerical model, and ending with rigorous sensitivity studies of the implemented model. Work on computational performance is necessary, so as to render program execution fast enough to ensure that the large sets of numerical experiments required run in reasonable times.
Study of the seismic behavior of piping systems using mechanical models of different degrees of fidelity
Piping systems are part of the equipment to which particular attention is paid as part of the safety review or design of nuclear installations. They are designed in accordance with codes, standards and regulations to withstand loads that occur or may occur over the life of a facility. These systems must therefore be designed to withstand accidental loads such as earthquakes. Feedback shows that piping systems generally behave well in the event of an earthquake. When failures are observed, they are more likely to be due to significant anchor movement, brittle materials, unwelded joints, corrosion, piping support failures, or seismic interactions. In practice, to be able to estimate the beyond design seismic behavior and the associated failure risks, the engineer can implement numerical models involving varying degrees of refinement depending on needs. This study consists of taking stock of the numerical modeling capabilities of piping systems under earthquake. For reasons of computational burden, global modeling based on beam elements is often favored, considering simplified material laws such as bilinear material laws with kinematic hardening. We know the “theoretical” limits of these models but it is difficult to have clear ideas about their real limits of applicability depending on the level of loading and the damage targeted. To make this assessment, we propose to interpret, using different numerical models involving different degrees of fidelity, the results of the experimental campaign carried out by the BARC and which was used for the MECOS benchmark (METallic COmponent margins under high Seismic loads).
Development, metrological validation and outdoor testing of a multitrack Raman/FO measurement unit dedicated to the safety of future cryogenic liquid hydrogen dispensing stations
Context: The domestic and industrial use of liquid hydrogen as the fuel of the future requires the definition of a suitable safety code. At present, tank separation criteria have been defined by anticipation using a conservative approach. It is therefore necessary to carry out full-scale experiments ("pool spreading") in order to provide input for calculation codes and build relevant standards. These experiments require the implementation of instrumentation adapted to the measurement of all gases present in free space (O2, N2, H2O, H2) in order to establish a measurement of partial pressures during each test, correlated with the other means of measurement in place (thermometry, catharometry, PIV, BOS, etc.).
Mission: In the context of an ANR-PEPR project (ESKHYMO) managed by CEA Liten, a Raman/FO Multitrack spectrometric measurement unit will be developed jointly by CEA List and CEA DES on the basis of an existing device. Raman measurement is multi-elemental, multi-track (a single measurement unit for several probes), non-explosive, and delivers a self-standardized measurement to a reference species (usually nitrogen at atmospheric pressure). The Raman/FO measurement unit comprises a laser, a spectrometer associated with a scientific CCD camera, and a fiber-optic circuit for remote measurement. The design of the Raman/FO probes will also be based on an existing CEA product, which will be miniaturized for deployment in field conditions. Four Raman/FO probes will be produced and then calibrated in air (climatic chamber) and hydrogen (shock tube or vacuum chamber) at CEA DES DM2S in Saclay. Finally, the final device will be deployed on the test site for multi-gas measurements during spraying experiments, in partnership with Air Liquide and accrediting bodies (INERIS).
Skills: Optics, laser, fiber optics, spectrometry
Seismic behavior of an overhead crane
Overhead cranes are part of the equipment in industrial installations to which special attention must be paid. They are generally located in the upper part of buildings and are potentially subject to significant levels of acceleration in the event of an earthquake, due to the amplification induced by the supporting structure. Consequently, they are potentially subjected to significant forces and can be the source of significant forces on the supporting structure. This study is a continuation of two previous test campaigns carried out on the Azalée shaking table of the EMSI laboratory, on a mock-up of an overhead crane. It aims to provide validated numerical models of this kind of equipment. Two lines of research are considered. The first axis aims to complement the “historical” test campaigns with static tests to justify the adjustment of the numerical models. The second axis consists of exploiting, by comparison tests/calculations, all of the tests that were carried out as part of a previous test campaign for statistical analysis purposes.
Earthquake effect on underground facilities
The Industrial Centre for Geological Disposal (Cigeo) is a project for a deep geological disposal facility for radioactive waste to be built in France. These wastes will be put in sealed packages in tunnels designed at 500 meters depth. The seals are made of a bentonite/sand mixture which has a high swelling capacity and a low water permeability. As a part of the long-term safety demonstration of the repository, it must be demonstrated that the sealing structures can fulfill their functions under seismic loads over their entire lifetime. In order to guarantee this future nuclear waste repository, CEA and Andra are collaborating to work on the potential scientific and engineering challenges involved.
The responses of underground repository to earthquake events are complex due to the spatially and temporally evolving hydro-mechanical properties of the surrounding media and the structure itself. Accurate modeling of the behavior, therefore, requires a coupled multiphysics numerical code to efficiently model the seismic responses for these underground repositories within their estimated lifespan of 100 thousand years.
The research will therefore, propose a performance assessment for sequential and parallel finite element numerical modeling for earthquake analysis of deep underground facilities. Then perform a synthetic data sampling to account for material uncertainties and based on the obtained results in the previous assessment, run a sensitivity analysis using a FEM or a metamodeling process. Finally, the results and knowledge gained within the span of this project will be processed and interpreted to provide responses for industrial needs.
Natural convection at high Ra numbers for nuclear safety: 2nd year
Thermal exchanges at very high Rayleigh numbers (Ra) exist on geophysical scale, at civil engineering scale and increasingly in industrial applications and here particularly in the energy sector. At this point, we mention the cooling of solar panels or the heat removal from nuclear power plants under accidental conditions. In fact, the passive safety concept of Small Modular Reactors (SMR) is based on the transfer of residual heat from the reactor to a water pool in which the reactor is placed. Since the outer reactor vessel is very high, heat exchange occurs by natural convection at Rayleigh numbers (Ra) between 1010 and 1016. Reliable heat transfer correlations exist to date only up to about Ra < 1012 with very high uncertainties in the extrapolation to higher Ra. Understanding the heat transfer at very high Ra is thus of fundamental and practical interest. The associated challenges are twofold:
• Numerical challenges: CFD codes cannot model turbulent heat transfer at very high Ra with sufficient accuracy and appropriate calculation time. Improved physical and numerical models are required, which use high performance computing (HPC) capabilities.
• Experimental challenges: Detailed experiments are essential for code validation. Since experiments in water require impractical huge dimensions, cryogenic experiments with helium are planned at CEA, based on the interesting physical properties of this fluid in the range of 5 K (high thermal expansion associated to low viscosity and thermal conduction).
Development and application of Inverse Uncertainty Quantification methods in thermal-hydraulics within the new OECD/NEA activity ATRIUM
Within the Best Estimate Plus Uncertainty methodologies (BEPU) for the safety analysis of the Nuclear Power Plants (NPPs), one of the crucial issue is to quantify the input uncertainties associated to the physical models in the code. Such a quantification consists of assessing the probability distribution of the input parameters needed for the uncertainty propagation through a comparison between simulations and experimental data. It is usually referred to as Inverse Uncertainty Quantification (IUQ).
In this framework, the Service of Thermal-hydraulics and Fluid dynamics (STMF) at CEA-Saclay has proposed a new international project within the OECD/NEA WGAMA working group. It is called ATRIUM (Application Tests for Realization of Inverse Uncertainty quantification and validation Methodologies in thermal-hydraulics). Its main objectives are to perform a benchmark on relevant Inverse Uncertainty Quantification (IUQ) exercises, to prove the applicability of the SAPIUM guideline and to promote best practices for IUQ in thermal-hydraulics. It is proposed to quantify the uncertainties associated to some physical phenomena relevant during a Loss Of Coolant Accident (LOCA) in a nuclear reactor. Two main IUQ exercises with increasing complexity are planned. The first one is about the critical flow at the break and the second one is related to the post-CHF heat transfer phenomena. A particular attention will be dedicated to the evaluation of the adequacy of the experimental databases for extrapolation to the study of a LOCA in a full-scale reactor. Finally, the obtained input model uncertainties will be propagated on a suitable Integral Effect Test (IET) to validate their application in experiments at a larger scale and possibly justify the extrapolation to reactor scale.
Experimentation and numerical simulation of lithium battery thermal runaway
In the current Energy transition context, the lithium battery is an essential technology to address the strong challenge of the electrical energy storage. However, Li battery severe solicitations/loadings can lead to a thermal runaway phenomenon, which can cause an outbreak of fire, even an explosive combustion of the cell or of the whole battery pack. If this phenomenon is well known, the research and development dedicated to the battery safety is emerging and must be consolidated. The post-doctorate global objective is to develop a numerical modelling and simulation strategy for thermal runaway occurring when a Li battery is subjected to mechanical/thermal/electrical abuse, in order to gain an understanding of the phenomenon, estimate the thermal spreading risk as a result of gas combustion, or study the runaway mechanical consequences (fluid structure interaction). This strategy relies on physical testing campaigns carried out as part of the post-doctorate, and on numerical tools developed by CEA (EUROPLEXUS, Cast3M). The work will be organised into three main content areas: Understanding and modelling of the phenomena on the basis of experimental tests (shock tube, abusive tests), Development of a numerical model representative of identified phenomena, Modelling including fluid-structure interaction (case deformation due to pressure increase).