Simulation of reactive gas-liquid multi-phase flows
The objective of this postdoctoral position is to develop and implement a simulation method for the simulation of a
sodium spray fire. Two key points need to be adressed. First, one needs to propose a proper representation of the sodium
droplets (dispersed phase) from their generation by a jet (separate phase) fragmentation to their behavior (motion,
oxidation, combustion) in the air atmosphere. This requires to derive a flow model that accounts for multiple components
with multiple interface topology regimes (dispersed and separate). Second, one needs to develop a robust discretization
strategy for this complex flow model.
The numerical work will be implemented in a new numerical tool to perform simulations of sodium spray fires developed at CEA. This tool is based on the canoP. Canop is a library designed for solving computational fluid dynamics problems using a cell-based
Adaptive Mesh Refinement (AMR) approach and parallel calculation.
Advanced modeling of thermal turbulent flows
For several decades, numerical simulations in fluid mechanics have significantly contributed to the design and maintenance of industrial installations. Turbulence modeling, a key area at the intersection of research and industry, has seen substantial advancements in both LES and RANS approaches. Since the early 2010s, hybrid methods that combine RANS and LES techniques have emerged to leverage the advantages of each, necessitating proficiency in both modeling types. The TrioCFD code developed at STMF, although capable of handling these models, has not seen adequate investment in modern approaches. To incorporate hybrid models, it is essential to update and enhance the current models. The proposed task is to identify the most relevant models for industrial applications, restructure the software to accommodate these models, and validate their performance.
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.