Analysis of Lattice-Boltzmann Methods

Contents

1 Introduction to lattice Boltzmann methods and their analysis

1.1 Initiation to lattice-Boltzmann methods

1.1.1 A brief introduction to kinetic theory

1.1.2 A primer of lattice-Boltzmann methods

1.2 Translational invariance and dimensional reduction

1.3 An abstract framework of numerical analysis

2 Scalings and singular limits on the basis of the D1P2 model

2.1 Hyperbolic versus parabolic scaling

2.2 A singularly perturbed initial value problem

2.2.1 The Fourier coefficient functions

2.2.2 Solution of the perturbed problem and convergence

2.2.3 Uniform convergence and convergence rate

2.3 Two-scale expansion and resolution of the initial layer

3 Analysis of a D1P3 lattice-Boltzmann equation

3.1 Energy estimate and stability

3.2 Regular expansion and consistency

3.3 Smoothness conditions and convergence

3.4 Initial conditions and irregular expansions

3.5 A glimpse of boundary conditions

4 Consistency of a D1P3 lattice-Boltzmann algorithm

4.1 Formal expansion

4.2 Consistency and asymptotic similarity

4.3 Construction of consistent population functions

4.4 Initial behavior

5 Long-term behavior of an advective lattice-Boltzmann scheme

5.1 Regular expansion

5.1.1 Analysis of the update rule

5.1.2 Smooth initialization and consistency

5.2 Multiscale Expansion

5.2.1 A numeric test to detect different time scales

5.2.2 Additional quadratic time scale

5.2.3 Emergence of a cubic time scale

6 Stability investigations around the D1P2 model

6.1 Basics concerning shift matrices

6.2 LB advection-diffusion scheme with periodic boundary conditions

6.2.1 An ℓ∞-stability result

6.2.2 The spectral limit set of the evolution matrices

6.2.3 Asymptotics and symmetry of eigenvalues

6.3 LB advection scheme with periodic boundary conditions

6.3.1 The CFL-condition and stability

6.3.2 Stability in the ℓ2-norm

6.3.3 Multiscale expansion and stability

6.4 LB diffusion scheme with bounce-back type boundary conditions

6.4.1 Evolution matrices and their spectra

6.4.2 Computing eigenbases

6.5 Towards the D1P3 scheme & Concluding remarks

7 Asymptotic analysis of a numeric boundary layer

7.1 Some remarks about interpolation and difference stencils

7.2 Model problem: 1D Poisson equation with Dirichlet BC

7.3 Discretization of Dirichlet boundary conditions

7.4 Stability of extrapolation schemes

7.5 Damping property of discrete inverse operators

7.6 Asymptotic expansions and convergence

7.7 Numeric experiments

Research Objectives and Topics

This thesis aims to provide a rigorous mathematical understanding of Lattice-Boltzmann methods, which are numerical schemes for solving evolutionary partial differential equations (PDEs). The central research question explores how Lattice-Boltzmann algorithms approximate PDEs and how numerical phenomena such as initial layers, boundary layers, and multiple time scales affect their stability and consistency.

• Rigorous analysis of Lattice-Boltzmann methods using regular and irregular asymptotic expansions (Hilbert expansion).
• Investigation of singular limits in scaled finite-velocity Boltzmann equations.
• Stability analysis of various algorithms, including spectral analysis and the role of the CFL-condition.
• Examination of numerical boundary layers, particularly for the one-dimensional Poisson equation.
• Bridging the gap between mesoscopic kinetic theory and macroscopic fluid dynamics.

Auszug aus dem Buch

Summary of Chapters

Chapter 1: Provides fundamental background on kinetic theory, the Boltzmann equation, and an introduction to the Lattice-Boltzmann method's mathematical framework.

Chapter 2: Analyzes the D1P2 Lattice-Boltzmann model, focusing on the distinction between hyperbolic and parabolic scaling and the investigation of initial layers.

Chapter 3: Presents an analysis of a D1P3 Lattice-Boltzmann equation, utilizing energy estimates and asymptotic expansions to establish convergence and stability.

Chapter 4: Conducts a formal consistency analysis of a D1P3 Lattice-Boltzmann algorithm discretizing the advection-diffusion equation, highlighting the recursive hierarchy of evolution equations.

Chapter 5: Discusses the long-term behavior of an advective Lattice-Boltzmann scheme, identifying different time scales using multiscale expansions.

Chapter 6: Investigates the stability of D1P2 models through spectral analysis of evolution matrices and examines the role of boundary conditions.

Chapter 7: Performs an asymptotic analysis of numerical boundary layers using the Poisson equation as a model problem, emphasizing the damping properties of discrete inverse operators.

Keywords

Lattice-Boltzmann methods, kinetic theory, asymptotic expansion, numerical stability, partial differential equations, initial layers, boundary layers, convergence analysis, singular limits, Hilbert expansion, Chapman-Enskog expansion, finite difference schemes, spectral analysis, grid refinement, computational fluid dynamics.

Frequently Asked Questions

What is the fundamental focus of this research?

The work focuses on the mathematical analysis of Lattice-Boltzmann algorithms to understand their consistency, convergence, and stability when solving evolutionary partial differential equations.

Which scientific methods are primarily employed?

The research relies heavily on asymptotic analysis, specifically using regular and irregular (Hilbert) expansions, multiscale expansions, and spectral analysis of discrete evolution operators.

What are the central thematic areas covered in this document?

The document covers kinetic theory, the derivation of Lattice-Boltzmann methods from Boltzmann equations, singular limits, the impact of initial and boundary layers, and the stability of these numerical schemes.

What is the primary research goal?

The main objective is to provide a rigorous mathematical foundation for Lattice-Boltzmann methods, moving beyond purely plausibility-based arguments to improve the reliability of these simulations in practical applications.

How is the consistency of the schemes analyzed?

Consistency is established by utilizing regular asymptotic expansions to derive how the discrete Lattice-Boltzmann equations approximate the macroscopic target equations, such as the Navier-Stokes or advection-diffusion equations.

Which specific "numerical phenomena" are identified and analyzed?

The thesis investigates numerical phenomena such as initial layers, boundary layers, and the presence of multiple time scales that can interfere with the accuracy and stability of the algorithms.

How do different velocity models (D1P2, D1P3, D2P9) relate to each other?

The thesis demonstrates that higher-dimensional or more complex models, such as the D2P9, can be reduced to simpler one-dimensional models like D1P3 or D1P2 under certain symmetry and initialization conditions.

What is the significance of the "damping property" mentioned in the later chapters?

The damping property is a crucial condition required for analyzing the stability of discrete inverse operators, particularly when dealing with singular residuals in boundary layer analysis.