Smooth initialization of lattice Boltzmann schemes

Authors:
Pieter Van Leemput;Martin Rheinländer;Michael Junk
Affiliations:
Department of Computer Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium;Fachbereich Mathematik und Statistik, Universität Konstanz, Postfach D194, 78457 Konstanz, Germany;Fachbereich Mathematik und Statistik, Universität Konstanz, Postfach D194, 78457 Konstanz, Germany
Venue:
Computers & Mathematics with Applications
Year:
2009

Citing 5
Cited 1

Scientific Computing with Ordinary Differential Equations

Scientific Computing with Ordinary Differential Equations
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems
Constraint-Defined Manifolds: a Legacy Code Approach to Low-Dimensional Computation

Journal of Scientific Computing
DNS and LES of decaying isotropic turbulence with and without frame rotation using lattice Boltzmann method

Journal of Computational Physics
An Efficient Newton-Krylov Implementation of the Constrained Runs Scheme for Initializing on a Slow Manifold

Journal of Scientific Computing

An Efficient Newton-Krylov Implementation of the Constrained Runs Scheme for Initializing on a Slow Manifold

Journal of Scientific Computing

Quantified Score

Hi-index	0.09

Visualization

Abstract

Lattice Boltzmann methods are paradigmatic discrete evolutions with incomplete initial conditions. This is due to the fact that the variables of the (mesoscopic) method outnumber the variables of the (macroscopic) problem to be solved. In such situations, most initializations which are compatible with the given macroscopic data lead to solutions with oscillatory or steep initial layers. In order to reduce such initial effects, we present a general approach to construct initial values which are compatible with the partial information available, and which guarantee a smooth start of the evolution. Our smoothness condition prescribes the unknown initial values as the polynomial backward extrapolation of the values obtained from a few time steps. Specifically for constant and linear extrapolation, we study the consistency, stability and accuracy of the approach in the case of a lattice Boltzmann method for one-dimensional advection. Moreover, the applicability of a simple iteration scheme as the solution method is investigated.