Home Global regularity for systems with p-structure depending on the symmetric gradient
Article Open Access

Global regularity for systems with p-structure depending on the symmetric gradient

  • Luigi C. Berselli ORCID logo EMAIL logo and Michael Růžička
Published/Copyright: October 11, 2018

Abstract

In this paper we study on smooth bounded domains the global regularity (up to the boundary) for weak solutions to systems having p-structure depending only on the symmetric part of the gradient.

MSC 2010: 76A05; 35D35; 35Q35

1 Introduction

In this paper we study regularity of weak solutions to the boundary value problem

(1.1) { - div 𝐒 ( 𝐃𝐮 ) = 𝐟 in  Ω , 𝐮 = 𝟎 on  Ω ,

where 𝐃𝐮 := 1 2 ( 𝐮 + 𝐮 ) denotes the symmetric part of the gradient 𝐮 and Ω 3 is a bounded domain with a C 2 , 1 boundary Ω .[1] Our interest in this system comes from the p-Stokes system

(1.2) { - div 𝐒 ( 𝐃𝐮 ) + π = 𝐟 in  Ω , div 𝐮 = 0 in  Ω , 𝐮 = 𝟎 on  Ω .

In both problems the typical example for 𝐒 we have in mind is

𝐒 ( 𝐃𝐮 ) = μ ( δ + | 𝐃𝐮 | ) p - 2 𝐃𝐮 ,

where p ( 1 , 2 ] , δ 0 and μ > 0 . In previous investigations of (1.2), only suboptimal results for the regularity up to the boundary have been proved. Here we mean suboptimal in the sense that the results are weaker than the results known for p-Laplacian systems, cf. [1, 13, 14]. Clearly, system (1.1) is obtained from (1.2) by dropping the divergence constraint and the resulting pressure gradient. Thus, system (1.1) lies in between system (1.2) and p-Laplacian systems, which depend on the full gradient 𝐮 .

We would like to stress that system (1.1) is of its own independent interest, since it is studied within plasticity theory, when formulated in the framework of deformation theory (cf. [11, 24]). In this context the unknown is the displacement vector field 𝐮 = ( u 1 , u 2 , u 3 ) , while the external body force 𝐟 = ( f 1 , f 2 , f 3 ) is given. The stress tensor 𝐒 , which is the tensor of small elasto-plastic deformations, depends only on 𝐃𝐮 . Physical interpretation and discussion of both systems (1.1) and (1.2) and the underlying models can be found, e.g., in [5, 11, 15, 19, 20].

We study global regularity properties of weak solutions to (1.1) in sufficiently smooth and bounded domains Ω; we obtain, for all p ( 1 , 2 ] , the optimal result, namely, that 𝐅 ( 𝐃𝐮 ) belongs to W 1 , 2 ( Ω ) , where the nonlinear tensor-valued function 𝐅 is defined in (2.5). This result has been proved near a flat boundary in [24] and is the same result as for p-Laplacian systems (cf. [1, 13, 14]). The situation is quite different for (1.2). There the optimal result, i.e., 𝐅 ( 𝐃𝐮 ) W 1 , 2 ( Ω ) , is only known for

  1. two-dimensional bounded domains (cf. [16], where even the p-Navier–Stokes system is treated);

  2. the space-periodic problem in d , d 2 , which follows immediately from interior estimates, i.e., 𝐅 ( 𝐃𝐮 ) W loc 1 , 2 ( Ω ) , which are known in all dimensions, and the periodicity of the solution;

  3. if the no-slip boundary condition is replaced by perfect slip boundary conditions (cf. [17]);

  4. in the case of small 𝐟 (cf. [6]).

We also observe that the above results for the p-Stokes system (apart those in the space periodic setting) require the stress tensor to be non-degenerate, that is, δ > 0 . In the case of homogeneous Dirichlet boundary conditions and three- and higher-dimensional bounded, sufficiently smooth domains only suboptimal results are known. To our knowledge, the state of the art for general data is that 𝐅 ( 𝐃𝐮 ) W loc 1 , 2 ( Ω ) , tangential derivatives of 𝐅 ( 𝐃𝐮 ) near the boundary belong to L 2 , while the normal derivative of 𝐅 ( 𝐃𝐮 ) near the boundary belongs to some L q , where q = q ( p ) < 2 (cf. [2, 4] and the discussion therein). We would also like to mention a result for another system between (1.2) and p-Laplacian systems, namely, if (1.2) is considered with 𝐒 depending on the full velocity gradient 𝐮 . In this case, it is proved in [7] that 𝐮 W 2 , r ( 3 ) W 0 1 , p ( 3 ) for some r > 3 , provided p < 2 is very close to 2.

In the present paper we extend the optimal regularity result for (1.1) of Seregin and Shilkin [24] in the case of a flat boundary to the general case of bounded sufficiently smooth domains and to possibly degenerate stress tensors, that is, the case δ = 0 . The precise result we prove is the following:

Theorem 1.1.

Let the tensor field S in (1.1) have ( p , δ ) -structure (cf. Definition 2.1) for some p ( 1 , 2 ] and δ [ 0 , ) , and let F be the associated tensor field to S defined in (2.5). Let Ω R 3 be a bounded domain with C 2 , 1 boundary, and let f L p ( Ω ) . Then the unique weak solution u W 0 1 , p ( Ω ) of problem (1.1) satisfies

Ω | 𝐅 ( 𝐃𝐮 ) | 2 d 𝐱 c ,

where c denotes a positive function which is non-decreasing in f p and δ, and which depends on the domain through its measure | Ω | and the C 2 , 1 -norms of the local description of Ω . In particular, the above estimate implies that u W 2 , 3 p p + 1 ( Ω ) .

2 Notations and preliminaries

In this section we introduce the notation we will use, state the precise assumptions on the extra stress tensor 𝐒 , and formulate the main results of the paper.

2.1 Notation

We use c , C to denote generic constants which may change from line to line, but they are independent of the crucial quantities. Moreover, we write f g if and only if there exists constants c , C > 0 such that c f g C f . In some cases we need to specify the dependence on certain parameters, and consequently we denote by c ( ) a positive function which is non-decreasing with respect to all its arguments. In particular, we denote by c ( δ - 1 ) a possibly critical dependence on the parameter δ as δ 0 , while c ( δ ) only indicates that the constant c depends on δ and will satisfy c ( δ ) c ( δ 0 ) for all δ δ 0 .

We use standard Lebesgue spaces ( L p ( Ω ) , p ) and Sobolev spaces ( W k , p ( Ω ) , k , p ) , where Ω 3 is a sufficiently smooth bounded domain. The space W 0 1 , p ( Ω ) is the closure of the compactly supported, smooth functions C 0 ( Ω ) in W 1 , p ( Ω ) . Thanks to the Poincaré inequality, we equip W 0 1 , p ( Ω ) with the gradient norm p . When dealing with functions defined only on some open subset G Ω , we denote the norm in L p ( G ) by p , G . As usual, we use the symbol to denote weak convergence, and to denote strong convergence. The symbol spt f denotes the support of the function f. We do not distinguish between scalar, vector-valued or tensor-valued function spaces. However, we denote vectors by boldface lower-case letter as, e.g., 𝐮 and tensors by boldface upper case letters as, e.g., 𝐒 . For vectors 𝐮 , 𝐯 3 , we denote 𝐮 𝑠 𝐯 := 1 2 ( 𝐮 𝐯 + ( 𝐮 𝐯 ) ) , where the standard tensor product 𝐮 𝐯 3 × 3 is defined as ( 𝐮 𝐯 ) i j := u i v j . The scalar product of vectors is denoted by 𝐮 𝐯 = i = 1 3 u i v i and the scalar product of tensors is denoted by 𝐀 𝐁 := i , j = 1 3 A i j B i j .

Greek lower-case letters take only the values 1 , 2 , while Latin lower-case ones take the values 1 , 2 , 3 . We use the summation convention over repeated indices only for Greek lower-case letters, but not for Latin lower-case ones.

2.2 ( p , δ ) -structure

We now define what it means that a tensor field 𝐒 has ( p , δ ) -structure, see [8, 23]. For a tensor 𝐏 3 × 3 , we denote its symmetric part by 𝐏 sym := 1 2 ( 𝐏 + 𝐏 ) sym 3 × 3 := { 𝐏 3 × 3 𝐏 = 𝐏 } . We use the notation | 𝐏 | 2 = 𝐏 𝐏 .

It is convenient to define for t 0 a special N-function[2] ϕ ( ) = ϕ p , δ ( ) , for p ( 1 , ) , δ 0 , by

(2.1) ϕ ( t ) := 0 t ( δ + s ) p - 2 s 𝑑 s .

The function ϕ satisfies, uniformly in t and independently of δ, the important equivalences

(2.2) ϕ ′′ ( t ) t ϕ ( t ) ,
ϕ ( t ) t ϕ ( t ) ,
(2.3) t p + δ p ϕ ( t ) + δ p .

We use the convention that if ϕ ′′ ( 0 ) does not exist, the left-hand side in (2.2) is continuously extended by zero for t = 0 . We define the shifted N-functions { ϕ a } a 0 (cf. [8, 9, 23]), for t 0 , by

ϕ a ( t ) := 0 t ϕ ( a + s ) s a + s 𝑑 s

Note that the family { ϕ a } a 0 satisfies the Δ 2 -condition uniformly with respect to a 0 , i.e., ϕ a ( 2 t ) c ( p ) ϕ a ( t ) holds for all t 0 .

Definition 2.1 ( ( p , δ ) -structure).

We say that a tensor field 𝐒 : 3 × 3 sym 3 × 3 , belonging to C 0 ( 3 × 3 , sym 3 × 3 ) C 1 ( 3 × 3 { 𝟎 } , sym 3 × 3 ) , satisfying 𝐒 ( 𝐏 ) = 𝐒 ( 𝐏 sym ) and 𝐒 ( 𝟎 ) = 𝟎 , possesses ( p , δ ) -structure if for some p ( 1 , ) , δ [ 0 , ) , and the N-function ϕ = ϕ p , δ (cf. (2.1)), there exist constants κ 0 , κ 1 > 0 such that

(2.4) i , j , k , l = 1 3 k l S i j ( 𝐏 ) Q i j Q k l κ 0 ϕ ′′ ( | 𝐏 sym | ) | 𝐐 sym | 2 , | k l S i j ( 𝐏 ) | κ 1 ϕ ′′ ( | 𝐏 sym | )

for all 𝐏 , 𝐐 3 × 3 , with 𝐏 sym 𝟎 , and all i , j , k , l = 1 , 2 , 3 . The constants κ 0 , κ 1 and p are called the characteristics of 𝐒 .

Remark 2.2.

  1. Assume that 𝐒 has ( p , δ ) -structure for some δ [ 0 , δ 0 ] . Then, if not otherwise stated, the constants in the estimates depend only on the characteristics of 𝐒 and on δ 0 , but they are independent of δ.

  2. An important example of a tensor field 𝐒 having ( p , δ ) -structure is given by 𝐒 ( 𝐏 ) = ϕ ( | 𝐏 sym | ) | 𝐏 sym | - 1 𝐏 sym . In this case, the characteristics of 𝐒 , namely, κ 0 and κ 1 , depend only on p and are independent of δ 0 .

  3. For a tensor field 𝐒 with ( p , δ ) -structure, we have k l S i j ( 𝐏 ) = k l S j i ( 𝐏 ) for all i , j , k , l = 1 , 2 , 3 and all 𝐏 3 × 3 , due to its symmetry. Moreover, from 𝐒 ( 𝐏 ) = 𝐒 ( 𝐏 sym ) , it follows that k l S i j ( 𝐏 ) = 1 2 k l S i j ( 𝐏 sym ) + 1 2 l k S i j ( 𝐏 sym ) for all i , j , k , l = 1 , 2 , 3 and all 𝐏 3 × 3 , and consequently k l S i j ( 𝐏 ) = l k S i j ( 𝐏 ) for all i , j , k , l = 1 , 2 , 3 and all 𝐏 sym 3 × 3 .

To a tensor field 𝐒 with ( p , δ ) -structure, we associate the tensor field 𝐅 : 3 × 3 sym 3 × 3 defined through

(2.5) 𝐅 ( 𝐏 ) := ( δ + | 𝐏 sym | ) p - 2 2 𝐏 sym .

The connection between 𝐒 , 𝐅 , and { ϕ a } a 0 is best explained in the following proposition (cf. [8, 23]).

Proposition 2.3.

Let S have ( p , δ ) -structure, and let F be as defined in (2.5). Then

( 𝐒 ( 𝐏 ) - 𝐒 ( 𝐐 ) ) ( 𝐏 - 𝐐 ) | 𝐅 ( 𝐏 ) - 𝐅 ( 𝐐 ) | 2 ϕ | 𝐏 sym | ( | 𝐏 sym - 𝐐 sym | )
ϕ ′′ ( | 𝐏 sym | + | 𝐏 sym - 𝐐 sym | ) | 𝐏 sym - 𝐐 sym | 2 ,
(2.6) | 𝐒 ( 𝐏 ) - 𝐒 ( 𝐐 ) | ϕ | 𝐏 sym | ( | 𝐏 sym - 𝐐 sym | )

uniformly in P , Q R 3 × 3 . Moreover, uniformly in Q R 3 × 3 ,

𝐒 ( 𝐐 ) 𝐐 | 𝐅 ( 𝐐 ) | 2 ϕ ( | 𝐐 sym | ) .

The constants depend only on the characteristics of S .

For a detailed discussion of the properties of 𝐒 and 𝐅 and their relation to Orlicz spaces and N-functions, we refer the reader to [23, 3]. Since in the following we shall insert into 𝐒 and 𝐅 only symmetric tensors, we can drop in the above formulas the superscript “ sym ” and restrict the admitted tensors to symmetric ones.

We recall that the following equivalence, which is proved in [3, Lemma 3.8], is valid for all smooth enough symmetric tensor fields 𝐐 sym 3 × 3 :

(2.7) | i 𝐅 ( 𝐐 ) | 2 ϕ ′′ ( | 𝐐 | ) | i 𝐐 | 2 .

The proof of this equivalence is based on Proposition 2.3. This proposition and the theory of divided differences also imply (cf. [4, Equation (2.26)]) that

(2.8) | τ 𝐅 ( 𝐐 ) | 2 ϕ ′′ ( | 𝐐 | ) | τ 𝐐 | 2

for all smooth enough symmetric tensor fields 𝐐 sym 3 × 3 .

A crucial observation in [24] is that the quantities in (2.7) are also equivalent to several further quantities. To formulate this precisely, we introduce, for i = 1 , 2 , 3 and for sufficiently smooth symmetric tensor fields 𝐐 , the quantity

(2.9) 𝒫 i ( 𝐐 ) := i 𝐒 ( 𝐐 ) i 𝐐 = k , l , m , n = 1 3 k l S m n ( 𝐐 ) i Q k l i Q m n .

Recall that in the definition of 𝒫 i ( 𝐐 ) there is no summation convention over the repeated Latin lower-case index i in i 𝐒 ( 𝐐 ) i 𝐐 . Note that if 𝐒 has ( p , δ ) -structure, then 𝒫 i ( 𝐐 ) 0 for i = 1 , 2 , 3 . The following important equivalences hold, first proved in [24].

Proposition 2.4.

Assume that S has ( p , δ ) -structure. Then the following equivalences are valid, for all smooth enough symmetric tensor fields Q and all i = 1 , 2 , 3 :

(2.10) 𝒫 i ( 𝐐 ) ϕ ′′ ( | 𝐐 | ) | i 𝐐 | 2 | i 𝐅 ( 𝐐 ) | 2 ,
(2.11) 𝒫 i ( 𝐐 ) | i 𝐒 ( 𝐐 ) | 2 ϕ ′′ ( | 𝐐 | ) ,

with the constants only depending on the characteristics of S .

Proof.

The assertions are proved in [24] using a different notation. For the convenience of the reader, we sketch the proof here. The equivalences in (2.10) follow from (2.7), (2.9) and the fact that 𝐒 has ( p , δ ) -structure. Furthermore, using (2.10), we have

| 𝒫 i ( 𝐐 ) | 2 | i 𝐒 ( 𝐐 ) | 2 | i 𝐐 | 2 c | i 𝐒 ( 𝐐 ) | 2 𝒫 i ( 𝐐 ) ϕ ′′ ( | 𝐐 | ) ,

which proves one inequality of (2.11). The other one follows from

| i 𝐒 ( 𝐐 ) | 2 c k , l = 1 3 | k l 𝐒 ( 𝐐 ) i Q k l | 2 c ( ϕ ′′ ( | 𝐐 | ) ) 2 | i 𝐐 | 2 c ϕ ′′ ( | 𝐐 | ) 𝒫 i ( 𝐐 ) ,

where we used (2.4) and (2.10). ∎

2.3 Existence of weak solutions

In this section we define weak solutions of (1.1), recall the main results of existence and uniqueness and discuss a perturbed problem, which is used to justify the computations that follow. From now on, we restrict ourselves to the case p 2 .

Definition 2.5.

We say that 𝐮 W 0 1 , p ( Ω ) is a weak solution to (1.1) if for all 𝐯 W 0 1 , p ( Ω ) ,

Ω 𝐒 ( 𝐃𝐮 ) 𝐃𝐯 𝑑 𝐱 = Ω 𝐟 𝐯 𝑑 𝐱 .

We have the following standard result.

Proposition 2.6.

Let the tensor field S in (1.1) have ( p , δ ) -structure for some p ( 1 , 2 ] and δ [ 0 , ) . Let Ω R 3 be a bounded domain with C 2 , 1 boundary, and let f L p ( Ω ) . Then there exists a unique weak solution u to (1.1) such that

Ω ϕ ( | 𝐃𝐮 | ) 𝑑 𝐱 c ( 𝐟 p , δ ) .

Proof.

The assertions follow directly from the assumptions by using the theory of monotone operators. ∎

In order to justify some of the following computations, we find it convenient to consider a perturbed problem, where we add to the tensor field 𝐒 with ( p , δ ) -structure a linear perturbation. Using again the theory of monotone operators one can easily prove the following proposition.

Proposition 2.7.

Let the tensor field S in (1.1) have ( p , δ ) -structure for some p ( 1 , 2 ] and δ [ 0 , ) , and let f L p ( Ω ) be given. Then there exists a unique weak solution u ε W 0 1 , 2 ( Ω ) of the problem

(2.12) { - div 𝐒 ε ( 𝐃𝐮 ε ) = 𝐟 in  Ω , 𝐮 ε = 𝟎 on  Ω ,

where

𝐒 ε ( 𝐐 ) := ε 𝐐 + 𝐒 ( 𝐐 ) , with  ε > 0 ,

i.e., u ε satisfies, for all v W 0 1 , 2 ( Ω ) ,

Ω 𝐒 ε ( 𝐃𝐮 ε ) 𝐃𝐯 𝑑 𝐱 = Ω 𝐟 𝐯 𝑑 𝐱 .

The solution u ε satisfies the estimate

(2.13) Ω ε | 𝐮 ε | 2 + ϕ ( | 𝐃𝐮 ε | ) d 𝐱 c ( 𝐟 p , δ ) .

Remark 2.8.

In fact, one could already prove more at this point. Namely, that for ε 0 , the unique solution 𝐮 ε converges to the unique weak solution 𝐮 of the unperturbed problem (1.1). Let us sketch the argument only, since later we get the same result with different easier arguments. From (2.13) and the properties of 𝐒 follows that

𝐮 ε 𝐮 in  W 0 1 , p ( Ω ) ,
𝐒 ( 𝐃𝐮 ε ) 𝝌 in  L p ( Ω ) .

Passing to the limit in the weak formulation of the perturbed problem, we get

Ω 𝝌 𝐃𝐯 𝑑 𝐱 = Ω 𝐟 𝐯 𝑑 𝐱 for all  𝐯 W 0 1 , p ( Ω ) .

One can not show directly that lim ε 0 Ω ε 𝐃𝐮 ε ( 𝐃𝐮 ε - 𝐃𝐮 ) 𝑑 𝐱 = 0 , since 𝐃𝐮 belongs to L p ( Ω ) only. Instead one uses the Lipschitz truncation method (cf. [10, 22]). Denoting by 𝐯 ε , j the Lipschitz truncation of ξ ( 𝐮 ε - 𝐮 ) , where ξ C 0 ( Ω ) is a localization, one can show, using the ideas from [10, 22], that

lim sup ε 0 | Ω ( 𝐒 ( 𝐃𝐮 ε ) - 𝐒 ( 𝐃𝐮 ) ) 𝐃𝐯 ε , j d 𝐱 | = 0 ,

which implies 𝐃𝐮 ε 𝐃𝐮 almost everywhere in Ω. Consequently, we have 𝝌 = 𝐒 ( 𝐃𝐮 ) , since weak and a.e. limits coincide.

2.4 Description and properties of the boundary

We assume that the boundary Ω is of class C 2 , 1 , that is, for each point P Ω , there are local coordinates such that, in these coordinates, we have P = 0 and Ω is locally described by a C 2 , 1 -function, i.e., there exist R P , R P ( 0 , ) , r P ( 0 , 1 ) and a C 2 , 1 -function a P : B R P 2 ( 0 ) B R P 1 ( 0 ) such that

  1. 𝐱 Ω ( B R P 2 ( 0 ) × B R P 1 ( 0 ) ) x 3 = a P ( x 1 , x 2 ) ,

  2. Ω P := { ( x , x 3 ) x = ( x 1 , x 2 ) B R P 2 ( 0 ) , a P ( x ) < x 3 < a P ( x ) + R P } Ω ,

  3. a P ( 0 ) = 𝟎 , and for all x = ( x 1 , x 2 ) B R P 2 ( 0 ) , | a P ( x ) | < r P ,

where B r k ( 0 ) denotes the k-dimensional open ball with center 0 and radius r > 0 . Note that r P can be made arbitrarily small if we make R P small enough. In the sequel, we will also use, for 0 < λ < 1 , the following scaled open sets:

λ Ω P := { ( x , x 3 ) x = ( x 1 , x 2 ) B λ R P 2 ( 0 ) , a P ( x ) < x 3 < a P ( x ) + λ R P } Ω P .

To localize near to Ω Ω P for P Ω , we fix smooth functions ξ P : 3 such that

  1. χ 1 2 Ω P ( 𝐱 ) ξ P ( 𝐱 ) χ 3 4 Ω P ( 𝐱 ) ,

where χ A ( 𝐱 ) is the indicator function of the measurable set A. For the remaining interior estimate, we localize by a smooth function 0 ξ 00 1 , with spt ξ 00 Ω 00 , where Ω 00 Ω is an open set such that dist ( Ω 00 , Ω ) > 0 . Since the boundary Ω is compact, we can use an appropriate finite sub-covering which, together with the interior estimate, yields the global estimate.

Let us introduce the tangential derivatives near the boundary. To simplify the notation, we fix P Ω , h ( 0 , R P 16 ) , and simply write ξ := ξ P , a := a P . We use the standard notation 𝐱 = ( x , x 3 ) and denote by 𝐞 i , i = 1 , 2 , 3 , the canonical orthonormal basis in 3 . In the following, lower-case Greek letters take the values 1 , 2 . For a function g, with spt g spt ξ , we define, for α = 1 , 2 ,

g τ ( x , x 3 ) = g τ α ( x , x 3 ) := g ( x + h 𝐞 α , x 3 + a ( x + h 𝐞 α ) - a ( x ) ) ,

and if Δ + g := g τ - g , we define tangential divided differences by d + g := h - 1 Δ + g . It holds that, if g W 1 , 1 ( Ω ) , then we have for α = 1 , 2

(2.14) d + g τ g = τ α g := α g + α a 3 g as  h 0 ,

almost everywhere in spt ξ (cf. [18, Section 3]). Conversely, uniform L q -bounds for d + g imply that τ g belongs to L q ( spt ξ ) . For simplicity, we denote a := ( 1 a , 2 a , 0 ) . The following variant of integration by parts will be often used.

Lemma 2.9.

Let spt g spt f spt ξ , and let h be small enough. Then

Ω f g - τ 𝑑 𝐱 = Ω f τ g 𝑑 𝐱 .

Consequently, Ω f d + g 𝑑 x = Ω ( d - f ) g 𝑑 x . Moreover, if in addition f and g are smooth enough and at least one vanishes on Ω , then

Ω f τ g d 𝐱 = - Ω ( τ f ) g 𝑑 𝐱 .

3 Proof of the main result

In the proof of the main result we use finite differences to show estimates in the interior and in tangential directions near the boundary, and for calculations involving directly derivatives in “normal” directions near the boundary. In order to justify that all occurring quantities are well posed, we perform the estimate for the approximate system (2.12).

The first intermediate step is the following result for the approximate problem.

Proposition 3.1.

Let the tensor field S in (1.1) have ( p , δ ) -structure for some p ( 1 , 2 ] and δ ( 0 , ) , and let F be the associated tensor field to S . Let Ω R 3 be a bounded domain with C 2 , 1 boundary and let f L p ( Ω ) . Then the unique weak solution u ε W 0 1 , 2 ( Ω ) of the approximate problem (2.12) satisfies[3]

(3.1) { Ω ε ξ 0 2 | 2 𝐮 ε | 2 + ξ 0 2 | 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 c ( 𝐟 p , ξ 0 2 , , δ ) , Ω ε ξ P 2 | τ 𝐃𝐮 ε | 2 + ξ P 2 | τ 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 c ( 𝐟 p , ξ P 2 , , a P C 2 , 1 , δ ) .

Here ξ 0 is a cut-off function with support in the interior of Ω, while for arbitrary P Ω the function ξ P is a cut-off function with support near to the boundary Ω , as defined in Section 2.4. The tangential derivative τ is defined locally in Ω P by (2.14). Moreover, there exists a constant C 1 > 0 such that[4]

(3.2) Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 + ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 c ( 𝐟 p , ξ 2 , , a C 2 , 1 , δ - 1 , ε - 1 , C 1 ) ,

provided that in the local description of the boundary, we have r P < C 1 in (b3).

In particular, these estimates imply that u ε W 2 , 2 ( Ω ) and that (2.12) holds pointwise a.e. in Ω.

The two estimates in (3.1) are uniform with respect to ε and could be also proved directly for problem (1.1). However, the third estimate (3.2) depends on ε and it is needed to justify all subsequent steps, which will give the proof of an estimate uniformly in ε , by using a different technique.

Proof of Proposition 3.1.

The proof of estimate (3.1) is very similar, being in fact a simplification (due to the fact that there is no pressure term involved) to the proof of the results in [4, Theorems 2.28 and 2.29]. On the other hand, the proof of (3.2) is different from the one in [4] due to the missing divergence constraint. In fact, it adapts techniques known from nonlinear elliptic systems. For the convenience of the reader, we recall the main steps here.

Fix P Ω and use in Ω P

𝐯 = d - ( ξ 2 d + ( 𝐮 ε | 1 2 Ω P ) ) ,

where ξ := ξ P , a := a P , and h ( 0 , R P 16 ) , as a test function in the weak formulation of (2.12). This yields

Ω ξ 2 d + 𝐒 ε ( 𝐃𝐮 ε ) d + 𝐃𝐮 ε 𝑑 𝐱 = - Ω 𝐒 ε ( 𝐃𝐮 ε ) ( ξ 2 d + 3 𝐮 ε - ( ξ - τ d - ξ + ξ d - ξ ) 3 𝐮 ε ) 𝑠 d - a d 𝐱
- Ω 𝐒 ε ( 𝐃𝐮 ε ) ξ 2 ( 3 𝐮 ε ) τ 𝑠 d - d + a - 𝐒 ε ( 𝐃𝐮 ε ) d - ( 2 ξ ξ 𝑠 d + 𝐮 ε ) d 𝐱
+ Ω 𝐒 ε ( ( 𝐃𝐮 ε ) τ ) ( 2 ξ 3 ξ d + 𝐮 ε + ξ 2 d + 3 𝐮 ε ) 𝑠 d + a d 𝐱
+ Ω 𝐟 d - ( ξ 2 d + 𝐮 ε ) d 𝐱 = : j = 1 8 I j .

From the assumption on 𝐒 , Proposition 2.3 and [4, Lemma 3.11], we have the following estimate:

Ω ε ξ 2 | d + 𝐮 ε | 2 + ε ξ 2 | d + 𝐮 ε | 2 + | d + 𝐅 ( 𝐃𝐮 ε ) | 2 + ϕ ( ξ | d + 𝐮 | ) + ϕ ( ξ | d + 𝐮 | ) d 𝐱
c Ω ξ 2 d + 𝐒 ε ( 𝐃𝐮 ε ) d + 𝐃𝐮 ε 𝑑 𝐱 + c ( ξ 1 , , a C 1 , 1 ) Ω spt ξ ϕ ( | 𝐮 ε | ) 𝑑 𝐱 .

The terms I 1 I 7 are estimated exactly as in [4, Equations (3.17)–(3.22)], while I 8 is estimated as the term I 15 in [4, (4.20)]. Thus, we get

Ω ε ξ 2 | d + 𝐮 ε | 2 + ε ξ 2 | d + 𝐮 ε | 2 + ξ 2 | d + 𝐅 ( 𝐃𝐮 ε ) | 2 + ϕ ( ξ | d + 𝐮 ε | ) + ϕ ( ξ | d + 𝐮 ε | ) d 𝐱
c ( 𝐟 p , ξ 2 , , a C 2 , 1 , δ ) .

This proves the second estimate in (3.1) by standard arguments. The first estimate in (3.1) is proved in the same way with many simplifications, since we work in the interior where the method works for all directions. This estimate implies that 𝐮 ε W loc 2 , 2 ( Ω ) and that the system (2.12) is well-defined point-wise a.e. in Ω.

To estimate the derivatives in the x 3 direction, we use equation (2.12) and it is at this point that we have changes with respect to the results in [4]. In fact, as usual in elliptic problems, we have to recover the partial derivatives with respect to x 3 by using the information on the tangential ones. In this problem the main difficulty is that the leading order term is nonlinear and depends on the symmetric part of the gradient. Thus, we have to exploit the properties of ( p , δ ) -structure of the tensor 𝐒 (cf. Definition 2.1). Denoting, for i = 1 , 2 , 3 ,[5] 𝔣 i := - f i - γ σ S i 3 ( 𝐃𝐮 ε ) 3 D γ σ 𝐮 ε - k , l = 1 3 k l S i β ( 𝐃𝐮 ε ) β D k l 𝐮 ε , we can re-write the equations in (2.12) as follows:

k = 1 3 k 3 S i 3 ( 𝐃𝐮 ε ) 3 D k 3 𝐮 ε + 3 α S i 3 ( 𝐃𝐮 ε ) 3 D 3 α 𝐮 ε = 𝔣 i a.e. in  Ω .

Contrary to the corresponding equality [4, Equation (3.49)], here we use directly all the equations in (1.1), and not only the first two. Now we multiply these equations not by 3 D i 3 𝐮 ε as expected, but by 3 D ^ i 3 𝐮 ε , where D ^ α β 𝐮 ε = 0 for α , β = 1 , 2 , D ^ α 3 𝐮 ε = D ^ 3 α 𝐮 ε = 2 D α 3 𝐮 ε for α = 1 , 2 , D ^ 33 𝐮 ε = D 33 𝐮 ε . Summing over i = 1 , 2 , 3 , we get, by using the symmetries in Remark 2.2 (iii), that

4 α 3 S β 3 ε ( 𝐃𝐮 ε ) 3 D α 3 𝐮 ε 3 D β 3 𝐮 ε + 2 α 3 S 33 ε ( 𝐃𝐮 ε ) 3 D α 3 𝐮 ε 3 D 33 𝐮 ε
(3.3) + 2 33 S β 3 ε ( 𝐃𝐮 ε ) 3 D 33 𝐮 ε 3 D β 3 𝐮 ε + 33 S 33 ε ( 𝐃𝐮 ε ) 3 D 33 𝐮 ε 3 D 33 𝐮 ε = i = 1 3 𝔣 i 3 D ^ i 3 𝐮 ε a.e. in  Ω .

To obtain a lower bound for the left-hand side, we observe that the terms on the left-hand side of (3.3) containing 𝐒 are equal to

i , j , k , l = 1 3 k l S i j ( 𝐃𝐮 ε ) Q i j Q k l

if we choose 𝐐 = 3 𝐃 ¯ 𝐮 ε , where D ¯ α β 𝐮 ε = 0 for α , β = 1 , 2 , D ¯ α 3 𝐮 ε = D ¯ 3 α 𝐮 ε = D α 3 𝐮 ε for α = 1 , 2 , and D ¯ 33 𝐮 ε = D 33 𝐮 ε . Thus, it follows from the coercivity estimate in (2.4) that these terms are bounded from below by κ 0 ϕ ′′ ( | 𝐃𝐮 ε | ) | 3 𝐃 ¯ 𝐮 ε | 2 . Similarly, we see that the remaining terms on the left-hand side of (3.3) are equal to ε | 3 𝐃 ¯ 𝐮 ε | 2 . Denoting 𝔟 i := 3 D i 3 𝐮 ε , i = 1 , 2 , 3 , we see that | 𝔟 | | 𝐃 ^ 𝐮 ε | | 𝐃 ¯ 𝐮 ε | . Consequently, we get from (3.3) the estimate

( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | 𝖇 | | 𝖋 | a.e. in  Ω .

By straightforward manipulations (cf. [4, Sections 3.2 and 4.2]), we can estimate the right-hand side as follows:

| 𝖋 | c ( | 𝐟 | + ( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) ( | τ 𝐮 ε | + a | 2 𝐮 ε | ) ) .

Note that we can deduce from 𝖇 information about 𝔟 ~ i := 33 2 u ε i , i = 1 , 2 , 3 , because

| 𝖇 | 2 | 𝖇 ~ | - | τ 𝐮 ε | - a | 2 𝐮 ε |

holds a.e. in Ω P . This and the last two inequalities imply

( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | 𝖇 ~ | c ( | 𝐟 | + ( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) ( | τ 𝐮 ε | + a | 2 𝐮 ε | ) ) a.e. in  Ω P .

Adding on both sides, for α = 1 , 2 and i , k = 1 , 2 , 3 , the term

( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | α i u ε k | ,

and using on the right-hand side the definition of the tangential derivative (cf. (2.14)), we finally arrive at

( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | 2 𝐮 ε | c ( | 𝐟 | + ( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) ( | τ 𝐮 ε | + a | 2 𝐮 ε | ) ) ,

which is valid a.e. in Ω P . Note that the constant c only depends on the characteristics of 𝐒 . Next, we can choose the open sets Ω P in such a way that a P ( x ) , Ω P is small enough, so that we can absorb the last term from the right-hand side, which yields

( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | 2 𝐮 ε | c ( | 𝐟 | + ( ε + ϕ ′′ ( | 𝐃𝐮 ε | ) ) | τ 𝐮 ε | ) a.e. in  Ω P ,

where again the constant c only depends on the characteristics of 𝐒 . By neglecting the second term on the left-hand side (which is non-negative), raising the remaining inequality to the power 2, and using that 𝐒 has ( p , δ ) -structure for p < 2 , we obtain

Ω ε ξ P 2 | 2 𝐮 ε | 2 d 𝐱 c Ω | 𝐟 | 2 d 𝐱 + ( ε + δ 2 ( p - 2 ) ) ε ( ε Ω ξ P 2 | τ 𝐮 ε | 2 d 𝐱 ) .

The already proven results on tangential derivatives and Korn’s inequality imply that the last integral from the right-hand side is finite. Thus, the properties of the covering imply the estimate in (3.2). ∎

3.1 Improved estimates for normal derivatives

In the proof of (3.2), we used system (2.12) and obtained an estimate that is not uniform with respect to ε . In this section, by following the ideas in [24], we proceed differently and estimate 𝒫 3 in terms of the quantities occurring in (3.1). The main technical step of the paper is the proof of the following result.

Proposition 3.2.

Let the hypotheses in Theorem 1.1 be satisfied with δ > 0 , and let the local description a P of the boundary and the localization function ξ P satisfy (b1)–(b3) and ( 1) (cf. Section 2.4). Then there exists a constant C 2 > 0 such that the weak solution u ε W 0 1 , 2 ( Ω ) of the approximate problem (2.12) satisfies,[6] for every P Ω ,

Ω ε ξ P 2 | 3 𝐃𝐮 ε | 2 + ξ P 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 C ( 𝐟 p , ξ P 2 , , a P C 2 , 1 , δ , C 2 ) ,

provided r P < C 2 in (b3).

Proof.

Let us fix an arbitrary point P Ω and a local description a = a P of the boundary and the localization function ξ = ξ P satisfying (b1)–(b3) and ( 1). In the following we denote by C constants that depend only on the characteristics of 𝐒 . First we observe that, by the results of Proposition 2.4, there exists a constant C 0 , depending only on the characteristics of 𝐒 , such that

1 C 0 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝒫 3 ( 𝐃𝐮 ε ) a.e. in  Ω .

Thus, using also the symmetry of 𝐃𝐮 ε and 𝐒 , we get

Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 + 1 C 0 ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 Ω ξ 2 ( ε 3 𝐃𝐮 ε + 3 𝐒 ( 𝐃𝐮 ε ) ) 3 𝐃𝐮 ε d 𝐱
= Ω i , j = 1 3 ξ 2 ( ε 3 D i j 𝐮 ε + 3 S i j ( 𝐃𝐮 ε ) ) 3 j u i d 𝐱
= Ω ξ 2 ( ε 3 D α β 𝐮 ε + 3 S α β ( 𝐃𝐮 ε ) ) 3 D α β 𝐮 ε d 𝐱
+ Ω ξ 2 ( ε 3 D 3 α 𝐮 ε + 3 S 3 α ( 𝐃𝐮 ε ) ) α D 33 𝐮 ε d 𝐱
+ Ω j = 1 3 ξ 2 3 ( ε D j 3 𝐮 ε + S j 3 ( 𝐃𝐮 ε ) ) 3 2 u ε j d 𝐱
= : I 1 + I 2 + I 3 .

To estimate I 2 , we multiply and divide by the quantity ϕ ′′ ( | 𝐃𝐮 ε | ) 0 , and use Young’s inequality and Proposition 2.4. This yields that, for all λ > 0 , there exists c λ > 0 such that

| I 2 | α = 1 2 Ω ξ 2 | 3 𝐒 ( 𝐃𝐮 ε ) | | α 𝐃𝐮 ε | ϕ ′′ ( | 𝐃𝐮 ε | ) ϕ ′′ ( | 𝐃𝐮 ε | ) d 𝐱 + λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 α = 1 2 Ω ε ξ 2 | α 𝐃𝐮 ε | 2 𝑑 𝐱
λ Ω ξ 2 | 3 𝐒 ( 𝐃𝐮 ε ) | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱 + c λ - 1 α = 1 2 Ω ξ 2 ϕ ′′ ( | 𝐃𝐮 ε | ) | α 𝐃𝐮 ε | 2 𝑑 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 α = 1 2 Ω ε ξ 2 | α 𝐃𝐮 ε | 2 𝑑 𝐱
C λ Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 a = 1 2 Ω ξ 2 | α 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 α = 1 2 Ω ε ξ 2 | α 𝐃𝐮 ε | 2 𝑑 𝐱 .

Here and in the following we denote by c λ - 1 constants that may depend on the characteristics of 𝐒 and on λ - 1 , while C denotes constants that may depend on the characteristics of 𝐒 only.

To treat the third integral I 3 , we proceed as follows: We use the following well-known algebraic identity, valid for smooth enough vectors 𝐯 and i , j , k = 1 , 2 , 3 :

(3.4) j k v i = j D i k 𝐯 + k D i j 𝐯 - i D j k 𝐯 ,

and equations (2.12) point-wise, which can be written, for j = 1 , 2 , 3 , as

3 ( ε D j 3 𝐮 ε + S j 3 ( 𝐃𝐮 ε ) ) = - f j - β ( ε D j β 𝐮 ε + S j β ( 𝐃𝐮 ε ) ) a.e. in  Ω .

This is possible due to Proposition 3.1. Hence, we obtain

| I 3 | j = 1 3 | Ω ξ 2 ( - f j - β S j β ( 𝐃𝐮 ε ) - ε β D j β 𝐮 ε ) ( 2 3 D j 3 𝐮 ε - j D 33 𝐮 ε ) d 𝐱 | .

The right-hand side can be estimated in a way similar to I 2 . This yields that, for all λ > 0 , there exists c λ > 0 such that

| I 3 | Ω ξ 2 ( | 𝐟 | + β = 1 2 | β 𝐒 ( 𝐃𝐮 ε ) | ) ( 2 | 3 𝐃𝐮 ε | + α = 1 2 | α 𝐃𝐮 ε | ) ϕ ′′ ( | 𝐃𝐮 ε | ) ϕ ′′ ( | 𝐃𝐮 ε | ) d 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱
λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱
+ c λ - 1 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 Ω ξ 2 | 𝐟 | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱
λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱
+ c λ - 1 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 ( 𝐟 p p + 𝐃𝐮 ε p p + δ p ) .

Observe that we used p 2 to estimate the term involving 𝐟 .

To estimate I 1 , we employ the algebraic identity (3.4) to split the integral as follows:

I 1 = Ω ξ 2 ( ε 3 D α β 𝐮 ε + 3 S α β ( 𝐃𝐮 ε ) ) ( α D 3 β 𝐮 ε + β D 3 α 𝐮 ε ) 𝑑 𝐱 - Ω ξ 2 ( ε 3 D α β 𝐮 ε + 3 S α β ( 𝐃𝐮 ε ) ) β α u ε 3 d 𝐱
= : A + B .

The first term is estimated in a way similar to I 2 , yielding, for all λ > 0 ,

| A | C λ Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱 .

To estimate B we observe that by the definition of the tangential derivative, we have

α β u ε 3 = α τ β u ε 3 - ( α β a ) D 33 𝐮 ε - ( β a ) α D 33 𝐮 ε ,

and consequently the term B can be split into the following three terms:

- Ω ξ 2 ( ε 3 D α β 𝐮 ε + 3 S α β ( 𝐃𝐮 ε ) ) ( α τ β u ε 3 - ( α β a ) D 33 𝐮 ε - ( β a ) α D 33 𝐮 ε ) d 𝐱 = : B 1 + B 2 + B 3 .

We estimate B 2 as follows:

| B 2 | Ω ξ 2 | 3 𝐒 ( 𝐃𝐮 ε ) | | 2 a | | 𝐃𝐮 ε | ϕ ′′ ( | 𝐃𝐮 ε | ) ϕ ′′ ( | 𝐃𝐮 ε | ) + ε ξ 2 | 3 𝐃𝐮 ε | | 2 a | | 𝐃𝐮 ε | d 𝐱
λ Ω ξ 2 | 3 𝐒 ( 𝐃𝐮 ε ) | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱 + c λ - 1 2 a 2 Ω ξ 2 | 𝐃𝐮 ε | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 2 a 2 Ω ε ξ 2 | 𝐃𝐮 ε | 2 𝑑 𝐱
λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 2 a 2 ρ ϕ ( | 𝐃𝐮 ε | ) + 1 8 Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + 2 ε 2 a 2 𝐃𝐮 ε 2 2 .

The term B 3 is estimated in a way similar to I 2 , yielding, for all λ > 0 ,

| B 3 | λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 a 2 β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱
+ λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 | a 2 β = 1 2 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱 .

Concerning the term B 1 , we would like to perform some integration by parts, which is one of the crucial observations we are adapting from [24]. Neglecting the localization ξ in B 1 , we would like to use that

Ω 3 S α β ε ( 𝐃𝐮 ε ) α τ β u ε 3 d 𝐱 = Ω α S α β ε ( 𝐃𝐮 ε ) 3 τ β u ε 3 d 𝐱 .

This formula can be justified by using an appropriate approximation that exists for 𝐮 ε W 0 1 , 2 ( Ω ) W 2 , 2 ( Ω ) , since τ 𝐮 ε = 𝟎 on Ω . More precisely, to treat the term B 1 , we use that the solution 𝐮 ε of (2.12) belongs to W 0 1 , 2 ( Ω ) W 2 , 2 ( Ω ) . Thus, τ ( 𝐮 ε | Ω P ) = 𝟎 on Ω P Ω , hence ξ P τ ( 𝐮 ε 3 ) = 𝟎 on Ω . This implies that we can find a sequence { ( 𝓢 n , 𝓤 n ) } C ( Ω ) × C 0 ( Ω ) such that ( 𝓢 n , 𝓤 n ) ( 𝐒 ε , τ 𝐮 ε ) in W 1 , 2 ( Ω ) × W 0 1 , 2 ( Ω ) , and perform calculations with ( 𝓢 n , 𝓤 n ) , showing then that all formulas of integration by parts are valid. The passage to the limit as n + is done only in the last step. For simplicity, we drop the details of this well-known argument (sketched also in [24]), and we write directly formulas without this smooth approximation. Thus, performing several integrations by parts, we get

Ω ξ 2 3 S α β ( 𝐃𝐮 ε ) α τ β u ε 3 d 𝐱 = Ω ( α ξ 2 ) S α β ( 𝐃𝐮 ε ) 3 τ β u ε 3 d 𝐱 - Ω ( 3 ξ 2 ) S α β ( 𝐃𝐮 ε ) α τ β u ε 3 d 𝐱
+ Ω ξ 2 α S α β ( 𝐃𝐮 ε ) 3 τ β u ε 3 d 𝐱

and

ε Ω ξ 2 3 D α β 𝐮 ε α τ β u ε 3 d 𝐱 = ε Ω ( α ξ 2 ) D α β 𝐮 ε 3 τ β u ε 3 d 𝐱 - ε Ω ( 3 ξ 2 ) D α β 𝐮 ε α τ β u ε 3 d 𝐱
+ ε Ω ξ 2 α D α β 𝐮 ε 3 τ β u ε 3 d 𝐱 .

This shows that

B 1 = Ω 2 ξ α ξ S α β ( 𝐃𝐮 ε ) 3 τ β u ε 3 d 𝐱 - Ω 2 ξ 3 ξ S α β ( 𝐃𝐮 ε ) α τ β u ε 3 d 𝐱
+ Ω ξ 2 α S α β ( 𝐃𝐮 ε ) 3 τ β u ε 3 d 𝐱 + ε Ω 2 ξ α ξ D α β 𝐮 ε 3 τ β u ε 3 d 𝐱
- ε Ω 2 ξ 3 ξ D α β 𝐮 ε α τ β u ε 3 d 𝐱 + ε Ω ξ 2 α D α β 𝐮 ε 3 τ β u ε 3 d 𝐱
= : B 1 , 1 + B 1 , 2 + B 1 , 3 + B 1 , 4 + B 1 , 5 + B 1 , 6 .

To estimate B 1 , 1 , B 1 , 3 , B 1 , 4 , B 1 , 6 , we observe that

3 τ β u ε 3 = τ β 3 u ε 3 = τ β D 33 𝐮 ε .

By using Young’s inequality, the growth properties of 𝐒 in (2.6) and (2.8), we get

| B 1 , 1 | ξ 2 Ω | 𝐒 ( 𝐃𝐮 ε ) | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱 + C β = 1 2 Ω ξ 2 ϕ ′′ ( | 𝐃𝐮 ε | ) | τ β 𝐃𝐮 ε | 2 𝑑 𝐱
ξ 2 ρ ϕ ( | 𝐃𝐮 ε | ) + C β = 1 2 Ω ξ 2 | τ β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱

and

| B 1 , 3 | β = 1 2 Ω ξ 2 | β 𝐒 α β ( 𝐃𝐮 ε ) | 2 ϕ ′′ ( | 𝐃𝐮 ε | ) 𝑑 𝐱 + β = 1 2 Ω ξ 2 ϕ ′′ ( | 𝐃𝐮 ε | ) | τ β 𝐃𝐮 ε | 2 𝑑 𝐱
C β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 + ξ 2 | τ β 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱 .

Similarly, we get

| B 1 , 4 | C ε ξ 2 𝐃𝐮 ε 2 2 + C β = 1 2 Ω ε ξ 2 | τ β 𝐃𝐮 ε | 2 𝑑 𝐱

and

| B 1 , 6 | C ε β = 1 2 Ω ξ 2 | β 𝐃𝐮 ε | 2 + ξ 2 | τ β 𝐃𝐮 ε | 2 d 𝐱 .

To estimate B 1 , 2 and B 15 , we observe that, using the algebraic identity (3.4) and the definition of the tangential derivative,

α τ β u ε 3 = α ( β 𝐮 ε 3 + β a 3 u ε 3 )
= α β u ε 3 + α β a D 33 𝐮 ε + β a α D 33 𝐮 ε
= α D β 3 𝐮 ε + β D α 3 𝐮 ε - 3 D α β 𝐮 ε + α β a D 33 𝐮 ε + β a α D 33 𝐮 ε .

Hence, by substituting and again the same inequalities as before, we arrive to the following estimates:

| B 1 , 2 | λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + C ( 1 + a 2 ) β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 ( 1 + 2 a ) ξ 2 ρ ϕ ( | 𝐃𝐮 ε | ) ,
| B 1 , 5 | λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 ( 1 + a 2 ) β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱 + c λ - 1 ( 1 + 2 a ) ξ 2 ε 𝐃𝐮 ε 2 2 .

Collecting all estimates and using that a r P 1 , we finally obtain

Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 + 1 C 0 ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱
λ Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + λ C Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱
    + c λ - 1 β = 1 2 Ω ξ 2 | β 𝐅 ( 𝐃𝐮 ε ) | 2 + ξ 2 | τ β 𝐅 ( 𝐮 ε ) | 2 d 𝐱 + c λ - 1 β = 1 2 Ω ε ξ 2 | β 𝐃𝐮 ε | 2 𝑑 𝐱
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) ( 𝐟 p p + ρ ϕ ( | 𝐃𝐮 ε | ) + ρ ϕ ( δ ) )
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) 𝐃𝐮 ε 2 2 .

The quantities that are bounded uniformly in L 2 ( Ω P ) are the tangential derivatives of ε 𝐃𝐮 ε and of 𝐅 ( 𝐃𝐮 ε ) . By definition, we have

β 𝐃𝐮 ε = τ β 𝐃𝐮 ε - β a 3 𝐃𝐮 ε ,
β 𝐅 ( 𝐃𝐮 ε ) = τ β 𝐅 ( 𝐃𝐮 ε ) - β a 3 𝐅 ( 𝐃𝐮 ε ) ,

and if we substitute, we obtain

Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 + 1 C 0 ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱
( λ + 4 a 2 ) Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 𝑑 𝐱 + ( λ C + c λ - 1 a 2 ) Ω ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 𝑑 𝐱
    + c λ - 1 β = 1 2 Ω ξ 2 | τ β 𝐅 ( 𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ε ξ 2 | τ β 𝐃𝐮 ε | 2 𝑑 𝐱
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) ( 𝐟 p p + ρ ϕ ( | 𝐃𝐮 ε | ) + ρ ϕ ( δ ) )
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) 𝐃𝐮 ε 2 2 .

By choosing first λ > 0 small enough such that λ C < 4 - 1 C 0 and then choosing in the local description of the boundary R = R P small enough such that c λ a < 4 - 1 C 0 , we can absorb the first two terms from the right-hand side into the left-hand side to obtain

Ω ε ξ 2 | 3 𝐃𝐮 ε | 2 + 1 C 0 ξ 2 | 3 𝐅 ( 𝐃𝐮 ε ) | 2 d 𝐱
c λ - 1 β = 1 2 Ω ξ 2 | τ β 𝐅 ( 𝐮 ε ) | 2 𝑑 𝐱 + c λ - 1 β = 1 2 Ω ε ξ 2 | τ β 𝐃𝐮 ε | 2 𝑑 𝐱
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) ( 𝐟 p p + ρ ϕ ( | 𝐃𝐮 ε | ) + ρ ϕ ( δ ) )
    + c λ - 1 ( 1 + 2 a 2 + ( 1 + 2 a 2 ) ξ 2 ) 𝐃𝐮 ε 2 2 ,

where now c λ depends on the fixed paramater λ, the characteristics of 𝐒 and on C 2 . The right-hand side is bounded uniformly with respect to ε > 0 , due to Proposition 3.1, proving the assertion of the proposition. ∎

Choosing now an appropriate finite covering of the boundary (for the details, see also [4]), Propositions 3.13.2 yield the following result.

Theorem 3.3.

Let the hypotheses in Theorem 1.1 with δ > 0 be satisfied. Then[7]

ε 𝐃𝐮 ε 2 2 + 𝐅 ( 𝐃𝐮 ε ) 2 2 C ( 𝐟 p , δ , Ω ) .

3.2 Passage to the limit

Once this has been proved, by means of appropriate limiting process, we can show that the estimate is inherited by 𝐮 = lim ε 0 𝐮 ε , since 𝐮 is the unique solution to the boundary value problem (1.1). We can now give the proof of the main result.

Proof of Theorem 1.1.

Let us firstly assume that δ > 0 . From Propositions 2.3 and 2.7, and Theorem 3.3, we know that 𝐅 ( 𝐃𝐮 ε ) is uniformly bounded with respect to ε in W 1 , 2 ( Ω ) . This also implies (cf. [3, Lemma 4.4]) that 𝐮 ε is uniformly bounded with respect to ε in W 2 , p ( Ω ) . The properties of 𝐒 and Proposition 2.7 also yield that 𝐒 ( 𝐃𝐮 ε ) is uniformly bounded with respect to ε in L p ( Ω ) . Thus, there exists a subsequence { ε n } (which converges to 0 as n + ) , 𝐮 W 2 , p ( Ω ) , 𝐅 ~ W 1 , 2 ( Ω ) , and 𝝌 L p ( Ω ) such that

𝐮 ε n 𝐮 in  W 2 , p ( Ω ) W 0 1 , p ( Ω ) ,
𝐃𝐮 ε n 𝐃𝐮 a.e. in  Ω ,
𝐅 ( 𝐃𝐮 ε n ) 𝐅 ~ in  W 1 , 2 ( Ω ) ,
𝐒 ( 𝐃𝐮 ε n ) 𝝌 in  L p ( Ω ) .

The continuity of 𝐒 and 𝐅 , and the classical result stating that the weak limit and the a.e. limit in Lebesgue spaces coincide (cf. [12]) imply that

𝐅 ~ = 𝐅 ( 𝐃𝐮 ) and 𝝌 = 𝐒 ( 𝐃𝐮 ) .

These results enable us to pass to the limit in the weak formulation of the perturbed problem (2.12), which yields

Ω 𝐒 ( 𝐃𝐮 ) 𝐃𝐯 𝑑 𝐱 = Ω 𝐟 𝐯 𝑑 𝐱 for all  𝐯 C 0 ( Ω ) ,

where we also used that lim ε n 0 Ω ε n 𝐃𝐮 ε n 𝐃𝐯 𝑑 𝐱 = 0 . By density, we thus know that 𝐮 is the unique weak solution of problem (1.1). Finally, the lower semi-continuity of the norm implies that

Ω | 𝐅 ( 𝐃𝐮 ) | 2 d 𝐱 lim inf ε n 0 Ω | 𝐅 ( 𝐃𝐮 ε n ) | 2 d 𝐱 c .

Note that in [3, Section 4] it is shown that

𝐮 W 2 , 3 p p + 1 ( Ω ) p c ( 𝐅 ( 𝐃𝐮 ) 2 2 + δ p ) ,

which implies the Sobolev regularity stated in Theorem 1.1. This finishes the proof in the case δ > 0 .

Let us now assume that δ = 0 . Propositions 3.1 and 3.2 are valid only for δ > 0 and thus cannot be used directly for the case that 𝐒 has ( p , δ ) -structure with δ = 0 . However, it is proved in [3, Section 3.1] that for any stress tensor with ( p , 0 ) -structure 𝐒 , there exist stress tensors 𝐒 κ , having ( p , κ ) -structure with κ > 0 , approximating 𝐒 in an appropriate way.[8] Thus, we approximate (2.12) by the system

{ - div 𝐒 ε , κ ( 𝐃𝐮 ε , κ ) = 𝐟 in  Ω , 𝐮 = 𝟎 on  Ω ,

where

𝐒 ε , κ ( 𝐐 ) := ε 𝐐 + 𝐒 κ ( 𝐐 ) , with  ε > 0 , κ ( 0 , 1 ) .

For fixed κ > 0 , we can use the above theory and the fact that the estimates are uniformly in κ to pass to the limit as ε 0 . Thus, we obtain that for all κ ( 0 , 1 ) , there exists a unique 𝐮 κ W 0 1 , p ( Ω ) satisfying, for all 𝐯 W 0 1 , p ( Ω ) ,

Ω 𝐒 κ ( 𝐃𝐮 κ ) 𝐃𝐯 𝑑 𝐱 = Ω 𝐟 𝐯 𝑑 𝐱

and

(3.5) Ω | 𝐅 κ ( 𝐃𝐮 κ ) | 2 + | 𝐅 κ ( 𝐃𝐮 κ ) | 2 d 𝐱 c ( 𝐟 p , Ω ) ,

where the constant is independent of κ ( 0 , 1 ) and 𝐅 κ : 3 × 3 sym 3 × 3 is defined through

𝐅 κ ( 𝐏 ) := ( κ + | 𝐏 sym | ) p - 2 2 𝐏 sym .

Now we can proceed as in [3]. Indeed, from (3.5) and the properties of ϕ p , κ (in particular (2.3)), it follows that 𝐅 κ ( 𝐃𝐮 κ ) is uniformly bounded in W 1 , 2 ( Ω ) , that 𝐮 κ is uniformly bounded in W 0 1 , p ( Ω ) and that 𝐒 κ ( 𝐃𝐮 κ ) is uniformly bounded in L p ( Ω ) . Thus, there exist 𝐀 W 1 , 2 ( Ω ) , 𝐮 W 0 1 , p ( Ω ) , 𝝌 L p ( Ω ) , and a subsequence { κ n } , with κ n 0 , such that

𝐅 ( 𝐃𝐮 κ n ) 𝐀 in  W 1 , 2 ( Ω ) ,
𝐅 κ n ( 𝐃𝐮 κ n ) 𝐀 in  L 2 ( Ω )  and a.e. in  Ω ,
𝐮 κ n 𝐮 in  W 0 1 , p ( Ω ) ,
𝐒 κ ( 𝐃𝐮 κ ) 𝝌 in  L p ( Ω ) .

Setting 𝐁 := ( 𝐅 0 ) - 1 ( 𝐀 ) , it follows from [3, Lemma 3.23] that

𝐃𝐮 κ n = ( 𝐅 κ n ) - 1 ( 𝐅 κ n ( 𝐃𝐮 κ n ) ) ( 𝐅 0 ) - 1 ( 𝐀 ) = 𝐁 a.e. in  Ω .

Since weak and a.e. limit coincide, we obtain that

𝐃𝐮 κ n 𝐃𝐮 = 𝐁 a.e. in  Ω .

From [3, Lemma 3.16] and [3, Corollary 3.22], it now follows that

𝐅 ( 𝐃𝐮 κ n ) 𝐅 0 ( 𝐃𝐮 ) in  W 1 , 2 ( Ω ) ,
𝐒 κ n ( 𝐃𝐮 κ n ) 𝐒 ( 𝐃𝐮 ) a.e. in  Ω .

Since weak and a.e. limit coincide, we obtain that

𝐃𝐮 = 𝝌 a.e. in  Ω .

Now we can finish the proof in the same way as in the case δ > 0 . ∎

Funding statement: The research that led to the present paper was partially supported by a grant of the group GNAMPA of INdAM.

References

[1] E. Acerbi and N. Fusco, Regularity for minimizers of nonquadratic functionals: The case 1 < p < 2 , J. Math. Anal. Appl. 140 (1989), no. 1, 115–135. 10.1016/0022-247X(89)90098-XSearch in Google Scholar

[2] H. Beirão da Veiga, P. Kaplický and M. Růžička, Boundary regularity of shear thickening flows, J. Math. Fluid Mech. 13 (2011), no. 3, 387–404. 10.1007/s00021-010-0025-ySearch in Google Scholar

[3] L. C. Berselli, L. Diening and M. Růžička, Existence of strong solutions for incompressible fluids with shear dependent viscosities, J. Math. Fluid Mech. 12 (2010), no. 1, 101–132. 10.1007/s00021-008-0277-ySearch in Google Scholar

[4] L. C. Berselli and M. Růžička, Global regularity properties of steady shear thinning flows, J. Math. Anal. Appl. 450 (2017), no. 2, 839–871. 10.1016/j.jmaa.2017.01.016Search in Google Scholar

[5] R. B. Bird, R. C. Armstrong and O. Hassager, Dynamic of Polymer Liquids, 2nd ed., John Wiley, New York, 1987. Search in Google Scholar

[6] F. Crispo and C. R. Grisanti, On the existence, uniqueness and C 1 , γ ( Ω ¯ ) W 2 , 2 ( Ω ) regularity for a class of shear-thinning fluids, J. Math. Fluid Mech. 10 (2008), no. 4, 455–487. 10.1007/s00021-008-0282-1Search in Google Scholar

[7] F. Crispo and P. Maremonti, A high regularity result of solutions to modified p-Stokes equations, Nonlinear Anal. 118 (2015), 97–129. 10.1016/j.na.2014.10.017Search in Google Scholar

[8] L. Diening and F. Ettwein, Fractional estimates for non-differentiable elliptic systems with general growth, Forum Math. 20 (2008), no. 3, 523–556. 10.1515/FORUM.2008.027Search in Google Scholar

[9] L. Diening and C. Kreuzer, Linear convergence of an adaptive finite element method for the p-Laplacian equation, SIAM J. Numer. Anal. 46 (2008), no. 2, 614–638. 10.1137/070681508Search in Google Scholar

[10] L. Diening, J. Málek and M. Steinhauer, On Lipschitz truncations of Sobolev functions (with variable exponent) and their selected applications, ESAIM Control Optim. Calc. Var. 14 (2008), no. 2, 211–232. 10.1051/cocv:2007049Search in Google Scholar

[11] M. Fuchs and G. Seregin, Variational Methods for Problems from Plasticity Theory and for Generalized Newtonian Fluids, Lecture Notes in Math. 1749, Springer, Berlin, 2000. 10.1007/BFb0103751Search in Google Scholar

[12] H. Gajewski, K. Gröger and K. Zacharias, Nichtlineare Operatorgleichungen und Operatordifferentialgleichungen, Akademie, Berlin, 1974. 10.1002/mana.19750672207Search in Google Scholar

[13] M. Giaquinta and G. Modica, Remarks on the regularity of the minimizers of certain degenerate functionals, Manuscripta Math. 57 (1986), no. 1, 55–99. 10.1007/BF01172492Search in Google Scholar

[14] E. Giusti, Direct Methods in the Calculus of Variations, World Scientific, River Edge, 2003. 10.1142/5002Search in Google Scholar

[15] L. M. Kachanov, Foundations of the Theory of Plasticity, revised ed., North-Holland Ser. Appl. Math. Mech. 12, North-Holland, Amsterdam, 1971, Search in Google Scholar

[16] P. Kaplický, J. Málek and J. Stará, C 1 , α -solutions to a class of nonlinear fluids in two dimensions—stationary Dirichlet problem, Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. Inst. Steklov. (POMI) 259 (1999), 89–121, 297. Search in Google Scholar

[17] P. Kaplický and J. Tichý, Boundary regularity of flows under perfect slip boundary conditions, Cent. Eur. J. Math. 11 (2013), no. 7, 1243–1263. 10.2478/s11533-013-0232-xSearch in Google Scholar

[18] J. Málek, J. Nečas and M. Růžička, On weak solutions to a class of non-Newtonian incompressible fluids in bounded three-dimensional domains: The case p 2 , Adv. Differential Equations 6 (2001), no. 3, 257–302. 10.57262/ade/1357141212Search in Google Scholar

[19] J. Málek and K. R. Rajagopal, Mathematical issues concerning the Navier–Stokes equations and some of its generalizations, Evolutionary Equations. Vol. II, Handb. Differ. Equ., Elsevier/North-Holland, Amsterdam (2005), 371–459. 10.1016/S1874-5717(06)80008-3Search in Google Scholar

[20] J. Nečas and I. Hlaváček, Mathematical Theory of Elastic and Elasto-plastic Bodies: An Introduction, Stud. Appl. Mech. 3, Elsevier Scientific, Amsterdam, 1980. Search in Google Scholar

[21] M. M. Rao and Z. D. Ren, Theory of Orlicz Spaces, Monogr. Textb. Pure Appl. Math. 146, Marcel Dekker, New York, 1991. Search in Google Scholar

[22] M. Růžička, Analysis of generalized Newtonian fluids, Topics in Mathematical Fluid Mechanics, Lecture Notes in Math. 2073, Springer, Heidelberg (2013), 199–238. 10.1007/978-3-642-36297-2_4Search in Google Scholar

[23] M. Růžička and L. Diening, Non-Newtonian fluids and function spaces, NAFSA 8—Nonlinear Analysis, Function Spaces and Applications. Vol. 8, Czech Academy of Sciences, Prague (2007), 94–143. Search in Google Scholar

[24] G. A. Seregin and T. N. Shilkin, Regularity of minimizers of some variational problems in plasticity theory, Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. Inst. Steklov. (POMI) 243 (1997), 270–298, 342–343. 10.1007/BF02673602Search in Google Scholar

Received: 2018-04-16
Revised: 2018-07-18
Accepted: 2018-07-18
Published Online: 2018-10-11

© 2020 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Frontmatter
  2. On the moving plane method for boundary blow-up solutions to semilinear elliptic equations
  3. Regularity of solutions of the parabolic normalized p-Laplace equation
  4. Cahn–Hilliard equation on the boundary with bulk condition of Allen–Cahn type
  5. Blow-up solutions for fully nonlinear equations: Existence, asymptotic estimates and uniqueness
  6. Radon measure-valued solutions of first order scalar conservation laws
  7. Ground state solutions for a semilinear elliptic problem with critical-subcritical growth
  8. Generalized solutions of variational problems and applications
  9. Existence and non-existence results for Kirchhoff-type problems with convolution nonlinearity
  10. Nonlinear Sherman-type inequalities
  11. Global regularity for systems with p-structure depending on the symmetric gradient
  12. Homogenization of a net of periodic critically scaled boundary obstacles related to reverse osmosis “nano-composite” membranes
  13. Noncoercive resonant (p,2)-equations with concave terms
  14. Evolutionary quasi-variational and variational inequalities with constraints on the derivatives
  15. Sharp estimates on the first Dirichlet eigenvalue of nonlinear elliptic operators via maximum principle
  16. Localization and multiplicity in the homogenization of nonlinear problems
  17. Remarks on a nonlinear nonlocal operator in Orlicz spaces
  18. A Picone identity for variable exponent operators and applications
  19. On the weakly degenerate Allen-Cahn equation
  20. Continuity results for parametric nonlinear singular Dirichlet problems
  21. Construction of type I blowup solutions for a higher order semilinear parabolic equation
  22. Singularly perturbed Choquard equations with nonlinearity satisfying Berestycki-Lions assumptions
  23. Comparison results for nonlinear divergence structure elliptic PDE’s
  24. Constant sign and nodal solutions for parametric (p, 2)-equations
  25. Monotonicity formulas for coupled elliptic gradient systems with applications
  26. Berestycki-Lions conditions on ground state solutions for a Nonlinear Schrödinger equation with variable potentials
  27. A class of semipositone p-Laplacian problems with a critical growth reaction term
  28. The role of superlinear damping in the construction of solutions to drift-diffusion problems with initial data in L1
  29. Reconstruction of Tesla micro-valve using topological sensitivity analysis
  30. Lewy-Stampacchia’s inequality for a pseudomonotone parabolic problem
  31. Global well-posedness of nonlinear wave equation with weak and strong damping terms and logarithmic source term
  32. Regularity Criteria for Navier-Stokes Equations with Slip Boundary Conditions on Non-flat Boundaries via Two Velocity Components
  33. Homoclinics for singular strong force Lagrangian systems
  34. A constructive method for convex solutions of a class of nonlinear Black-Scholes equations
  35. On a class of nonlocal nonlinear Schrödinger equations with potential well
  36. Superlinear Schrödinger–Kirchhoff type problems involving the fractional p–Laplacian and critical exponent
  37. Regularity for minimizers for functionals of double phase with variable exponents
  38. Boundary blow-up solutions to the Monge-Ampère equation: Sharp conditions and asymptotic behavior
  39. Homogenisation with error estimates of attractors for damped semi-linear anisotropic wave equations
  40. A-priori bounds for quasilinear problems in critical dimension
  41. Critical growth elliptic problems involving Hardy-Littlewood-Sobolev critical exponent in non-contractible domains
  42. On the Sobolev space of functions with derivative of logarithmic order
  43. On a logarithmic Hartree equation
  44. Critical elliptic systems involving multiple strongly–coupled Hardy–type terms
  45. Sharp conditions of global existence for nonlinear Schrödinger equation with a harmonic potential
  46. Existence for (p, q) critical systems in the Heisenberg group
  47. Periodic traveling fronts for partially degenerate reaction-diffusion systems with bistable and time-periodic nonlinearity
  48. Some hemivariational inequalities in the Euclidean space
  49. Existence of standing waves for quasi-linear Schrödinger equations on Tn
  50. Periodic solutions for second order differential equations with indefinite singularities
  51. On the Hölder continuity for a class of vectorial problems
  52. Bifurcations of nontrivial solutions of a cubic Helmholtz system
  53. On the exact multiplicity of stable ground states of non-Lipschitz semilinear elliptic equations for some classes of starshaped sets
  54. Sign-changing multi-bump solutions for the Chern-Simons-Schrödinger equations in ℝ2
  55. Positive solutions for diffusive Logistic equation with refuge
  56. Null controllability for a degenerate population model in divergence form via Carleman estimates
  57. Eigenvalues for a class of singular problems involving p(x)-Biharmonic operator and q(x)-Hardy potential
  58. On the convergence analysis of a time dependent elliptic equation with discontinuous coefficients
  59. Multiplicity and concentration results for magnetic relativistic Schrödinger equations
  60. Solvability of an infinite system of nonlinear integral equations of Volterra-Hammerstein type
  61. The superposition operator in the space of functions continuous and converging at infinity on the real half-axis
  62. Estimates by gap potentials of free homotopy decompositions of critical Sobolev maps
  63. Pseudo almost periodic solutions for a class of differential equation with delays depending on state
  64. Normalized multi-bump solutions for saturable Schrödinger equations
  65. Some inequalities and superposition operator in the space of regulated functions
  66. Area Integral Characterization of Hardy space H1L related to Degenerate Schrödinger Operators
  67. Bifurcation of time-periodic solutions for the incompressible flow of nematic liquid crystals in three dimension
  68. Morrey estimates for a class of elliptic equations with drift term
  69. A singularity as a break point for the multiplicity of solutions to quasilinear elliptic problems
  70. Global and non global solutions for a class of coupled parabolic systems
  71. On the analysis of a geometrically selective turbulence model
  72. Multiplicity of positive solutions for quasilinear elliptic equations involving critical nonlinearity
  73. Lack of smoothing for bounded solutions of a semilinear parabolic equation
  74. Gradient estimates for the fundamental solution of Lévy type operator
  75. π/4-tangentiality of solutions for one-dimensional Minkowski-curvature problems
  76. On the existence and multiplicity of solutions to fractional Lane-Emden elliptic systems involving measures
  77. Anisotropic problems with unbalanced growth
  78. On a fractional thin film equation
  79. Minimum action solutions of nonhomogeneous Schrödinger equations
  80. Global existence and blow-up of weak solutions for a class of fractional p-Laplacian evolution equations
  81. Optimal rearrangement problem and normalized obstacle problem in the fractional setting
  82. A few problems connected with invariant measures of Markov maps - verification of some claims and opinions that circulate in the literature
Downloaded on 5.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/anona-2018-0090/html
Scroll to top button