Modélisation des écoulements multiphasiques réactifs
Modélisation et simulation de la propagation des feux de forêts
Thermodynamique des mélanges
Thermodynamics, Numerical Waves, Interfaces, Combustion Team
Présentation
The TONIC team is developing an activity of modeling of strongly multi-scale phenomena. It covers in particular multiphase and/or reactive flows, from the scale of the isolated injector (a few mm) to the scale of a fully developed forest fire (several hectares).
Adapted numerical methods are developed in parallel, in particular for soil imaging (detection of slicks by acoustic analysis), or for the modeling of radiative transfers.
In parallel to these multi-scale developments, analytical work is carried out to support the construction of models. An important research effort is devoted to the modeling of the thermodynamics of multiphase mixtures (thermochemical equilibrium calculations, complex thermodynamic closures), or to the development of reduced kinetic models for combustion.
Gabriel Meletti, Stéphane Abide, Stephane Viazzo, Jezabel Curbelo, Uwe Harlander. Wave-like spirals and spontaneous oscillations in strato-rotational flows. Journal of Fluid Mechanics, 2026, 1035, pp.A38. ⟨10.1017/jfm.2026.11559⟩. ⟨hal-05636292⟩ Plus de détails...
This study investigates the dynamics of strato-rotational instability (SRI) in a stratified, rotating fluid, focusing on the interaction between axial modes and spiral components. Through numerical analysis, we find that SRI induces oscillatory behaviours that change the mean flow, leading to the selective activation of distinct axial wavenumbers associated with upward and downward propagating spiral modes. These results suggest wave–mean flow interactions. The use of Radon transforms (RTs) allowed us to separate these spiral components, showing that each upward and downward component was individually modulated, but out of phase with each other. Inspired by the RT findings, a simplified toy model was developed to interpret the spiral pattern changes linked to amplitude modulations. The model considers two wave-like spirals propagating in opposite axial directions, linearly interacting. By incorporating out-of-phase individual spiral modulations, the model reproduces the observed spiral pattern transitions, offering a straightforward interpretation of the underlying physical processes. To explore the mechanism of individual spiral modulations, we consider a quasi-biennial oscillation (QBO)-like framework derived from the Navier–Stokes aligns in a rotating frame. These findings contribute to a better understanding of low-frequency SRI dynamics and may offer insights into similar phenomena in geophysical and astrophysical contexts.
Gabriel Meletti, Stéphane Abide, Stephane Viazzo, Jezabel Curbelo, Uwe Harlander. Wave-like spirals and spontaneous oscillations in strato-rotational flows. Journal of Fluid Mechanics, 2026, 1035, pp.A38. ⟨10.1017/jfm.2026.11559⟩. ⟨hal-05636292⟩
Jinhua Lu, Thomas Gregorczyk, Song Zhao, Pierre Boivin. Phase-field-based recursive regularized multiphase lattice Boltzmann model with a consistent pressure scheme. International Journal of Multiphase Flow, 2026, 195, pp.105500. ⟨10.1016/j.ijmultiphaseflow.2025.105500⟩. ⟨hal-05344425⟩ Plus de détails...
Multiphase lattice Boltzmann models with enhanced stability and no deviation terms.
• Consistent pressure scheme decoupled from density and viscosity variations.
• The proposed model shows superior numerical stablity and accuracy.
Jinhua Lu, Thomas Gregorczyk, Song Zhao, Pierre Boivin. Phase-field-based recursive regularized multiphase lattice Boltzmann model with a consistent pressure scheme. International Journal of Multiphase Flow, 2026, 195, pp.105500. ⟨10.1016/j.ijmultiphaseflow.2025.105500⟩. ⟨hal-05344425⟩
Benoît Péden, Pierre Boivin, Nicolas Odier. Large-Eddy Simulation of a 3D airblast injector using a diffuse interface four-equation model: effects of evaporation and combustion. Combustion and Flame, 2026, 285, pp.114771. ⟨10.1016/j.combustflame.2026.114771⟩. ⟨hal-05557933⟩ Plus de détails...
This work presents Large-Eddy Simulations of a three-dimensional airblast-type injector using a diffuse-interface Multi-Fluid approach. A four-equation model is employed, including a consistent phase transition solver and a thermodynamic closure suitable for evaporating and reacting flows. The influence of evaporation and combustion on the spray and flow dynamics is investigated through a comparative analysis of cold, evaporative, and reactive configurations. The method is first validated against reference results and known behavior for similar injector geometries. It is shown that the addition of evaporation significantly alters the liquid fuel distribution, particularly in the inner recirculation zone, while combustion further modifies both liquid and gaseous fuel fields due to temperatureinduced evaporation and fuel consumption. The reacting case exhibits typical flame features, including hollow cone structures and localized high-temperature zones near stoichiometric mixture fractions. These phenomena align well with expected flame behavior under airblast conditions. Phase transition and combustion also have a notable impact on the velocity field, with increased expansion and stronger recirculation induced by heat release. The proposed model captures these effects in a unified framework. Finally, the present multi-physics approach enables consistent and efficient simulation of multiphase, reactive sprays, providing physical insight into the coupled interaction between atomization, evaporation, and combustion. The method shows good numerical performance on the 3D injector, with a reduced computational time of 2.1 × 10 -5 s.mpi/node/it, which has no overcost compared to the Lagrangian reference model. The fully explicit treatment of the equation of state (NASG) ensures excellent robustness on complex geometries, while avoiding the iterative procedure required by cubic-type EoS. These numerical properties make the DIM suitable for industrial LES configurations involving evaporation and combustion, and further model development.
Benoît Péden, Pierre Boivin, Nicolas Odier. Large-Eddy Simulation of a 3D airblast injector using a diffuse interface four-equation model: effects of evaporation and combustion. Combustion and Flame, 2026, 285, pp.114771. ⟨10.1016/j.combustflame.2026.114771⟩. ⟨hal-05557933⟩
A. Fayet, Stéphane Mimouni, Luc Favre, Catherine Colin, Pierre Boivin, et al.. Modeling and numerical simulation of boiling flows: application and dataset release of the DEBORA experiment. International Journal of Heat and Mass Transfer, In press, ⟨10.2139/ssrn.5718803⟩. ⟨hal-05571906⟩ Plus de détails...
Computational Fluid Dynamics (CFD) is widely used in nuclear engineering for safety related studies or for new design investigations. Co-developed by EDF, CEA, ASNR and Framatome, the NEPTUNE_CFD code is specialized in nuclear thermal-hydraulic applications allowing the simulation of two-phase flows on complex geometries. Recently, a new Heat Flux Partitioning (HFP) model has been proposed by Favre et al. [1] for a thorough description of the boiling phenomena, including, among others, the effect of wall sliding bubbles. However, an excessive increase in computation time follows the subsequent modeling improvement. This paper presents an optimization of the bubble sliding calculation returning to a reasonable computation time compatible with industrial applications. The newly developed model is then validated using NEPTUNE_CFD and compared to the DEBORA experimental data, featuring R12 coolant boiling flow within a characteristic non-dimensional scope of a Pressurized Water Reactor (PWR). The improvement in the wall temperature calculation is demonstrated by several simulations implementing the new HFP model. To support the community's validation and benchmarking efforts, the complete DEBORA experimental dataset is made publicly available for the first time as part of this work, provided under a CC BY 4.0 license. This contribution advances both modeling capabilities and data availability, promoting transparency and reproducibility in multiphase CFD for nuclear applications
A. Fayet, Stéphane Mimouni, Luc Favre, Catherine Colin, Pierre Boivin, et al.. Modeling and numerical simulation of boiling flows: application and dataset release of the DEBORA experiment. International Journal of Heat and Mass Transfer, In press, ⟨10.2139/ssrn.5718803⟩. ⟨hal-05571906⟩
Journal: International Journal of Heat and Mass Transfer
Jian Cardenas, Song Zhao, Isabelle Raspo, Guillaume Chiavassa, Pierre Boivin. Lattice-Boltzmann methods for supercritical fluids flows. Journal of Supercritical Fluids, 2026, 230, pp.106838. ⟨10.1016/j.supflu.2025.106838⟩. ⟨hal-05608205⟩ Plus de détails...
This article develops two algorithms for the thermal Lattice Boltzmann Method (LBM) to simulate supercritical fluid dynamics using real-fluid thermodynamics. The first algorithm employs a compressible formulation (LBM-C), while the second is based on a Low Mach Number approximation (LBM-LMN) to enhance computational efficiency without compromising physical fidelity. Both approaches use a conservative scheme compatible with cubic equations of state (EOS), enabling accurate representation of non-ideal behaviors under supercritical conditions and ensuring mass, momentum, and energy conservation. Validation on canonical supercritical flow benchmarks demonstrates that the LBM-LMN approach achieves accuracy comparable to the compressible formulation while reducing the overall computational time by a factor of about 15, as quantified by the RTTS parameter.
Jian Cardenas, Song Zhao, Isabelle Raspo, Guillaume Chiavassa, Pierre Boivin. Lattice-Boltzmann methods for supercritical fluids flows. Journal of Supercritical Fluids, 2026, 230, pp.106838. ⟨10.1016/j.supflu.2025.106838⟩. ⟨hal-05608205⟩
8 juin 2026
- Modeling of supercritical CO2 flow and heat transfer applied to an innovative thermal machine structure / Soutenance de thèse Jian Cardenas Cabezas
Doctorant : Jian CARDENAS CABEZAS
Date et lieu : e lundi 8 juin 2026 à 12h30 dans l’Amphithéâtre N°3, Centrale Méditerranée ; 38 rue Frédéric Joliot Curie,13013 Marseille
Abstract: Improving the energy efficiency of industrial systems is a major challenge in the energy transition. A significant share of primary energy is dissipated as waste heat, particularly at low and medium temperatures, where recovery solutions remain limited. In this context, supercritical carbon dioxide (sCO2 ) cycles are attracting growing interest because of their compactness, their high efficiency potential, and the favorable properties of the fluid near the critical point. However, their deployment is still hindered by technological barriers and by the difficulty of accurately modeling the real fluid. This thesis investigates an innovative concept of thermally driven compression, called SHREC, developed with the industrial partner CIXTEN. Unlike conventional Brayton cycles, which rely on an electrically driven mechanical compressor, the SHREC device directly converts thermal energy, ideally recovered from waste heat, into compression work. By reducing the need for mechanical compression, it aims to lower auxiliary electricity consumption and improve the overall efficiency of sCO2 cycles. The main scientific difficulty lies in modeling supercritical CO2. Near the critical point, small variations in temperature or density induce strong nonlinear changes in density, heat capacity, compressibility, and transport properties. Under such conditions, ideal-gas assumptions or constant-property approximations become unsuitable. In addition, the SHREC device operates at low Mach number and in complex geometries, which imposes stringent requirements in terms of stability, mass conservation, and thermodynamic consistency. To address these challenges, a multi-scale modeling approach was developed. At the system scale, a zero-dimensional (0D) thermodynamic model, based on a cubic equation of state, was established to describe the behavior of the supercritical fluid. It allows rapid prediction of overall performance and enables cycle analysis. Experimental tests on a reduced-scale prototype were carried out to measure the evolution of pressure, temperature, produced power, coefficient of performance, and exergy destruction. The results show that the main irreversibilities are located in the expander-compressor unit and in the main heat exchangers. At the local scale, a three-dimensional numerical framework based on the Lattice Boltzmann Method was developed and adapted to real-fluid conditions. A compressible low-Mach-number formulation was implemented and coupled with a cubic equation of state to reproduce supercritical thermodynamic effects. Complex geometries were handled using the Immersed Boundary Method. The approach was first validated on supercritical jets from the literature, and then applied to the SHREC device to resolve transient temperature and pressure fields. The comparison between CFD simulations and the 0D model shows that the assumption of instantaneous pressure equilibrium between the chambers is not strictly satisfied. The 0D model does not capture local pressure imbalances, but it reproduces the overall compression and expansion dynamics correctly, whereas the CFD simulations provide a detailed description of the internal flow and heat transfer mechanisms. This work shows that reliable investigation of sCO2 systems requires a consistent combination of real-fluid thermodynamics, reduced-order modeling, numerical simulations, and experimental validation. It contributes to the development of thermally driven compression technologies for waste heat recovery and lays the foundations for compact systems for next-generation energy cycles.
M. RIBERT Guillaume INSA Rouen Normandie - Examinateur
Mme RASPO Isabelle CNRS, M2P2 - Examinatrice
M. MELDI Marcello ENSAM - Président du Jury
M. FERRASSE Jean-Henry Aix Marseille Université, M2P2 - Co-encadrant
M. FAVIER Julien Aix Marseille Université, M2P2 - Co-Directeur de thèse
M. SCHMITT Thomas CNRS, EM2C -Rapporteur
M. SILVA Gonçalo Universidade de Évora, Portugal - Rapporteur
6 février 2026
- Study of turbulent transport of energetic particles in nuclear fusion plasmas nuclear fusion plasmas by trajectory simulations and artificial intelligence techniques / PhD Defense Benoît Clavier
Doctorant : Benoît CLAVIER
Date et lieu : le vendredi 6 février à 14h00, M2P2 - salle Labus, Centrale Méditerranée
Abstract: This thesis studies the turbulent transport of charged particles in magnetized fusion plasmas by combining reduced turbulence models, numerical trajectory simulations, and data-driven approaches based on artificial intelligence. After presenting the physical framework of radial transport in a tokamak and the Hasegawa–Wakatani model, Eulerian and Lagrangian diagnostics are developed to obtain reference transport measurements. The work then analyzes the transport of test particles in different turbulent regimes, highlighting the limitations of certain classical approximations and the complexity of energetic particle dynamics. The study is extended to a more realistic three-dimensional ion-temperature-gradient (ITG) turbulence, allowing scaling laws for radial diffusion to be established. Finally, a synthetic turbulence generation model based on a Convolutional Variational Autoencoder (CVAE) coupled with a dynamic model is proposed to efficiently reproduce turbulence and accelerate transport studies, illustrating the potential of data-driven approaches for future research in plasma physics.
Jury
David ZARZOSO-FERNANDEZ - Chargé de recherche - CNRS M2P2 - Directeur de thèse
Emmanuel FRéNOD - Professeur des universités - Université Bretagne Sud - Co-directeur de thèse
Victor TRIBALDOS - Professeur des universités - Universidad Carlos III de Madrid - Rapporteur
Julien LE SOMMER - Directeur de recherche - CNRS, IGE Grenoble - Examinateur
Maxime LESUR - Professeur des universités - Université de Lorraine - Institut Jean Lamour - Rapporteur
Mitra FOULADIRAD - Professeure des universités - Centrale Méditerranée - Président
