Researcher profile

Pierre Sagaut

Pierre Sagaut contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2021arXiv

Hydrodynamic limits and numerical errors of isothermal lattice Boltzmann schemes

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 pre-factor. It is at the origin of over-dissipation and severe instabilities in the low-viscosity regime.

preprint2020arXiv

A linear stability analysis of compressible hybrid lattice Boltzmann methods

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, formulation (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 analysed. 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_{ref}$ and (2) introducing a controllable numerical dissipation $σ$ 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.

preprint2020arXiv

Analysis and reduction of spurious noise generated at grid refinement interfaces with the lattice Boltzmann method

The present study focuses on the unphysical effects induced by the use of non-uniform grids in the lattice Boltzmann method. In particular, the convection of vortical structures across a grid refinement interface is likely to generate spurious noise that may impact the whole computation domain. This issue becomes critical in the case of aeroacoustic simulations, where accurate pressure estimations are of paramount importance. The purpose of this article is to identify the issues occurring at the interface and to propose possible solutions yielding significant improvements for aeroacoustic simulations. More specifically, this study highlights the critical involvement of non-physical modes in the generation of spurious vorticity and acoustics. The identification of these modes is made possible thanks to linear stability analyses performed in the fluid core, and non-hydrodynamic sensors specifically developed to systematically emphasize them during a simulation. Investigations seeking pure acoustic waves and sheared flows allow for isolating the contribution of each mode. An important result is that spurious wave generation is intrinsically due to the change in the grid resolution (i.e. aliasing) independently of the details of the grid transition algorithm. Finally, the solution proposed to minimize spurious wave amplitude consists of choosing an appropriate collision model in the fluid core so as to cancel the non-hydrodynamic mode contribution regardless the grid coupling algorithm. Results are validated on a convected vortex and on a turbulent flow around a cylinder where a huge reduction of both spurious noise and vorticity are obtained.

preprint2020arXiv

Lattice Boltzmann method for computational aeroacoustics on non-uniform meshes: a direct grid coupling approach

The present study proposes a highly accurate lattice Boltzmann direct coupling cell-vertex 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~\textit{et al. 2020}) on the deleterious effect of non-hydrodynamic modes crossing mesh transitions, which can be addressed using an appropriate collision model. The Hybrid Recursive Regularized model is then used for this study. The grid coupling algorithm is assessed and compared to a widely-used cell-vertex algorithm on an acoustic pulse test case, a convected vortex and a turbulent circular cylinder wake flow at high Reynolds number.

preprint2020arXiv

Lattice Boltzmann Method for wave propagation in elastic solids with a regular lattice: Theoretical analysis and validation

The von Neumann stability analysis along with a Chapman-Enskog analysis is proposed for a single-relaxation-time lattice Boltzmann Method (LBM) for wave propagation in isotropic linear elastic solids, using a regular D2Q9 lattice. Different boundary conditions are considered: periodic, free surface, rigid interface. An original absorbing layer model is proposed to prevent spurious wave reflection at domain boundaries. The present method is assessed considering several test cases. First, a spatial Gaussian force modulated in time by a Ricker wavelet is used as a source. Comparisons are made with results obtained using a classical Fourier spectral method. Both P and S waves are shown to be very accurately predicted. The case of Rayleigh surface waves is then addressed to check the accuracy of the method.

preprint2020arXiv

Structural stability of Lattice Boltzmann schemes

The goal of this work is to determine classes of traveling solitary wave solutions for Lattice Boltzmann schemes by means of an hyperbolic ansatz. It is shown that spurious solitary waves can occur in finite-difference solutions of nonlinear wave equation. The occurence of such a spurious solitary wave, which exhibits a very long life time, results in a non-vanishing numerical error for arbitrary time in unbounded numerical domain. Such a behavior is referred here to have a structural instability of the scheme, since the space of solutions spanned by the numerical scheme encompasses types of solutions (solitary waves in the present case) that are not solutions of the original continuous equations. This paper extends our previous work about classical schemes to Lattice Boltzmann schemes.

preprint2013arXiv

A stochastic view of isotropic turbulence decay

A stochastic EDQNM approach is used to investigate self-similar decaying isotropic turbulence at high Reynolds number ($400 \leq Re_λ\leq 10^4$). The realistic energy spectrum functional form recently proposed by Meyers & Meneveau is generalised by considering some of the model constants as random parameters, since they escape measure in most experimental set-ups. The induced uncertainty on the solution is investigated building response surfaces for decay power-law exponents of usual physical quantities. Large-scale uncertainties are considered, the emphasis being put on Saffman and Batchelor turbulence. The sensitivity of the solution to initial spectrum uncertainties is quantified through probability density functions of the decay exponents. It is observed that initial spectrum shape at very large scales governs the long-time evolution, even at high Reynolds number, a parameter which is not explicitly taken into account in many theoretical works. Therefore, a universal asymptotic behavior in which kinetic energy decays as $t^{-1}$ is not detected. But this decay law is observed at finite Reynolds number with low probability for some initial conditions.

preprint2013arXiv

Quantification of errors in large-eddy simulations of a spatially-evolving mixing layer

A stochastic approach based on generalized Polynomial Chaos (gPC) is used to quantify the error in Large-Eddy Simulation (LES) of a spatially-evolving mixing layer flow and its sensitivity to different simulation parameters, viz. the grid stretching in the streamwise and lateral directions and the subgrid scale model constant ($C_S$). The error is evaluated with respect to the results of a highly resolved LES (HRLES) and for different quantities of interest, namely the mean streamwise velocity, the momentum thickness and the shear stress. A typical feature of the considered spatially evolving flow is the progressive transition from a laminar regime, highly dependent on the inlet conditions, to a fully-developed turbulent one. Therefore the computational domain is divided in two different zones (\textit{inlet dependent} and \textit{fully turbulent}) and the gPC error analysis is carried out for these two zones separately. An optimization of the parameters is also carried out for both these zones. For all the considered quantities, the results point out that the error is mainly governed by the value of the $C_S$ constant. At the end of the inlet-dependent zone a strong coupling between the normal stretching ratio and the $C_S$ value is observed. The error sensitivity to the parameter values is significantly larger in the inlet-dependent upstream region; however, low error values can be obtained in this region for all the considered physical quantities by an ad-hoc tuning of the parameters. Conversely, in the turbulent regime the error is globally lower and less sensitive to the parameter variations, but it is more difficult to find a set of parameter values leading to optimal results for all the analyzed physical quantities.

preprint2012arXiv

Analysis of the absorbing layers for the weakly-compressible lattice Boltzmann schemes

It has been demonstrated that Lattice Boltzmann schemes (LBSs) are very efficient for Computational AeroAcoustics (CAA). In order to handle the issue of absorbing acoustic boundary conditions for LBS, three kinds of damping terms are proposed and added into the right hand sides of the governing equations of LBS. From the classical theory, these terms play an important role to absorb and minimize the acoustic wave reflections from computational boundaries. Meanwhile, the corresponding macroscopic equations with the damping terms are recovered for analyzing the macroscopic behaviors of the these damping terms and determining the critical absorbing strength. Further, in order to detect the dissipation and dispersion behaviors, the linearized LBS with the damping terms is derived and analyzed. The dispersive and dissipative properties are explored in the wave-number spaces via the Von Neumann analysis. The related damping strength critical values and the optimal absorbing term are addressed. Finally, some benchmark problems are implemented to assess the theoretical results.

preprint2011arXiv

Optimal low-dispersion low-dissipation LBM schemes for computational aeroacoustics

Lattice Boltmzmann Methods (LBM) have been proved to be very effective methods for computational aeroacoustics (CAA), which have been used to capture the dynamics of weak acoustic fluctuations. In this paper, we propose a strategy to reduce the dispersive and disspative errors of the two-dimensional (2D) multi-relaxation-time lattice Boltzmann method (MRT-LBM). By presenting an effective algorithm, we obtain a uniform form of the linearized Navier-Stokes equations corresponding to the MRT-LBM in wave-number space. Using the matrix perturbation theory and the equivalent modified equation approach for finite difference methods, we propose a class of minimization problems to optimize the free-parameters in the MRT-LBM. We obtain this way a dispersion-relation-preserving LBM (DRP-LBM) to circumvent the minimized dispersion error of the MRT-LBM. The dissipation relation precision is also improved.And the stability of the MRT-LBM with the small bulk viscosity is guaranteed. Von Neuman analysis of the linearized MRT-LBM is performed to validate the optimized dispersion/dissipation relations considering monochromatic wave solutions. Meanwhile, dispersion and dissipation errors of the optimized MRT-LBM are quantitatively compared with the original MRT-LBM . Finally, some numerical simulations are carried out to assess the new optimized MRT-LBM schemes.