In this paper, we prove the existence of solutions for an elliptic system. More precisely, we combine the potential theory with the sub-super solution method and use the properties of the well-known Kato class to justify our existence results. The novelty of our study is that we consider either the bounded or the exterior domain; Also, the nonlinearities may be singular near the boundary. Some examples are presented to validate our main results.
Citation: Abdeljabbar Ghanmi, Hadeel Z. Alzumi, Noureddine Zeddini. A sub-super solution method to continuous weak solutions for a semilinear elliptic boundary value problems on bounded and unbounded domains[J]. Electronic Research Archive, 2024, 32(6): 3742-3757. doi: 10.3934/era.2024170
[1] | S. Bandyopadhyay, M. Chhetri, B. B. Delgado, N. Mavinga, R. Pardo . Maximal and minimal weak solutions for elliptic problems with nonlinearity on the boundary. Electronic Research Archive, 2022, 30(6): 2121-2137. doi: 10.3934/era.2022107 |
[2] | Chenghua Gao, Enming Yang, Huijuan Li . Solutions to a discrete resonance problem with eigenparameter-dependent boundary conditions. Electronic Research Archive, 2024, 32(3): 1692-1707. doi: 10.3934/era.2024077 |
[3] | Yijun Chen, Yaning Xie . A kernel-free boundary integral method for reaction-diffusion equations. Electronic Research Archive, 2025, 33(2): 556-581. doi: 10.3934/era.2025026 |
[4] | Limin Guo, Weihua Wang, Cheng Li, Jingbo Zhao, Dandan Min . Existence results for a class of nonlinear singular $ p $-Laplacian Hadamard fractional differential equations. Electronic Research Archive, 2024, 32(2): 928-944. doi: 10.3934/era.2024045 |
[5] | Youngmok Jeon, Dongwook Shin . Immersed hybrid difference methods for elliptic boundary value problems by artificial interface conditions. Electronic Research Archive, 2021, 29(5): 3361-3382. doi: 10.3934/era.2021043 |
[6] | Meng Wang, Naiwei Liu . Qualitative analysis and traveling wave solutions of a predator-prey model with time delay and stage structure. Electronic Research Archive, 2024, 32(4): 2665-2698. doi: 10.3934/era.2024121 |
[7] | Qingcong Song, Xinan Hao . Positive solutions for fractional iterative functional differential equation with a convection term. Electronic Research Archive, 2023, 31(4): 1863-1875. doi: 10.3934/era.2023096 |
[8] | Yong Zhou, Jia Wei He, Ahmed Alsaedi, Bashir Ahmad . The well-posedness for semilinear time fractional wave equations on $ \mathbb R^N $. Electronic Research Archive, 2022, 30(8): 2981-3003. doi: 10.3934/era.2022151 |
[9] | Jon Johnsen . Well-posed final value problems and Duhamel's formula for coercive Lax–Milgram operators. Electronic Research Archive, 2019, 27(0): 20-36. doi: 10.3934/era.2019008 |
[10] | Massimo Grossi . On the number of critical points of solutions of semilinear elliptic equations. Electronic Research Archive, 2021, 29(6): 4215-4228. doi: 10.3934/era.2021080 |
In this paper, we prove the existence of solutions for an elliptic system. More precisely, we combine the potential theory with the sub-super solution method and use the properties of the well-known Kato class to justify our existence results. The novelty of our study is that we consider either the bounded or the exterior domain; Also, the nonlinearities may be singular near the boundary. Some examples are presented to validate our main results.
In recent years, semilinear elliptic boundary value problems have attracted more attention. This is due to their importance in several fields such as chemical reactions, population evolution, and pattern formation, see, for instance, [1] for other related applications. Due to their importance, several researchers have concentrated on the development of problems involving semilinear elliptic operators. In this paper, we will continue in this direction, so we fix a C1,1 domain Λ in RN, (N≥3), with a nonempty compact boundary ∂Λ, and we consider the following problem:
{−Δu=ψ(x,u) in Λ,u=f, on ∂Λ, | (1.1) |
and in the case when Λ is unbounded, we assume the following supplementary condition
lim|x|→∞u(x)=c. | (1.2) |
Also, we consider the following semilinear elliptic system
{−Δu=ψ1(x,u,v) in Λ,−Δv=ψ2(x,u,v) in Λ,u=f1,v=f2, on ∂Λ, | (1.3) |
and in the case when Λ is unbounded, we assume the following supplementary condition
lim|x|→∞u(x)=c1 and lim|x|→∞v(x)=c2, | (1.4) |
where c, c1, and c2 are real numbers, and the functions f, f1, and f2 are continuous on ∂Λ. Problems like (1.1) or like (1.3) are extensively studied by several authors and by different methods. For interested readers, we refer to the works of, Akô [2] (Schauder's estimates in the Banach space of Hölder continuous functions), Alsaedi et al. [3] (combination of the Karamata regular variation theory with a related comparison principle), Amann [4] (fixed point index), Clément and Sweers [5] (monotone iterative and sub-super solution methods), Cui [6] (sub-super solution method and Sobolev-Morrey's inequality), Keller [7] (strong Maximum principle), Montenegro and Ponce [8] (method of sub-supersolutions combined with Schauder's fixed point theorem), Montenegro and Suárez [9] (adequate sub-super solution method for singular systems), [10] (iterative method combined with Schauder's type and Sobolev inequalities), Noussair and Swanson [11] (Atkinson's theorem), Ogata [12] (monotone arguments combined with iterative methods), Rǎdulecu and Repovš [13] (monotone argument combined with Variational method). Several works in the literature treat problems like (1.1) or problems like (1.3) in the case where the nonlinearities ψ, ψ1, and ψ2 are continuous concerning x (in the case of classical solutions) or Caratheodory functions (in the case of weak solutions in the Sobolev spaces); Moreover, these problems are generally considered in regular bounded domains. Our goal in this paper is to ensure the validity of the sub-super solution method for continuous distributional solutions in the cases where f, ψ1, and ψ2 may be singular concerning x near the boundary. This will be done using a fixed-point argument based on the compacity property of a class of potential functions defined in Λ. To be more precise, these singularities are related to the Kato class K(Λ) which is introduced and studied by Bachar et al. [14] for the exterior domain, and by Mâagli and Zribi [15] for the bounded domain. This class is used to find solutions for several semilinear elliptic problems. We cite, for example, the papers of Alsaedi et al. [16], Bachar et al. [14], Ghanmi et al. [17], Mâagli and Zribi [15], and Zeddini and Sari [18]. Our results for (1.1) apply to prove the existence of continuous solutions for singular nonlinearities such that |ψ(x,u)|≤1(δ(x)λ)|h(u)| on a C1,1-bounded domain with δ(x)=dist(x,∂Λ) and λ<2, and in the case where Λ is a C1,1-exterior domain, our results for (1.1) apply to nonlinearities satisfying |ψ(x,u)|≤1(1+|x|)μ−λ(δ(x)λ)|h(u)| with λ<2<μ. Our results concerning sub-super-solution methods involving nonlinearities that may be singular near the boundary are new and have not been discussed before.
In this paper, we continue to study such problems using the Kato class. The novelty of our study is that we consider either the bounded or the exterior domains; moreover, the nonlinearities used in our problems can be singular, which means more complicated manipulation of our study. More precisely, we transform our problem to an equivalent integral equation, and after that, we define an associated operator, which is (using the properties of the Kato class) relatively compact. Finally, we prove that the fixed points of the associated operator are weak solutions for the studied problem.
Before giving the main results of this paper, we assume the following hypothesis:
(H1) ψ∈B(Λ×R), such that for almost every x∈Λ, the function t→ψ(x,t) is continuous. Moreover, for all C>0, there exists a function pC∈K(Λ) such that for any t∈ and x∈Λ, we have
|ψ(x,t)|≤pC(x),∀(x,t)∈Λ×[−C,C], |
where B(D×R) and K(Λ) are introduced in Section 2.
(H2) For i=1,2, the function ψi∈B(Λ×R and the map (s,t)→ψi(x,s,t) are continuous on R×R for almost every x∈Λ.
(H3) For every C>0 there exists a nonnegative function pi,C(i=1,2), such that
|ψi(x,s,t)|≤pi,C(x),∀(x,s,t)∈Λ×[−C,C]×[−C,C]. |
Our main results from this work are the following theorems:
Theorem 1.1. Assume that the function f is continuous on ∂Λ and the hypothesis (H1) is satisfied. If (1.1) and (1.2) have a continuous sub solution u_ and a continuous super solution ¯u with u_≤¯u in Λ, then problems (1.1) and (1.2) admit a continuous weak solution u∈C(Λ∪∂∞Λ) satisfying
u_≤u≤¯u,in ¯Λ. |
Theorem 1.2. Assume that the functions fi(i=1,2) are continuous on ∂Λ and the hypotheses (H2)and(H3) are satisfied. If (1.3) and (1.4) have a double pair of continuous sub-super solution (u_,v_) and (¯u,¯v), then, problems (1.3) and (1.4) admit a continuous weak solution (u,v)∈C(Λ∪∂∞Λ)×C(Λ∪∂∞Λ) satisfying in addition
u_≤u≤¯u,andv_≤v≤¯v,in ¯Λ, |
where ∂∞Λ and the notion of sub-super solution are introduced in Section 2.
Next, in Section 2, we introduce some notations and we present several properties of the Kato class K(Λ). Section 3 is devoted to the proofs of our main results. As applications of our main results, four examples are presented in Section 4 to validate the above theorems.
In this section, we begin by giving some notations, which will be used later in Section 3. We denote by B(Λ) the set of all Borel measurable functions in Λ, by B+(Λ) the subset of all nonnegative functions of B(Λ), and by Bb(Λ) the subset of the bounded ones. Also, we denote by C0(Λ) the set of continuous functions in Λ that tend to zero near ∂Λ and satisfy, in addition, lim|x|→∞u(x)=0 in the case of unbounded domain. Now, let us denote by C(¯Λ) the subset of B(Λ) composed by continuous functions in ¯Λ in the case when Λ is bounded and in the case when Λ is unbounded, C(¯Λ∪{∞}) will denotes the subset of B(Λ) composed of continuous functions in ¯Λ for which the limit as |x|→∞ exists and is finite.
∂∞Λ will denotes ∂Λ in the case when Λ is bounded and ∂Λ∪{∞} in the case when Λ is unbounded. Consequently, C(Λ∪∂∞Λ) will denote C(¯Λ∪{∞}) if Λ is unbounded and C(¯Λ) if Λ is bounded. Also, we denote by D(Λ) the set of all C∞-functions in Λ with compact support in Λ. Now, we recall that the supremum norm is defined for u∈C(Λ∪∂∞D) by
‖u‖∞=supx∈¯Λ|u(x)|, |
furthermore it is well known that the normed space (C0(Λ),‖.‖∞) is also a Banach space. The Green function of the Dirichlet Laplacian in Λ will be denoted by GΛ: moreover, for a given function p in B+(Λ), the Green potential Vp of a function p is defined on Λ as follows:
Vp(x)=∫ΛGΛ(x,y)p(y)dy. |
It is well known (see [19] p.52) that if p∈L1loc(Λ) is such that Vp∈L1loc(Λ), then (in the sense of distributions) we have
Δ(Vp)=−p in Λ. | (2.1) |
Let f be a nonnegative continuous function on ∂Λ, then, HΛf will denotes the unique solution in C2(Λ)∩C(¯Λ) of the following problem
{Δu=0 in Λu=f on ∂Λ, |
and in the case of an unbounded domain, the above problem is subject to the following condition:
lim|x|→∞u(x)=0. |
Hereafter, ξ will denote the following function
ξ=1−HΛ1. |
It is not difficult to see that the function ξ is harmonic and equal to zero at the boundary of Λ; moreover, if Λ is unbounded, then lim|x|→∞ξ(x)=1.
Since we use the potential theory, it is natural to define the Kato class, which is defined in the following definition.
Definition 2.1. (See [14,15].) A function p∈B(Λ) is said to be in the Kato class K(Λ) if we have
limσ→0supx∈Λ∫Λ∩D(x,σ)ρ(y)ρ(x)GΛ(x,y)|p(y)|dy=0, |
in the case when Λ is bounded, and in addition, in the case of an unbounded domain, we have
limC→∞supx∈Λ∫Λ∩{|y|≥C}ρ(y)ρ(x)GΛ(x,y)|p(y)|dy=0, |
where D(x,α) is the open ball with center x and radius α, δ(x)=d(x,∂Λ), and ρ(x)=min(1,δ(x)).
We note that if Λ is bounded, then we will use the following elementary inequality:
11+dδ(x)≤ρ(x)≤δ(x), |
where d is the diameter of Λ. This means that we can replace ρ by δ in Definition 2.1. Moreover, Zeddini and Sari [18] proved that in the case when Λ is a C1,1-bounded domain, then, this definition is equivalent to
limσ→0(sup(x,y)∈Λ×Λ∫Λ∩(D(x,σ)∪D(y,σ))GΛ(x,z)GΛ(z,y)GΛ(x,y)|p(z)|dz)=0. |
The Kato class K(Λ) is very important in the manipulation of the Green potential, it is quite rich and it contains several functions, as shown in the following example.
1) If Λ is bounded, then the function x↦1(δ(x))λ is in K(Λ) if and only if λ∈(−∞,2).
2) In the case when Λ is the open unit ball, then a radial function p is in K(Λ) if and only if ∫10r(1−r)|p(r)|dr<∞.
3) If Λ is an exterior domain. Then the function x→1(|x|+1)μ−λδ(x)λ is in K(Λ) if and only if λ<2<μ.
4) If Λ is the exterior of the unit closed ball, a radial function p is in K(Λ) if and only if ∫∞1(r−1)|p(r)|dr<∞.
Remark 2.2. If Λ is a C1,1 exterior domain and a function q is nontrivial and nonnegative in K(Λ). Then its Green potential, Vq, is positive in Λ. Indeed, q∈L1loc(Λ), moreover, there exists a compact subset F of Λ such that
0<∫Fq(y)dy<∞. |
Without loss of generality, we can assume that 0∉Λ. Then it has been proved in Bachar et al. [14], that there exists C>0 such that
Cδ(x)|x|N−1δ(y)|y|N−1≤GΛ(x,y),∀(x,y)∈Λ2. |
Hence, for every x∈Λ, we have
Vq(x)=∫ΛGΛ(x,y)q(y)dy≥Cδ(x)|x|N−1∫Fδ(y)|y|N−1q(y)dy≥Cδ(x)|x|N−1infz∈F(δ(z)|z|N−1)∫Fq(y)dy>0. |
Next, to prove the existence of solutions, we use the sub-super solution method, so in the following, we define such a notion.
Definition 2.3. A function u_∈C(Λ∪∂∞Λ) is said to be a continuous sub-solution of the problems (1.1) and (1.2) if the following statements are true:
(i) limx⟶ξ∈∂Λu_(x)≤f(ξ), and in the case of an unbounded domain, it satisfies, in addition, lim|x|⟶∞u_(x)≤c.
(ii) For any nonnegative function φ∈D(Λ), we get
∫Λu_(x)Δφ(x)+ψ(x,u_(x))φ(x)dx≥0. |
The definition of a super-solution of the problems (1.1) and (1.2) is obtained similarly by reversing the inequality in the last definition.
An analog definition is adopted for the problems (1.3) and (1.4). This type of definition is utilized by Gfaifia et al. [20], and Pao [21]. Next, we recall this general definition.
Definition 2.4. ((u_,v_), (¯u,¯v))∈(C(D∪∂∞Λ))2×(C(D∪∂∞Λ))2 is said to be a sub-super solution of problems (1.3) and (1.4) if the following statements hold:
(i) For all x∈Λ, we have
u_(x)≤¯u(x),and v_(x)≤¯v(x). |
(ii) We have
{limx⟶ξ∈∂Λu_(x)≤f1(ξ)≤limx⟶ξ∈∂Λ¯u(x)limx⟶ξ∈∂Λv_(x)≤f2(ξ)≤limx⟶ξ∈∂Λ¯v(x), |
and if in addition Λ is unbounded, then we get
{lim|x|⟶∞u_(x)≤c1≤lim|x|⟶∞¯u(x),lim|x|⟶∞v_(x)≤c2≤lim|x|⟶∞¯v(x). |
(iii) For any nonnegative function φ∈D(Λ) and any (u,v)∈[u_,¯u]×[v_,¯v], we have
{∫Du_(x)Δφ(x)+ψ1(x,u_(x),v(x))φ(x)dx≥0,∫Dv_(x)Δφ(x)+ψ2(x,u(x),v_(x))φ(x)dx≥0, |
moreover, the last inequalities hold by reversing the inequality and replacing u_ and v_ by ¯u and ¯v, respectively.
In this section, we present the proofs of our main results (Theorems 1.1 and 1.2). Firstly, let us introduce the following key result, which can be found in Mâagli and Zribi [15] in the case when Λ is bounded, and in Bachar et al. [14] in the case when Λ is unbounded.
Proposition 3.1. For a given function p in K(Λ). The following statements hold:
1) Vp is a continuous functions in Λ and tends to zero on ∂Λ.
2) The family of functions
{Vq such that |q|≤|p|}, |
is equicontinuous in Λ∪∂∞Λ.
3) p∈L1loc(Λ).
We note that from the well-known Ascoli's theorem and (2) in the previous proposition, we deduce that {Vq such that |q|≤|p|} is relatively compact in C0(Λ).
In order to prove Theorem 1.1, we begin by defining the auxiliary function ˜ψ on Λ×R by
˜ψ(x,u(x))={ψ(x,u_(x)) if u(x)<u_(x),ψ(x,u(x)) if u_(x)≤u(x)≤¯u(x),ψ(x,¯u(x)) if ¯u(x)<u(x). | (3.1) |
Now, let us consider the following associated problem
{−Δu=˜ψ(x,u) in Λ,u=f, on ∂Λlim|x|→∞u(x)=c. | (3.2) |
We begin by proving that problem (3.2) admits a weak solution u in C(Λ∪∂∞Λ). Since the function t↦ψ(x,t) is continuous, the function t↦˜ψ(x,t) is also continuous. On the other hand, if we put C=‖u_‖∞+‖¯u‖∞, then from hypothesis (H1) there exists a nonnegative function pC∈K(Λ) such that
|ψ(x,t)|≤pC(x),∀(x,t)∈Λ×[−C,C]. |
Therefore, we deduce that
|˜ψ(x,u(x))|≤pC(x),∀(x,u)∈Λ×C(Λ∪∂∞Λ). |
So, Proposition 3.1 implies that
˜ψ(.,cξ(.)+HΛϕ(.)+v(.))∈K(Λ), |
moreover, for each v∈C0(Λ) the family
{V(˜ψ(.,cξ+HΛϕ+v)):v∈C0(Λ)}, |
is relatively compact in (C0(Λ),‖.‖∞).
Now we define the operator T:C0(Λ)→C0(Λ) by
Tv(x)=V(˜ψ(.,cξ+HΛϕ+v))(x). |
It is clear from the above information that T(C0(Λ)) is relatively compact in (C0(Λ),‖.‖∞). On the other hand, the operator T is continuous. Indeed, let {vn} be a sequence in C0(Λ) that converges uniformly to v∈C0(Λ). Since we have
|˜ψ(y,cξ+HΛϕ+vn(y))|≤pM(y),∀(y,n)∈Λ×N, |
then, by combining the dominated convergence theorem and Proposition 3.1 with the continuity of ψ concerning the second variable, we deduce
{limn→∞Tvn(x)=Tv(x),∀x∈Λ,lim‖x‖→∞Tvn(x)=0, uniformly in n. |
From the equicontinuity of T(C0(Λ)), we deduce that the pointwise convergence implies uniform convergence. This fact implies that Tvn converges uniformly to Tv in C0(Λ). Which implies that T is continuous. By combining this fact and the fact that T(C0(Λ)) is relatively compact in (C0(Λ),‖.‖∞) with the Schauder fixed point theorem, we deduce that T has a fixed point v∈C0(Λ). Now if we put u=cξ+HΛf+v, then we get u∈C(Λ∪∂∞Λ) and
{limx→σ∈∂Λu(x)=f(σ),lim|x|→∞u(x)=c, |
moreover, u satisfies the following integral equation:
u=cξ+HΛf+V(˜ψ(.,u)) in Λ. |
Next, let us prove that the function u satisfies the following inequality
u_(x)≤u(x)≤¯u(x),∀x∈¯Λ. |
To do this, we proceed by contradiction, and we suppose that there exists x0∈Λ such that u(x0)<u_(x0). In this case, we define the following set:
E={x∈Λ such that u(x)<u_(x)}. |
Then, from our assumption, the set E is nonempty; moreover, from the fact that the functions u and u_ are continuous, we see that E is an open set in Λ. Moreover, we have u−u_=0 on ∂E. On the other hand, from Eq (3.1), we can see that for all x∈E, we have ˜ψ(x,u(x))=˜ψ(x,u_(x)). Hence, for any nonnegative function φ in D(E), we have
∫E(u(x)−u_(x))(−Δφ)(x)dx≥∫E(˜ψ(x,u(x))−˜ψ(x,u_(x)))φ(x)dx=0. |
The above information shows that u−u_ is a continuous superharmonic function in E with a boundary value equal to zero. On the other hand, from the definition of the sub-solution, in the case when E is unbounded, we have
lim|x|→∞,x∈E(u(x)−u_(x))=c−lim|x|→∞,x∈Eu_(x)≥0. |
Hence the maximum principle (see [22, p.397–398]) can be applied, and we conclude that u−u_≥0 in E. This contradicts the fact that E is nonempty. So u_≤u in ¯Λ. Similar arguments can be used to show that u≤¯u in ¯Λ. Now, the fact that u_(x)≤u(x)≤¯u(x), implies that for any x∈Λ we have ˜ψ(x,u(x))=ψ(x,u(x)). Finally, since u is a solution for the problem (3.1) and since ˜ψ(x,u(x))=ψ(x,u(x)), then u is also a solution for the problem (1.1). The proof is now finished.
We suppose that problems (1.3) and (1.4) admit a double pair of continuous sub-supersolutions (u_,v_) and (¯u,¯v). As in the proof of Theorem 1.1, for y∈C(¯Λ∪{∞}) and x∈¯Λ, we define
θ1(y)(x)={u_(x) if y(x)<u_(x),y(x) if u_(x)≤y(x)≤¯u(x),¯u(x) if ¯u(x)<y(x), |
and
θ2(y)(x)={v_(x) if y(x)<v_(x),y(x) if v_(x)≤y(x)≤¯v(x),¯u(x) if ¯v(x)<y(x). |
Clearly
θ1(y)∈[u_,¯u],and θ2(y)∈[v_,¯v]. |
On the other hand, by the fact that the functions u_,v_,¯u, and ¯v are in C(¯Λ∪{∞}), we deduce the existence of a positive constant C, such that for any x∈¯Λ, we have
{−C≤u_(x)≤¯u(x)≤C,−C≤v_(x)≤¯v(x)≤C. |
So, from (H3), there exist two nonnegative functions p1,C, p2,C in K(Λ) such that for all (x,y)∈Λ×C(¯Λ∪{∞}), we have
{|ψ1(x,θ1(y)(x),θ2(y)(x))|≤p1,C(x),|ψ2(x,θ1(y)(x),θ2(y)(x))|≤p2,C(x). |
Hence, from Proposition 3.1, we deduce that the map ω→ψi(ω,θ1(y)(ω),θ2(y)(ω)) is in K(Λ), moreover, for i∈{1,2}, the set
{V(ψi(.,θ1(y),θ2(y))):y∈C(¯Λ∪{∞})}, |
is relatively compact in (C0(Λ),‖.‖∞).
Now, we consider the Banach space C0(Λ)×C0(Λ), which is equipped with the following norm
‖(χ1,χ2)‖=‖χ1‖∞+‖χ2‖∞, |
and we define the operator T:C0(Λ)×C0(Λ)→C0(Λ)×C0(Λ) by
T(χ1,χ2)=(T1(χ1,χ2),T2(χ1,χ2)), |
where T1(χ1,χ2) and T2(χ1,χ2) are defined by
{T1(χ1,χ2)=Vψ1(.,θ1[c1ξ+HΛf1+χ1],θ2[c2ξ+HΛf2+χ2])(x),T2(χ1,χ2)=Vψ2(.,θ1[c1ξ+HΛf1+χ1],θ2[c2ξ+HΛf2+χ2])(x). |
We note that T1(χ1,χ2) and T2(χ1,χ2) are the unique pair of solutions to the following problem:
{−Δy=ψ1(x,θ1[c1ξ+HΛf1+χ1],θ2[c2ξ+HΛf2+χ2]) in Λ,−Δz=ψ1(x,θ1[c1ξ+HΛf1+χ1],θ2[c2ξ+HΛf2+χ2]) in Λ,y=0,z=0, on ∂Λlim|x|→∞y(x)=0 and lim|x|→∞z(x)=0. |
The same arguments as in the proof of Theorem 1.1, show that the set T(C0(Λ)×C0(Λ)) is relatively compact in C0(Λ)×C0(Λ), and moreover, the operator T is continuous. So, the Schauder fixed point theorem implies that T has a fixed point (χ1,χ2)∈C0(Λ)×C0(Λ).
Now, if we put
{u=c1ξ+HΛf1+χ1,v=c2ξ+HΛf2+χ2, |
then, u,v∈C(¯Λ∪{∞}), moreover, we have
{limx→σ∈∂Λu(x)=f1(σ),limx→σ∈∂Λv(x)=f2(σ), |
and
{lim|x|→∞u(x)=c1,lim|x|→∞v(x)=c2. |
Also, we have
{u(x)=c1ξ(x)+HΛf1(x)+Vψ1(.,θ1(u),θ2(v))(x),v(x)=c2ξ(x)+HΛf2(x)+Vψ2(.,θ1(u),θ2(v))(x). |
Finally, we will prove that u_≤u≤¯u and v_≤v≤¯v in ¯Λ. Since the proofs are similar for the cases u_≤u, u≤¯u, v_≤v, and v≤¯v, we will only prove that u_≤u. By contradiction, we assume that this is not true, so the set E1={x∈Λ such that u(x)<u_(x)} is nonempty. Since u and u_ are continuous, then E1 is open; u−u_=0 on ∂E1, and in the case when E1 is unbounded, we have
lim|x|→∞,x∈E1(u(x)−u_(x))=c−lim|x|→∞,x∈E1u_(x)≥0. |
Moreover since θ1(u)(x)=u_(x) for every x∈E1 and θ2(v)∈[v_,¯v], then for every φ∈C∞c(E1) with φ≥0 we have
∫E1[u(x)−u_(x)]Δφ(x)dx=∫E1u(x)Δφ(x)dx−∫E1u_(x)Δφ(x)dx=−∫E1ψ1(x,θ1(u)(x),θ2(v)(x))φ(x)dx−∫E1u_(x)Δφ(x)dx≥−∫E1ψ1(x,u_(x),θ2(v)(x))φ(x)dx+∫E1ψ1(x,u_(x),θ2(v)(x))φ(x)dx=0. |
The above information shows that u−u_ is a continuous super-harmonic function in E1 and satisfies u−u_≥0 on ∂E1 and lim|x|→∞,x∈E1u(x)−u_(x)≥0, if E1 is unbounded. Hence, from the maximum principle [22, p.397–398] we deduce that u−u_≥0 in E1. This contradicts the definition of E1 and so E1 is empty. Which proves that u_≤u in Λ. Consequently, we conclude that θ1(u)=u, θ2(v)=v, and (u,v) is a continuous weak solution of problem (1.3). The proof of Theorem 1.2 is now completed.
In this section, we present several applications of the main results. These applications approve and validate the main results of this paper. Since the case of the unbounded domain is more general than the bounded domain, throughout this section, we assume that Λ is a C1,1 exterior domain, and the letter i will denote the integer 1 or 2.
Application 1. Let α>0 and a be a nontrivial nonnegative function in B(K(Λ)), and consider the following problem:
{−Δu=λa(x)(1−uα) in Λ,u=0, on ∂Λ,lim|x|→∞u(x)=0. | (4.1) |
By the fact that a∈K(Λ), we deduce that Va∈C0(Λ). So if we put ¯u=λVa, then it is not difficult to see that for λ>0, the function ¯u is a continuous super-solution of problem (4.1). On the other hand, if we define the function f on [0,1‖Va‖∞) by
f(t)=t1−‖Va‖α∞tα. |
Then it is not difficult to see that the function f is differentiable and increasing on (0,1‖Va‖∞), and since f(0)=0 and limt→1‖Va‖∞f(t)=∞, then the function f:[0,1‖Va‖∞)→[0,∞) is a bijection.
Now, if λ>0, then there exists ε∈(0,1‖Va‖∞) such that 0<ε<f(ε)<λ, which implies that
ε<λ(1−εα‖Va‖α∞)<λ(1−(εVa)α). |
So, εa≤λa(1−(εVa)α) in Λ, and the function u_=εVa becomes a continuous weak sub-solution of (4.1) satisfying u_≤¯u. Hence, it follows from Theorem 1.1 that problem (4.1) has a continuous weak solution u satisfying εVa≤u≤λVa.
Application 2. Let α>0 and let a and b be two nontrivial nonnegative functions in B(K(Λ)), and we consider the following problem
{−Δu=λ(a(x)+b(x)uα) in Λ,u=0, on ∂Λ,lim|x|→∞u(x)=0. | (4.2) |
We begin by remarking that, similar to Application 1, for any λ>0, the function u_=λVa is a continuous sub-solution of (4.2). So to use the main theorems of this paper, we will find a positive continuous weak supersolution of (4.2).
Next, if we define the function g on [0,∞) by
g(t)=t1+‖V(a+b)‖α∞tα. |
Then a simple calculation shows that g is differentiable on (0,∞) and
g′(t)=1+(1−α)‖V(a+b)‖α∞tα(1+‖V(a+b)‖α∞tα)2. | (4.3) |
To discuss the monotonicity of the function g, we distinguish two cases.
Case 1: In this case, we consider the sublinear case, which means that α∈(0,1). In this case, we see from (4.3) that g is increasing, and so, it is a bijection from [0,∞) into [0,∞). Hence, for any λ>0, there exists M∈(0,∞) such that 0<λ<g(M)<M. Now, it is not difficult to see that the function ¯u=MV(a+b) is a continuous super-solution of (4.2), which satisfies in addition u_≤¯u. Finally, Theorem 1.1 implies that (4.2) has a continuous weak solution u; moreover, we have
λV(a)≤u≤MV(a+b). |
Case 2: In this case, we consider the super-solution case, which is the case when α≥1. In this case, easily, the function g is increasing on [0,t0] and decreasing on to [t0,∞) and satisfies g(0)=0 and limt→∞g(t)=0, where
t0=1(α−1)1α‖V(a+b)‖∞. |
Put
λ∗=g(t0)=(α−1)1−1αα‖V(a+b)‖∞. |
Then, for any λ in (0,λ∗], we see that λ<t0. So the function ¯u=t0V(a+b) is a continuous super-solution to problem (4.2), which satisfies in addition u_≤¯u. Again, Theorem 1.1 implies that (4.2) has a continuous weak solution satisfying
λVa≤u≤t0V(a+b). |
Application 3. For a positive real numbers α1,α2,β1,β2, and a nontrivial nonnegative functions a, b in B(K(Λ)), we consider the following system:
{−Δu=λa(x)(1−uα1vβ1) in Λ,−Δv=μb(x)(1−uα2vβ2) in Λ,u=v=0, on ∂Λ,lim|x|→∞u(x)=lim|x|→∞v(x)=0=0. | (4.4) |
We begin by remarking that for every λ>0, μ>0, the function
fi:[0,1‖Va‖αiαi+βi∞‖Vb‖βiαi+βi∞)→[0,∞), |
defined by:
fi(t)=t1−‖Va‖αi∞‖Vb‖βi∞tαi+βi, |
is an increasing bijection. So, there exists ε>0 such that
0<ε<fi(ε)<min(λ,μ). |
From which we deduce that
εa≤λa(1−(εVa)α1(εVb)β1), and εb≤μb(1−(εVa)α2(εVb)β2). |
This allows us to prove that the double pair (u_,v_)=(εVa,εVb), (¯u,¯v)=(λVa,μVb) is a continuous sub-supersolution of system (4.4). Hence, it follows from Theorem 1.2 that (4.4) has a continuous weak solution (u,v) satisfying εVa≤u≤λVa and εVb≤v≤μVb.
Application 4. For a positive real numbers α1,α2,β1,β2, and a nontrivial nonnegative functions a, b in B(K(Λ)), we consider the following system:
{−Δu=λ(a1(x)+b1(x)uα1vβ1) in Λ,−Δv=μ(a2(x)+b2(x)uα2vβ2) in Λ,u=v=0, on ∂Λ,lim|x|→∞u(x)=lim|x|→∞v(x)=0. | (4.5) |
We define the function gi on [0,∞) by
gi(t)=t1+‖V(a1+b1)‖αi∞‖V(a2+b2)‖βi∞tαi+βi. |
Clearly, gi is differentiable on [0,∞) and
g′i(t)=1+(1−αi−βi)‖V(a1+b1)‖αi∞‖V(a2+b2)‖βi∞tαi+βi(1+‖V(a1+b1)‖αi∞‖V(a2+b2)‖βi∞tαi+βi)2. |
Thus we will discuss four cases.
Case 1: 0<αi+βi<1. In this case, each gi is an increasing bijection from [0,∞) to [0,∞). Hence, for every 0<λ and 0<μ, there exists M∈(0,∞) such that 0<max(λ,μ)<min(g1(M),g2(M))<M. The double pair of continuous functions (u_,v_)=(λVa1,μVa2) and (¯u,¯v)=(MV(a1+b1),MV(a2+b2)) is a continuous sub-supersolution of system (4.5). Hence, it follows from Theorem 1.2 that (4.5) has a continuous weak solution (u,v) satisfying λVa1≤u≤MV(a1+b1) and μVa2≤v≤MV(a2+b2).
Case 2: αi+βi≥1. If we put
ti=1(αi+βi−1)1αi+βi‖V(a1+b1)‖αiαi+βi∞‖V(a2+b2)‖βiαi+βi∞. |
Then we see that gi is increasing on [0,ti] and decreasing on [ti,∞) and satisfies gi(0)=0 and limt→∞gi(t)=0. So, if we take λ∗=g1(t1) and μ∗=g2(t2). Then for every 0<λ≤λ∗ and any 0<μ≤μ∗, the double pair
(u_,v_)=(λVa1,μVa2)and(¯u,¯v)=(t1V(a1+b1),t2V(a2+b2)), |
is a continuous sub-supersolution of the system (4.5). Hence, it follows from Theorem 1.2 that (4.5) has a continuous weak solution (u,v) satisfying λVa1≤u≤t1V(a1+b1) and μVa2≤v≤t2V(a2+b2).
Case 3: 0<α1+β1<1 and α2+β2≥1. In this case, we obtain that for every 0<λ and 0<μ≤μ∗=g2(t2) the double pair (u_,v_)=(λVa1,μVa2) and (¯u,¯v)=(MV(a1+b1),t2V(a2+b2)) is a continuous sub-supersolution to system (4.5), where M is chosen so that 0<λ<g1(M). Hence, it follows from Theorem 1.2 that (4.5) has a continuous weak solution (u,v); moreover, we have
λVa1≤u≤MV(a1+b1)andμVa2≤v≤t2V(a2+b2). |
Case 4: 0<α2+β2<1, and α1+β1≥1. Inspired by cases 2 and 3, we can prove that for every 0<λ≤λ∗=g1(t1), and for every μ>0, the double pair
(u_,v_)=(λVa1,μVa2)and(¯u,¯v)=(t1V(a1+b1),MV(a2+b2)), |
is a continuous sub-supersolution of system (4.5), where M is chosen so that 0<μ<g2(M). Hence, it follows from Theorem 1.2 that (4.5) has a continuous weak solution (u,v) satisfying in addition
λVa1≤u≤t1V(a1+b1)andμVa2≤v≤MV(a2+b2). |
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This work was funded by the University of Jeddah, Saudi Arabia, under grant No. UJ-22-DR-56. The authors, therefore, acknowledge with thanks the University's technical and financial support.
The authors declare there is no conflict of interest.
[1] |
H. Zou, A priori estimates for a semilinear elliptic system without variational structure and their applications, Math. Ann., 323 (2002), 713–735. https://doi.org/10.1007/s002080200324 doi: 10.1007/s002080200324
![]() |
[2] |
K. Akô, On the Dirichlet problem for quasi-linear elliptic differential equations of the second order, J. Math. Soc. Jpn., 13 (1961), 45–62. https://doi.org/10.2969/jmsj/01310045 doi: 10.2969/jmsj/01310045
![]() |
[3] |
R. S. Alsaedi, H. Mâagli, V. D. Rădulescu, N. Zeddini, Asymptotic behavior of positive large solutions of quasilinear logistic problems, Electron. J. Qual. Theory Differ. Equations, 28 (2015), 1–15. https://doi.org/10.14232/ejqtde.2015.1.28 doi: 10.14232/ejqtde.2015.1.28
![]() |
[4] |
H. Amann, On the existence of positive solutions of nonlinear elliptic boundary value problems, Indiana Univ. Math. J., 21 (1972), 125–146. https://doi.org/10.1512/iumj.1972.21.21012 doi: 10.1512/iumj.1972.21.21012
![]() |
[5] | P. Clément, G. Sweers, Getting a solution between sub- and supersolutions without monotone iteration, Rend. Ist. Mat. Univ. Trieste, 19 (1987), 189–194. Available from: http://hdl.handle.net/10077/4937. |
[6] |
S. Cui, Existence and nonexistence of positive solutions for singular semilinear elliptic boundary value problems, Nonlinear Anal. Theory Methods Appl., 41 (2000), 149–176. https://doi.org/10.1016/S0362-546X(98)00271-5 doi: 10.1016/S0362-546X(98)00271-5
![]() |
[7] |
H. B. Keller, Elliptic boundary value problems suggested by nonlinear diffusion processes, Arch. Ration. Mech. Anal., 35 (1969), 363–381. https://doi.org/10.1007/BF00247683 doi: 10.1007/BF00247683
![]() |
[8] |
M. Montenegro, A. C. Ponce, The sub-supersolution method for weak solutions, Proc. Am. Math. Soc., 136 (2008), 2429–2438. https://doi.org/10.1090/S0002-9939-08-09231-9 doi: 10.1090/S0002-9939-08-09231-9
![]() |
[9] |
M. Montenegro, A. Suárez, Existence of a positive solution for a singular system, Proc. R. Soc. Edinburgh Sect. A: Math., 140 (2010), 435–447. https://doi.org/10.1017/S0308210509000705 doi: 10.1017/S0308210509000705
![]() |
[10] |
E. S. Noussair, On the existence of solutions of nonlinear elliptic boundary value problems, J. Differ. Equations, 34 (1979), 482–495. https://doi.org/10.1016/0022-0396(79)90032-9 doi: 10.1016/0022-0396(79)90032-9
![]() |
[11] | E. S. Noussair, C. A. Swanson, Positive solutions of semilinear Schrödinger equations in exterior domains, Indiana Univ. Math. J., 28 (1979), 993–1003. Available from: https://www.jstor.org/stable/24892551. |
[12] | A. Ogata, On bounded positive solutions of nonlinear elliptic boundary value problems in an exterior domain, Funkcialaj Ekvacioj, 17 (1974), 207–222. |
[13] |
V. D. Rădulescu, D. Repovš, Combined effects in nonlinear problems arising in the study of anisotropic continuous media, Nonlinear Anal. Theory Methods Appl., 75 (2012), 1524–1530. https://doi.org/10.1016/j.na.2011.01.037 doi: 10.1016/j.na.2011.01.037
![]() |
[14] |
I. Bachar, H. Mâagli, N. Zeddini, Estimates on the green function and existence of positive solutions of nonlinear singular elliptic equations, Commun. Contemp. Math., 5 (2003), 401–434. https://doi.org/10.1142/S0219199703001038 doi: 10.1142/S0219199703001038
![]() |
[15] |
H. Mâagli, M. Zribi, On a new Kato class and singular solutions of a nonlinear elliptic equation in bounded domains of Rn, Positivity, 9 (2005), 667–686. https://doi.org/10.1007/s11117-005-2782-z doi: 10.1007/s11117-005-2782-z
![]() |
[16] |
R. S. Alsaedi, H. Mâagli, N. Zeddini, Positive solutions for some competitive elliptic systems, Math. Slovaca, 64 (2014), 61–72. https://doi.org/10.2478/s12175-013-0187-1 doi: 10.2478/s12175-013-0187-1
![]() |
[17] |
A. Ghanmi, H. Mâagli, S. Turki, N. Zeddini, Existence of positive bounded solutions for some nonlinear elliptic systems, J. Math. Anal. Appl., 352 (2009), 440–448. https://doi.org/10.1016/j.jmaa.2008.04.029 doi: 10.1016/j.jmaa.2008.04.029
![]() |
[18] |
N. Zeddini, R. S. Sari, Existence of positive continuous weak solutions for some semilinear elliptic eigenvalue problems, Opuscula Math., 42 (2022), 489–519. https://doi.org/10.7494/OpMath.2022.42.3.489 doi: 10.7494/OpMath.2022.42.3.489
![]() |
[19] | K. L. Chung, Z. Zhao, From Brownian motion to Schrödinger's equation, Bull. Am. Math. Soc., 39 (2001), 109–111. Available from: https://www.ams.org/journals/bull/2002-39-01/S0273-0979-01-00925-9/S0273-0979-01-00925-9.pdf. |
[20] |
R. Guefaifia, G. C. G. dos Santos, T. Bouali, R. Jan, S. Boulaaras, A. Alharbi, Sub-super solutions method combined with schauder's fixed point for existence of positive weak solutions for anisotropic non-local elliptic systems, Mathematics, 10 (2022), 4479. https://doi.org/10.3390/math10234479 doi: 10.3390/math10234479
![]() |
[21] | C. V. Pao, Nonlinear Parabolic and Elliptic Equations, Plenum Press, New York, 1992. |
[22] | R. Dautray, J. L. Lions, Mathematical Analysis and Numerical Methods for Science and Technology-Volume 1-Physical Origins and Classical Methods, Springer-Verlag, 1999. |
1. | Abdujabar Rasulov, Gulnora Raimova, D. Rudoy, A.N. Altybaev, M. Petkovich, N. Miletic, Application of branching processes to solve boundary value problems for the semi-linear Helmholtz equation, 2024, 583, 2267-1242, 07014, 10.1051/e3sconf/202458307014 |