N. Fedorczak, C. Arnas, L. Cappelli, L. Colas, Y. Corre, et al.. Survey of tungsten gross erosion from main plasma facing components in WEST during a L-mode high fluence campaign. Nuclear Materials and Energy, 2024, 41 (4), pp.101758. ⟨10.1016/j.nme.2024.101758⟩. ⟨cea-04816563⟩ Plus de détails...
An initial high fluence campaign was performed in WEST, in 2023, on the newly installed actively cooled tungsten divertor composed of ITER-grade monoblocks. The campaign consisted in the repetition of a 60 s long Deuterium L-mode pulse in attached divertor conditions, cumulating over 10000s of plasma exposure. A maximum deuterium fluence of approximately 5 ⋅ 10 26 m -2 was reached in the outer strike point region, representative of a few high performance ITER pulses. Gross tungsten erosion inferred from visible spectroscopy shows that the most eroded plasma facing component is the inner divertor target with rates ten times larger than on the outer divertor target. The outer midplane tungsten bumpers, located a few centimeters from the plasma, show gross erosion rates two times lower than at the outer divertor. We conclude that the outer midplane bumpers have a negligible contribution to the long range tungsten migration and deposition onto the lower divertor. The cumulated gross erosion rate on the inner divertor translates in an effective gross erosion thickness of about 20 μm, while it is about 2 μm for the outer divertor. Strikingly, these orderings coincide with the thickness of deposits found locally on the divertor: the exposed surfaces of high field side monoblocks are covered with several tens of μm tungsten deposits, while on the lower field side, few μm thin tungsten deposits are only found on the magnetically shadowed parts of monoblocks. The strong impact of those deposits on WEST operation, namely perturbation of surface temperature measurement with infra-red thermography, and the emission of flakes causing radiative perturbation of the confined plasma, calls for anticipating similar issues in ITER. In particular, the start of research operation shall consider the definition of a divertor erosion budget in order to anticipate the formation of deleterious deposits.
N. Fedorczak, C. Arnas, L. Cappelli, L. Colas, Y. Corre, et al.. Survey of tungsten gross erosion from main plasma facing components in WEST during a L-mode high fluence campaign. Nuclear Materials and Energy, 2024, 41 (4), pp.101758. ⟨10.1016/j.nme.2024.101758⟩. ⟨cea-04816563⟩
Stefano Di Genova, Alberto Gallo, Luca Cappelli, Nicolas Fedorczak, Hugo Bufferand, et al.. Global analysis of tungsten migration in WEST discharges using numerical modelling. Nuclear Fusion, 2024, ⟨10.1088/1741-4326/ad82f9⟩. ⟨hal-04739577⟩ Plus de détails...
Plasma discharges in the tungsten (W) Environment Steady-state Tokamak (WEST) are strongly impacted by W contamination. In WEST experiments, due to W contamination, the power radiated in the plasma (PRad) is on average, around 50% of the total power injected into the plasma (PTOT). Furthermore, this radiated power fraction (fRad) is almost insensitive to plasma conditions. The causes behind this experimental trend are not fully understood. In this contribution, a 3D numerical model is used to analyze the W migration in the WEST boundary plasma in different plasma scenarios. The WEST experimental database is sampled to obtain a scan of simulation input parameters. These parameters mimic the WEST plasma conditions over a chosen experimental campaign. The simulation results are compared to WEST diagnostics measurements (reflectometry, Langmuir probes, and visible spectroscopy) to verify that the simulated plasma conditions are representative of the WEST database. The W contamination trend is analysed: the W density (nW) strongly decreases when the radial distance between the separatrix and WEST antennas (Radial Outer Gap, ROG) increases. On the other hand, at a given ROG, nW increases proportionally with the power entering the scrape-off layer (PSOL). PRad is estimated with a simple 0D model. For a fixed ROG, fRad is not sensitive to plasma conditions. These trends are qualitatively and, at times, quantitatively comparable to what is observed in WEST experiments: the simulated trends are related to the poorly screened W influx caused by the erosion of the main chamber Plasma-Facing Components (PFCs). Thus, this numerical analysis suggests a possible interpretation of WEST experimental trends.
Stefano Di Genova, Alberto Gallo, Luca Cappelli, Nicolas Fedorczak, Hugo Bufferand, et al.. Global analysis of tungsten migration in WEST discharges using numerical modelling. Nuclear Fusion, 2024, ⟨10.1088/1741-4326/ad82f9⟩. ⟨hal-04739577⟩
M. Scotto D’abusco, I. Kudashev, G. Giorgiani, Anna Glasser, F. Schwander, et al.. First integrated core-edge fluid simulation of ITER’s Limiter-Divertor transition with SolEdge-HDG. Nuclear Materials and Energy, 2024, pp.101750. ⟨10.1016/j.nme.2024.101750⟩. ⟨hal-04720290⟩ Plus de détails...
This work explores the Limiter-Divertor transition (L-D) during the current ramp-up of ITER's Q=10 baseline plasma scenario at various central line-integrated density n_(li) values. The analysis, based on transport simulations performed with the latest version of SoleEdge-HDG, focuses on the time evolution of heat and ion particle fluxes, revealing regions of elevated temperature on the inner wall and plasma-facing components (PFCs) despite moderate loads.The investigation also delves into the effects of perpendicular convection flux terms on density build-up, comparing different formulations and their interplay with auxiliary heating sources. Furthermore, the paper shows the impact of taking into account the evolution of the parallel neutral momentum on plasma and neutral density at the targets in the context of an ITER steady-state scenario.
M. Scotto D’abusco, I. Kudashev, G. Giorgiani, Anna Glasser, F. Schwander, et al.. First integrated core-edge fluid simulation of ITER’s Limiter-Divertor transition with SolEdge-HDG. Nuclear Materials and Energy, 2024, pp.101750. ⟨10.1016/j.nme.2024.101750⟩. ⟨hal-04720290⟩
Q. Mao, Umberto d'Ortona, J. Favier. Hydrodynamic coupling of a cilia–mucus system in Herschel–Bulkley flows. Journal of Fluid Mechanics, 2024, 994, pp.A8. ⟨10.1017/jfm.2024.600⟩. ⟨hal-04735292⟩ Plus de détails...
The yield stress and shear thinning properties of mucus are identified as critical for ciliary coordination and mucus transport in human airways. We use here numerical simulations to explore the hydrodynamic coupling of cilia and mucus with these two properties using the Herschel–Bulkley model, in a lattice Boltzmann solver for the fluid flow. Three mucus flow regimes, i.e. a poorly organized regime, a swirly regime, and a fully unidirectional regime, are observed and analysed by parametric studies. We systematically investigate the effects of ciliary density, interaction length, Bingham number and flow index on the mucus flow regime formation. The underlying mechanism of the regime formation is analysed in detail by examining the variation of two physical quantities (polarization and integral length) and the evolution of the flow velocity, viscosity and shear-rate fields. Mucus viscosity is found to be the dominant parameter influencing the regime formation when enhancing the yield stress and shear thinning properties. The present model is able to reproduce the solid body rotation observed in experiments (Loiseau et al. , Nat. Phys. , vol. 16, 2020, pp. 1158–1164). A more precise prediction can be achieved by incorporating non-Newtonian properties into the modelling of mucus as proposed by Gsell et al. ( Sci. Rep. , vol. 10, 2020, 8405).
Q. Mao, Umberto d'Ortona, J. Favier. Hydrodynamic coupling of a cilia–mucus system in Herschel–Bulkley flows. Journal of Fluid Mechanics, 2024, 994, pp.A8. ⟨10.1017/jfm.2024.600⟩. ⟨hal-04735292⟩
L Cappelli, N Fedorczak, E Serre. Semi-analytical modelling of prompt redeposition in a steady-state plasma. Nuclear Fusion, 2024, 64 (10), pp.106028. ⟨10.1088/1741-4326/ad6c5e⟩. ⟨hal-04685227⟩ Plus de détails...
A steady-state, 1D semi-analytical model for prompt redeposition based on the separation between redeposition caused by the electric field in the sheath and redeposition related to gyromotion is here described. The model allows for the estimation of not only the fraction of promptly redeposited flux but also the energy and angular distribution of the non-promptly redeposited population, along with their average charge state. Thus, the temperature and mean parallel-to-B velocity of the non-promptly redeposited flux are also available. The semi-analytical model was validated against equivalent Monte Carlo simulations across a broad range of input parameters. In this paper the eroded material under exam was tungsten (W) for which the code demonstrated consistent agreement with respect to numerical results, within its defined validity limits. The model can theoretically provide a solution for any material, temperature and electron density profile in the sheath, monotonic potential drop profile, and sputtered particles energy and angular distribution at the wall. As such, this code emerges as a potential tool for addressing the boundary redeposition phenomenon in fluid impurity transport simulations.
L Cappelli, N Fedorczak, E Serre. Semi-analytical modelling of prompt redeposition in a steady-state plasma. Nuclear Fusion, 2024, 64 (10), pp.106028. ⟨10.1088/1741-4326/ad6c5e⟩. ⟨hal-04685227⟩
R Varennes, G Dif-Pradalier, P Ghendrih, V Grandgirard, O Panico, et al.. Turbulent relaxation patterns in SOL plasma. Plasma Physics and Controlled Fusion, 2024, 66 (10), pp.105008. ⟨10.1088/1361-6587/ad705c⟩. ⟨hal-04797031⟩ Plus de détails...
Abstract Relaxations of localized over-density in a plane transverse to the magnetic field are numerically investigated under the effect of drift-wave and interchange drives in SOL conditions. Such a controlled departure from thermodynamic equilibrium allows the investigation of fundamental processes at play in cross-field transport. Interchange instabilities generate ballistic outward radial flux with low amplitude zonal flow patterns, whereas drift-wave instabilities result in symmetric radial flux with large amplitude zonal flow patterns. When both instabilities are considered, the combined effects tend to favor drift-waves, leading to a weaker outward flux with larger zonal flow patterns.
R Varennes, G Dif-Pradalier, P Ghendrih, V Grandgirard, O Panico, et al.. Turbulent relaxation patterns in SOL plasma. Plasma Physics and Controlled Fusion, 2024, 66 (10), pp.105008. ⟨10.1088/1361-6587/ad705c⟩. ⟨hal-04797031⟩
I Kudashev, M Scotto D’abusco, A Glasser, E Serre, F Schwander, et al.. Global particle buildup simulations with gas puff scan: application to WEST discharge. Frontiers in Physics, 2024, 12, ⟨10.3389/fphy.2024.1407534⟩. ⟨hal-04703430⟩ Plus de détails...
This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheathlimited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.
I Kudashev, M Scotto D’abusco, A Glasser, E Serre, F Schwander, et al.. Global particle buildup simulations with gas puff scan: application to WEST discharge. Frontiers in Physics, 2024, 12, ⟨10.3389/fphy.2024.1407534⟩. ⟨hal-04703430⟩
Elena Alekseenko, A.A. Sukhinov, B. Roux. Modeling of multi-fractional suspended particle pathways in a shallow water basin under influence of strong winds. Regional Studies in Marine Science, 2024, 73, pp.103477. ⟨10.1016/j.rsma.2024.103477⟩. ⟨hal-04515082⟩ Plus de détails...
In this study, we investigate the complex dynamics of multi-fractional suspended particle transport in a shallow water basin subjected to strong wind conditions. Our research focuses on understanding the interplay between wind-induced advection and particle settlement, and its implications for sediment redistribution. Through our analysis, we reveal the distinct behaviors of different sediment fractions. Clay particles, constituting the lowest fraction in sediment cores, remain suspended throughout the simulation due to their low settlement velocity, with relatively stable concentrations. Conversely, the dominant fraction, medium silt, is suspended during intense wind events but quickly settles to the bed due to its higher settling velocity. Wind stress exceeding 0.05 Pa triggers particulate matter erosion, leading to its presence in the water column. Additionally, we explore the 2D distribution of sediment characteristics, including thickness, dry density, and mud fraction, to identify areas prone to erosion and deposition. Our findings demonstrate that coastal areas of the Taganrog Bay experienced significant erosion following strong wind events, exhibiting the thinnest sediment thickness and the highest dry bulk density. Deposition areas, characterized by thicker sediment layers and lower dry density, were often found in proximity to erosion zones, indicating the influence of particle resuspension and settlement processes. Furthermore, we analyze the implications of our findings on the vulnerability of specific regions to erosion and deposition. The central part of the sea contains moderately thicker sediment layers with a moderately high mud fraction, representing a zone of fine sediment accumulation. These fine sediments, including fine silt and clay, remain suspended for longer durations and are redistributed over greater distances by currents. Overall, our study provides valuable understanding into the multi-fractional suspended particle pathways and their interaction with strong winds in shallow water basins. The results contribute to a better understanding of sediment dynamics, which has implications for coastal management, environmental monitoring, and the preservation of benthic ecosystems.
Elena Alekseenko, A.A. Sukhinov, B. Roux. Modeling of multi-fractional suspended particle pathways in a shallow water basin under influence of strong winds. Regional Studies in Marine Science, 2024, 73, pp.103477. ⟨10.1016/j.rsma.2024.103477⟩. ⟨hal-04515082⟩
Alix Limoges, Jacques Piazzola, Christophe Yohia, Quentin Rodier, William Bruch, et al.. Study of the Atmospheric Transport of Sea-Spray Aerosols in a Coastal Zone Using a High-Resolution Model. Atmosphere, 2024, 15 (6), pp.702. ⟨10.3390/atmos15060702⟩. ⟨hal-04779475⟩ Plus de détails...
Fine-scale models for the transport of marine aerosols are of great interest for the study of micro-climates and air quality in areas of complex topography, such as in urbanized coastal areas. To this end, the MIO laboratory implemented the Meso-NH model in its LES version over the northwest Mediterranean coastal zone using a recent sea-spray source function. Simulated meteorological parameters and aerosol concentrations are compared to experimental data acquired in the Mediterranean coastal zone in spring 2008 on board the R/V Atalante. Key findings indicate that the large eddy simulation (LES) mode closely matches with the experimental data, enabling an in-depth analysis of the numerical model ability to predict variations in aerosol concentrations. These variations are influenced by different wind directions, which lead to various fetch distances typical of coastal zones.
Alix Limoges, Jacques Piazzola, Christophe Yohia, Quentin Rodier, William Bruch, et al.. Study of the Atmospheric Transport of Sea-Spray Aerosols in a Coastal Zone Using a High-Resolution Model. Atmosphere, 2024, 15 (6), pp.702. ⟨10.3390/atmos15060702⟩. ⟨hal-04779475⟩
Louis Lamérand, Didier Auroux, Philippe Ghendrih, Francesca Rapetti, Eric Serre. Inverse problem for determining free parameters of a reduced turbulent transport model for tokamak plasma. Advances in Computational Mathematics, 2024, 50 (3), pp.39. ⟨10.1007/s10444-024-10135-6⟩. ⟨hal-04569449⟩ Plus de détails...
Two-dimensional transport codes for the simulation of tokamak plasma are reduced version of full 3D fluid models where plasma turbulence has been smoothed out by averaging. One of the main issues nowadays in such reduced models is the accurate modelling of transverse transport fluxes resulting from the averaging of stresses due to fluctuations. Transverse fluxes are assumed driven by local gradients, and characterised by ad hoc diffusion coefficients (turbulent eddy viscosity), adjusted by hand in order to match numerical solutions with experimental measurements. However, these coefficients vary substantially depending on the machine used, type of experiment and even the location inside the device, reducing drastically the predictive capabilities of these codes for a new configuration. To mitigate this issue, we recently proposed an innovative path for fusion plasma simulations by adding two supplementary transport equations to the mean-flow system for turbulence characteristic variables (here the turbulent kinetic energy k and its dissipation rate ) to estimate the turbulent eddy viscosity. The remaining free parameters are more driven by the underlying transport physics and hence vary much less between machines and between locations in the plasma. In this paper, as a proof of concept, we explore, on the basis of digital twin experiments, the efficiency of the assimilation of data to fix these free parameters involved in the transverse turbulent transport models in the set of averaged equations in 2D.
Louis Lamérand, Didier Auroux, Philippe Ghendrih, Francesca Rapetti, Eric Serre. Inverse problem for determining free parameters of a reduced turbulent transport model for tokamak plasma. Advances in Computational Mathematics, 2024, 50 (3), pp.39. ⟨10.1007/s10444-024-10135-6⟩. ⟨hal-04569449⟩
Franck Corset, Mitra Fouladirad, Christian Paroissin. Imperfect and worse than old maintenances for a gamma degradation process. Applied Stochastic Models in Business and Industry, 2024, ENBIS 2022, 40 (3), pp.620-639. ⟨10.1002/asmb.2849⟩. ⟨hal-04462980⟩ Plus de détails...
This article considers a condition‐based maintenance for a system subject to deterioration. The deterioration is modeled by a non‐homogeneous gamma process, more precisely the gamma process and the preventive maintenance are imperfect or worse than old. The corrective maintenance actions are as good as new. The maintenance efficiency or non‐efficiency parameters as well as the deterioration parameters are considered to be unknown. The monitoring data under consideration give indirect information on the maintenance parameters. Therefore, an expected maximum algorithm is applied for parameter estimation.
Franck Corset, Mitra Fouladirad, Christian Paroissin. Imperfect and worse than old maintenances for a gamma degradation process. Applied Stochastic Models in Business and Industry, 2024, ENBIS 2022, 40 (3), pp.620-639. ⟨10.1002/asmb.2849⟩. ⟨hal-04462980⟩
Journal: Applied Stochastic Models in Business and Industry
Uwe Ehrenstein. Generalization to differential–algebraic equations of Lyapunov–Schmidt type reduction at Hopf bifurcations. Communications in Nonlinear Science and Numerical Simulation, 2024, 131, pp.107833. ⟨10.1016/j.cnsns.2024.107833⟩. ⟨hal-04408097⟩ Plus de détails...
The Lyapunov-Schmidt procedure, a well-known and powerful tool for the local reduction of nonlinear systems at bifurcation points or for ordinary differential equations (ODEs) at Hopf bifurcations, is extended to the context of strangeness-free differential-algebraic equations (DAEs), by generalizing the comprehensive presentation of the method for ODEs provided in the classical textbook by Golubitsky and Schaeffer [Applied mathematical sciences, {\bf 51}, Springer (1985)]. The appropriate setting in the context of DAEs at Hopf bifurcations is first detailed, introducing suitable operators and addressing the question of appropriate numerical algorithms for their construction as well. The different steps of the reduction procedure are carefully reinterpreted in the light of the DAE context and detailed formulas are provided for systematic and rational construction of the bifurcating local periodic solution, whose stability is shown, likely to the ODE context, to be predicted by the reduced equations. As an illustrative example, a classical DAE model for an electric power system is considered, exhibiting both supercritical and subcritical Hopf bifurcations, demonstrating the prediction capability of the reduced system with regard to the global dynamics.
Uwe Ehrenstein. Generalization to differential–algebraic equations of Lyapunov–Schmidt type reduction at Hopf bifurcations. Communications in Nonlinear Science and Numerical Simulation, 2024, 131, pp.107833. ⟨10.1016/j.cnsns.2024.107833⟩. ⟨hal-04408097⟩
Journal: Communications in Nonlinear Science and Numerical Simulation
Jingtao Ma, Lincheng Xu, Jérôme Jacob, Eric Serre, Pierre Sagaut. An averaged mass correction scheme for the simulation of high subsonic turbulent internal flows using a lattice Boltzmann method. Physics of Fluids, 2024, 36 (3), ⟨10.1063/5.0192360⟩. ⟨hal-04514161⟩ Plus de détails...
This paper addresses the simulation of internal high-speed turbulent compressible flows using lattice Boltzmann method (LBM) when it is coupled with the immersed boundary method for non-body-fitted meshes. The focus is made here on the mass leakage issue. The recent LBM pressure-based algorithm [Farag et al. Phys. Fluids 32, 066106 (2020)] has shown its superiority on classical density-based algorithm to simulate high-speed compressible flows. Following our previous theoretical work on incompressible flows [Xu et al. Phys. Fluids 34, 065113 (2022)], we propose an averaged mass correction technique to mitigate mass leakage when simulating high-Mach-number compressible flows. It is adapted to deal here with a density, which is decoupled from the zero-moment definition. The simulations focus on two generic but canonical configurations of more complex industrial devices, the straight channel at different angles of inclination at Mach numbers (Ma) ranging from 0.2 to 0.8, and the National Aeronautics and Space Administration Glenn S-duct at Ma = 0.6. The present results show that mass leakage can be a critical issue for the accuracy of the solution and that the proposed correction technique effectively mitigates it and leads to significant improvements in the prediction of the solution.
Jingtao Ma, Lincheng Xu, Jérôme Jacob, Eric Serre, Pierre Sagaut. An averaged mass correction scheme for the simulation of high subsonic turbulent internal flows using a lattice Boltzmann method. Physics of Fluids, 2024, 36 (3), ⟨10.1063/5.0192360⟩. ⟨hal-04514161⟩
Raffael Düll, Hugo Bufferand, Eric Serre, Guido Ciraolo, Virginia Quadri, et al.. Introducing electromagnetic effects in Soledge3X. Contributions to Plasma Physics, 2024, pp.e202300147. ⟨10.1002/ctpp.202300147⟩. ⟨hal-04474339⟩ Plus de détails...
In the pedestal region, electromagnetic effects affect the evolution of micro‐instabilities and plasma turbulence. The transport code Soledge3X developed by the CEA offers an efficient framework for turbulent 3D simulation on an electrostatic model with a fixed magnetic field. The physical accuracy of the model is improved with electromagnetic induction, driven by the local value of the parallel component of the electromagnetic vector potential , known from Ampère's law. It is solved implicitly in a coupled system with the vorticity equation on the electric potential . The consequence is a basic electromagnetic behavior in the form of shear Alfvén waves. A finite electron mass prevents unphysical speeds but requires solving for the time evolution of the parallel current density in the generalized Ohm's law. This term can be analytically included with little computational overhead in the system on and and improves its numerical condition, facilitating the iterative solving procedure. Simulations on a periodic slab case let us observe the predicted bifurcation of the wave propagation speed between the Alfvén wave and the electron thermal wave speeds for varying perpendicular wavenumbers. The first results on a circular geometry with a limiter attest to the feasibility of turbulent electromagnetic scenarios.
Raffael Düll, Hugo Bufferand, Eric Serre, Guido Ciraolo, Virginia Quadri, et al.. Introducing electromagnetic effects in Soledge3X. Contributions to Plasma Physics, 2024, pp.e202300147. ⟨10.1002/ctpp.202300147⟩. ⟨hal-04474339⟩
Frédéric Schwander, Eric Serre, Hugo Bufferand, Guido Ciraolo, Philippe Ghendrih. Global fluid simulations of edge plasma turbulence in tokamaks: a review. Computers and Fluids, 2024, 270, pp.106141. ⟨10.1016/j.compfluid.2023.106141⟩. ⟨hal-04352255⟩ Plus de détails...
With ITER, the largest tokamak ever built, and the growing number of fusion energy startups in the world, the need for numerical simulations has never been more crucial to progress towards the successful operation of fusion reactors. From fundamental plasma physics to engineering, a hierarchy of models exists from high-fidelity (gyro-)kinetic models in (5D) 6D to 0D fluid transport models. In this paper, we review the state-of-the-art of 3D turbulence fluid simulations in edge tokamak configurations. The widely used drift-reduced Braginskii equations are introduced together with the dedicated boundary conditions modelling plasma wall interactions. If until recently most of the models were focused on electrostatic turbulence driven by interchange-like instabilities, in recent years models have incorporated electromagnetic effects allowing fluctuations of the magnetic field. Specific features of the edge plasma configurations, which make these equations specially challenging to resolve and stressful for the numerical methods, are detailed. In particular, the strong anisotropy of the flow as well as the complex geometric characteristics lead to the development of dedicated discretization schemes and meshing, which are implemented in state-of-the-art codes reviewed here. It appears that the latter can be differentiated by their mesh construction as well by the manner in which they handle parallel gradients (aligned or not along the magnetic field). The review shows that no consensus on the optimal combination between meshing and discretization schemes, if it exists, has been found. Finally, examples of recent achievements show that 3D turbulence simulations of medium-sized tokamaks are currently achievable, but that ITERsize tokamaks and thermonuclear plasmas still require significant progress.
Frédéric Schwander, Eric Serre, Hugo Bufferand, Guido Ciraolo, Philippe Ghendrih. Global fluid simulations of edge plasma turbulence in tokamaks: a review. Computers and Fluids, 2024, 270, pp.106141. ⟨10.1016/j.compfluid.2023.106141⟩. ⟨hal-04352255⟩
Ivan Kudashev, Anna Medvedeva Glasser, Manuel Scotto D’abusco, Eric Serre. Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG. IEEE Transactions on Plasma Science, 2024, pp.1-6. ⟨10.1109/tps.2024.3384031⟩. ⟨hal-04552343⟩ Plus de détails...
Plasma–wall interaction is one of the key research topics on the way to controlled fusion. To study the best operational designs with reduced heat and particle fluxes onto tokamak plasma facing components (PFCs) comprehensive plasma simulations are required. A recent implementation of a hybridized discontinuous Galerkin scheme into a new version of SolEdge code has the advantage of using magnetic equilibrium-free mesh. This allows us to conduct pioneering 2-D transport simulations of a full discharge in the WEST tokamak. In this work, we implemented plasma transport coefficients as functions of coordinate in the poloidal plane and neutral diffusion as a function of neutral mean free path. Moreover, the perpendicular convection flux terms were added to the code. Using the new features, a few test cases were investigated. The influence of nonconstant transport coefficients on the simulated particle and heat fluxes onto the WEST tokamak PFCs are demonstrated.
Ivan Kudashev, Anna Medvedeva Glasser, Manuel Scotto D’abusco, Eric Serre. Impact of Variable Perpendicular Transport Coefficients in WEST Simulations Using SolEdge-HDG. IEEE Transactions on Plasma Science, 2024, pp.1-6. ⟨10.1109/tps.2024.3384031⟩. ⟨hal-04552343⟩
L. Cappelli, J. Guterl, N. Fedorczak, D.L. Rudakov, G. Sinclair, et al.. Model validation of tungsten erosion and redeposition properties using biased tungsten samples on DiMES. Nuclear Materials and Energy, 2023, 37, pp.101551. ⟨10.1016/j.nme.2023.101551⟩. ⟨hal-04546765⟩ Plus de détails...
An experiment was performed in the DIII-D lower divertor to validate numerical SOL tungsten (W) impurity erosion and redeposition simulations against experimental data. The net and gross erosion of W were calculated as a function of the voltage (or bias) applied to the exposed material. Five samples were inserted into the DIII-D lower divertor using the Divertor Material Evaluation System (DiMES) manipulator and exposed to constant L-mode attached plasma conditions. Each sample was partially coated with W. During plasma shots, samples were biased with respect to the machine vessel ground, ranging from −60 V to 25 V. The ERO2.0 code was used to numerically simulate the experiment aiming to compare the numerical results with experimental measures. A good agreement is found between estimated and measured tungsten erosion at least for negative biases.
L. Cappelli, J. Guterl, N. Fedorczak, D.L. Rudakov, G. Sinclair, et al.. Model validation of tungsten erosion and redeposition properties using biased tungsten samples on DiMES. Nuclear Materials and Energy, 2023, 37, pp.101551. ⟨10.1016/j.nme.2023.101551⟩. ⟨hal-04546765⟩
Hassan Hachem, Hai Canh Vu, Mitra Fouladirad. Different methods for RUL prediction considering sensor degradation. Reliability Engineering and System Safety, 2023, 243, pp.109897. ⟨10.1016/j.ress.2023.109897⟩. ⟨hal-04543787⟩ Plus de détails...
Predicting the Remaining Useful Lifetime (RUL) of a system has become one of the primary goals of engineering and reliability researchers. RUL prediction is based on the measurement data collected from sensors (e.g. vibration data, temperature data). The collected data is may be inaccurate owing to sensor problems. These problems are often ignored or modeled by a Gaussian noise in most previous work. However, due to various operation circumstances and the aging impact, the sensor itself will ultimately deteriorate and its performance will deteriorate. The Gaussian noise with a constant mean is then not appropriate to fully capture the sensor degradation. In this context, this study focuses on predicting the RUL considering the sensor degradation. For this purpose, a joint model of sensor degradation and system degradation is firstly developed. In this model, the sensor degradation is modeled by Wiener and Gamma processes instead of Gaussian noise. Then, different estimation methods based on the particle filter, a popular model-based technique, were proposed to predict the RUL based on the joint degradation model. To study the performances of our methods, numerical analyzes were carried out. The obtained results confirm the performance and advantages of the proposed methods.
Hassan Hachem, Hai Canh Vu, Mitra Fouladirad. Different methods for RUL prediction considering sensor degradation. Reliability Engineering and System Safety, 2023, 243, pp.109897. ⟨10.1016/j.ress.2023.109897⟩. ⟨hal-04543787⟩
Journal: Reliability Engineering and System Safety
Jingqi Zhang, Mitra Fouladirad, Nikolaos Limnios. A Semi-Markov Model with Geometric Renewal Processes. Methodology and Computing in Applied Probability, 2023, 25 (4), pp.85. ⟨10.1007/s11009-023-10060-z⟩. ⟨hal-04543367⟩ Plus de détails...
We consider a repairable system modeled by a semi-Markov process (SMP), where we include a geometric renewal process for system degradation upon repair, and replacement strategies for non-repairable failure or upon N repairs. First Pérez-Ocón and Torres-Castro studied this system (Pérez-Ocón and Torres-Castro in Appl Stoch Model Bus Ind 18(2):157–170, 2002) and proposed availability calculation using the Laplace Transform. In our work, we consider an extended state space for up and down times separately. This allows us to leverage the standard theory for SMP to obtain all reliability related measurements such as reliability, availability (point and steady-state), mean times and rate of occurrence of failures of the system with general initial law. We proceed with a convolution algebra, which allows us to obtain final closed form formulas for the above measurements. Finally, numerical examples are given to illustrate the methodology.
Jingqi Zhang, Mitra Fouladirad, Nikolaos Limnios. A Semi-Markov Model with Geometric Renewal Processes. Methodology and Computing in Applied Probability, 2023, 25 (4), pp.85. ⟨10.1007/s11009-023-10060-z⟩. ⟨hal-04543367⟩
Journal: Methodology and Computing in Applied Probability
Alex Karagrigoriou, Ioannis Mavrogiannis, Georgia Papasotiriou, Ilia Vonta. An Exponentiality Test of Fit Based on a Tail Characterization against Heavy and Light-Tailed Alternatives. Risks, 2023, 11 (10), pp.169. ⟨10.3390/risks11100169⟩. ⟨hal-04543899⟩ Plus de détails...
Log-concavity and log-convexity play a key role in various scientific fields, especially in those where the distinction between exponential and non-exponential distributions is necessary for inferential purposes. In the present study, we introduce a testing procedure for the tail part of a distribution which can be used for the distinction between exponential and non-exponential distributions. The conspiracy and catastrophe principles are initially used to establish a characterization of (the tail part of) the exponential distribution, which is one of the main contributions of the present work, leading the way for the construction of the new test of fit. The proposed test and its implementation are thoroughly discussed, and an extended simulation study has been undertaken to clarify issues related to its implementation and explore the extent of its capabilities. A real data case is also investigated.
Alex Karagrigoriou, Ioannis Mavrogiannis, Georgia Papasotiriou, Ilia Vonta. An Exponentiality Test of Fit Based on a Tail Characterization against Heavy and Light-Tailed Alternatives. Risks, 2023, 11 (10), pp.169. ⟨10.3390/risks11100169⟩. ⟨hal-04543899⟩
Mathis Pasquier, Stéphane Jay, Jérôme Jacob, Pierre Sagaut. A Lattice-Boltzmann-Based Modelling Chain for Traffic-Related Atmospheric Pollutant Dispersion at the Local Urban Scale. Building and Environment, 2023, 242, pp.110562. ⟨10.1016/j.buildenv.2023.110562⟩. ⟨hal-04190005⟩ Plus de détails...
Urban traffic-related air pollution is a major source of environmental and health damage and is difficult to quantify due to its inherent physical complexity. We construct a CFD-based simulation framework coupling an efficient numerical method for turbulent fluid flows with a microscopic traffic model and an emissions model to simulate road transport pollutant dispersion at the urban microscale. We improve the open-source Lattice-Boltzmann based CFD software OpenLB to overcome its original stability deficiencies for high Reynolds number flows. A stable recursive regularization procedure with a double distribution function approach is proposed to solve an advection diffusion equation for passive scalar transport at high Reynolds number. The code is successfully validated on three reference cases of increasing complexity and the traffic model SUMO along with a physical engine emissions model are coupled with OpenLB to simulate traffic-induced pollution from a road network in a realistic complex geometry. Transient flow features are analysed and the time-averaged concentration levels in different neighbourhoods of the considered geometry are evaluated: high concentration levels are observed close to the streets but also inside specific building infrastructures due to complex wind dynamics. Analyses of altitudinal concentration variations show that flow recirculations located close to traffic lights can drive pollutant over the buildings and increase concentration levels inside inner courtyards. Time-averaged concentration maps are constructed using both spatially uniform and non-uniform line sources and it is shown that using uniform sources leads to up to 20% local concentration overestimations inside the urban canopy.
Mathis Pasquier, Stéphane Jay, Jérôme Jacob, Pierre Sagaut. A Lattice-Boltzmann-Based Modelling Chain for Traffic-Related Atmospheric Pollutant Dispersion at the Local Urban Scale. Building and Environment, 2023, 242, pp.110562. ⟨10.1016/j.buildenv.2023.110562⟩. ⟨hal-04190005⟩
Isabelle Cheylan, Tom Fringand, Jérôme Jacob, Julien Favier. Analysis of the immersed boundary method for turbulent fluid-structure interaction with Lattice Boltzmann method. Journal of Computational Physics, 2023, 492, pp.112418. ⟨10.1016/j.jcp.2023.112418⟩. ⟨hal-04543923⟩ Plus de détails...
L. Cappelli, N. Fedorczak, J. P. Gunn, S. Di Genova, J. Guterl, et al.. Study of the erosion and redeposition of W considering the kinetic energy distribution of incident ions through a semi-analytical model. Plasma Physics and Controlled Fusion, 2023, 65 (9), pp.095001. ⟨10.1088/1361-6587/ace282⟩. ⟨hal-04190861⟩ Plus de détails...
In today’s nuclear fusion devices, erosion of high-Z metallic plasma-facing materials (PFMs) is mainly caused by physical sputtering. That is, by the exchange of energy between plasma ions and the atoms in the walls. In most of the numerical codes currently in use impinging plasma is approximated as a fluid. By averaging the incident particles’ energy distribution the high-energy population of the eroded material is underestimated. For heavy materials such as W, high-energy eroded particles tend to ionize far from the wall and they are less affected by the sheath electric field hence, not being attracted back to the wall, they have a higher chance to contaminate the core plasma. This could in turn result in an underestimation of the net erosion sources. In this work, a semi-analytical model was developed to include the energy distribution of the incident particles. Then, by Monte Carlo method, the net erosion of tungsten from a smooth PFM was calculated. The results show that the kinetic description in energy is important only for incident particles ionized once. For instance, it is particularly important for plasma ions such as Deuterium. It is seen that Deuterium contribution to the W net sources is not always negligible if compared to light impurities or to tungsten self-sputtering in the range of plasma parameters tested. Finally, results show that the difference between the fluid and kinetic models becomes more pronounced for high-screening plasma conditions.
L. Cappelli, N. Fedorczak, J. P. Gunn, S. Di Genova, J. Guterl, et al.. Study of the erosion and redeposition of W considering the kinetic energy distribution of incident ions through a semi-analytical model. Plasma Physics and Controlled Fusion, 2023, 65 (9), pp.095001. ⟨10.1088/1361-6587/ace282⟩. ⟨hal-04190861⟩
Chenglei Wang, Simon Gsell, Umberto d'Ortona, Julien Favier. Generalized-Newtonian fluid transport by an instability-driven filament. Journal of Fluid Mechanics, 2023, 965, pp.A6. ⟨10.1017/jfm.2023.381⟩. ⟨hal-04544010⟩ Plus de détails...
Cilia are micro-scale hair-like organelles. They can exhibit self-sustained oscillations which play crucial roles in flow transport or locomotion. Recent studies have shown that these oscillations can spontaneously emerge from dynamic instability triggered by internal stresses via a Hopf bifurcation. However, the flow transport induced by an instability-driven cilium still remains unclear, especially when the fluid is non-Newtonian. This study aims at bridging these gaps. Specifically, the cilium is modelled as an elastic filament, and its internal actuation is represented by a constant follower force imposed at its tip. Three generalized Newtonian behaviours are considered, i.e. the shear-thinning, Newtonian and shear-thickening behaviours. Effects of four key factors, including the filament zero-stress shape, Reynolds number ( $Re$ ), follower-force magnitude and fluid rheology, on the filament dynamics, fluid dynamics and flow transport are explored through direct numerical simulation at $Re$ of 0.04 to 5 and through a scaling analysis at $Re \approx 0$ . The results reveal that even though it is expected that inertia vanishes at $Re \ll 1$ , inertial forces do alter the filament dynamics and deteriorate the flow transport at $Re\ge 0.04$ . Regardless of $Re$ , the flow transport can be improved when the flow is shear thinning or when the follower force increases. Furthermore, a linear stability analysis is performed, and the variation of the filament beating frequency, which is closely correlated with the filament dynamics and flow transport, can be predicted.
Chenglei Wang, Simon Gsell, Umberto d'Ortona, Julien Favier. Generalized-Newtonian fluid transport by an instability-driven filament. Journal of Fluid Mechanics, 2023, 965, pp.A6. ⟨10.1017/jfm.2023.381⟩. ⟨hal-04544010⟩
E.V. Kuidjo Kuidjo, M.G. Rodio, R. Abgrall, P. Sagaut. Comparison of bubbles interaction mechanisms of two-group Interfacial Area Transport Equation model. International Journal of Multiphase Flow, 2023, 163, pp.104399. ⟨10.1016/j.ijmultiphaseflow.2023.104399⟩. ⟨cea-04483130⟩ Plus de détails...
This work deals with 3D simulations of complex bubbly, cap-bubbly and churn regimes exhibiting bubbles of different shapes and with broad bubble size distribution. The first contribution of this work is to investigate and compare several bubble interaction mechanisms of coalescence and fragmentation for the 2-Group Interfacial Area Transport Equation (IATE) model. For two of these models, this is the first time their performances are assessed within a CFD code. The second contribution is to propose and assess a novel model of fragmentation and coalescence. Finally, a validation versus experimental data on three different configurations and three different regimes is performed. In particular, the interaction mechanisms are analysed for one specific regime. The implementation of the two-group IATE model has been systematically performed in the 3D NEPTUNE CFD code.
E.V. Kuidjo Kuidjo, M.G. Rodio, R. Abgrall, P. Sagaut. Comparison of bubbles interaction mechanisms of two-group Interfacial Area Transport Equation model. International Journal of Multiphase Flow, 2023, 163, pp.104399. ⟨10.1016/j.ijmultiphaseflow.2023.104399⟩. ⟨cea-04483130⟩
Shang-Gui Cai, Jérôme Jacob, Pierre Sagaut. Immersed boundary based near-wall modeling for large eddy simulation of turbulent wall-bounded flow. Computers and Fluids, 2023, 259, pp.105893. ⟨10.1016/j.compfluid.2023.105893⟩. ⟨hal-04543910⟩ Plus de détails...
Shang-Gui Cai, Jérôme Jacob, Pierre Sagaut. Immersed boundary based near-wall modeling for large eddy simulation of turbulent wall-bounded flow. Computers and Fluids, 2023, 259, pp.105893. ⟨10.1016/j.compfluid.2023.105893⟩. ⟨hal-04543910⟩
E.V. Kuidjo Kuidjo, M.G. Rodio, R. Abgrall, P. Sagaut. Comparison of bubbles interaction mechanisms of two-group Interfacial Area Transport Equation model. International Journal of Multiphase Flow, 2023, 163, pp.104399. ⟨10.1016/j.ijmultiphaseflow.2023.104399⟩. ⟨hal-04543708⟩ Plus de détails...
This work deals with 3D simulations of complex bubbly, cap-bubbly and churn regimes exhibiting bubbles of different shapes and with broad bubble size distribution. The first contribution of this work is to investigate and compare several bubble interaction mechanisms of coalescence and fragmentation for the 2-Group Interfacial Area Transport Equation (IATE) model. For two of these models, this is the first time their performances are assessed within a CFD code. The second contribution is to propose and assess a novel model of fragmentation and coalescence. Finally, a validation versus experimental data on three different configurations and three different regimes is performed. In particular, the interaction mechanisms are analysed for one specific regime. The implementation of the two-group IATE model has been systematically performed in the 3D NEPTUNE CFD code.
E.V. Kuidjo Kuidjo, M.G. Rodio, R. Abgrall, P. Sagaut. Comparison of bubbles interaction mechanisms of two-group Interfacial Area Transport Equation model. International Journal of Multiphase Flow, 2023, 163, pp.104399. ⟨10.1016/j.ijmultiphaseflow.2023.104399⟩. ⟨hal-04543708⟩
Minh Nguyen, Jean-François Boussuge, Pierre Sagaut, Juan-Carlos Larroya-Huguet. Large eddy simulation of a row of impinging jets with upstream crossflow using the lattice Boltzmann method. International Journal of Heat and Mass Transfer, 2023, 212, pp.124256. ⟨10.1016/j.ijheatmasstransfer.2023.124256⟩. ⟨hal-04546755⟩ Plus de détails...
Large eddy simulations of a row of seven jets emerging from a perforated pipe and impinging on a flat heated plate were carried out using a compressible hybrid thermal lattice Boltzmann solver. The average Reynolds number of the emerging jets was , with an exit-to-plate distance of 3 jet diameters. Two levels of upstream crossflow were simulated: one with weak cross flow, with a velocity ratio of , and one with strong cross flow, with a velocity ratio of . The flow field and heat transfer statistics were validated against experimental PIV and infrared thermography data from a recent study, showing good agreement. The effect of varying the Mach number on the wall heat transfer was subsequently tested. It was found that an increased Mach number lead to an increased value of the Nusselt number near the stagnation points.
Minh Nguyen, Jean-François Boussuge, Pierre Sagaut, Juan-Carlos Larroya-Huguet. Large eddy simulation of a row of impinging jets with upstream crossflow using the lattice Boltzmann method. International Journal of Heat and Mass Transfer, 2023, 212, pp.124256. ⟨10.1016/j.ijheatmasstransfer.2023.124256⟩. ⟨hal-04546755⟩
Journal: International Journal of Heat and Mass Transfer
G. Farag, P. Boivin, P. Sagaut. Linear interaction approximation for shock/disturbance interaction in a Noble–Abel stiffened gas. Shock Waves, 2023, ⟨10.1007/s00193-023-01131-8⟩. ⟨hal-04097657⟩ Plus de détails...
When departure from the ideal gas equation of state is considered, the Noble-Abel stiffened gas model is an appealing and versatile candidate due to its simple form. The Linear Interaction Approximation formalism is extended to consider non-ideal gas effects introduced by this equation of state. Kovásznay decomposition and adequate definition of the energy of disturbances are provided in the context of this equation of state. Changes with respect to ideal gas are investigated on transfer functions, critical angle and compression factor. Those differences yield concrete effects on the damping and transfer of fluctuations across shock waves. Those changes are further illustrated by considering the interaction of an entropy spot with a Mach 3 stationary shock wave.
G. Farag, P. Boivin, P. Sagaut. Linear interaction approximation for shock/disturbance interaction in a Noble–Abel stiffened gas. Shock Waves, 2023, ⟨10.1007/s00193-023-01131-8⟩. ⟨hal-04097657⟩
Jérémie Labasse, Uwe Ehrenstein, Guillaume Fasse, Frédéric Hauville. Thrust scaling for a large-amplitude heaving and pitching foil with application to cycloidal propulsion. Ocean Engineering, 2023, 275, pp.114169. ⟨10.1016/j.oceaneng.2023.114169⟩. ⟨hal-04032117⟩ Plus de détails...
A numerical solution procedure using the mesh-superposition approach, known as the Chimera method, together with the OpenFOAM toolbox environment is used to compute the forces generated by large amplitude heaving and pitching foil. The possibility of fitting thrust prediction laws, based on classical potential flow theories, with the numerically computed forces is explored, for a Reynolds number of 5 10 4. It is shown, first for a pure heaving motion and subsequently by adding a harmonic pitching motion, that theoretical scaling may be fitted to numerical time-averaged thrust data, even in the case of large amplitude motions. The thrust-prediction law is shown to still apply to pitching-rotating motions, such as those of blades in cycloidal propulsion devices, the mean pressure correction due to the additional surging motion being small. The synchronized rotation-pitching of three foils typical of a cross-flow propeller configuration is addressed as well. The numerical global thrust results are shown to be in general agreement with the theoretical prediction, but also with blade-embedded load cell measurements for an experimental device developed by the French Naval Academy Research Institute.
Jérémie Labasse, Uwe Ehrenstein, Guillaume Fasse, Frédéric Hauville. Thrust scaling for a large-amplitude heaving and pitching foil with application to cycloidal propulsion. Ocean Engineering, 2023, 275, pp.114169. ⟨10.1016/j.oceaneng.2023.114169⟩. ⟨hal-04032117⟩
Hasan Misaii, Mitra Fouladirad, Firoozeh Haghighi. Optimal task-driven time-dependent covariate-based maintenance policy. Journal of Computational and Applied Mathematics, 2023, pp.115315. ⟨10.1016/j.cam.2023.115315⟩. ⟨hal-04162527⟩ Plus de détails...
Alexandre Di-Marco, Jerome Jacob, Pierre Sagaut. Unsteady characteristics of pressure and swirl distortion in helicopter intake: A lattice Boltzmann method approach. Aerospace Science and Technology, 2023, 138, pp.108333. ⟨10.1016/j.ast.2023.108333⟩. ⟨hal-04543807⟩ Plus de détails...
The present paper focuses on the analysis of the unsteady characteristics of pressure and swirl distortion and their relationship inside a complex air intake, including a plenum chamber, using Large-Eddy Simulation based on the lattice Boltzmann method (LBM). The steady-state analysis of the pressure field and swirl angle in the Aerodynamic Interface Plane (AIP) shows a complex pressure and swirl distortion pattern, as well as a high level of unsteadiness. Additionally, two vortical structures with equal intensity located on the top and bottom sides of the AIP plane were identified. The unsteady pressure distortion analysis highlighted fluctuations between two positions of the maximum circumferential pressure distortion in the AIP plane. One of these positions presents mean and peak distortion values lower than the other position. Furthermore, an alternating vortex pattern between two swirl patterns in the AIP plane is also identified. Finally, the relationship between pressure and swirl unsteadiness was investigated, and a link between the swirl distortion pattern and pressure distortion was identified. The identification of this relationship will allow geometrical changes in the air intake, leading to a reduction of the average circumferential distortion level, as well as a reduction of the maximum peak level achievable during a flight phase.
Alexandre Di-Marco, Jerome Jacob, Pierre Sagaut. Unsteady characteristics of pressure and swirl distortion in helicopter intake: A lattice Boltzmann method approach. Aerospace Science and Technology, 2023, 138, pp.108333. ⟨10.1016/j.ast.2023.108333⟩. ⟨hal-04543807⟩
Thomas Gianoli, Jean‐françois Boussuge, Pierre Sagaut, Jérôme de Laborderie. Development and validation of Navier–Stokes characteristic boundary conditions applied to turbomachinery simulations using the lattice Boltzmann method. International Journal for Numerical Methods in Fluids, 2023, 95 (4), pp.528-556. ⟨10.1002/fld.5160⟩. ⟨hal-04063964⟩ Plus de détails...
This article reports a procedure to implement as well as to validate non-reflecting boundary conditions applied for turbomachinery simulations, using Navier-Stokes characteristic boundary conditions in a compressible lattice Boltzmann solver. The implementation of both an inlet condition imposing total pressure, total temperature, and flow angles, as well as an outlet condition imposing a static pressure profile that allows the simulation to reach a simplified radial equilibrium, is described within the context of a lattice Boltzmann approach. The treatment at the boundaries relies on the characteristic methodology to derive conditions which are non-reflecting in terms of acoustics and is also compatible with turbulence injection at the inlet. These properties are evaluated on test cases of increasing complexity, ranging from a simple 2D periodic domain to an S-duct stage with turbulence injection.
Thomas Gianoli, Jean‐françois Boussuge, Pierre Sagaut, Jérôme de Laborderie. Development and validation of Navier–Stokes characteristic boundary conditions applied to turbomachinery simulations using the lattice Boltzmann method. International Journal for Numerical Methods in Fluids, 2023, 95 (4), pp.528-556. ⟨10.1002/fld.5160⟩. ⟨hal-04063964⟩
Journal: International Journal for Numerical Methods in Fluids
Imran Afgan, Yacine Kahil, Sofiane Benhamadouche, Mohamed Ali, Ahmed Alkaabi, et al.. Cross flow over two heated cylinders in tandem arrangements at subcritical Reynolds number using large eddy simulations. International Journal of Heat and Fluid Flow, 2023, 100, pp.109115. ⟨10.1016/j.ijheatfluidflow.2023.109115⟩. ⟨hal-04546846⟩ Plus de détails...
This study analyses the heat transfer and flow characteristics of cross-flow over two heated infinite cylinders in a tandem (in-line) configuration. Non-isothermal Large Eddy Simulations (LES) using the dynamic Smagorinsky model were conducted at a fixed Reynolds number of 3, 000 (based on the free stream velocity and the cylinder diameter). A range of cylinder gap ratios (1.0 ≤ L/D ≤ 5.0) was investigated (in increments of 0.25) with two different Prandtl numbers Pr = 0.1 and 1.0. Results show that the flow structures vary according to the order of the patterns: (i) Extended body regime: without attachment for low L/D (1.0 1.25) where cylinders behave as a single bluff body with top–bottom vortex shedding, (ii) Shear layer reattachment regime: with reattachment for moderate L/D (1.5 3.75) where the detached shear layer from the upstream cylinder reattaches to the down- stream cylinder, and (iii) Co-shedding regime: for high gap ratios (3.75 ≤ L/D ≤ 5.0) a phenomenon called “jumping”, where the two cylinders behave as isolated bluff bodies. Furthermore, it was observed that the average Nusselt number of both cylinders experience a drastic variation at a critical spacing ratio (between 3.75 ≤ L/D ≤ 4.0). For L/D ≤ 3.0, the average Nusselt number of the upstream cylinder was found to be higher than that of the downstream one. However, for spacing ratios L/D > 3.0, the average Nusselt number was similar for both cylinders. For the downstream cylinder, the maximum Nusselt number was located at the separation angle and was found to be independent of the spacing ratio
Imran Afgan, Yacine Kahil, Sofiane Benhamadouche, Mohamed Ali, Ahmed Alkaabi, et al.. Cross flow over two heated cylinders in tandem arrangements at subcritical Reynolds number using large eddy simulations. International Journal of Heat and Fluid Flow, 2023, 100, pp.109115. ⟨10.1016/j.ijheatfluidflow.2023.109115⟩. ⟨hal-04546846⟩
Journal: International Journal of Heat and Fluid Flow
Richard M Lueptow, Rainer Hollerbach, Eric Serre. Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper: part 1. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2243), ⟨10.1098/rsta.2022.0140⟩. ⟨hal-03989060⟩ Plus de détails...
In 1923, the Philosophical Transactions published G. I. Taylor’s seminal paper on the stability of what we now call Taylor–Couette flow. In the century since the paper was published, Taylor’s ground-breaking linear stability analysis of fluid flow between two rotating cylinders has had an enormous impact on the field of fluid mechanics. The paper’s influence has extended to general rotating flows, geophysical flows and astrophysical flows, not to mention its significance in firmly establishing several foundational concepts in fluid mechanics that are now broadly accepted. This two-part issue includes review articles and research articles spanning a broad range of contemporary research areas, all rooted in Taylor’s landmark paper. This article is part of the theme issue ‘Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper (part 1)’.
Richard M Lueptow, Rainer Hollerbach, Eric Serre. Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper: part 1. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2243), ⟨10.1098/rsta.2022.0140⟩. ⟨hal-03989060⟩
Journal: Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Rainer Hollerbach, Richard M Lueptow, Eric Serre. Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper: part 2. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2246), ⟨10.1098/rsta.2022.0359⟩. ⟨hal-04057399⟩ Plus de détails...
In 1923, the Philosophical Transactions published G. I. Taylor’s seminal paper on the stability of what we now call Taylor–Couette flow. In the century since the paper was published, Taylor’s ground-breaking linear stability analysis of fluid flow between two rotating cylinders has had an enormous impact on the field of fluid mechanics. The paper’s influence has extended to general rotating flows, geophysical flows and astrophysical flows, not to mention its significance in firmly establishing several foundational concepts in fluid mechanics that are now broadly accepted. This two-part issue includes review articles and research articles spanning a broad range of contemporary research areas, all rooted in Taylor’s landmark paper. This article is part of the theme issue ‘Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper (part 2)’.
Rainer Hollerbach, Richard M Lueptow, Eric Serre. Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper: part 2. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2246), ⟨10.1098/rsta.2022.0359⟩. ⟨hal-04057399⟩
Journal: Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
S. Di Genova, G. Ciraolo, A. Gallo, J. Romazanov, N. Fedorczak, et al.. First 3D modeling of tungsten erosion and migration in WEST discharges adopting a toroidally non-symmetric wall geometry. Nuclear Materials and Energy, 2023, 34, pp.101340. ⟨10.1016/j.nme.2022.101340⟩. ⟨hal-03988791⟩ Plus de détails...
Numerical analyses are a key tool to investigate tungsten (W) sources and contamination in W Environment steady-state tokamak (WEST) plasma discharges. Modelling activity was performed in order to study W erosion and migration at WEST plasma-facing components (PFCs), using for the first time a toroidally asymmetric wall geometry provided by toroidally localized objects representing WEST outer limiter or antennae. 3D non-axisymmetric SOLEDGE transport simulations were performed with simplifying assumptions (pure Deuterium plasma, fluid model for neutrals) to reproduce WEST boundary plasma, and used as background for ERO2.0 simulations modelling W erosion, re-deposition, and migration. On the sides of the toroidally localized limiters/antennae, two thin W stripes were considered in order to model WEST W antennae protections. Simulations suggest antennae protections contribution to dominate W contamination in the considered simulations settings, highlighting the need of further analyses with different configurations using this kind of tools.
S. Di Genova, G. Ciraolo, A. Gallo, J. Romazanov, N. Fedorczak, et al.. First 3D modeling of tungsten erosion and migration in WEST discharges adopting a toroidally non-symmetric wall geometry. Nuclear Materials and Energy, 2023, 34, pp.101340. ⟨10.1016/j.nme.2022.101340⟩. ⟨hal-03988791⟩
Ivan Kudashev, Anna Medvedeva, Nicolas Fedorszak, David Zarzoso, Manuel Scotto d’ Abusco, et al.. Development of a set of synthetic diagnostics for the WEST tokamak to confront 2D transport simulations and experimental data. Journal of Instrumentation, 2023, 18 (02), pp.C02058. ⟨10.1088/1748-0221/18/02/C02058⟩. ⟨hal-04010344⟩ Plus de détails...
Significant scientific effort has been focused on optimizing the scenarios and plasma parameters for tokamak operations. The lack of comprehensive understanding of underlying physical processes leads to simplifications used both in plasma simulation codes and for diagnostics, which is also complicated by the harsh plasma environment. One of the main tools to couple, check and verify these assumptions are the synthetic diagnostics. In this work we demonstrate current results of the development of the set of synthetic diagnostics for the WEST tokamak to couple experimental data with the SolEdge3X-HDG 2D transport code.
Ivan Kudashev, Anna Medvedeva, Nicolas Fedorszak, David Zarzoso, Manuel Scotto d’ Abusco, et al.. Development of a set of synthetic diagnostics for the WEST tokamak to confront 2D transport simulations and experimental data. Journal of Instrumentation, 2023, 18 (02), pp.C02058. ⟨10.1088/1748-0221/18/02/C02058⟩. ⟨hal-04010344⟩
Thomas Cartier-Michaud, Philippe Ghendrih, Virginie Grandgirard, Eric Serre. Verification and accuracy check of simulations with PoPe and iPoPe. Journal of Computational Physics, 2023, 474, pp.111759. ⟨10.1016/j.jcp.2022.111759⟩. ⟨hal-03871954⟩ Plus de détails...
The theoretical background of the PoPe and iPoPe verification scheme is presented. Verification is performed using the simulation output of production runs. The computing overhead is estimated to be at most 10%. PoPe or iPoPe calculations can be done offline provided the necessary data is stored, for example additional time slices, or online where iPoPe is more effective. The computing overhead is mostly that of storing the necessary data. The numerical error is determined and split into a part proportional to the operators, which are combined to form the equations to be solved, thus modifying their control parameters, completed by a residual error orthogonal to these operators. The accuracy of the numerical solution is determined by this modification of the control parameters. The PoPe and iPoPe methods are illustrated in this paper with simulations of a simple mechanical system with chaotic trajectories evolving into a strange attractor with sensitivity to initial conditions. We show that the accuracy depends on the particular simulation both because the properties of the numerical solution depend on the values of the control parameter, and because the target accuracy will depend on the problem that is addressed. One shows that for a case close to bifurcations between different states, the accuracy is determined by the level of detail of the bifurcation phenomena one aims at describing. A unique verification index, the PoPe index, is proposed to characterise the accuracy, and consequently the verification, of each production run. The PoPe output allows one to step beyond verification and analyse for example the numerical scheme efficiency. For the chosen example at fixed PoPe index, therefore at fixed numerical error, one finds that the higher order integration scheme, comparing order 4 to order 2 Runge-Kutta time stepping, reduces the computation cost by a factor 4.
Thomas Cartier-Michaud, Philippe Ghendrih, Virginie Grandgirard, Eric Serre. Verification and accuracy check of simulations with PoPe and iPoPe. Journal of Computational Physics, 2023, 474, pp.111759. ⟨10.1016/j.jcp.2022.111759⟩. ⟨hal-03871954⟩
Gauthier Wissocq, Said Taileb, Song Zhao, Pierre Boivin. A hybrid lattice Boltzmann method for gaseous detonations. Journal of Computational Physics, 2023, 494, pp.112525. ⟨10.1016/j.jcp.2023.112525⟩. ⟨hal-04244340⟩ Plus de détails...
This article is dedicated to the construction of a robust and accurate numerical scheme based on the lattice Boltzmann method (LBM) for simulations of gaseous detonations. This objective is achieved through careful construction of a fully conservative hybrid lattice Boltzmann scheme tailored for multi-species reactive flows. The core concept is to retain LBM low dissipation properties for acoustic and vortical modes by using the collide and stream algorithm for the particle distribution function, while transporting entropic and species modes via a specifically designed finite-volume scheme. The proposed method is first evaluated on common academic cases, demonstrating its ability to accurately simulate multi-species compressible and reactive flows with discontinuities: the convection of inert species, a Sod shock tube with two ideal gases and a steady one-dimensional inviscid detonation wave. Subsequently, the potential of this novel approach is demonstrated in one- and two-dimensional inviscid unsteady gaseous detonations, highlighting its ability to accurately recover detonation structures and associated instabilities for high activation energies. To the authors' knowledge, this study is the first successful simulation of detonation cellular structures capitalizing on the LBM collide and stream algorithm.
Gauthier Wissocq, Said Taileb, Song Zhao, Pierre Boivin. A hybrid lattice Boltzmann method for gaseous detonations. Journal of Computational Physics, 2023, 494, pp.112525. ⟨10.1016/j.jcp.2023.112525⟩. ⟨hal-04244340⟩
Jingtao Ma, Qiuxiang Huang, Yi Zhu, Yuan-Qing Xu, Fang-Bao Tian. Effects of fluid rheology on dynamics of a capsule through a microchannel constriction. Physics of Fluids, 2023, 35 (9), pp.091901. ⟨10.1063/5.0165614⟩. ⟨hal-04546805⟩ Plus de détails...
This paper numerically investigates the impact of fluid rheology on the behaviors of a spherical capsule through a microchannel constriction. Different flow scenarios are considered: a Newtonian capsule in a viscoelastic matrix, a Newtonian capsule in a Newtonian matrix, and a viscoelastic capsule in a Newtonian matrix. The results demonstrate that the capsule's lengths undergo oscillations during the passage through the constriction, with three stages of evolution. When approaching the constriction, the capsule respectively experiences increase and decrease in its length and height. While within or exiting the constriction, the length of the capsule continuously decreases, and the height generally increases. As the capsule moves away from the constriction, the capsule relaxes to different profiles in different flows. Detailed analysis on the effects of the fluid viscoelasticity on the capsule's lengths in different stages is provided. In addition, the behaviors of a red blood cell passing through a microchannel constriction are also examined. This study sheds light on the complex behaviors of a spherical capsule and red blood cell in microchannel constriction, emphasizing the significant influence of fluid rheology on their deformation and shape changes.
Jingtao Ma, Qiuxiang Huang, Yi Zhu, Yuan-Qing Xu, Fang-Bao Tian. Effects of fluid rheology on dynamics of a capsule through a microchannel constriction. Physics of Fluids, 2023, 35 (9), pp.091901. ⟨10.1063/5.0165614⟩. ⟨hal-04546805⟩
Franck Corset, Mitra Fouladirad, Christian Paroissin. Imperfect condition-based maintenance for a gamma degradation process in presence of unknown parameters. Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability, 2023, 237 (3), pp.546-561. ⟨10.1177/1748006X221134132⟩. ⟨hal-03842177⟩ Plus de détails...
A system subject to degradation is considered. The degradation is modelled by a gamma process. A condition-based maintenance policy with perfect corrective and an imperfect preventive actions is proposed. The maintenance cost is derived considering a Markov-renewal process. The statistical inference of the degradation and maintenance parameters by the maximum likelihood method is investigated. A sensibility analysis to different parameters is carried out and the perspectives are detailed.
Franck Corset, Mitra Fouladirad, Christian Paroissin. Imperfect condition-based maintenance for a gamma degradation process in presence of unknown parameters. Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability, 2023, 237 (3), pp.546-561. ⟨10.1177/1748006X221134132⟩. ⟨hal-03842177⟩
Journal: Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability
D. Martinand, E. Serre, B. Viaud. Instabilities and routes to turbulence in rotating disc boundary layers and cavities. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2243), pp.20220135. ⟨10.1098/rsta.2022.0135⟩. ⟨hal-03989074⟩ Plus de détails...
Studied for more than a century, first in the field of geophysics, flows over rotating discs present a great diversity of complex instability behaviours that are not yet fully understood. While the primary instabilities are now well characterized experimentally, theoretically and numerically, their role in the transition mechanisms to turbulence remains an open question that still challenges the scientific community. This article brings together the main results of the literature related to the instabilities over rotating discs, but also in the connected problem of rotating cavities, and reviews the main scenarios currently assumed to describe the flow breakdown to turbulence. A particular focus is on more recent studies of generic flows in rotating cavities bounded by two coaxial rotating discs, that occur in many industrial systems, the performances of which and their improvement are linked to a better understanding of these mechanisms. This article is part of the theme issue ‘Taylor–Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (part 1)’.
D. Martinand, E. Serre, B. Viaud. Instabilities and routes to turbulence in rotating disc boundary layers and cavities. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2023, 381 (2243), pp.20220135. ⟨10.1098/rsta.2022.0135⟩. ⟨hal-03989074⟩
Journal: Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Pierre Sagaut, V.K. Suman, P. Sundaram, M.K. Rajpoot, Y.G. Bhumkar, et al.. Global spectral analysis: Review of numerical methods. Computers and Fluids, 2023, 261, pp.105915. ⟨10.1016/j.compfluid.2023.105915⟩. ⟨hal-04546492⟩ Plus de détails...
The design and analysis of numerical methods are usually guided by the following: (a) von Neumann analysis using Fourier series expansion of unknowns, (b) the modified differential equation approach, and (c) a more generalized approach that analyzes numerical methods globally, using Fourier–Laplace transform to treat the total or disturbance quantities in terms of waves. This is termed as the global spectral analysis (GSA). GSA can easily handle non-periodic problems, by invoking wave properties of the field through the correct numerical dispersion relation, which is central to the design and analysis. This has transcended dimensionality of the problem, while incorporating various physical processes e.g. by studying convection, diffusion and reaction as the prototypical elements involved in defining the physics of the problem. Although this is used for fluid dynamical problems, it can also explain many multi-physics and multi-scale problems. This review describes this powerful tool of scientific computing, with new results originating from GSA: (i) providing a common framework to analyze both hyperbolic and dispersive wave problems; (ii) analyze numerical methods by comparing physical and numerical dispersion relation, which leads to the new class of dispersion relation preserving (DRP) schemes; (iii) developing error dynamics as a distinct tool, identifying sources of numerical errors involving both the truncation and round-off error. Such studies of error dynamics provide the epistemic tool of analysis rather than an aleatoric tool, which depends on uncertainty quantification for high performance computing (HPC). One of the central themes of GSA covers the recent advances in understanding numerical phenomenon like focusing, which defied analysis so far. An application of GSA shown here for the objective evaluation of the so-called DNS by pseudo-spectral method for spatial discretization along with time integration by two-stage Runge–Kutta method is performed. GSA clearly shows that this should not qualify as DNS for multiple reasons. A new design of HPC methods for peta- and exa-flop computing tools necessary for parallel computing by compact schemes are also described.
Pierre Sagaut, V.K. Suman, P. Sundaram, M.K. Rajpoot, Y.G. Bhumkar, et al.. Global spectral analysis: Review of numerical methods. Computers and Fluids, 2023, 261, pp.105915. ⟨10.1016/j.compfluid.2023.105915⟩. ⟨hal-04546492⟩
Jacob Johnston, Sarah Dischinger, Mostafa Nassr, Ji Yeon Lee, Pedram Bigdelou, et al.. A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers. Journal of Membrane Science, 2023, 675, pp.121508. ⟨10.1016/j.memsci.2023.121508⟩. ⟨hal-04546833⟩ Plus de détails...
Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. We also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.
Jacob Johnston, Sarah Dischinger, Mostafa Nassr, Ji Yeon Lee, Pedram Bigdelou, et al.. A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers. Journal of Membrane Science, 2023, 675, pp.121508. ⟨10.1016/j.memsci.2023.121508⟩. ⟨hal-04546833⟩
Hasan Misaii, Firoozeh Haghighi, Mitra Fouladirad. Opportunistic perfect preventive maintenance policy in presence of masked data. Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability, 2022, 236 (6), pp.1024-1036. ⟨10.1177/1748006X211058936⟩. ⟨hal-04064544⟩ Plus de détails...
In this paper, the maintenance optimization problem of multi-component system is considered. It is assumed that the exact cause of system failure might be masked. That is, the exact cause of failure is unknown, and we only know that it belongs to a set called mask set. Both opportunistic perfect preventive maintenance (OPPM) and perfect corrective maintenance are considered. Threshold of OPPM and inter-inspection interval are considered as decision parameters which are optimized using long-run cost rate criteria. The applicability of the proposed maintenance policy is investigated using an illustrative example.
Hasan Misaii, Firoozeh Haghighi, Mitra Fouladirad. Opportunistic perfect preventive maintenance policy in presence of masked data. Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability, 2022, 236 (6), pp.1024-1036. ⟨10.1177/1748006X211058936⟩. ⟨hal-04064544⟩
Journal: Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability
Philippe Ghendrih, Guilhem Dif-Pradalier, Olivier Panico, Yanick Sarazin, Hugo Bufferand, et al.. Role of avalanche transport in competing drift wave and interchange turbulence. Journal of Physics: Conference Series, 2022, 2397, pp.012018. ⟨10.1088/1742-6596/2397/1/012018⟩. ⟨hal-04498387⟩ Plus de détails...
We complete the 2D 2-fields turbulence model previously used with an interchange-like instability by slightly modifying the parallel loss terms to drive drift wave instabilities. We show that the instability driven by temperature fluctuations of the sheath losses is identical to that of the drift wave turbulence. The linear analysis is performed and used to select control parameters that yield identical maximum growth rates for the interchange alone and drift wave alone instability. Combining the two instabilities doubles the maximum growth rate. The non-linear simulations are used to analyse the SOL width. The simulations allow one to identify a low field side SOL region where interchange and drift wave are unstable and a high field side SOL region where only the drift wave is unstable. The SOL profiles appear exponential in the region close to the source but depart from a simple exponential fall-off in the far SOL. The low field side SOL width is found to be larger in the interchange alone case, slightly smaller when both instabilities are present and finally narrower when only the drift waves. For the high field side SOL, without interchange, the drift wave SOL width is observed to be identical to that on the low field side and larger than that when both instabilities at play. The Sherwood dimensionless parameter, ratio of convective particle flux divided by the diffusive particle flux, is used to compare the efficiency of turbulent transport. The profiles of the Sherwood parameter for time and flux surface averaged transport indicate that turbulent transport is dominant close to the separatrix but is less effective towards the far SOL. The Sherwood parameter evolution, determined with the flux-surface averaged transport, indicates that outward avalanche transport with corrugations governs the case with interchange only. When combining the two instabilities, outward avalanche transport is less pronounced and inward avalanche transport is observed, reducing the overall turbulent transport efficiency. The avalanche transport with drift waves only compared to interchange only is found to be inhibited.
Philippe Ghendrih, Guilhem Dif-Pradalier, Olivier Panico, Yanick Sarazin, Hugo Bufferand, et al.. Role of avalanche transport in competing drift wave and interchange turbulence. Journal of Physics: Conference Series, 2022, 2397, pp.012018. ⟨10.1088/1742-6596/2397/1/012018⟩. ⟨hal-04498387⟩
Georis Billo, Michel Belliard, Pierre Sagaut. Comparison of several interpolation methods to reconstruct field data in the vicinity of a finite element immersed boundary. Computers & Mathematics with Applications, 2022, 123, pp.123-135. ⟨10.1016/j.camwa.2022.08.002⟩. ⟨hal-04064030⟩ Plus de détails...
Thermal-hydraulics safety requirements for the second and third generation of nuclear reactors led to the development of innovative passive safety systems. In particular, new devices must be developed involving numerical simulations for turbulent two-phase flows around complex geometries. To reduce the time-consuming mesh generation phase when testing various geometries, we use a fictitious domain approach. More specifically, we choose the Penalized Direct Forcing Method to take into account inflow obstacles. Following a recent work, involving the resolution of the one-phase incompressible Navier-Stokes equations using a projection scheme and the Finite Element Method, this paper focuses on different techniques to recover data from the discrete immersed boundary and different ways to achieve order 2 in space via linear interpolation. Indeed, we investigate two data reconstruction approaches (one based on various weighted averaging, the other based on optimization) and compare their results for cylindrical and NACA0012 airfoil shapes: they provide similar accuracy but the weighting is much faster in terms of execution time. We also investigate three different interpolation types: unidirectional, multi-directional and a new hybrid between the two. The Taylor-Couette flow and the flow around a circular cylinder are used to carry out mesh convergence studies. Globally, order 2 in space is numerically assessed in both 2 and ∞ norms for all the interpolation types, which is consistent with theoretical expectations-even if the space convergence order is a bit higher for the multi-directional approach. For the flow around a circular cylinder, the values of aerodynamic coefficients and Strouhal number are in good agreement with the literature, especially when using directional interpolation. Finally, an industrial case, representative of passive safety systems, is presented to assess the robustness and capability of the method. The simulations tend to show that, here again, the directional interpolation offers the best behavior when dealing with complex geometries and relatively coarse meshes.
Georis Billo, Michel Belliard, Pierre Sagaut. Comparison of several interpolation methods to reconstruct field data in the vicinity of a finite element immersed boundary. Computers & Mathematics with Applications, 2022, 123, pp.123-135. ⟨10.1016/j.camwa.2022.08.002⟩. ⟨hal-04064030⟩
Journal: Computers & Mathematics with Applications
Ivan Kudashev, Anna Medvedeva, Manuel Scotto D’abusco, Nicolas Fedorszak, Stefano Di Genova, et al.. Development of a set of synthetic diagnostics for the confrontation between 2D transport simulations and WEST tokamak experimental data. Applied Sciences, 2022, 12 (19), pp.9807. ⟨10.3390/app12199807⟩. ⟨hal-03982630⟩ Plus de détails...
Transport codes are frequently used for describing fusion plasmas with the aim to prepare tokamak operations. Considering novel codes, such as SolEdge3X-HDG, synthetic diagnostics are a common technique used to validate new models and confront them with experimental data. The purpose of this study is to develop a set of synthetic diagnostics, starting from bolometer and visible cameras for the WEST tokamak, in order to compare the code results with the experimental data. This research is done in the framework of Raysect and Cherab Python libraries. This allows us to process various synthetic diagnostics in the same fashion in terms of 3D ray tracing with volume emitters developed specifically for fusion plasmas. We were able to implement the WEST tokamak model and the design of bolometer and visible cameras. Synthetic signals, based on full-discharge WEST plasma simulation, were used for to compare the SolEdge3X-HDG output plasma with experimental data. The study also considers the optical properties of the plasma-facing components (PFCs) and their influence on the performance of diagnostics. The paper shows a unified approach to synthetic diagnostic design, which will be further extended to cover the remaining diagnostics on the WEST tokamak.
Ivan Kudashev, Anna Medvedeva, Manuel Scotto D’abusco, Nicolas Fedorszak, Stefano Di Genova, et al.. Development of a set of synthetic diagnostics for the confrontation between 2D transport simulations and WEST tokamak experimental data. Applied Sciences, 2022, 12 (19), pp.9807. ⟨10.3390/app12199807⟩. ⟨hal-03982630⟩
Shang-Gui Cai, Sajad Mozaffari, Jérôme Jacob, Pierre Sagaut. Application of immersed boundary based turbulence wall modeling to the Ahmed body aerodynamics. Physics of Fluids, 2022, 34 (9), pp.095106. ⟨10.1063/5.0098232⟩. ⟨hal-04065468⟩ Plus de détails...
This paper applies a recently developed immersed boundary-turbulence wall modeling approach to turbulent flows over a generic car geometry, known as the Ahmed body, under massive flow separation within a lattice Boltzmann solver. Although the immersed boundary method combined with hierarchical Cartesian grid offers high flexibility in automatic grid generation around complex geometries, the near-wall solution is significantly deteriorated compared to the body-fitted simulation, especially when coupled to wall models for turbulent flows at high Reynolds number. Enhanced wall treatments have been proposed in the literature and validated for attached flow configurations. In this work, the Ahmed body with a slant surface of angle 35 is considered where the flow separates massively over the slant surface and the vertical base. The large eddy simulation is performed with a Reynolds stress constraint near-wall. The eddy viscosity is computed dynamically by taking into account the actually resolved Reynolds stresses. It approaches the mixing length eddy viscosity in attached boundary layers and returns to the subgrid eddy viscosity in detached boundary layers. An explicit equilibrium wall model has also been proposed to accelerate the calculation. Comparison with the no-slip boundary condition on the separated surfaces shows that the near-wall treatments with the equilibrium wall model operate reasonably well on both attached and detached boundary layers.
Shang-Gui Cai, Sajad Mozaffari, Jérôme Jacob, Pierre Sagaut. Application of immersed boundary based turbulence wall modeling to the Ahmed body aerodynamics. Physics of Fluids, 2022, 34 (9), pp.095106. ⟨10.1063/5.0098232⟩. ⟨hal-04065468⟩
Fatemeh Safaei, Jafar Ahmadi, Mitra Fouladirad. Optimal N-policy for the maintenance of k-out-of-n systems with dynamic minor repairs considering second-hand component income. International Journal of Production Research, 2022, pp.1-18. ⟨10.1080/00207543.2022.2120107⟩. ⟨hal-04064518⟩ Plus de détails...
This paper proposes a maintenance policy for a k-out-of-n: F system operating to fulfil several jobs without interruptions. The system is replaced when the Nth job completes or at the kth failure, whichever occurs first. So, there are some non-failed components when a replacement is done. These components can be sold as secondhand products to continue working in other systems for a while. The price of them is also considered in the proposed maintenance policy. An optimal maintenance policy is studied to minimise long-run average cost under the constraint of relative mean operating time. Comprehensive numerical studies are done to assess the effect of the model parameters on the optimal solutions. Also, to demonstrate the applicability of the proposed plan, a data set related to wind turbine generator failures is considered as a case study.
Fatemeh Safaei, Jafar Ahmadi, Mitra Fouladirad. Optimal N-policy for the maintenance of k-out-of-n systems with dynamic minor repairs considering second-hand component income. International Journal of Production Research, 2022, pp.1-18. ⟨10.1080/00207543.2022.2120107⟩. ⟨hal-04064518⟩
Journal: International Journal of Production Research
D Auroux, P Ghendrih, L Lamerand, F Rapetti, E Serre. Asymptotic behaviour, non-local dynamics and data assimilation tailoring of the reduced κ − ε model to address turbulent transport of fusion plasmas. Physics of Plasmas, 2022. ⟨hal-03811621⟩ Plus de détails...
The high-dimensional and multiscale nature of fusion plasma flows require the development of reduced models to be implemented in numerical codes capable of capturing the main features of turbulent transport in a sufficiently short time to be useful during tokamak operation. This paper goes further in the analysis of the dynamics of the κ − ε model based on the turbulent kinetic energy κ and its dissipation rate ε [Baschetti et al., Nuc. Fus 61, 106020 (2021)] to improve the predictability of the transverse turbulent transport in simulation codes. Present 1D results show further capabilities with respect to current models (based on constant effective perpendicular diffusion) and on the standard quasi-linear approach. The nonlinear dependence of D in κ and ε estimated from two additional transport equations allow to introduce some non-locality in the transport model. This is illustrated by the existence of parameter ranges with turbulence spreading. The paper also addresses another issue related to the uncertainties on the inherent free parameters of such reduced model. The study proposes a new approach in the fusion community based on a variational data assimilation involving the minimisation of a cost function defined as the distance between the reference data and the calculated values. The results are good, and show the ability of the data assimilation to reduce uncertainties on the free parameters, which remains a critical point to ensure the total reliability of such an approach.
D Auroux, P Ghendrih, L Lamerand, F Rapetti, E Serre. Asymptotic behaviour, non-local dynamics and data assimilation tailoring of the reduced κ − ε model to address turbulent transport of fusion plasmas. Physics of Plasmas, 2022. ⟨hal-03811621⟩
Tuan Dung Nguyen, Christophe Besse, François Rogier. High-order Scharfetter-Gummel-based schemes and applications to gas discharge modeling. Journal of Computational Physics, 2022, 461, pp.111196. ⟨10.1016/j.jcp.2022.111196⟩. ⟨hal-03352814⟩ Plus de détails...
A generalized Scharfetter-Gummel method is proposed to construct the numerical flux for one-dimensional drift-diffusion equations. Instead of taking a constant approximation of the flux as Scharfetter and Gummel did in [1], we consider a p-degree polynomial with p ≥ 1. The high order moments of the approximating flux function serve as intermediaries to bring numerical correction to the Scharfetter-Gummel flux, that the other end turns out to be the solution derivatives. Therefore, local solution reconstructions are required. The resulting schemes are high order and discretize at the same time the convective and diffusive fluxes without having to employ separately different methods to do so. The new schemes with p = 1 and p = 2 are employed to simulate atmospheric pressure discharge where they are applied to the continuity equations for electrons and ions, and solved simultaneously with Poisson's equation. Numerical results indicate that our method are robust and highly accurate.
Tuan Dung Nguyen, Christophe Besse, François Rogier. High-order Scharfetter-Gummel-based schemes and applications to gas discharge modeling. Journal of Computational Physics, 2022, 461, pp.111196. ⟨10.1016/j.jcp.2022.111196⟩. ⟨hal-03352814⟩
B. Bugeat, J.-Ch. Robinet, J.-C. Chassaing, P. Sagaut. Low-frequency resolvent analysis of the laminar oblique shock wave/boundary layer interaction. Journal of Fluid Mechanics, 2022, 942, pp.A43. ⟨10.1017/jfm.2022.390⟩. ⟨hal-04064039⟩ Plus de détails...
Resolvent analysis is used to study the low-frequency behaviour of the laminar oblique shock wave/boundary layer interaction (SWBLI). It is shown that the computed optimal gain, which can be seen as a transfer function of the system, follows a first-order low-pass filter equation, recovering the results of Touber & Sandham (J. Fluid Mech., vol. 671, 2011, pp. 417-465). This behaviour is understood as proceeding from the excitation of a single stable, steady global mode whose damping rate sets the time scale of the filter. Different Mach and Reynolds numbers are studied, covering different recirculation lengths L. This damping rate is found to scale as 1/L, leading to a constant Strouhal number St L as observed in the literature. It is associated with a breathing motion of the recirculation bubble. This analysis furthermore supports the idea that the low-frequency dynamics of the SWBLI is a forced dynamics, in which background perturbations continuously excite the flow. The investigation is then carried out for three-dimensional perturbations for which two regimes are identified. At low wavenumbers of the order of L, a modal mechanism similar to that of two-dimensional perturbations is found and exhibits larger values of the optimal gain. At larger wavenumbers, of the order of the boundary layer thickness, the growth of streaks, which results from a non-modal mechanism, is detected. No interaction with the recirculation region is observed. Based on these results, the potential prevalence of three-dimensional effects in the low-frequency dynamics of the SWBLI is discussed.
B. Bugeat, J.-Ch. Robinet, J.-C. Chassaing, P. Sagaut. Low-frequency resolvent analysis of the laminar oblique shock wave/boundary layer interaction. Journal of Fluid Mechanics, 2022, 942, pp.A43. ⟨10.1017/jfm.2022.390⟩. ⟨hal-04064039⟩
Mitra Fouladirad, Massimiliano Giorgio, Gianpaolo Pulcini. A transformed gamma process for bounded degradation phenomena. Quality and Reliability Engineering International, 2022, 39 (2), pp.546 - 564. ⟨10.1002/qre.3167⟩. ⟨hal-04063942⟩ Plus de détails...
Most of the stochastic models adopted to describe the evolution over time of degradation phenomena of technological units assume that their degradation level can increase indeterminately. However, these degradation phenomena are typically subjected to obvious bounds, if only because technological units have finite size. In fact, very often, this inconsistency does not significantly affect the effectiveness of unbounded degradation models, since degrading units are usually assumed to fail when their degradation level exceeds a failure threshold that is much smaller than the obvious bounds. Nevertheless, in some cases, due to the very nature of the underlying degradation mechanism, less obvious bounds could exist, which are not necessarily far from the failure thresholds. The question that arises is whether the use of a bounded degradation model, in this latter type of experimental situations, could be beneficial. For this purpose, since a bounded degradation process should necessarily have dependent increments, in this paper we investigate the potential of a new bounded transformed gamma (TG) process to adequately describe bounded degradation phenomena and predict their future evolution. Differently from other existing gamma process based bounded degradation models, here the upper bound is treated as an unknown parameter that has to be estimated from the available degradation data. A numerical example is presented where the parameters of the proposed model are estimated from simulated data. Then the model is applied to a set of wear measures of cylinder liners that equip a diesel engine for marine propulsion, which have also stimulated this study. Model parameters are estimated by using the maximum likelihood (ML) method. The fitting ability of the proposed new bounded process is compared to that of an unbounded gamma process, which was previously adopted to analyze the same liner wear data. Obtained results are critically discussed in the paper.
Mitra Fouladirad, Massimiliano Giorgio, Gianpaolo Pulcini. A transformed gamma process for bounded degradation phenomena. Quality and Reliability Engineering International, 2022, 39 (2), pp.546 - 564. ⟨10.1002/qre.3167⟩. ⟨hal-04063942⟩
Journal: Quality and Reliability Engineering International
Marwa Belhaj Salem, Mitra Fouladirad, Estelle Deloux. Variance Gamma process as degradation model for prognosis and imperfect maintenance of centrifugal pumps. Reliability Engineering and System Safety, 2022, 223, pp.108417. ⟨10.1016/j.ress.2022.108417⟩. ⟨hal-04063946⟩ Plus de détails...
This paper considers the degradation modelling of a non-monotonous health indicator. A Variance Gamma process is proposed for the degradation modelling and its calibration in presence of data is discussed. The remaining useful lifetime estimation based on this model is considered and its sensitivity to parameters estimates is analysed. The model is applied to real data of leakage rate of centrifugal pump. Eventually, an imperfect maintenance policy is proposed and optimized.
Marwa Belhaj Salem, Mitra Fouladirad, Estelle Deloux. Variance Gamma process as degradation model for prognosis and imperfect maintenance of centrifugal pumps. Reliability Engineering and System Safety, 2022, 223, pp.108417. ⟨10.1016/j.ress.2022.108417⟩. ⟨hal-04063946⟩
Journal: Reliability Engineering and System Safety
Hoda Kamranfar, Kambiz Ahmadi, Mitra Fouladirad. On inference in a class of exponential distribution under imperfect maintenance. Communications in Statistics - Simulation and Computation, 2022, pp.1-27. ⟨10.1080/03610918.2022.2103567⟩. ⟨hal-04064552⟩ Plus de détails...
This paper deals with statistical inference for lifetime data in presence of imperfect maintenance. For the maintenance model, the Sheu and Griffith model is considered. The lifetime distribution belongs to exponential distribution class. The maximum likelihood estimation procedure of the model parameters is discussed, and confidence intervals are provided using the asymptotic likelihood theory and bootstrap approach. Based on conjugate and discrete priors, Bayesian estimators of the model parameters are developed under symmetric and asymmetric loss functions. The proposed methodologies are applied to simulated data and sensitivity analysis to different parameters and data characteristics is carried out. The effect of model misspecification is also assessed within this class of distributions through a Monte Carlo simulation study. Finally, two datasets are analyzed for demonstrative aims.
Hoda Kamranfar, Kambiz Ahmadi, Mitra Fouladirad. On inference in a class of exponential distribution under imperfect maintenance. Communications in Statistics - Simulation and Computation, 2022, pp.1-27. ⟨10.1080/03610918.2022.2103567⟩. ⟨hal-04064552⟩
Journal: Communications in Statistics - Simulation and Computation
J. Etminan, H. Kamranfar, M. Chahkandi, M. Fouladirad. Analysis of time-to-failure data for a repairable system subject to degradation. Journal of Computational and Applied Mathematics, 2022, 408, pp.114098. ⟨10.1016/j.cam.2022.114098⟩. ⟨hal-04063947⟩ Plus de détails...
In this paper, a gradually deteriorating system with imperfect repair is considered. The deterioration is modeled by a stationary stochastic process. The system fails once the deterioration level exceeds a given threshold L. At failure, an imperfect repair is performed and the deterioration level is reduced to a fixed value r, say. The system can be repaired n − 1 times and will be replaced after the nth failure. The article aims to estimate the parameters of the proposed deterioration process based on the observed failures. To this end we consider the Wiener and Gamma processes which are the most common used stochastic process models. In Wiener process, an explicit expression for the estimators is obtained. Birnbaum-Saunders approximation is extended to estimate the parameters in Gamma process. An optimal replacement policy is also discussed. Finally, a Monte-Carlo simulation is conducted to investigate the performance of estimators.
J. Etminan, H. Kamranfar, M. Chahkandi, M. Fouladirad. Analysis of time-to-failure data for a repairable system subject to degradation. Journal of Computational and Applied Mathematics, 2022, 408, pp.114098. ⟨10.1016/j.cam.2022.114098⟩. ⟨hal-04063947⟩
Journal: Journal of Computational and Applied Mathematics
Olivier Mesdjian, Chenglei Wang, Simon Gsell, Umberto D’ortona, Julien Favier, et al.. Longitudinal to Transverse Metachronal Wave Transitions in an In Vitro Model of Ciliated Bronchial Epithelium. Physical Review Letters, 2022, 129 (3), pp.038101. ⟨10.1103/PhysRevLett.129.038101⟩. ⟨hal-03741505⟩ Plus de détails...
Olivier Mesdjian, Chenglei Wang, Simon Gsell, Umberto D’ortona, Julien Favier, et al.. Longitudinal to Transverse Metachronal Wave Transitions in an In Vitro Model of Ciliated Bronchial Epithelium. Physical Review Letters, 2022, 129 (3), pp.038101. ⟨10.1103/PhysRevLett.129.038101⟩. ⟨hal-03741505⟩
G. Piraccini, F. Schwander, E. Serre, G. Giorgiani, M. Scotto d'Abusco. Spatial adaptivity in SOLEDGE3X‐HDG for edge plasma simulations in versatile magnetic and reactor geometries. Contributions to Plasma Physics, 2022, 62 (5-6), ⟨10.1002/ctpp.202100185⟩. ⟨hal-04063959⟩ Plus de détails...
With the ultimate goal to predict plasmas heat and particle fluxes in ITER opera- tion, more efforts are required to deal with realistic magnetic configurations and tokamak geometries. In an attempt to achieve this goal, we propose an adaptive mesh refinement method added to a fluid solver based on a high-order hybrid discontinuous Galerkin (HDG) method. Based on unstructured meshes, this magnetic equilibrium free numerical scheme has shown promising and encour- aging features to solve 2D/3D transport reduced Braginski fluid equations. To improve its numerical efficiency, a mesh refinement based on h-adpativity is investigated. We describe here an adaptive refinement strategy on a reduced edge particle transport model based on electron density and parallel momentum. This strategy is illustrated in realistic tokamak wall geometry. Computations performed show potential gains in the required number of degrees of freedom against benchmark computations with uniform meshes, along with the poten- tial to give an automated, goal-oriented, mesh generation technique for edge transport simulations in 2D
G. Piraccini, F. Schwander, E. Serre, G. Giorgiani, M. Scotto d'Abusco. Spatial adaptivity in SOLEDGE3X‐HDG for edge plasma simulations in versatile magnetic and reactor geometries. Contributions to Plasma Physics, 2022, 62 (5-6), ⟨10.1002/ctpp.202100185⟩. ⟨hal-04063959⟩
Hasan Misaii, Firoozeh Haghighi, Mitra Fouladirad. Optimal shock‐based maintenance policy for a system in a dynamic environment. Applied Stochastic Models in Business and Industry, 2022, 38 (6), pp.918 - 934. ⟨10.1002/asmb.2686⟩. ⟨hal-04063939⟩ Plus de détails...
In this article, we consider a single-unit system that operates in a dynamic environment and is subject to shocks. Shocks only affect the system (nonlethal shock) and do not cause sudden failure, and arrive according to a counting process. Both homogeneous and nonhomogeneous Poisson processes are considered for shocks arrival modeling. In order to model the dynamic environment and consider shock effects, a multiplicative failure rate model is proposed. Both corrective maintenance and shock-based preventive maintenance are considered, and two policies are proposed. In the first proposed policy, the system is replaced by a new one upon a failure or based on the predetermined number of shocks, whichever comes first, while the second proposed policy extends the first one by considering an imperfect preventive repair at each inspection time. The inspection times are periodic and the interinspection interval is considered as a decision variable. The proposed policies are optimized according to long-run cost rate criteria. Numerical examples illustrate the applicability and efficiency of the proposed policies.
Hasan Misaii, Firoozeh Haghighi, Mitra Fouladirad. Optimal shock‐based maintenance policy for a system in a dynamic environment. Applied Stochastic Models in Business and Industry, 2022, 38 (6), pp.918 - 934. ⟨10.1002/asmb.2686⟩. ⟨hal-04063939⟩
Journal: Applied Stochastic Models in Business and Industry
M. Nguyen, J. Boussuge, P. Sagaut, J. Larroya-Huguet. Large eddy simulation of a thermal impinging jet using the lattice Boltzmann method. Physics of Fluids, 2022, 34 (5), pp.055115. ⟨10.1063/5.0088410⟩. ⟨hal-03669901⟩ Plus de détails...
A compressible Hybrid Lattice Boltzmann Method solver is used to perform a wall-resolved Large eddy simulation of an isothermal axisymmetric jet issuing from a pipe and impinging on a heated flat plate at a Reynolds number of 23 000, a Mach number of 0.1, and an impingement distance of two jet diameters. The jet flow field statistics, Nusselt number profile (including the secondary peak), and shear stress profile were well reproduced. The azimuthal coherence of the primary vortical structures was relatively low, leading to no discernible temporal periodicity of the azimuthally averaged Nusselt number at the location of the secondary peak. While local unsteady near-wall flow separation was observed in the wall jet, this flow separation did not exhibit azimuthal coherence and was not found to be the only cause of the thermal spots blue, which lead to the secondary peak in the Nusselt number, as stream-wise oriented structures also played a significant role in increasing the local heat transfer.
M. Nguyen, J. Boussuge, P. Sagaut, J. Larroya-Huguet. Large eddy simulation of a thermal impinging jet using the lattice Boltzmann method. Physics of Fluids, 2022, 34 (5), pp.055115. ⟨10.1063/5.0088410⟩. ⟨hal-03669901⟩
Antoine Galko, Simon Gsell, Umberto d'Ortona, Laurent Morin, Julien Favier. Pulsated Herschel-Bulkley flows in two-dimensional channels: A model for mucus clearance devices. Physical Review Fluids, 2022, 7 (5), pp.053301. ⟨10.1103/PhysRevFluids.7.053301⟩. ⟨hal-03863329⟩ Plus de détails...
Guanxiong Wang, Song Zhao, Pierre Boivin, Eric Serre, Pierre Sagaut. A new hybrid lattice-Boltzmann method for thermal flow simulations in low-Mach number approximation. Physics of Fluids, 2022, 34 (4), pp.046114. ⟨10.1063/5.0091517⟩. ⟨hal-03796386⟩ Plus de détails...
A new low-Mach algorithm for the thermal lattice Boltzmann method (LBM) is proposed aiming at reducing the computational cost of thermal flow simulations in the low Mach number limit. The well-known low Mach number approximation is adopted to accelerate the simulations by enlarging the time step through re-scaling the psuedoacoustic speed to the same order of the fluid motion velocity. This specific process is inspired by the similarity between the artificial compressibility method and the isothermal LBM and is further extended to its thermal counterpart. It must be emphasized that such low-Mach acceleration strategy is in a general form, thus can be easily applied to other compressible LB methods. The present method overcomes the drawback of the classical pressure gradient scaling method due to the pressure gradient changing. The new algorithm is validated by various well-documented academic test cases in laminar [one dimensional gravity column, 2D (two dimensional) rising thermal bubble, and 2D differentially heated square cavity] and turbulent [3D (three dimensional) Taylor–Green vortex and 3D heated cylinder] regimes. All the results show excellent agreement with the reference data and high computational efficiency.
Guanxiong Wang, Song Zhao, Pierre Boivin, Eric Serre, Pierre Sagaut. A new hybrid lattice-Boltzmann method for thermal flow simulations in low-Mach number approximation. Physics of Fluids, 2022, 34 (4), pp.046114. ⟨10.1063/5.0091517⟩. ⟨hal-03796386⟩
H Bufferand, J Balbin, S Baschetti, J Bucalossi, G Ciraolo, et al.. Implementation of multi-component Zhdanov closure in SOLEDGE3X. Plasma Physics and Controlled Fusion, 2022, 64 (5), pp.055001. ⟨10.1088/1361-6587/ac4fac⟩. ⟨hal-03981783⟩ Plus de détails...
The multi-component fluid closure derived by Zhdanov (2002 Transport Processes in Multicomponent Plasma (London: Taylor and Francis)) is implemented in the fluid code SOLEDGE3X-EIRENE to deal with arbitrary edge plasma composition. The closure assumes no distinction between species such as light versus heavy species separation. The work of Zhdanov is rewritten in a matricial form in order to clearly link friction forces and heat fluxes to the different species velocities and temperature gradients.
H Bufferand, J Balbin, S Baschetti, J Bucalossi, G Ciraolo, et al.. Implementation of multi-component Zhdanov closure in SOLEDGE3X. Plasma Physics and Controlled Fusion, 2022, 64 (5), pp.055001. ⟨10.1088/1361-6587/ac4fac⟩. ⟨hal-03981783⟩
Rouae Ben Dhia, Nils Tilton, Denis Martinand. Impact of osmotic pressure on the stability of Taylor vortices. Journal of Fluid Mechanics, 2022, 933, pp.A51. ⟨10.1017/jfm.2021.1101⟩. ⟨hal-03533753⟩ Plus de détails...
We use linear stability analysis and direct numerical simulations to investigate the coupling between centrifugal instabilities, solute transport and osmotic pressure in a Taylor–Couette configuration that models rotating dynamic filtration devices. The geometry consists of a Taylor–Couette cell with a superimposed radial throughflow of solvent across two semi-permeable cylinders. Both cylinders totally reject the solute, inducing the build-up of a concentration boundary layer. The solute retroacts on the velocity field via the osmotic pressure associated with the concentration differences across the semi-permeable cylinders. Our results show that the presence of osmotic pressure strongly alters the dynamics of the centrifugal instabilities and substantially reduces the critical conditions above which Taylor vortices are observed. It is also found that this enhancement of the hydrodynamic instabilities eventually plateaus as the osmotic pressure is further increased. We propose a mechanism to explain how osmosis and instabilities cooperate and develop an analytical criterion to bound the parameter range for which osmosis fosters the hydrodynamic instabilities.
Rouae Ben Dhia, Nils Tilton, Denis Martinand. Impact of osmotic pressure on the stability of Taylor vortices. Journal of Fluid Mechanics, 2022, 933, pp.A51. ⟨10.1017/jfm.2021.1101⟩. ⟨hal-03533753⟩
M Raghunathan, Y Marandet, H Bufferand, G Ciraolo, Ph Ghendrih, et al.. Multi-temperature Generalized Zhdanov Closure for Scrape-Off Layer/Edge Applications. Plasma Physics and Controlled Fusion, 2022. ⟨hal-03810579⟩ Plus de détails...
The derivation of the multi-temperature generalized Zhdanov closure is provided starting from the most general form of the left hand side of the moment averaged kinetic equation with the Sonine-Hermite polynomial ansatz for an arbitrary number of moments. The process of arriving at the reduced higher-order moment equations, with its assumptions and approximations, is explicitly outlined. The generalized multi-species, multi-temperature coefficients from the authors' previous article are used to compute values of higher order moments such as heat flux in terms of the lower order moments. Transport coefficients and the friction and thermal forces for magnetic confinement fusion relevant cases with the generalized coefficients are compared to the scheme with the single-temperature coefficients previously provided by Zhdanov et al. It is found that the 21N-moment multi-temperature coefficients are adequate for most cases relevant to fusion. Furthermore, the 21N-moment scheme is also tested against the trace approximation to determine the range of validity of the trace approximation with respect to fusion relevant plasmas. Possible refinements to the closure scheme are illustrated as well, in order to account for quantities which might be significant in certain schemes such as the drift approximation.
M Raghunathan, Y Marandet, H Bufferand, G Ciraolo, Ph Ghendrih, et al.. Multi-temperature Generalized Zhdanov Closure for Scrape-Off Layer/Edge Applications. Plasma Physics and Controlled Fusion, 2022. ⟨hal-03810579⟩
X Litaudon, F Jenko, D Borba, D Borodin, B J Braams, et al.. EUROfusion-theory and advanced simulation coordination (E-TASC): programme and the role of high performance computing. Plasma Physics and Controlled Fusion, 2022, 64 (3), pp.034005. ⟨10.1088/1361-6587/ac44e4⟩. ⟨hal-03562886⟩ Plus de détails...
Abstract This paper is a written summary of an overview oral presentation given at the 1st Spanish Fusion High Performance Computer (HPC) Workshop that took place on the 27 November 2020 as an online event. Given that over the next few years ITER 24 24 ITER (‘The Way’ in Latin) is the world’s largest tokamak under construction in the south of France: a magnetic fusion device that has been designed to prove the feasibility of fusion as a large-scale and carbon-free source of energy ( https://www.iter.org/ ). will move to its operation phase and the European-DEMO design will be significantly advanced, the EUROfusion consortium has initiated a coordination effort in theory and advanced simulation to address some of the challenges of the fusion research in Horizon EUROPE (2021–2027), i.e. the next EU Framework Programme for Research and Technological Development. This initiative has been called E-TASC, which stands for EUROfusion-Theory and Advanced Simulation Coordination. The general and guiding principles of E-TASC are summarized in this paper. In addition, an overview of the scientific results obtained in the pilot phase (2019–2020) of E-TASC are provided while highlighting the importance of the required progress in computational methods and HPC techniques. In the initial phase, five pilot theory and simulation tasks were initiated: towards a validated predictive capability of the low to high transition and pedestal physics; runaway electrons in tokamak disruptions in the presence of massive material injection; fast code for the calculation of neoclassical toroidal viscosity in stellarators and tokamaks; development of a neutral gas kinetics modular code; European edge and boundary code for reactor-relevant devices. In this paper, we report on recent progress made by each of these projects.
X Litaudon, F Jenko, D Borba, D Borodin, B J Braams, et al.. EUROfusion-theory and advanced simulation coordination (E-TASC): programme and the role of high performance computing. Plasma Physics and Controlled Fusion, 2022, 64 (3), pp.034005. ⟨10.1088/1361-6587/ac44e4⟩. ⟨hal-03562886⟩
Gauthier Wissocq, Pierre Sagaut. Hydrodynamic limits and numerical errors of isothermal lattice Boltzmann schemes. Journal of Computational Physics, 2022, 450, pp.110858. ⟨10.1016/j.jcp.2021.110858⟩. ⟨hal-04064045⟩ Plus de détails...
With the aim of better understanding the numerical properties of the lattice Boltzmann method (LBM), a general methodology is proposed to derive its hydrodynamic limits in the discrete setting. It relies on a Taylor expansion in the limit of low Knudsen numbers. With a single asymptotic analysis, two kinds of deviations with the Navier-Stokes (NS) equations are explicitly evidenced: consistency errors, inherited from the kinetic description of the LBM, and numerical errors attributed to its space and time discretization. The methodology is applied to the Bhatnagar-Gross-Krook (BGK), the regularized and the multiple relaxation time (MRT) collision models in the isothermal framework. Deviation terms are systematically confronted to linear analyses in order to validate their expressions, interpret them and provide explanations for their numerical properties. The low dissipation of the BGK model is then related to a particular pattern of its error terms in the Taylor expansion. Similarly, dissipation properties of the regularized and MRT models are explained by a phenomenon referred to as hyperviscous degeneracy. The latter consists in an unexpected resurgence of high-order Knudsen effects induced by a large numerical prefactor. It is at the origin of over-dissipation and severe instabilities in the low-viscosity regime.
Gauthier Wissocq, Pierre Sagaut. Hydrodynamic limits and numerical errors of isothermal lattice Boltzmann schemes. Journal of Computational Physics, 2022, 450, pp.110858. ⟨10.1016/j.jcp.2021.110858⟩. ⟨hal-04064045⟩
D. Galassi, C. Theiler, T. Body, F. Manke, P. Micheletti, et al.. Validation of edge turbulence codes in a magnetic X-point scenario in TORPEX. Physics of Plasmas, 2022, 29 (1), pp.012501. ⟨10.1063/5.0064522⟩. ⟨hal-03566373⟩ Plus de détails...
D. Galassi, C. Theiler, T. Body, F. Manke, P. Micheletti, et al.. Validation of edge turbulence codes in a magnetic X-point scenario in TORPEX. Physics of Plasmas, 2022, 29 (1), pp.012501. ⟨10.1063/5.0064522⟩. ⟨hal-03566373⟩
M Scotto d'Abusco, G Giorgiani, J F Artaud, H Bufferand, G Ciraolo, et al.. Core-edge 2D fluid modeling of full tokamak discharge with varying magnetic equilibrium: from WEST start-up to ramp-down. Nuclear Fusion, 2022, ⟨10.1088/1741-4326/ac47ad⟩. ⟨hal-03509800⟩ Plus de détails...
In the present work we investigate for the first time the 2D fluid transport of the plasma in WEST during an entire discharge from the start-up to the ramp-down (shot #54487). The evolution of density profile, electron and ion temperatures together with the experimental magnetic equilibrium, total current and gas-puff rate is investigated. Comparisons with the interferometry diagnostic show a remarkable overall qualitative agreement during the discharge that can be quantitative at some locations in the plasma core. If at the onset of the X-points during the ramp-up the electron heat flux is dominant at the target, present results show that the ion heat flux becomes dominant during the stationary phase of the discharge. Using a simple model for erosion, present results assess the tungsten sputtering due to deuterium ions during the start-up and ramp-up phases of the discharge and confirms the need to consider full discharge simulation to accurately treat the W source of contamination. This work also demonstrates the interest of developing magnetic equilibrium free solver including efficient time integration to step toward predictive capabilities in the future for fusion operation.
M Scotto d'Abusco, G Giorgiani, J F Artaud, H Bufferand, G Ciraolo, et al.. Core-edge 2D fluid modeling of full tokamak discharge with varying magnetic equilibrium: from WEST start-up to ramp-down. Nuclear Fusion, 2022, ⟨10.1088/1741-4326/ac47ad⟩. ⟨hal-03509800⟩
Guanxiong Wang, Song Zhao, Pierre Boivin, Eric Serre, Pierre Sagaut. A new hybrid Lattice-Boltzmann method for thermal flow simulations in low-Mach number approximation. Physics of Fluids, 2022, Physics of fluids, 34 (046114). ⟨hal-03636905⟩ Plus de détails...
Guanxiong Wang, Song Zhao, Pierre Boivin, Eric Serre, Pierre Sagaut. A new hybrid Lattice-Boltzmann method for thermal flow simulations in low-Mach number approximation. Physics of Fluids, 2022, Physics of fluids, 34 (046114). ⟨hal-03636905⟩
D.S. Oliveira, T. Body, D. Galassi, C. Theiler, E. Laribi, et al.. Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case. Nuclear Fusion, 2022, 62 (9), pp.096001. ⟨10.1088/1741-4326/ac4cde⟩. ⟨hal-03740877⟩ Plus de détails...
Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21 , contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.
D.S. Oliveira, T. Body, D. Galassi, C. Theiler, E. Laribi, et al.. Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case. Nuclear Fusion, 2022, 62 (9), pp.096001. ⟨10.1088/1741-4326/ac4cde⟩. ⟨hal-03740877⟩
Gauthier Wissocq, Thomas Coratger, Gabriel Farag, Song Zhao, Pierre Boivin, et al.. Restoring the conservativity of characteristic-based segregated models: application to the hybrid lattice Boltzmann method. Physics of Fluids, 2022, 34 (4), pp.046102. ⟨10.1063/5.0083377⟩. ⟨hal-03627520⟩ Plus de détails...
A general methodology is introduced to build conservative numerical models for fluid simulations based on segregated schemes, where mass, momentum and energy equations are solved by different methods. It is here especially designed for developing new numerical discretizations of the total energy equation, adapted to a thermal coupling with the lattice Boltzmann method (LBM). The proposed methodology is based on a linear equivalence with standard discretizations of the entropy equation, which, as a characteristic variable of the Euler system, allows efficiently decoupling the energy equation with the LBM. To this extent, any LBM scheme is equivalently written under a finite-volume formulation involving fluxes, which are further included in the total energy equation as numerical corrections. The viscous heat production is implicitly considered thanks to the knowledge of the LBM momentum flux. Three models are subsequently derived: a first-order upwind, a Lax-Wendroff and a third-order Godunov-type schemes. They are assessed on standard academic test cases: a Couette flow, entropy spot and vortex convections, a Sod shock tube, several two-dimensional Riemann problems and a shock-vortex interaction. Three key features are then exhibited: 1) the models are conservative by construction, recovering correct jump relations across shock waves, 2) the stability and accuracy of entropy modes can be explicitly controlled, 3) the low dissipation of the LBM for isentropic phenomena is preserved.
Gauthier Wissocq, Thomas Coratger, Gabriel Farag, Song Zhao, Pierre Boivin, et al.. Restoring the conservativity of characteristic-based segregated models: application to the hybrid lattice Boltzmann method. Physics of Fluids, 2022, 34 (4), pp.046102. ⟨10.1063/5.0083377⟩. ⟨hal-03627520⟩
Lincheng Xu, Eric Serre, Pierre Sagaut. A theoretical analysis of mass leakage at boundaries within the lattice Boltzmann method. Physics of Fluids, 2022, Physics of fluids, 34 (065113). ⟨hal-03683744⟩ Plus de détails...
Mass leakage at boundaries can be a critical issue for reliability of the lattice Boltzmann (LB) method based on Cartesian grids. Despite numerous work based on the LB method, the intrinsic macroscopic mechanisms causing mass leakage are still not fully charac- terised, but are essential to improve the mass conservation of LB simulations. In this paper, an original theoretical investigation of mass leakage at boundaries is proposed within the general LB framework. It is demonstrated that the mass leakage originates from the in- trinsic deficiency of the wall-cut LB links at boundary nodes in recovering macroscopic momenta. From a mesoscopic-level definition, i.e. the net loss of distribution functions during the streaming process, the local mass leakage at individual boundary nodes and its averaged value along smooth boundaries are mathematically expressed using macroscopic variables. The local mass leakage is shown to be dominated by terms proportional to the tangential momentum component. In contrast, the averaged mass leakage is shown to be contributed from various terms including the boundary curvature, the tangential momen- tum, and the gradients of density, momentum and momentum flux. Meanwhile, amplitude of the averaged mass leakage is theoretically estimated to be proportional to the local grid spacing, based on which a first-order accurate correction scheme is proposed. In addition, both the local and averaged mass leakage are demonstrated to be significantly dependent on boundary orientation with respect to the grid. The proposed theoretical analysis is assessed by performing numerical experiments. Two-dimensional weakly compressible flows through straight and curved moving channels are considered to estimate each term appearing in the theoretical analysis. The numerical results are in very good agreement with the proposed analysis, and the proposed mass correction scheme based on the av- eraged mass leakage effectively cures the mass leakage problems in the considered test cases.
Lincheng Xu, Eric Serre, Pierre Sagaut. A theoretical analysis of mass leakage at boundaries within the lattice Boltzmann method. Physics of Fluids, 2022, Physics of fluids, 34 (065113). ⟨hal-03683744⟩
Giacomo Piraccini, Frédéric Schwander, Eric Serre, Giorgio Giorgiani, Manuel Scotto D’abusco. Spatial adaptivity in SOLEDGE3X‐HDG for edge plasma simulations in versatile magnetic and reactor geometries. Contributions to Plasma Physics, 2022, 18th International Workshop on Plasma Edge Theory in Fusion Devices September 13‐15, 2021, organized by the EPFL Swiss Plasma Center, 62 (5-6), ⟨10.1002/ctpp.202100185⟩. ⟨hal-04489377⟩ Plus de détails...
With the ultimate goal to predict plasmas heat and particle fluxes in ITER operation, more efforts are required to deal with realistic magnetic configurations and tokamak geometries. In an attempt to achieve this goal, we propose an adaptive mesh refinement method added to a fluid solver based on a high‐order hybrid discontinuous Galerkin (HDG) method. Based on unstructured meshes, this magnetic equilibrium free numerical scheme has shown promising and encouraging features to solve 2D/3D transport reduced Braginski fluid equations. To improve its numerical efficiency, a mesh refinement based on h‐adpativity is investigated. We describe here an adaptive refinement strategy on a reduced edge particle transport model based on electron density and parallel momentum. This strategy is illustrated in realistic tokamak wall geometry. Computations performed show potential gains in the required number of degrees of freedom against benchmark computations with uniform meshes, along with the potential to give an automated, goal‐oriented, mesh generation technique for edge transport simulations in 2D.
Giacomo Piraccini, Frédéric Schwander, Eric Serre, Giorgio Giorgiani, Manuel Scotto D’abusco. Spatial adaptivity in SOLEDGE3X‐HDG for edge plasma simulations in versatile magnetic and reactor geometries. Contributions to Plasma Physics, 2022, 18th International Workshop on Plasma Edge Theory in Fusion Devices September 13‐15, 2021, organized by the EPFL Swiss Plasma Center, 62 (5-6), ⟨10.1002/ctpp.202100185⟩. ⟨hal-04489377⟩
Giacomo Piraccini, Marcello Capasso, Manuel Scotto d'Abusco, Giorgio Giorgiani, Frédéric Schwander, et al.. Recent upgrades in a 2D turbulent transport solver based on a hybrid discontinuous Galerkin method for the simulation of fusion plasma in tokamak. Fluids, 2022, ⟨10.3390/fluids7020063⟩. ⟨hal-03562497⟩ Plus de détails...
The simulation of fusion plasmas in realistic magnetic configurations and tokamak ge-1 ometries still requires the development of advanced numerical algorithms owing to the complexity 2 of the problem. In this context, we propose a Hybrid Discontinuous Galerkin (HDG) method to 3 solve 2D transport fluid equations in realistic magnetic and tokamak wall geometries. This high-4 order solver can handle magnetic equilibrium free structured and unstructured meshes allowing a 5 much more accurate discretization of the plasma facing components than current solvers based 6 on magnetic field aligned methods associated to finite-differences (volumes) discretization. In 7 addition, the method allows to handle realistic magnetic equilibrium, eventually non steady, a 8 critical point in the modelling of full discharges including ramp up and ramp down phases. In 9 this paper, we introduce the HDG algorithm with a special focus on recent developments related 10 to the treatment of the cross-field diffusive terms, and to an adaptive mesh refinement technique 11 improving the numerical efficiency and robustness of the scheme. The updated solver is verified 12 with a manufactured solution method, and numerical tests are provided to illustrate the new 13 capabilities of the code.
Giacomo Piraccini, Marcello Capasso, Manuel Scotto d'Abusco, Giorgio Giorgiani, Frédéric Schwander, et al.. Recent upgrades in a 2D turbulent transport solver based on a hybrid discontinuous Galerkin method for the simulation of fusion plasma in tokamak. Fluids, 2022, ⟨10.3390/fluids7020063⟩. ⟨hal-03562497⟩
Florian Renard, Gauthier Wissocq, Jean-François Boussuge, Pierre Sagaut. A linear stability analysis of compressible hybrid lattice Boltzmann methods. Journal of Computational Physics, 2021, 446, pp.110649. ⟨10.1016/j.jcp.2021.110649⟩. ⟨hal-03514639⟩ Plus de détails...
An original spectral study of the compressible hybrid lattice Boltzmann method (HLBM) on standard lattice is proposed. In this framework, the mass and momentum equations are addressed using the lattice Boltzmann method (LBM), while finite difference (FD) schemes solve an energy equation. Both systems are coupled with each other thanks to an ideal gas equation of state. This work aims at answering some questions regarding the numerical stability of such models, which strongly depends on the choice of numerical parameters. To this extent, several one- and two-dimensional HLBM classes based on different energy variables, formulations (primitive or conservative), collision terms and numerical schemes are scrutinized. Once appropriate corrective terms introduced, it is shown that all continuous HLBM classes recover the Navier-Stokes-Fourier behavior in the linear approximation. However, striking differences arise between HLBM classes when their discrete counterparts are analyzed. Multiple instability mechanisms arising at relatively high Mach number are pointed out and two exhaustive stabilization strategies are introduced: (1) decreasing the time step by changing the reference temperature T-r and (2) introducing a controllable numerical dissipation a via the collision operator. A complete parametric study reveals that only HLBM classes based on the primitive and conservative entropy equations are found usable for compressible applications. Finally, an innovative study of the macroscopic modal composition of the entropy classes is conducted. Through this study, two original phenomena, referred to as shear-to-entropy and entropy-to-shear transfers, are highlighted and confirmed on standard two-dimensional test cases. (C) 2021 Elsevier Inc. All rights reserved.
Florian Renard, Gauthier Wissocq, Jean-François Boussuge, Pierre Sagaut. A linear stability analysis of compressible hybrid lattice Boltzmann methods. Journal of Computational Physics, 2021, 446, pp.110649. ⟨10.1016/j.jcp.2021.110649⟩. ⟨hal-03514639⟩
Thomas Astoul, Gauthier Wissocq, Jean-François Boussuge, Alois Sengissen, Pierre Sagaut. Lattice Boltzmann method for computational aeroacoustics on non-uniform meshes: A direct grid coupling approach. Journal of Computational Physics, 2021, 447, pp.110667. ⟨10.1016/j.jcp.2021.110667⟩. ⟨hal-03514616⟩ Plus de détails...
The present study proposes an accurate lattice Boltzmann direct coupling algorithm, well suited for industrial purposes, making it highly valuable for aeroacoustic applications. It is indeed known that the convection of vortical structures across a grid refinement interface, where cell size is abruptly doubled, is likely to generate spurious noise that may corrupt the solution over the whole computational domain. This issue becomes critical in the case of aeroacoustic simulations, where accurate pressure estimations are of paramount importance. Consequently, any interfering noise that may pollute the acoustic predictions must be reduced. The proposed grid refinement algorithm differs from conventionally used ones, in which an overlapping mesh layer is considered. Instead, it provides a direct connection allowing a tighter link between fine and coarse grids, especially with the use of a coherent equilibrium function shared by both grids. Moreover, the direct coupling makes the algorithm more local and prevents the duplication of points, which might be detrimental for massive parallelization. This work follows our first study (Astoul et al. 2020 [1]) on the deleterious effect of non-hydrodynamic modes crossing mesh transitions, which can be addressed using an appropriate collision model: the hybrid recursive regularization. The grid coupling algorithm is assessed in the context of computational aeroacoustics and compared to a widely-used cell-vertex algorithm. The validation benchmark includes the simulation of (1) an acoustic pulse, (2) a vortex transport by a mean flow, and finally, (3) a turbulent circular cylinder wake flow at high Reynolds number. In the end, the proposed approach is proven to drastically reduced the spurious noise generated at grid interfaces, hence, paving the way for accurate and efficient aeroacoustic simulations based on lattice Boltzmann methods. (C) 2021 Elsevier Inc. All rights reserved.
Thomas Astoul, Gauthier Wissocq, Jean-François Boussuge, Alois Sengissen, Pierre Sagaut. Lattice Boltzmann method for computational aeroacoustics on non-uniform meshes: A direct grid coupling approach. Journal of Computational Physics, 2021, 447, pp.110667. ⟨10.1016/j.jcp.2021.110667⟩. ⟨hal-03514616⟩
