Appendix B

The AImage-Stokes system

The aim of this section is to present regularity results for the (non-stationary) AImage-Stokes system depending on the right hand side (in divergence form). Let us fix for this section a bounded domain GR3Image with C2Image-boundary and a time interval (0,T)Image. Moreover, let A:R3×3R3×3Image be an elliptic tensor.

B.1 The stationary problem

The AImage-Stokes problem (in the pressure-free formulation) with right hand side fL1(G)Image reads as: find vLD0,div(G)Image such that

GA(ε(v),ε(φ))dx=Gfφdxfor allφC0,div(G).

Image (B.1.1)

The right hand side can also be given in divergence form, i.e.

GA(ε(v),ε(φ))dx=GF:φdxfor allφC0,div(G)

Image (B.1.2)

for FL1(G)Image. For certain purposes it is convenient to discuss the problem with a fixed divergence. To be precise for gL1(G)Image we are seeking for a function vLD0(G)Image with divv=gImage satisfying (B.1.1) or (B.1.2). We have the following LqImage-estimates.

Lemma B.1.1

Let GR3Image be a bounded C2Image-domain and 1<q<Image.

a) Let fLq(G)Image and gW1,q(G)Image with Ggdx=0Image. Then there is a unique solution wW2,qW1,q0(G)Image to (B.1.1) such that divv=gImage and

B|2w|qdxcG|f|qdx+cG|g|qdx,

Image

where c only depends on AImage and q.

b) Let FLq(G)Image and gLq0(G)Image. Then there is a unique solution wW1,q0(G)Image to (B.1.2) such that divv=gImage and

G|w|qdxcG|F|qdx+cG|g|qdx,

Image

where c only depends on AImage and q.

In case A=IImage both parts follow from [10], Thm. 4.1. However, the main tool in [10] is the theory from [6,7] where very general linear systems are investigated. Hence it is clear that the results also hold in case of an arbitrary elliptic tensor AImage.

Corollary B.1.1

Let the assumptions of Lemma B.1.1 be satisfied. Assume further that f=0Image and

g=divg+g0

Image (B.1.3)

with gW1,q(G)Image and g0Lq0(G)Image with sptg0GImage. Then we have

G|w|qdxcG|g|qdx+cG|gN|qdH2+cL3(sptg0)βG|g0|qdx

Image

for some β>0Image.

Proof

We follow the lines of [133, Thm. 2.4]. Let ψLq(G)Image be arbitrary. In accordance to Lemma B.1.1 a) there is a unique solution uW2,qW1,q0,div(G)Image to the AImage-Stokes problem with right-hand-side ψ such that

G|2u|qdxcG|ψ|qdx.

Image

By De Rahm's theorem there is a unique ϑLq(G)Image such that

GA(ε(u),ε(φ))dx=Gϑdivφ+Gψφdxfor all φW1,q0(G).

Image

We can conclude that ϑW1,q(G)Image with

G|ϑ|qdxc(G|2u|qdx+G|ψ|qdx)cG|ψ|qdx.

Image (B.1.4)

We proceed by

Gwψdx=GA(ε(u),ε(w))dxGϑdivwdx=Gϑdivwdx=GϑdivgdxGg0ϑdx=GϑgdxGgNϑdH2Gg0ϑdx.

Image

We finish the proof by estimating this integrals separately using (B.1.4). For the first one we obtain

|Gϑgdx|(G|ϑ|qdx)1q(G|g|qdx)1q(G|ψ|qdx)1q(G|g|qdx)1q.

Image

Similarly, we estimate using the trace theorem

|GϑgNdH2|(G|ϑ|qdH2)1q(G|gN|qdH2)1q(G|ϑ|qdx)1q(G|gN|qdH2)1q(G|ψ|qdx)1q(G|gN|qdH2)1q.

Image

In case q<3Image we obtain for the third integral

|Gg0ϑdx|(G|ϑ|3q3qdx)3qdq(G|g0|3q4q3dx)4q33q(G|ϑ|qdx)qL3(sptg0)β(G|g0|qdx)1q(G|ψ|qdx)1qL3(sptg0)β(G|g0|qdx)1q,

Image

where β=13Image. If q3Image we can modify the proof by replacing 3q3qImage with an arbitrary number r>qImage and choose β=3rImage. As ψ is arbitrary, the claim follows.  □

B.2 The non-stationary problem

Now we turn to the parabolic problem and the first result is a local LqImage-estimate for weak solutions. In case of the AImage-heat system this follows from the continuity of the corresponding semigroup (see [132]). It is also known for the non-stationary Stokes-system (see [133] and [88]) but not in our setting.

Theorem B.2.48

Let fLq(Q0)Image for some q>2Image, where Q0:=(0,T)×R3Image and let vL2((0,T);W1,20,div(R3))Image be the unique weak solution to

Q0vtφdxdtQ0A(ε(v),ε(φ))dx=Q0fφdxdt

Image (B.2.5)

for all φC0,div([0,T)×R3)Image. Then we have 2vLqloc(Q0)Image and there holds

T0B|2v|qdxdtcBQ0|f|qdxdt

Image

for all balls BR3Image.

Proof

The main ingredient is the proof of the following auxiliary result which has been used in a similar version in [42].

i) We start with interior estimates. Let

Qr:=Qr(x0,t0):=(t0r2,t0+r2)×Br(x0)

Image

be a parabolic cylinder such that 4QrQ0Image. We claim the following: There is a constant N1>0Image such that for every ε>0Image there is δ=δ(ε)>0Image such that (since fL2(0,T;L2loc(R3))Image the standard interior regularity theory implies 2vL2(0,T;L2loc(R3))Image and tvL2(0,T;L2loc(R3))Image)

L4(Qr{M(|f|)2>N21})εL4(Qr)Qr{M(|2v|2)>1}{M(|f|2)>δ2}.

Image (B.2.6)

Let us assume for simplicity that r=1Image. In fact, we will establish (B.2.6) by showing

Q1{M(|2v|2)1}{M(|f|2)δ2}L4(Q1{M(|v|)2>N21})<εL4(Q1)

Image (B.2.7)

and applying a simple scaling argument. In order to show (B.2.7) we compare v with a solution to a homogeneous problem on

Q4=(t4042,t40+42)×B4(x40)Q0

Image

which is smooth in the interior. So let us define h as the unique solution to

{thdivA(ε(h))=πhin Q4,divh=0in Q4,h=von I4×B4,h(t40,)=v(t40,)in B4.

Image (B.2.8)

We test the difference of both equations with vhImage. This yields by the ellipticity of AImage

suptI4B4|v(t)h(t)|2dx+Q4|ε(v)ε(h)|2dxdtcQ4|f|2dxdt+cQ4|vh|2dxdt.

Image

An application of Korn's inequality and Gronwall's lemma implies

suptI4B4|v(t)h(t)|2dx+Q4|vh|2dxdtcQ4|f|2dxdt.

Image (B.2.9)

First we insert t(vh)Image which yields similarly

Q4|t(vh)|2dx+suptI4B4|vh|2dxcQ4|f|2dxdt.

Image (B.2.10)

We can introduce the pressure terms πv,πhL2(I4,L20(B4))Image in the equations for v and h and show

Q4|πvπh|2dxcQ4|f|2dxdt.

Image (B.2.11)

Estimate (B.2.11) can be shown by using the Bogovskiĭ  operator, see Section 2.1. Setting Bog=BogB4Image we gain due to (B.2.10) for any φC0(Q4)Image that

Q4(πvπh)φdxdt=Q4(πvπh)divBog(φφB4)dxdt=Q4A(ε(vh),ε(Bog(φφB4)))dxdt+Q4fBog(φφB4))dxdtQ4t(vh)Bog(φφB4)dxdtc((vh)2+f2+t(vh)2)Bog(φ(φ)B4)2c(Q4|f|2dxdt)12(Q4|φ|2dxdt)12.

Image

Now we choose a cut off function ηC0(B4)Image with 0η1Image and η1Image on B3Image. We insert γ(η2γ(vh))Image in the equation for vhImage and sum over γ{1,2,3}Image. We gain

suptI4B4η2|(vh)|2dx+Q4η2|ε(v)ε(h)|2dxdtcQ4fγ(η2γ(hh))dxdt+c(η)Q4|vh|2dxdt+cQ4(ππh)γ(η2γ(hh))dxdt.

Image

We estimate the term involving f by

Q4fγ(η2γ(vh))dxdtc(κ)Q4|f|2dxdt+κQ4η2|2v2h|2dxdt+c(η)Q4|vh|2dxdt,

Image

where κ>0Image is arbitrary. The term involving ππhImage can be estimated in the same fashion. Choosing κ>0Image small enough and using the inequality |2u|c|ε(u)|Image as well as (B.2.9)(B.2.11) shows

suptI4B3|(vh)(t)|2dx+I4B3|2v2h|2dxdtcQ4|f|2dxdt.

Image (B.2.12)

Now, let us assume that (B.2.7)1 holds. Then there is a point (t0,x0)Q1Image such that

Qσ(t0,x0)|2v|2dxdt1,Qσ(t0,x0)|f|2dxdtδ2

Image (B.2.13)

for all σ>0Image. Since Q4Q6(t0,x0)Image we have

Q4|2v|2dxdtc,Q4|f|2dxdtcδ2.

Image (B.2.14)

As h is smooth we know that

N20:=supQ3|2h|2<.

Image (B.2.15)

From this we aim to conclude that

Q1{M(|2v|2)>N21}Q1{M(χQ3|2v2h|2)>N20}

Image (B.2.16)

for N21:=max{4N20,25}Image. To establish (B.2.16) suppose that

(t,x)Q1{M(χQ3|2v2h|2)N20}.

Image (B.2.17)

If σ2Image we have Qσ(t,x)Q3Image and gain by (B.2.15)

Qσ(t,x)|2v|2dxdt2Qσ(t,x)χQ3|2v2h|2dxdt+2Qσ(t,x)χQ3|2h|2dxdt4N20.

Image

If σ2Image we have by (B.2.13)

Qσ(t,x)|2v|2dxdt25Q2σ(t0,x0)|2v|2dxdt25.

Image

Combining the both cases yields (B.2.16). This implies together with the continuity of the maximal function on L2Image, (B.2.12) and (B.2.14)

L4(Q1{M(|2v|2)>N21})Ld+1(Q1{M(χQ3|2v2h|2)>N20})cN20Q3|2v2h|2dxdtcN20Q3|f|2dxdtcN20δ2=εL4(Q1),

Image

choosing δ:=c1/2N0εImage. So we have shown (B.2.7) which yields (B.2.6) by a scaling argument.

If (B.2.6)1 holds then we have

L4(Qr{M(|2v|)2>N21})εL4(Qr{M(|2v|2)>1}{M(|f|2)>δ2})ε(L4(Qr{M(|2v|2)>1})+L4(Qr{M(|f|2)>δ2})).

Image

Multiplying the equation for v by some small number ϱ=ϱ(fq,2v2)Image we can assume that

L4(Qr{M(|2v|2)>N21})<ε.

Image (B.2.18)

By induction we can establish that

L4(Qr{M(|2v|2)>N2k1})εkL4(Qr{M(|2v|2)>1}+cki=1εiL4(Qr{M(|f|2)>δ2N2(ki)1}).

Image

In the induction step one has to introduce v1:=vN1Image which is a solution to the AImage-Stokes problem with right hand side f1:=fN1Image. Now we will show 2vLq(Qr)Image. For this we use the equivalence for 1q0<Image

k=1Lq0kL4(Qr{(t,x):|u(t,x)|>θLk})<uLq0(Qr),

Image

which holds for any measurable function u, see [44]. Here L>1Image and θ>0Image are arbitrary. So we aim to prove that

k=1(N21)q2kL4(Qr{M(|2v|2)>N2k1})<

Image (B.2.19)

to conclude M(|2v|2)Lq2(Qr)Image and hence 2vLq(Qr)Image. Since fLq(Qr)Image we have M(|f|2)Lq/2(Qr)Image (recall (5.2.4)) and we have

k=1(N21)q2kL4(Qr{M(|f|2)>δ2N2k1})<.

Image

We obtain

k=1(N1)qkL4(Qr{M(|2v|2)>N2k1})ck=1(N1)qkεkL4(Qr{M(|2v|2)>1})+ck=1(N1)qkki=1εiL4(Qr{M(|f|2)>δ2N2(ki)1})crk=1(εN1)qk+ci=1εi(N1)qik=i(N1)q(ki)L4(Qr{M(|f|2)>δ2N2(ki)1})crk=1(εNq1)k.

Image

If we choose εNq1<1Image the sum in (B.2.19) is converging and we have 2vLq(Qr)Image. Since the mapping f2vImage is linear we gain the desired estimate

Qr|2v|qdxdtcrQ0|f|qdxdt.

Image (B.2.20)

ii) Now let Q1Image be a cylinder such that 4Q1(,0]×R3Image. Moreover, assume that Q1Q0Image. We consider the solution ˜hImage to

{t˜hdivA(ε(˜h))=π˜hin ˜Q4,div˜h=0in ˜Q4,˜h=von ˜I4×B4,˜h(t40,)=0in B4,

Image (B.2.21)

where ˜Im:=Im(0,T)Image and ˜Qm:=˜Im×BmImage. We can establish a variant of (B.2.12) on ˜Q4Image. Now we have sup˜Q3|2˜h|2<Image due to the smooth initial datum of ˜hImage (recall that v(0,)=0Image a.e.). So we can finish the proof as before and gain 2vLq(Q1)Image. This implies again (B.2.20).

iii) The situation 4Q1[T,)×R3Image is uncritical again and we can assume that ii) and iii) do not occur for the same cylinder (by choosing sufficiently small cubes).

Covering the set (0,T)×BImage by smaller cylinders and combing i)–iii) yield the desired estimate.  □

Corollary B.2.1

Under the assumptions of Theorem B.2.48 we have for all balls BR3Image the following estimates for some constant cBImage which does not depend on T.

a) The following holds

T0B(|vT|q+|vT|q+|2v|q)dxdtcBQ0|f|qdxdt.

Image

b) We have tvLq(0,T;Lqloc(R3))Image together with

T0B|tv|qdxdtcBQ0|f|qdxdt.

Image

c) There is πLq((0,T),W1,qloc(R3))Image such that

Q0vtφdxdtQ0A(ε(v),ε(φ))dx=Q0πdivφdxdt+Q0fφdxdt

Image

for all φC0([0,T)×R3)Image.

d) The following holds

T0B(|πT|q+|π|q)dxdtcBQ0|f|qdxdt.

Image

Proof

The estimate in a) is a simple scaling argument. Having a solution v defined on (0,T)×R3Image we gain a solution ˆvImage on (0,1)×R3Image by setting

ˆv(s,x):=1Tv(Ts,Tx).

Image

Now we apply Theorem B.2.48 to ˆvImage. The constant which appears is independent of T. Transforming back to v yields the claimed inequality.

b) For φC0(Q0)Image with φ(t,x)=τ(t)ψ(x)Image where φC0(G)Image (GR3Image a bounded Lipschitz domain) we have

Q0vtφdxdt=T0tτR3v(ψdiv+Ψ)dxdt=T0tτR3vψdivdxdt=Q0vtφdivdxdt

Image

where ψdiv:=ψΔ1GdivψImage and Ψ:=Δ1GdivψImage. (Here Δ1GImage is the solution operator to the Laplace equation with zero boundary datum on ∂G.) Here we took into account Ψ|G=0Image as well as divv=0Image. Using 2vL2(0,T;L2loc(R3))Image we proceed by

Q0vtφdxdt=T0G(fdivAε(v))φdivdxdtc(T0G(|2v|q+|f|q)dxdt)1q(T0G|φdiv|qdxdt)1qc(Q0|f|qdxdt)1q(T0G|φ|qdxdt)1q.

Image

In the last step we used the estimate from Theorem B.2.48 and continuity of Δ1GdivImage on Lq(G)Image. Duality implies tvLq(0,T;Lqloc(R3))Image and we can introduce the pressure function πLq(0,T;Lqloc(R3))Image as claimed in b) by De Rahm's Theorem. Using the equation for v and the estimates in a) and b) we gain

Q|π|qdxdtcQ0|f|qdxdt.

Image

The estimate for π in d) follows again by scaling.  □

Corollary B.2.2

Let fLq(Q+0)Image for some q>2Image where Q+0:=(0,T)×R3+Image, R3+=R3[x3>0]Image, and let vL2((0,T);W1,20,div(R3+))Image with v|x3=0=0Image be the unique weak solution to

Q+0vtφdxdtQ+0A(ε(v),ε(φ))dx=Q+0fφdxdt

Image (B.2.22)

for all φC0,div([0,T)×R3+)Image. Then the results from Theorem B.2.48 and Corollary B.2.1 hold for v for all half balls B+(z)=B(z)[x3>0]R3Image with z3=0Image.

Proof

We will show a variant of the LqImage-estimate from Theorem B.2.48 on half balls B+Image, i.e.

T0B+|2v|qdxdtcBT0Q+|f|qdxdt.

Image (B.2.23)

From this we can follow estimates in the fashion of Corollary B.2.1 as done there. In order to establish (B.2.23) we will proceed as in the proof of Theorem B.2.48 replacing all balls with half ball. So let Q1R4Image such that 4Q1Q0Image (the other situation can be shown along the modifications indicated at the end of the proof of Theorem B.2.48). Moreover, assume that Q1=I1×B1(z)Image where z3=0Image. We compare v with the unique solution h+Image to

{th+divA(ε(h+))=πh+in Q+4,divh+=0in Q+4,h+=von I4×B+4,h+(0,)=v(0,)in B+4.

Image (B.2.24)

We gain a version of the estimate (B.2.9) and (B.2.10) on half-balls. In fact, there holds

suptI4B+4|vh+|2dx+Q+4|vh+|2dxdtcQ+4|f|2dxdt,

Image (B.2.25)

Q+4|t(vh+)|2dx+suptI4B+4|vh+|2dxcQ+4|f|2dxdt.

Image (B.2.26)

We can introduce the pressure terms πv,πh+L2(I4,L20(B+4))Image in the equations for v and h+Image and show

Q+4|πvπh+|2dxcQ+4|f|2dxdt.

Image (B.2.27)

This can be done as in the proof of (B.2.11) using the Bogovskiĭ  operator on B+4Image. Estimate (B.2.27) can be shown by using the Bogovskiĭ  operator introduced in [24]. Now we insert γ(η2γ(vh))Image for γ{1,2}Image in the equation for vhImage. Here we choose ηC0(B4)Image with 0η1Image and η1Image on B3Image. This yields together with (B.2.25)(B.2.27)

Q+3|˜(vh)|2dxcQ+4|f|2dxdt,

Image (B.2.28)

where ˜:=(1,2)Image. Finally, the only term which is missing is 23(vh)Image. On account of div(vh)=0Image we have (cf. [21])

|23(vh)|c(|˜(πvπh+)|+|˜(vh)|+|f|).

Image (B.2.29)

So we have to estimate derivatives of the pressure. In fact we have

Q+3|˜(πvπh+)|2dxcQ+4|f|2dxdt.

Image (B.2.30)

We can show this similarly to the proof of (B.2.27) replacing φ by γφImage and using (B.2.28). Combining (B.2.28)(B.2.30) implies

Q+3|2v2h|2dxcQ+4|f|2dxdt.

Image

Moreover, we know supQ+3|2h+|2<Image. Note that h+=0Image on Q3[x3=0]Image. This allows to show 2vLq(Q+1)Image as in the proof of Theorem B.2.48.  □

Corollary B.2.3

Let fLq(Qν,ξ0)Image for some q>2Image where Qν,ξ0:=(0,T)×R3ν,ξImage, R3ν,ξ=R3[(xξ)ν>0]Image for some ν,ξR3Image. Let vL2(0,T;W1,20,div(R3ν,ξ))Image with v|(xξ)ν=0=0Image be the unique weak solution to

Qν,ξ0vtφdxdtQν,ξ0A(ε(v),ε(φ))dx=Qν,ξ0fφdxdt

Image (B.2.31)

for all φC0,div([0,T)×R3ν,ξ)Image. Then the results from Theorem B.2.48 and Corollary B.2.1 hold for v for all half balls Bν,ξ(z)=B(z)[(xξ)ν>0]R3Image with (zξ)ν=0Image.

Proof

The proof follows easily from Corollary B.2.2 by rotation of the coordinate system. There is an orthogonal matrix VR3×3Image such that the mapping z=V(xξ)Image transforms R3ν,ξImage to R3+Image. We define

˜v(t,x)=V1v(t,V(xξ)),˜f(t,x)=V1f(t,V(xξ)),

Image

as well as the bilinear form ˜AImage by

˜A(ζ,ξ):=A(VζV1,VξV1),ζ,ξR3×3.

Image

Note that the ellipticity constants of AImage and ˜AImage coincide as V is an orthogonal matrix. Now it is easy to see that ˜vL2(0,T;W1,20,div(R3+))Image satisfies ˜v|x3=0=0Image and is a solution to (B.2.22) with right hand side ˜fImage and bilinear form ˜AImage. Hence Corollary B.2.3.  □

Theorem B.2.49

Let Q:=(0,T)×GImage with a bounded domain GR3Image having a C2Image-boundary. Let fLq(Q)Image for some q>2Image. Then there is a unique weak solution vL(0,T;L2(G))Lq(0,T;W1,q0,div(G))Image to

QvtφdxdtQA(ε(v),ε(φ))dx=Qfφdxdt

Image (B.2.32)

for all φC0,div([0,T)×G)Image such that 2vLq(Q)Image. Moreover, we have

Q|2v|qdxdtcQ|f|qdxdt.

Image

Proof

Due to the local LqImage-theory for the whole space problem and the half-space problem which follow from Corollary B.2.1 and Corollary B.2.2 (with the right scaling in T) the proof is similar to [133], Thm. 4.1, in the case A=IImage. Note that LqImage-estimates for the stationary problem on bounded domains with given divergence are stated in Lemma B.1.1. We want to invert the operator

L:YLq(I;Lqdiv(G)),vPq(tvdivAε(v)).

Image

The space YImage is given by

Y:=Lq(I;W1,q0,divW2,q(G))W1,q(I;Lq(G)){v(0,)=0,vG=0}

Image

and PqImage is the Helmholtz projection from Lq(G)Image into Lqdiv(G)Image. The latter one is defined by

Lqdiv(G):=C0,div(G)q.

Image

The Helmholtz-projection PquImage of a function uLq(G)Image can be defined as Pqu:=uhImage, where h is the solution to the Neumann-problem

{Δh=divuonG,NB(hu)=0onG.

Image

We will try to find an operator R:Lq(I;Lqdiv(G))YImage such that

LR=I+τ

Image (B.2.33)

with τ<1Image. The range of LRImage (which then is Lq(0,T;Lqdiv(G))Image) is contained in the range of LImage. So LImage is onto.

Let GkImage, 0=1,...,NImage be a covering of G such that G0GImage and GkImage covers a (small) boundary strip of G. There are local coordinates

z=Zk(yk)=(yk1,yk2,yk3Fk(yk1,yk2)),yk=(yk1,yk2)Bλ(ξk)

Image

where FkImage is a C2Image-function and 0<λλ0<1Image. Here ξkImage denotes the center point of Sk=GGkImage and νkImage the outer unit normal of SkImage at the point ξkImage. In this coordinate system we have a flat boundary which is contained in the plane {(xξk)νk=0}Image. We consider a decomposition of unity (ζk)Nk=0C0(R3)Image with respect to GkImage such that sptζkGkImage. We can assume that |lζk|cλlImage for l=1,2Image and that the multiplicity of the covering of G by the domains GkImage does not depend on λ.

Furthermore, M0fImage and MfImage is the extension of a function f (by zero) to the whole space or the half space respectively. Note that if fLqdiv(G)Image, then M0fLqdiv(R3)Image. Finally, we denote by UkfImage (for fLq(R3+)Image) and U0fImage (for fLq(R3)Image) the solution on the half space (corresponding to the plane {(xξk)νk=0}Image) and the whole space respectively (see Corollary B.2.1 and Corollary B.2.3). By QkfImage we denote the pressure corresponding to UkfImage. Now we define the operators

R0f:=ζ0U0M0f+Nk=1ζkZ1kUkMZkf,P0f:=Nk=1Z1kσQkMZkf.

Image

The idea is in the interior to extend the force to the whole space, compute the whole space solution and localize again. At the boundary it is more tricky since we have to flatten the problem before considering the half space problem (and of course we have to transform back after solving it). Due to the involved cut-off functions R0fImage is in general not divergence-free. This will be corrected in the following way: we set

Rf=R0f+R1f,Pf=P0f+P1f.

Image

Here R1f=wImage and P1f=sImage, where (w,s)Image is the unique solution to the stationary problem

divAε(w)+s=0,divw=divR0f,w|G=0,

Image

cf. Lemma B.1.1. Now we clearly have RfYImage. We need to establish (B.2.33). We abbreviate

uk:=Z1kUkMZkf,πk:=Z1kQkMZkf.

Image

We define ˜k=x+ZkFkz3Image and accordingly ˜divkImage as well as ˜εkImage. We gain on QTkImage

tuk˜divkA(˜εk(uk))+˜kπk=Z1kZkf=f.

Image

There holds

tRfdivAε(Rf)+Pf=f+Sf+tR1f,

Image (B.2.34)

Sf=Aε(U0M0f)ζ0divA(ζ0U0M0f)Nk=1Aε(uk)ζkNk=1divA(ζkuk)+Nk=1ζkπkNk=1ζk(˜k)πk+Nk=1ζk(˜divkA(˜εk())divA(ε()))uk.

Image (B.2.35)

From (B.2.34) it follows

LRf=f+PqSf+PqtR1f

Image

i.e. (B.2.33) with τ=PqS+PqtR1Image. We need to estimate the norms of the operators SImage and tR1Image. If we choose T small enough the first two terms in (B.2.35) are small in accordance to Corollary B.2.1. The same is true for the first three sums as a consequence of Corollary B.2.3. All together we have

5i=1Tiqδ(λ,T)fq

Image

with δ(λ,T)0Image for T0Image (and any fixed λ). Note that δ(λ,T)Image does not depend on N on account of the localization. We will argue similarly for the next three sums assuming that the gradients of the FkImage are small (meaning that λ is small). Here, we gain

T6q+T7qκ(λ)fq

Image

with κ(λ)0Image for λ0Image. Note that κ(λ)Image does neither depend on T nor on N. By choosing first λ small enough such that κ(λ)18Image and then T small enough such that δ(λ,T)18Image we can follow

Sfq14fq.

Image (B.2.36)

Now we are going to show the same for tR1Image. In order to achieve this we consider the function w=tR1fImage which is the solution to

divAε(tw)+s=0,divw=divtR0f,w|G=0.

Image

We have the identity

divtR0f=Nk=0ζktuk=Nk=0˜kζktuk+Nk=0(˜k)ζktuk=Nk=1˜kζk(˜divkA(˜εk(uk))˜kπk)+Nk=0(˜k)ζkf+Nk=1(˜k)ζktuk=:T1+T2+T3,

Image

where u0=U0M0fImage and ˜0=Image. We use the formula

T1=Nk=1˜divk(˜kζkA(˜εk(uk))˜kζkπk)Nk=1(˜2kζkA(˜εk(uk))˜Δkζkπk)=:T11+T21,

Image

as well as

(˜k)gk=ZkFkz3gk=Zkz3(Fkgk)=νk(ZkFkgk)=νk(ZkFkgk)=div(νkZ1k(Fkgk)).

Image

Now we can apply Corollary B.1.1 with divg=T11Image and g0=T21+T2+T3Image. By Corollary B.2.2 there are constants δ(λ,T),δ(λ,T)Image with δ(λ,T)0Image and δ(λ,T)0Image for T0Image (and any fixed λ) such that

gqδ(λ,T)fq,g0qc(1+δ(λ,T))fq.

Image

As a consequence of Corollary B.1.1 this yields

tR1fqc(λβ+δ(λ,T)+δ(λ,T))fq.

Image

Choosing first λ and then T small enough we gain

tR1fq14fq.

Image (B.2.37)

Combining (B.2.36) and (B.2.37) implies τ12Image. Hence LImage is onto (recall (B.2.33)). This means we have shown the claim for T sufficiently small, say T=T01Image. It is easy to extend it to the whole interval. Let (v,π)Image be the solution on [0,T0]Image. We extend it in an even manner to the interval [0,2T0]Image. On the interval [T0,2T0]Image we define (v,π)Image as solution to the AImage-Stokes system with right-hand-side

f(t,x)=f(t,x)f(x,2T0t)+2tv(x,2T0t).

Image

If we set v=0,π=0Image on [0,T0]Image then (v+v,π+π)Image is the solution on [0,2T0]Image. This can be repeated to construct the solution on [0,T]Image.  □

B.3 The non-stationary problem in divergence form

In order to treat problems with right hand side in divergence form we consider the AImage-Stokes operator

Aq:=PqdivA(ε()).

Image

The AImage-Stokes operator AqImage enjoys the same properties than the Stokes operator AqImage (see for instance [87]).

For the AImage-Stokes operator it holds D(Aq)=W1,q0,divW2,q(G)Image, where D denotes the domain, and

u2,qc1Aquqc2u2,q,uD(Aq),

Image (B.3.38)

GAquwdx=GuAqwdxuD(Aq),wD(Aq).

Image (B.3.39)

Inequality (B.3.38) is a consequence of Lemma B.1.1 a) and the continuity of PqImage.

Since AqImage is positive its root A12qImage is well-defined with D(A12q)=W1,q0,div(G)Image and

u1,qc1A12quqc2u1,q,uD(A12q),

Image (B.3.40)

GA12quwdx=GuA12qwdxuD(A12q),wD(A12q).

Image (B.3.41)

Finally, the inverse operator A12q:Lqdiv(G)W1,q0,div(G)Image is defined and it holds

A12quqcuq,uD(A12q),

Image (B.3.42)

GA12quwdx=GuA12qwdxuD(A12q),wD(A12q).

Image (B.3.43)

From the definition of the square root of a positive self-adjoint operator follows also that

A12q:W2,qW1,q0,div(G)W1,q0,div(G),A12q:W1,q0,div(G)W2,qW1,q0,div(G),

Image

together with

A12quqcu2,q,uW2,qW1,q0,div(G),

Image (B.3.44)

A12quqcuq,uW1,q0,div(G).

Image (B.3.45)

Finally we state the main result of this section.

Theorem B.3.50

Let Q:=(0,T)×GImage with a bounded domain GR3Image having a C2Image-boundary. Let FLq(Q)Image, where q(1,)Image. There is a unique solution wLq(0,T;W1,q0,div(G))Image to

QwtφdxdtQA(ε(w),ε(φ))dxdt=QF:φ

Image (B.3.46)

for all φC0,div([0,T)×G)Image. Moreover we have

Q|w|qdxdtcQ|F|qdxdt,

Image

where c only depends on AImage and q.

Proof

Let us first assume that q>2Image. Then Theorem B.2.49 applies. We set f:=A12qdivFImage which is defined via the duality

GA12qdivFφdx=GF:A12qφdx,φC0,div(G),

Image

using (B.3.43). So we gain fLqdiv(G)Image with

fqcFq.

Image (B.3.47)

We define ˜wLq(0,T;W1,q0,div(G))Image as the unique solution to

Q˜wtφdxdtQA(ε(˜w),ε(φ))dxdt=Qfφdxdt

Image (B.3.48)

for all φC0,div([0,T)×G)Image. Theorem B.2.49 yields ˜wLq(0,T;W2,q(G))Image and

˜w2,qcfq.

Image (B.3.49)

We want to return to the original problem and set w:=A12q˜wImage thus we have wLq(0,T;W1,q0,div(G))Image. Since A12q:W1,q0,divW2,q(G)W1,q0,div(G)Image we can replace φ by A12qφImage in (B.3.48). This implies using (B.3.41) and the definition of f

Qwtφdxdt+QdivA(ε(˜w)):A12qφdxdt=QF:φ

Image

for all φC0,div(Q)Image. On account of A12qφW1,q0,div(G)Image and A12q˜wW1,q0,div(G)Image we gain due to (B.3.41)

QdivA(ε(˜w)):A12qφdxdt=QAq˜w:A12qφdxdt=QA12q˜w:Aqφdxdt=QwdivA(ε(φ))dxdt=QA(ε(w),ε(φ))dxdt

Image

using (B.3.41) and wW1,q0,div(G)Image. This shows that w is the unique solution to (6.2.27). Moreover, we obtain the desired regularity estimate via

Q|w|qdxdtcQ|A12qw|qdxdt=cQ|Aq˜w|qdxdtcQ|2˜w|qcQ|f|qdxdtcQ|F|qdxdt

Image

as a consequence of (B.3.40), the definition of w, (B.3.38), (B.3.49), and (B.3.47). A simple scaling argument shows that the inequality is independent of the diameter of I and B. So we have shown the claim for q>2Image.

The case q=2Image follows easily from a priori estimates and Korn's inequality. So let us assume that q<2Image. Duality arguments show that

1qQ|w|qdxdt=supGLq(Q)[Qw:Gdxdt1qQ|G|qdxdt].

Image

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

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

Image (B.3.50)

for all ξC0,div((0,T]×G)Image. This is a backward parabolic equation with end datum zero. We have that tzGLq(0,T;W1,qdiv(G))Image such that test-functions can be chosen from the space Lq(0,T;W1,q0,div(G))Image. Due to q>2Image the first part of the proof (applied to ˜z˜G(t,)=zG(Tt,)Image, where ˜G(t,)=G(Tt,)Image) yields

Q|zG|qdxdtcQ|G|qdxdt.

Image

This and w(0,)=0Image implies (using w as a test-function in (B.3.50))

Q|w|qdxdtcsupGLq(Q)[QA(ε(w),ε(zG))dxdtQtzGwdxdtQ|zG|qdxdt]csupξC0,div(Q)[QA(ε(w),ε(ξ))dxdtQwtξdxdtQ|ξ|qdxdt].

Image

The equation for w and Young's inequality finally give

Q|w|qdxdtcsupξC0,div(Q)[QF:ξdxdtQ|ξ|qdxdt]cQ|F|qdxdt

Image

and hence the claim.  □

References

[6] S. Agmon, A. Douglis, L. Nirenberg, Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions. I, Commun. Pure Appl. Math. 1959;12:623–727.

[7] S. Agmon, A. Douglis, L. Nirenberg, Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions. II, Commun. Pure Appl. Math. 1964;17:35–92.

[10] C. Amrouche, V. Girault, Decomposition of vector spaces and application to the Stokes problem in arbitrary dimension, Czechoslov. Math. J. 1994;44:109–140.

[21] H. Beirão da Veiga, P. Kaplický, M. Růžička, Boundary Regularity of Shear Thickening Flows, J. Math. Fluid Mech. 2011;13:387–404.

[24] M.E. Bogovskiĭ, Solutions of some problems of vector analysis, associated with the operators div and grad, In: Theory of Cubature Formulas and the Application of Functional Analysis to Problems of Mathematical Physics. Novosibirsk: Akad. Nauk SSSR Sibirsk. Otdel. Inst. Mat.; 1980:5–40 149 (in Russian).

[42] S.-S. Byun, L. Wang, LpImage estimates for parabolic equations in Reifenberg domains, J. Funct. Anal. 2005;223:44–85.

[44] L.A. Caffarelli, X. Cabré, Fully Nonlinear Elliptic Equations. Colloq. Publ. – Am. Math. Soc.. Providence, RI: American Mathematical Society; 1995;vol. 43 vi+104 pp.

[87] G.P. Galdi, C.G. Simader, H. Sohr, A class of solutions to stationary Stokes and Navier–Stokes equations with boundary data in W1/q,qImage, Math. Ann. 2005;331:41–74.

[88] M. Giaquinta, Multiple Integrals in the Calculus of Variations and Nonlinear Elliptic Systems. Ann. Math. Stud.. Princeton, New Jersey: Princeton University Press; 1983;vol. 105.

[132] Z. Shen, Resolvent Estimates in LpImage for eliptic systems in Lipschitz domains, J. Funct. Anal. 1995;133:224–251.

[133] V.A. Solonnikov, Estimates for solutions of nonstationary Navier–Stokes equations, J. Sov. Math. 1977;8:467–528.

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

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