Semilinear problems with poly-Laplace type operators

Abstract

The paper deals with semilinear operator equations involving iterates of a strongly monotone symmetric linear operator. In particular there are consider semilinear polyharmonic equations subject to the Navier boundary conditions. A careful analysis is made on the energetic spaces associated to such problems and a number of existence results are obtained by using a fixed point approach.

Authors

Radu Precup
Babes-Bolyai University, Cluj-Napoca, Romania
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy

Keywords

polyharmonic equation, iterates of symmetric linear operators, energetic space

Paper coordinates

R. Precup, Semilinear problems with poly-Laplace type operators, Proceedings of the Romanian Academy Series A, 23 (2022) no. 4, pp. 319-328.

PDF

About this paper

Print ISSN

1454-9069

Online ISSN

1454-9069

google scholar link

Paper (preprint) in HTML form

Semilinear problems with poly-Laplace type operators

Semilinear problems with poly-Laplace type operators

Radu Precup
Faculty of Mathematics and Computer Science and Institute of
Advanced Studies in Science and Technology,
Babeş-Bolyai University, 400084 Cluj-Napoca, Romania &
Tiberiu Popoviciu  Institute of Numerical Analysis, Romanian
Academy, P.O. Box 68-1, 400110 Cluj-Napoca, Romania
Abstract

The paper deals with semilinear operator equations involving iterates of a strongly monotone symmetric linear operator. In particular there are consider semilinear polyharmonic equations subject to the Navier boundary conditions. A careful analysis is made on the energetic spaces associated to such problems and a number of existence results are obtained by using a fixed point approach.


Mathematics Classification: 35J40, 46E35, 47J25

Key words: Polyharmonic equation, iterates of symmetric linear operators, energetic space.

1 Introduction

There is known the bi-Laplace equation Δ2u=0 whose solutions are called biharmonic functions. The equation arises as a model for the elastic equilibrium in the theory of elasticity. Also there are known its generalizations, the poly-Laplace equations Δpu=0, p>2, whose solutions are said to be polyharmonic of order p (see, [9] and [10]). The operator Δ2=ΔΔ is referred as the bi-Laplacian and Δp=Δ(Δp1) is said to be the Laplacian of order p. The non-homogeneous versions of these equations are

Δpu=h,

and when considered in a domain Ω, in case p=2, there has been added the boundary condition

u=dudν=0on Ω, (1.1)

where ν is the unit normal vector to the boundary, or the boundary condition

u=Δu=0on Ω (1.2)

(see [8]). For p>2, condition (1.1) can be generalized following Lauricella [7] as follows

u=dudν==dp1udνp1=0on Ω,

and (1.2), as suggested by Riquier [13], by

u=Δu=Δ2u==Δp1u=0on Ω.

For the classical theory of polyharmonic functions we refer the reader to the volume [9] which brings together the entire contribution of Miron Nicolescu to this field and which allows obtaining information on contributions originating from old, less accessible publications.

Modern theory has introduced the concept of weak solution and Sobolev spaces as natural framework for the study of these operators and of the associated semilinear problems. Thus the problem

{Δ2u=f(x,u,u,Δu)in Ωu=uν=0on Ω (1.3)

involving the natural boundary condition uν=0can be naturally addressed in H02(Ω) endowed with the equivalent norm |Δu|L2(Ω). Other studies (see, e.g., [1], [2], [4] and [11]) have aimed to treat problems of type (1.3) under the boundary conditions u=Δu=0 on Ω (called Navier boundary conditions [5]) by looking for solutions in the space H2(Ω)H01(Ω) with norm |Δu|L2(Ω). The problem is that the condition ”Δu=0on Ω” not being a natural boundary condition it does not follow from the variational formulation of the problem. This is the reason to restrict the study to a subspace of functions in order to give a meaning to the equality Δu=0 on the boundary. This will be one of our main goals in this work. Roughly speaking we suggest that the iterative nature of the differential operator to reflect on its energetic space and consequently on some basic inequalities. We lead this discussion more generally considering instead of Laplacian a strongly monotone symmetric linear operator A. Thus our results will concern semilinear operator equations of the form

Apu=h+F(u,Au,,Ap1u),

where Ap is the p-th iterate of A, defined recursively by Ap=AAp1. The whole approach is based on the theory of the energetic space XA associated to A. There are thus obtained existence results for the problem

{Apu=h+F(u,Au,,Ap1u)u,Au,,Ap1uXA (1.4)

where hXA is given and F is on the position of a perturbation of h. In particular, we obtain results for semilinear poly-Laplace equations.

2 Preliminaries

In this section we recall the notion of energetic space (see [14]) and some related results.

2.1 The energetic space

Let X be a real Hilbert space with the inner product (,)X and the norm ||X. Let Y be a linear subspace of X and A:YX be a strongly monotone symmetric linear operator, that is, a linear operator satisfying

(Au,v)X=(u,Av)Xfor all u,vY,
(Au,u)Xc2|u|X2for all uY (2.1)

and some constant c>0. Then, endowed with the energetic inner product

(u,v)A:=(Au,v)X(u,vY)

and the energetic norm

|u|A=(Au,u)X1/2(uY),

Y becomes a pre-Hilbert space. Its completion (see, e.g., [6, Section I 4.3]) denoted by XA is called the energetic space of A. In virtue of (2.1), any Cauchy sequence in the energetic norm is also a Cauchy sequence in the norm ||X. This allows us to see XA as a subset of the original complete space X, and the elements of XA as limits in X of Cauchy sequences from (Y,||A). Furthemore, the energetic inner product and norm can be extended from Y to XA by

(u,v)XA:=limk(uk,vk)A,|u|XA:=lim|uk|A,

where (uk) and (vk) are Cauchy sequences in (Y,||A) that converge in X to u and v, respectively.

2.2 Abstract Poincaré’s inequality

Inequality (2.1) can be extended by density from Y to XA showing that

|u|XAc|u|Xfor all uXA. (2.2)

Thus c is an embedding constant for the continuous inclusion XAX. We call this inequality Poincaré’s inequality.

If the embedding XAX is compact, then there is a largest embedding constant c and the inequality is reached. Indeed, if we denote

λ:=inf{|u|XA2:uXA,|u|X=1},

then λc2 and if we take any minimizing sequence (uk), that is

ukXA,|uk|X=1,|uk|XA2λ,

then using the compactness of the embedding XAX and passing eventually to a subsequence we can assume that uku in X, for some uX. Furthermore, from the identity

|ukum|XA2+|uk+um|XA2=2(|uk|XA2+|um|XA2),

since |uk+um|XA2λ|uk+um|X2, we deduce

|ukum|XA22(|uk|XA2+|um|XA2)λ|uk+um|X20

as k,m. Hence (uk) is a Cauchy sequence in XA. Let vXA be such that ukv in XA. Then ukv in X too, and the uniqueness of the limit implies that v=u. Consequently, |uk|XA|u|XA, that is λ=|u|XA2. Thus the infimum λ is reached and λ is the best constant c in (2.2). Thus, in case that the embedding XAX is compact, Poincaré’s inequality reads as follows:

|u|X1λ|u|XAfor all uXA.

2.3 The dual of the energetic space

Having XAX, for the dual spaces we have XXA and if, based on Riesz’ theorem, we assume the identification X=X, then one has

XAXXA.

In addition, from (2.2) we also have

|u|Xc|u|XAfor all uX. (2.3)

Indeed, if uX, then  for any vX, one hasu,v=(u,v)X, where by  , we mean the value of a linear functional at a given element. Then

|u|XA=supvXA{0}|u,v||v|XA=supvXA{0}|(u,v)X||v|XAsupvXA{0}|u|X|v|X|v|XA1c|u|X.

Notice in case that the embedding XAX is compact, so is the embedding XXA and in (2.3) we may take as in (2.2) the best constant c=λ.

2.4 Extension of operator A

Clearly we can define the linear operator A~:XAXA by

A~u,v=(u,v)XAfor all u,vXA.

In particular, if u,vY, then since (u,v)XA=(Au,v)X, one has A~u,v=(Au,v)X, which by the density of Y into XA can be extended to all vXA. Thus the functionals A~u and (Au,)X act identically in XA. The last one is a continuous linear functional on X which in virtue of Riesz’s representation theorem is identified with Au. In this sense, as continuous linear functionals on XA, one has A~u=Au, and therefore A~ can be seen as an extension of A from Y to XA. It is common to use the same symbol A for the extension A~. Thus A:XAXA.

2.5 The inverse of operator A

In the previous subsection we have that for every uXA there is a unique element denoted AuXA with

Au,v=(u,v)XAfor all vXA. (2.4)

Conversely, for every hXA by Riesz’s theorem, there is a unique element uXA with

h,v=(u,v)XAfor all vXA.

Clearly Au=h and thus u=A1h. Hence the inverse of A is the operator A1:XAXA defined by

(A1h,v)XA=h,vfor all vXA.

The two linear operators A and A1 are isometries between XA and XA. Indeed, letting v=u in (2.4) gives

|u|XA2=Au,u|Au|XA|u|XA,

whence |u|XA|Au|XA. The converse inequality comes from

|Au|XA=supvXA{0}|Au,v||v|XA=supvXA{0}|(u,v)XA||v|XA|u|XA.

Hence

|Au|XA=|u|XA(uXA),|A1h|XA=|h|XA(hXA).

2.6 Weak solutions to linear operator equations

Consider the operator equation associated to A,

Au=h.

By a (strong) solution we mean an element uY such that Au=h. Obviously this is possible if hX. By a weak solution we mean an element uXA satisfying the identity

(u,v)XA=h,vfor all vXA.

When speaking about weak solutions we may assume more generally that hXA. In view of the previous subsection, for each hXA, the equation has a unique weak solution, namely u=A1h.

Note that looking for weak solutions to a semilinear equation

Au=Φ(u),

where Φ:XAXA is any mapping, reduces to solving the fixed point equation

u=A1Φ(u),uXA.

3 Semilinear operator equations involving iterates of a symmetric linear operator

We now come back to problem (1.4), where A is a linear operator as in Introduction.

3.1
Functional framework

Looking at the required conditions on the elements u,Au,,Ap1u to belong to the energetic space XA of the operator A, we may seek solutions in the space

H:=A(p1)(XA).

Here Ak=A1(A(k1)) for k=2,,p1. Since A1:XAXA and XAXA one has

H=A(p1)(XA)A(p2)(XA)A1(XA)XA. (3.1)

We endow H with the inner product and norm

(u,v)H:=(Ap1u,Ap1v)XA,|u|H:=|Ap1u|XA.

Note that the functional ||H is indeed a norm on H since if for some uH, one has |u|H=0, then Ap1u=0, whence Ap2u=0 and so on until we obtain u=0.

Lemma 3.1

The space H endowed with the inner product (,)H is a Hilbert space which continuously embeds in XA.

Proof. Let (uk) be any Cauchy sequence in H. Then (Ap1uk) is Cauchy in XA, so convergent in XA to some vXA. Since the embedding XAXA is continuous, we then have Ap1ukv in XA. Next, the continuity of A1 from XA to XA implies Ap2ukA1v in XA. Repeating the above reasoning we arrive to the conclusion that ukA(p1)v in XA, that is uku:=A(p1)v in H. This proves that (H,||H) is complete.   

Knowing the operator A1 from XA to XA and the inclusions (3.1) we immediately can see that for every hXA there is a unique uH, namely u=Aph,   which solves the non-homogeneous equation

Apu=h.

Consequently, solving a semi-linear equation of the form

Apu=Φ(u),

where Φ:HXA is any mapping, is equivalent to the fixed point equation

u=ApΦ(u),uH

for the operator ApΦ:HH.

3.2 Existence and uniqueness under a Lipschitz condition

Using Banach contraction principle we obtain the following result on problem (1.4).

Theorem 3.2

Let F:XpX satisfy

|F(u)F(v)|Xi=1pai|uivi|X (3.2)

for all u=(u1,,up),v=(v1,,vp)Xp and some nonnegative constants ai,i=1,,p. If

θ:=i=1paic2(p+1i)<1, (3.3)

then problem (1.4) has a unique solution uH.

Proof. Problem (1.4) is equivalent to the fixed point equation

u=Ap(h+F(u,Au,,Ap1u)),uH.

Using (3.2) and Poincaré’s inequality (2.2), for any u,vH, we have

|ApF(u,Au,,Ap1u)ApF(v,Av,,Ap1v)|H
= |A1(F(u,Au,,Ap1u)F(v,Av,,Ap1v))|XA
= |F(u,Au,,Ap1u)F(v,Av,,Ap1v)|XA
1c|F(u,Au,,Ap1u)F(v,Av,,Ap1v)|X
1ci=1pai|Ai1(uv)|X.

Furthermore, for ant wH, one has

|Ap1w|X 1c|Ap1w|XA=1c|w|H,
|Ap2w|X 1c|Ap2w|XA=1c|Ap1w|XA1c2|Ap1w|X1c3|w|H,
|Ap3w|X 1c|Ap3w|XA=1c|Ap2w|XA1c2|Ap2w|X1c5|w|H.

Repeating the above estimations for p4,, 0, we obtain

|Ai1w|X1c2(pi)+1|w|H,i=1,,p. (3.5)

Then

|ApF(u,Au,,Ap1u)ApF(v,Av,,Ap1v)|H
(i=1paic2(p+1i))|uv|H
= θ|uv|H,

which in view of (3.3) shows that the operator

N:=Ap(h+F(u,Au,,Ap1u))

is a contraction on H. The conclusion now follows from Banach contraction principle.   

3.3 Existence under a linear growth condition

If instead of the Lipschitz condition (3.2) we only have a linear growth condition on F and we assume that the embedding XAX is compact, then we can still prove the existence of at least one solution by using Schauder’s fixed point theorem.

Theorem 3.3

Assume that the embedding XAX is compact and that F:XpX is continuous and satisfies

|F(u)|XC+i=1pai|ui|X (3.7)

for all u=(u1,,up)Xp and some C>0 and nonnegative constants ai,i=1,,p. If condition (3.3) holds, then problem (1.4) has at least one solution uH with

|u|HCc1+|h|XA1θ.

Proof. As above we now have

|ApF(u,Au,,Ap1u)|H 1c|F(u,Au,,Ap1u)|X
Cc1+θ|u|H.

Since θ<1, R=(Cc1+|h|XA)/(1θ)>0 and N is a self mapping of the closed ball BR of H centered at the origin and of radius R. On the other hand N0(u):=ApF(u,Au,,Ap1u) can be decomposed as

N=A(p1)A1JFJ0P,

where

P : HXAp,Pu=(u,Au,,Ap1u);J0:XApXp,Ju=u;
F : XpX;J:XXA,Ju=u;A1:XAXA;A(p1):XAH.

All these operators are continuous and bounded (send bounded sets to bounded sets) and J0 is compact. As a result their composition N0 is completely continuous. Now the conclusion follows from Schauder’s fixed point theorem applied to N in the ball BR.   

3.4 Existence via a priori bounds

We may replace the growth condition on F by a sign type condition as shows the following theorem.

Theorem 3.4

Assume that F:XApXA is completely continuous and satisfies

F(v),vpα|vp|X2 (3.8)

for all v=(v1,,vp)XAp and some α[0,c2). Then problem (1.4) has at least one solution uH. Moreover, any solution uH of the problem satisfies

|u|H|h|XA/(1αc2). (3.9)

Proof. Using a similar reasoning as in the previous proof we can show that the operator ApF:HH is completely continuous. We now prove that the set of all possible solutions of the equations

u=μAp(h+F(u,Au,,Ap1u))

for μ[0,1] is bounded as (3.9) shows. Indeed, if u is such a solution, then

|u|H2 = μ(Aph,u)H+μ(ApF(u,Au,,Ap1u),u)H
= μ(A1h,Ap1u)XA+μ(A1F,Ap1u)XA
|h|XA|u|H+μ(A1F,Ap1u)XA.

Next since v:=(u,Au,,Ap1u)XAp, based on (3.8), one has

(A1F,Ap1u)XA=F,Ap1uα|Ap1u|X2αc2|u|H2.

It follows that

|u|H2αc2|u|H2+|h|XA|u|H,

whence (3.9). The existence of a solution is guaranteed by the Leray-Schauder principle.   

4 Semilinear problems with poly-Laplace operators

The results established in Section 3 can be easily applied to problems involving poly-Laplace operators, more exactly to the problem

{Δpu=h+f(x,u,Δu,,Δp1u)in Ωu=Δu==Δp1u=0on Ω. (4.1)

Here Ωn is bounded open, X=L2(Ω), A=Δ, XA=H01(Ω),XA=H1(Ω), hH1(Ω)andf:Ω×p. Hence

H=(Δ)(p1)H01(Ω),(u,v)H=ΩΔp1uΔp1v,|u|H=|Δp1u|H01(Ω).

Also the compactness of the imbedding XAX holds and the imbedding constant in Poincaré’s inequality is c=λ1, where λ1 is the first eigenvalue of the Dirichlet problem for Laplacian (for the theory of elliptic problems, see, e.g., [3] or [12]).

In this case, F is the superposition operator

F(u,(Δ)u,,(Δ)p1u)(x)=f(x,u(x),Δu(x),,Δp1u(x))

(xΩ,uH).

Theorem 3.2 yields the following result.

Corollary 4.1

Let f satisfy the Carathéodory conditions, f(,0)L2(Ω) and

|f(x,u)f(x,v)|i=1pai|uivi| (4.2)

for all u=(u1,,up),v=(v1,,vp)p and some nonnegative constants ai,i=1,,p. If θ<1(θ being given by (3.3) with c=λ1), then problem (4.1) has a unique solution u(Δ)(p1)H01(Ω).

Proof. According to the main theorem abut Nemytskii’s superposition operator, F maps L2(Ω;p) to L2(Ω). In addition, for any u,vL2(Ω;p), from (4.2) we find

|F(u)F(v)|L2(Ω)i=1pai|uivi|L2(Ω).

Thus Theorem 3.2 is applicable and gives the result.   

Theorem 3.3 yields the following result.

Corollary 4.2

Let f satisfy the Carathéodory conditions and

|f(x,u)|ψ(x)+i=1pai|ui| (4.3)

for all u=(u1,,up)p, a.e. xΩ, some nonnegative constants ai,i=1,,p and a function ψL2(Ω). If θ<1(θ being given by (3.3) with c=λ1), then problem (4.1) has at least one solution u(Δ)(p1)H01(Ω) with

|Δp1u|H01(Ω)|ψ|L2(Ω)/λ1+|h|H1(Ω)1θ.

Using Theorem 3.4 we obtain the following result.

Corollary 4.3

Let f satisfy the Carathéodory conditions and

|f(x,u)|ψ(x)+i=1pai|ui|q (4.4)

for all u=(u1,,up)p, a.e. xΩ, some nonnegative constants ai,i=1,,p, a number 1q<21=2/(2)=(n+2)/(n2)(n3)) and a function ψLq0(Ω), where q0((2),2/q]. In addition assume that

vpf(x,v)αvp2

for every vp and some α[0,λ1). Then problem (4.1) has at least one solution u(Δ)(p1)H01(Ω). Moreover, any solution u(Δ)(p1)H01(Ω) satisfies

|Δp1u|H01(Ω)|h|H1(Ω)/(1αλ1).

Proof. Let q1=q0q. Clearly q1[1,2].Hence the embedding H01(Ω)Lq1(Ω) is continuous, while since q0>(2), the embedding Lq0(Ω)H1(Ω) is compact. In addition since q=q1/q0, from (4.4) we have that Nemytskii’s superposition operator Nf is well-defined, continuous and bounded from Lq1(Ω)p to Lq0(Ω). Then our operator F(u)=f(,u()) can be decomposed as F=JNfP, where

P : H01(Ω)pLq1(Ω)p,Pu=u;
Nf : Lq1(Ω)pLq0(Ω),Nf(v)(x)=f(x,v(x));
J : Lq0(Ω)H1(Ω),J1u=u.

Since J is compact one deduces that F:H01(Ω)pH1(Ω) is completely continuous.

We now check condition (3.8). For vH01(Ω)p, one has

F(v),vp=JNfP(v),vp=Ωvp(x)f(x,v(x))αΩvp(x)2=α|vp|L2(Ω)2.

Hence the assumptions of Theorem 3.4 are fulfilled and the conclusion follows.   

In contrast with the general case of equations involving iterates of a linear operator A, the case of the Laplace operator is a special one due to the representation of the Laplacian Δ=as a composition of two differential operators, the gradient and the divergence. This particularity allows nonlinear terms of semilinear equations also to depend on gradient. Thus, instead of problem (4.1) we can consider more generally the problem

{Δpu=h+f(x,u,Δu,,Δp1u,u,Δu,,Δp1u)in Ωu=Δu==Δp1u=0on Ω. (4.5)

Then looking to extend to this problem the results in Corollaries 4.1 and 4.2, the expression of constant θ in (3.3) should be completed by terms involving odd powers of 1/c (1/λ1). For example, if f:Ω×p×np is such that

|f(x,u,𝐮)f(x,v,𝐯)|i=1p(ai|uivi|+bi|𝐮i𝐯i|)

for all u,vp and 𝐮,𝐯(n)p (where applied to vectors from n, notation || stands for the Euclidian norm), then trying to follow the estimation made for (3.2) we arrive to the final sum

1ci=1p(ai|Δi1(uv)|L2(Ω)+bi|Δi1(uv)|L2(Ω;n)).

According to (3.5) we have

|Δi1w|L2(Ω)1c2(pi)+1|w|H,i=1,,p,

which help in the estimation

|Δi1w|L2(Ω;n)=|Δi1w|H01(Ω)=|Δiw|H1(Ω)1c|Δiw|L2(Ω)1c2(pi)|w|H.

Then the analogue of (3.2) for the new operator

N0(u):=Δpf(,u,Δu,,Δp1u,u,Δu,,Δp1u),

is the estimate

|N0(u)N0(v)|Hi=1p(aic2(pi)+2+bic2(pi)+1)|uv|H.

Thus the contraction condition guaranteeing the existence and uniqueness of the solution of (4.5) is now

θ~:=i=1p(aic2(pi)+2+bic2(pi)+1)<1.

An analogue result to Corollary 4.2 can be established under the growth condition on f,

|f(x,u,𝐮)|C+i=1p(ai|ui|+bi|𝐮i|)

and the same condition θ~<1 on the constants ai and bi.

References

  • [1] F. Bernis, J. Garcia-Azorebo and I. Peral, Existence and multiplicity of nontrivial solutions in semilinear critical problems of fourth order, Adv. Differential Equations 1 (1996), 210–240.
  • [2] M. Bhakta, Solutions to semilinear elliptic PDE’s with biharmonic operator and singular potential, Electronic J. Differential Equations, 2016 (2016), 261, 1–17.
  • [3] H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer, New York, 2011.
  • [4] X. Cheng, Z. Feng and L. Wei, Existence and multiplicity of nontrivial solutions for a semilinear biharmonic equation with weight functions, Discrete Cont. Dyn. Syst. Ser. S, 14 (2021), 3067–3083.
  • [5] F. Gazzola, H.-C. Grunau and G. Sweers, Polyharmonic Boundary Value Problems, Springer, Berlin, 2009.
  • [6] L.V. Kantorovich and G.P. Akilov, Functional Analysis, Pergamon Press, Oxford, 1982.
  • [7] G. Lauricella, Integrazione dell’equazione Δ2(Δ2u)=0 in un campo di forma circolare, Atti Accad. Sci. Torino Cl. Sci. Fis. Mat. Natur. 31 (1895-96), p. 1010.
  • [8] E. Mathieu, Mémoire sur l’équation aux différences partielles du quatrième ordre ΔΔu=0 et sur l’équilibre d’élasticité d’un corps solide, J. Math. Pures Appl. 2e série, 14 (1869), 378–421.
  • [9] M. Nicolescu, Opera matematică. Funcţii poliarmonice, Ed. Academiei, Bucureşti, 1980.
  • [10] M. Nicolescu, Les Fonctions Polyharmoniques, Actualité Sci. 331, Paris, Herman, 1936.
  • [11] M. Pérez-Llanos and A. Primo, Semilinear biharmonic problems with a singular term, J. Differential Equations 257 (2014), 3200–3225.
  • [12] R. Precup, Linear and Semilinear Partial Differential Equations, De Gruyter, Berlin, 2013.
  • [13] Ch. Riquier, Sur quelques problèmes relatifs à l’équation aux dérivées partielles Δnu=0, J. Math. Pures Appl. (9) 5 (1926), 297–394.
  • [14] E. Zeidler, Applied Functional Analysis: Applications to Mathematical Physics, Springer, New York, 1995.

[1] F. BERNIS, J. GARCIA-AZOREBO, I. PERAL, Existence and multiplicity of nontrivial solutions in semilinear critical problems of fourth order, Adv. Differential Equations, 1, pp. 210–240, 1996.
[2] M. BHAKTA, Solutions to semilinear elliptic PDE’s with biharmonic operator and singular potential, Electronic J. Differential Equations, 2016, 261, pp. 1–17, 2016.
[3] H. BREZIS, Functional analysis, Sobolev spaces and partial differential equations, Springer, 2011.
[4] X. CHENG, Z. FENG, L. WEI, Existence and multiplicity of nontrivial solutions for a semilinear biharmonic equation with weight functions, Discrete Cont. Dyn. Syst. Ser. S, 14, pp. 3067–3083, 2021.
[5] F. GAZZOLA, H.-C. GRUNAU, G. SWEERS, Polyharmonic boundary value problems, Springer, 2009.
[6] L.V. KANTOROVICH, G.P. AKILOV, Functional analysis, Pergamon Press, 1982.
[7] G. LAURICELLA, Integrazione dell’equazione ∆2(∆2u) = 0 in un campo di forma circolare, Atti Accad. Sci. Torino Cl. Sci. Fis. Mat. Natur., 31, pp. 1010, 1895-96.
[8] E. MATHIEU, Memoire sur l’ equation aux diff erences partielles du quatrieme ordre  ∆∆u = 0 et sur l’equilibre d’ elasticite d’un  corps solide, J. Math. Pures Appl. 2e serie, 14, pp. 378–421, 1869.
[9] M. NICOLESCU, Opera matematica. Functii poliarmonice , Ed. Academiei, Bucuresti, 1980.
[10] M. NICOLESCU, Les fonctions polyharmoniques, Actualite Sci.  331, Herman, 1936.
[11] M. PEREZ-LLANOS, A. PRIMO, Semilinear biharmonic problems with a singular term, J. Differential Equations, 257, pp. 3200–3225, 2014.
[12] R. PRECUP, Linear and semilinear partial differential equations, De Gruyter, 2013.
[13] CH. RIQUIER, Sur quelques problemes relatifs a l’ equation aux derivees partielles  ∆nu = 0, J. Math. Pures Appl., 5, 9, pp. 297–394, 1926.
[14] E. ZEIDLER, Applied functional analysis: applications to mathematical physics, Springer, 1995

2022

Related Posts