Skip to main content

L p -estimates for a transmission problem of mixed elliptic-parabolic type

Abstract

We consider the situation when an elliptic problem in a subdomain Ω 1 of an n-dimensional bounded domain Ω is coupled via inhomogeneous canonical transmission conditions to a parabolic problem in Ω Ω 1 . In particular, we can treat elliptic-parabolic equations in bounded domains with discontinuous coefficients. Using Fourier multiplier techniques, we prove an a priori estimate for strong solutions to the equations in L p -Sobolev spaces.

MSC:35B45, 35M12.

1 Introduction

In the present paper we prove a priori estimates in L p -Sobolev spaces for the solution of a transmission problem of elliptic-parabolic type with discontinuous coefficients. More precisely, we consider a bounded domain Ω R n which is divided into two subdomains Ω 1 , Ω 2 separated by a closed contour ΓΩ and a boundary value problem of the form

A ( x , D ) u = f 1 in  Ω 1 , ( A ( x , D ) λ ) u = f 2 in  Ω 2 , C ( x , D ) u = h on  Ω .
(1.1)

Here A(x,D) is a differential operator of order 2m, C(x,D) is a column of boundary operators C 1 ,, C m , and λ is a complex parameter. We assume f k L p ( Ω k ) and are looking for a solution u W p 2 m (Ω). The top-order coefficients of the operator A(x,D) are assumed to be continuous up to the boundary in each subdomain Ω k but may have jumps across the interface Γ. The condition u W p 2 m (Ω) leads to the canonical transmission conditions along Γ, given by

[ [ ν j 1 u ] ] =0(j=1,,2m),
(1.2)

where [[ ν j 1 u]] stands for the jump of the (j1)th normal derivative of u along the interface Γ. Generalizing (1.2), we will consider inhomogeneous transmission conditions of the form

B(x,D)u=g,
(1.3)

where B= ( B 1 , , B 2 m ) with

B j (x,D)u:= ν j 1 u 1 ν j 1 u 2 (j=1,,2m).

Here we have set u k :=u | Ω k for k=1,2.

The aim of the paper is to prove uniform a priori estimates for the solutions of (1.1) and (1.3) under suitable ellipticity and smoothness assumptions on A and C, see Section 2 below for the precise formulations. To give an idea of our results, let us for the moment assume that f 1 = f 2 =0 and h=0 in (1.1) and (1.3). In classical elliptic theory, in the case of an uncoupled system we would expect a uniform a priori estimate of the form

u 1 W p 2 m ( Ω 1 ) C ( j = 1 m g j W p 2 m j + 1 1 / p ( Γ ) + u 1 L p ( Ω 1 ) ) .

On the other hand, the classical parabolic (in the sense of parameter-elliptic) a priori estimate would read

u 2 W p 2 m ( Ω 2 ) C j = 1 m g j W p 2 m j + 1 1 / p ( Γ ) .

Here W p s := W p s + | λ | s / 2 m L p is the typical parameter-dependent norm appearing in parabolic theory. Concerning the coupled system (1.1), (1.3), the question arises if we still have similar estimates for u 1 and u 2 . We will see below that this is true in some sense. More precisely, we will obtain

u 1 W p 2 m ( Ω 1 ) + u 2 W p m ( Ω 2 ) C ( j = 1 2 m g j W p 2 m j + 1 1 / p ( Γ ) + u 1 L p ( Ω 1 ) ) , | λ | 1 / 2 u 1 W p m ( Ω 1 ) + u 2 W p m ( Ω 2 ) C ( j = 1 2 m g j W p 2 m j + 1 1 / p ( Γ ) + | λ | 1 / 2 u 1 L p ( Ω 1 ) ) .

This can be seen as a mixture of elliptic and parabolic a priori estimates. Note that we do not reach the full order 2m with respect to u 2 in the first inequality and not the full power |λ| with respect to u 1 in the second inequality. The general result for f0 and h0 and the precise formulation are stated in Section 2 below.

Applications of problem (1.1), (1.3) (in its parabolic form, i.e., the parameter λ being replaced by the time derivative) can be found, e.g., in [1], including the heat equation in a domain with vanishing thermal capacity in some subdomain and a model of an electric field generated by a current in a partially non-conducting domain. On the other hand, the problem under consideration is closely related to spectral problems with indefinite weight functions of the form

( A ( x , D ) λ ω ( x ) ) u=fin Ω,C(x,D)u=0on Ω.

Here ω is a weight function which may change sign and may vanish on a set of positive measure. Such spectral problems have been investigated, e.g., in a series of papers by Faierman (see [24]) and by Pyatkov [5, 6], see also [7] and the references therein. In particular, in the paper [4] a Calderón method of reduction to the boundary was applied to deal with the case where ω vanishes on a set Ω 1 of positive measure. For this, unique solvability of the Dirichlet boundary value problem in Ω 1 had to be assumed. Transmission problems of purely parabolic type (where the parameter λ is present in each subdomain) and L p -a priori estimates for their solution were considered in [8] and [23]. Transmission problems in L p were also studied, with the same methods as in the present note, by Shibata and Shimizu in [9].

A standard approach to treat transmission problems is to use (locally) a reflection technique in one subdomain resulting in a system of differential operators which are coupled by the transmission conditions. A general theory of parameter-dependent systems can be found in a series of papers by Volevich and his co-authors (see [10] and the references therein). Here the so-called Newton polygon method leads to uniform a priori estimates for the solution. However, in the present case the Newton polygon is of trapezoidal form and thus not regular. Therefore, the Newton polygon approach cannot be applied to the transmission problem (1.1). On the other side, the resulting system is not parameter-elliptic in the classical sense ([11]) and is not covered by the standard parameter-elliptic theory. We also note the connection to singularly perturbed problems where a similar Newton polygon structure appears, cf. [12]. The analysis of the elliptic-parabolic system below also serves as a starting point for more general (and nonlinear) elliptic-parabolic systems as, for instance, appearing in lithium battery models (see [13]). A detailed investigation of the nonlinear elliptic-parabolic lithium battery model and solvability in L p -Sobolev spaces can be found in the second author’s thesis [14]. In [15] and [16] mathematical models for lithium battery systems can be found which lead to inhomogeneous transmission conditions.

In Section 2 we will state the precise assumptions and the main result of the present paper. The boundary value problem is analyzed by a localization method and the investigation of the model problem in the half-space. An explicit description of the solution of the model problem (in terms of Fourier multipliers) and resulting estimates can be found in Section 3. Finally, the proof of the main a priori estimate is given in Section 4.

2 Statement of the problem and main result

Let 1<p<, nN, k N 0 :={0,1,2,}, and Ω R n be open. By ( L p (Ω), 0 , p , Ω ) and ( W p k (Ω), k , p , Ω ) we denote the Lebesgue and Sobolev spaces on Ω with their standard norms. We will further make use of the seminorms

| u | k , p , Ω := | α | = k D α u 0 , p , Ω ( k N 0 , u W p k ( Ω ) ) ,

where we used the standard notation D α := ( i ) | α | ( x 1 ) α 1 ( x n ) α n . For real non-integer s>0 let W p s (Ω):= B p p s (Ω) denote the Besov space on Ω with its standard norm. Besides the standard norms, for the treatment of parameter-elliptic problems the following parameter-dependent norms will be convenient: Let θ(0,π] and let λ Σ ¯ θ be a complex parameter, varying in a closed sector Σ ¯ θ with vertex at 0 where Σ θ :={zC{0}:|arg(z)|<θ}. Then for mN and k{0,1,,2m}, we define

u k , p , Ω := u k , p , Ω + | λ | k 2 m u 0 , p , Ω ( u W p k ( Ω ) ) .
(2.1)

On the boundary, we will consider parameter-dependent trace norms given by

u k 1 / p , p , Γ := u k 1 / p , p , Γ + | λ | k 1 / p 2 m u 0 , p , Γ ( u W p k 1 / p ( Γ ) ) .

By Fu we denote the Fourier transform of u, and ( F u)( ξ , x n ) stands for the partial Fourier transform with respect to the first n1 variables x :=( x 1 ,, x n 1 ).

Let Ω R n be a bounded domain with boundary Ω of class C 2 m 1 , 1 , and let Γ be a closed C 2 m 1 , 1 Jordan contour in Ω, having no points with Ω in common. Denote by Ω 1 and Ω 2 the resulting subdomains such that Ω 1 Ω 2 =, Ω 1 =Γ, and Ω ¯ = Ω 1 Ω 2 ¯ . Note that, due to our assumptions, there is no contact point of Γ and Ω. We define u i :=u | Ω i and will consider the differential operators A 1 (x,D)=A(x,D) for x Ω 1 and A ˜ 2 (x,D)λ=A(x,D)λ for x Ω 2 . Slightly generalizing the form of equation (1.1), we consider differential operators of even order 2m of the following structure:

A 1 (x,D)= | α | 2 m a α ( 1 ) (x) D α and A ˜ 2 (x,D,λ)= | α | + k 2 m a α k ( 2 ) (x) λ k / 2 m D α

with mN and λ Σ ¯ θ for some θ[0,π). Furthermore, let the boundary operators C j of order 0 m j 2m1 be of the form

C j (x,D)= | γ | m j c j γ (x) D γ ,

being defined on Ω. We will write for short (A, C 1 ,, C m ) when we refer to the boundary value problem (1.1).

Assumption 2.1

  1. (1)

    Smoothness assumptions on the coefficients. We assume

    a α ( 1 ) { C ( Ω 1 ¯ ) ( | α | = 2 m ) , L ( Ω 1 ) ( | α | < 2 m ) , a α k ( 2 ) { C ( Ω 2 ¯ ) ( | α | = 2 m ) , L ( Ω 2 ) ( | α | < 2 m )

    for the coefficients of the differential operators and c j γ C 2 m m j 1 , 1 (Ω) for the coefficients of the boundary operators.

  2. (2a)

    Ellipticity of A 1 . For the principal symbol A 1 0 (x,ξ):= | α | = 2 m a α ( 1 ) (x) ξ α , we have A 1 0 (x,ξ)0 (x Ω ¯ 1 , ξ R n {0}).

  3. (2b)

    Ellipticity with parameter of the boundary value problem ( A ˜ 2 , C 1 ,, C m ). The principal symbol of A ˜ 2 satisfies

    A ˜ 2 0 (x,ξ,λ):= | α | + k = 2 m a α k ( 2 ) (x) λ k / 2 m ξ α 0

    for all x Ω ¯ 2 and all (ξ,λ)( R n × Σ ¯ θ ){(0,0)}, and the Shapiro-Lopatinskii condition is satisfied for ( A ˜ 2 , C 1 ,, C m ) at each point x 0 Ω. If C j 0 (x,D)= | γ | = m j c j γ (x) D γ denotes the principal symbol of the boundary operator, this condition reads as follows: For x 0 Ω let the boundary value problem ( A ˜ 2 , C 1 ,, C m ) be rewritten in local coordinates associated with x 0 , i.e. in coordinates resulting from the original ones by rotation and translation such that the positive x n -axis coincides with the direction of the inner normal vector. Then for all ( ξ ,λ)( R n 1 × Σ ¯ θ ){(0,0)} and h j C, the ODE problem on the halfline

    A ˜ 2 0 ( x 0 , ξ , D n , λ ) v ( x n ) = 0 in  ( 0 , ) , C j 0 ( x 0 , ξ , D n ) v ( x n ) = h j at  x n = 0 ( j = 1 , , m ) , v ( x n ) 0 ( x n )

    admits a unique solution. Here, D n :=i x n .

  4. (3)

    Assumptions on the data. We assume f 1 L p ( Ω 1 ), f 2 L p ( Ω 2 ), g j W p 2 m j + 1 1 / p (Γ) for j=1,,2m, and h j W p 2 m m j 1 / p (Ω) for j=1,,m.

  5. (4)

    In addition, we assume proper ellipticity, i.e. the polynomials A 1 0 (x, ξ ,t) and A ˜ 2 0 (x, ξ ,t,λ)C[t] of order 2m from conditions (2a) and (2b) have exactly m roots in each half-plane C ± :={zC:±Imz>0} for all x Ω ¯ 1 and x Ω ¯ 2 , respectively, and for all ξ R n 1 {0} and λ Σ ¯ θ . Proper ellipticity allows a decomposition of the form A 1 0 (x, ξ ,t)= A 1 + (x, ξ ,t) A 1 (x, ξ ,t) with

    A 1 + ( x , ξ , t ) : = j = 1 m ( t τ j ( x , ξ ) ) and A 1 ( x , ξ , t ) : = j = m + 1 2 m ( t τ j ( x , ξ ) ) ,
    (2.2)

    where τ j denote the roots in C + (jm) and C (j>m), respectively. A similar decomposition with an additional dependence on λ also holds for A ˜ 2 0 . We remark that proper ellipticity holds automatically if n3.

Under these assumptions, we consider the inhomogeneous transmission boundary value problem

A 1 ( x , D ) u 1 = f 1 in  Ω 1 , A ˜ 2 ( x , D , λ ) u 2 = f 2 in  Ω 2 , B j ( x , D ) u = g j on  Γ ( j = 1 , , 2 m ) , C j ( x , D ) u 2 = h j on  Ω ( j = 1 , , m ) .
(2.3)

Here we have set B j (x,D)u:= ν j 1 u 1 ν j 1 u 2 where ν denotes the derivative in direction of the outer normal with respect to Ω 2 . Our main result is the following a priori estimate for solutions to (2.3). Here, a solution of (2.3) is defined as a pair ( u 1 , u 2 ) belonging to the Sobolev space W p 2 m ( Ω 1 )× W p 2 m ( Ω 2 ) for which the system (2.3) is satisfied as equality of L p -functions.

Theorem 2.2 (A priori estimate for the transmission boundary value problem)

Let Assumption  2.1 be satisfied and let u=( u 1 , u 2 ) W p 2 m ( Ω 1 )× W p 2 m ( Ω 2 ) be a solution to the transmission problem (2.3). Then there exists λ 0 >0 such that for all λ Σ ¯ θ with |λ| λ 0 the following estimates hold:

u 1 2 m , p , Ω 1 + u 2 m , p , Ω 2 + | λ | 1 / 2 u 2 0 , p , Ω 2 C ( f 1 0 , p , Ω 1 + f 2 0 , p , Ω 2 + u 1 0 , p , Ω 1 + j = 1 2 m g j 2 m j + 1 1 / p , p , Γ + j = 1 m h j 2 m m j 1 / p , p , Ω ) ,
(2.4)
u 1 2 m , p , Ω 1 + | λ | 1 / 2 u 1 m , p , Ω 1 + u 2 2 m , p , Ω 2 C ( | λ | 1 / 2 f 1 0 , p , Ω 1 + f 2 0 , p , Ω 2 + | λ | 1 / 2 u 1 0 , p , Ω 1 + j = 1 2 m g j 2 m j + 1 1 / p , p , Γ + j = 1 m h j 2 m m j 1 / p , p , Ω ) .
(2.5)

Note that with respect to g, inequality (2.4) is of elliptic type and (2.5) is of parameter-elliptic type. Due to the fact that the boundary operators C j act on u 2 , we have parameter-elliptic norms with respect to h j in both inequalities.

Remark 2.3 Our main task will be to study the problem for constant coefficient operators A 1 (D) and A ˜ 2 (D,λ) in the half-spaces R ± n without lower order terms. This simplification can be justified by performing a localization procedure, using a finite covering Ω ¯ k = 1 N U k with appropriate open sets U k , a corresponding partition of unity and perturbation results. For a detailed explanation of the localization procedure, we refer to [8], pp.151-153, but here we briefly list the types of local problems one has to deal with. If U k ¯ Ω i , one faces a local elliptic (i=1) or parameter-elliptic (i=2) operator in the whole space. For these situations, the estimates for u i are well-known results, see [17], Theorem 5.3.2, for the elliptic and [8], Proposition 2.5, for the parameter-elliptic case. If U k Ω, the local problem is a standard boundary value problem in the half-space and the desired estimate is contained in [8], Proposition 2.6. It remains to consider the case where U k intersects both Ω 1 and Ω 2 , and in the sequel we restrict our considerations to the corresponding local model problem. This reads

A 1 ( D ) u 1 = f 1 in  R + n , A ˜ 2 ( D , λ ) u ˜ 2 = f ˜ 2 in  R n , D n j 1 ( u 1 u 2 ) = g j on  R n 1 ( j = 1 , , 2 m ) .
(2.6)

The reflection τ n : R n R n , x( x , x n ) will be useful to treat problem (2.6). Therefore, we will use the notation A 2 (ξ,λ):= A ˜ 2 ( τ n (ξ),λ)= A ˜ 2 ( ξ , ξ n ,λ) for the symbol of the reflected operator, which is parameter-elliptic in R + n . We set u 2 (x):= u ˜ 2 ( τ n (x)) and f 2 (x):= f ˜ 2 ( τ n (x)).

By this substitution, we may rewrite (2.6) as a system in the half-space R + n :

A ( D , λ ) u = f  in  R + n , D n j 1 ( u 1 + ( 1 ) j u 2 ) = g j ( j = 1 , , 2 m )  on  R n 1 .
(2.7)

Here we have set

A ( D , λ ) : = ( A 1 ( D ) 0 0 A 2 ( D , λ ) ) , u : = ( u 1 u 2 ) , f : = ( f 1 f 2 ) .

Remark 2.4 We see that the determinant of the principal symbol det( A 0 (ξ,λ))=det(A(ξ,λ))= A 1 (ξ) A 2 (ξ,λ) vanishes at the points (0,λ) R n ×(0,). Hence the standard theory for parameter-elliptic systems is not applicable in this case. Due to continuity and homogeneity of the principal symbols we have the estimate

| A 1 ( ξ ) A 2 ( ξ , λ ) | C | ξ | 2 m ( | λ | + | ξ | 2 m )
(2.8)

with a constant C>0. Operators whose principal symbols allow an estimate of the form (2.8) are also called N-elliptic with parameter. Here the ‘N’ stands for the Newton polygon which is related to the principal symbol. In the case of (2.8), the Newton polygon is not regular, and therefore this equation is not covered by the results on N-ellipticity as in [10].

Remark 2.5 The boundary conditions in (2.6) are called canonical transmission conditions. In the case g j =0, they are equivalent to the condition U W p 2 m ( R n ) for

U ( x , x n ) :={ u 1 ( x , x n ) ( x n 0 ) , u ˜ 2 ( x , x n ) ( x n < 0 ) .

Note that in (2.6) the number of conditions equals the order of the operator, in contrast to boundary value problems. We will show in Lemma 3.1 below that the ODE system corresponding to the transmission problem (2.6) is uniquely solvable. This is an analogue of the Dirichlet boundary conditions which are absolutely elliptic, i.e., for every properly elliptic operator the Dirichlet boundary value problem satisfies the Shapiro-Lopatinskii condition.

3 Fundamental solutions and solution operators

To represent the solution in terms of fundamental solutions, we start with the observation that the ODE system obtained from (2.6) by partial Fourier transform is uniquely solvable. This is the analogue of the Shapiro-Lopatinskii condition for transmission problems. For detailed discussions of this condition for boundary value problems, we refer to [18], Section 6.2 and [19], Chapter 11. The assertion of the following lemma is formulated for our situation of one elliptic and one parameter-elliptic operator but of course it also holds in the cases when both operators are of the same type.

To simplify our notation, we define q:= λ 1 / 2 m and consider the differential operator A ˜ 2 (D,q)= | α | + k 2 m a α k ( 2 ) q k D α with q Σ ¯ := Σ ¯ θ / ( 2 m ) .

Lemma 3.1 Suppose the operators A 1 (x,D) and A ˜ 2 (x,D,q) are elliptic and parameter-elliptic in Σ ¯ , respectively. Fix x 0 Ω, ξ R n 1 {0}, q Σ ¯ and let h j C (j=1,,2m). Then the ODE problem

A 1 0 ( x 0 , ξ , D n ) u 1 = 0 ( x n > 0 ) , A ˜ 2 0 ( x 0 , ξ , D n , q ) u ˜ 2 = 0 ( x n < 0 ) , D n j 1 ( u 1 u ˜ 2 ) | x n = 0 = h j ( j = 1 , , 2 m ) , u 1 ( x n ) 0 ( x n ) , u ˜ 2 ( x n ) 0 ( x n )
(3.1)

admits a unique solution.

Proof In the sequel, we do not write down the dependence of the polynomials and their roots on x 0 explicitly and fix ξ R n 1 {0} as well as q Σ ¯ . We decompose A 1 0 ( ξ ,t) and A ˜ 0 2 ( ξ ,t,q) as indicated in (2.2) into A i ± 0 ( ξ ,t). Let M 1 denote the m-dimensional space of stable solutions to

A 1 0 ( ξ , D n ) v=0( x n >0),v0( x n )

and let M 2 denote the m-dimensional space of stable solutions to

A ˜ 2 0 ( ξ , D n , q ) w=0( x n <0),w0( x n ).

Let B 1 :={ v 1 ,, v m } and B 2 :={ w 1 ,, w m } be a basis of M 1 and M 2 , respectively. Then B:= B 1 B 2 is obviously a subset of the 2m-dimensional space of solutions to the equation

P ( ξ , D n , q ) u( x n ):= A 1 + 0 ( ξ , D n ) A 2 0 ( ξ , D n , q ) u( x n )=0on R
(3.2)

and B is linearly independent: Suppose there are nontrivial α j , β j C (j=1,,m) with

j = 1 m α j v j = j = 1 m β j w j .

Then (3.2) would possess a solution which is bounded on the entire real line, which contradicts the fact that the polynomial P( ξ ,t,q) has only roots with nonzero imaginary part. Hence B is a fundamental system to (3.2) and the determinant of the Wronskian matrix W( x n ) is nonzero:

detW( x n )=det( v 1 ( x n ) w m ( x n ) D n 2 m 1 v 1 ( x n ) D n 2 m 1 w m ( x n ) )0( x n R).
(3.3)

Now suppose that (v,w) is a solution to (3.1). Then there exist constants α i , β i C for i=1,,m, such that

v( x n )= j = 1 m α j v j ( x n )andw( x n )= j = 1 m β j w j ( x n ).

If we plug in this approach into the transmission conditions, we obtain the system of linear equations to determine α j and β j :

( v 1 ( 0 ) w m ( 0 ) D n 2 m 1 v 1 ( 0 ) D n 2 m 1 w m ( 0 ) )( α 1 β m )=( h 1 h 2 m ).

From (3.3) it now follows that the coefficients exist and are uniquely determined, which proves the assertion. □

From now on, we restrict ourselves to the model problem (2.7) which is the only non-standard step in the proof of the main theorem, see Remark 2.3. We first consider the case f=0 in (2.7), i.e. we study

A ( D , q ) u = 0 in  R + n , B ( D n ) u = g on  R n 1 .
(3.4)

Here u= ( u 1 , u 2 ) , g= ( g 1 , , g 2 m ) ,

A(D,q)=( A 1 ( D ) 0 0 A 2 ( D , q ) ),B( D n )=( B ( 1 , 1 ) ( D n ) B ( 1 , 2 ) ( D n ) B ( 2 , 1 ) ( D n ) B ( 2 , 2 ) ( D n ) )

with

B ( 1 , 1 ) ( D n ) : = ( D n j 1 ) j = 1 , , m , B ( 1 , 2 ) ( D n ) : = ( ( 1 ) j D n j 1 ) j = 1 , , m , B ( 2 , 1 ) ( D n ) : = ( D n j 1 ) j = m + 1 , , 2 m , B ( 2 , 2 ) ( D n ) : = ( ( 1 ) j D n j 1 ) j = m + 1 , , 2 m .

Note that B ( 1 , 1 ) ( D n ) and B ( 2 , 1 ) ( D n ) are also called generalized Dirichlet and Neumann conditions, respectively.

Due to Lemma 3.1, the ODE system corresponding to (3.4) is uniquely solvable. The main step in the proof of Theorem 2.2 will be to find a priori estimates for the fundamental solutions of this ODE system. In the following, I k stands for the (k×k)-dimensional unit matrix.

Definition 3.2 The fundamental solution

ω: ( R n 1 { 0 } ) ×(0,)× Σ ¯ C 2 × 2 m , ( ξ , x n , q ) ω ( ξ , x n , q )

is defined as the unique solution of the ODE system (in x n )

A ( ξ , D n , q ) ω ( ξ , x n , q ) = 0 ( x n > 0 ) , B ( D n ) ω ( ξ , x n , q ) | x n = 0 = I 2 m , ω ( ξ , x n , q ) 0 ( x n ) .

Following an idea of Leonid Volevich [20], we represent the solutions in a specific way. For this, we consider the elliptic boundary value problem ( A 1 (D), B ( 1 , 1 ) ( D n )) and the parameter-elliptic boundary value problem ( A 2 (D,q), B ( 2 , 2 ) ( D n )) separately. It is well known that the (generalized) Dirichlet and Neumann boundary conditions are absolutely elliptic, hence the Shapiro-Lopatinskii condition holds for both subproblems. We will call the canonical basis for these boundary value problems the basic solutions Y ( 1 ) and Y ( 2 ) . More precisely, we define the following.

Definition 3.3 We define the basic solution

Y ( 1 ) : ( R n 1 { 0 } ) ×(0,) C 1 × m , ( ξ , x n ) Y ( 1 ) ( ξ , x n ) ,

as the unique solution of the ODE system

A 1 ( ξ , D n ) Y ( 1 ) ( ξ , x n ) = 0 ( x n > 0 ) , B ( 1 , 1 ) ( D n ) Y ( 1 ) ( ξ , x n ) | x n = 0 = I m , Y ( 1 ) ( ξ , x n ) 0 ( x n ) .
(3.5)

Analogously, the basic solution

Y ( 2 ) : ( R n 1 { 0 } ) ×(0,)× Σ ¯ C 1 × m , ( ξ , x n , q ) Y ( 2 ) ( ξ , x n , q ) ,

is defined as the unique solution of the ODE system

A 2 ( ξ , D n , q ) Y ( 2 ) ( ξ , x n , q ) = 0 ( x n > 0 ) , B ( 2 , 2 ) ( D n ) Y ( 2 ) ( ξ , x n , q ) | x n = 0 = I m , Y ( 2 ) ( ξ , x n , q ) 0 ( x n ) .
(3.6)

We set

Y ( ξ , x n , q ) = ( Y k ( j ) ( ξ , x n , q ) ) j = 1 , 2 k = 1 , , 2 m :=( Y ( 1 ) ( ξ , x n ) 0 0 Y ( 2 ) ( ξ , x n , q ) ).

The advantage of the basic solutions Y ( 1 ) , Y ( 2 ) lies in the fact that classical (parameter-)elliptic estimates are easily available for them. We have to compare these solutions with the fundamental solution ω. Let j{1,,2m}. As the function ω 1 j is a solution of A 1 ( ξ , D n ) ω 1 j =0 ( x n >0), it can be written as a linear combination of the basic solutions. Therefore, we can write

ω 1 j ( ξ , x n , q ) = k = 1 m Y k ( 1 ) ( ξ , x n , q ) ψ k j ( ξ , q )

with unknown coefficients ψ k j . The analogous representation holds for ω 2 j . In matrix notation, we obtain

ω ( ξ , x n , q ) =Y ( ξ , x n , q ) Ψ ( ξ , q )

with Ψ( ξ ,q)= ( ψ k j ( ξ , q ) ) k , j = 1 , , 2 m . By the definition of the fundamental solution, we have

I 2 m = B ( D n ) ω ( ξ , x n , q ) | x n = 0 = B ( D n ) Y ( ξ , x n , q ) | x n = 0 Ψ ( ξ , q ) .

Therefore,

Ψ ( ξ , q ) = ( I m B ( 1 , 2 ) ( D n ) Y ( 2 ) ( ξ , x n , q ) | x n = 0 B ( 2 , 1 ) ( D n ) Y ( 1 ) ( ξ , x n ) | x n = 0 I m ) 1 .
(3.7)

Remark 3.4 Due to the unique solvability of equations (3.5) and (3.6), we have for ( ξ ,q)( R n 1 {0})× Σ ¯ the following scaling properties for all r>0:

Y ( 1 ) ( ξ r , r x n ) = Y ( 1 ) ( ξ , x n ) Δ 1 ( r ) , Y ( 2 ) ( ξ r , r x n , q r ) = Y ( 2 ) ( ξ , x n , q ) Δ 2 ( r ) ,

where we used the abbreviations

Δ 1 ( r ) : = diag ( 1 , r , , r m 1 ) , Δ 2 ( r ) : = diag ( r m , , r 2 m 1 ) = r m Δ 1 ( r ) .

(See also (3.10) below for an explicit representation of Y ( 1 ) and Y ( 2 ) .) We will apply this with r:=| ξ | for Y ( 1 ) and r:=| ξ |+|q| for Y ( 2 ) . Note that these scaling properties also yield the identities

B ( 2 , 1 ) ( D n ) Y ( 1 ) ( ξ , 0 ) = Δ 2 ( r ) B ( 2 , 1 ) ( D n ) Y ( 1 ) ( ξ r , 0 ) Δ 1 ( r ) 1 , B ( 1 , 2 ) ( D n ) Y ( 2 ) ( ξ , 0 , q ) = Δ 1 ( r ) B ( 1 , 2 ) ( D n ) Y ( 2 ) ( ξ r , 0 , q r ) Δ 2 ( r ) 1 .
(3.8)

We summarize the representation of the solution in form of solution operators:

Lemma 3.5 Let g j = 1 2 m W p 2 m j + 1 1 / p ( R n 1 ), and let u W p 2 m ( R + n ) be a solution of (3.4). Let g ˜ j = 1 2 m W p 2 m j + 1 ( R + n ) be an extension of g to the half-space. Then u has the form

u= T 1 g ˜ + T 2 ( n g ˜ ),

where n := x n and where the solution operators T 1 and T 2 are given by

( T 1 φ ) ( x , x n ) = 0 ( F ) 1 ( n Y ) ( ξ , x n + y n , q ) Ψ ( ξ , q ) ( F φ ) ( ξ , y n ) d y n , ( T 2 φ ) ( x , x n ) = 0 ( F ) 1 Y ( ξ , x n + y n , q ) Ψ ( ξ , q ) ( F φ ) ( ξ , y n ) d y n .

Here the basic solution Y( ξ , x n ,q) is defined in Definition  3.3, and the coefficient matrix Ψ( ξ ,q) is defined in (3.7).

Proof By definition of the fundamental solution, we have u= ( F ) 1 ω(, x n )Fg. Writing this in the form

u= 0 y n [ ( F ) 1 ω ( , x n + y n ) ( F g ˜ ) ( , y n ) ] d y n ,

(‘Volevich trick’) and noting that ω( ξ , x n ,q)=Y( ξ , x n ,q)Ψ( ξ ,q), we obtain the above representation. □

Our proofs are based on the Fourier multiplier concept, see, e.g., [18]. Here a function m L ( R n ) is called an L p -Fourier multiplier if T m :S( R n ) L ( R n ), f F 1 mF (being defined on the Schwartz space S( R n )) extends to a continuous mapping T m L( L p ( R n )). We will apply Michlin’s theorem to prove the Fourier multiplier property. For this, we introduce the notion of a Michlin function.

Definition 3.6 Let M:( R n 1 {0})× Σ ¯ C k × be a matrix-valued function. Then we call M a Michlin function if M(,q) C [ n 2 ] + 1 ( R n 1 {0}) for all q Σ ¯ and if there exists a constant C>0, independent of q, γ , and ξ , such that

| ξ | | γ | | ξ γ M ( ξ , q ) | C ( ξ R n 1 { 0 } , q Σ ¯ , γ N 0 n 1  with  | γ | [ n 2 ] + 1 ) .

Remark 3.7 (a) Michlin’s theorem (see [17], Section 2.2.4) states that every Michlin function is an L p -Fourier multiplier for all p(1,).

  1. (b)

    By the product rule one immediately sees that the product of Michlin functions is a Michlin function, too.

  2. (c)

    Let M:( R n 1 {0})× Σ ¯ C k × k be a Michlin function, and let M( ξ ,q) be invertible for all ξ and q. If the norm of the inverse matrix is bounded by a constant independent of ξ and q, then also ( ξ ,q)M ( ξ , q ) 1 is a Michlin function. This follows iteratively noting that

    ξ j ξ j M ( ξ , q ) 1 =M ( ξ , q ) 1 ( ξ j ξ j M ( ξ , q ) ) M ( ξ , q ) 1 .

Now we will show that the basic solution Y as well as the coefficient matrix Ψ satisfy uniform estimates. Here and in the following, C stands for a generic constant which may vary from inequality to inequality but is independent of the variables appearing in the inequality. We will scale the functions with | ξ | and with

ρ:=ρ ( ξ , q ) := | ξ | +|q|.
(3.9)

Lemma 3.8 (a) For all N 0 and all x n >0, the function

M 1 ( ) ( ξ , x n , q ) := x n ( | ξ | 0 0 ρ ) n + 1 Y ( ξ , x n , q ) ( Δ 1 ( | ξ | ) 0 0 Δ 2 ( ρ ) )

is a Michlin function with constant independent of x n (0,).

  1. (b)

    The functions

    C 1 ( ξ , q ) : = Δ 1 ( ρ ) 1 ( B ( 1 , 2 ) ( D n ) Y ( 2 ) ) ( ξ , 0 , q ) Δ 2 ( ρ ) , C 2 ( ξ ) : = Δ 2 ( | ξ | ) 1 ( B ( 2 , 1 ) ( D n ) Y ( 1 ) ) ( ξ , 0 ) Δ 1 ( | ξ | )

are Michlin functions.

Proof We use an explicit description of the basic solutions. According to [21], Section 1, there exist polynomials (with respect to τ) N 1 ( ξ ,τ),, N m ( ξ ,τ) and N m + 1 ( ξ ,τ,q),, N 2 m ( ξ ,τ,q) such that

1 2 π i γ 1 N k ( ξ , τ ) A 1 + ( ξ , τ ) τ j 1 d τ = δ j k ( j , k = 1 , , m ) , 1 2 π i γ 2 N k ( ξ , τ , q ) A 2 + ( ξ , τ , q ) τ j 1 d τ = δ j k ( j , k = m + 1 , , 2 m )

with δ j k being the Kronecker delta symbol. Here γ 1 = γ 1 ( ξ ) is a smooth closed contour in the upper half-plane C + , depending on ξ and enclosing the m roots of the polynomial A 1 ( ξ ,) with positive imaginary part, while γ 2 = γ 2 ( ξ ,q) is a smooth closed contour in C + depending on ( ξ ,q) and enclosing the m roots of A 2 ( ξ ,,q) in C + . Moreover, N k is positively homogeneous in its arguments of degree mk for k=1,,2m while A 1 + and A 2 + are positively homogeneous in their arguments of degree m.

This leads to the following representation for the basic solutions Y ( 1 ) = ( Y k ( 1 ) ) k = 1 , , m and Y ( 2 ) = ( Y k ( 2 ) ) k = m + 1 , , 2 m :

Y k ( 1 ) = 1 2 π i γ 1 N k ( ξ , τ ) A 1 + ( ξ , τ ) e i x n τ d τ ( k = 1 , , m ) , Y k ( 2 ) = 1 2 π i γ 2 N k ( ξ , τ , q ) A 2 + ( ξ , τ , q ) e i x n τ d τ ( k = m + 1 , , 2 m ) .
(3.10)

To prove part (a), we will show that for all j N 0

x n | ξ | k j n j Y k ( 1 ) ( ξ , x n ) and x n ρ k j n j Y k ( 2 ) ( ξ , x n , q )
(3.11)

are Michlin functions. Setting j:=+1 and noting the definitions of Δ 1 and Δ 2 , this immediately implies (a). Similarly, to show (b) we have to prove that

| ξ | k j 1 n j Y k ( 1 ) ( ξ , 0 ) and ρ k j 1 n j Y k ( 2 ) ( ξ , 0 , q )
(3.12)

are Michlin functions. We will restrict ourselves to Y k ( 2 ) , the result for Y k ( 1 ) follows in the same way.

For j N 0 and k{m+1,,2m}, we substitute ττ/ρ in the integral representation (3.10) and obtain

ξ γ [ x n ρ k j n j Y k ( 2 ) ( ξ , x n , q ) ] = 1 2 π i γ 2 ( ξ , q ) ξ γ [ ρ k j N k ( ξ , τ , q ) A 2 + ( ξ , τ , q ) ] τ j x n e i x n τ d τ = 1 2 π i γ 2 ( ξ / ρ , q / ρ ) ξ γ [ ρ k j N k ( ξ , ρ τ , q ) A 2 + ( ξ , ρ τ , q ) ] ( ρ τ ) j x n e i ρ x n τ ρ d τ = 1 2 π i γ ˜ 2 ρ j ξ γ [ H k ( ξ , ρ τ , q ) ] τ j ( ρ x n ) e i ρ x n τ d τ

with H k ( ξ ,ρ,τ):= ρ k j N k ( ξ ,ρτ,q)/ A 2 + ( ξ ,ρτ,q). Note for the first equality that it is not necessary to differentiate the contour γ 2 ( ξ ,q) because it may be chosen locally independent of ξ . In the last equality, we replaced the contour γ 2 ( ξ ρ , q ρ ) by a fixed contour γ ˜ 2 which is possible by a compactness argument.

Due to the properties of N k and A 2 + , the function H k is homogeneous of degree −j in its arguments. Therefore, ξ γ H k is homogeneous of degree j| γ | in its arguments, and we obtain

ξ γ [ H k ( ξ , ρ τ , q ) ] = ρ j | γ | ( ξ γ H k ) ( ξ ρ , τ , q ρ ) .

From the fact that γ ˜ 2 may be chosen in C + and the elementary inequality t e t 1 (t0) we get

| ( ρ x n ) e i ρ x n τ | =(ρ x n ) e ρ x n Im τ 1 Im τ C

for τ γ ˜ 2 . Inserting this and the homogeneity of H k into the above representation, we see

| ξ γ [ x n ρ k j n j Y k ( 2 ) ( ξ , x n , q ) ] | C ρ j ρ j | γ | C | ξ | | γ |

which shows (3.11). In the same way, for the proof of (3.12) we set x n =0 in the above integral representation and obtain

| ξ γ [ ρ k j 1 n j Y k ( 2 ) ( ξ , 0 , q ) ] | = | 1 2 π i γ ˜ 2 ρ j ξ γ [ H k ( ξ , ρ τ , q ) ] τ j d τ | C ρ j ρ j | γ | C | ξ | | γ | .

This finishes the proof of (3.11) and (3.12) for Y k ( 2 ) . For Y k ( 1 ) , we use the substitution ττ/| ξ | in the integral representation. As indicated above, (a) and (b) are immediate consequences of (3.11) and (3.12), respectively. □

The last lemma in connection with the following result is the essential step for the proof of the a priori estimates from the main theorem.

Lemma 3.9 The functions

M 2 ( ξ , q ) : = ( Δ 1 ( | ξ | ) 1 0 0 | ξ | m Δ 1 ( ρ ) 1 ) Ψ ( ξ , q ) ( Δ 1 ( | ξ | ) 0 0 | ξ | m Δ 1 ( ρ ) ) , M ˜ 2 ( ξ , q ) : = ( | ξ | 1 I m 0 0 ρ 1 I m ) M 2 ( ξ , q ) ( | ξ | I m 0 0 ρ I m )

are Michlin functions.

Proof By Lemma 3.8(b), we have

Ψ ( ξ , q ) = ( I m Δ 1 ( ρ ) C 1 ( ξ , q ) Δ 2 ( ρ ) 1 Δ 2 ( | ξ | ) C 2 ( ξ ) Δ 1 ( | ξ | ) 1 I m ) 1

with Michlin functions C 1 and C 2 . For M 2 we obtain

M 2 ( ξ , q ) = ( I m ( | ξ | ρ ) m Δ 1 ( ρ | ξ | ) C 1 ( ξ , q ) Δ 1 ( | ξ | ρ ) C 2 ( ξ ) I m ) 1 .
(3.13)

By a homogeneity argument we see that Δ 1 (| ξ |/ρ) and ( | ξ | / ρ ) m Δ 1 (ρ/| ξ |) are Michlin functions, and therefore the matrix on the right-hand side of (3.13) is a Michlin function. In order to apply Remark 3.7(c), we have to show that the norm of M 2 ( ξ ,q) is uniformly bounded.

For this, we write M 2 ( ξ ,q) in the form of a Schur complement: For an invertible block matrix, we have

( I m A ( 1 , 2 ) A ( 2 , 1 ) I m ) 1 =( I m + A ( 1 , 2 ) S 1 A ( 2 , 1 ) A ( 1 , 2 ) S 1 S 1 A ( 2 , 1 ) S 1 )

with S:= I m A ( 2 , 1 ) A ( 1 , 2 ) . Applied to the matrix M 2 , we obtain

M 2 ( ξ , q ) =( I m + ( | ξ | ρ ) m Δ 1 ( ρ | ξ | ) C 1 S 1 Δ 1 ( | ξ | ρ ) C 2 ( | ξ | ρ ) m Δ 1 ( ρ | ξ | ) C 1 S 1 S 1 Δ 1 ( | ξ | ρ ) C 2 S 1 )
(3.14)

with

S ( ξ , q ) := I m ( | ξ | ρ ) m Δ 1 ( | ξ | ρ ) C 2 ( ξ ) Δ 1 ( ρ | ξ | ) C 1 ( ξ , q ) .
(3.15)

By (3.8), the matrices C 1 and C 2 and, consequently, the matrix M 2 are homogeneous of degree 0 in their arguments. Thus we can write S in the form

S ( ξ , q ) =S ( ξ | ξ | , q | ξ | ) ( ξ R n 1 { 0 } , q Σ ¯ ) .

We set η := ξ /| ξ | and

Λ:= ρ | ξ | = | ξ | + | q | | ξ | =1+ | q | | ξ |

and write S as

S ( ξ , q ) = I m Λ m Δ 1 ( 1 Λ ) C 2 ( η ) Δ 1 (Λ) C 1 ( η , q | ξ | ) .

The matrices C 1 , C 2 , and Δ 1 (1/Λ) are bounded for all ξ R n 1 {0} and q Σ ¯ . By | Δ 1 (Λ)|C Λ m 1 for all Λ1, we see that there exists a Λ 0 >1 such that

| Λ m Δ 1 ( 1 Λ ) C 2 ( η ) Δ 1 ( Λ ) C 1 ( η , q | ξ | ) | 1 2

holds for all ξ R n 1 {0} and q Σ ¯ with |q| Λ 0 | ξ |. For these ξ and q, a Neumann series argument shows that the norm of S 1 ( ξ ,q) is bounded by 2.

For ξ R n 1 {0} and q Σ ¯ with |q| Λ 0 | ξ |, the tuple ( ξ /| ξ |,q/| ξ |) belongs to the compact set {( η , q ˜ ):| η |=1, q ˜ Σ ¯ ,| q ˜ | Λ 0 }. Now we use the fact that for all ξ R n 1 {0} and q Σ ¯ , the matrix B( D n )Y( ξ ,0,q) is invertible, and therefore the matrix on the right-hand side of (3.13) is invertible, too. This yields the invertibility of S, and by continuity the inverse matrix S 1 ( ξ ,q) is bounded for these ξ and q.

Therefore, we have seen that | S 1 ( ξ ,q)|C holds for all ξ R n 1 {0} and q Σ ¯ . From the explicit description of M 2 ( ξ ,q) in (3.14) and the uniform boundedness of the other coefficients in (3.14), we see that | M 2 ( ξ ,q)|C holds for all ξ and q. By Remark 3.7(c), M 2 is a Michlin function.

The above proof also shows that the modification M ˜ 2 is a Michlin function. Note that S( ξ ,q) remains unchanged and that we obtain an additional factor ρ/| ξ | in the right upper corner which does not affect the boundedness. □

4 Proof of the a priori estimate

In this section, we will investigate the mapping properties of the solution operators T 1 , T 2 introduced in Lemma 3.5. As above, let ρ:=| ξ |+|q|. In the following, we will use the abbreviations D :=i( x 1 x n 1 ) and L( D ,q):= ( F ) 1 L( ξ ,q) F . Based on Lemma 3.8 and 3.9 and on the continuity of the Hilbert transform, it is not difficult to obtain the following result.

Lemma 4.1 (a) Let

L 1 ( ξ , q ) : = ( ρ m | ξ | m 0 0 ρ 2 m ) , L 2 ( ξ , q ) : = ( ρ m | ξ | m Δ 1 ( | ξ | ) 1 0 0 ρ m Δ 1 ( ρ ) 1 ) .

Then for all φS ( R + n ) 2 m and all N 0 we have

L 1 ( D , q ) ( | D | 0 0 ( | D | + | q | ) ) n T 1 φ L p ( R + n ) C L 2 ( D , q ) φ L p ( R + n ) .

The same holds when L 1 and L 2 are replaced by L 1 ( 0 ) := | ξ | m ρ m L 1 and L 2 ( 0 ) := | ξ | m ρ m L 2 , respectively.

  1. (b)

    Let

    L ˜ 2 ( ξ , q ) :=( ρ m | ξ | m 1 Δ 1 ( | ξ | ) 1 0 0 ρ m 1 Δ 1 ( ρ ) 1 ).

Then for all φS ( R + n ) 2 m and all N 0 we have

L 1 ( D , q ) ( | D | 0 0 ( | D | + | q | ) ) n T 2 φ L p ( R + n ) C L ˜ 2 ( D , q ) φ L p ( R + n ) .

Proof (a) For fixed N 0 , let φ ˜ := L 2 ( D ,q)φ and

u ˜ := L 1 ( D , q ) ( | D | 0 0 ( | D | + | q | ) ) n T 1 φ.

We have to show that u ˜ L p ( R + n ) C φ ˜ L p ( R + n ) . For this, we write

L 1 ( ξ , q ) ( | ξ | 0 0 ρ ) n + 1 Y ( ξ , x n , q ) Ψ ( ξ , q ) L 2 ( ξ , q ) 1 = x n 1 M 1 ( ) ( ξ , x n , q ) ( ρ m | ξ | m I 2 m ) M 2 ( ξ , q ) ( ρ m | ξ | m I 2 m ) = x n 1 M 1 ( ) ( ξ , x n , q ) M 2 ( ξ , q ) .

Inserting this into the definition of the solution operator, we obtain

u ˜ = 0 1 x n + y n ( F ) 1 M 1 ( ) ( ξ , x n + y n , q ) M 2 ( ξ , q ) F φ ˜ ( ξ , y n ) d y n .

Therefore,

u ˜ L p ( R + n ) p = 0 F u ˜ ( , x n ) L p ( R n 1 ) p d x n 0 [ 0 1 x n + y n ( F ) 1 M 1 ( ) ( ξ , x n + y n , q ) × M 2 ( ξ , q ) F φ ˜ ( ξ , y n ) L p ( R n 1 ) d y n ] p d x n C 0 [ 0 1 x n + y n φ ˜ ( , x n ) L p ( R n 1 ) d y n ] p d x n C 0 φ ˜ ( , y n ) L p ( R n 1 ) p d y n = C φ ˜ L p ( R + n ) p .

Here we used the fact that M 1 ( ) and M 2 are Michlin functions and therefore Fourier multipliers and that the (one-sided) Hilbert transform

ϕHϕ,(Hϕ)( x n ):= 0 ϕ ( y n ) x n + y n d y n

induces a bounded operator in L p ((0,)) for every p(1,).

This shows the first statement in (a). Obviously, the uniform estimate also holds in the case when L 1 and L 2 are multiplied with the same factor, as this factor cancels out.

The proof of (b) follows exactly in the same way with M 2 being replaced by M ˜ 2 from Lemma 3.5. □

The next result shows the key estimate for the solution of (3.4).

Theorem 4.2 Let u ( W p 2 m ( R + n ) ) 2 be a solution of A(D,q)u=0, B( D n )u=g with g j = 1 2 m W p 2 m j + 1 1 / p ( R n 1 ). Let g ˜ j = 1 2 m W p 2 m j + 1 ( R + n ) be an extension of g to the half-space. Let q 0 >0. Then for all q Σ ¯ with |q| q 0 , the inequalities

| u 1 | 2 m , p , R + n + | u 2 | m , p , R + n + | q | m | u 2 | 0 , p , R + n C j = 1 2 m g ˜ j 2 m j + 1 , p , R + n , | u 1 | 2 m , p , R + n + | q | m | u 1 | m , p , R + n + u 2 2 m , p , R + n C j = 1 2 m g ˜ j 2 m j + 1 , p , R + n

hold.

Proof In this proof, we will write := L p ( R + n ) . We use the equivalences

| φ | k , p , R + n = 0 k | D | k n φ ,φ k , p , R + n = 0 k ( | D | + | q | ) k n φ

which can easily be seen by a Michlin type argument. With this, we get

| u 1 | 2 m , p , R + n + | q | m | u 1 | m , p , R + n + u 2 2 m , p , R + n = 0 2 m | D | 2 m n u 1 + = 0 m | q | m | D | m n u 1 + = 0 2 m ( | D | + | q | ) 2 m n u 2 C = 0 2 m L 1 ( D , q ) ( | D | 0 0 ( | D | + | q | ) ) n u

and

L 2 ( D , q ) g ˜ j = 1 m ( | D | + | q | ) m | D | m j + 1 g ˜ j + j = m + 1 2 m ( | D | + | q | ) 2 m j + 1 g ˜ j j = 1 m | g ˜ j | 2 m j + 1 , p , R + n + | q | m j = 1 m | g ˜ j | m j + 1 , p , R + n + j = m + 1 2 m g ˜ j 2 m j + 1 , p , R + n .

Let u ( W p 2 m ( R + n ) ) 2 be a solution of A(D,q)u=0, B( D n )u=g, and let g ˜ be an extension of g. By a density argument, we may assume that u ( S ( R + n ) ) 2 . By Lemma 3.5, we have u= T 1 g ˜ + T 2 ( n g ˜ ). Applying Lemma 4.1, we get

| u 1 | 2 m , p , R + n + | q | m | u 1 | m , p , R + n + u 2 2 m , p , R + n C [ j = 1 m ( | g ˜ j | 2 m j + 1 , p , R + n + | q | m | g ˜ j | m j + 1 , p , R + n )
(4.1)
+ j = m + 1 2 m g ˜ j 2 m j + 1 , p , R + n + j = 1 m | n g ˜ j | 2 m j , p , R + n + | q | m | n g ˜ j | m j , p , R + n + j = m + 1 2 m n g ˜ j 2 m j , p , R + n ] C [ j = 1 m ( g ˜ j 2 m j + 1 , p , R + n + | q | m g ˜ j m j + 1 , p , R + n )
(4.2)
+ j = m + 1 2 m g ˜ j 2 m j + 1 , p , R + n ].
(4.3)

Inserting the inequality

g ˜ j 2 m j + 1 , p , R + n + | q | m g ˜ j m j + 1 , p , R + n C g ˜ j 2 m j + 1 , p , R + n

into the right-hand side, we see that

| u 1 | 2 m , p , R + n + | q | m | u 1 | m , p , R + n + u 2 2 m , p , R + n C j = 1 2 m g ˜ j 2 m j + 1 , p , R + n .
(4.4)

On the other hand, inserting the inequality

g ˜ j 2 m j + 1 , p , R + n + | q | m g ˜ j m j + 1 , p , R + n C | q | m g ˜ j 2 m j + 1 , p , R + n

(which holds for all |q| q 0 with a constant C depending on q 0 ), we get in particular

u 2 2 m , p , R + n C | q | m j = 1 2 m g ˜ j 2 m j + 1 , p , R + n .

Dividing by | q | m , we see that this implies

| u 2 | m , p , R + n + | q | m u 2 C j = 1 2 m g ˜ j 2 m j + 1 , p , R + n .
(4.5)

In the same way as above, we can apply Lemma 4.1 with L 1 ( 0 ) and L 2 ( 0 ) instead of L 1 and L 2 , respectively. We see that

| u 1 | 2 m , p , R + n = 0 2 m | D | 2 m n u 1 C = 0 2 m L 1 ( 0 ) ( D , q ) ( | D | 0 0 ( | D | + | q | ) ) n u C ( L 2 ( 0 ) ( D , q ) g ˜ + L ˜ 2 ( 0 ) ( D , q ) n g ˜ ) C ( j = 1 m | D | 2 m j + 1 g ˜ j + j = m + 1 2 m | D | m ( | D | + | q | ) m j + 1 g ˜ j + j = 1 m | D | 2 m j n g ˜ j + j = m + 1 2 m | D | m ( | D | + | q | ) m j n g ˜ j ) .

With the inequality

| D | m ( | D | + | q | ) m j + 1 g ˜ j C | D | 2 m j + 1 g ˜ j (j=m+1,,2m)

this gives

| u 1 | 2 m , p , R + n C j = 1 2 m g ˜ j 2 m j + 1 , p , R + n .

This and equations (4.4) and (4.5) yield the statements of the theorem. □

Now we can consider the problem A(D,q)u=f, B( D n )u=g in the half-space. As mentioned in Remark 2.3, this finishes the proof of the main theorem.

Theorem 4.3 Let u ( W p 2 m ( R + n ) ) 2 be a solution of A(D,q)u=f, B( D n )u=g with f ( L p ( R + n ) ) 2 and g j = 1 2 m W p 2 m j + 1 1 / p ( R n 1 ). Let q 0 >0. Then for all q Σ ¯ with |q| q 0 the following a priori estimates hold:

u 1 2 m , p , R + n + u 2 m , p , R + n + | q | m u 2 0 , p , R + n C ( f 1 0 , p , R + n + f 2 0 , p , R + n + j = 1 2 m g j 2 m j + 1 1 / p , p , R n 1 + u 1 0 , p , R + n ) ,
(4.6)
u 1 2 m , p , R + n + | q | m u 1 m , p , R + n + u 2 2 m , p , R + n C ( | q | m f 1 0 , p , R + n + f 2 0 , p , R + n + j = 1 2 m g j 2 m j + 1 1 / p , p , R n 1 + | q | m u 1 0 , p , R + n ) .
(4.7)

Proof (i) We start the proof with some preliminary remarks. Let r + :φφ | R + n be the restriction operator from R n to R + n . Then r + is a retraction from W p k ( R n ) to W p k ( R + n ) for every k N 0 , and there exists a co-retraction (independent of k), i.e. a total extension operator e + L( W p k ( R + n ), W p k ( R n )) satisfying r + e + = id W p k ( R + n ) for all k (see [22], Theorem 5.21).

For every j{0,,k1}, the trace operator to the boundary γ j u:= n j u | R n 1 is a bounded operator from W p k ( R + n ) to W p k 1 / p ( R n 1 ). This holds both with respect to the parameter-independent norms and the parameter-dependent norms . For the latter, we refer to [8], Proposition 2.2. There exists a parameter-dependent extension operator E q L( W p k 1 / p ( R n 1 ), W p k ( R + n )) which satisfies γ 0 E q = id W p k 1 / p ( R n 1 ) and whose operator norm with respect to the parameter-dependent norms is bounded by a constant independent of q for all q Σ ¯ with |q| q 0 (see, e.g., [8], Proposition 2.3). In particular, we will consider E 1 which is a parameter-independent continuous extension operator.

Let ψ C ( R n ) with 0ψ1, ψ(ξ)=0 for |ξ|1 and ψ(ξ)=1 for |ξ|2. Then a simple application of Michlin’s theorem shows that R 1 (D):= F 1 ψ(ξ) A 1 1 (ξ)F induces a bounded linear operator R 1 (D)L( W p k ( R n ), W p k + 2 m ( R n )) for all k N 0 . Due to the compact support of 1ψ, the related operator (1ψ)(D) belongs to L( L p ( R n ), W p k ( R n )) for all k N 0 . Note that A 1 (D) and ψ(D) commute due to A 1 (D)ψ(D)= F 1 A 1 (ξ)ψ(ξ)F.

  1. (ii)

    Let u ( W p 2 m ( R + n ) ) 2 be a solution of A(D)u=f, B(D)u=g, and let q Σ ¯ with |q| q 0 . We define f ˜ 1 := A 1 (D) e + u 1 . Then f ˜ 1 L p ( R n ) and r + f ˜ 1 = A 1 (D) r + e + u 1 = A 1 (D) u 1 = f 1 . For

    v 1 := r + [ ( 1 ψ ) ( D ) e + u 1 + R 1 ( D ) f ˜ 1 ] ,

we obtain v 1 W p 2 m ( R + n ) and

A 1 ( D ) v 1 = r + A 1 ( D ) ( 1 ψ ) ( D ) e + u 1 + r + A 1 ( D ) R 1 ( D ) f ˜ 1 = r + ( 1 ψ ) ( D ) A 1 ( D ) e + u 1 + r + ( A 1 ψ A 1 1 ) ( D ) f ˜ 1 = r + ( 1 ψ ) ( D ) f ˜ 1 + r + ψ ( D ) f ˜ 1 = r + f ˜ 1 = f 1 .

By the continuity of the involved operators, we have

v 1 2 m , p , R + n C ( u 1 0 , p , R + n + f 1 0 , p , R + n ) .
(4.8)
  1. (iii)

    Similarly, we set v 2 := r + A 2 ( D , q ) 1 e + f 2 . It is well known (or easily seen by Michlin’s theorem) that v 2 W p 2 m ( R + n ) with A 2 (D,q) v 2 = f 2 and

    v 2 2 m , p , R + n C f 2 0 , p , R + n .
    (4.9)
  2. (iv)

    We define v:= ( v 1 , v 2 ) ( W p 2 m ( R + n ) ) 2 and w:=uv. Then w is a solution of A(D)w=0, B(D)w=gB(D)v. Applying the parameter-independent extension operator E 1 to every component of g, we define g ˜ := E 1 g j = 1 2 m W p 2 m j + 1 ( R + n ). An extension h ˜ of B(D)v is given by omitting the trace to the boundary. Note that h ˜ j = n j 1 v 1 ± n j 1 v 2 .

For the left-hand side of (4.6), we remark that for w= ( w 1 , w 2 ) we have

w 1 2 m , p , R + n C ( | w 1 | 2 m , p , R + n + w 1 0 , p , R + n ) .

By Theorem 4.3, we obtain

w 1 2 m , p , R + n + w 2 m , p , R + n + | q | m w 2 0 , p , R + n C ( f 0 , p , R + n + j = 1 2 m g ˜ j + n j 1 ( v 1 ± v 2 ) 2 m j + 1 , p , R + n + w 1 0 , p , R + n ) .
(4.10)

From (4.8) we see that

n j 1 v 1 2 m j + 1 , p , R + n C v 1 2 m , p , R + n C ( u 1 0 , p , R + n + f 1 0 , p , R + n ) .

For v 2 we obtain n j 1 v 2 2 m j + 1 , p , R + n C f 2 0 , p , R + n in the same way from (4.9). Inserting this into (4.10), we obtain the first inequality (4.6) of the theorem.

  1. (v)

    The proof of (4.7) follows the same lines. However, here we start with the refined estimate (4.3). For the left-hand side of (4.7), we note that

    u 1 2 m , p , R + n + | q | m u 1 m , p , R + n C ( | u 1 | 2 m , p , R + n + | q | m | u 1 | m , p , R p n + | q | m u 1 0 , p , R + n ) .

Now we define g ˜ := E q g with the parameter-dependent extension operator E q from part (i). Then the term on the right-hand side of (4.3) equals

j = 1 m g ˜ j + n j 1 ( v 1 ± v 2 ) 2 m j + 1 , p , R + n + | q | m j = 1 m g ˜ j + n j 1 ( v 1 ± v 2 ) m j + 1 , p , R + n + j = m + 1 2 m g ˜ j + n j 1 ( v 1 ± v 2 ) 2 m j + 1 , p , R + n .
(4.11)

For j=1,,m, we can estimate

g ˜ j 2 m j + 1 , p , R + n + | q | m g ˜ j m j + 1 , p , R + n C g ˜ j 2 m j + 1 , p , R + n C g j 2 m j + 1 1 / p , p , R n 1 .

Concerning the terms involving v 1 , we use

n j 1 v 1 2 m j + 1 , p , R + n + | q | m n j 1 v 1 m j + 1 , p , R + n v 1 2 m , p , R + n + | q | m v 1 m , p , R + n C | q | m v 1 2 m , p , R + n C | q | m ( u 1 0 , p , R + n + f 1 0 , p , R + n )

for j=1,,m and

n j 1 v 1 2 m j + 1 , p , R + n C | q | m v 1 2 m , p , R + n C | q | m ( u 1 0 , p , R + n + f 1 0 , p , R + n )

for j=m+1,,2m. Finally, the terms involving v 2 can be estimated by

n j 1 v 2 2 m j + 1 , p , R + n C v 2 2 m , p , R + n C f 2 0 , p , R + n .

So we see that all terms in (4.11) can be estimated by the right-hand side of (4.7), and the proof of (4.7) is finished. □

Remark 4.4 (a) The estimate (2.5) does not imply uniqueness of a solution to (2.3) because the elliptic part u 1 of the solution appears in a norm of lower order on the right-hand side of the estimate. Nevertheless, in bounded domains such estimates give rise to the Fredholm property of a corresponding solution operator.

  1. (b)

    For g=0 and f 1 =0, we obtain in particular

    |λ| u 2 0 , p , R + n C ( f 2 0 , p , R + n + | λ | 1 / 2 u 1 0 , p , R + n )

from (4.7). This is the basis for resolvent estimates and spectral properties of the corresponding L p -realization in the case where the Dirichlet problem for A 1 (x,D) in Ω 1 is invertible. Here we have a connection to eigenvalue problems with weights and the Calderón method as studied in, e.g., [4].

References

  1. Gebauer B: Sensitivity analysis of a parabolic-elliptic problem. Q. Appl. Math. 2007, 65: 591-604.

    Article  MathSciNet  Google Scholar 

  2. Faierman M: A transmission problem for elliptic equations involving a parameter and a weight. Glas. Mat. 2000, 35: 89-109.

    MathSciNet  Google Scholar 

  3. Faierman M: A transmission problem involving a parameter and a weight. Math. Nachr. 2001, 229: 73-89. 10.1002/1522-2616(200109)229:1<73::AID-MANA73>3.0.CO;2-P

    Article  MathSciNet  Google Scholar 

  4. Faierman M: A boundary problem for an elliptic system of differential equations involving a discontinuous weight. Math. Nachr. 2009, 282: 408-421. 10.1002/mana.200710744

    Article  MathSciNet  Google Scholar 

  5. Pyatkov SG, Abasheeva NL: Solvability of boundary value problems for operator-differential equations of mixed type. A degenerate case. Sib. Mat. Zh. 2002, 43: 678-693. English transl. in Sib. Math. J. 43, 549-561 (2002)

    Article  MathSciNet  Google Scholar 

  6. Pyatkov SG: Indefinite elliptic spectral problems. Sib. Mat. Zh. 1998, 39: 409-426. English transl. in Sib. Math. J. 39, 358-372 (1998)

    Article  MathSciNet  Google Scholar 

  7. Behrndt J: Spectral theory of elliptic differential operators with indefinite weights. Proc. R. Soc. Edinb. A 2013, 143: 21-38. 10.1017/S0308210511000965

    Article  MathSciNet  Google Scholar 

  8. Agranovich MS, Denk R, Faierman M: Weakly smooth nonselfadjoint spectral elliptic boundary problems. Math. Top. 14. In Spectral Theory, Microlocal Analysis, Singular Manifolds. Akademie Verlag, Berlin; 1997:138-199.

    Google Scholar 

  9. Shibata Y, Shimizu S:Maximal L p - L q regularity for the two-phase Stokes equations; model problems. J. Differ. Equ. 2011, 251: 373-419. 10.1016/j.jde.2011.04.005

    Article  MathSciNet  Google Scholar 

  10. Denk R, Mennicken R, Volevich L: Boundary value problems for a class of elliptic operator pencils. Integral Equ. Oper. Theory 2000, 38: 410-436. 10.1007/BF01228606

    Article  MathSciNet  Google Scholar 

  11. Agranovich MS, Vishik MI: Elliptic boundary value problems depending on a parameter. Sov. Math. Dokl. 1963, 4: 325-329.

    Google Scholar 

  12. Denk R, Volevich L: A priori estimates for a singularly perturbed mixed-order boundary value problem. Russ. J. Math. Phys. 2000, 7: 288-318.

    MathSciNet  Google Scholar 

  13. Wu J, Xu J, Zou H: On the well posedness of mathematical model for Li-ion battery systems. Methods Appl. Anal. 2006, 13: 275-298.

    MathSciNet  Google Scholar 

  14. Seger, T: Elliptic-parabolic systems with applications to lithium-ion battery models. PhD Thesis, University of Konstanz (2013)

  15. Iliev O, Latz A, Zausch J, et al.: Modeling of species and charge transport in Li-ion batteries based on non-equilibrium thermodynamics. Lecture Notes in Computer Science 6046. In Numerical Methods and Applications Edited by: Dimov I. 2011, 329-337.

    Google Scholar 

  16. Dreyer W, Gaberšček M, Guhlke C, Huth R, Jamnik J: Phase transition in a rechargeable lithium battery. Eur. J. Appl. Math. 2011, 22: 267-290. 10.1017/S0956792511000052

    Article  Google Scholar 

  17. Triebel H North-Holland Mathematical Library 18. In Interpolation Theory, Function Spaces, Differential Operators. North-Holland, Amsterdam; 1978.

    Google Scholar 

  18. Denk R, Hieber M, Prüss J: -Boundedness, Fourier multipliers and problems of elliptic and parabolic type. Mem. Am. Math. Soc. 2003., 166: Article ID 788

    Google Scholar 

  19. Wloka J: Partial Differential Equations. Cambridge University Press, Cambridge; 1987.

    Book  Google Scholar 

  20. Volevich, LR: Private communication (2004)

  21. Agmon S, Douglis A, Nirenberg L: Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions. I. Commun. Pure Appl. Math. 1959, 12: 623-727. 10.1002/cpa.3160120405

    Article  MathSciNet  Google Scholar 

  22. Adams RA, Fournier JJF Pure and Applied Mathematics 140. In Sobolev Spaces. 2nd edition. Academic Press, New York; 2003.

    Google Scholar 

  23. Denk R, Faierman M, Möller M: An elliptic boundary problem for a system involving a discontinuous weight. Manuscr. Math. 2002, 108: 289-317. 10.1007/s002290200264

    Article  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Robert Denk.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

Both authors contributed equally to this work.

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

Denk, R., Seger, T. L p -estimates for a transmission problem of mixed elliptic-parabolic type. Bound Value Probl 2014, 22 (2014). https://doi.org/10.1186/1687-2770-2014-22

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-22

Keywords