This paper studies the regularity of solutions to mixed boundary value problems of two-dimensional linear second-order elliptic equations in an angular region. When the corresponding angle in the angular region tends to 0, we consider the regularity influence of the solution of the equation under the assumption that the boundary conditions, inhomogeneous terms and coefficients have sufficient regularity. The study of mixed boundary value problems for elliptic equations in angular regions has important theoretical significance and practical application value. This type of problem has been systematically discussed in the study of elliptic equations, see [1] - [9] . Recently, in the study of hypersonic flow, we encountered a type of small-angle boundary value problem for elliptic equations. The regularity of its solution is of great significance.

This article considers the simplest type of mixed boundary value problems. Assume that the origin O = ( 0 , 0 ) is a corner point, and consider the angular region:

Ω ε : = { ( x , y ) | 0 < y < ε x , x > 0 } (2.1)

Its boundary consists of the following two rays:

l 0 : = { ( x , y ) | y = 0 , x ≥ 0 } , (2.2)

and

l ε : = { ( x , y ) | y = ε x , x ≥ 0 } , (2.3)

where 0 < ε ≪ 1 .

Assume that:

B : = { ( x , y ) | x 2 + y 2 ≤ 1 } . (2.4)

We consider the regularity of the solution to the following boundary value problem for the elliptic equation at the origin O:

u ( x , ε x ) = f ( x )     on     l ε ∩ B ; (2.5)

b u x + u y = 0     on     l 0 ∩ B ; (2.6)

Δ u = u x x + u y y = 0     in     Ω ε ∩ B ; (2.7)

where b is a fixed constant, f ∈ C m , α ( [ 0 , 1 ] ) .

Regarding this problem, we have the results of the Lieberman’s weighted Hölder norm corner point elliptic equation problem [6] . i.e. u ∈ C 1 , β ( Ω ε ∩ B ) , in particular u ∈ C 1 , β ( O ) , where β is determined by ε and b. [10] proved the case of second-order regularity. Inspired by its results, we will prove the following regularity theorem by constructing a similar barrier function.

Theorem 1. For any natural number m > 2 , there exists a sufficiently small positive number ε 0 such that when 0 < ε < ε 0 , we have u ∈ C m , α ( O ) . In fact, we will prove that there exists a polynomial of degree m, denoted as P m ( x , y ) such that the following inequality holds:

| u − P m | ≤ C r m + α ,

where C is a constant determined by ‖ f ‖ C m , α ( [ 0 , 1 ] ) , b and ‖ u ‖ L ∞ ( Ω ε ∩ ∂ B ) .

Consider the higher-order boundary conditions corresponding to (2.5)-(2.7). Taking the k-order derivative of x along l ε , we have:

∑ i = 0 k C k i ε k − i ∂ x i ∂ y k − i u ( x , ε x ) = f ( k ) ( x )     on     l ε ∩ B . (3.1)

Taking the ( k − 1 ) -order derivative of x along l 0 , we have:

b ∂ x k u ( x , 0 ) + ∂ x k − 1 ∂ y u ( x , 0 ) = 0     on     l 0 ∩ B . (3.2)

Taking the partial derivative ∂ x i ∂ y k − i − 2 of the function, where i = 0 , ⋯ , k − 2 , we have:

∂ x i + 2 ∂ y k − i − 2 u ( x , y ) + ∂ x i ∂ y k − i u ( x , y ) = 0     in     Ω ε ∩ B . (3.3)

We consider a polynomial of degree m, denoted as P m , which is defined as:

P m ( x , y ) : = ∑ k = 0 m ∑ i = 0 k P i , k − i x i y k − i .

We have the following theorem:

Theorem 2. There exists a unique polynomial of degree m, denoted as P m ( x , y ) , which satisfies for all k = 0 , ⋯ , m ,

∑ i = 0 k C k i ε k − i ∂ x i ∂ y k − i P m ( 0 , 0 ) = f ( k ) ( 0 ) ; (3.4)

b ∂ x k P m ( 0 , 0 ) + ∂ x k − 1 ∂ y P m ( 0 , 0 ) = 0 ; (3.5)

and

∂ x i + 2 ∂ y k − i − 2 P m ( 0 , 0 ) + ∂ x i ∂ y k − i P m ( 0 , 0 ) = 0 (3.6)

where i = 0 , ⋯ , k − 2 .

Proof. Obviously, P m ( x , y ) is determined by ∂ x i ∂ y k − i P m ( 0 , 0 ) , ( 0 ≤ i + j ≤ m ) . We note that for a given 0 ≤ k ≤ m , (3.4)-(3.6) form a linear system of equations of ( k + 1 ) -order with respect to ∂ x i ∂ y k − i P m ( 0 , 0 ) . According to Cramer’s law, we only need to verify that its coefficient determinant is not zero.

In fact, let ε = 0 , the matrix form corresponding to (3.4)-(3.6) is:

( 1 0 0 0 ⋯ 0 0 0 b 1 0 0 ⋯ 0 0 0 1 0 1 0 ⋯ 0 0 0 0 1 0 1 ⋯ 0 0 0 ⋮ ⋮ ⋮ ⋮ ⋱ ⋮ ⋮ ⋮ 0 0 0 0 ⋯ 1 0 1 ) ( ∂ x k P m ( 0 , 0 ) ∂ x k − 1 ∂ y P m ( 0 , 0 ) ∂ x k − 2 ∂ y 2 P m ( 0 , 0 ) ∂ x k − 3 ∂ y 3 P m ( 0 , 0 ) ⋮ ∂ y k P m ( 0 , 0 ) ) = ( f ( k ) ( 0 ) 0 0 0 ⋮ 0 )

Obviously, its coefficient determinant value is 1, and the equation has a unique solution. We note that when ε is very small, the coefficient determinant of (3.4)-(3.6) is a small perturbation of the above equation, so there must be a unique solution. The theorem is proved.

Proof of Theorem 1. Thus, for all k = 0 , ⋯ , m , we have:

∂ l ε k ( u − P m ) ( 0 , 0 ) = 0.

And thus, there is a boundary estimate:

| u − P m | ( x , ε x ) ≤ ‖ f ‖ C m , α | x | m + α ≤ 2 ‖ f ‖ C m , α r m + α . (3.7)

At the same time, for all k = 0 , ⋯ , m − 1 , we have:

∂ l 0 k ( ( b ∂ x + ∂ y ) ( u − P m ) ) ( 0 , 0 ) = 0.

And thus, there is a boundary estimate:

| ( b ∂ x + ∂ y ) ( u − P m ) ( x , 0 ) | ≤ M r m ≤ M r m − 1 + α . (3.8)

Considering the function Δ ( u − P m ) in Ω ε , then for all i = 0 , ⋯ , k − 2 and k = 2 , ⋯ , m , we have:

∂ x i ∂ y k − i − 2 ( Δ ( u − P m ) ) ( 0 , 0 ) = 0.

Thereby, we obtain an equation estimate:

| Δ ( u − P m ) ( x , y ) | ≤ M r m − 1 ≤ M r m − 2 + α . (3.9)

When considering polar coordinates:

{ r = x 2 + y 2 , θ = arctan y x ,

in this way, we have some differential relationships:

{ r x = x x 2 + y 2 = cos θ , θ x = − y x 2 + y 2 = − sin θ r , r y = y x 2 + y 2 = sin θ , θ y = x x 2 + y 2 = cos θ r , r x x = 1 x 2 + y 2 − x 2 ( x 2 + y 2 ) 3 / 2 = 1 r − cos 2 θ r , θ x x = 2 x y ( x 2 + y 2 ) 2 = 2 sin θ cos θ r 2 , r y y = 1 x 2 + y 2 − y 2 ( x 2 + y 2 ) 3 / 2 = 1 r − sin 2 θ r , θ y y = − 2 x y ( x 2 + y 2 ) 2 = − 2 sin θ cos θ r 2 , r x y = − x y ( x 2 + y 2 ) 3 / 2 = sin θ cos θ r , θ x y = − 1 x 2 + y 2 + 2 y 2 ( x 2 + y 2 ) 2 = cos 2 θ − sin 2 θ r 2 .

We assume that the barrier function has the form:

U = r m + α sin ( κ θ + φ 0 ) . (3.10)

where κ and φ 0 are undetermined constants. In this way, we first have:

U | l ε = U | θ = arctan ε = r m + α sin ( κ arctan ε + φ 0 ) . (3.11)

Then, due to

{ U x = ( m + α ) x r r m + α − 1 sin ( κ θ + φ 0 ) − κ y r r m + α − 1 cos ( κ θ + φ 0 ) ; U y = ( m + α ) y r r m + α − 1 sin ( κ θ + φ 0 ) + κ x r r m + α − 1 cos ( κ θ + φ 0 ) .

We have the expression for the boundary oblique derivative:

b U x + U y | l 0 = b U x + U y | θ = 0 = b ( m + α ) r m + α − 1 sin φ 0 + κ r m + α − 1 cos φ 0 = r m + α − 1 ( b ( m + α ) sin φ 0 + κ cos φ 0 ) = κ cos φ 0 r m + α − 1 ( b ( m + α ) κ tan φ 0 + 1 ) . (3.12)

Finally, applying the Laplace operator Δ, we have:

Δ U = U r r + U r r + U θ θ r 2 = ( ( m + α ) 2 − κ 2 ) r m + α − 2 sin ( κ θ + φ 0 ) . (3.13)

Now, we set κ and φ 0 . On the one hand, let κ > 0 be sufficiently large and satisfy:

( m + α ) 2 − κ 2 ≤ − 1 ; (3.14)

and

| b ( m + α ) κ | ≤ 1 2 . (3.15)

On the other hand, let:

φ 0 = 3 π 4 . (3.16)

It can be seen from the above selection that κ and φ 0 here only depend on m, α and b. In this way, for the κ and φ 0 that have been determined above, we can make ε 0 sufficiently small such that when ε < ε 0 , we have:

κ arctan ε + φ 0 = κ arctan ε + 3 π 4 ∈ [ 3 π 4 , 5 π 6 ] . (3.17)

Let’s summarize the properties following the above settings. First, from (3.11) and (3.17), we have:

U | θ = arctan ε = r m + α sin ( κ arctan ε + φ 0 ) ≥ r m + α 2 . (3.18)

Secondly, from (3.14), (3.15) and (3.16), we have:

b U x + U y | θ = 0 = κ cos φ 0 r m + α − 1 ( b ( m + α ) κ tan φ 0 + 1 ) ≤ − r m + α − 1 4 . (3.19)

Finally, from (3.14), (3.16) and (3.17),

Δ U = ( ( m + α ) 2 − κ 2 ) r m + α − 2 sin ( κ θ + φ 0 ) ≤ − r m + α − 2 2 . (3.20)

Here, θ ∈ [ 0 , arctan ε ] . Furthermore,

U | r = 1 = r m + α sin ( κ θ + φ 0 ) | r = 1 = sin ( κ θ + φ 0 ) ≥ 1 2 . (3.21)

Thus, for a positive constant C, we consider C U ± ( u − P m ) . Combining (3.7)-(3.9) and (3.18)-(3.21), by choosing a sufficiently large C, we can get:

C U ± ( u − P m ) > 0     on     l ε ∩ B ;

( b ∂ x + ∂ y ) ( C U ± ( u − P m ) ) < 0     on     l 0 ∩ B ;

C U ± ( u − P m ) > 0     on     Ω ε ∩ ∂ B ;

Δ ( C U ± ( u − P m ) ) < 0     in     Ω ε ∩ B .

According to the extreme value principle, we have:

C U ± ( u − P m ) > 0     in     Ω ε ∩ B ,

that is

| u − P m | ≤ C U ≤ C r m + α     in     Ω ε ∩ B .

The theorem is proved.

This paper studies the regularity of solutions to mixed boundary value problems of second-order linear elliptic equations with small angles in a two-dimensional region. We proved that under the assumption of sufficient regularity of boundary conditions, inhomogeneous terms and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order. Our proof process mainly used the construction of the barrier function, the theory of Schauder estimation, and the regularity conclusion of the solution to the oblique derivative problems in Lipschitz domains. The results can also be extended to high dimensions and high orders. For details, please refer to Lieberman’s work [6] [11] , which has a more systematic explanation.

Sincerely thank Dian Hu for his professional education and careful guidance.

The author declares no conflicts of interest regarding the publication of this paper.

Wu, M.Y. (2024) The Regularity of Solutions to Mixed Bound-ary Value Problems of Second-Order Elliptic Equations with Small Angles. Journal of Applied Mathematics and Physics, 12, 1043-1049. https://doi.org/10.4236/jamp.2024.124064