• Keine Ergebnisse gefunden

Analysis of a D1P3 lattice-Boltzmann equation

3.1 Energy estimate and stability

Multiplying equation (3.2) by ǫ2 and sending it to zero seems to avoid the trouble with infinite or undefined terms. However, under the hypothesis that fǫ behaves decently (e.g. staying bounded), the differential equation becomes in the limit

0 =J0f(t, x) =hf,1iw−f,

which is dramatically different from (3.2). In contrast to the initial value problem having a unique solution, this purely algebraic limit equation permits an infinite multitude of solutions. Any function of the formf(t, x) =w(t, x)w, wherew:XT → Ris someL-periodic function respecting the initial conditionw(0, x) =v0(x), solves the limit problem.

What may not be grasped at first glance is the fact that w is not at all arbitrary but the solution of (3.4). Clearly, this reveals the singular nature of the equation.

In order to disentangle the confusing situation, we will make the ansatz, that fǫ depends in a regular manner on ǫ. This means, that we try to approximate the solution by a truncated power series in ǫ. Amazingly, this procedure provides a lot of insight, although only a special class of initial values can be captured4.

3.1 Energy estimate and stability

Prior to the asymptotic consistency analysis, we establish a stability result. This states grosso modo, that the solution of the lattice-Boltzmann equation depends continuously on its initial value and a possible source term. The result is formulated in a way directly adapted to our purpose, guided by the following idea:

Consider an exact solution fǫ and an approximate solution ˆfǫ of the lattice-Boltz-mann equation. We can interpret the approximate solution as an exact solution with an additional source term corresponding to the residual. If stability holds, the two functions will stay close together during their evolution, provided they are initially almost equal and the residual is small. So asymptotic similarity can be inferred from stability, as soon as consistency is ensured, which means that the residual (i.e. the additional source driving ˆfǫ) vanishes withǫ.

A crucial role in the stability proof comes up to the following lemma which relies on the special structure of the equilibrium operator. In the sequel let W :F → F be the linear operator5 corresponding to a multiplication with the weight function w and letB =W1 its inverse.

Lemma 3.1. BJǫ is a symmetric (self-adjoint) operator with respect to the stan-dard scalar product in F. There exist orthogonal projectors Πk, k∈ {1,2,3}, with mutually disjoint ranges such that BJǫ permits the following representation:

BJǫ =−P3

k=1

λkΠk=−P3

k=1

λkΠkΠk (3.5)

4It turns out that the zeroth order term of the expansion must be proportional tow. Therefore the initialization (3.1) has been chosen so restrictive. Since fǫ(0, x) ∈ F, we have fǫ(0, x) span{w,s,1θ s2}. That is why the general initialization should be given byfǫ(0, x) = a(x)w+ b(x)s+c(x)`1

θ s2´

with some functionsa, b, c:X →R. Moreoverǫ-dependent initializations are imaginable, e.g. in the form of truncated power series.

5Iff= (f,f0,f+)thenW = diag(1,θ−1θ ,1).

Here the −λk’s (withλk≥0) denote the non-positiveeigenvalues.

Proof: SinceBJǫ =ω(BEǫ−B) and due to the symmetry ofB, the first part of the assertion is proved, if we verify, thatBEǫis symmetric too. Havingw=w·1=W·1, we find

BEǫ = (1−ǫ2τ c)h·,1i B W| {z }

=I

·1= (1−ǫ2τ c)h·,1i1.

So, BEǫ is just a multiple of the rank-one orthogonal projector h1,11 ih·,1i1, that is symmetric according to proposition 3.1. This implies that BEǫ has only a single non-zero eigenvalue given by 1−ǫ2τ c with the associated eigenvector is 1∈ F. By thespectral theoremforsymmetricoperators in Euclidean6vector spacesBJǫhas an orthonormal basis of eigenvectors. Let Πk fork∈ {1,2,3}denote theorthogonal projector7 on thek’th eigenvector. Since Πjkrefer to different eigenvectors ifj6= k we get ΠjΠk = 0. Using again proposition 3.1, the orthogonality of the projector Πk entails its symmetry, i.e. Πk = Πk, and thus Πk = Π2k = ΠkΠk. Designating with−λk the eigenvalue ofBJǫbelonging to Πk, the spectral decomposition ofBJǫ in (3.5) follows.

Let us finally check that all eigenvalues of BJǫ are non-positive. Considering first the case c= 0, we find

BJǫ =BJ0 =ωB(h·,1iw−I) =ω h·,1i1−B

1 1 1

1 1 1 1 1 1

· ·

· θ θ1 ·

· ·

Taking recourse to the explicit matrix representation of BJ0, it can be verified that BJ0 has the following non-positive eigenvalues:

0, −2ωθ, −ω(2θ+θθ1 −3) withθ >1.

HenceBJ0 is negative semi-definite. As the projector−ǫ2τ ch·,1i1 is negative semi-definite as well (recall τ >0 and c≥0) we conclude that BJǫ =BJ0−ǫ2τ ch·,1i1 is also negative semi-definite. Hence positive eigenvalues are excluded.

Proposition 3.1. Let P be a projector in a vector space V, i.e. P = P2. Then the following two statements are equivalent:

i) P is orthogonal, i.e. v, wV hP v,(IP)wi=h(IP)v, P wi= 0 ii) P is symmetric, i.e. P=P.

Proof:

i)ii): Expanding the assumption yields: hP v, wi=hP v, P wi=hv, P wi. This implies: 0 =hP v, wi − hv, P wi=h(PP)v, wi.

Since this holds for allv, wV, we conclude: PP= 0 P=P. ii)i): From the hypothesis we infer: 0 =PP2=PPP =P(IP).

Hence we find for allv, wV: 0 =hv, P(IP)wi=hP v,(IP)wi.

6By anEuclidean vector spacewe understand a finite-dimensional, real vector space equipped with a scalar product.

7Ifv6= 0 is a vector in an Euclidean vector space, then hv,vi1 , vivis theorthogonal projection onv.

3.1. Energy estimate and stability 131

Let us make the notion of an approximate solution more precise by recalling the definition ofconsistency. Anǫ-dependent function ˆfǫ ∈ Cper1 (XT,F) is said to be con-sistent of orderα to the lattice-Boltzmann equation, if there exists anǫ-dependent functionrǫ∈ Cper(XT,F) – called the residue function– such that:

i) ˆfǫ is an (exact) solution of the lattice-Boltzmann equation with the additional source termǫαrǫ:

tˆfǫ − ǫ1S∂xˆfǫ = ǫ2Jǫˆfǫ + gw + ǫαrǫ. (3.6) ii) rǫ is uniformly bounded with respect toǫ:

∃C >0, ∀ǫ∈(0,1] :

Now we can formulate the central stability result of this chapter:

Theorem 3.1. Let f ∈ Cper1 (XT;F) be a solution of the lattice-Boltzmann equation.

Suppose that ˆfǫ ∈ Cper1 (XT;F) is consistent to the lattice-Boltzmann equation of order α >0 and

kfǫ(0,·)−ˆfǫ(0,·)kL2(X;F)< K0ǫα (3.8) for some constant K0. Then there exists a constantK >0independent of ǫ∈(0,1]

such that

sup

t[0,T]kfǫ(t,·)−ˆfǫ(t,·)kL2(X;F)< Kǫα .

Proof: For convenience let us introduce thedeviationdǫ :=fǫ−ˆfǫ. As our equilib-rium operator is linear, the evolution equation fordǫis simply derived by subtracting (3.6) from (3.2).

tdǫ+1ǫS∂xdǫ= ǫ12 Jǫdǫ−ǫαrǫ (3.9) Moreover note, that B is a positive operator in F, whose operator norm satisfies 1≤min{2θ,θθ1} ≤ kBk ≤max{2θ,θθ1}ifθ >1 andkBk= 2 forθ= 1. Therefore the norm defined by ||| · |||2L2(X;F) :=R

XhB·,·i is equivalent to the standard norm in L2(X;F) and we have k·k2L2(X;F)≤ ||| · |||2L2(X;F) ≤ kBk k·k2L2(X;F).

Now the proof is done by an energy estimate of the deviation. This is obtained by taking the scalar product of equation (3.9) with Bdǫ and integrating over the spatial domain X. Let us consider the second integral on the left hand side, which is given in more detail by:

Since fǫ,ˆfǫ and consequently dǫ are L-periodic, i.e. dǫ(t,0) = dǫ(t, L), the integral vanishes. Remarking that hBdǫ, Jǫdǫi = hdǫ, BJǫdǫi, because B is diagonal as a multiplication operator and a fortiori symmetric, we can exploit lemma 3.1. So equation (3.10) becomes

The expression resulting from the integral over the collision term (the first one of the right hand side) is obviously negative. Thus we can estimate applying additionally the Cauchy-Schwarz inequality: for arbitrary reals a, b, ǫ. Finally, we apply Gronwall’s inequality8:

|||dǫ(t,·)|||2L2(X;F) ≤ exp(t)|||dǫ(0, .)|||2L2(X;F)exp(t)kBk Z t

0 krǫ(σ,·)k2L2(X;F)dσ By the initial condition (3.8) we have|||dǫ(0,·)|||L2(X;F)≤ kBkK0ǫα. From the defini-tion of consistency we know, that the integral is dominated by some constantC >0 for 0≤t≤T. So we end up with

Due to its linearity and the constancy of the parameters ω and c, its form will

8Gronwall’s lemma in differential form claims, that the non-negative function η ∈ C1([0, T]) satisfies the estimate

if the differential inequality dtdη(t) φ(t)η(t) +ψ(t) holds true with the nonnegative functions φ, ψ∈ L1([0, T]). For further reference see [17] page 624.