Skip to main content

Optimal partial regularity of second-order parabolic systems under natural growth condition

Abstract

We consider the regularity for weak solutions of second-order nonlinear parabolic systems under a natural growth condition when m>2, and obtain a general criterion for a weak solution to be regular in the neighborhood of a given point. In particular, we get the optimal regularity by the method of A-caloric approximation introduced by Duzaar and Mingione.

1 Introduction

Electrorheological fluids are special viscous liquids, that are characterized by their ability to undergo significant changes in their mechanical properties when an electric field is applied. This property can be exploited in technological applications, e.g., actuators, clutches, shock absorbers, and rehabilitation equipment to name a few [1].

A model was developed for these liquids within the framework of rational mechanics [2, 3]; it takes into account the complex interactions between the electro-magnetic fields and the moving liquid. If the fluid is assumed to be incompressible, it turns out that the relevant equations of the model are the system

div(E+P)=0,
(1.1)
curlE=0,
(1.2)
ρ 0 v t divS+ ρ 0 [v]v+ϕ= ρ 0 f+[E]P,
(1.3)
divv=0,
(1.4)

where E is the electric field, P is the polarization, ρ 0 is the density, v is the velocity, S is the extra stress, ϕ is the pressure, and f is the mechanical force. In fact, in a model capable of explaining many of the observed phenomena, the extra stress has the form

S = α 21 ( ( 1 + | D | 2 ) p 1 2 1 ) E E + ( α 31 + α 33 | E | 2 ) ( 1 + | D | 2 ) p 2 2 D + α 51 ( 1 + | D | 2 ) p 2 2 ( D E E + E D E ) ,
(1.5)

where α i j are material constants, and where the material function p depends on the strength of the electric field | E | 2 and satisfies

1< p p ( | E | 2 ) p 0 <.
(1.6)

Since the material function p, which essentially determines S, depends on the magnitude of the electric field | E | 2 , we have to deal with an elliptic or parabolic system of partial differential equations with the so-called non-standard growth conditions, i.e., the elliptic operator S satisfies

S(D,E)D c 0 ( 1 + | E | 2 ) ( 1 + | D | 2 ) p 2 2 | D | 2 ,
(1.7)
| S ( D , E ) | c 1 ( 1 + | D | 2 ) p 0 1 2 | E | 2 .
(1.8)

Equality (1.5) of electrorheological fluids with the conditions (1.7) and (1.8) encouraged us to considered the partial regularity of a more simple and standard model as the following:

u t i α = 1 n D α A i α (z,u,Du)= B i (z,u,Du),i=1,2,,N,
(1.9)

where Ω R n is a bounded domain and T>0, z=(x,t) with xΩ, 0<tT, denote a point in Q T =Ω×(T,0). Let u(z)=( u 1 (z), u 2 (z),, u N (z)) be a vector-valued function defined in Q T . Denote by Du the gradient of u, i.e., Du= { D α u i } i = 1 , , N ; α = 1 , , n . m>2 is a real number.

In order to define the weak solution of (1.9), one needs to impose some regularity conditions and constructer conditions to A i α and B i . For a vector field A i α : Q T × R N × R n N , we shall denote the coefficients by A i α (z,u,p)= A i α (x,t,u,p) if z=(x,t), u R N and p R n N . We assume that the functions (z,u,p) A i α (z,u,p); (z,u,p) A i α p β j (z,u,p) are continuous in Q T × R N × R n N and that the following growth and ellipticity conditions are satisfied:

(H1) There exists a constant L such that

| A i α ( z , u , p ) | L ( 1 + | p | ) m 2 for all z Q T ,u R n  and p R n N .

(H2) A i α (z,u,p) are differentiable functions in p and there exists a constant L such that

| A i α p β i (z,u,p)|L ( 1 + | p | 2 ) m 2 2 for all z Q T ,u R n  and p R n N .

(H3) A i α is uniformly strongly elliptic, that is, for some λ>0, we have

( A i α p β i ( z , u , p ) p ˜ α i ) p ˜ β j λ | p ˜ | 2 ( 1 + | p | 2 ) m 2 2 for all z Q T ,u R n  and p, p ˜ R n N ,

where λ>0 and 1L<. Now we shall specify the regularity assumptions on A i α (z,u,p) with respect to the ‘coefficient’ (z,u) and assume that the function (z,u) A i α ( z , u , p ) 1 + | p | is Hölder continuous with respect to the parabolic metric | x x 0 | 2 + | t t 0 | with Hölder exponent β(0,1) but not necessarily uniformly Hölder continuous; namely we shall assume that:

(H4) There exists a constant L such that

| A i α (z,u,p) A i α ( z 0 , u 0 ,p)|Lθ ( | u | + | u 0 | , | x x 0 | + | t t 0 | + | u u 0 | ) ( 1 + | p | ) m 2

for any z=(x,t) and z 0 =( x 0 , t 0 ) in Q T . u and u 0 in R n and for all p R n N , where θ(y,s)=min{1, K ˜ (y) s β }, K ˜ :[0,)(1,) is a given non-decreasing function. Note that θ is concave in the argument. This is the standard way to prescribe (non-uniform) Hölder continuity of the function A i α (z,u,p). We find it a bit difficult to handle, therefore, in many points of the paper, we shall use:

(H4′) For β(0,1) and K:[0,)[L,) monotone nondecreasing such that

| A i α ( z , u , p ) A i α ( z 0 , u 0 , p ) | K ( | u | ) ( | x x 0 | + | t t 0 | + | u u 0 | ) β ( 1 + | p | ) m 2 ,

valid for any z=(x,t) and z 0 =( x 0 , t 0 ) in Q T , u and u 0 in R n and p R n N .

(H5) There exist constants a and b such that

| B i ( z , u , p ) | a | p | m +b, sup Q T |u|=V,2aV<λ.

Finally, we remark a trial consequence of the continuity of A i α p β j ; this implies the existence of a function ω:[0,)×[0,)[0,) with ω(t,0)=0 for all t such that tω(t,s) is nondecreasing for fixed s, s ω m (t,s) is concave and nondecreasing for fixed t, and such that

( H 6 ) | A i α p β j ( x , t , u , p ) A i α p β j ( x 0 , t 0 , u 0 , p 0 ) | L ( 1 + | p | 2 + | p 0 | 2 ) m 2 2 × ω ( | u | + | p | , | x x 0 | 2 + | t t 0 | + | u u 0 | 2 + | p p 0 | 2 )

for any z=(x,t) and z 0 =( x 0 , t 0 ) in Q T , any u, u 0 in R n and p, p 0 R n N whenever |u|+|p|+|u u 0 |+|p p 0 |M.

From (H2) and (H3) we immediately deduce the following:

| A i α (z,u,p) A i α (z,u,q)|L ( 1 + | p | 2 + | q | 2 ) m 2 2 |pq|,
(1.10)
( A i α ( z , u , p ) A i α ( z , u , q ) ) (pq)λ ( 1 + | p | 2 + | q | 2 ) m 2 2 | p q | 2
(1.11)

for all z Q T , u R N and p,q R n N .

Definition 1.1 By a weak solution of (1.9) under the assumptions (H1)-(H5), we mean a vector-valued function u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) such that

Q T ( A i α ( z , u , D u ) D α φ i u i φ t i ) dz= Q T B i (z,u,Du) φ i dz
(1.12)

for all φ C 0 ( Q T , R N ).

In [4] Duzaar and Mingione considered the partial regularity of homogeneous systems of (1.9) with m2 under the natural growth condition. In this paper, we extend their results to the case of m>2. We have to overcome the difficulty of m>2. Motivated by the works of Duzaar [4, 5], Chen and Tan [69] and Tan [10], we use the technique of ‘A-caloric approximation’ to establish the optimal partial regularity of nonlinear parabolic systems (1.9). In fact, the use of the ‘A-caloric approximation lemma’ allows optimal regularity, without the use of Reverse-Hölder inequalities and (parabolic) Gehring’s lemma. The method is based on an approximation result that we called the ‘A-caloric approximation lemma’. This is the parabolic analogue of the classical harmonic approximation lemma of De Giorgi [11, 12] and allows to approximate functions with solutions to parabolic systems with constant coefficients in a similar way as the classical harmonic approximation lemma does with harmonic functions. And we can obtain the following theorem.

Theorem 1.1 Let u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) be a weak solution to system (1.9) under the assumptions (H1)-(H4) and the natural growth condition (H5) and denote by Q 0 the set of regularity points of u in Q T :

Q 0 = { z Q T : D u C β , β / 2 ( O , R n N ) , O Q T is a neighborhood of z } .

Then Q 0 is an open subset with full measure, and therefore

Du C β , β / 2 ( Q 0 , R n N ) ,| Q T Q 0 |=0.

At the end of the section, we summarize some notions which we will be used in this paper. For x 0 R n , t 0 R, we denote B( x 0 ,R)={x R n :|x x 0 |<R}, Q(( x 0 , t 0 ),R)=B( x 0 ,R)×( t 0 R 2 , t 0 ). If v is an integrable function in Q( z 0 ,ρ)= Q ρ ( z 0 )= B ρ ( x 0 )×( t 0 ρ 2 , t 0 ), z 0 =( x 0 , t 0 ), we will denote its average by , where α n denotes the volume of the unit ball in R n . We remark that in the following, when not crucial, the ‘center’ of the cylinder will be often unspecified, e.g., Q ρ ( z 0 )= Q ρ ; the same convention will be adopted for balls in R n therefore denoting B( x 0 ,ρ)= B ρ ( x 0 ). Finally, in the rest of the paper, the symbol C will denote a positive, finite constant that may vary from line to line; the relevant dependencies will be specified.

2 The A-caloric approximation technique and preliminaries

In this section we introduce the A-caloric approximation lemma [4] and some preliminaries. Recall a strongly elliptic bilinear form A i α on R n N with an ellipticity constant λ>0, and upper bound Λ>0 means that λ | p ˜ | 2 A i α ( p ˜ , p ˜ ), A i α (p, p ˜ )Λ|p|| p ˜ |, p, p ˜ R n N , we define A-caloric approximation function.

Definition 2.1 We shall say that a function h L 2 (1,0; W 1 , 2 ( B ρ , R N )) is A-caloric on Q ρ if it satisfies

Q ρ ( h i φ t i A i α ( D h , D α φ i ) ) dz=0for all φ C 0 ( Q ρ , R N ) .

Remark 2.2 Obviously, when A( p ˜ , p ˜ ) | p ˜ | 2 for every p ˜ R n N , then an A-caloric function is just a caloric function h t h0.

Lemma 2.3 (A-caloric approximation lemma)

There exists a positive function δ(n,N,λ,Λ,ε)1 with the following property: Whenever A is a bilinear form on R n N , which is strongly ellipticity constant λ>0 and upper bound Λ, ε is a positive number, and u L 2 (1,0; W 1 , 2 (B, R N )) with

Q ( | u | 2 + | D u | 2 ) dz1,
(2.1)

is approximatively A-caloric in the sense that

| Q ( u φ t A ( D u , D φ ) ) dz|δ sup Q |Dφ| for all φ C 0 ( Q , R N ) ,
(2.2)

then there exists an A-caloric function h such that

Q ( | h | 2 + | D h | 2 ) dz1,and Q | u h | 2 dzε.
(2.3)

Actually, we could have directly applied Theorem 5 of [13] with the choice X= W 1 , 2 (B, R N ), B= L 2 (B, R N ), R= W l , 2 (B, R N ), F= ( v k ) k N , p=2 to conclude that ( v k ) k N is relatively compact in L 2 ( Q T , R N )= L 2 (1,0; L 2 (B, R N )).

Lemma 2.4 There exists a positive function δ(n,N,λ,Λ,ε)1 with the following property: Whenever A is a bilinear form on R n N which is strongly ellipticity constant λ>0 and upper bound Λ, ε is a positive number, and u L 2 ( t 0 ρ 2 , t 0 ; W 1 , 2 ( B ρ ( x 0 ), R N )) with

(2.4)

is approximatively A-caloric in the sense that

(2.5)

then there exists h L 2 ( t 0 ρ 2 , t 0 ; W 1 , 2 ( B ρ ( x 0 ), R N )) A-caloric on Q ρ ( z 0 ) such that

(2.6)

For u L 2 ( Q ρ ( z 0 ), R N ) we denote by l z 0 , ρ the unique affine function (in space) l(z)=l(x) minimizing , amongst all affine functions a(z)=a(x) which are independent of t. To get an explicit formula for l z 0 , ρ , we note that such a unique minimum point exists and takes the form l z 0 , ρ (x)= ξ z 0 , ρ + ν z 0 , ρ (x x 0 ), where ν z 0 , ρ R n N . A straightforward computation yields that , for any affine function a(x)=ξ+ν(x x 0 ) with ξ R N and ν R n N . This implies in particular that and .

For convenience we recall from [14] the following.

Lemma 2.5 Let u L 2 ( Q ρ ( z 0 ), R N ), 0<θ<1, and l z 0 , ρ respectively l z 0 , θ ρ the unique affine functions minimizing respectively . Then there holds

Moreover, if Du L 2 ( Q ρ ( z 0 ), R n N ), we have

3 Caccioppoli second inequality

In this section we prove Caccioppoli’s second inequality.

Theorem 3.1 (Caccioppoli second inequality)

Let u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) be a weak solution to (1.9) under the assumptions (H1)-(H4) and the natural growth condition (H5). Then, for any M>0, any affine function l(z)=l(x) independent of t and satisfying |l( z 0 )|+|Dl|M, and any Q ρ ( z 0 ) Q T with 0<ρ<R1, we have

Proof We take the test function φ= η 2 ξ 2 (ul), where η(x) C 0 1 ( B R ( x 0 )) is a cut-off function in space such that 0η1, η1 in B ρ ( x 0 ), |Dη| 1 ( R ρ ) . While ξ C 1 (R) is a cut-off function in time such that, with 0<σ<ρ being arbitrary,

{ ξ 1 , on  ( t 0 ρ 2 , t 0 σ 2 ) , ξ 0 , on  ( , t 0 R 2 ) ( t 0 , ) , 0 ξ 1 , on  R , ξ t 0 , on  ( t 0 ρ 2 , ) , | ξ t | 1 | R ρ | 2 , on  ( t 0 R 2 , t 0 ρ 2 ) .

Thus, we obtain

Q R ( z 0 ) A i α ( z , u , D u ) D ( u l ) i ξ 2 η 2 d z = 2 Q R ( z 0 ) A i α ( z , u , D u ) ξ 2 η η ( u l ) i d z + Q R ( z 0 ) u i t φ i d z + Q R ( z 0 ) B i ( z , u , D u ) φ i d z .

We further have

Q R ( z 0 ) A i α ( z , u , D l ) D α ( u l ) i ξ 2 η 2 d z = 2 Q R ( z 0 ) A i α ( z , u , D l ) ξ 2 η η ( u l ) i d z Q R ( z 0 ) A i α ( z , u , D l ) D α φ i d z

and

0= Q R ( z 0 ) A i α ( z 0 , l ( z 0 ) , D l ) D α φ i dz.

Adding these equations and using l t 0, we deduce

Q R ( z 0 ) [ A i α ( z , u , D u ) A i α ( z , u , D l ) ] D ( u l ) ξ 2 η 2 d z = 2 Q R ( z 0 ) [ A i α ( z , u , D u ) A i α ( z , u , D l ) ] ξ 2 η η ( u l ) d z Q R ( z 0 ) [ A i α ( z , u , D l ) A i α ( z , l , D l ) ] D α φ i d z Q R ( z 0 ) [ A i α ( z , l , D l ) A i α ( z 0 , l ( z 0 ) , D l ) ] D α φ i d z + Q R ( z 0 ) ( u l ) i t φ i d z + Q R ( z 0 ) B i ( z , u , D u ) φ i d z I + II + III + IV + V .
(3.1)

By (1.10) and Young’s inequality, we have

I ε Q R ( z 0 ) ( 1 + | D l | 2 ) m 2 2 | D u D l | 2 ξ 2 η 2 d z + ε ( m 1 ) C m Q R ( z 0 ) ξ m | η | m | u l | m d z + C 2 ε Q R ( z 0 ) ( 1 + | D l | 2 ) m 2 2 ξ 2 | η | 2 | u l | 2 d z + ε Q R ( z 0 ) | D u D l | m ξ m m 1 η m m 1 d z .
(3.2)

By the condition (H4′) and Young’s inequality, we can get

II ε Q R ( z 0 ) ξ 2 η 2 | D u D l | 2 d z + ( 1 ε + 1 ) 1 | R ρ | 2 Q R ( z 0 ) ξ 2 η | u l | 2 d z + ( 1 ε + 2 2 1 β ) [ K ( | l | ) ( 1 + | D l | ) m 2 ] 2 1 β α n R n + 2 + 2 β 1 β .
(3.3)

Similarly, we can estimate III as follows:

III ε Q R ( z 0 ) ξ 2 η 2 | D u D l | 2 d z + Q R ( z 0 ) ξ 2 | η | 2 | u l | 2 d z + ( 1 ε + 4 ) [ K ( | l | ) ( 1 + | D l | ) m 2 + β ] 2 α n R n + 2 + 2 β .
(3.4)

Using the fact that ξ0 on (, t 0 R 2 )( t 0 ,), taking into account that ξ ξ t 0 for t> t 0 ρ 2 and | ξ t | 1 | R ρ | 2 , we infer

IV = Q R ( z 0 ) ( u l ) i t φ i d z = Q R ( z 0 ) | u l | 2 η 2 t ( ξ 2 ) d z + 1 2 Q R ( z 0 ) ξ 2 η 2 t | u l | 2 d z = 1 2 Q R ( z 0 ) | u l | 2 η 2 t ( ξ 2 ) d z = Q R ( z 0 ) | u l | 2 η 2 ξ ξ t d z 1 | R ρ | 2 Q R ( z 0 ) | u l | 2 d z ,
(3.5)

and for μ positive to be fixed later, we have

V = Q R ( z 0 ) a | D u | m ξ 2 η 2 | u l | d z + Q R ( z 0 ) ( | u l | R ρ ξ η ) ( ξ η b ( R ρ ) ) d z Q R ( z 0 ) a [ ( 1 + μ ) | D u D l | m + ( 1 + 1 μ ) | D l | m ] ξ 2 η 2 | u l | d z + 1 2 ε Q R ( z 0 ) ( | u l | R ρ ξ η ) 2 d z + ε 2 Q R ( z 0 ) ξ 2 η 2 b 2 R 2 d z a V ( 1 + μ ) Q R ( z 0 ) ξ 2 η 2 | D u D l | m d z + 1 ε Q R ( z 0 ) ( | u l | R ρ ξ η ) 2 d z + ε 2 [ a 2 ( 1 + 1 μ ) 2 | D l | 2 m + b 2 ] α n R n + 4 .
(3.6)

By (1.11) we have

Q R ( z 0 ) [ A i α ( z , u , D u ) A i α ( z , u , D l ) ] D α ( u l ) i ξ 2 η 2 d z λ Q R ( z 0 ) ( 1 + | D u | 2 + | D l | 2 ) m 2 2 | D u D l | 2 ξ 2 η 2 d z λ Q R ( z 0 ) [ ( 1 + | D l | 2 ) m 2 2 | D u D l | 2 ξ 2 η 2 + | D u D l | m ξ 2 η 2 ] d z .
(3.7)

Combining (3.2)-(3.7) in (3.1) and noting that R 2 β 1 β R 2 β (R1), that 2 2 1 β >4, that [ K ( | l | ) ( 1 + | D l | ) m 2 + β ] 2 [ K ( | l | ) ( 1 + | D l | ) m 2 ] 2 1 β (for K1), choosing ε sufficiently small and taking into account that 2aVλ, that ξ1 for t[ t 0 ρ 2 , t 0 σ 2 ], that η1 on B ρ ( x 0 ), we infer that

t 0 ρ 2 t 0 σ 2 B ρ ( x 0 ) [ ( 1 + | D l | 2 ) m 2 2 | D u D l | 2 + | D u D l | m ] d x d t C 1 [ ( 1 + | D l | 2 ) m 2 2 Q R ( z 0 ) 1 | R ρ | 2 | u l | 2 d z + Q R ( z 0 ) 1 | R ρ | m | u l | m d z ] + C 2 [ K ( | l | ) ( 1 + | D l | ) m 2 ] 2 1 β α n R n + 2 + 2 β + C 3 [ a 2 ( 1 + 1 μ ) 2 | D l | 2 m + b 2 ] α n R n + 4 .

Then the desired result follows by taking the limit σ0. □

4 The proof of the main theorem

The next lemma is a prerequisite for applying the A-caloric approximation technique.

Lemma 4.1 Let u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) be a weak solution to (1.9) under the assumptions (H1)-(H6). Then for any M>0, we have

for any Q ρ ( z 0 ) Q T and φ C 0 ( Q ρ ( z 0 ), R N ) with ρ1 and any affine function l(z)=l(x) independent of time, satisfying |l( z 0 )|+|Dl|M. Here C Eu = C Eu (M,L,m) and we write

Proof Without loss of generality, we can assume that sup Q ρ ( z 0 ) |Dφ|1. From (1.12) and the fact that and , we deduce

In turn, we split the first integral as follows:

I= 1 | Q ρ ( z 0 ) | s 1 ()dz+ 1 | Q ρ ( z 0 ) | s 2 ()dz= I 1 + I 2 ,

and s 1 = Q ρ ( z 0 ){z:|DuDl|1}, s 2 = Q ρ ( z 0 ){z:|DuDl|>1}.

We proceed estimating the two resulting pieces. As for I 1 , using (H6), the fact that s ω m (t,s) is concave and Jensen’s inequality (note that m 1 m > 1 2 ), we get

To estimate I 2 , we preliminarily observe that, using Hölder inequality,

| s 2 | s 2 | D u D l | d z ( s 2 d z ) 1 2 ( s 2 | D u D l | 2 d z ) 1 2 | s 2 | ( Q ρ | D u D l | 2 d z ) 1 2 ,

and therefore

Similarly, we also have

Using (H1), (H2) and the previous inequality, we then conclude the estimate of I 2 as follows:

Combining the estimates found for I 1 and I 2 , we have

|I|L ( 1 + M 2 ) m 2 2 ω(M+1,Φ) Φ +2L ( 1 + M ) m 2 Φ.

For the remaining pieces, using (H4′), we deduce

Here we have used that K1 and the assumption that ρ1. Using again (H4′) and Young’s inequality, we estimate

and

Noting the definition of H and combining the estimates just found for I, II, III and IV, we obtain

A simple scaling argument yields the result for general φ. □

The next lemma is a standard estimate for weak solutions to linear parabolic systems with constant coefficients [15], Lemma 5.1.

Lemma 4.2 Let h L 2 ( t 0 ρ 2 , t 0 ; W 1 , 2 ( B ρ ( x 0 ), R N )) be a weak solution in Q ρ ( z 0 )= B ρ ( x 0 )×( t 0 ρ 2 , t 0 ) of the following linear parabolic system with constant coefficients:

where the coefficients A i α satisfy A i α (p,p)λ | p | 2 , A i α (p, p ˜ )L|p|| p ˜ | for any p, p ˜ R n N . Then h is smooth in Q ρ ( z 0 ) and there exists a constant C pa = C pa (n,N,L/λ)1 such that

ψ ˜ ( z 0 ,θρ) C pa θ 2 ψ ˜ ( z 0 ,ρ),0<θ<1.

Here we write

In the following we consider a weak solution u of the nonlinear parabolic system (1.9) on a fixed sub-cylinder Q ρ ( z 0 ) Q T and ρ1.

Lemma 4.3 Given M>0 and 0<β<α<1, there exist θ(0, 1 2 ) and δ(0,1] depending only on n, N, λ, L, β, α and m such that if

ω ( M + 1 , Ψ ˜ ( z 0 , ρ , l z 0 , ρ ) ) + Ψ ˜ ( z 0 , ρ , l z 0 , ρ ) δ 2 ,

on Q ρ ( z 0 ) Q T for some 0<ρ1 and such if

| l z 0 , ρ ( z 0 ) | +|D l z 0 , ρ |M,

then

Ψ ˜ ( z 0 ,θρ, l z 0 , θ ρ ) θ 2 α Ψ ˜ ( z 0 ,ρ, l z 0 , ρ )+ C 6 ρ 2 β H 2 (M)

for

Ψ ˜ ( z 0 ,ρ, l z 0 , ρ )=Ψ( z 0 ,ρ, l z 0 , ρ )+ H 2 (M) ρ 2 β .

Proof Given M>1. And we shall always consider ρ1. We first want to apply Lemma 4.1 on Q ρ / 2 ( z 0 ) to ul, where l(z)=l(x) is an affine function independent of t satisfying |l( z 0 )|+|Dl|M. We observe that Ψ has the following property:

(4.1)

From Caccioppoli’s second inequality, we infer

Φ( z 0 ,ρ/2,l) C cac [ 2 m Ψ ( z 0 , ρ , l ) + 2 H ( M ) ρ 2 β ] = C ˜ cac Ψ ˜ ( z 0 ,ρ,l).
(4.2)

From Lemma 4.1 we therefore get, for any φ C 0 ( Q ρ / 2 ( z 0 ), R N ), that

(4.3)

where C ˜ Eu = C ˜ Eu (L,M,m).

For given ε>0 to be specified later, we let δ=δ(n,N,λ,L,ε)(0,1] to be constant from Lemma 2.3. Define γ= C ˜ Eu Ψ ( z 0 , ρ ) + 4 δ 2 H 2 ( M ) ρ 2 β and w= γ 1 (ul).

Then from (4.3) we deduce that, for all φ C 0 ( Q ρ / 2 ( z 0 ), R N ), the following holds:

(4.4)

Moreover, we estimate, using Caccioppoli’s second inequality, (4.1) and (4.2),

(4.5)

provided we have chosen C ˜ Eu 1 large enough.

Assuming the smallness condition,

ω ( M + 1 , Ψ ˜ ( z 0 , ρ , l z 0 , ρ ) ) + Ψ ˜ ( z 0 , ρ , l z 0 , ρ ) δ 2 ,
(4.6)

satisfied. Then (4.4) and (4.5) allow us to apply Lemma 2.4, i.e., they yield the existence of h L 2 ( t 0 ρ 2 , t 0 ; W 1 , 2 ( B ρ ( x 0 ), R N )) solving the A i α p β j -heat equation on Q ρ / 2 ( z 0 ) and satisfying

(4.7)

and

(4.8)

From Lemma 4.2 we recall that h satisfies, for any 0<θ<1, the a priori estimate (note that C pa = C pa (n,N,λ,L)1)

Here we have used that , and and (4.7). Combining the previous estimate with (4.8), we deduce

(4.9)

Recalling back (ul) via w= u l γ , we arrive at

(4.10)

Next we use the minimizing property of l z 0 , θ ρ / 2

(4.11)

At the same time, from (4.11), we can see that: For 2mn+2 (n3), we have 2<m< m , where

m ={ m ( n + 2 ) n m + 2 if  n + 2 > m , m > m if  m = n + 2

with 1 m < 1 m < 1 2 . Therefore we can find s[0,1] such that 1 m = 1 s 2 + s m .

Using Sobolev’s, Caccioppoli’s and Young’s inequalities together with (4.11), we have

(4.12)

Using Lemma 2.5, Caccioppoli’s inequality, (4.4), (4.6), (4.12) and Young’s inequality, we obtain

(4.13)

From (4.12) and (4.13), we conclude

(4.14)

provided γ 2 ( m m ) / m θ 2 + ( n + m ) m / m and we fixed ε= θ n + 6 . That it is to say,

(4.15)

Combining (4.11) and (4.15) yields the desired estimate

Ψ( z 0 ,θρ/2, l z 0 , θ ρ / 2 ) C 5 θ 2 ( Ψ ( z 0 , ρ , l ) + 4 δ 2 H 2 ( M ) ρ 2 β )
(4.16)

for C 5 = C 4 +12 C pa . Given β<α<1, we choose 0<θ<1 such that 2 2 α C 5 θ 2 θ 2 α with θ=θ(n,m,N,λ,L,α,β). This also fixes the constants ε=ε(n,m,N,λ,L,α,β) and δ=δ(n,m,N,λ,L,α,β)(0,1]. Thus we have shown Lemma 4.3. □

In the following, we want to iterate Lemma 4.3. That is,

Lemma 4.4 For M>1 and Q ρ ( z 0 ) Q T , suppose that the conditions

( i ) | l z 0 , ρ | + | ( D l ) z 0 , ρ | M ; ( ii ) ρ ρ 0 ( M ) ; ( iii ) Ψ ˜ ( ρ ) Ψ ˜ 0 ( M )

are satisfied. Then, for every jN, we have

Ψ ˜ ( z 0 , θ j ρ , l z 0 , θ ρ ) θ 2 α j Ψ ˜ ( z 0 ,ρ, l z 0 , ρ )+ C 6 (M) ( θ j ρ ) 2 β H 2 (M)

and

| l z 0 , θ j ρ |+ | ( D l ) z 0 , θ j ρ | 2M.

Moreover, the limit

Γ z 0 = lim j ( D u ) z 0 , θ j ρ / 2

exists, and the estimate

is valid for a constant C=C(n,N,λ,α,L,β,M,m).

Proof For fixed z 0 we shall denote l z 0 , ρ l ρ . For given M>1 (and β<α<1), we determine δ=δ(2M), θ=θ(2M) and C 6 = C 6 (2M) according to Lemma 4.3. Then we can find Ψ ˜ 0 (M)>0 sufficiently small such that

ω ( M + 1 , 2 Ψ ˜ 0 ( M ) ) + Ψ ˜ 0 ( M ) δ 2
(4.17)

and

Ψ ˜ 0 (M) M 2 θ n + 4 ( 1 θ α ) 2 4 ( n + 2 ) 2 .
(4.18)

Given this, we can also find ρ 0 (M)(0,1] so small that, writing

C 7 (M)= C 6 ( 2 M ) θ 2 β θ 2 α ,

we have

C 7 (M) ρ 0 ( M ) 2 β H 2 (M)min { δ 2 16 , Ψ ˜ 0 ( M ) , M 2 θ n + 4 ( 1 θ β ) 2 4 ( n + 2 ) 2 } .
(4.19)

Now, suppose that the conditions (i), (ii) and (iii) are satisfied on Q ρ ( z 0 ) Q T . Then, for j=1,2,3, , we shall show

( I ) j Ψ ˜ ( z 0 , θ j ρ , l z 0 , θ ρ ) θ 2 α j Ψ ˜ ( z 0 , ρ , l z 0 , ρ ) + C 7 ( M ) ( θ j ρ ) 2 β H 2 ( M ) , ( II ) j | l z 0 , θ j ρ ( z 0 ) | + | ( D l ) z 0 , θ j ρ | 2 M .

Note first that ( I ) j combined with (ii), (iii) and (4.19) yields

( I ) j Ψ ˜ ( θ j ρ ) 2 Ψ ˜ 0 (M).

Moreover, we have ρ ρ 0 (M)1 and | l z 0 , ρ |+| ( D l ) z 0 , ρ |M. There we can apply Lemma 4.3 to conclude that ( I ) 1 holds. Furthermore, using Lemma 2.5, (iii) and (4.18), we deduce

i.e., ( II ) 1 holds. We now assume that ( I ) ι and ( II ) ι for ι=1,2,,j1 hold. We can apply Lemma 4.3 to calculate

Ψ ˜ ( θ j ρ ) θ 2 α j Ψ ˜ ( ρ ) + C 6 ( 2 M ) ( θ j ρ ) 2 β θ 2 β ι = 0 j 1 θ 2 ( α β ) ι θ 2 α j Ψ ˜ ( ρ ) + C 6 ( 2 M ) θ 2 β θ 2 α ( θ j ρ ) 2 β = θ 2 α j Ψ ˜ ( ρ ) + C 7 ( 2 M ) ( θ j ρ ) 2 β ,

showing ( I ) j . To show ( II ) j we estimate

Here we have used in turn Lemma 2.5, the definition of Ψ( θ ι 1 ρ) and ( I ) ι for ι=1,2,,j1.

Since |D l θ j ρ |2M. We are in a position to apply Theorem 3.1. We obtain

Φ ( θ j ρ / 2 , ( D u ) θ j ρ / 2 ) Φ ( θ j ρ , ( D u ) θ j ρ ) C cac ( 2 M ) Ψ ˜ ( θ j ρ ) C cac ( 2 M ) ( θ 2 α j Ψ ˜ ( ρ ) + C 7 ( 2 M ) ( θ j ρ ) 2 β ) .
(4.20)

We now consider 0<rρ/2. We fix kN{0} with θ k + 1 ρ/2<r θ k ρ/2. Then the previous estimate implies

Next, we show that ( ( D u ) θ j ρ / 2 ) j N is a Cauchy sequence in R n N . For K>j we deduce

This proves the claim. Therefore the limit Γ z 0 = lim j ( D u ) θ j ρ / 2 R n N exists and from the previous estimate, we infer (taking the limit k)

| ( D u ) θ j ρ / 2 Γ z 0 | C 8 (M) θ 2 α j Ψ ˜ ( ρ ) + ( θ j ρ ) 2 β .

Combining this with (4.20), we arrive at

For 0<rρ/2, we find kN{0} with θ k + 1 ρ/2<r< θ k ρ/2. Then the previous estimate implies

This proves the assertion of the lemma. □

An immediate consequence of the previous lemma and of isomorphism theorem of Campanato-Da Prato [16] is the following result.

Theorem 4.1 (Description of regularity points)

Let u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) be a weak solution to the system (1.9) under the assumptions (H1)-(H3) and (H4′), (H5), and denote by Σ the singular set of u. Then Σ Σ 0 Σ 2 , where

and

Σ 2 = { z 0 Q T : lim sup ρ 0 ( | ( u ) z 0 , ρ | + | ( D u ) z 0 , ρ | ) = } .

At last, we have the following.

Theorem 4.2 (Almost everywhere regularity)

Let u L m (T,0; W 1 , m (Ω, R N )) L ( Q T ; R N ) be a weak solution to the system (1.9) under the assumptions (H1)-(H3) and (H4), (H5), and denote by Σ the singular set of u. Then Σ Σ 1 Σ 2 , where Σ 2 is as in Theorem  4.1 and

Proof We start taking a point z 0 ( x 0 , t 0 ) Q T such that

(4.21)

and

sup ρ > 0 | ( u ) z 0 , ρ | + sup ρ > 0 | ( D u ) z 0 , ρ | M<.
(4.22)

The proof is complete if we show that such points are regularity points.

Step 1: a comparison estimate. Consider the unique weak solution v L m ( t 0 4 ρ 2 , t 0 ; W 1 , m ( B ρ ( x 0 ), R N )) of the initial boundary value problem

{ Q 2 ρ ( z 0 ) ( v i φ t i A i α ( z 0 , ( u ) z 0 , ρ , D v ) D α φ i ) d z = 0 , φ C 0 ( Q 2 ρ ( z 0 ) , R N ) , v = u , on  B 2 ρ ( x 0 ) × { t 0 4 ρ 2 } B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , t 0 ) .

Then the difference uv satisfies

Q 2 ρ ( z 0 ) [ ( u v ) i φ t i ( A i α ( z , u , D u ) A i α ( z 0 , ( u ) z 0 , ρ , D v ) ) D α φ i ] d z + Q 2 ρ ( z 0 ) B i ( z , u , D u ) φ i d z = 0

for every φ C 0 ( Q 2 ρ ( z 0 ), R N ). We now choose φ=χ(t) ( u v ) i with χ1 for (,s), χ0 on (s+ε,), and χ(t)=(s+εt)/ε for sts+ε, where [s,s+ε]( t 0 4 ρ 2 , t 0 ). Then

1 2 Q 2 ρ ( z 0 ) t ( | u v | 2 χ ) d z + 1 2 Q 2 ρ ( z 0 ) | u v | 2 t χ d z Q 2 ρ ( z 0 ) ( A i α ( z , u , D u ) A i α ( z 0 , ( u ) z 0 , ρ , D v ) ) ( D u i D v i ) χ d z + Q 2 ρ ( z 0 ) B i ( z , u , D u ) χ ( t ) ( u v ) i d z = 0 .

Letting ε0, we easily obtain that for a.e. s( t 0 4 ρ 2 , t 0 )

1 2 u ( , s ) v ( , s ) L 2 ( B 2 ρ ( x 0 ) ) 2 + B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , t 0 ) ( A i α ( z 0 , ( u ) z 0 , ρ , D u ) A i α ( z 0 , ( u ) z 0 , ρ , D v ) ) D ( u v ) i χ d z = B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) ( A i α ( z 0 , ( u ) z 0 , ρ , D u ) A i α ( z , u , D u ) ) D ( u v ) i d z + B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) B i ( z , u , D u ) ( u v ) i d z = 0 .

The second term of the left-hand side of the previous equation can be estimated by the use of monotonicity, i.e., (H3). We therefore obtain

1 2 u ( , s ) v ( , s ) L 2 ( B 2 ρ ( x 0 ) ) 2 + λ ( 1 + | D v | 2 + | D u | 2 ) m 2 2 B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) | D u D v | 2 d z B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) ( A i α ( z 0 , ( u ) z 0 , ρ , D u ) A i α ( z , u , D u ) ) D ( u v ) i d z + B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) B i ( z , u , D u ) ( u v ) i d z = I + II .
(4.23)

To estimate the right-hand side, we use (H4) which easily yields

| A i α ( z 0 , ( u ) z 0 , ρ , D u ) A i α ( z , u , D u ) | L θ ( 2 | ( u ) z 0 , 2 ρ | + | u ( u ) z 0 , 2 ρ | , 4 ρ + | u ( u ) z 0 , 2 ρ | ) ( 1 + | D u | ) m 2 .

Using the previous estimate, Young’s inequality and the fact that θ1, we have

| I | λ 2 ( 1 + | D v | 2 + | D u | 2 ) m 2 2 B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) | D u D v | 2 d z + 2 L 2 λ Q 2 ρ ( z 0 ) θ ( 2 | ( u ) z 0 , 2 ρ | + | u ( u ) z 0 , 2 ρ | , 4 ρ + | u ( u ) z 0 , 2 ρ | ) ( 1 + | D u | m ) d z .

Having combined the previous estimate with (4.23), we arrive at

1 2 u ( , s ) v ( , s ) L 2 ( B 2 ρ ( x 0 ) ) 2 + λ 2 ( 1 + | D v | 2 + | D u | 2 ) m 2 2 B 2 ρ ( x 0 ) × ( t 0 4 ρ 2 , s ) | D u D v | 2 d z 2 L 2 λ Q 2 ρ ( z 0 ) θ ( 2 | ( u ) z 0 , 2 ρ | + | u ( u ) z 0 , 2 ρ | , 4 ρ + | ( u ) z 0 , 2 ρ | ) ( 1 + | D u | m ) d z + II = 2 L 2 λ III + II .
(4.24)

We shall provide on estimate for III. We denote , σ>0.

If we let A t σ A t ={z Q 2 ρ ( z 0 ):|u ( u ) z 0 , 2 ρ |t}, then

| A t | 1 t Q 2 ρ ( z 0 ) | u ( u ) z 0 , 2 ρ | dz | Q 2 ρ | t ε 2 ρ .
(4.25)

We now split III

III= A t ()dz+ Q 2 ρ ( z 0 ) A t ()dz=IV+V

and estimate IV and V. We have, using that θ1, (4.25) and (4.22)

IV 2 m 1 Q 2 ρ ( z 0 ) | D u ( D u ) z 0 , 2 ρ | m d z + ( 1 + 2 m 1 | ( D u ) z 0 , 2 ρ | m ) | A t | 2 m 1 Q 2 ρ ( z 0 ) | D u ( D u ) z 0 , 2 ρ | m d z + 2 m 1 ( 1 + 2 M 2 ) | Q 2 ρ | t ε 2 ρ .

From the definition of θ, we have

V4K(2M+t) ( ρ + t ) β Q 2 ρ ( z 0 ) ( 1 + | D u | m ) dz.

Noting that sup Q T u=V, we have

II 2 m 1 a V Q 2 ρ ( z 0 ) | D u ( D u ) z 0 , 2 ρ | m d z + Q 2 ρ ( z 0 ) | u v | 2 ρ 2 d z + ( 2 m 1 a M m + b ) α n ρ n + 2 ρ 2 .

We now choose the parameter t carefully, i.e., t= ε 2 ρ and let ε suitably small. Then connecting the previous estimates for II, III, IV and V to (4.24), we easily have the estimate we were interested in, that is,

(4.26)

In particular, we see that

(4.27)

We observe that, as a consequence of (4.21) and (4.22), we have that

lim inf ρ 0 S(ρ)=0.
(4.28)

Step 2: A Poincare-type inequality. Let us define

v ˜ =v ( D u ) z 0 , 2 ρ (x x 0 ).

Therefore v ˜ solves

Q ρ ( z 0 ) ( v ˜ i φ t i A ˜ i α ( D v ˜ ) D α φ i ) dz=0,φ C 0 ( Q 2 ρ ( z 0 ) , R N ) ,

where A ˜ i α (p)= A i α ( z 0 , ( u ) z 0 , 2 ρ , ( D v ) z 0 , 2 ρ +p) for every p R n N . From [17], Theorem 3.1, we conclude that v ˜ W 1 , 2 ( t 0 ρ 2 , t 0 ; W 1 , 2 ( B 2 ρ ( x 0 ) , R N )) and that

In view of the previous estimate, using the Poincare inequality for v and (4.26), we find

where C=C(n,λ,L).

Finally, by comparison, we get the Poincare inequality for u via (4.26) and the previous estimate

(4.29)

for a constant C=C(n,λ,L).

Step 3: Conclusion. From the previous estimate and (4.28), the assertion readily follows. Indeed if z 0 Q T satisfies (4.21) and (4.22), then we have

therefore z 0 is a regular point in view of Theorem 4.1. □

References

  1. Wittle M: Computer simulations of an electrorheological fluid. J. Non-Newton. Fluid Mech. 1993, 37: 233–263.

    Article  Google Scholar 

  2. Rajagopal KR, Ruzicka M: On the modeling of electrorheological materials. Mech. Res. Commun. 1996, 23: 401–407. 10.1016/0093-6413(96)00038-9

    Article  MATH  Google Scholar 

  3. Rajagopal KR, Ruzicka M: Mathematical modeling of electrorheological materials. Contin. Mech. Thermodyn. 2001, 13: 59–78. 10.1007/s001610100034

    Article  MATH  Google Scholar 

  4. Duzaar F, Mingione G: Second order parabolic systems, optimal regularity, and singular sets of solutions. Ann. Inst. Henri Poincaré, Anal. Non Linéaire 2005, 22: 705–751. 10.1016/j.anihpc.2004.10.011

    Article  MathSciNet  MATH  Google Scholar 

  5. Duzaar F, Grotowski JF: Optimal interior partial regularity for nonlinear elliptic systems: the method of A -harmonic approximation. Manuscr. Math. 2000, 103: 267–298. 10.1007/s002290070007

    Article  MathSciNet  MATH  Google Scholar 

  6. Chen S, Tan Z: The method of p -harmonic approximation and optimal interior partial regularity for energy minimizing p -harmonic maps under the controllable growth condition. Sci. China Ser. A 2007, 50(1):105–115.

    Article  MathSciNet  MATH  Google Scholar 

  7. Chen S, Tan Z: The method of A -harmonic approximation and optimal interior partial regularity for nonlinear elliptic systems under the controllable growth condition. J. Math. Anal. Appl. 2007, 335: 20–42. 10.1016/j.jmaa.2007.01.042

    Article  MathSciNet  MATH  Google Scholar 

  8. Chen S, Tan Z: Optimal interior partial regularity for nonlinear elliptic systems. Discrete Contin. Dyn. Syst., Ser. A 2010, 27(3):981–993.

    Article  MathSciNet  MATH  Google Scholar 

  9. Chen S, Tan Z:Optimal interior partial regularity for nonlinear elliptic systems for the case 1<m<2 under natural growth condition. J. Inequal. Appl. 2010., 2010: Article ID 680714

    Google Scholar 

  10. Zhong T: C 1 , α -partial regularity of nonlinear parabolic systems. J. Partial Differ. Equ. 1992, 5: 23–34.

    MathSciNet  Google Scholar 

  11. De Giorgi, E: Frontiere orientate di misura minima, Sem. Scuola Normale Superiore Pisa (1960–1961)

  12. De Giorgi E: Un esempio di estremali discontinue per un problema variazionale di tipo ellittico. Boll. Unione Mat. Ital. 1968, 1(4):135–137.

    MathSciNet  MATH  Google Scholar 

  13. Simon L:Compact sets in the space L p (0,T;B). Ann. Mat. Pura Appl. 1987, 146(4):65–96.

    MathSciNet  MATH  Google Scholar 

  14. Kronz M: Partial regularity results for minimizers of quasiconvex functionals of higher order. Ann. Inst. Henri Poincaré, Anal. Non Linéaire 2002, 19: 81–112. 10.1016/S0294-1449(01)00072-5

    Article  MathSciNet  MATH  Google Scholar 

  15. Campanato S:Equazioni paraboliche del secondo ordine e spazi L 2 , θ (Ω,δ). Ann. Mat. Pura Appl. 1996, 73(4):55–102.

    MathSciNet  Google Scholar 

  16. Da Prato G:Spazi L 2 , θ (Ω,δ) e loro proprieta. Ann. Mat. Pura Appl. 1965, 69(4):383–392.

    Article  MathSciNet  Google Scholar 

  17. Campanato S: On the nonlinear parabolic systems in divergence form Hölder continuity and partial Hölder continuity of the solutions. Ann. Mat. Pura Appl. 1984, 137(4):83–122.

    Article  MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

Supported by the National Natural Science Foundation of China (Nos: 11201415, 11271305), the Natural Science Foundation of Fujian Province (2012J01027) and the Training Programme Foundation for Excellent Youth Researching Talents of Fujian’s Universities (JA12205).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zhong Tan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SC participated in design of the study and drafted the manuscript. ZT participated in conceived of the study and the amendment of the paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Chen, S., Tan, Z. Optimal partial regularity of second-order parabolic systems under natural growth condition. Bound Value Probl 2013, 152 (2013). https://doi.org/10.1186/1687-2770-2013-152

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-152

Keywords