Chapter 6

Solenoidal Lipschitz truncation

Abstract

In this chapter we present the solenoidal Lipschitz truncation for non-stationary problems: we show how to construct a Lipschitz truncation which preserves the divergence-free character of a given Sobolev function. As a matter of fact, it suffices to have distributional time-derivatives in the sense of divergence-free test-functions. After this, we present the AImage-Stokes approximation for non-stationary problems. It aims at approximating almost solutions to the non-stationary AImage-Stokes system by exact solutions. Thanks to the solenoidal Lipschitz truncation this can be done on the level of gradients.

Keywords

Solenoidal Lipschitz truncation; Divergence-free constraint; Parabolic PDEs; Inverse curl-operator; AImage-Stokes approximation; Almost solutions

In this chapter we develop a non-stationary counterpart of the solenoidal Lipschitz truncation from Chapter 3. Here, the main difficulty is to handle problems connected with the distributional time derivative of the function we aim to truncate. Let us be a little bit more precise. Let Q0=I0×B0R×R3Image be a space time cylinder and σ(1,)Image. Let uLσ(I0,W1,σdiv(B0))Image and GLσ(Q0)Image satisfy

Q0tuξdxdt=Q0G:ξdxdtfor allξC0,div(Q0).

Image (6.0.1)

The main purpose of the solenoidal Lipschitz truncation is to avoid the appearance of the pressure function. Hence we start in (6.0.1) with an equation on the level of divergence-free test-functions. Unfortunately, this is not enough information on the time derivative for a Poincaré-type inequality as in (5.2.6). Hence the approach from [65] as explained in Section 5.2 will not give LImage-estimates for the gradient of the truncation, cf. the proof of Lemma 5.2.3. Our aim is to construct a truncation which preserves the properties from [65] and is, in addition, divergence-free.

We will show that there is a truncation uλImage of u with roughly the following properties (see Theorem 6.1.25 for a precise formulation).

(a) uλL(Q0)Image with uλcλImage and divuλ=0Image.

(b) uλ=uImage a.e. outside a suitable set OλImage.

(c) There holds

|tu,uλu|+χOαλuλppcλp|Oλ|δ(λ),

Image

with δ(λ)0Image if λImage.

In the following we sketch the construction on a heuristic level. In fact, the rigorous approach which we shall present in the next section requires a series of localization arguments, so it is quite technical. Let us start with a function

uL(I0;L2(B0))Lp(I0;W1,p0,div(B0))

Image

with tu=divHImage in Ddiv(B0)Image, where HLσ(B)Image for some σ>1Image. We define

w:=curl1uL(I0;W1,2(B0))Lp(I0;W2,pdiv(B0)).

Image

It follows that tΔw=curldivHImage in D(B0)Image. Also we can obtain an information about the time derivative of w as a distribution acting on all test-functions. However, we do not have control about a possible harmonic part of w. Hence we decompose w into a harmonic and anti-harmonic part. To do this we define, pointwise in time,

w(t)=z(t)+h(t),

Image

where z(t)ΔW2,p0(B0)Image and Δh(t)=0Image. This decomposition is based on a singular integral operator which is continuous on LpImage-spaces such that

z,wL(I0;W1,2(B0))Lp(I0;W2,p(B0)).

Image (6.0.2)

Moreover, we have

tΔz=tw=curldivH

Image

in D(B0)Image. As z is anti-harmonic by construction this yields

tzσcHσ.

Image

In fact, tzImage is a measurable function. Now, we truncate z to zλImage with an approach similar to (5.2.8). This truncation satisfies with 2zλcλImage as well as zλ=zImage in OλImage, where Oλ=Oλ(M(2z);M(tz))Image. Finally, we set

uλ:=curlzλ+curlh.

Image

Obviously, we have divuλ=0Image. Due to (6.0.2) and the properties of harmonic functions we have hL(I0;Wk,2(B0))Image for any kNImage (at least locally in space). Hence uλImage has the same regularity as curlzλImage. In particular, uλImage is bounded.

In Section 6.2 we develop the AImage-Stokes approximation for non-stationary problems, see [29]. This is, on the one hand, a non-stationary variant of the AImage-Stokes approximation from Section 3.3. On the other hand it is a fluid-mechanical counterpart of the AImage-caloric approximation from [68] which is concerned with the AImage-heat equation.

6.1 Solenoidal truncation – evolutionary case

In this section we examine solenoidal functions, whose time derivative is only a distribution acting on solenoidal test-functions. Let uLσ(I0,W1,σdiv(B0))Image be such that (6.0.1) holds for some GLσ(Q0)Image. So the time derivative is only well defined via the duality with solenoidal test functions. The goal of this section is to construct a solenoidal truncation uλImage of u which preserves the properties of the truncation in [65].

First we extend our function u in a suitable way to the whole space and then apply the inverse curl operator. Let γC0(B0)Image with χ12B0γχB0Image, where B0Image is a ball. Let C0Image denote the annulus B012B0Image. Then according to Theorem 2.1.6 (with A(t)=B(t)=tqImage) there exists a Bogovskiĭ operator BogC0:C0,(C0)C0(C0)Image which is bounded from Lq(C0)W1,q0(C0)Image for all q(1,)Image, and such that divBogC0=IdImage. Define

˜u:=γuBogC0(div(γu))=γuBogC0(γu).

Image

Then div˜u=0Image on I0×B0Image and ˜u(t)W1,σ0(B0)Image, so we can extend ˜uImage by zero in space to ˜uLσ(I0,W1,σdiv(R3))Image. Since ˜u=uImage on I0×12B0Image, we have

Q0t˜uξdxdt=Q0G:ξdxdtfor allξC0,div(12Q0).

Image (6.1.3)

Now, we define, pointwise in time,

w:=curl1(˜u)=curl1(γuBogC0(γu)).

Image

Overall, we get the following lemma.

Lemma 6.1.1

We have

curlw=˜u=uin12Q0divw=0inR3

Image

and

w(t)Ls(R3)cs˜u(t)La(B0)w(t)Ls(R3)cs˜u(t)Ls(B0)2w(t)Ls(R3)cs˜u(t)Ls(B0),

Image

for a=max{1,3s3+s}Image, tI0Image and s(1,)Image.

Let us derive from (6.0.1) the equation for w. For ψC0(12Q0)Image we have

Q0tucurlψdxdt=Q0G:curlψdxdt.

Image

We use u=curlwImage and partial integration to show that

Q0twcurlcurlψdxdt=Q0G:curlψdxdt.

Image

Now, because

Q0wtdivψdxdt=Q0divwtdivψdxdt=0

Image

and curlcurlψ=Δψ+divψImage we obtain

Q0wtΔψdxdt=Q0G:curlψdxdt

Image (6.1.4)

for every ψC0(12Q0)Image. We can rewrite this as

Q0wtΔψdxdt=Q0H:2ψdxdt,

Image (6.1.5)

with |G||H|Image pointwise. In particular, in the sense of distributions we have

tΔw=curldivG=divdivH.

Image

So in passing from u to w we got a system valid for all test functions ψC0(Q0)Image. However, we only have control of tΔwImage, so that the time derivative of the harmonic part of w cannot be seen. Hence, a parabolic Poincaré inequality for w still does not hold; i.e. twImage is not controlled! In order to remove this harmonic invariance we will replace w by some function z such that tΔw=tΔzImage. This will imply that tzImage can be controlled by H. To define z conveniently we need some auxiliary results.

For a ball BR3Image and a function fLs(B)Image we define Δ2BΔfImage as the weak solution FW2,s0(B)Image of

BΔFΔφdx=BfΔφdxfor allφC0(B).

Image (6.1.6)

Then fΔ(Δ2BΔf)Image is harmonic on BImage.

According to [117] and Lemma 2.1 of [140] we have the following variational estimate.

Lemma 6.1.2

Let s(1,)Image. Then for all gW2,s0(B)Image we have

2gscssupφC0(B)2φs1BΔgΔφdx.

Image

This implies the following two corollaries.

Corollary 6.1.1

Let s(1,)Image. Then

B|2(Δ2BΔf)|sdxcssupφC0(B)2φs1BfΔφdxcsB|f|sdx

Image

for fLs(B)Image, where csImage is independent of the ball BImage.

Proof

The claim follows by Lemma 6.1.2,

BΔ(Δ2BΔf)Δφdx=BfΔφdx,

Image

and Hölder's inequality.  □

Corollary 6.1.2

Let s(1,)Image. Then

23B|3(Δ2BΔf)|sdxcsB|f|sdxforfW1,s(B),23B|4(Δ2BΔf)|sdxcsB|2f|sdxforfW2,s(B),

Image

where csImage is independent of the ball BImage.

Proof

The claim follows from Corollary 6.1.1 by standard interior regularity theory (difference quotients, localization and Poincaré's inequality).  □

For VLs(B)Image we define Δ2BdivdivVImage as the weak solution FW2,s0(B)Image of

BΔFΔφdx=BV2φdxfor allφC0(B).

Image

Similar to Corollary 6.1.1 we get the following result.

Corollary 6.1.3

Let s(1,)Image. Then

B|2(Δ2BdivdivV)|sdxcsB|V|sdx

Image

for VLs(B)Image, where csImage is independent of the ball BImage.

The next lemma shows the wanted control of the time derivative.

Lemma 6.1.3

For a cube Q=I×BQ0Image let zQ:=ΔΔ2BΔwImage. Then for s(1,)Image we have

Q|zQ|sdxcsQ|w|sdxQ|zQr|sdx+23Q|zQ|sdxcsQ|w|sdxQ|zQ(r)2|sdx+23Q|zQr|sdx+23Q|2zQ|sdxcsQ|2w|sdx,Q|tzQ|sdxdtcsQ|H|sdxdt,

Image

where r:=rBImage.

Proof

The estimate of zQImage in terms of w follows directly from Corollary 6.1.1 and integration over time. The estimate of zQImage in terms of wImage and 2wImage follows from this by Poincaré's inequality, using the fact that we can subtract a linear polynomial from w without changing the definition of zQImage. The other estimate for zQImage and 2zQImage follow analogously from Corollary 6.1.2.

For all ρC0(I)Image and φC0(B)Image it follows from (6.1.5) that

IBwΔφdxtρdt=IBH:2φdxρdt.

Image

Let dhtImage denote the forward difference quotient in time with step size h. We use ρ(t):=tht˜ρ(τ)dτImage with ˜ρC0(I)Image and h sufficiently small. Then tρ=dht˜ρImage and

IBwΔφdxdht˜ρdt=IBH:2φdxtht˜ρ(τ)dτdt.

Image

This implies that

IBdhtwΔφdx˜ρdt=IBt+htH(τ)dτ:2φdx˜ρdt.

Image

Since this is valid for all choices of ˜ρImage we have

BdhtwΔφdx=Bt+htH(τ)dτ:2φdx.

Image

Since dhtzQ=dht(ΔΔ2BΔw)=ΔΔ2BΔ(dhtw)Image, it follows by Corollary 6.1.1 that

(B|dhtzQ|sdx)1scsupφC0(B)2φs1BdhtwΔφdx=supφC0(B)2φs1(Bt+htH(τ)dτ:2φdx)c(Bt+ht|H(τ)|sdτdx)1s.

Image

Integrating over time and passing to the limit h0Image yields

(Q|tzQ|sdxdt)1sc(Q|H|sdxdt)1s

Image

which finishes the proof.  □

Defining z(t):=z12Q0(t)=ΔΔ212B0Δw(t)Image for t12I0Image, we then have

Q0ztΔψdxdt=Q0wtΔψdxdt=Q0H:2ψdxdt,

Image (6.1.7)

for all ψC0(12Q0)Image. Since the function Δ212B0w(t)W2,s0(12B0)Image, we can extend it by zero to a function in W2,s(R3)Image. In this sense it is natural to extend z(t)Image by zero to a function in Ls(R3)Image.

Note that Lemma 6.1.3 enables us to control tzImage by H in Ls(12Q0)Image.

Lemma 6.1.4

We have

z(t)Ls(13B0)cs˜u(t)L3ss+3(B0)z(t)Ls(13B0)cs˜u(t)Ls(B0)2z(t)Ls(13B0)cs˜u(t)Ls(B0),

Image

for tIImage and s(1,)Image.

Proof

This follows from Corollary 6.1.1, Corollary 6.1.2 and Lemma 6.1.1.  □

For λ,α>0Image and σ>1Image we define

Oαλ:={Mασ(χ13Q0|2z|)>λ}{αMασ(χ13Q0|tz|)>λ}.

Image (6.1.8)

Later we will choose α=λ2pImage and σ smaller than the integrability exponent of tzImage. We want to redefine z on OαλImage. The first step is to cover OαλImage by well selected cubes. By the lower-semi-continuity property of the maximal functions the set OαλImage is open. We assume in the following that OαλImage is non-empty. (In the case that OαλImage is empty, we do not need to truncate at all.) We cover OαλImage by an α-parabolic Whitney covering {Qi}Image with partition of unity in accordance with Lemmas 5.2.1 and 5.2.2.

Due to property (PW3) we have that 16Qj(Rd+1Oαλ)Image. Thus, the definition of OαλImage implies that

(16Qj|2z|σχ13Q0dxdt)1σλ,

Image (6.1.9)

α(16Qj|tz|σχ13Q0dxdt)1σλ.

Image (6.1.10)

Lemma 6.1.5

Assume that there exists c0>0Image such that λp|Oαλ|c0Image with p>2dd+2Image. Then the following holds:

if λλ0=λ0(c0)Image, α=λ2pImage and Qi14Q0Image, then Qi13Q0Image and Qj13Q0Image for all jAiImage.

Proof

Let Qi14Q0Image. We claim that Qi724Q013Q0Image. Let si:=αr2iImage. It suffices to show that ri,si0Image as λImage. Because QiOαλImage and by assumption, we find

λ2rd+2i=λpsirdicλp|Qi|cλp|Oαλ|cc0.

Image (6.1.11)

This implies that ri(cc0λ2)1d+20Image as λImage. Moreover, ri=s12iα12Image and (6.1.11) imply

cc0λpsirdi=λpsd+22iα32=λd+22pdsd+22i.

Image

If p>2dd+2Image, then λImage implies si0Image as desired.

The claim on jAiImage follows from the fact that QiImage and QjImage have comparable size and that 724Q0Image is strictly contained in 13Q0Image.  □

Let us show that the assumption λp|Oαλ|c0Image from Lemma 6.1.5 is satisfied in our situation. To do this we assume from now on that

α:=λ2p

Image (6.1.12)

and that σ<min{p,p}Image.

Lemma 6.1.6

Let c0:=2zpLp(13Q0)+tzpLp(13Q0)Image. Then

λp|Oαλ|c0.

Image

Proof

It follows from the weak-type estimate of MασImage in (5.2.5), if σ<min{p,p}Image, then

|Oαλ|cλp2zpLp(13Q0)+c(λα1)ptzpLp(13Q0)=cλp(2zpLp(13Q0)+tzpLp(13Q0)).

Image

 □

In the following we choose λ0Image such that the conclusion of Lemma 6.1.5 is valid and assume λλ0Image. Without loss of generality we can assume further that

λ0(13Q0|2z|σdxdt)1σ+r20(13Q0|z|σdxdt)1σ.

Image (6.1.13)

We define

I:={i:Qi14Q0}.

Image

Then Lemma 6.1.5 implies that Qi13Q0Image (and Qj13Q0Image for jAiImage) for all iIImage. For each iIImage we define local approximation ziImage for z on QiImage by

zi:=Π0IiΠ1Bi(z),

Image (6.1.14)

where Π1Bi(z)Image is the first order averaged Taylor polynomial [37,63] with respect to space and Π0IiImage is the zero order averaged Taylor polynomial in time. Note that this definition implies the Poincaré-type inequality.

Lemma 6.1.7

For all jNImage and 1s<Image, if 2z,tzLs(14Q0)Image, then

Qj|zzjr2j|sdxdt+Qj|(zzj)rj|sdxdtcQj|2z|sdxdt+cαsQj|tz|sdxdt.

Image

Proof

The estimate is a consequence of Fubini's Theorem, Poincaré estimates and the properties of the averaged Taylor polynomials see Lemma 3.1 of [63]. We find

Qj|zzjr2j|sdxdtcQj|zΠ1Bj(z)r2j|sdxdt+cBjIj|Π1Bj(z)Π0IjΠ1Bj(z)r2j|sdxdtcQj|2z|sdxdt+cαIjBj|tΠ1Bj(z)|sdxdt.

Image

Now the continuity of Π1BjImage on LsImage gives the estimate. Similarly we find (since all norms for polynomials are equivalent)

Qj|(zzj)rj|sdxdtcQj|(zΠ1Bj(z))rj|sdxdt+cQj|(Π1Bj(z)Π0IjΠ1Bj(z))rj|sdxdtcQj|(zΠ1Bj(z))rj|sdxdt+cQj|Π1Bj(z)Π0IjΠ1Bj(z)r2j|sdxdtcQj|2z|sdxdt+cαQj|tz|sdxdt.

Image

 □

We can now define our truncation zαλImage for λλ0Image on 14Q0Image by

zαλ:=ziIφi(zzi).

Image (6.1.15)

It suffices to sum over i with Qi14Q0Image.

Since the φiImage are locally finite, this sum is pointwise well-defined. We will see later that the sum converges also in other topologies. Using iIφi=1Image on 14Q0Image, we can also write zαλImage in the form

zαλ={zin14Q0Oαλ,iIφiziin14Q0Oαλ.

Image (6.1.16)

In the following we describe some properties of the truncation (e.g. 2zαλL(14Q0)Image).

Lemma 6.1.8

For all jNImage and all kNImage with QjQkImage we have

(a) Qj|2z|dxdt+αQj|tz|dxdtcλImage.

(b) zjzkL(Qj)cQj|zzj|dxdt+cQk|zzk|dxdtImage.

(c) zjzkL(Qj)cr2jλImage.

Proof

Part (a) follows from Qj16QjImage and 16QjOλImage, so

16Qj(|2z|+α|tz|)χ13Q0dxdt(16Qj(|2z|+α|tz|)σχ13Q0dxdt)1σcλ.

Image

Part (b) follows from the geometric property of the QjImage. If QjQkImage, then |QjQk|cmax{|Qj|,|Qk|}Image. This and the norm equivalence for linear polynomials imply

zjzkL(Qj)cQjQk|zjzk|dxdtcQj|zjz|dx+cQk|zzk|dx.

Image

Finally, (c) is a consequence of Lemma 6.1.7, (a) and (b).  □

Next, we prove the stability of the truncation.

Lemma 6.1.9

Let 1<s<Image and zLs(R;W2,s(R3))Image. Then

zαλLs(14Q0)czLs(13Q0),zαλLs(14Q0)czLs(13Q0)+cαr0tzLs(13Q0),2zαλLs(14Q0)+αtzαλLs(13Q0)c2zLs(13Q0)+cαtzLs(13Q0).

Image

Moreover, the sum in (6.1.15) converges in Ls(14I0,W2,s(14B0))Image.

Proof

We first show that the sum in (6.1.15) converges absolutely in Ls(14Q0)Image:

14Q0|zzαλ|sdxiIQi|zzi|sdxdtciIQi|z|sdxdtc13Q0|z|sdt,

Image

where we used continuity of the mapping zziImage in Ls(Qi)Image, (PP1) and the finite intersection property of QiImage (PW4). We start by showing the estimate for the second derivatives

Oαλ|2(zzαλ)|sdxdt=|iIQi2(φi(zzi))dxdt|ciIQi|2z|s+|(zzi)ri|s+|zzir2i|sdxdt.

Image

For the time derivative we find (since ziImage is constant in time), that

Oαλ|t(zzαλ)|sdxdt=|iIQit(φi(zzi))dxdt|ciIQi|tz|s+|zziαr2i|sdxdt.

Image (6.1.17)

Using Lemma 6.1.7 and the finite intersection of the QiImage shows that

14Q0|2(zzαλ)|s+αs|t(zzαλ)|sdxdtciIQi|2z|s+αs|tz|sdxdtc13Q0|2z|s+αs|tz|sdxdt.

Image

The estimate of the gradient is analogous, since

Oαλ|(zzαλ)|sdxdtiI|(zzi)|s+|zziri|sdxdt.

Image

 □

The truncation zαλImage has better regularity properties than z. Indeed, zImage is Lipschitz.

Lemma 6.1.10

For λ>λ0Image we have

2zαλL(14Q0)+r10zαλL(14Q0)+r20zαλL(14Q0)+αtzαλL(14Q0)cλ.

Image

Proof

If (t,x)QiImage, then

|2zαλ(t,x)|=|jAi2(φjzj)(t,x)|jAi|2(φj(zjzi))(t,x)|

Image

because {φj}Image is a partition of unity. Now we find (since all norms on polynomials are equivalent, #AjcImage and Lemma 6.1.8) that

|2zαλ(t,x)|cjAizizjL(Qi)r2icλ.

Image

Concerning the time derivative for (t,x)QiImage since ziImage is constant in time we find that

|tzαλ(t,x)|=|tjAi(φjzj)(t,x)|jAi|t(φj)(zjzi)(t,x)|jAizizjL(Qi)αr2icλα.

Image

The zero order term is estimated by Poincaré's inequality; first in time and then in space

r20zαλL(14I0;L(14B0))cαtzαλL(14Q0)+cr20zαλL1(14I0;L(14B0))cλ+c2zαλL(14Q0)+cr20zαλL1(14Q0).

Image

This implies, by the norm equivalence of polynomials, Jensen's inequality, Lemma 6.1.9 and (6.1.13),

r10zαλL(14Q0)+r20zαλL(14Q0)cλ+r20zLσ(13Q0)cλ.

Image

 □

The next lemma will control the time error we get when we apply the truncation as a test function.

Lemma 6.1.11

For all ζC0(14Q0)Image with 2ζcImage and λλ0Image,

|14Q0t(zzαλ)Δ(ζzαλ)dxdt|cα1λ2|Oαλ|,

Image

where the constant c is independent of α and λ.

Proof

We use Hölder's inequality and Lemma 6.1.10 to derive

(I):=|14Q0t(zzαλ)Δ(ζzαλ)dxdt|iI(Qi|t(φi(zzi)|σdxdt)1σ(Qi|Δ(ζzαλ)|σdxdt)1σcλiI|Qi|(Qi|t(φi(zzi)|σdxdt)1σ.

Image

Combining this with (6.1.17), (6.1.9) and (6.1.10) yields

(I)cλiI|Qi|(α1(Qj|2z|σdxdt)1σ+(Qj|tz|σdxdt)1α)cα1λ2iI|Qi|cα1λ2|Oαλ|,

Image

using the local finiteness of the {Qi}Image in the final step.  □

Theorem 6.1.24

Let 1<p<Image with p,p>σImage. Let wmImage and HmImage satisfy tΔwm=divdivHmImage in the sense of distributions D(12Q0)Image, see (6.1.5). Further assume that wmImage is a weak null sequence in Lp(12I0;W2,p(12B0))Image and a strong null sequence in Lσ(12Q0)Image. Further, assume that Hm=H1m+H2mImage such that H1mImage is a weak null sequence in Lp(Q0)Image and H2mImage converges strongly to zero in Lσ(Q0)Image. Define zm:=ΔΔ212B0ΔwmImage pointwise in time on 12I0Image. Then there is a double sequence (λm,k)R+Image and k0NImage such that

(a) 22kλm,k22k+1Image

such that the double sequence zm,k:=zαm,kλm,kImage with αm,k:=λ2pm,kImage satisfies the following properties for all kk0Image

(b) {zm,kz}Om,k:=Oαm,kλm,kImage,

(c) 2zm,kL(14Q0)cλm,kImage,

(d) zm,k0Image and zm,k0Image in L(14Q0)Image for mImage and k fixed,

(e) 2zm,k0Image in L(14Q0)Image for mImage and k fixed,

(f) We have for all ζC0(14Q0)Image

|(t(zmzm,k))Δ(ζzm,k)dxdt|cλpm,k|Om,k|,

Image

(g) limsupmλpm,k|Om,k|c2ksupm(2zmp+cH1m1p1p)Image.

Proof

Let us assume that λm,kImage satisfies (a). We will choose the precise values of λm,kImage later. Due to Lemma 6.1.3 we have zm0Image in Lp(14I0;W2,p(14B0))Image; this is due to the fact that the operator wΔΔ212B0Δw=zImage is linear and continuous in Lp(14I0;W2,p(14B0))Image. Then the properties (b) and (c) follow from Lemma 6.1.10. Moreover, Corollary 6.1.1 ensures that the strong convergence in Lσ(12Q0)Image transfers from wmImage to zmImage. By Lemma 6.1.9 we get the same for zm,kImage and that the sequence 2zm,kImage is, for fixed k and s, bounded in Ls(14Q0)Image. The combination of these convergence properties implies (by interpolation) (d). Moreover, the boundedness of 2zm,kImage in Ls(14Q0)Image implies the weak convergence of a subsequence. Since (d) ensures that the limit is zero, we get, by the usual arguments, weak convergence of the whole sequence. This proves (e). Moreover, (f) follows by Lemma 6.1.11 and the choice of αm,kImage.

It remains to choose 22kλm,k22k+1Image such that (g) holds. We use the decomposition

tzm=ΔΔ212B0divdivHm=ΔΔ212B0divdivH1m+ΔΔ212B0divdivH2m=:h1m+h2m.

Image

We decompose

Om,k={Mαm,kσ(χ13Q0|2zm,k|)>λm,k}{αm,kMαm,kσ(χ13Q0|tzm|)>λm,k}{Mαm,kσ(χ13Q0|2zm,k|)>λm,k}{αm,kMαm,kσ(χ13Q0|h1m|)>12λm,k}{αm,kMαm,kσ(χ13Q0|h2m|)>12λm,k}=:IIIIII.

Image

Define

gm:=2Mαm,kσ(χ13Q0|2zm|)+(2Mαm,kσ(χ13Q0|h1m|))1p1.

Image

Then by the boundedness of MσImage on LpImage and LpImage (using p,p>σImage), as well as Corollary 6.1.3, we have

gmp2Mαm,kσ(χ13Q0|2zm|)p+(2Mαm,kσ(χ13Q0|h1m|))1p1p=2Mαm,kσ(χ13Q0|2zm|)p+2Mαm,kσ(χ13Q0|h1m|)1p1pc2zmLp(13Q0)+ch1m1p1Lp(12Q0)c2zmLp(13Q0)+cH1m1p1Lp(12Q0).

Image

Let K:=supm(2zmp+ch1m1p1p)Image. In particular, gmpKImage uniformly in k. Note that

III={Mαm,kσ(χ13Q0|2zm,k|)>λm,k}{(Mαm,kσ(χ13Q0|h1m|))1p1>λm,k}{2Mαm,kσ(χ13Q0|2zm,k|)+(2Mαm,kσ(χ13Q0|h1m|))1p1>λm,k}={gm>λm,k}.

Image

We estimate

Rd+1|gm|pdx=Rd+101ptp1χ{|gm|>t}dtdxRd+1kN1p2kχ{|gm|>2k+1}dxjN2j+11k=2j1p2kp|{|gm|>2k+1}|.

Image

For fixed m,jImage the sum over k involves 2jImage summands and not all of them can be large. Consequently there exists λm,k{22k+1,,22k+1}Image, such that

λpm,k|{|gm|>λm,k}|c2kKp

Image

uniformly in m and k, and hence

λpm,k|III|λpm,k|{gm>λm,k}|c2kKp.

Image (6.1.18)

On the other hand, from the weak-LσImage estimate for Mαm,kσImage we see that

limsupm(λpm,k|III|)=limsupm(λpm,k|{αm,kMαm,kσ(χ13Q0|h2m|)>12λm,k}|)limsupm(cλpm,kh2mσLσ(13Q0)(αm,k/λm,k)σ).

Image

Since 22k+1λm,k22k+1Image, αm,k=λ2pm,kImage and h2m0Image in Lσ(12Q0)Image (which is a consequence of H2m0Image in Lσ(12Q0)Image and Corollary 6.1.3), it follows that

limsupm(λpm,k|III|)=0.

Image

This and (6.1.18) prove (g).  □

Theorem 6.1.25

Let 1<p<Image with p,p>σImage. Let umImage and GmImage satisfy tum=divGmImage in the sense of distributions Ddiv(Q0)Image. Assume that umImage is a weak null sequence in Lp(I0;W1,p(B0))Image and a strong null sequence in Lσ(Q0)Image and bounded in L(I0,T;Lσ(B0))Image. Further assume that Gm=G1m+G2mImage such that G1mImage is a weak null sequence in Lp(Q0)Image and G2mImage converges strongly to zero in Lσ(Q0)Image. Then there is a double sequence (λm,k)R+Image and k0NImage with

(a) 22kλm,k22k+1Image

such that the double sequences um,k:=uαm,kλm,kL1(Q0)Image, αm,k:=λ2pm,kImage and Om,k:=Oαm,kλm,kImage (defined in Theorem 6.1.24) satisfy the following properties for all kk0Image.

(b) um,kLs(14I0;W1,s0,div(16B0))Image for all s<Image and supp(um,k)16Q0Image.

(c) um,k=umImage a.e. on 18Q0Om,kImage.

(d) um,kL(14Q0)cλm,kImage.

(e) um,k0Image in L(14Q0)Image for mImage and k fixed.

(f) um,k0Image in L(14Q0)Image for mImage and k fixed.

(g) limsupmλpm,k|Om,k|c2k.Image

(h) limsupm|Gm:um,kdxdt|cλpm,k|Om,k|Image.

Proof

We define, pointwise in time on I0Image,

˜um:=γumBogB012B0(γum),wm:=curl1˜um,zm:=ΔΔ212Q0Δwm,

Image

where γC0(Q0)Image with χ12Q0γχQ0Image. Then we apply Theorem 6.1.24 to the sequence zmImage. Finally, let

um,k:=curl(ζzm,k)+curl(ζ(wmzm)),

Image (6.1.19)

where ζC0(16Q0)Image with χ18Q0ζχ16Q0Image. This means on 18Q0Image we have

um,k=um+curl(zm,kzm).

Image

Note that curl(wmzm)Image is harmonic (in space) on 12Q0Image and bounded in time, due to the assumption that umImage is bounded uniformly in L(I0;Lσ(B0))Image, which transfers to wmImage and zmImage by Lemma 6.1.1 and 6.1.4. This allows us to estimate the higher order spaces derivatives on 14Q0Image by lower order ones on 12Q0Image. This, (6.1.19) and Theorem 6.1.24 immediately imply all the claimed properties except (h).

The claim of (g) follows exactly as (g) of Theorem 6.1.24.

Let us prove (h). It follows by simple density arguments that um,kImage is an admissible test function for the equation tum=divGmImage. We thus obtain

Gm:um,kdxdt=tumum,kdxdt=(tcurlwm)curl(ζzm,k)dxdt+(tcurlwm)curl(ζ(wmzm))dxdt=(tzm)Δ(ζzm,k)dxdt(tzm)Δ(ζ(wmzm))dxdt=:T1+T2.

Image

Here we took into account curlcurlwm=ΔwmImage (due to divwm=0Image) and Δwm=ΔzmImage. By assumption GmImage is bounded in Lσ(Q0)Image. Using regularity properties of harmonic functions (for wmzmImage) as well as Lemma 6.1.3 and Lemma 6.1.1 we gain (after choosing a subsequence)

(Q0|Δ(ζ(wmzm))|σdxdt)1σcr20(14Q0|wmzm|3σ3+σdxdt)3+σ3σcr20(12Q0|wm|3σ3+σdxdt)3+σ3σcr30(Q0|˜um|σdxdt)1σ0asm.

Image

Since, additionally, tzmImage is uniformly bounded in Lσ(12Q0)Image by Lemma 6.1.3, we have T20Image as mImage. Furthermore, there holds

T1=(t(zmzm,k))Δ(ζzm,k)dxdt+(tzm,k)Δ(ζzm,k)dxdt=:T1,1+T1,2,

Image

where the first term can be bounded using Theorem 6.1.24 (f). So it remains to show that

T1,2:=(tzm,k)Δ(ζzm,k)dxdt0asm.

Image

We have

T1,2=12t(|zm,k|2)ζdxdt+(tzm,k)div(ζzm,k)dxdt=12|zm,k|2tζdxdt+(tzm,k)div(ζzm,k)dxdt.

Image

The first term is estimated by Theorem 6.1.24(d). For the second we use Lemma 6.1.9 and Lemma 6.1.3 (s=σImage) to find

|tzm,k||div(ζzm,k)|dxdtc(13Q0|Gm|σ+|2zm|σdxdt)1σ(13Q0|zm,k|σ+|zm,k|σdxdt)1σ.

Image

Now because GmImage and 2zmImage are uniformly bounded in Lσ(12Q0)Image we find by Theorem 6.1.24 (d), that

limmT1,2=0,

Image

which proves the claim of (h).  □

The following corollary is useful in the application of the solenoidal Lipschitz truncation.

Corollary 6.1.4

Let all assumptions of Theorem 6.1.25 be satisfied with ζC0(16Q0)Image with χ18Q0ζχ16Q0Image as in the proof of Theorem 6.1.25. If additionally umImage is uniformly bounded in L(I0,Lσ(B0))Image, then for every KLp(16Q0)Image

limsupm|((G1m+K):um)ζχOm,kdxdt|c2k/p.

Image

Proof

It follows from (f), (g) and (h) of Theorem 6.1.25 that

limsupm|(Gm+K):um,kdxdt|cλpm,k|Om,k|c2k.

Image (6.1.20)

Recall that um,k=curl(ζzm,k)+curl(ζ(wmzm))Image. So, by Theorem 6.1.24, we have zm,k,zm,k0Image in L(14Q0)Image as mImage with k fixed. Since umImage is a strong null sequence in Lσ(Q0)Image and is bounded in L(I0,Lσ(B0))Image we see that um0Image strongly in Ls(I0,Lσ(B0))Image for any s(1,)Image. By continuity of the Bogovskiĭ operator (see Theorem 2.1.6 with A(t)=B(t)=tσImage) we have the same convergence for ˜umImage. Now, Lemma 6.1.1 implies wm=curl1˜um0Image in Ls(I0,W1,σ(R3))Image. Using zm:=ΔΔ212Q0ΔwmImage and Corollary 6.1.1 we also have zm0Image in Ls(I0,W1,σ(R3))Image. Since zmwmImage is harmonic on 14Q0Image, we have zmwm0Image in Ls(I0,W2,s(16B0))Image. These convergences imply that

um,k=ζcurlzm,k+am,k,

Image

with am,k0Image in Ls(16Q0)Image as mImage with k fixed. This, the boundedness of GmImage in Lσ(16Q0)Image, KLp(16Q0)Image and (6.1.20) imply (using s>σImage)

limsupm|((Gm+K):curlzm,k)ζdxdt|c2k.

Image

Since Gm=G1m+G2mImage, G2m0Image in Lσ(16Q0)Image and zm,k0Image in Lσ(16Qz)Image for mImage and k fixed, we have

limsupm|((G1m+K):curlzm,k)ζdxdt|c2k.

Image (6.1.21)

The boundedness of G1mImage and K in Lp(16Q0)Image and Theorem 6.1.24 and (g) prove

limsupm|((G1m+K):curlzm,k)ζχOm,kdxdt|c2k/p.

Image

This, together with (6.1.21) and zm,k=zmImage in Om,kImage yield

limsupm|((G1m+K):curlzm)ζχOm,kdxdt|c2k/p.

Image

Recall that zmwm0Image in Ls(I0,W2,s(16B0))Image for any s(1,)Image. This and the boundedness of G1mImage in Lp(Q0)Image allows us to replace zmImage in the previous integral by wmImage. Now curlwm=umImage proves the claim.  □

The next corollary follows by combining Lemma 6.1.9, Lemma 6.1.10, Theorem 6.1.24 (g) (with α=1Image) and the continuity of curl1Image with a scaling procedure.

Corollary 6.1.5

For some σ>0Image let uLσ(I0;W1,σdiv(B0))L(I;Lσ(B0))Image with tu=divHImage in Ddiv(Q0)Image for some HLσ(Q0)Image. Then for every m01Image and γ>0Image there exist λ[2m0γ,22m0γ]Image and a function uλImage with the following properties.

(a) It holds uλL(I0,W1,0,div(B0))Image with uλcλImage.

(b) We have

λσLd+1(12Q0{uλu})|Q0|cm0(Q0rσ0|u|σ+|u|σdxdt+Q0|H|σdxdt).

Image

(c) It holds

Q0|uλ|σdxdtc(Q0|u|σ+Q0rσ0|H|σdxdt),Q0|uλ|σdxdtc(Q0rσ0|u|σ+|u|σdxdt+Q0|H|σdxdt).

Image

(d) We have t(uuλ)Lσ(12I0,W1,σ(12B0))Image and

12Q0(uuλ)tφdxdtc(κ)12Q0χ{uλu}|φ|σdxdt+κ(Q0rσ0|u|σ+|u|σ+|H|σdxdt)

Image

for all φC0(12Q0)Image and all κ>0Image.

Proof

We apply the arguments used in the proof of Theorem 6.1.25 to the constant sequence u with the choice α=1Image. So we set

uλ:=curl(ζzλ)+curl(ζ(wz)),

Image (6.1.22)

where ζC0(16Q0)Image with χ18Q0ζχ16Q0Image. Hence, in 18Q0Image we have

uλ=u+curl(zλz).

Image

We immediately obtain the claim of (a). As a consequence of the Lemmas 6.1.1, 6.1.4 and 6.1.9 we obtain the inequalities

Q|uλ|σdxdtc(Q0|u|σ+Qrσ0|tz|σdxdt),Q0|uλ|σdxdtc(Q0rσ0|u|σ+|u|σdxdt+Q|tz|σdxdt),

Image

claimed in c). Finally we can replace tzImage by H on account of Lemma 6.1.3.

It remains to find good levels. We define for some s(1,max{σ,σ})Image

Oλ:={Ms(χ13Q0|2z|)>λ}{Ms(χ13Q0|H|)>λ},g:=Ms(χ13Q0|2z|)+Ms(χ13Q0|H|).

Image

It now follows from the continuity of MsImage in (5.2.4), together with Lemmas 6.1.1 and 6.1.4

Rd+1|g|σdxc(Q0rσ0|u|σ+|u|σdxdt+Q0|H|σdxdt).

Image (6.1.23)

Furthermore, the following holds for every m0NImage and every γ>0Image

Rd+1|g|σdx=Rd+101σtσ1χ{|g|>t}dtdxRd+12m01m=m01σ(2mγ)σχ{|g|>γ2m+1}dx.

Image

So, there is m1{m0,...,2m01}Image such that

Rd+1(2m1γ)σχ{|g|>γ2m1+1}dxcm0Rd+1|g|σdx.

Image

Setting λ=γ2m1+1Image yields

λσ|13Q0{|g|>λ}|cm0Rd+1|g|σdx.

Image

Combining this with (6.1.23) gives the estimate in b) due to the definition of OλImage.

Finally, we prove d). We have uλu=curl(zλz)Image in 18Q0Image such that Lemma 6.1.3 and 6.1.9 imply t(uλu)Lσ(18I0,W1,σ(18B0))Image. Moreover, we have for φC0(18Q0)Image

18Q0(uuλ)tφdxdt=18Q0χOλt(zzλ)curlφdxdtc(κ)12Q0χOλ|φ|σdxdt+κ(Q0|t(zλz)|σdxdt),

Image

as a consequence of Young's inequality. Applying Lemmas 6.1.3, 6.1.4 and 6.1.9 yields

Q0|t(zλz)|σdxdtc(Q0rσ0|u|σ+|u|σdxdt+Q0|H|σdxdt).

Image

So we have shown the estimate claimed in (d) on 18Q0Image.

A simple scaling argument allows us the obtain all the estimates on 12Q0Image.  □

Remark 6.1.10

The higher dimensional case

For general dimensions, the solenoidal Lipschitz truncation is best understood in terms of differential forms. We start with ˜uImage as given in (6.1.3). Now, we have to find w such that curlw=˜uImage and divw=0Image. Let us define the 1-form α in RdImage associated to the vector field ˜uImage by α:=i˜uidxiImage. Then we need to find a 2-form G such that dG=αImage and dG=0Image, where d is the outer derivative and dImage its adjoint by the scalar product for k-forms. Similar to w=curl1˜u=curlΔ1˜uImage we get G by G:=dΔ1αImage. Since we are on the whole space, Δ1Image can be constructed by mollification with c|x|2dImage. Thus, we have

G(x)=(dΔ1α)(x)=di(Rdui(y)|xy|d2dy)dxi.

Image

Let us explain how to substitute the equation tΔw=curldivGImage, see (6.1.4). Instead of test functions ψ with divψ=0Image we use the associated 1-forms β=iψidxiImage with dβ=0Image. Thus there exists a 2-form γ with dγ=0Image. Then

t˜u,ψ=tα,β=tdG,dγ=tddG,γ=tΔG,γ,

Image

where we used Δ=dd+ddImage and dα=0Image in the last step. Note that −Δ applied to the form G is the same as −Δ applied to the vector field of all components of G. Now we define w as the associated vector field (with (d2)Image components) of G and we arrive again at an equation for tΔwImage. This concludes the construction; the rest can be done exactly as for dimension three. The restriction p>65Image used in this section will change to 2dd+2Image.

6.2 AImage-Stokes approximation – evolutionary case

By AImage we denote a symmetric, elliptic tensor, i.e.

c0|τ|2A(τ,τ)c1|τ|2for allτRd×d.

Image (6.2.24)

We set |A|:=c1/c0Image. Let BRdImage be a ball and J=(t0,t1)Image a bounded interval. We set Q=J×BImage. For a function wL1(Q)Image with twLq(J;W1,q(B))Image we introduce the unique function HwLq0(Q)Image with

Qwtφdxdt=QHw:φdxdt

Image

for all φC0,div(Q)Image. We begin with a variational inequality for the non-stationary AImage-Stokes system.

Lemma 6.2.1

Suppose that (6.2.24) holds and that q>1Image. There holds for every uCw([t0,t1];L1(B))Lq(J;W1,q(B))Image with u(t0,)=0Image a.e.

Q|u|qdxdtcsupξC0,div(Q)[Q(A(ε(u),ε(ξ))utξ)dxdtQ(|ξ|q+|Hξ|q)dxdt],

Image

where c only depends on AImage, q and d.

Proof

Duality arguments show that

1qQ|u|qdxdt=supGLq(Q)[Qu:Gdxdt1qQ|G|qdxdt].

Image

For a given GLq(Q)Image let zGImage be the unique Lq(J;W1,q0,div(B))Image-solution to

Qztξdxdt+QA(ε(z),ε(ξ))dxdt=QG:ξdxdt

Image (6.2.25)

for all ξC0,div((t0,t1]×B)Image. This is a backward parabolic equation with end datum zero. We have that tzGLq(J;W1,qdiv(B))Image, so test-functions can be chosen from the space Lq(J;W1,q0,div(B))Image. Due to Theorem B.3.50 (which can be applied to ˜z˜G(t,)=zG(t1t,)Image, where ˜G(t,)=G(t1t,)Image) this solution satisfies

Q|zG|qdxdt+Q|HzG|qdxdtcQ|G|qdxdt.

Image

In other words, the mapping Lq(B)GzGLq(J;W1,q0,div(B))Image is continuous. This and u(t0,)=0Image yield (using u as a test-function in (6.2.25))

Q|u|qdxdtcsupGLq(Q)[QA(ε(u),ε(zG))dxdtQtzGudxdtQ(|zH|q+|HzG|q)dxdt]csupξC0,div(Q)[QA(ε(u),ε(ξ))dxdtQutξdxdtQ(|ξ|q+|Hξ|q)dxdt],

Image

which yields the claim.  □

Let us now state the AImage-Stokes approximation. In the following let BImage be a ball with radius r and J an interval with length 2r2Image. Let ˜QImage denote either Q=J×BImage or 2Q. We use similar notations for ˜JImage and ˜BImage.

Theorem 6.2.26

Suppose that (6.2.24) holds. Let vLqs(2˜J;W1,qsdiv(2˜B))Image, q,s>1Image, be an almost AImage-Stokes solution in the sense that

|2Qvtξdxdt2QA(ε(v),ε(ξ))dxdt|δ2˜Q|ε(v)|dxdtξ

Image (6.2.26)

for all ξC0,div(2Q)Image and some small δ>0Image. Then the unique solution wLq(J;W1,q0,div(B))Image to

QwtξdxdtQA(ε(w),ε(ξ))dxdt=QvtξdxdtQA(ε(v),ε(ξ))dxdt

Image (6.2.27)

for all ξC0,div([t0,t1)×B)Image satisfies

Q|wr|qdxdt+Q|w|qdxdtκ(2˜Q|v|qsdxdt)1s.

Image

It holds κ=κ(q,s,δ)Image and limδ0κ(q,s,δ)=0Image. The function h:=vwImage is called the AImage-Stokes approximation of v.

Remark 6.2.11

From the proof of Theorem 6.2.26 we have the following stability result, choosing p=qs=qImage.

Q|wr|pdxdt+Q|w|pdxdtc2˜Q|v|pdxdt.

Image

Indeed κ stays bounded if s1Image.

Proof

Let w be defined as in (6.2.27). Combining Poincaré's inequality with Lemma 6.2.1 and (6.2.27) shows that

Q|wr|qdxdt+Q|w|qdxdtcsupξC0,div(Q)[QA(ε(v),ε(ξ))dxdtQvtξdxdtQ(|ξ|q+|Hξ|q)dxdt].

Image (6.2.28)

In the following let us fix ξC0,div(Q)Image. Let

γ:=(Q|ξ|qdxdt+Q|Hξ|qdxdt)1q,

Image

and m0NImage, m01Image. Due to Corollary 6.1.5, applied with σ=qImage, we find λ[2m0γ,22m0γ]Image and ξλL(4J;W1,0,div(4B))Image such that

ξλL(4Q)cλ,

Image (6.2.29)

λqLd+1(2Q{ξλξ})|Q|cm0(Q|ξ|qdxdt+Q|Hξ|qdxdt),

Image (6.2.30)

4Q|ξλ|qdxdtc(Q|ξ|qdxdt+Qrq|Hξ|qdxdt),

Image (6.2.31)

4Q|ξλ|qdxdtc(Q|ξ|qdxdt+Q|Hξ|qdxdt).

Image (6.2.32)

Note that ξ can be extended by 0 to 4Q thus the equation

tξ=divBogB(tξ)=:divHξ

Image

holds on 4Q by the properties of BogBImage (since HξImage can be extended as well). For the properties of the Bogovskiĭoperator we refer to Section 2.1, in particular Theorem 2.1.6. Corollary 6.1.5 (d) implies that t(ξξλ)Lq(2J,W1,q(2B))Image and

2Jt(ξξλ),φdtc(κ)2Qχ{ξξλ}|φ|qdxdt+κ(Q|ξ|q+|Hξ|qdxdt),

Image (6.2.33)

for all φW1,q0(2Q)Image. For ηC0(2Q)Image with η1Image on Q, |kη|crkImage and |tk1η|cr(k+1)Image (k=1,2Image) we see that

QA(ε(v),ε(ξ))dxdtQvtξdxdt=2d+22QA(ε(v),ε(ηξBog2BB(ηξ)))dxdt2Qvt(ηξ...)dxdt=2d+2(2QA(ε(v),ε(ηξλBog2BB(ηξλ)))dxdt+2Qtv(ηξλ...)dxdt)+2d+22QA(ε(v),ε(η(ξξλ)Bog2BB(η(ξξλ))))dxdt+2d+22Qtv(η(ξξλ)Bog2BB(η(ξξλ)))dxdt=:2d+2(I+II+III).

Image

Note that the time-derivative of v exists in the W1,divImage-sense as a consequence of (6.2.26), so all terms are well-defined by the properties of ξλImage. We have the following inequality on account of the continuity properties of ∇Bog on LpImage-spaces, (6.2.31), (6.2.32) and Poincaré's inequality (we set ˜ξλ:=ξξλImage):

2Q|Ψλ|qdxdt:=2Q|(η˜ξλ)Bog2BB(η˜ξλ)|qdxdtc2Q|˜ξλ|qdxdt+c2Q|˜ξλr|qdxdtcQ|ξ|qdxdt+cQ|ξr|qdxdt+cQ|Hξ|qdxdtcQ|ξ|qdxdt+cQ|Hξ|qdxdt.

Image (6.2.34)

Young's inequality, with an appropriate choice of ε>0Image, together with (6.2.31) and (6.2.32), implies that

IIc(ε)2Q|ε(v)|qχ{ξξλ}dxdt+ε2Q|Ψλ|qdxdtc2Q|ε(v)|qχ{ξξλ}dxdt+13Q|ξ|q+|Hξ|qdxdt=:II1+II2,

Image

where c depends on AImage, q and qImage. Hölder's inequality now yields

II1c(2Q|v|qsdxdt)1s(Ld+1(2Q{ξλξ})|Q|)11s.

Image

If follows from (6.2.30), by the choice of γ and λγImage that

Ld+1(2Q{ξλξ})|Q|cγqm0λqcm0.

Image (6.2.35)

Thus

II1c(2Q|v|qsdxdt)1s(cm0)11s.

Image

We choose m0Image sufficiently large that

II1κ3(2Q|v|qsdxdt)1s.

Image

Since t(ξξλ)Lq(2J,W1,q(2B))Image we can write III as

III=2Qvtη(ξξλ)dxdt+2Qηvt(ξξλ)dxdt2QvBog2BB(tη(ξξλ))dxdt2QBog2BB(v)ηt(ξξλ)dxdt=:III1+III2+III3+III4.

Image

The Bogovskiĭ  operator is continuous from L2L2Image. Hence its dual (in the sense of L2Image-duality) is continuous from L2L2Image. Therefore Bog2BB(v)Image is well-defined. We consider the four terms separately. For the first one we have

III1c2I2BBχ{ξλξ}|vr||ξξλr|dxdtc(ε)2I2BB|vr|qχ{ξλξ}dxdt+ε2Q|ξξλr|qdxdt=:c(ε)III11+εIII12.

Image

Poincaré's inequality and Young's inequality yield

III11c(2I2BB|vr|qsdxdt)1s(Ld+1(2Q{ξλξ})|Q|)11sc(2Q|v|qsdxdt)1s(Ld+1(2Q{ξλξ})|Q|)11s.

Image

Arguing as for the term II1Image implies that

III11κ12(2Q|v|qsdxdt)1s.

Image

Moreover, we gain from (6.2.31) and Poincaré's inequality

III12cQ|ξr|qdxdt+cQ|Hξ|qdxdtcQ|ξ|qdxdt+cQ|Hξ|qdxdt,

Image

and finally

III1κ12(2Q|v|qsdxdt)1s+112(Q|ξ|q+|Hξ|qdxdt).

Image

The formulation in (6.2.26) does not change if we subtract terms which are constant in space from v (note that (tξ)=0Image for every t due to tξ(t,)C0,div(B)Image). So we can assume that

2BBv(t)dx=0for a.e.t2J.

Image (6.2.36)

As a consequence of (6.2.33), (6.2.36) and Poincaré's inequality we obtain similarly as for III1Image

III2c(ε)2Qχ{ξξλ}|(ηv)|qdxdt+ε(Q|ξ|q+|Hξ|qdxdt)κ12(2Q|v|qsdxdt)1s+112(Q|ξ|q+|Hξ|qdxdt).

Image

Taking into account continuity properties of the Bogovskiĭ  operator from LqW1,q0Image we can estimate III3Image via (we use again (6.2.36) and Poincaré's inequality)

III3(2Q2BB|vr|qsdxdt)1qs×(2Qr(qs)|Bog2BB(tη(ξξλ)))|(qs)dxdt)1(qs)c(2Q|v|qsdxdt)1qs(2Qr2(qs)|tη(ξξλ)|(qs)dxdt)1(qs)c(2Q|v|qsdxdt)1qs(2Qχ{ξξλ}|ξξλr|(qs)dxdt)1(qs).

Image

Hence, from Young's inequality for every ε>0Image, we deduce that

III3εq12(2Q|v|qsdxdt)1s+cεq(2Qχ{ξξλ}|ξξλr|(qs)dxdt)q(qs)=:εq12III31+cεqIII32.

Image

It follows due to Hölder's inequality, Poincaré's inequality, (6.2.30), (6.2.32) and (6.2.35) for m0Image large enough

III32(Ld+1(2Q{ξλξ})|Q|)11s(2Q|ξξλr|qdxdt)c(Ld+1(2Q{ξλξ})|Q|)11s(4Q|ξξλ4r|qdxdt)c(Ld+1(2Q{ξλξ})|Q|)11s(4Q(|ξ|q+|ξλ|q)dxdt)κ12c(Q|ξ|q+|Hξ|qdxdt).

Image

Choosing ε:=κ1/qImage implies

III3κ12(2Q|v|qsdxdt)1s+112(Q|ξ|q+|Hξ|qdxdt).

Image

By (6.2.33) and (6.2.35) we have for m0Image large enough

III4c2Qχ{ξξλ}|(ηBog2BB(v))|qdxdt+112(Q|ξ|q+|Hξ|qdxdt)ε(2Q|(ηBog2BB(v))|sqdxdt)1s+112(Q|ξ|q+|Hξ|qdxdt)=:εIII41+112III42.

Image

Due to the continuity of Bog(div())Image on LpImage for any 1<p<Image (see [85, III.3, Theorem 3.3] and Theorem 2.1.7 for the Bogovskiĭ  operator and negative norms) we have continuity of BogImage as well. This, Poincaré's inequality (note that Bog2BB(v)Lp0(2BB)Image) and (6.2.36) yield

III41c(2Q|Bog2BB(v)r2|sqdxdt+2Q|Bog2BB(v)r|sqdxdt)1sc(2Q|Bog2BB(v)r|sqdxdt)1sc(2I2BB|vr|sqdxdt)1sc(2Q|v|sqdxdt)1s,

Image

and hence for ε:=κ/12cImage,

III4κ12(2Q|v|qsdxdt)1s+112(Q|ξ|q+|Hξ|qdxdt).

Image

Combining the estimates for III1ImageIII4Image, we see that

IIIκ3(2Q|v|qsdxdt)1s+13(Q|ξ|q+|Hξ|qdxdt).

Image

Since v is an almost AImage-Stokes solution and ξλcλc2m0γImage we have

|I|δ2˜Q|v|dxdtξλ,2Qδ(2˜Q|v|qsdxdt)1qsc2m0γ.

Image

We apply Young's inequality and Jensen's inequality to give

|I|δ2m0c(2˜Q|v|qsdxdt)1s+δ2m0cγqδ2m0c(2˜Q|v|qsdxdt)1s+δ2m0c(Q|ξ|qdxdt+Q|Hξ|qdxdt).

Image

Now, we choose δ>0Image so small such that δ2m0cκ/3Image. Thus

|I|κ3(2˜Q|v|qsdxdt)1s+13(Q|ξ|qdxdt+Q|Hξ|qdxdt).

Image

Combining the estimates for I, II and III we have established

2QA(ε(v),ε(ξ))dxdtQvtξdxdtκ(2˜Q|v|qsdxdt)1s+Q|ξ|qdxdt+Q|Hξ|qdxdt.

Image

Inserting this in (6.2.28) shows the claim.  □

References

[29] D. Breit, Existence theory for generalized Newtonian fluids, In: Recent Advances in Partial Differential Equations and Applications. Contemp. Math.. Providence, RI: Amer. Math. Soc.; 2016;vol. 666:99–110.

[37] S.C. Brenner, L.R. Scott, The Mathematical Theory of Finite Element Methods. Texts Appl. Math.. New York: Springer-Verlag; 1994;vol. 15.

[63] L. Diening, M. Růžička, Interpolation operators in Orlicz Sobolev spaces, Numer. Math. 2007;107(1):107–129.

[65] L. Diening, M. Růžička, J. Wolf, Existence of weak solutions for unsteady motions of generalized Newtonian fluids, Ann. Sc. Norm. Super. Pisa, Cl. Sci. (5) 2010;IX:1–46.

[68] F. Duzaar, G. Mingione, Second order parabolic systems, optimal regularity, and singular sets of solutions, Ann. Inst. Henri Poincaré, Anal. Non Linéaire 2005;22:705–751.

[85] G.P. Galdi, An introduction to the mathematical theory of the Navier–Stokes equations, vol. I. Springer Tracts Nat. Philos.. Berlin–New York: Springer; 1994;vol. 38.

[117] R. Müller, Das schwache Dirichletproblem in LqImage für den Bipotentialoperator in beschränkten Gebieten und in Außengebieten. [PhD thesis] Universität Bayreuth, Bayreuth / Bayreuth. Math. Schr. 49; 1994/1995:115–211.

[140] J. Wolf, Existence of weak solutions to the equations of nonstationary motion of non-Newtonian fluids with shear-dependent viscosity, J. Math. Fluid Mech. 2007;9:104–138.

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
18.218.151.44