/Publications /Publication of team Instabilities, Turbulence & Control
Téléphone: +33(0)491 11 85 47
Publication of team Instabilities, Turbulence & Control
2023
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⟩
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⟩
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-04064988⟩ 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-04064988⟩
Journal: Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability
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⟩
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
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⟩
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
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
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
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⟩
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
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⟩
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
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⟩
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
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
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...
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⟩
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⟩
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⟩
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⟩
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⟩
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, In press. ⟨hal-03842177⟩ Plus de détails...
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, In press. ⟨hal-03842177⟩
Journal: Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability
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⟩
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⟩
T. Coratger, G. Farag, S. Zhao, Pierre Boivin, P. Sagaut. Large-eddy lattice-Boltzmann modeling of transonic flows. Physics of Fluids, 2021, 33 (11), pp.115112. ⟨10.1063/5.0064944⟩. ⟨hal-03424286⟩ Plus de détails...
T. Coratger, G. Farag, S. Zhao, Pierre Boivin, P. Sagaut. Large-eddy lattice-Boltzmann modeling of transonic flows. Physics of Fluids, 2021, 33 (11), pp.115112. ⟨10.1063/5.0064944⟩. ⟨hal-03424286⟩
Felix Marlow, Jérôme Jacob, Pierre Sagaut. A multidisciplinary model coupling Lattice-Boltzmann-based CFD and a Social Force Model for the simulation of pollutant dispersion in evacuation situations. Building and Environment, 2021, 205, pp.108212. ⟨10.1016/j.buildenv.2021.108212⟩. ⟨hal-03514660⟩ Plus de détails...
In closed rooms with limited convection human motion can considerably affect the airflow and thus the dispersion of pollutant. However, in Computational Fluid Dynamics (CFD) simulations on air quality and safety for human beings this effect is generally not considered, which is mainly due to a lack of a well-founded and detailed estimation of the human behavior and the high computational cost of taking into account moving objects in CFD meshes. This work addresses this issue by coupling multidisciplinary methods to allow for a more realistic simulation of pollutant dispersion by taking into account the influence of human movements. A Social Force Model predicts trajectory and speed of each person moving in a complex environment. A lattice Boltzmann-based CFD tool provides a Large Eddy Simulation of the unsteady turbulent airflow with pollutant dispersion and thermal effects. And an Actuator Line Model supplies the CFD tool with body forces that mimic the impact of moving objects on the airflow, thus, avoiding computationally expensive dynamic meshing. The capability of the coupled model is demonstrated on three realistic evacuation scenarios with various pollutant sources and a wide range of scales (dimension from 10 to 100 m, occupation from 10 to 6000 persons). The results allow to access instantaneous environmental parameters like pollutant concentration for each person during the course of the evacuation, enabling the assessment of associated health risks.
Felix Marlow, Jérôme Jacob, Pierre Sagaut. A multidisciplinary model coupling Lattice-Boltzmann-based CFD and a Social Force Model for the simulation of pollutant dispersion in evacuation situations. Building and Environment, 2021, 205, pp.108212. ⟨10.1016/j.buildenv.2021.108212⟩. ⟨hal-03514660⟩
Guanxiong Wang, Lincheng Xu, Eric Serre, Pierre Sagaut. Large temperature difference heat dominated flow simulations using a pressure-based lattice Boltzmann method with mass correction. Physics of Fluids, 2021, 33 (11), pp.116107. ⟨10.1063/5.0073178⟩. ⟨hal-03438869⟩ Plus de détails...
This paper addresses simulation of heat dominated compressible flows in a closed cavity using a pressure-based lattice Boltzmann (LB) method, in which thermal effects are modeled by applying a pressure-featured zero-order moment of distribution functions. A focus is made on the conservation of mass at boundary nodes, which is a challenging issue that significantly complicated by the density-decoupled zero-order moment here. The mass leakage at boundary nodes is mathematically quantified, which enables an efficient local mass correction scheme. The performance of this solver is assessed by simulating buoyancy-driven flows in a closed deferentially heated cavity with large temperature differences (non-Boussinesq) at Rayleigh numbers ranging from 103 to 107. Simulations show that mass leakage at solid walls in such configurations is a critical issue to obtain reliable solutions, and it eventually leads to simulations overflow when the cavity is inclined. The proposed mass correction scheme is, however, shown to be effective to control the mass leakage and get accurate solutions. Thus, associated with the proposed mass conservation scheme, the pressure-based LB method becomes reliable to study natural convection dominated flows at large temperature differences in closed geometries with mesh aligned boundaries or not
Guanxiong Wang, Lincheng Xu, Eric Serre, Pierre Sagaut. Large temperature difference heat dominated flow simulations using a pressure-based lattice Boltzmann method with mass correction. Physics of Fluids, 2021, 33 (11), pp.116107. ⟨10.1063/5.0073178⟩. ⟨hal-03438869⟩
Felix Marlow, Jérôme Jacob, Pierre Sagaut. A multidisciplinary model coupling Lattice-Boltzmann-based CFD and a Social Force Model for the simulation of pollutant dispersion in evacuation situations. Building and Environment, 2021, 205, pp.108212. ⟨10.1016/j.buildenv.2021.108212⟩. ⟨hal-03597658⟩ Plus de détails...
In closed rooms with limited convection human motion can considerably affect the airflow and thus the dispersion of pollutant. However, in Computational Fluid Dynamics (CFD) simulations on air quality and safety for human beings this effect is generally not considered, which is mainly due to a lack of a well-founded and detailed estimation of the human behavior and the high computational cost of taking into account moving objects in CFD meshes. This work addresses this issue by coupling multidisciplinary methods to allow for a more realistic simulation of pollutant dispersion by taking into account the influence of human movements. A Social Force Model predicts trajectory and speed of each person moving in a complex environment. A lattice Boltzmann-based CFD tool provides a Large Eddy Simulation of the unsteady turbulent airflow with pollutant dispersion and thermal effects. And an Actuator Line Model supplies the CFD tool with body forces that mimic the impact of moving objects on the airflow, thus, avoiding computationally expensive dynamic meshing. The capability of the coupled model is demonstrated on three realistic evacuation scenarios with various pollutant sources and a wide range of scales (dimension from 10 to 100 m, occupation from 10 to 6000 persons). The results allow to access instantaneous environmental parameters like pollutant concentration for each person during the course of the evacuation, enabling the assessment of associated health risks.
Felix Marlow, Jérôme Jacob, Pierre Sagaut. A multidisciplinary model coupling Lattice-Boltzmann-based CFD and a Social Force Model for the simulation of pollutant dispersion in evacuation situations. Building and Environment, 2021, 205, pp.108212. ⟨10.1016/j.buildenv.2021.108212⟩. ⟨hal-03597658⟩
H. Bufferand, J. Bucalossi, G. Ciraolo, G. Falchetto, A. Gallo, et al.. Progress in edge plasma turbulence modelling hierarchy of models from 2D transport application to 3D fluid simulations in realistic tokamak geometry. Nuclear Fusion, 2021, 61 (11), pp.116052. ⟨10.1088/1741-4326/ac2873⟩. ⟨hal-03377162⟩ Plus de détails...
This contribution presents the recent effort at CEA and French federation for Fusion to simulate edge plasma transport with the new code SOLEDGE3X. The latter can be used both as a 2D transport code or as a 3D turbulence code. It makes possible simulating edge plasma up to the first wall including the complex wall geometry. It also includes neutral recycling and impurity sputtering, seeding and transport. In order to improve turbulence description in transport simulation, a reduced model for turbulence intensity prediction has been derived and implemented, based on "kepsilon" like models from the neutral fluid community. Applications to a JET L-mode detached plasma and to a WEST plasma are used as illustration of the code abilities
H. Bufferand, J. Bucalossi, G. Ciraolo, G. Falchetto, A. Gallo, et al.. Progress in edge plasma turbulence modelling hierarchy of models from 2D transport application to 3D fluid simulations in realistic tokamak geometry. Nuclear Fusion, 2021, 61 (11), pp.116052. ⟨10.1088/1741-4326/ac2873⟩. ⟨hal-03377162⟩
Georis Billo, Michel Belliard, Pierre Sagaut. A Finite Element Penalized Direct Forcing Immersed Boundary Method for infinitely thin obstacles in a dilatable flow. Computers & Mathematics with Applications, 2021, 99, pp.292-304. ⟨10.1016/j.camwa.2021.08.005⟩. ⟨hal-03514671⟩ Plus de détails...
In the framework of the development of new passive safety systems for the second and third generations of nuclear reactors, the numerical simulations, involving complex turbulent two-phase flows around thin or massive inflow obstacles, are privileged tools to model, optimize and assess new design shapes. In order to match industrial demands, computational fluid dynamics tools must be the fastest, most accurate and most robust possible. To face this issue, we have chosen to solve the Navier-Stokes equations using a projection scheme for a mixture fluid coupled with an Immersed Boundary (IB) approach: the penalized direct forcing method - a technique whose characteristics inherit from both penalty and immersed boundary methods - adapted to infinitely thin obstacles and to a Finite Element (FE) formulation. Various IB conditions (slip, no-slip or Neumann) for the velocity on the IB can be managed by imposing Dirichlet values in the vicinity of the thin obstacles. To deal with these imposed Dirichlet velocities, we investigated two variants: one in which we use the obstacle velocity and another one in which we use linear interpolations based on discrete geometrical properties of the IB (barycenters and normal vectors) and the FE basis functions. This last variant is motivated by an increase of the accuracy/computation time ratio for coarse meshes. As a first step, concerning academic test cases for one-phase dilatable-fluid laminar flows, the results obtained via those two variants are in good agreement with analytical and experimental data. Moreover, when compared to each other, the linear interpolation variant increases the spatial order of convergence as expected. An industrial test case illustrates the advantages and drawbacks of this approach. In a shortcoming second step, to face two-phase turbulent fluid simulations, some methodology modifications will be considered such as adapting the projection scheme to low-compressible fluid and immersed wall-law boundary conditions.
Georis Billo, Michel Belliard, Pierre Sagaut. A Finite Element Penalized Direct Forcing Immersed Boundary Method for infinitely thin obstacles in a dilatable flow. Computers & Mathematics with Applications, 2021, 99, pp.292-304. ⟨10.1016/j.camwa.2021.08.005⟩. ⟨hal-03514671⟩
Journal: Computers & Mathematics with Applications
Georis Billo, Michel Belliard, Pierre Sagaut. A Finite Element Penalized Direct Forcing Immersed Boundary Method for infinitely thin obstacles in a dilatable flow. Computers & Mathematics with Applications, 2021, 99, pp.292-304. ⟨10.1016/j.camwa.2021.08.005⟩. ⟨hal-03596009⟩ Plus de détails...
In the framework of the development of new passive safety systems for the second and third generations of nuclear reactors, the numerical simulations, involving complex turbulent two-phase flows around thin or massive inflow obstacles, are privileged tools to model, optimize and assess new design shapes. In order to match industrial demands, computational fluid dynamics tools must be the fastest, most accurate and most robust possible. To face this issue, we have chosen to solve the Navier-Stokes equations using a projection scheme for a mixture fluid coupled with an Immersed Boundary (IB) approach: the penalized direct forcing method-a technique whose characteristics inherit from both penalty and immersed boundary methods-adapted to infinitely thin obstacles and to a Finite Element (FE) formulation. Various IB conditions (slip, no-slip or Neumann) for the velocity on the IB can be managed by imposing Dirichlet values in the vicinity of the thin obstacles. To deal with these imposed Dirichlet velocities, we investigated two variants: one in which we use the obstacle velocity and another one in which we use linear interpolations based on discrete geometrical properties of the IB (barycenters and normal vectors) and the FE basis functions. This last variant is motivated by an increase of the accuracy/computation time ratio for coarse meshes. As a first step, concerning academic test cases for one-phase dilatable-fluid laminar flows, the results obtained via those two variants are in good agreement with analytical and experimental data. Moreover, when compared to each other, the linear interpolation variant increases the spatial order of convergence as expected. An industrial test case illustrates the advantages and drawbacks of this approach. In a shortcoming second step, to face two-phase turbulent fluid simulations, some methodology modifications will be considered such as adapting the projection scheme to low-compressible fluid and immersed wall-law boundary conditions.
Georis Billo, Michel Belliard, Pierre Sagaut. A Finite Element Penalized Direct Forcing Immersed Boundary Method for infinitely thin obstacles in a dilatable flow. Computers & Mathematics with Applications, 2021, 99, pp.292-304. ⟨10.1016/j.camwa.2021.08.005⟩. ⟨hal-03596009⟩
Journal: Computers & Mathematics with Applications
Pierre Magnico. Molecular dynamics study on water and hydroxide transfer mechanisms in PSU-g-alkyl-TMA membranes at low hydration: Effect of side chain length. International Journal of Hydrogen Energy, 2021, 46 (68), pp.33915-33933. ⟨10.1016/j.ijhydene.2021.07.081⟩. ⟨hal-03358207⟩ Plus de détails...
Molecular dynamics simulations with anion exchange membranes (alkyl trimethyl ammonium grafted onto polysulfone) are performed to investigate the influence of the spacer length on the transport properties, on the molecular exchange mechanisms between the functional group and the aqueous phase and on the hydrogen bond network. This is especially insightful that in this work the hydration number is small. In this condition the aqueous phase must be thought as an assembly of small clusters. The results show an unexpected dependence of the water and hydroxide (OH) diffusivity on the temperature and the water uptake. The distribution of the cluster size bonded to OH explain partially the OH diffusivity. “Hopping” and “caging” motions are observed with the self-part of the Van Hove functions even at high temperature. The characteristic time of the survival probability correlation function around the functional groups is a decreasing function of the alkyl length.
Pierre Magnico. Molecular dynamics study on water and hydroxide transfer mechanisms in PSU-g-alkyl-TMA membranes at low hydration: Effect of side chain length. International Journal of Hydrogen Energy, 2021, 46 (68), pp.33915-33933. ⟨10.1016/j.ijhydene.2021.07.081⟩. ⟨hal-03358207⟩
S. Baschetti, H. Bufferand, G. Ciraolo, Ph Ghendrih, E. Serre, et al.. Self-consistent cross-field transport model for core and edge plasma transport. Nuclear Fusion, 2021, 61 (10), pp.106020. ⟨10.1088/1741-4326/ac1e60⟩. ⟨hal-03380310⟩ Plus de détails...
A two-equation model to self-consistently determine cross-field fluxes in the edge and scrape-off layer region of diverted plasma is used to complete 2D mean-field edge transport description of plasma wall interaction. Inspired by the Reynolds Average Navier-Stokes simulations for neutral fluids, this model is based on the local evolution of the turbulent kinetic energy κ and its dissipation rate ε. These two equations are algebraically derived for RANS modeling and are very slightly modified and adapted to describe self-consistent plasma turbulent transport. The general features of the model are discussed and bridged to the well-known predatorprey and quasilinear models commonly used to investigate plasma transport. Specific closures are proposed based on the interchange turbulence. Results of the 1D model are confronted to experimental evidence by analyzing the computed SOL width and comparing the results to the existing scaling law for L-mode plasmas. Introducing a dependence on the shear of large scale flows, typically the zonal flows, 1D simulations can exhibit an H-mode like transition when increasing the input power, generating an increased stored energy thanks to an interface barrier located at the separatrix. Further 2D plasma-wall interaction simulations for WEST are analyzed that show a good match with the experimental profiles, as well as a ballooned transport driving turbulent transport in the divertor SOL and nearly no transport in the private flux region. The SOL width of WEST is also recovered. These results show the remarkable capability of the κ-ε model to capture key aspects of the physics of turbulent transport throughout the plasma knowing that a unique scalar free parameter is available to tune cross field transport in the whole 2D cross section of the plasma.
S. Baschetti, H. Bufferand, G. Ciraolo, Ph Ghendrih, E. Serre, et al.. Self-consistent cross-field transport model for core and edge plasma transport. Nuclear Fusion, 2021, 61 (10), pp.106020. ⟨10.1088/1741-4326/ac1e60⟩. ⟨hal-03380310⟩
S. Di Genova, A. Gallo, N. Fedorczak, H. Yang, G. Ciraolo, et al.. Modelling of tungsten contamination and screening in WEST plasma discharges. Nuclear Fusion, 2021, 61 (10), pp.106019. ⟨10.1088/1741-4326/ac2026⟩. ⟨hal-03380329⟩ Plus de détails...
The WEST experiment is currently operating with tungsten plasma-facing components and testing ITER-like divertor monoblocks. In order to support WEST experiments interpretation, numerical analyses were carried out. Starting from WEST experimental data, realistic background plasma conditions were reproduced through SolEdge-EIRENE and used as input for ERO2.0 simulations to investigate tungsten migration. Tungsten contamination due to the different plasma-facing components was modelled under different plasma conditions, highlighting a non-negligible contribution of tungsten coming from the tokamak main chamber. Tungsten penetration factor was computed and used as an indication for tungsten screening by the background plasma at the different tokamak plasma-facing components. Simulations showed the main chamber components to be very weakly screened. Light impurities charge was showed to influence not only tungsten sputtering, but also its probability to enter the confined plasma. Simulations results indicated that even when the tungsten source is not heavily influenced by self-sputtering, contamination of the confined plasma can be strongly impacted by it in low density background plasma conditions. Finally, a one-to-one comparison between tungsten visible spectroscopy at the lower divertor from experimental data and from synthetic diagnostics was performed, showing that it is possible to reproduce a realistic lower divertor signal following experimental evidence on light impurities asymmetry between the targets.
S. Di Genova, A. Gallo, N. Fedorczak, H. Yang, G. Ciraolo, et al.. Modelling of tungsten contamination and screening in WEST plasma discharges. Nuclear Fusion, 2021, 61 (10), pp.106019. ⟨10.1088/1741-4326/ac2026⟩. ⟨hal-03380329⟩
Johan Degrigny, Shang-Gui Cai, Jean-François Boussuge, Pierre Sagaut. Improved wall model treatment for aerodynamic flows in LBM. Computers and Fluids, 2021, 227, pp.105041. ⟨10.1016/j.compfluid.2021.105041⟩. ⟨hal-03326170⟩ Plus de détails...
The article deals with an improved treatment of wall models for the simulation of turbulent flows in the framework of Immersed Wall Boundaries on Cartesian grids. The emphasis is put on the implementa-tion in a Lattice-Boltzmann Method solver without loss of generality, since the proposed approach can be used in Navier-Stokes-based solvers in a straightforward way. The proposed improved wall model im-plementation relies on the combination of several key elements, namely i) the removal of grid points too close to the solid surface and ii) an original computation of wall normal velocity gradient and iii) the interpolation scheme. The new method is successfully assessed considering URANS simulations focusing on steady solutions of the Zero Pressure Gradient turbulent flat plate boundary layer and the turbulent flow around a NACA0012 airfoil at several angles of attack.
Johan Degrigny, Shang-Gui Cai, Jean-François Boussuge, Pierre Sagaut. Improved wall model treatment for aerodynamic flows in LBM. Computers and Fluids, 2021, 227, pp.105041. ⟨10.1016/j.compfluid.2021.105041⟩. ⟨hal-03326170⟩
Isabelle Cheylan, Julien Favier, Pierre Sagaut. Immersed boundary conditions for moving objects in turbulent flows with the lattice-Boltzmann method. Physics of Fluids, 2021, 33 (9), pp.095101. ⟨10.1063/5.0062575⟩. ⟨hal-03597108⟩ Plus de détails...
An immersed boundary method is coupled to a turbulent wall model and Large Eddy Simulation, within the Lattice-Boltzmann framework. The method is able to handle arbitrarily moving objects immersed in a high Reynolds number flow and to accurately capture the shear layer and near wall effects. We perform a thorough numerical study which validates the numerical method on a set of test-cases of increasing complexity, in order to demonstrate the application of this method to industrial conditions. The robustness and accuracy of the method are assessed first in a static laminar configuration, then in a mobile laminar case, and finally in a static and oscillating turbulent simulation. In all cases, the proposed method shows good results compared to the available data in the literature.
Isabelle Cheylan, Julien Favier, Pierre Sagaut. Immersed boundary conditions for moving objects in turbulent flows with the lattice-Boltzmann method. Physics of Fluids, 2021, 33 (9), pp.095101. ⟨10.1063/5.0062575⟩. ⟨hal-03597108⟩
Isabelle Cheylan, Julien Favier, Pierre Sagaut. Immersed boundary conditions for moving objects in turbulent flows with the lattice-Boltzmann method. Physics of Fluids, 2021, 33 (9), pp.095101. ⟨10.1063/5.0062575⟩. ⟨hal-03514710⟩ Plus de détails...
An immersed boundary method is coupled to a turbulent wall model and Large Eddy Simulation, within the Lattice-Boltzmann framework. The method is able to handle arbitrarily moving objects immersed in a high Reynolds number flow and to accurately capture the shear layer and near wall effects. We perform a thorough numerical study which validates the numerical method on a set of test-cases of increasing complexity, in order to demonstrate the application of this method to industrial conditions. The robustness and accuracy of the method are assessed first in a static laminar configuration, then in a mobile laminar case, and finally in a static and oscillating turbulent simulation. In all cases, the proposed method shows good results compared to the available data in the literature.
Isabelle Cheylan, Julien Favier, Pierre Sagaut. Immersed boundary conditions for moving objects in turbulent flows with the lattice-Boltzmann method. Physics of Fluids, 2021, 33 (9), pp.095101. ⟨10.1063/5.0062575⟩. ⟨hal-03514710⟩
M Raghunathan, Y Marandet, H Bufferand, G Ciraolo, Ph Ghendrih, et al.. Generalized Collisional Fluid Theory for Multi-Component, Multi-Temperature Plasma Using The Linearized Boltzmann Collision Operator for Scrape-Off Layer/Edge Applications. Plasma Physics and Controlled Fusion, 2021, 63 (6), pp.064005. ⟨10.1088/1361-6587/abf670⟩. ⟨hal-03384547⟩ Plus de détails...
M Raghunathan, Y Marandet, H Bufferand, G Ciraolo, Ph Ghendrih, et al.. Generalized Collisional Fluid Theory for Multi-Component, Multi-Temperature Plasma Using The Linearized Boltzmann Collision Operator for Scrape-Off Layer/Edge Applications. Plasma Physics and Controlled Fusion, 2021, 63 (6), pp.064005. ⟨10.1088/1361-6587/abf670⟩. ⟨hal-03384547⟩
M. Bahlali, H. Yoo, Julien Favier, Pierre Sagaut. A lattice Boltzmann direct coupling overset approach for the moving boundary problem. Physics of Fluids, 2021, 33 (5), pp.053607. ⟨10.1063/5.0044994⟩. ⟨hal-03326151⟩ Plus de détails...
We propose a new direct coupling scheme based on the overset technique to tackle moving boundary problems within the lattice Boltzmann framework. The scheme is based on the interpolation of distribution functions rather than moments, that is, macroscopic variables, and includes an additional hypothesis ensuring mass and momentum conservation at the interface nodes between fixed and moving grids. The method is assessed considering four test cases and considering both the vortical and the acoustic fields. It is shown that the direct coupling method results are in very good agreement with reference results on a configuration without any moving subdomain. Moreover, it is demonstrated that the direct coupling method provides an improvement of the accuracy of the lattice Boltzmann overset algorithm for aeroacoustics. In particular, a convected vortex test case is studied and reveals that the direct coupling approach leads to a better ability to conserve the vortex structure over time, as well as a reduction in spurious acoustic distorsions at the fixed/moving interface.
M. Bahlali, H. Yoo, Julien Favier, Pierre Sagaut. A lattice Boltzmann direct coupling overset approach for the moving boundary problem. Physics of Fluids, 2021, 33 (5), pp.053607. ⟨10.1063/5.0044994⟩. ⟨hal-03326151⟩
H. Yoo, M. Bahlali, Julien Favier, Pierre Sagaut. A hybrid recursive regularized lattice Boltzmann model with overset grids for rotating geometries. Physics of Fluids, 2021, 33 (5), pp.057113. ⟨10.1063/5.0045524⟩. ⟨hal-03597721⟩ Plus de détails...
Simulating rotating geometries in fluid flows for industrial applications remains a challenging task for general fluid solvers and in particular for the lattice Boltzmann method (LBM) due to inherent stability and accuracy problems. This work proposes an original method based on the widely used overset grids (or Chimera grids) while being integrated with a recent and optimized LBM collision operator, the hybrid recursive regularized model (HRR). The overset grids are used to actualize the rotating geometries where both the rotating and fixed meshes exist simultaneously. In the rotating mesh, the fictitious forces generated from its non-inertial rotating reference frame are taken into account by using a second order discrete forcing term. The fixed and rotating grids communicate with each other through the interpolation of the macroscopic variables. Meanwhile, the HRR collision model is selected to enhance the stability and accuracy properties of the LBM simulations by filtering out redundant higher order non-equilibrium tensors. The robustness of the overset HRR algorithm is assessed on different configurations, undergoing mid-to-high Reynolds number flows, and the method successfully demonstrates its robustness while exhibiting the second order accuracy.
H. Yoo, M. Bahlali, Julien Favier, Pierre Sagaut. A hybrid recursive regularized lattice Boltzmann model with overset grids for rotating geometries. Physics of Fluids, 2021, 33 (5), pp.057113. ⟨10.1063/5.0045524⟩. ⟨hal-03597721⟩
H. Yoo, M. Bahlali, Julien Favier, Pierre Sagaut. A hybrid recursive regularized lattice Boltzmann model with overset grids for rotating geometries. Physics of Fluids, 2021, 33 (5), pp.057113. ⟨10.1063/5.0045524⟩. ⟨hal-03326134⟩ Plus de détails...
Simulating rotating geometries in fluid flows for industrial applications remains a challenging task for general fluid solvers and in particular for the lattice Boltzmann method (LBM) due to inherent stability and accuracy problems. This work proposes an original method based on the widely used overset grids (or Chimera grids) while being integrated with a recent and optimized LBM collision operator, the hybrid recursive regularized model (HRR). The overset grids are used to actualize the rotating geometries where both the rotating and fixed meshes exist simultaneously. In the rotating mesh, the fictitious forces generated from its non-inertial rotating reference frame are taken into account by using a second order discrete forcing term. The fixed and rotating grids communicate with each other through the interpolation of the macroscopic variables. Meanwhile, the HRR collision model is selected to enhance the stability and accuracy properties of the LBM simulations by filtering out redundant higher order non-equilibrium tensors. The robustness of the overset HRR algorithm is assessed on different configurations, undergoing mid-to-high Reynolds number flows, and the method successfully demonstrates its robustness while exhibiting the second order accuracy.
H. Yoo, M. Bahlali, Julien Favier, Pierre Sagaut. A hybrid recursive regularized lattice Boltzmann model with overset grids for rotating geometries. Physics of Fluids, 2021, 33 (5), pp.057113. ⟨10.1063/5.0045524⟩. ⟨hal-03326134⟩
Elisa Buffa, Jérôme Jacob, Pierre Sagaut. Lattice-Boltzmann-based large-eddy simulation of high-rise building aerodynamics with inlet turbulence reconstruction. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 212, pp.104560. ⟨10.1016/j.jweia.2021.104560⟩. ⟨hal-03596056⟩ Plus de détails...
Boltzmann-based Large-Eddy Simulation approach for wind load prediction on high-rise building is proposed and validated. An extension of the original incompressible Synthetic Eddy Method to reconstruct inlet turbulence is proposed within the Lattice-Boltzmann framework, including a low-noise frozen density variant. Extensive successful comparisons with experimental data are carried out, for both quantities defined on the building surface and in its wake. A detailed sensitivity analysis of the results with respect to inlet turbulence reconstruction, boundary conditions at the building surface and grid resolution is also provided. An almost unique set of comparisons with experimental data is presented, including mean and rms values, spectra, but also peak values of pressure at the building surface.
Elisa Buffa, Jérôme Jacob, Pierre Sagaut. Lattice-Boltzmann-based large-eddy simulation of high-rise building aerodynamics with inlet turbulence reconstruction. Journal of Wind Engineering and Industrial Aerodynamics, 2021, 212, pp.104560. ⟨10.1016/j.jweia.2021.104560⟩. ⟨hal-03596056⟩
Journal: Journal of Wind Engineering and Industrial Aerodynamics
Shang-Gui Cai, Pierre Sagaut. Explicit wall models for large eddy simulation. Physics of Fluids, 2021, 33 (4), pp.041703. ⟨10.1063/5.0048563⟩. ⟨hal-03597083⟩ Plus de détails...
Algebraic explicit wall models covering the entire inner region of the turbulent boundary layer are proposed to reduce the computational effort for large eddy simulation of wall-bounded turbulent flows. The proposed formulas are given in closed forms with either logarithmicor power-function-based laws of the wall, allowing straightforward evaluation of the friction velocity on near wall grids independent of their locations in the turbulent boundary layer. The performance of the proposed models is demonstrated by the wall modeled large eddy simulation of a turbulent plane channel flow.
Shang-Gui Cai, Pierre Sagaut. Explicit wall models for large eddy simulation. Physics of Fluids, 2021, 33 (4), pp.041703. ⟨10.1063/5.0048563⟩. ⟨hal-03597083⟩
Jérémie Janin, Fabien Duval, Christophe Friess, Pierre Sagaut. A new linear forcing method for isotropic turbulence with controlled integral length scale. Physics of Fluids, 2021, 33 (4), pp.045127. ⟨10.1063/5.0045818⟩. ⟨hal-03326165⟩ Plus de détails...
Turbulence is a common feature to all flows that surround us. Despite its ubiquity, particularly in industrial flows, it is very difficult to provide a mathematical framework to the generation of turbulent eddies. Several techniques have been proposed which are able to reproduce the main features of turbulent flows, such as realistic pressure and velocity fluctuations, exhibiting proper space- and time-correlations. These techniques are usually first evaluated upon sustained homogeneous isotropic turbulence by introducing body forces to the Navier-Stokes equations. Among these techniques, Lundgren suggested a successful forcing, applied in physical space. The latter approach unfortunately lacks predicting the integral length scale of turbulence. The present study provides a forcing method based on a reconstruction approach which consists in building fluctuations with a turbulent synthetic velocity field based on a prescribed energy spectrum model. The proposed approach is assessed by performing large-eddy simulations of a sustained homogeneous isotropic turbulence in a triply periodic box of size L = 2pi. Properties of the new forcing technique are discussed, drawing on both spatial and time correlations and also on the shape of energy spectrum together with the level of resolved turbulent kinetic energy. A special attention is put on the control of resolved turbulent energy. In this framework, an efficient selective forcing technique is derived, making use of spectral space features. The results show that the proposed approach allows to drive efficiently the resolved kinetic energy towards its target value while preserving the integral length scale independent of the domain size. It is observed that the resulting longitudinal length scale is overestimated by 13%, while the two-time correlations are recovered when using stochastic frequencies.
Jérémie Janin, Fabien Duval, Christophe Friess, Pierre Sagaut. A new linear forcing method for isotropic turbulence with controlled integral length scale. Physics of Fluids, 2021, 33 (4), pp.045127. ⟨10.1063/5.0045818⟩. ⟨hal-03326165⟩
E. Laribi, E. Serre, P. Tamain, H. Yang. Impact of negative triangularity on edge plasma transport and turbulence in TOKAM3X simulations. Nuclear Materials and Energy, 2021, pp.101012. ⟨10.1016/j.nme.2021.101012⟩. ⟨hal-03214958⟩ Plus de détails...
The impact of triangularity on edge plasma transport and turbulence is addressed from full 3D turbulence simulations performed with TOKAM3X. Flux driven fluid simulations are run on analytical magnetic equilibria generated with positive and negative triangularity δ in a bottom limiter configuration. The conservation of the energy is assured by the increase of the bottom limiter radial position from δ > 0 to δ < 0. Changing the triangularity impacts both the plasma equilibrium and the turbulence. In particular, negative triangularity leads to a reduction of the density and electron temperature decay lengths in agreement with the literature. Concerning the turbulence, in all the simulations, it remains ballooned with an enhanced level of fluctuations at low field side in comparison to the high field one. Moreover, no clear trend is visible on the relative level of fluctuations of both density and electron temperature in the CFR whereas an enhancement (resp. reduction) is visible in the scrape-off layer at the low field side midplane for the negative (resp. positive) triangularity simulations. This behaviour differs from TCV and DIII-D measurements which show the benefit of negative triangularity in terms of turbulence reduction and increased confinement. However, no conclusion is drawn from our preliminary study concerning the impact of triangularity on the turbulent transport. Change in triangularity impacts many simulation control parameters, as in the experiments, and that the analysis of its impact alone on the dynamics of the plasma is not obvious in this configuration.
E. Laribi, E. Serre, P. Tamain, H. Yang. Impact of negative triangularity on edge plasma transport and turbulence in TOKAM3X simulations. Nuclear Materials and Energy, 2021, pp.101012. ⟨10.1016/j.nme.2021.101012⟩. ⟨hal-03214958⟩
E. Laribi, E. Serre, P. Tamain, H. Yang. Impact of negative triangularity on edge plasma transport and turbulence in TOKAM3X simulations. Nuclear Materials and Energy, 2021, pp.101012. ⟨10.1016/j.nme.2021.101012⟩. ⟨hal-03214958⟩ Plus de détails...
The impact of triangularity on edge plasma transport and turbulence is addressed from full 3D turbulence simulations performed with TOKAM3X. Flux driven fluid simulations are run on analytical magnetic equilibria generated with positive and negative triangularity δ in a bottom limiter configuration. The conservation of the energy is assured by the increase of the bottom limiter radial position from δ > 0 to δ < 0. Changing the triangularity impacts both the plasma equilibrium and the turbulence. In particular, negative triangularity leads to a reduction of the density and electron temperature decay lengths in agreement with the literature. Concerning the turbulence, in all the simulations, it remains ballooned with an enhanced level of fluctuations at low field side in comparison to the high field one. Moreover, no clear trend is visible on the relative level of fluctuations of both density and electron temperature in the CFR whereas an enhancement (resp. reduction) is visible in the scrape-off layer at the low field side midplane for the negative (resp. positive) triangularity simulations. This behaviour differs from TCV and DIII-D measurements which show the benefit of negative triangularity in terms of turbulence reduction and increased confinement. However, no conclusion is drawn from our preliminary study concerning the impact of triangularity on the turbulent transport. Change in triangularity impacts many simulation control parameters, as in the experiments, and that the analysis of its impact alone on the dynamics of the plasma is not obvious in this configuration.
E. Laribi, E. Serre, P. Tamain, H. Yang. Impact of negative triangularity on edge plasma transport and turbulence in TOKAM3X simulations. Nuclear Materials and Energy, 2021, pp.101012. ⟨10.1016/j.nme.2021.101012⟩. ⟨hal-03214958⟩
S. Guo, Y. Feng, Pierre Sagaut. On the use of conservative formulation of energy equation in hybrid compressible lattice Boltzmann method. Computers and Fluids, 2021, 219, pp.104866. ⟨10.1016/j.compfluid.2021.104866⟩. ⟨hal-03326128⟩ Plus de détails...
Effect of density variations on mass conservation properties is widely recognized in the lattice Boltzmann method (LBM), thus non-conservative form of scalar transport equation was commonly adopted within the framework of hybrid LBM. Focusing on the compressible hybrid LBM, mass conservation and its effect on energy conservation equation are studied in this paper. Starting from the analysis on mass conservation law recovered by LBM, the consistency between conservative and non-conservative formulations of energy conservation equation based on various thermodynamic variables and lattice Boltzmann equation is addressed. Driven by the theoretical analysis, a set of modified consistent energy equations in entropy and internal energy form is derived to reduce the error terms and improve the consistency. The theoretical analysis and modified energy equations are intensively evaluated by several numerical test cases, e.g., the isentropic vortex convection, three-dimensional compressible Taylor-Green vortex and shock-vortex interaction.
S. Guo, Y. Feng, Pierre Sagaut. On the use of conservative formulation of energy equation in hybrid compressible lattice Boltzmann method. Computers and Fluids, 2021, 219, pp.104866. ⟨10.1016/j.compfluid.2021.104866⟩. ⟨hal-03326128⟩
S. Guo, Y. Feng, Pierre Sagaut. On the use of conservative formulation of energy equation in hybrid compressible lattice Boltzmann method. Computers and Fluids, 2021, 219, pp.104866. ⟨10.1016/j.compfluid.2021.104866⟩. ⟨hal-03597478⟩ Plus de détails...
Effect of density variations on mass conservation properties is widely recognized in the lattice Boltzmann method (LBM), thus non-conservative form of scalar transport equation was commonly adopted within the framework of hybrid LBM. Focusing on the compressible hybrid LBM, mass conservation and its effect on energy conservation equation are studied in this paper. Starting from the analysis on mass conservation law recovered by LBM, the consistency between conservative and non-conservative formulations of energy conservation equation based on various thermodynamic variables and lattice Boltzmann equation is addressed. Driven by the theoretical analysis, a set of modified consistent energy equations in entropy and internal energy form is derived to reduce the error terms and improve the consistency. The theoretical analysis and modified energy equations are intensively evaluated by several numerical test cases, e.g., the isentropic vortex convection, three-dimensional compressible Taylor-Green vortex and shock-vortex interaction.
S. Guo, Y. Feng, Pierre Sagaut. On the use of conservative formulation of energy equation in hybrid compressible lattice Boltzmann method. Computers and Fluids, 2021, 219, pp.104866. ⟨10.1016/j.compfluid.2021.104866⟩. ⟨hal-03597478⟩
Florian Renard, Yongliang Feng, Jean-François Boussuge, Pierre Sagaut. Improved compressible hybrid lattice Boltzmann method on standard lattice for subsonic and supersonic flows. Computers and Fluids, 2021, 219, pp.104867. ⟨10.1016/j.compfluid.2021.104867⟩. ⟨hal-03326159⟩ Plus de détails...
A D2Q9 Hybrid Lattice Boltzmann Method (HLBM) is proposed for the simulation of both compressible subsonic and supersonic flows. This HLBM is an extension of the model of Feng et al. [1], which has been found, via different test cases, to be unstable for supersonic regimes. To circumvent this limitation, we propose:: (1) a new discretization of the lattice closure correction term that makes possible the simulation of supersonic flows, (2) a corrected viscous stress tensor that takes into account polyatomic gases, and (3) a novel discretization of the viscous heat production term fitting with the regularized formalism. The result is a hybrid method that resolves the mass and momentum equations with an LBM algorithm, and resolves the entropy-based energy equation with a finite volume method. This approach fully recovers the physics of the Navier-Stokes-Fourier equations with the ideal gas equation of state, and is valid from subsonic to supersonic regimes. It is then successfully assessed with both smooth flows and flows involving shocks. The proposed model is shown to be an efficient, accurate, and robust alternative to classic Navier-Stokes methods for the simulation of compressible flows.
Florian Renard, Yongliang Feng, Jean-François Boussuge, Pierre Sagaut. Improved compressible hybrid lattice Boltzmann method on standard lattice for subsonic and supersonic flows. Computers and Fluids, 2021, 219, pp.104867. ⟨10.1016/j.compfluid.2021.104867⟩. ⟨hal-03326159⟩
Jérémie Janin, Fabien Duval, Christophe Friess, Pierre Sagaut. A new linear forcing method for isotropic turbulence with controlled integral length scale. Physics of Fluids, American Institute of Physics, 2021, 33 (4), pp.045127. ⟨10.1063/5.0045818⟩. ⟨hal-03326165⟩ Plus de détails...
Turbulence is a common feature to all flows that surround us. Despite its ubiquity, particularly in industrial flows, it is very difficult to provide a mathematical framework for the generation of turbulent eddies. Several methods have been proposed which are able to reproduce realistic features for velocity fluctuations, exhibiting proper space- and time-correlations. Focusing on physical space forcing, these methods are usually first evaluated upon sustained homogeneous isotropic turbulence by introducing a body force to the Navier-Stokes equations. Since the pioneering work of Lundgren, these techniques usually experience difficulties in predicting the integral length scale. The present study provides a forcing through a reconstruction approach which consists in building velocity fluctuations with a prescribed energy spectrum model. The proposed approach is assessed by performing large-eddy simulations of a sustained homogeneous isotropic turbulence in a triply periodic box. Properties of this forcing technique are discussed, drawing on both spatial and time correlations and also on the shape of energy spectrum together with the level of resolved turbulent kinetic energy. A special attention is put on the control of resolved turbulent energy. In this framework, an efficient selective forcing technique is derived, making use of spectral space features. The results show that the proposed approach allows to drive efficiently the resolved kinetic energy toward its target value while preserving the integral length scale independent of the domain size. It is observed that the resulting longitudinal length scale is overestimated by 13%, while the two-time correlations are recovered when using stochastic frequencies.
Jérémie Janin, Fabien Duval, Christophe Friess, Pierre Sagaut. A new linear forcing method for isotropic turbulence with controlled integral length scale. Physics of Fluids, American Institute of Physics, 2021, 33 (4), pp.045127. ⟨10.1063/5.0045818⟩. ⟨hal-03326165⟩
Shang-Gui Cai, Johan Degrigny, Jean-François Boussuge, Pierre Sagaut. Coupling of turbulence wall models and immersed boundaries on Cartesian grids. Journal of Computational Physics, 2021, 429, pp.109995. ⟨10.1016/j.jcp.2020.109995⟩. ⟨hal-03597064⟩ Plus de détails...
An improved coupling of immersed boundary method and turbulence wall models on Cartesian grids is proposed, for producing smooth wall surface pressure and skin friction at high Reynolds numbers. Spurious oscillations are frequently observed on these quantities with most immersed boundary wall modeling methods, especially for the skin friction which is found to be very sensitive to the solid surface's position and orientation against the Cartesian grids. The problem originates from the irregularity of the wall distance on the stair-step grid boundaries where the immersed boundary conditions are applied. To reduce this directional error, several modifications are presented to enhance the near wall solution. First, the commonly used interpolation for the flow velocity is replaced by one for the friction velocity, which has much less variation near wall. The concept of using a fictitious point to retrieve flow fields in the wall normal direction is abandoned and the interpolation is performed in the wall parallel plane with existing fluid points. Secondly, the velocity gradients at the approximated boundary are computed with advanced schemes and the normal gradient of the tangential velocity is reconstructed from the wall laws. To further protect the near wall solution, the normal velocity gradient and the working viscosity from the Spalart-Allmaras turbulence model are enforced by their theoretical solutions in the interior fluid close to the wall. Additionally, various post-processing algorithms for reconstructing wall surface quantities and force integrations are investigated. Other related factors are also discussed for their effects on the results. The validity of present method has been demonstrated through numerical benchmark tests on a flat plate at zero pressure gradient, both aligned and inclined with respect to the grid, as well as aerodynamic cases of NACA 23012 airfoil and NASA trap wing.
Shang-Gui Cai, Johan Degrigny, Jean-François Boussuge, Pierre Sagaut. Coupling of turbulence wall models and immersed boundaries on Cartesian grids. Journal of Computational Physics, 2021, 429, pp.109995. ⟨10.1016/j.jcp.2020.109995⟩. ⟨hal-03597064⟩
The boundary slip error resulting from the interpolation/spreading non-reciprocity of the direct-forcing immersedboundary method is analyzed based on a simple and generic theoretical framework. In explicit implementations, the slip error scales with the Courant number, as predicted by the analysis and confirmed by lattice-Boltzmann simulation results. Using an analytical approximation of the non-reciprocity error, the immersed-boundary force can be corrected in order to prevent boundary slip and flow penetration. This a priori correction leads to a major improvement of the no-slip condition while avoiding any additional computational time or implementation effort.
Simon Gsell, Julien Favier. Direct-forcing immersed-boundary method: a simple correction preventing boundary slip error. Journal of Computational Physics, 2021, 435, pp.110265. ⟨10.1016/j.jcp.2021.110265⟩. ⟨hal-03425864⟩
Jincheng Lou, Jacob Johnston, Tzahi Y. Cath, Denis Martinand, Nils Tilton. Computational fluid dynamics simulations of unsteady mixing in spacer-filled direct contact membrane distillation channels. Journal of Membrane Science, 2021, 622, pp.118931. ⟨10.1016/j.memsci.2020.118931⟩. ⟨hal-03597607⟩ Plus de détails...
Direct contact membrane distillation (DCMD) is a promising means of concentrating brines to their saturation limit. During that process, membrane spacers play a key role in temperature polarization, concentration polarization, and mineral scaling. These interactions are not well understood, because they are difficult to study experimentally and numerically, and the flow regimes are not fully charted. We consequently develop a tailored in-house CFD code that simulates unsteady two-dimensional heat and mass transport in plate-and-frame DCMD systems with cylindrical spacers. The code uses a combination of finite-volume methods in space, projection methods in time, and recent advances in immersed boundary methods for the spacer surfaces. Using the code, we explore how the transition to unsteady laminar vortex shedding affects polarization and permeate production of DCMD systems. We show that the impact of spacers can be explained by examining the various steady and unsteady vortical flow structures generated in the bulk and near the membranes. Overall, we show that though unsteady vortex structures tend to mix temperature polarization layers with the bulk, they are not similarly able to mix the concentration layers. Rather, vortical structures tend to create regions of preferential salt accumulation. In the vortex shedding regime, the net result is that spacers often increase vapor production at the expense of increasing the risk of mineral scaling.
Jincheng Lou, Jacob Johnston, Tzahi Y. Cath, Denis Martinand, Nils Tilton. Computational fluid dynamics simulations of unsteady mixing in spacer-filled direct contact membrane distillation channels. Journal of Membrane Science, 2021, 622, pp.118931. ⟨10.1016/j.memsci.2020.118931⟩. ⟨hal-03597607⟩
Yongliang Feng, Johann Miranda‐fuentes, Shaolong Guo, Jérôme Jacob, Pierre Sagaut. ProLB: A Lattice Boltzmann Solver of Large‐Eddy Simulation for Atmospheric Boundary Layer Flows. Journal of Advances in Modeling Earth Systems, 2021, 13 (3), pp.e2020MS002107. ⟨10.1029/2020MS002107⟩. ⟨hal-03326123⟩ Plus de détails...
A large-eddy simulation tool is developed for simulating the dynamics of atmospheric boundary layers (ABLs) using lattice Boltzmann method (LBM), which is an alternative approach for computational fluid dynamics and proved to be very well suited for the simulation of low-Mach flows. The equations of motion are coupled with the global complex physical models considering the coupling among several mechanisms, namely basic hydro-thermodynamics and body forces related to stratification, Coriolis force, canopy effects, humidity transport, and condensation. Mass and momentum equations are recovered by an efficient streaming, collision, and forcing process within the framework of LBM while the governing equations of temperature, liquid, and vapor water fraction are solved using a finite volume method. The implementation of wall models for ABL, subgrid models, and interaction terms related to multiphysic phenomena (e.g., stratification, condensation) is described, implemented, and assessed in this study. An immersed boundary approach is used to handle flows in complex configurations, with application to flows in realistic urban areas. Applications to both wind engineering and atmospheric pollutant dispersion are illustrated.
Yongliang Feng, Johann Miranda‐fuentes, Shaolong Guo, Jérôme Jacob, Pierre Sagaut. ProLB: A Lattice Boltzmann Solver of Large‐Eddy Simulation for Atmospheric Boundary Layer Flows. Journal of Advances in Modeling Earth Systems, 2021, 13 (3), pp.e2020MS002107. ⟨10.1029/2020MS002107⟩. ⟨hal-03326123⟩
Journal: Journal of Advances in Modeling Earth Systems
Shang-Gui Cai, Johan Degrigny, Jean-François Boussuge, Pierre Sagaut. Coupling of turbulence wall models and immersed boundaries on Cartesian grids. Journal of Computational Physics, 2021, 429, pp.109995. ⟨10.1016/j.jcp.2020.109995⟩. ⟨hal-03326140⟩ Plus de détails...
An improved coupling of immersed boundary method and turbulence wall models on Cartesian grids is proposed, for producing smooth wall surface pressure and skin friction at high Reynolds numbers. Spurious oscillations are frequently observed on these quantities with most immersed boundary wall modeling methods, especially for the skin friction which is found to be very sensitive to the solid surface's position and orientation against the Cartesian grids. The problem originates from the irregularity of the wall distance on the stair-step grid boundaries where the immersed boundary conditions are applied. To reduce this directional error, several modifications are presented to enhance the near wall solution. First, the commonly used interpolation for the flow velocity is replaced by one for the friction velocity, which has much less variation near wall. The concept of using a fictitious point to retrieve flow fields in the wall normal direction is abandoned and the interpolation is performed in the wall parallel plane with existing fluid points. Secondly, the velocity gradients at the approximated boundary are computed with advanced schemes and the normal gradient of the tangential velocity is reconstructed from the wall laws. To further protect the near wall solution, the normal velocity gradient and the working viscosity from the Spalart-Allmaras turbulence model are enforced by their theoretical solutions in the interior fluid close to the wall. Additionally, various post-processing algorithms for reconstructing wall surface quantities and force integrations are investigated. Other related factors are also discussed for their effects on the results. The validity of present method has been demonstrated through numerical benchmark tests on a flat plate at zero pressure gradient, both aligned and inclined with respect to the grid, as well as aerodynamic cases of NACA 23012 airfoil and NASA trap wing.
Shang-Gui Cai, Johan Degrigny, Jean-François Boussuge, Pierre Sagaut. Coupling of turbulence wall models and immersed boundaries on Cartesian grids. Journal of Computational Physics, 2021, 429, pp.109995. ⟨10.1016/j.jcp.2020.109995⟩. ⟨hal-03326140⟩
Jincheng Lou, Jacob Johnston, Tzahi Cath, Denis Martinand, Nils Tilton. Computational fluid dynamics simulations of unsteady mixing in spacer-filled direct contact membrane distillation channels. Journal of Membrane Science, 2021, 622, pp.118931. ⟨10.1016/j.memsci.2020.118931⟩. ⟨hal-03515318⟩ Plus de détails...
Direct contact membrane distillation (DCMD) is a promising means of concentrating brines to their saturation limit. During that process, membrane spacers play a key role in temperature polarization, concentration polarization, and mineral scaling. These interactions are not well understood, because they are difficult to study experimentally and numerically, and the flow regimes are not fully charted. We consequently develop a tailored in-house CFD code that simulates unsteady two-dimensional heat and mass transport in plate-and-frame DCMD systems with cylindrical spacers. The code uses a combination of finite-volume methods in space, projection methods in time, and recent advances in immersed boundary methods for the spacer surfaces. Using the code, we explore how the transition to unsteady laminar vortex shedding affects polarization and permeate production of DCMD systems. We show that the impact of spacers can be explained by examining the various steady and unsteady vortical flow structures generated in the bulk and near the membranes. Overall, we show that though unsteady vortex structures tend to mix temperature polarization layers with the bulk, they are not similarly able to mix the concentration layers. Rather, vortical structures tend to create regions of preferential salt accumulation. In the vortex shedding regime, the net result is that spacers often increase vapor production at the expense of increasing the risk of mineral scaling.
Jincheng Lou, Jacob Johnston, Tzahi Cath, Denis Martinand, Nils Tilton. Computational fluid dynamics simulations of unsteady mixing in spacer-filled direct contact membrane distillation channels. Journal of Membrane Science, 2021, 622, pp.118931. ⟨10.1016/j.memsci.2020.118931⟩. ⟨hal-03515318⟩
Y. Feng, J. Miranda-Fuentes, Jérôme Jacob, Pierre Sagaut. Hybrid lattice Boltzmann model for atmospheric flows under anelastic approximation. Physics of Fluids, 2021, 33 (3), pp.036607. ⟨10.1063/5.0039516⟩. ⟨hal-03597258⟩ Plus de détails...
Lattice Boltzmann (LB) method for atmospheric dynamics is developed by considering the characteristics of the anelastic approximation. After introducing reference base state values in atmospheric flows, an LB model, with an external force term, has been constructed in anelastic framework. In the proposed anelastic LB model, mass and momentum conservation equations are solved by the LB method with a regularization procedure, and temperature field or scalar transport is simulated by finite volume method. The derived macroscopic governing equations from the anelastic model are analyzed and discussed in Chapman-Enskog asymptotic expansion. The anelastic LB model is assessed considering three benchmarks including a non-hydrostatic atmospheric inviscid convection, two-dimensional density currents, and inertia-gravity waves in stably stratified atmospheric layer. The validations demonstrate that the anelastic extension of the LB method can simulate atmospheric flows effectively and accurately. Besides, the proposed model offers a unified framework for both Boussinesq approximation and anelastic approximation, which is largely free of characteristic depth of atmospheric flows.
Y. Feng, J. Miranda-Fuentes, Jérôme Jacob, Pierre Sagaut. Hybrid lattice Boltzmann model for atmospheric flows under anelastic approximation. Physics of Fluids, 2021, 33 (3), pp.036607. ⟨10.1063/5.0039516⟩. ⟨hal-03597258⟩
Y. Feng, J. Miranda-Fuentes, Jérôme Jacob, Pierre Sagaut. Hybrid lattice Boltzmann model for atmospheric flows under anelastic approximation. Physics of Fluids, 2021, 33 (3), pp.036607. ⟨10.1063/5.0039516⟩. ⟨hal-03326143⟩ Plus de détails...
Lattice Boltzmann (LB) method for atmospheric dynamics is developed by considering the characteristics of the anelastic approximation. After introducing reference base state values in atmospheric flows, an LB model, with an external force term, has been constructed in anelastic framework. In the proposed anelastic LB model, mass and momentum conservation equations are solved by the LB method with a regularization procedure, and temperature field or scalar transport is simulated by finite volume method. The derived macroscopic governing equations from the anelastic model are analyzed and discussed in Chapman-Enskog asymptotic expansion. The anelastic LB model is assessed considering three benchmarks including a non-hydrostatic atmospheric inviscid convection, two-dimensional density currents, and inertia-gravity waves in stably stratified atmospheric layer. The validations demonstrate that the anelastic extension of the LB method can simulate atmospheric flows effectively and accurately. Besides, the proposed model offers a unified framework for both Boussinesq approximation and anelastic approximation, which is largely free of characteristic depth of atmospheric flows.
Y. Feng, J. Miranda-Fuentes, Jérôme Jacob, Pierre Sagaut. Hybrid lattice Boltzmann model for atmospheric flows under anelastic approximation. Physics of Fluids, 2021, 33 (3), pp.036607. ⟨10.1063/5.0039516⟩. ⟨hal-03326143⟩
Isabelle Cheylan, Song Zhao, Pierre Boivin, Pierre Sagaut. Compressible pressure-based Lattice-Boltzmann applied to humid air with phase change. Applied Thermal Engineering, 2021, pp.116868. ⟨10.1016/j.applthermaleng.2021.116868⟩. ⟨hal-03180596⟩ Plus de détails...
A new compressible pressure-based Lattice Boltzmann Method is proposed to simulate humid air flows with phase change. The variable density and compressible effects are fully resolved, effectively lifting the Boussinesq approximation commonly used, e.g. for meteorological flows. Previous studies indicate that the Boussinesq assumption can lead to errors up to 25%, but the model remains common, for compressible models often suffer from a lack of stability. In order to overcome this issue, a new pressure-based solver is proposed, exhibiting excellent stability properties. Mass and momentum conservation equations are solved by a hybrid recursive regularized Lattice-Boltzmann approach, whereas the enthalpy and species conservation equations are solved using a finite volume method. The solver is based on a pressure-based method coupled with a predictor-corrector algorithm, and incorporates a humid equation of state, as well as a specific boundary condition treatment for phase change. In particular, boundary conditions that handle mass leakage are also proposed and validated. Three test cases are investigated in order to validate this new approach: the Rayleigh-Bénard instability applied to humid air, the atmospheric rising of a condensing moist bubble, and finally the evaporation of a thin liquid film in a vertical channel. Results indicate that the proposed pressure-based Lattice-Boltzmann model is stable and accurate on all cases.
Isabelle Cheylan, Song Zhao, Pierre Boivin, Pierre Sagaut. Compressible pressure-based Lattice-Boltzmann applied to humid air with phase change. Applied Thermal Engineering, 2021, pp.116868. ⟨10.1016/j.applthermaleng.2021.116868⟩. ⟨hal-03180596⟩
