Research

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

Shuhong Chen1 and Zhong Tan2*

Author Affiliations

1 School of Mathematics and Statistics, Minnan Normal University, Zhangzhou, Fujian, 363000, China

2 School of Mathematical Science, Xiamen University, Xiamen, Fujian, 361005, China

For all author emails, please log on.

Boundary Value Problems 2013, 2013:152  doi:10.1186/1687-2770-2013-152

 Received: 17 August 2012 Accepted: 10 May 2013 Published: 26 June 2013

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

### 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.

##### Keywords:
nonlinear parabolic systems; natural growth condition; A-caloric approximation; optimal partial regularity

### 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)

curl E = 0 , (1.2)

ρ 0 v t div S + ρ 0 [ v ] v + ϕ = ρ 0 f + [ E ] P , (1.3)

div v = 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 , D u ) = B i ( z , u , D u ) , i = 1 , 2 , , N , (1.9)

where Ω R n is a bounded domain and T > 0 , z = ( x , t ) with x Ω , 0 < t T , 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., D u = { 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 1 L < . 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 , 2 a V < λ .

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 | p q | , (1.10)

( A i α ( z , u , p ) A i α ( z , u , q ) ) ( p q ) λ ( 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 ) d z = Q T B i ( z , u , D u ) φ i d z (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 m 2 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 [6-9] 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.1Let 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 ofuin 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

D u 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 ) ) d z = 0 for 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 h 0 .

Lemma 2.3 (A-caloric approximation lemma)

There exists a positive function δ ( n , N , λ , Λ , ε ) 1 with the following property: WheneverAis 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 ) d z 1 , (2.1)

is approximativelyA-caloric in the sense that

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

then there exists anA-caloric functionhsuch that

Q ( | h | 2 + | D h | 2 ) d z 1 , and Q | u h | 2 d z ε . (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.4There exists a positive function δ ( n , N , λ , Λ , ε ) 1 with the following property: WheneverAis 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 approximativelyA-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.5Let u L 2 ( Q ρ ( z 0 ) , R N ) , 0 < θ < 1 , and l z 0 , ρ respectively l z 0 , θ ρ the unique affine functions minimizingrespectively. Then there holds

Moreover, if D u 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 oftand satisfying | l ( z 0 ) | + | D l | M , and any Q ρ ( z 0 ) Q T with 0 < ρ < R 1 , we have

Proof We take the test function φ = η 2 ξ 2 ( u l ) , 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 d z .

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 β ( R 1 ), that 2 2 1 β > 4 , that [ K ( | l | ) ( 1 + | D l | ) m 2 + β ] 2 [ K ( | l | ) ( 1 + | D l | ) m 2 ] 2 1 β (for K 1 ), choosing ε sufficiently small and taking into account that 2 a V λ , 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.1Let 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 ) | + | D l | 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 ( ) d z + 1 | Q ρ ( z 0 ) | s 2 ( ) d z = I 1 + I 2 ,

and s 1 = Q ρ ( z 0 ) { z : | D u D l | 1 } , s 2 = Q ρ ( z 0 ) { z : | D u D l | > 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 , Φ ) Φ + 2 L ( 1 + M ) m 2 Φ .

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

Here we have used that K 1 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.2Let 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 . Thenhis 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.3Given M > 0 and 0 < β < α < 1 , there exist θ ( 0 , 1 2 ) and δ ( 0 , 1 ] depending only onn, N, λ, L, β, αandmsuch 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 u l , where l ( z ) = l ( x ) is an affine function independent of t satisfying | l ( z 0 ) | + | D l | 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 ( u l ) .

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 ( u l ) 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 2 m n + 2 ( n 3 ), 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.4For 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 j N , 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 ρ | 2 M .

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 δ = δ ( 2 M ) , θ = θ ( 2 M ) and C 6 = C 6 ( 2 M ) 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 , , j 1 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 , , j 1 .

Since | D l θ j ρ | 2 M . 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 k N { 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 k N { 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 ofu. 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 ofu. 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 u v 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 s t s + ε , 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 ρ | d z | Q 2 ρ | t ε 2 ρ . (4.25)

We now split III

III = A t ( ) d z + Q 2 ρ ( z 0 ) A t ( ) d z = 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

V 4 K ( 2 M + t ) ( ρ + t ) β Q 2 ρ ( z 0 ) ( 1 + | D u | m ) d z .

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 ) d z = 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. □

### 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.

### 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).

### References

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

2. Rajagopal, KR, Ruzicka, M: On the modeling of electrorheological materials. Mech. Res. Commun.. 23, 401–407 (1996). Publisher Full Text

3. Rajagopal, KR, Ruzicka, M: Mathematical modeling of electrorheological materials. Contin. Mech. Thermodyn.. 13, 59–78 (2001). Publisher Full Text

4. Duzaar, F, Mingione, G: Second order parabolic systems, optimal regularity, and singular sets of solutions. Ann. Inst. Henri Poincaré, Anal. Non Linéaire. 22, 705–751 (2005). Publisher Full Text

5. Duzaar, F, Grotowski, JF: Optimal interior partial regularity for nonlinear elliptic systems: the method of A-harmonic approximation. Manuscr. Math.. 103, 267–298 (2000). Publisher Full Text

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. 50(1), 105–115 (2007)

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.. 335, 20–42 (2007). Publisher Full Text

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

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, Article ID 680714 (2010)

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

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.. 1(4), 135–137 (1968)

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

14. Kronz, M: Partial regularity results for minimizers of quasiconvex functionals of higher order. Ann. Inst. Henri Poincaré, Anal. Non Linéaire. 19, 81–112 (2002). Publisher Full Text

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

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

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.. 137(4), 83–122 (1984)