Prediction of elastic wave dispersion effects using a semi-analytical model under high-frequency approximation
Ultrasonic testing (UT) methods are a fundamental component of non-destructive testing (NDT). They are widely used to inspect mechanical components such as welds (in nuclear and petrochemical industries) and composite material structures (in aeronautics). To understand the physical phenomena involved in a given configuration, simulation is a valuable tool and sometimes an essential step in implementing the inspection process.
Modeling approaches fall into two main categories: purely numerical models based on finite elements (FE) and semi-analytical methods derived from high-frequency (HF) approximations, such as paraxial rays. While the latter are often favored for their computational efficiency, they introduce simplifications that can compromise the quantitative accuracy of results, particularly for phenomena like dispersion (variation in wave speed with frequency), which are common in certain industrial contexts.
This thesis project aims to enhance the paraxial ray approach by integrating models of dispersive interfaces (composite interplies, coupling layers), dispersive viscoelastic media, and a modal guided wave model. The goal is to develop a simulation tool capable of faithfully reproducing realistic inspection configurations, thereby improving the representativeness of the results.
Electromagnetic Signature Modeling and AI for Radar Object Recognition
This PhD thesis offers a unique opportunity to work at the crossroads of electromagnetics, numerical simulations, and artificial intelligence, contributing to the development of next-generation intelligent sensing and recognition systems. The intern will join the Antenna & Propagation Laboratory at CEA-LETI, Grenoble (France), a world-class research environment equipped with state-of-the-art tools for propagation channel characterization and modelling. A collaboration with the University of Bologna (Italy) is planned during the PhD.
This PhD thesis aims to develop advanced electromagnetic models of near-field radar backscattering, tailored to radar and Joint Communication and Sensing (JCAS) systems operating at mmWave and THz frequencies. The research will focus on the physics-based modeling of the radar signatures of extended objects, accounting for near-field effects, multistatic and multi-antenna configurations, as well as the influence of target materials and orientations. These models will be validated through electromagnetic simulations and dedicated measurement campaigns, and subsequently integrated into scene-level and multipath propagation simulation tools based on ray tracing. The resulting radar signatures will be exploited to train artificial intelligence algorithms for object recognition, material property inference, and radar imaging. In parallel, physics-assisted AI approaches will be investigated to accelerate electromagnetic simulations and reduce their computational complexity. The final objective of the thesis is to integrate radar backscattering-based information into a 3D Semantic Radio SLAM framework, in order to improve localization, mapping, and environmental understanding in complex or partially obstructed scenarios.
We are seeking a student at engineering school or Master’s level (MSc/M2), with a strong background in signal processing, electromagnetics, radar, or telecommunications. An interest in artificial intelligence, physics-based modeling, and numerical simulation is expected. Programming skills in Matlab and/or Python are appreciated, as well as the ability to work at the interface between theoretical models, simulations, and experimental validation. Scientific curiosity, autonomy, and strong motivation for research are essential.The application must include a CV, academic transcripts, and a motivation letter.
Hybrid CPU-GPU Preconditioning Strategies for Exascale Finite Element Simulations
Exascale supercomputers are based on heterogeneous architectures that combine CPUs and GPUs, making it necessary to redesign numerical algorithms to fully exploit all available resources. In large-scale finite element simulations, the solution of linear systems using iterative solvers and algebraic multigrid (AMG) preconditioners remains a major performance bottleneck.
The objective of this PhD is to study and develop hybrid preconditioning strategies adapted to such heterogeneous systems. The work will investigate how multilevel and AMG techniques can be structured to efficiently use both CPUs and GPUs, without restricting computations to a single type of processor. Particular attention will be paid to data distribution, task placement, and CPU–GPU interactions within multilevel solvers.
From a numerical point of view, the research will focus on the analysis and construction of multilevel operators, including grid hierarchies, intergrid transfer operators, and smoothing procedures on avalible GPU's and CPU's. The impact of these choices on convergence, spectral properties, and robustness of preconditioned iterative methods will be studied. Mathematical criteria guiding the design of efficient hybrid preconditioners will be investigated and validated on representative finite element problems, e.g., regional-scale earthquake analysis.
These developments will be coupled with domain decomposition and parallelization strategies adapted to heterogeneous architectures. Particular attention will be paid to CPU–GPU data transfers, memory usage, and the balance between compute-bound and memory-bound kernels. The interaction between numerical choices and hardware constraints, such as CPU and GPU memory hierarchies, will be designed and developed to ensure scalable and efficient implementations.
Control & optimization of fuel cell temperature
Proton exchange membrane fuel cells (PEMFC) represent a key technology for the development of clean and sustainable energy systems, particularly for heavy-duty transport applications where their energy density is very attractive. However, in order to represent a viable industrial alternative, a number of obstacles still need to be overcome, including operating costs and, above all, the durability of the systems under real-world conditions. Among the levers for action, optimizing operating conditions is a promising avenue for limiting the degradation phenomena occurring within the cell. The operating temperature is a particularly key parameter because it affects all aspects of the system, from the kinetics of degradation mechanisms to the thermal capacity that the system can dissipate, including the water balance within the fuel cell. Despite the influence of this parameter on durability, it is generally only optimized at the system level to achieve the best performance, the shortest possible response time and to limit the size of the thermal management system.
The aim of this thesis is to work on optimizing the temperature management of a fuel cell within a system, taking into account not only performance but also sustainability criteria. To do this, the impact of operating temperature on degradation mechanisms will be analyzed using various simulation tools already available at LITEN and the teams' fifteen years of experience in studying PEMFC fuel cell degradation. Various thermal architectures will be proposed and evaluated in conjunction with the work on temperature control optimization. The latter will be implemented on a real fuel cell system in order to demonstrate the relevance of the proposed solution using concrete experimental data.
Development of a new numerical scheme, based on T-coercivity, for discretizing the Navier-Stokes equations.
In the TrioCFD code, the discretization of the Navier-Stokes equations leads to a three-step algorithm (see Chorin'67, Temam'68): velocity prediction, pressure solution, velocity correction. If an implicit time discretization scheme is to be used, the pressure solution step is particularly costly. Thus, most simulations are performed using an explicit time scheme, for which the time step depends on the mesh size, which can be very restrictive. We would like to develop an implicit time discretization scheme using a stabilized formulation of the Navier-Stokes problem based on explicit T-coercivity (see Ciarlet-Jamelot'25). It would then be possible to solve an implicit scheme directly without a correction step, which could significantly improve the performance of the calculations. This would also allow the use of the P1-P0 finite element pair, which is frugal in terms of degrees of freedom but unstable for a classical formulation.
Development of automatic gamma spectrum analysis using a hybrid machine learning algorithm for the radiological characterization of nuclear facilities decommissioning.
The application of gamma spectrometry to radiological characterization in nuclear facility decommissioning, requires the development of specific algorithms for automatic gamma spectrum analysis. In particular, the classification of concrete waste according to its level of contamination, is a crucial issue for controlling decommissioning costs.
Within CEA/List, LNHB, in collaboration with CEA/DEDIP, has been involved for several years in the development of tools for the automatic analysis of low-statistics gamma spectra, which can be applied to scintillator detectors (NaI(Tl), plastics). In this context, an original approach based on a hybrid machine learning/statistics spectral unmixing algorithm has been developed for the identification and quantification of radionuclides in the presence of significant deformations in the measured spectrum, due in particular to interactions between the gamma emission from the radioactive source and its environment.
The proposed subject follows on from thesis work that led to the development of the hybrid algorithm with the aim of extending this approach to the radiological characterization of concrete surfaces. The candidate will be involved in the evolution of the hybrid machine learning/statistical algorithm for the characterization of concrete for classification as conventional waste. The work will include a feasibility study of modeling the deviations of the learned model to optimize the robustness of decision-making.
Novel architecture and signal processing for mobile optical telecommunications
Free-Space Optical Communications (FSO) rely on transmitting data via light between two distant points, eliminating the need for fibers or cables. This approach is particularly valuable when wired connections are impractical or prohibitively expensive.
However, these links are highly susceptible to atmospheric conditions—fog, rain, dust, and thermal turbulence—which attenuate or distort the light beam, significantly degrading communication quality. Current solutions remain costly and limited, both in terms of optical compensation hardware and signal processing algorithms.
Within this framework, the thesis aims to design high-performance, robust mobile optical links capable of adapting to dynamic and disturbed environments. The study will focus on leveraging Silicon-based Optical Phased Arrays (OPAs)—a technology derived from low-cost LiDAR systems—offering a promising path toward compact, integrated, and cost-effective architectures.
The primary focus of the research will be developing advanced algorithmic approaches for signal processing and compensation. The PhD candidate will be tasked with designing a dedicated simulation environment to evaluate and validate architectural choices and algorithmic strategies before practical experimentation.
The overarching goal is to propose an integrated, flexible, and reliable architecture that ensures uninterrupted optical communication in motion, with potential applications in aerospace, space, and terrestrial domains.
Topologic optimization of µLED's optical performance
The performance of micro-LEDs (µLEDs) is crucial for micro-displays, a field of expertise at the LITE laboratory within CEA-LETI. However, simulating these components is complex and computationally expensive due to the incoherent nature of light sources and the involved geometries. This limits the ability to effectively explore multi-parameter design spaces.
This thesis proposes to develop an innovative finite element method to accelerate simulations and enable the use of topological optimization. The goal is to produce non-intuitive designs that maximize performance while respecting industrial constraints.
The work is divided into three phases:
- Develop a fast and reliable simulation method by incorporating appropriate physical approximations for incoherent sources and significantly reducing computation times.
- Design a robust topological optimization framework that includes fabrication constraints to generate immediately realizable designs.
- Realize such a metasurface on an existing shortloop in the laboratory. This part is optional and will be tackled only if we manage to seize an Opportunity to finance the prototype, via the inclusion of the thésis inside the "metasurface
topics" of european or IPCEI projets in the lab .
The expected results include optimized designs for micro-displays with enhanced performance and a methodology that can be applied to other photonic devices and used by other laboratories from DOPT.
Investigation of polytopal methods apllied to CFD and optimized on GPU architecture
This research proposal focuses on the study and implementation of polytopal methods for solving the equations of fluid mechanics. These methods aim to handle the most general meshes possible, overcoming geometric constraints or those inherited from CAD operations such as extrusions or assemblies that introduce non-conformities. This work also falls within the scope of high-performance computing, addressing the increase in computational resources and, in particular, the development of massively parallel computing on GPUs.
The objective of this thesis is to build upon existing polytopal methods already implemented in the TRUST software, specifically the Compatible Discrete Operator (CDO) and Discontinuous Galerkin (DG) methods. The study will be extended to include convection operators and will investigate other methods from the literature, such as Hybrid High Order (HHO), Hybridizable Discontinuous Galerkin (HDG), and Virtual Element Method (VEM).
The main goals are to evaluate:
1. The numerical behavior of these different methods on the Stokes/Navier-Stokes equations;
2. The adaptability of these methods to heterogeneous architectures such as GPUs.
Robust multi-material topological optimization under manufacturability constraints applied to the design of superconducting magnets for high-field MRI
MRI scanners are invaluable tools for medicine and research, whose operation is based on exploiting the properties of atomic nuclei immersed in a very intense static magnetic field. In almost all MRI scanners, this field is generated by a superconducting electromagnet.
The design of electromagnets for MRI must meet very demanding requirements in terms of the homogeneity of the field produced. In addition, as the magnetic field becomes more intense, the forces exerted on the electromagnet increase, raising the issue of the mechanical strength of the windings. Finally, the “manufacturability” of the electromagnet imposes constraints on the shapes of acceptable solutions. The design of superconducting electromagnets for MRI therefore requires a meticulous effort to optimize the design, subject to constraints based on magneto-mechanical multiphysics modeling.
A new innovative multiphysics topological optimization methodology has been developed, based on a density method (SIMP) and a finite element code. This has made it possible to produce magnet designs that meet the constraints on the homogeneity of the magnetic field produced and on the mechanical strength of the windings. However, the solutions obtained are not feasible in practice, both in terms of the manufacturability of the coils (cable windings) and their integration with a supporting structure (coils held in place by a steel structure).
The objective of this thesis is to enhance the topological optimization method by formalizing and implementing manufacturing constraints related to the winding method, residual stresses resulting from pre-tensioning the cables during winding, and the presence of a structural material capable of absorbing the forces transmitted by the coils.