Study of the Velocity-Vorticity-Pression formulation for discretising the Navier-Stokes equations.
The incompressible Navier-Stokes equations are among the most widely used models to describe the flow of a Newtonian fluid (i.e. a fluid whose viscosity is independent of the external forces applied to the fluid). These equations model the fluid's velocity field and pressure field. The first of the two equations is none other than Newton's law, while the second derives from the conservation of mass in the case of an incompressible fluid (the divergence of velocity vanishes). The numerical approximation of these equations is a real challenge because of their three-dimensional and unsteady nature, the vanishing divergence constraint and the non-linearity of the convection term. Various discretisation methods exist, but for most of them, the mass conservation equation is not satisfied exactly. An alternative is to introduce the vorticity of the fluid as an additional unknown, equal to the curl of the velocity. The Navier-Stokes equations are then rewritten with three equations. The post-doc involves studying this formulation from a theoretical and numerical point of view and proposing an efficient algorithm for solving it, in the TrioCFD code.
Modeling of the MADISON fuel irradiation device for the future JHR reactor
The Jules Horowitz Reactor (RJH), currently under construction at CEA's Cadarache site, will irradiate materials and fuels in support of the French and international nuclear industry, as well as producing radioelements for medical use. To carry out its missions, the reactor will be equipped with numerous experimental devices. In particular, the MADISON device, currently under design, will irradiate 2 or 4 fuel samples under nominal stationary or operational transient conditions. The loop is representative of light-water reactor operating conditions, with single-phase and two-phase forced convection.
The main objective of the Post-Doc is to model the MADISON device and all associated heat exchanges precisely, in order to help determine the overall heat balance during the test and thus improve the accuracy of the linear power imposed on the samples. To this end, a coupled thermal model (describing the fuel rods and device structures) / CFD thermal-hydraulic model (describing the coolant) will be established using the NEPTUNE_CFD/SYRTHES code. The modeling will be validated based on results obtained from similar modeling carried out on the ISABELLE-1 and ADELINE single-rod devices in the OSIRIS and RJH reactors. The proposed approach fits in with the logic of developing digital twins of the RJH experimental devices.
Cascade of circulicity in compressible turbulence
In this post-doctorate, we propose to study the properties of the small scales of forced compressible homogeneous turbulence. More precisely, exact statistical relations similar to the Monin-Yaglom relation will be investigated. The idea, detailed in reference [1], is to understand how the transfer of circulicity is organized in the inertial range. Circulicity is a quantity associated with angular momentum and, by extension, with vortex motions. The analysis of its inertial properties allows to complete the description of the energy cascade already highlighted in previous works [2,3].
The objective of the post-doctorate is to carry out and exploit direct simulations of compressible homogeneous turbulence with forcing, in order to highlight the inertial properties of circulicity .
To this end, the post-doctoral student will be given access to the very large computing center (TGCC) as well as a code, Triclade, solving the compressible Navier-Stokes equations [4]. This code does not have a forcing mechanism and the first task will therefore be to add this functionality. Once this task has been accomplished, simulations will be carried out by varying the nature of the forcing and in particular the ratio between its solenoidal and dilatational components. These simulations will then be exploited by analyzing the transfer terms of circulicity.
[1] Soulard and Briard. Submitted to Phys. Rev. Fluids. Preprint at arXviv:2207.03761v1
[2] Aluie. Phys. Rev. Lett. 106(17):174502, 2011.
[3] Eyink and Drivas.Phys. Rev. X 8(1):011022, 2018.
[4] Thornber et al. Phys. Fluids 29:105107, 2017.
HPC simulations for PEM fuel cells
The goal is to improve TRUST-FC software -a joint development between LITEN and DES institutes at CEA- for detailed full 3D simulation of hydrogene PEM fuel cells and to run simulations on whole real bipolar plate geometries. Funded by AIDAS virtual lab (CEA/Forshungs Zentrum Juelich), a fully coupled electro-chemical, fluidic and thermal model has been built, based on CEA software TRUST. The model has been benchmarked against its FZJ counterpart (Open fuelcell, based on OpenFoam). The candidate will adapt the software and toolchain to larger and larger meshes up to billion cells meshes required to model a full bipolar plate. Besides, he will introduce two phase flow models in order to address the current technological challenges (local flooding or dryout). This ambitious project is actively supported by close collaboration with CEA/DES and FZJ.
ACCELERATING a DSN SWEEP KERNEL ALGORITHM FOR NEUTRONICS BY PORTING ON GPU.
In the framework of the Programmes Transversaux de Compétences (PTC or literally Cross-XXX Programme), the DES/ISAS/DM2S/SERMA/LLPR and the CEA-DIF are both working on the porting of deterministic neutron transport codes on GPU.
The DM2S within the Energies Direction (DES) is responsible for research and development activities on the numerical methods and codes for reactor physics, amongst which the APOLLO3® code. The neutronics laboratory of CEA-DIF is responsible for developing tools for deterministic methods in neutronics for the Simulation programme.
These two laboratories are actively preparing for the advent of new generation of supercomputers where GPU (Graphical Processing Units) will be predominant. Indeed, the underlying numerical problems to be solved along with the working methodology as well as the conclusions and experience which will be obtained from such studies may be rationalised between both laboratories. Thus, this work has given rise to this postdoctoral position which will be common to both teams. The postdoctoral researcher will be formally based at SERMA at CEA Saclay, with nevertheless regular meetings with the CEA-DIF scientists.
The postdoctoral research work is to study the acceleration of a toy model of a 3D discrete ordinates diamond-differencing sweep kernel (DSN) by porting the code on GPU. This work hinges on porting experiments which have previously been carried by both teams following two different approaches: a ‘’high-level’’ one based on the Kokkos framework for DES and a ‘’low-level’’ approach based on Cuda for CEA-DIF.