Next Article in Journal
A Comparison of a Priori Estimates of the Solutions of a Linear Fractional System with Distributed Delays and Application to the Stability Analysis
Next Article in Special Issue
Certain Unified Integrals Involving a Multivariate Mittag–Leffler Function
Previous Article in Journal
Coupled Optimal Results with an Application to Nonlinear Integral Equations
Previous Article in Special Issue
Pseudo-Lucas Functions of Fractional Degree and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Singular Integral Neumann Boundary Conditions for Semilinear Elliptic PDEs

1
Department of Mathematics, Anand International College of Engineering, Jaipur 303012, Rajasthan, India
2
Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman AE 346, United Arab Emirates
3
Centre for Mathematics and Natural Sciences, Leipzig University of Applied Sciences, PF 30 11 66, 04251 Leipzig, Germany
*
Author to whom correspondence should be addressed.
Submission received: 28 February 2021 / Revised: 15 April 2021 / Accepted: 21 April 2021 / Published: 24 April 2021
(This article belongs to the Special Issue Special Functions Associated with Fractional Calculus)

Abstract

:
In this article, we discuss semilinear elliptic partial differential equations with singular integral Neumann boundary conditions. Such boundary value problems occur in applications as mathematical models of nonlocal interaction between interior points and boundary points. Particularly, we are interested in the uniqueness of solutions to such problems. For the sublinear and subcritical case, we calculate, on the one hand, illustrative, rather explicit solutions in the one-dimensional case. On the other hand, we prove in the general case the existence and—via the strong solution of an integro-PDE with a kind of fractional divergence as a lower order term—uniqueness up to a constant.

1. Introduction

The aim of this article is to discuss (possibly singular) semilinear elliptic PDEs of the form
Δ u + c ( · , u ) = 0 in Ω with c ( x , u ) : = Ω b ( x , y , u ) d S ( y )
on bounded C 2 -domains Ω R N subject to (possibly singular) integral Neumann boundary conditions
u n = I u on Ω with ( I u ) ( y ) : = Ω b ( x , y , u ( x ) ) d x
for a (possibly singular) kernel b : Ω × Ω × R R . Such integral boundary value problems occur in applications as mathematical models of nonlocal interaction between interior points and boundary points [1] and have been abstractly studied in the non-singular case by [2]. The linear case of Poisson’s problem and (non-integral) singular Neumann boundary conditions has been discussed in [3]; see [4] for the general linear case and [5] for nonlinear boundary conditions involving a measure. The existence of positive solutions to semilinear singular elliptic problems has been studied in [6] for homogeneous Dirichlet boundary conditions; for (non-integral) homogeneous Neumann boundary conditions, see [7]. For the semilinear problem (1) with (possibly singular) integral Neumann boundary condition (2), existence of very weak solutions u L p ( Ω ) , 1 < p < , with zero average Ω u d x = 0 has been shown in [8] under the assumptions on b that
(A1)
b : Ω × Ω × R R is a Carathéodory function; i.e., b ( x , y , u ) is measurable w.r.t. ( x , y ) Ω × Ω for every u R and continuous w.r.t. u for a.e. ( x , y ) Ω × Ω .
(A2)
b is bounded by | b ( x , y , u ) | b 0 ( x , y ) + b 1 ( x , y ) | u | r with non-negative functions satisfying b 0 ( x , y ) | x y | s + N 1 q L q ( Ω × Ω ) and b 1 ( x , y ) | x y | s + N 1 q L p q p r q ( Ω × Ω ) for some 1 q < ( p ) * , 0 r < max ( p / q , 1 ) and 0 < s < 1 ,
I.e., in the subcritical (w.r.t. x , y ) and sublinear (w.r.t. u) case.
In this article, on the one hand, we discuss one-dimensional examples, which illustrate that without fixing the average Ω u d x the problem (1) and (2), has a one-dimensional continuum of solutions; on the other hand, we show that if the average Ω u d x is fixed, then problems (1) and (2), have a unique very weak solution under the additional assumption.
(A3)
b ( x , y , u ) is monotone w.r.t. u for a.e. ( x , y ) Ω × Ω with derivative 0 b u ( x , y , u ) b 2 ( x , y ) bounded by a function satisfying b 2 ( x , y ) | x y | s + N 1 q L p q p q ( Ω × Ω ) ; particularly, we require p < q < ( p ) * .
For example, b ( x , y , u ) : = b 0 ( x , y ) + b 2 ( x , y ) f ( u ) with f ( u ) : = u | u | 1 1 r | u | r 1 u | u | 1 for 0 r < 1 and functions b 0 , b 2 as in (A2), (A3), satisfies all three conditions. A main difference to previous results of other authors like [1,2] regarding PDEs with integral boundary conditions is that, here, the equation as well as the boundary condition may be singular.

2. Preliminaries

Let Ω R N be a bounded domain. For an exponent p [ 1 , ] , we denote the dual exponent by p : = p p 1 , and the Sobolev conjugate by p * : = N p N p for 1 p < N resp. Consider p * R as arbitrary large for p = N resp. Let p * : = for p > N . Particularly, the dual of the Banach space L p ( Ω ) of the p-integrable function u on Ω can be identified with L p ( Ω ) , the Sobolev embedding W 0 1 , p ( Ω ) L p * ( Ω ) holds for p [ 1 , ] , and the embedding W 0 1 , p ( Ω ) L q ( Ω ) is compact for 1 q < p * due to the Rellich–Kondrachov theorem. Hereby, W 1 , p ( Ω ) denotes the Sobolev space of functions u on Ω having a p-integrable weak gradient u , and functions in W 0 1 , p ( Ω ) additionally satisfy Dirichlet boundary conditions u = 0 on Ω in the sense of traces [9]. Similarly, W 2 , p ( Ω ) denotes the Sobolev space of functions u on Ω having p-integrable second-order derivatives.
Let us exemplify the difference between strong, weak, and very weak solutions of linear elliptic partial differential equations (PDEs) by considering Poisson’s equation Δ u = f for a right hand side (r.h.s.) f. If u W 2 , p ( Ω ) , then Δ u = div ( u ) = i = 1 N 2 u x i 2 is a p-integrable function; thus, for f L p ( Ω ) it makes sense to require Δ u = f almost everywhere (a.e.) in Ω , and in this case, u W 2 , p ( Ω ) is called a strong solution Poisson’s equation. Correspondingly, A : W 2 , p ( Ω ) L p ( Ω ) , A u : = Δ u , is called the strong realization of the negative Laplacian.
However, if the r.h.s. is merely a distribution f ( W 0 1 , p ( Ω ) ) * , then the existence of strong solutions cannot be guaranteed. However, weak solutions u W 0 1 , p ( Ω ) of Poisson’s equation subject to Dirichlet boundary conditions on Ω may exist in the sense that
A u , v : = Ω u · v d x = f , v
holds for every v W 0 1 , p ( Ω ) , and the operator A : W 0 1 , p ( Ω ) ( W 0 1 , p ( Ω ) ) * defined on the left hand side (l.h.s.) is called the weak realization of the negative Dirichlet–Laplacian. Note that the l.h.s. arises from Ω ( Δ u ) v d x via partial integration using v = 0 on Ω . If the r.h.s. is an even worse distribution f ( { v W 2 , p ( Ω ) | v = 0 on Ω } ) * , then another partial integration leads to the notion of a very weak solution.
Definition 1.
A function u L p ( Ω ) is called a very weak solution of Poisson’s equation subject to Dirichlet boundary condition, if
Ω u ( Δ v ) d x = f , v
is valid for every v W 2 , p ( Ω ) with v = 0 on Ω , and A : L p ( Ω ) ( { v W 2 , p ( Ω ) | v = 0 o n Ω } ) * defined by letting A u , v be the l.h.s. of (4) is called the very weak realization of the negative Dirichlet–Laplacian.
Note that during partial integration the term Ω u v n d x arises, but as the derivative of v in direction of the outer normal vector n along the boundary Ω can be arbitrary, for functions f, validity of (4) ensures formally that u satisfies Dirichlet boundary conditions. Particularly, if ( p ) * > N , then { v W 2 , p ( Ω ) | v = 0 on Ω } W 0 1 , ( p ) * ( Ω ) C 0 , α ( Ω ¯ ) holds with α = 1 N ( p ) * by Sobolev embeddings into Hölder spaces. Thus, if b : Ω × Ω R is a non-negative kernel not integrable over Ω × Ω , but satisfying b ( x , y ) | x y | α L 1 ( Ω × Ω ) , then for the possibly non-integrable function f ( x ) : = Ω b ( x , y ) d y on Ω it still makes sense to consider very weak solutions of (4), where the definition f , v : = Ω f ( x ) v ( x ) d x of the r.h.s. makes sense due to
Ω | f ( x ) v ( x ) | d x = Ω Ω b ( x , y ) | x y | α v ( x ) v ( y ) | x y | α d y d x Ω Ω b ( x , y ) | x y | α d y d x v C 0 , α ,
where v ( y ) = 0 for y Ω is used. This case is very similar to Brézis problem, where Poisson’s equation Δ u = f subject to Dirichlet boundary conditions u = 0 on Ω is considered for a measurable function f L 1 ( Ω ) satisfying Ω | f ( x ) | dist ( x , Ω ) d x < + . Brézis et al. [10] (see also [11]) have proved the existence and uniqueness of a very weak solution u L 1 ( Ω ) to this problem, and there is the notion of a very weak solution that prominently occurred for the first time. In this article, we are mainly interested in the very weak solution for the case of singular integral Neumann boundary conditions (instead of homogeneous Dirichlet boundary conditions) and semilinear (instead of linear) elliptic PDEs. To prove the existence and uniqueness, considering the above mentioned facts about Sobolev spaces, we used functional analytic methods like topological degree theory [12].

3. One-Dimensional Examples

In this section, we illustrate by one-dimensional examples—one example is linear but inhomogeneous; another one is sublinear but has vanishing r.h.s—that without fixing the average Ω u d x , the problem (1) and (2), has a one-dimensional continuum of solutions.

3.1. A Linear Example with Singular Right Hand Side

Example 1
(in [8]). In the one-dimensional case of the interval Ω : = ( 0 , 1 ) , where the boundary Ω = { 0 , 1 } consists of two points, and for the kernel b ( x , 0 , u ) : = a ( a 1 ) x s ( u 1 ) containing a parameter a R , b ( x , 1 , u ) : = 0 , which is affine linear w.r.t. u and has a singularity of order s > 0 as x 0 , Equation (1) reads as u ( x ) + a ( a 1 ) x s u ( x ) = a ( a 1 ) x s . In the special case s = 2 , this is an inhomogeneous Cauchy–Euler ODE, and every solution has the form
u ( x ) = C 1 x a + C 2 x 1 a + 1 , C 1 , C 2 R .
Note that u does not admit a classical derivative u ( 0 ) at x = 0 for parameters 1 a < 1 2 except for a = 0 , or a < 0 and C 1 = 0 . The integral Neumann boundary conditions (2) read in this example formally as u ( 0 ) = 0 1 a ( a 1 ) x 2 ( u ( x ) 1 ) d x and u ( 1 ) = 0 . In other words, the outflow through the boundary point 0 is given by the integral 0 1 a ( a 1 ) x 2 ( u ( x ) 1 ) d x ; i.e., there is a nonlocal dependence of the outflow on a singularly weighted average of u over the domain, while there is no flow through the boundary point 1. Note that for every solution u of the form (5), the integral 0 1 a ( a 1 ) x 2 ( u ( x ) 1 ) d x is infinite within the parameter range 1 a < 1 2 , except for a = 0 , or a < 0 and C 1 = 0 . Therefore, the integral Neumann boundary condition at the boundary point 0 reads formally as ± = ± and thus is singular; i.e., there is an explosive inflow or outflow. Nonetheless, Equation (1) together with the singular integral Neumann boundary condition (2) at x = 0 as well as the classical Neumann boundary condition u ( 1 ) = 0 at x = 1 (corresponding to b ( x , 1 , u ) : = 0 ) is satisfied for 1 < a < 1 2 by the integrable functions u ( x ) : = C ( ( 1 a ) x a a x 1 a ) + 1 for arbitrary C R in the very weak sense that
0 1 u ( x ) v ( x ) d x + 0 1 a ( a 1 ) x 2 ( u ( x ) 1 ) ( v ( x ) v ( 0 ) ) d x = 0
holds for every v C 2 ( [ 0 , 1 ] ) with v ( 0 ) = 0 = v ( 1 ) . Thus, the very weak solutions of this singular integral boundary value problem form a one-dimensional affine linear subspace. Finally, fixing the average 0 1 u ( x ) d x determines the constant C and hence the solution uniquely.

3.2. A Semilinear Example with Sublinear Kernel

Example 2.
Again, consider the one-dimensional case Ω : = ( 0 , 1 ) , now for the kernel b ( x , 0 , u ) : = 1 x s u r with signed power u r : = | u | r 1 u , b ( x , 1 , u ) : = 0 , which is sublinear for 0 < r < 1 and has a singularity of order s > 0 as x 0 . Then, Equation (1) reads as u ( x ) + 1 x s ( u ( x ) ) r = 0 and is called an Emden–Fowler or Thomas–Fermi equation. As this equation is a second-order ODE, its general solution contains two constants. The boundary condition u ( 1 ) = 0 relates these constants by one equation. However, as the integral Neumann boundary condition u ( 0 ) = 0 1 1 x s ( u ( x ) ) r d x is simply implied by integrating the ODE u ( x ) = 1 x s ( u ( x ) ) r over Ω, it does not provide another equation for the constants. Hence, there is a one-dimensional continuum of solutions.
It remains to show that these solutions satisfy a singular Neumann BC at x = 0 for u ( 1 ) = 0 and u ( 1 ) > 0 (and similarly for u ( 1 ) < 0 ). Elementarily, u ( x ) > 0 for 0 < x 1 implies due to u ( x ) = 1 x s ( u ( x ) ) r that u is positive near x; i.e., u is strictly convex near x. Hence, u ( x ) is strictly monotone increasing and, thus, due to u ( 1 ) = 0 negative for 0 < x < 1 , so that u ( x ) is strictly monotone decreasing and positive for 0 < x 1 . In the weak formulation 0 1 u v d x + 0 1 1 x s ( u ( x ) ) r ( v ( x ) v ( 0 ) ) d x = 0 , put v ( x ) : = x s 1 to obtain 0 1 1 x ( u ( x ) ) r d x = ( s 1 ) 0 1 1 x 2 s ( u ( x ) ) d x . As u ( x ) u ( 1 ) , we obtain + = u ( 1 ) r 0 1 1 x d x ( s 1 ) 0 1 1 x 2 s ( u ( x ) ) d x , and, hence, | u ( x ) | = u ( x ) tends to + as x 0 due to 2 s < 1 Therefore, for s > 1 and 0 < r < 1 , every solution u 0 of u ( x ) + 1 x s ( u ( x ) ) r = 0 on Ω = ( 0 , 1 ) to the boundary value u ( 1 ) = 0 satisfies a singular Neumann boundary condition at x = 0 .
In the special case of s = 2 , solutions can be more precisely described by u ( x ) : = y ( ln ( x ) ) , where y ( t ) solves the autonomous ODE y ( t ) + y ( t ) + ( y ( t ) ) r = 0 and t is related to x by t : = ln ( x ) 0 . This autonomous ODE can principally be solved by substituting z ( y ) : = y ( t ( y ) ) , which leads to the first-order ODE z ( y ) = 1 + | z ( y ) | r 1 . In the case r = 1 2 , we obtain by separation of variables the equation | z ( y ) | 2 | z ( y ) | r + 2 ln ( 1 + | z ( y ) | r ) = y C with a constant C R . As z ( y ( t ) ) = y ( t ) holds, y ( t ) satisfies the implicit ODE | y ( t ) | 2 | y ( t ) | r + 2 ln ( 1 + | y ( t ) | r ) = y ( t ) C . A phase plot of y is shown in the left of Figure 1 for C = 0 , and a numerical solution u with u ( 1 ) = 0 and u ( 1 ) = 1 showing singularity at x = 0 is plotted to the right in Figure 1.

4. Existence

In this section, for dimensions N > 1 , let us briefly extend the proof of existence of very weak solutions from [8] to functions with an arbitrary pregiven average.
Theorem 1.
Let Ω R N be a bounded C 2 -domain and let b : Ω × Ω × R R be a kernel satisfying (A1) and (A2) for some 1 < p < . Then, for every constant a R , there exists u L p ( Ω ) , satisfying Ω u ( x ) d x = a , which solves (1) and (2) in the very weak sense that
Ω u ( Δ v ) d x + Ω Ω b ( x , y , u ( x ) ) ( v ( x ) v ( y ) ) d S ( y ) d x = 0
is valid for every v W 2 , p ( Ω ) with v n = 0 on Ω .
In the following, denote by
  • D the subspace
    D : = { v W 2 , p ( Ω ) | v n = 0 on Ω } W 2 , p ( Ω )
    which has a compact embedding D W 1 , q ( Ω ) for q < ( p ) * .
  • W s , q ( Ω ; Ω ) , the mixed fractional Sobolev space of order s and exponent q, which can be considered a space of functions v on Ω ¯ such that v | Ω is q-integrable over Ω and v ( x ) v ( y ) | x y | s + ( N 1 ) / q is q-integrable over Ω × Ω (particularly, v | Ω L q ( Ω ) ), where functions on Ω ¯ are identified if they coincide a.e. on Ω and a.e. on Ω .
  • A : L p ( Ω ) D * ,
    A u , v : = Ω u ( Δ v ) d x ,
    the very weak realization of the negative Neumann–Laplacian Δ .
  • B : L p ( Ω ) ( W s , q ( Ω ; Ω ) ) * ,
    B u , v : = Ω Ω b ( x , y , u ( x ) ) ( v ( x ) v ( y ) ) d S ( y ) d x
    the realization of the nonlinearity and the integral Neumann boundary conditions, which may be singular.
Proof of Theorem 1.
Restrict A and B to the closed affine linear subspace U : = { u L p ( Ω ) | Ω u d x = a } of L p ( Ω ) (of codimension 1) so that A becomes injective. By (A1) and (A2), B : L p ( Ω ) ( W s , q ( Ω ; Ω ) ) * ( W 1 , q ( Ω ) ) * is a bounded continuous mapping which becomes compact when viewed as mapping into D * due to the compactness of D W 1 , q ( Ω ) . Thus, we can view A + B as a perturbation of A by a compact operator B. To conclude, regarding the existence of very weak solutions u L p ( Ω ) with Ω u d x = a , we may apply topological degree theory as in [8], but for this, we need to exclude the existence of solutions with an arbitrarily large L p -norm:
For u L p ( Ω ) with Ω u d x = a let v D be a strong solution of
Δ v = | u | p 2 u λ
subject to homogeneous Neumann boundary conditions v n = 0 on Ω , where the constant λ : = 1 | Ω | Ω | u | p 2 u d x is chosen such that the right hand side of (9) has zero integral. Thus, the compatibility condition is satisfied, and there exists a strong solution v of (9) subject to homogeneous Neumann boundary conditions, which is unique up to an additive constant. Moreover, such a solution v satisfies
[ v ] s , q , Ω , Ω C v q , Ω C | u | p 2 u λ p , Ω C ( u p , Ω p 1 + λ p , Ω ) C u p , Ω p 1
with a constant C < independent of u, as λ p , Ω = | λ | | Ω | 1 / p | Ω | 1 / p | Ω | 1 / p u p , Ω p 1 = u p , Ω p 1 holds by the definition of λ . Further, by the definition of v we have
A u , v = Ω | u | p d x λ Ω u d x u p , Ω p a | Ω | 1 / p u p , Ω p 1
due to λ 1 | Ω | Ω | u | p 1 d x and Jensen’s inequality 1 | Ω | Ω | u | p 1 d x p 1 | Ω | Ω | u | p d x . From (A2) we can conclude
| B u , v | ( C 1 + C 2 u p , Ω r ) [ v ] s , q , Ω , Ω
with 0 r < max ( p / q , 1 ) , and together with (10) both inequalities imply
( A + B ) u , v u p , Ω p a | Ω | 1 / p u p , Ω p 1 ( C 1 + C 2 u p , Ω r ) u p , Ω p 1 .
As the right hand side tends to infinity for u p , Ω due to r < 1 , the equation ( A + B ) u = 0 has no solution u L p ( Ω ) satisfying Ω u d x = a on the boundary of a sufficiently large ball around zero in L p ( Ω ) . Therefore, like A u = 0 , ( A + B ) u = 0 also has a solution u L p ( Ω ) with Ω u d x = a . □
Remark 1.
Existence even holds in the linear case r = 1 , p < q < ( p ) * , provided that C 2 < 1 or equivalently b 1 ( x , y ) | x y | s + N 1 q p q p q , Ω × Ω is sufficiently small.

5. Duality

To prepare the proof of uniqueness of very weak solutions to (1) and (2), let us consider in this section the affine linear case b ( x , y , u ) = b ˜ ( x , y ) u 1 | Ω | f ˜ ( x ) , where the sublinearity required in (A2) as well as the smallness condition mentioned in Remark 1 are not satisfied, but instead b ˜ satisfies the same condition as b 2 in (A3); i.e., b ˜ 0 and b ˜ ( x , y ) | x y | s + N 1 q L p q p q ( Ω × Ω ) and f ˜ L p ( Ω ) have average Ω f ˜ d x = 0 . Then problems (1) and (2) become linear and read (with w instead of u) as
Δ w + c ˜ ( · ) w = f ˜ in Ω with c ˜ ( x ) : = Ω b ˜ ( x , y ) d S ( y ) , w n ( y ) = Ω b ˜ ( x , y ) w ( x ) d x for y Ω .
The dual to this linear problem is
Δ v + B ˜ v = f λ in Ω , v n = 0 on Ω ,
where f L p ( Ω ) , λ R and
( B ˜ v ) ( x ) : = Ω b ˜ ( x , y ) ( v ( x ) v ( y ) ) d S ( y )
is a kind of fractional divergence of the function v with kernel b ˜ ( x , y ) .
Example 3.
In the one-dimensional nonsingular case b ˜ 1 2 , f 0 , (12) with (13) reads as v + v 1 2 ( v ( 0 ) + v ( 1 ) ) = 0 or equivalently as v + v + 0 1 ( 1 2 t ) v ( t ) v ( t ) d t = 0 , which is an integro-ODE. Differentiating v + v 1 2 ( v ( 0 ) + v ( 1 ) ) = 0 gives v + v = 0 with fundamental system v ( x ) = C 1 cosh ( x ) + C 2 sinh ( x ) + C 3 , but the constant functions v ( x ) = C 3 simply satisfy the homogeneous Neumann boundary conditions v ( 0 ) = 0 = v ( 1 ) ; i.e., we have uniqueness up to constants.
Lemma 1.
For b ˜ 0 satisfying b ˜ ( x , y ) | x y | s + N 1 q L p q p q ( Ω × Ω ) with p < q < ( p ) * and every f L p ( Ω ) , there exists a λ R such that (12) and (13) have, with the addition of constants, a unique strong solution v W 2 , p ( Ω ) .
Proof. 
With the strong realization A s : D L p , A s v : = Δ v , of the Neumann–Laplacian and the compact operator B ˜ : D W 1 , q ( Ω ) W s , q ( Ω ; Ω ) L p ( Ω ) , we seek a solution v D of ( A s + B ˜ ) v = f λ , i.e., a strong solution satisfying homogeneous Neumann boundary conditions. A problem closely related to (12) is
Δ v ˜ + B ˜ v ˜ = f λ ˜ in Ω , v ˜ n = Ω b ˜ ( x , y ) ( v ˜ ( x ) v ˜ ( y ) ) d x on Ω
with λ ˜ : = 1 | Ω | Ω f d x . In fact, note that the weak formulation of this problem is the identification of a v ˜ W 1 , 2 ( Ω ) such that
Ω v ˜ · w d x + Ω Ω b ˜ ( x , y ) ( v ˜ ( x ) v ˜ ( y ) ) ( w ( x ) w ( y ) ) d S ( y ) d x = Ω ( f λ ˜ ) w d x
for every w W 1 , 2 ( Ω ) . For b ˜ ( x , y ) 0 with b ˜ ( x , y ) | x y | 2 s + ( N 1 ) L ( Ω × Ω ) , the existence and uniqueness of the solution v ˜ W 1 , 2 ( Ω ) with Ω v ˜ d x = 0 follows Lax–Milgram and Poincáre inequality. In fact, the second term in (15) has the potential Φ ( v ) : = 1 2 Ω Ω b ˜ ( x , y ) | v ( x ) v ( y ) | 2 d S ( y ) d x , and due to b ˜ 0 , the whole operator on the l.h.s. of the weak formulation has a coercive and weakly lower semicontinuous potential on W 1 , 2 ( Ω ) / R .
Further, solving Poisson’s equation Δ v ¯ = λ ¯ subject to inhomogeneous Neumann BC v ¯ n = Ω b ˜ ( x , y ) ( v ˜ ( x ) v ˜ ( y ) ) d x on Ω , where λ ¯ : = Ω Ω b ˜ ( x , y ) ( v ˜ ( x ) v ˜ ( y ) ) d S ( y ) d x is such that a (strong) solution v ¯ exists. Then v : = v ˜ v ¯ satisfies Δ v + B ˜ v = f λ in Ω with λ : = λ ˜ + λ ¯ and v n = 0 on Ω , i.e., (12), in a weak sense. Finally, due to B ˜ v L p ( Ω ) and f L p ( Ω ) , the weak solution v is in fact a strong solution, and, by approximation, this also holds true under the regularity assumptions on b ˜ 0 in Lemma 1. □

6. Uniqueness

To prove the uniqueness of very weak solutions to (1) and (2), we additionally require b to satisfy (A3); i.e. the derivative b u is bounded by a function b 2 ( x , y ) independent of u satisfying b 2 ( x , y ) | x y | s + N 1 q L p q p q ( Ω × Ω ) . As a consequence, the difference quotient ( x , y , u , u ˜ ) b ( x , y , u ) b ( x , y , u ˜ ) u u ˜ is continuous w.r.t. ( u , u ˜ ) for a.e. ( x , y ) Ω × Ω and thus a Carathéodory function.
Theorem 2.
Let Ω R N be a bounded C 2 -domain and let b : Ω × Ω × R R be a kernel satisfying (A1)–(A3) for some 1 < p < . Then, for every constant a R , the very weak solution u L p ( Ω ) of (1) and (2) satisfying Ω u ( x ) d x = a is unique.
Proof. 
Let u , u ˜ L p ( Ω ) be two very weak solutions of (1) and (2), with identical averages Ω u d x = Ω u ˜ d x . Let v W 2 , p ( Ω ) be a strong solution of
Δ v + B ˜ v = | u u ˜ | p 2 ( u u ˜ ) λ in Ω , v n = 0 on Ω ,
where
B ˜ v ( x ) : = Ω b ˜ ( x , y ) ( v ( x ) v ( y ) ) d S ( y )
is a kind of fractional divergence of the function v with kernel b ˜ ( x , y ) : = b ( x , y , u ( x ) ) b ( x , y , u ˜ ( x ) ) u ( x ) u ˜ ( x ) , and λ R is chosen such that there exists a solution v satisfying homogeneous Neumann boundary conditions v n = 0 on Ω . Note that Lemma 1 guarantees the existence of such a strong solution. Now, test the difference between the equations solved by u and u ˜ with the test function v to obtain
Ω | u ( x ) u ˜ ( x ) | p d x Ω Ω b ˜ ( x , y ) ( u ( x ) u ˜ ( x ) ) ( v ( x ) v ( y ) ) d S ( y ) d x + Ω Ω ( b ( x , y , u ( x ) ) b ( x , y , u ˜ ( x ) ) ) ( v ( x ) v ( y ) ) d S ( y ) d x = 0
and due to b ˜ ( x , y ) ( u ( x ) u ˜ ( x ) ) = b ( x , y , u ( x ) ) b ( x , y , u ˜ ( x ) ) , Ω | u u ˜ | p d x = 0 . Hence, u = u ˜ holds. □

7. Conclusions

We were able to prove uniqueness up to a constant for very weak solutions to singular semilinear elliptic PDEs subject to singular integral Neumann boundary conditions. Our method to test for two very weak solutions of (1) and (2), with identical averages of the difference in the equations determined by the strong solution of a dual problem (which is an integro-PDE with a kind of fractional divergence as a lower-order term here), seems to be rather promising and may be applicable to show the uniqueness of very weak solutions for many other problems.
However, our results are still not optimal. For example, while we used an intermediate compact embedding W 2 , p ( Ω ) W 1 , ( p ) * ( Ω ) , it seems to be an open question for which ( s , q ) the direct embedding W 2 , p ( Ω ) W s , q ( Ω ; Ω ) is compact. In fact, the embedding and interpolation properties of mixed fractional Sobolev spaces have not been studied in depth in the literature to date. Therefore, although there are some first steps, as in the Appendix B of [13], for example, how singular a domain-boundary kernel b can be without destroying the existence of solutions is an open question.

Author Contributions

Conceptualization, J.M. and P.A.; methodology, J.M.; software, G.S.; validation, J.M., P.A. and G.S.; formal analysis, J.M.; investigation, J.M.; resources, P.A.; data writing—original draft preparation, J.M.; writing—review and editing, J.M., P.A. and G.S.; visualization, G.S.; project administration, P.A.; funding acquisition, P.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German DAAD/BMBF and the Indian DST for financial support within the DAAD-DST project PPP 2019 “Mathematical models of non-local interactions” grant number 57453186.

Acknowledgments

Jochen Merker is grateful to Jean-Michel Rakotoson, who suggested during the conference in Toulouse on the occasion of Peter Takáč’s birthday 2019 to prove uniqueness up to constants for solutions to semilinear elliptic PDEs with singular integral Neumann boundary conditions, and shared an idea of how to prove uniqueness under a kind of smallness condition.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Avalishvili, G.; Avalishvili, M.; Gordeziani, D. On a nonlocal problem with integral boundary conditions for a multidimensional elliptic equation. Appl. Math. Lett. 2011, 24, 566–571. [Google Scholar] [CrossRef] [Green Version]
  2. Shakhmurov, V. Abstract elliptic equations with integral boundary conditons. Chin. Ann. Math. B 2016, 37, 625–642. [Google Scholar] [CrossRef]
  3. Merker, J.; Rakotoson, J.-M. Very weak solutions of Poisson’s equation with singular data under Neumann boundary conditions. Calc. Var. Partial. Differ. Equ. 2015, 52, 705–726. [Google Scholar] [CrossRef]
  4. Merker, J. On a generalization of Sobolev-Slobodecki spaces and an associated fractional derivative. In Proceedings of the International Conference on Fractional Differentiation and Its Applications (ICFDA), Amman, Jordan, 16–18 July 2018. [Google Scholar]
  5. Oussama Boukarabila, Y.; Véron, L. Nonlinear boundary value problems relative to harmonic functions. Nonl. Anal. 2015, 201, 112090. [Google Scholar] [CrossRef]
  6. Lair, A.; Shaker, A.W. Classical and Weak Solutions of a Singular Semilinear Elliptic Problem. J. Math. Anal. Appl. 1997, 211, 371–385. [Google Scholar] [CrossRef] [Green Version]
  7. Chabrowski, J. On the Neumann problem with singular and superlinear nonlinearities. Commun. Appl. Anal. 2009, 13, 327–340. [Google Scholar]
  8. Merker, J. Mixed fractional Sobolev spaces and elliptic PDEs with singular integral boundary data. Math. Meth. Appl. Sci. 2021, 44, 2859–2867. [Google Scholar] [CrossRef] [Green Version]
  9. Adams, R.A.; Fournier, J.J.F. Sobolev Spaces; Elsevier: Amsterdam, The Netherlands, 2003. [Google Scholar]
  10. Brézis, H.; Cazenave, T.; Martel, Y.; Ramiandrisoa, A. Blow-up for ut − Δu = g(u) revisited. Adv. Differ. Equ. 1996, 1, 73–90. [Google Scholar]
  11. Brézis, H.; Cabré, X. Some simple nonlinear PDE’s without solutions. Boll. Unione Mat. Ital. 1998, 1, 223–262. [Google Scholar]
  12. Cho, Y.; Chen, Y.Q. Topological Degree Theory and Applications; Chapman and Hall/CRC: New York, NY, USA, 2006. [Google Scholar]
  13. Merker, J. Very weak solutions of linear elliptic PDEs with singular data and irregular coefficients. Differ. Equ. Appl. 2018, 10, 3–20. [Google Scholar] [CrossRef]
Figure 1. Left: phase portrait; right: solution curve.
Figure 1. Left: phase portrait; right: solution curve.
Axioms 10 00074 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Agarwal, P.; Merker, J.; Schuldt, G. Singular Integral Neumann Boundary Conditions for Semilinear Elliptic PDEs. Axioms 2021, 10, 74. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10020074

AMA Style

Agarwal P, Merker J, Schuldt G. Singular Integral Neumann Boundary Conditions for Semilinear Elliptic PDEs. Axioms. 2021; 10(2):74. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10020074

Chicago/Turabian Style

Agarwal, Praveen, Jochen Merker, and Gregor Schuldt. 2021. "Singular Integral Neumann Boundary Conditions for Semilinear Elliptic PDEs" Axioms 10, no. 2: 74. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10020074

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop