Stationary solutions of Fokker-Planck equations with nonlinear reaction terms in bounded domains

Abstract

Using an operator approach, we discuss stationary solutions to Fokker-Planck equations and systems with nonlinear reaction terms. The existence of solutions is obtained by using Banach, Schauder and Schaefer fixed point theorems, and for systems by means of Perov’s fixed point theorem. Using the Ekeland variational principle, it is proved that the unique solution of the problem minimizes the energy functional, and in case of a system that it is the Nash equilibrium of the energy functionals associated to the component equations.

Authors

Radu Precup
Department of Mathematics Babes-Bolyai University, Cluj-Napoca, Romania

Paola Rubbioni
Department of Mathematics and Computer Science, University of Perugia, Perugia, Italy

Keywords

Elliptic equation; Reaction-diffusion equation; Semi-linear Fokker-Planck equation; Fixed point; Variational method; Nash type equilibrium

Paper coordinates

R. Precup, P. Rubbioni, Stationary solutions of Fokker-Planck equations with nonlinear reaction terms in bounded domains, Potential Analysis, 57 (2022), 181–199, https://doi.org/10.1007/s11118-021-09911-6

PDF

??

About this paper

Journal

Potential Analysis

Publisher Name

Springer

Print ISSN

0926-2601

Online ISSN

1572-929X

google scholar link

[1] Barbu, V.: Generalized solutions to nonlinear Fokker-Planck equations. J. Differ. Equat. 261, 2446–2471 (2016), MathSciNet Article Google Scholar
[2] Bengfort, M., Malchow, H., Hilker, F.M.: The Fokker–Planck law of diffusion and pattern formation in heterogeneous environments. J. Math. Biol. 73, 683–704 (2016), MathSciNet Article Google Scholar
[3] Brezis, H.: Functional Analysis, Sobolev Spaces and Partial Differential Equations. Springer, New York (2011), MATH Google Scholar
[4] Bogachev, V.I., Krylov, N.V., Röckner, M., Shaposhnikov, S.V.: Fokker–Planck–Kolmogorov Equations. Amer. Math. Soc., Providence, Rhode Island (2015), Book Google Scholar
[5] Cardinali, T., Precup, R., Rubbioni, P.: A unified existence theory for evolution equations and systems under nonlocal conditions. J. Math. Anal. Appl. 432, 1039–1057 (2015), MathSciNet Article Google Scholar
[6] Cardinali, T., Precup, R., Rubbioni, P.: Heterogeneous vectorial fixed point theorems. Mediterr. J. Math. 14(2), Art. 83, 12 (2017), MathSciNet Article Google Scholar
[7] Chupin, L.: Fokker-planck equation in bounded domain. Ann. Inst. Fourier, Grenoble 60, 217–255 (2010), MathSciNet Article Google Scholar
[8] de Figueiredo, D.G.: Lectures on the Ekeland variational principle with applications and detours, Tata Institute of Fundamental Research, Bombay Lectures on Mathematics and Physics-Mathematics, vol. 81. Springer, Berlin (1989), Google Scholar
[9] Droniou, J.: Non-coercive elliptic problems. Potential Anal. 17, 181–203 (2002), MathSciNet Article Google Scholar
[10] Ekeland, I.: On the variational principle. J. Math. Anal. Appl. 47, 324–353 (1974), MathSciNet Article Google Scholar
[11] Franck, T.D.: Nonlinear Fokker-Planck Equations. Fundamentals and Applications. Springer, Berlin (2005), Google Scholar
[12] Jebelean, P., Precup, R.: Poincaré inequalities in reflexive cones. Appl. Math. Lett. 24, 359–363 (2011), MathSciNet Article Google Scholar
[13] Markowich, P.A., Villani, C.: On the trend to equilibrium for the Fokker-Planck equation: an interplay between physics and functional analysis. Matematica Contemporanea (SBM) 19, 1–31 (2000), MathSciNet MATH Google Scholar
[14] Mikhlin, S.G.: Linear Partial Differential Equations. Vysshaya Shkola, Moscow (1977). [in Russian], Google Scholar
[15] Okubo, A., Levin, S.A.: Diffusion and Ecological Problems: Modern Perspectives, Interdisciplinary Applied Mathematics, 2nd edn., vol. 14. Springer, New York (2001), Book Google Scholar
[16] Precup, R.: Nash-type equilibria and periodic solutions to nonvariational systems. Adv. Nonlinear Anal. 3, 197–207 (2014), MathSciNet MATH Google Scholar
[17] Precup, R.: Linear and Semilinear Partial Differential Equations. De Gruyter, Berlin (2013), MATH Google Scholar
[18] Precup, R.: The role of matrices that are convergent to zero in the study of semilinear operator systems. Math. Comput. Modelling 49(3-4), 703–708 (2009), MathSciNet Article Google Scholar
[19] Risken, H.: The Fokker-Planck Equation: Methods of Solution and Applications. Springer, Berlin (1996), Book Google Scholar
[20] Sattin, F.: Fick’s law and Fokker-Planck equation in inhomogeneous environments. Phys. Lett. A 372(22), 3941–3945 (2008), MathSciNet Article Google Scholar
[21] Zeidler, E.: Applied Functional Analysis: Applications to Mathematical Physics. Springer, Berlin (1995), Book Google Scholar

Paper (preprint) in HTML form

Stationary solutions of Fokker-Planck equations with nonlinear reaction terms in bounded domains

Stationary solutions of Fokker-Planck equations with nonlinear reaction terms in bounded domains

Radu Precup
Department of Mathematics, Babeş-Bolyai University, Cluj-Napoca, ROMANIA
E-mail address: r.precup@math.ubbcluj.ro
Paola Rubbioni
Department of Mathematics and Computer Science, University of Perugia, Perugia, ITALY
E-mail address: paola.rubbioni@unipg.it

Abstract: Using an operator approach, we discuss stationary solutions to Fokker-Planck equations and systems with nonlinear reaction terms. The existence of solutions is obtained by using Banach, Schauder and Schaefer fixed point theorems, and for systems by means of Perov’s fixed point theorem. Using the Ekeland variational principle, it is proved that the unique solution of the problem minimizes the energy functional, and in case of a system that it is the Nash equilibrium of the energy functionals associated to the component equations.

Key words: elliptic equation, reaction-diffusion equation, semi-linear Fokker-Planck equation, fixed point, variational method, Nash type equilibrium.

Mathematics Subject Classification: 35J60, 47H30, 47J05

1 Introduction


The Fokker-Planck equation arises as a mathematical model in many areas of physics and biology, mostly connected with the analysis of random phenomena (see, e.g., [1, 2, 4, 15, 19, 20]). It has the form

wt− div ⁢(𝐃⁢∇w+w⁢𝐅)=h,

where 𝐃=𝐃⁢(x) is a symmetric (diffusion) matrix, 𝐅=𝐅⁢(x) is a given vector field, h=h⁢(t,x) is the source term, and w=w⁢(t,x) is a probability distribution. It is the continuity equation

wt+div ⁢J=h,

for the flux density J=−𝐃⁢∇w−w⁢𝐅 involving both diffusion, by the term 𝐃⁢∇w, and drift, by w⁢𝐅. In case that 𝐃=D⁢I, where I is the identity matrix and D is a constant, the equation reads

wt−D⁢Δ⁢w−div ⁢(w⁢𝐅)=h.

We consider the semi-linear case, where the source term h is a reaction term h⁢(w) depending on the state, namely equations of the form

wt−D⁢Δ⁢w−div ⁢(w⁢𝐅)=h⁢(w).

The reaction term h⁢(w) can be very general and nonlocal both in time and space. In particular, the cases of the equations with memory and of integral-differential equations are included.  The reaction term can involve power-type nonlinearities, or rational functions simulating the saturation effect and making the equation a self-limiting model (see [2]).

In this paper, we consider only the stationary equation in a smooth bounded domain Ω⊂ℝd, with a potential flow 𝐅=−ε⁢D⁢∇H, where ε>0 and H is a sufficiently smooth function, let it be in C1⁢(Ω), having a number of properties as shown below, but which could be very irregular near the boundary. Thus we consider the following semi-linear problem

{−Δ⁢w+ε⁢div ⁢(w⁢∇H)=Φ⁢(w)in ⁢Ω,∫Ωw=1. (1.1)

Assuming that exp⁡H∈L∞⁢(Ω) and making the substitution w=(u+1)⁢ρ, with

ρ:=exp⁡(ε⁢H)/|exp⁡(ε⁢H)|L1⁢(Ω),

the above problem becomes

{−Δ⁢u−ε⁢∇u⋅∇H=Ψ⁢(u)in ⁢Ω∫Ωρ⁢u=0 (1.2)

where Ψ⁢(u)=1ρ⁢Φ⁢(ρ⁢(u+1)), leading to the study of renormalized solutions u in a special weighted space.

We also consider systems of such equations modeling the evolution of many randomly diffusing particles. This is the case of chemical reactions involving several reagents that react and diffuse simultaneously.

Compared to other approaches in the literature (see [4, 7, 9, 11, 13, 19]), our approach is essentially based on the theory of nonlinear operators and by this, the specificity of the subject is brought inside the unifying nonlinear functional analysis. We first consider the solution operator associated to the non-homogeneous problem, which is defined by using the general theory of positive-define self-adjoint linear operators, and next its composition with the nonlinear mapping giving the right-hand side (nonlinearity) of the semi-linear problem. Then, joint suitable properties of the solution operator and nonlinearity allow us to make use of several fixed point principles: Banach’s fixed point theorem, which guarantees the existence and uniqueness of the solution, and its property of being a global minimum of the energy functional; Schauder’s and Schaefer’s fixed point theorems, which not only guarantee the existence of a solution, but also give its localization in terms of the energetic norm.

Our approach to reaction-diffusion systems of Fokker-Planck equations is based on the vector method that uses matrices instead of constants, vector-valued norms and Perov’s fixed point theorem (for the vector approach to nonlinear systems, see [5, 6, 16, 18]). In this case, the obtained solution is a Nash equilibrium of the energy functionals associated to the equations of the system. The variational properties of solutions are obtained by means of Ekeland’s principle.

2 Preliminaries. Linear Fokker-Planck equations

Our approach to linear Fokker-Planck equations makes use of the variational theory of positive-define symmetric linear operators (see [14, Ch. 4], or [21, Ch. 5]). The application of this theory to linear Fokker-Planck equations is detailed in this section.

2.1 The Fokker-Planck operator

Consider the Banach normalized weighted spaces

Lρq={u:ρ1/q⁢u∈Lq⁢(Ω),∫Ωρ⁢u=0}(1≤q<+∞)

with norm

|u|Lρqq=∫Ωρ⁢|u|q.

For q=2, we endow Lρ2 with the inner product and norm

(u,v)ρ=∫Ωρ⁢u⁢v,|u|ρ=(∫Ωρ⁢u2)12.

Consider the linear operator in Lρ,2 defined by

ℒ⁢u=−Δ⁢u−ε⁢∇u⋅∇H

with the domain

D⁢(ℒ)={u∈C02⁢(Ω):∫Ωρ⁢u=0},

where C02⁢(Ω) is the space of all functions in C2⁢(Ω) with compact support included in Ω. For any u∈C02⁢(Ω), Δ⁢u∈C0⁢(Ω), and since H∈C1⁢(Ω), one has ∇H∈C⁢(Ω,ℝd). Hence ℒ⁢u∈C0⁢(Ω)⊂Lρ2, that is ℒ is well-defined. Also D⁢(ℒ) is dense in Lρ2. Indeed, if u∈Lρ2, then v:=ρ⁢u∈L2⁢(Ω) and in view of the density of C0∞⁢(Ω) into L2⁢(Ω), there exists in C0∞⁢(Ω) a sequence (vk) with vk→v in L2⁢(Ω). Let φk∈C0∞⁢(Ω) be such that φk→1 in L2⁢(Ω), and let

uk:=1ρ⁢(vk−ck⁢φk),

where ck=∫Ωρ⁢vk/∫Ωρ⁢φk, Clearly uk∈D⁢(ℒ). Also

|uk−u|ρ=|vk−ck⁢φk−v|L2⁢(Ω)≤|vk−v|L2⁢(Ω)+ck⁢|φk|L2⁢(Ω). (2.1)

Hence uk→u in Lρ2 if ck→0. To show this, first note that

0⁢<∫Ωρ≤∫Ωρ|⁢1−φk|+∫Ωρ⁢φk≤|ρ|L1⁢(Ω)⁢|1−φk|L2⁢(Ω)+∫Ωρ⁢φk

whence we have that the sequence (∫Ωρ⁢φk) is bounded from below by a positive number C. Then ck≤(1/C)⁢∫Ωρ⁢vk. Next, in view of ∫Ωρ⁢v=0, one has

|∫Ωρ⁢vk|=|∫Ωρ⁢(vk−v)|≤|ρ|L1⁢(Ω)⁢|vk−v|L2⁢(Ω)→0.

Hence ck→0 and from (2.1) it follows that uk→u in Lρ2. Therefore ℒ is densely defined on Lρ2.

The operator ℒ is symmetric. Indeed, since ∇ρ=ε⁢ρ⁢∇H, we have

(ℒ⁢u,v)ρ = −∫Ωρ⁢v⁢(Δ⁢u+ε⁢∇u⋅∇H)=∫Ω∇(ρ⁢v)⋅∇u−∫Ωv⁢∇u⋅∇ρ
= ∫Ωρ⁢∇u⋅∇v=(u,ℒ⁢v)ρ.

Finally note that

(ℒ⁢u,u)ρ=∫Ωρ⁢|∇u|2>0

for every u∈D⁢(ℒ)∖{0}, that is the operator ℒ is strictly positive.

2.2 The energetic space

We may endow D⁢(ℒ) with two inner products

⟨u,v⟩=(u,v)ρ+(ℒ⁢u,v)ρ,[u,v]=(ℒ⁢u,v)ρ

and the corresponding norms

‖u‖2=|u|ρ2+(ℒ⁢u,u)ρ,[u]2=(ℒ⁢u,u)ρ.

Let Eℒ (called the energetic space of ℒ) be the completion of the prehilbertian space (D(ℒ),∥⋅∥) and let us use the same notations ⟨⋅,⋅⟩, [⋅,⋅], ∥⋅∥ and [⋅] for the corresponding maps extended by density to Eℒ. Since |u|ρ≤‖u‖ for all u∈D⁢(ℒ), we have D⁢(ℒ)⊂Eℒ⊂Lρ2 with dense and continuous embeddings. Recall that, from the construction of the completion, any element u of Eℒ can be seen as the limit in Lρ2 of a sequence of functions from D⁢(ℒ) which is fundamental with respect the norm ∥⋅∥, and that this limit is common for all such sequences (uk),(vk) which are equivalent in the sense that ‖uk−vk‖→0. If (uk) is a fundamental sequence in D⁢(ℒ), then there exist v,vi∈L2⁢(Ω), i=1,⋅⋅,d such that

ρuk→v,ρ∂uk∂xi→vi(i=1,⋅⋅,d)in L2(Ω).

Thus, if we denote

u:=vρ,∂u∂xi:=viρ(i=1,⋅⋅,d),

then we may say that for every u,v∈Eℒ,

⟨u,v⟩ = ∫Ωρ⁢(u⁢v+∇u⋅∇v),[u,v]=∫Ωρ⁢∇u⋅∇v, (2.2)
‖u‖2 = ∫Ωρ⁢(u2+|∇u|2),[u]2=∫Ωρ⁢|∇u|2. (2.3)

Notice that the functional [⋅] is only a semi-norm on Eℒ. To make it a norm, equivalent to the norm ∥⋅∥ on Eℒ, we need a compactness assumption. To this aim, we state the following condition:

(Cq)

The embedding D⁢(ℒ)⊂Lρq is compact, i.e., any sequence of functions in D⁢(ℒ) which is bounded with respect to the norm ∥⋅∥ has a subsequence that converges in Lρq.

Clearly condition (Cq) implies that the embedding Eℒ⊂Lρq is also compact.

The next condition (H) gives an exact representation of the space Eℒ, and consequently, it is sufficient for (Cq) to hold for some values of q.

(H)

There exists a constant c>0 such that

|∇H|⁢δΩ≤cin ⁢Ω,

where δΩ gives the distance to the boundary ∂Ω, i.e.

δΩ⁢(x)=miny∈∂Ω⁡|x−y|⁢(x∈Ω).
Proposition 2.1

If condition (H) is satisfied for a constant c sufficiently small, then

Eℒ=1ρ⁢H0,ρ1⁢(Ω) (2.4)

where H0,ρ1⁢(Ω)={v∈H01⁢(Ω):∫Ωρ⁢v=0}, and [u], |ρ⁢u|H01⁢(Ω) and ‖u‖ are equivalent norms on Eℒ.

Proof. For any u∈D⁢(ℒ), one has

∇(ρ⁢u)=ρ⁢∇u+12⁢ρ⁢u⁢∇ρ=ρ⁢∇u+12⁢ρ⁢u⁢∇ρρ.

Since ρ−1⁢∇ρ=ε⁢∇H, from (H) and Hardy’s inequality [3], we can estimate the last addendum of the previous identity as follows

|12⁢ρ⁢u⁢ε⁢∇H|L2⁢(Ω)≤ε⁢c2⁢|ρ⁢uδΩ|L2⁢(Ω)≤ε⁢c~2⁢|∇(ρ⁢u)|L2⁢(Ω).

Consequently

(∫Ω|∇(ρ⁢u)|2)12≤(∫Ωρ⁢|∇u|2)12+ε⁢c~2⁢(∫Ω|∇(ρ⁢u)|2)12,
(∫Ωρ⁢|∇u|2)12≤(∫Ω|∇(ρ⁢u)|2)12+ε⁢c~2⁢(∫Ω|∇(ρ⁢u)|2)12,

whence, if c>0 is small enough that 1−ε⁢c~/2>0, we obtain

|ρ⁢u|H01⁢(Ω)≤c0⁢[u],[u]≤c1⁢|ρ⁢u|H01⁢(Ω), (2.5)

where c0=1/(1−ε⁢c~/2) and c1=1+ε⁢c~/2. As a result, for any sequence of functions uk∈D⁢(ℒ) which is fundamental with respect to the norm ∥⋅∥, the sequence ρ⁢uk is fundamental with respect to the norm |⋅|H01⁢(Ω), and conversely, if vk is fundamental with respect to the norm |⋅|H01⁢(Ω), then the sequence uk=vk/ρ is fundamental with respect to the norm ∥⋅∥. This proves (2.4), while (2.5) together with the continuous embedding (1/ρ)⁢H0,ρ1⁢(Ω)⊂Lρ2 shows that [u], |ρ⁢u|H01⁢(Ω) and ‖u‖ are equivalent norms on Eℒ.

Remark 2.1

The behavior on the boundary ∂Ω and the regularity of the solution w=(u+1)⁢ρ with u∈Eℒ strongly depend on the behavior and regularity of H. Thus, under assumption (H), if H∈H1⁢(Ω), then w∈W1,1⁢(Ω). Indeed, from (H) one has w=ρ+u⁢ρ=ρ+ρ⁢v, where v∈H01⁢(Ω). Now, since ρ was assumed in L∞⁢(Ω),

∂xiρ=ε⁢ρ⁢∂xiH∈L2⁢(Ω),

hence ρ∈H1⁢(Ω). Furthermore,

∂xi(ρ⁢v)=ρ⁢∂xiv+ε2⁢ρ⁢v⁢∂xiH∈L1⁢(Ω).

Therefore w∈ W1,1⁢(Ω). For an exhaustive discussion of regularity of solutions we refer the reader to [4, Chapter 1].

2.3 The Poincaré inequality

Assume that condition (C2) holds. The space Eℒ being reflexive (as a Hilbert space), one deduces from a result in paper [12] that

μ:=infu∈Eℒ,u≠0[u]2|u|ρ2>0

and the infimum is reached. From this, we have the Poincaré inequality

μ⁢|u|ρ2≤[u]2for ⁢u∈Eℒ

which ensures that [⋅] is a norm on Eℒ, equivalent to the norm ∥⋅∥. Let Eℒ′ be the dual of (Eℒ,[⋅]). If we identify Lρ2 to its dual, then we have

D⁢(ℒ)⊂Eℒ⊂Lρ2⊂Eℒ′ (2.6)

where the last embedding is compact too. For f∈Eℒ′ and u∈Eℒ, let (f,u) be the value of the linear functional f at u. In case that f∈Lρ2, one has (f,u)=(f,u)ρ.

Throughout the paper we assume that condition (C2) holds.

2.4 The solution operator

Returning to the operator ℒ, for a fixed f∈Eℒ′, we define the weak solution of the stationary problem

{ℒ⁢u=fin ⁢Ω∫Ωρ⁢u=0 (2.7)

as being a function u∈Eℒ such that for every v∈Eℒ, one has that [u,v]=(f,v). In particular, if f∈Lρ2, this identity becomes

∫Ωρ⁢∇u⋅∇v=∫Ωρ⁢f⁢vfor ⁢v∈Eℒ.

From Riesz’s representation theorem, since (f,⋅) is a continuous linear functional on (Eℒ,[⋅]), it follows that problem (2.7) has a unique weak solution uf. Thus we may define the solution operator

ℒ−1:Eℒ′→Eℒ,ℒ−1⁢f:=uf.

Recall that under condition (C2), the operator ℒ has a sequence of eigenvalues (λk) with 0<λ1=μ≤⋅⋅≤λk≤⋅⋅, λk→+∞, and correspondingly a sequence (ϕk) of eigenfunctions, which is orthonormal and complete in Lℒ2. Also the sequence (ϕk/λk) is orthonormal and complete in and (Eℒ,[⋅]). This yields the Fourier representation of the solution operator:

ℒ−1⁢f=∑[uf,ϕkλk]⁢ϕkλk=∑(f,ϕk)λk⁢ϕk,

where the series converges in Eℒ and Lρ2.

Also note that ℒ−1 is an isometry between Eℒ′ and Eℒ, i.e., [ℒ−1⁢f]=|f|Eℒ′ for every f∈Eℒ′, and that the exact Poincaré inequality

λ1⁢|u|ρ2≤[u]2for ⁢u∈Eℒ (2.8)

is accompanied by the Poincaré inequality for the dual, namely

λ1⁢|f|Eℒ′2≤|f|ρ2for ⁢f∈Lρ2. (2.9)

Indeed, if f∈Lρ2, then using (2.8) we have

|f|Eℒ′ = supu∈Eℒ,u≠0|(f,u)|[u]=supu∈Eℒ,u≠0|(f,u)ρ|[u]
≤ supu∈Eℒ,u≠0|f|ρ⁢|u|ρ[u]≤supu∈Eℒ,u≠01λ1⁢|f|ρ⁢[u][u]=|f|ρλ1,

that is (2.9).

2.5 The energy functional

According to the variational theory of positive-define symmetric linear operators, for each fixed f∈Eℒ′, the functional J:Eℒ→ℝ,

J⁢u=12⁢[u]2−(f,u)

is C1 and J′⁢u=ℒ⁢u−f, more exactly

(J′⁢u,v)=[u,v]−(f,v)for all ⁢u,v∈Eℒ.

Therefore the weak solution of problem (2.7) is the critical point of the energy functional J.

3 Semilinear Fokker-Planck equations

We now turn back to the semi-linear problems (1.1) and (1.2).

3.1 Existence and uniqueness via Banach’s fixed point theorem

Our first result is about the existence and uniqueness of the solution to the semilinear problem (1.2) and consequently to (1.1).

Let j0 and j be the canonical injections of the embeddings Eℒ⊂Lρ2 and Lρ2⊂Eℒ′, respectively.

Notice that problem (1.2) is equivalent with the fixed point equation u=ℒ−1⁢Ψ⁢(u) in Eℒ. In view of embeddings (2.6), we may discuss three cases:

  • •

    Ψ maps Eℒ into Eℒ′;

  • •

    Ψ maps Eℒ into Lρ2; here by the composition ℒ−1⁢Ψ we mean ℒ−1⁢j⁢Ψ;

  • •

    Ψ maps Lρ2 into Lρ2; here by ℒ−1⁢Ψ we mean ℒ−1⁢j⁢Ψ⁢j0.

Our first results are existence and uniqueness theorems, the first in terms of Ψ and the second in terms of Φ.

Theorem 3.1

Problems (1.2) and (1.1) have unique weak solution u∈ Eℒ and w=(u+1)⁢ρ∈L2⁢(Ω), respectively, if one of the following conditions holds:

(a)

Ψ:Eℒ→Eℒ′ and there is a constant 0≤a0<1 such that

|Ψ⁢(u)−Ψ⁢(v)|Eℒ′≤a0⁢[u−v]for ⁢u,v∈Eℒ.
(b)

Ψ:Eℒ→Lρ2 and there is a constant 0≤a1<λ1 such that

|Ψ⁢(u)−Ψ⁢(v)|ρ≤a1⁢[u−v]for ⁢u,v∈Eℒ.
(c)

Ψ:Lρ2→Lρ2 and there is a constant 0≤a2<λ1 such that

|Ψ⁢(u)−Ψ⁢(v)|ρ≤a2⁢|u−v|ρfor ⁢u,v∈Lρ2.

Proof. (a) Under condition (a), for any u,v∈Eℒ, one has

[ℒ−1⁢Ψ⁢u−ℒ−1⁢Ψ⁢v]=|Ψ⁢u−Ψ⁢v|Eℒ′≤a0⁢[u−v], (3.1)

and the conclusion follows from Banach’s contraction principle.

(b) In case that Ψ takes values in Lρ2, using the Poincaré inequality (2.9), we have

|Ψ⁢u−Ψ⁢v|Eℒ′≤1λ1⁢|Ψ⁢u−Ψ⁢v|ρ≤a1λ1⁢[u−v]

and the result follows from case (a) where a0=a1/λ1.

(c) This case reduces to (b) with a1=a2/λ1 since in virtue of (2.8), one has

|Ψ⁢(u)−Ψ⁢(v)|ρ≤a2⁢|u−v|ρ≤a2λ1⁢[u−v].

□

Theorem 3.2

Problems (1.2) and (1.1) have a unique weak solution u∈ Eℒ and w=(u+1)⁢ρ∈L2⁢(Ω), respectively, if one of the following conditions holds:

(d)

Φ:L2⁢(Ω)→L2⁢(Ω),Φ⁢(ρ)∈Lρ−12⁢(Ω), and there is a constant 0≤a<λ1 such that

|Φ⁢(u)⁢(x)−Φ⁢(v)⁢(x)|≤a⁢|u⁢(x)−v⁢(x)|for ⁢u,v∈L2⁢(Ω)⁢ and a.a. ⁢x∈Ω. (3.2)
(e)

Φ⁢(u)⁢(x)=f⁢(x,u⁢(x)), where f:Ω×ℝ→ℝ satisfies the Carathéodory conditions, f⁢(⋅,0)∈L2⁢(Ω), f⁢(⋅,ρ⁢(⋅))∈Lρ−12, and there exists 0≤a<λ1 such that

|f⁢(x,u)−f⁢(x,v)|≤a⁢|u−v|for all ⁢u,v∈ℝ⁢ and a.a. ⁢x∈Ω. (3.3)

Proof. (d) First, if u∈Lρ2, then ρ⁢u∈L2⁢(Ω) and

|ρ⁢Ψ⁢(u)|L2⁢(Ω) = |1ρ⁢Φ⁢((u+1)⁢ρ)|L2⁢(Ω)
≤ |1ρ⁢(Φ⁢((u+1)⁢ρ)−Φ⁢(ρ))|L2⁢(Ω)+|1ρ⁢Φ⁢(ρ)|L2⁢(Ω)
≤ a⁢|ρ⁢u|L2⁢(Ω)+|1ρ⁢Φ⁢(ρ)|L2⁢(Ω)=a⁢|u|ρ+|1ρ⁢Φ⁢(ρ)|L2⁢(Ω)<+∞.

Hence Ψ⁢(u)∈Lρ2. Also, for any u,v∈Lρ2, we have

|Ψ⁢(u)−Ψ⁢(v)|ρ = |1ρ⁢(Φ⁢((u+1)⁢ρ)−Φ⁢((v+1)⁢ρ))|ρ
= |1ρ⁢(Φ⁢((u+1)⁢ρ)−Φ⁢((v+1)⁢ρ))|L2⁢(Ω)
≤ a⁢|ρ⁢(u−v)|L2⁢(Ω)=a⁢|u−v|ρ.

Thus we are in case (c) of Theorem 3.1.

(e) Under the assumptions of f, the Nemytskii operator Φ maps L2⁢(Ω) into itself. In addition (3.3) immediately yields (3.2). Hence we are in case (d). □

3.2 Variational characterization of the solution

The next result gives a variational characterization of the solution guaranteed by the previous theorems.

Theorem 3.3

(j) Assume that Ψ is in case (a) of Theorem 3.1 and there is a C1-functional Θ:Eℒ→ℝ bounded from above on bounded sets and such that Ψ=Θ′. Then the unique solution u of problem (1.2) is the unique minimum point of the energy functional

J⁢v=12⁢[v]2−Θ⁢(v). (3.4)

(jj) Assume that Φ is in case (d) of Theorem 3.2 and there is a C1-functional θ:L2⁢(Ω)→ℝ bounded from above on bounded sets and such that Φ⁢(v)=ρ⁢θ′⁢(v−ρ) for v∈L2⁢(Ω). Then the unique solution u of problem (1.2) is the unique minimum point of the energy functional (3.4) for Θ⁢(v)=θ⁢(ρ⁢v), v∈Eℒ.

(jjj) Assume that f is in case (e) of Theorem 3.2 and in addition that there exists h∈L2⁢(Ω) and a small enough c≥0 such that

|f⁢(x,s)|≤ρ⁢(x)⁢(c⁢|s|+h⁢(x))for ⁢s∈ℝ⁢ and a.a. ⁢x∈Ω. (3.5)

Then the unique solution u of problem (1.2) is the unique minimum point of the energy functional (3.4) for

Θ⁢(v)=∫Ω(∫0ρ⁢v1ρ⁢f⁢(x,τ+ρ)⁢𝑑τ).

Proof. (j) Since J′=ℒ−Ψ, one has that the unique solution of (1.2) is the unique critical point of J. Let B be a closed ball of the space Eℒ with center at the origin and positive radius R≥|Ψ⁢(0)|Eℒ′/(1−a0). Then, by (3.1) we immediately see that

ℒ−1⁢Ψ⁢(B)⊂B. (3.6)

To prove that the solution u minimizes J we use the weak form of Ekeland’s variational principle [8, 10]. Note that the boundedness of Θ guarantees the functional J to be bounded from below on B. Indeed, for any v∈B, one has J⁢v≥−Θ⁢(v)≥c>−∞, where Θ⁢(v)≤c for every v∈B. Then, using the weak form of Ekeland’s variational principle, there is a minimizing sequence (uk) of elements in B such that

J⁢uk≤infBJ+1k (3.7)

and

J⁢uk≤J⁢v+1k⁢[v−vk] (3.8)

for all v∈B. For any fixed index k, choose

vt=uk−t⁢ℒ−1⁢J′⁢uk, 0<t<1.

Using ℒ−1⁢Ψ⁢uk=uk−ℒ−1⁢J′⁢uk, one has

vt=(1−t)⁢uk+t⁢ℒ−1⁢Ψ⁢uk.

Here, one has uk∈B, and by the invariance property (3.6), ℒ−1⁢Ψ⁢uk also belongs to B. Since B is convex, it follows that vt∈B for every t∈(0,1). Replacing v by vt into (3.15) and then dividing by t, yields

t−1⁢(J⁢uk−J⁢(uk−t⁢ℒ−1⁢J′⁢uk))≤1k⁢[ℒ−1⁢J′⁢uk],

whence letting t go to zero, one finds

(J′⁢uk,ℒ−1⁢J′⁢uk)≤1k⁢[ℒ−1⁢J′⁢uk],

that is

[ℒ−1⁢J′⁢uk]≤1k.

Hence

J⁢uk→infBJand ⁢ℒ−1⁢J′⁢uk→0in ⁢Eℒ. (3.9)

Let vk:=ℒ−1⁢J′⁢uk. From

[uk+p−uk]≤[vk+p−vk]+[ℒ−1⁢Ψ⁢uk+p−ℒ−1⁢Ψ⁢uk]

and the contraction condition (3.1), we obtain

[uk+p−uk]≤(1−a0)−1⁢[vk+p−vk].

The sequence (vk) is a Cauchy sequence (as a convergent sequence). Consequently (uk) is a Cauchy sequence, hence convergent to some u∈B. Now passing to the limit in (3.9) gives

J⁢u=infBJandJ′⁢u=0.

Do to the fact that R was taken arbitrary bigger than |Ψ⁢(0)|Eℒ′/(1−a0) and to the uniqueness of the critical point of J, we may conclude that u minimizes J on the whole space Eℒ, that is J⁢u=infEℒJ.

(jj) For u,v∈Eℒ, we have ρ⁢u, ρ⁢v∈L2⁢(Ω) and

Θ⁢(u+t⁢v)−Θ⁢(u) = θ⁢(ρ⁢(u+t⁢v))−θ⁢(ρ⁢u)
= t⁢(θ′⁢(ρ⁢u),ρ⁢v)L2⁢(Ω)+o⁢(|t|)
= t⁢(1ρ⁢Φ⁢(ρ⁢(u+1)),ρ⁢v)L2⁢(Ω)+o⁢(|t|)
= t⁢(1ρ⁢Φ⁢(ρ⁢(u+1)),v)ρ+o⁢(|t|)
= t⁢(Ψ⁢(u),v)ρ+o⁢(|t|).

It follows that Θ′⁢(u)=Ψ⁢(u). Hence we are in case (j).

(jjj) Under the assumptions on f, the function g⁢(x,s)=ρ⁢(x)−1⁢f⁢(x,s+ρ⁢(x)) satisfies the Carathéodory conditions and the growth inequality |g⁢(x,s)|≤c⁢|s|+h~⁢(x), with h~⁢(x)=h⁢(x)+c⁢ρ⁢(x), which makes the Nemytskii operator Ng associated to g to act in L2⁢(Ω). In addition the functional

θ⁢(v)=∫Ω(∫0vg⁢(x,τ)⁢𝑑τ)

is C1 on L2⁢(Ω) and θ′⁢(v)=Ng⁢(v) for v∈L2⁢(Ω). On the other hand, Φ⁢(v)=f⁢(⋅,v). Now it is easy to check the equality Φ⁢(v)=ρ⁢θ′⁢(v−ρ). Further, the functional J is bounded from below on B (where B is like in case (j)) provided c is small enough. In fact, thanks to (3.5), after some computations including the use of Poincaré’s inequality leading a constant c1, we obtain an estimate of the type

J⁢v≥(12−c⁢c1)⁢[v]2+c2⁢[v]+c3.

This shows that J is bounded from below if c<1/(2⁢c1). Therefore Ekeland’s principle applies as well in this case. □

The next theorem gives the variational characterization of the solution of problem (1.1) and it is a direct consequence of the previous result just by making the change of variable w=(u+1)⁢ρ.

Theorem 3.4

(k) Under the assumptions of Theorem 3.1, if in addition Φ⁢(v)=ρ⁢Θ′⁢(vρ−1), where Θ:Eℒ→ℝ is a C1-functional, then the unique solution w of problem (1.1) is the unique minimum point of the energy functional

J0⁢v=12⁢[vρ]2−Θ0⁢(v), (3.10)

where Θ0⁢(v)=Θ⁢(vρ−1).

(kk) Assume that Φ is in case (d) of Theorem 3.2 and there is a C1-functional θ:L2⁢(Ω)→ℝ such that Φ⁢(v)=ρ⁢θ′⁢(v−ρ) for v∈L2⁢(Ω). Then the unique solution w of problem (1.1) is the unique minimum point of the energy functional (3.10) for Θ0⁢(v)=θ⁢(v−ρ), v∈Eℒ.

(kkk) Under the assumptions of case (e) of Theorem 3.2 the unique solution w of problem (1.1) is the unique minimum point of the functional

J0⁢(v)=∫Ω(12⁢ρ⁢|∇(vρ)|2−∫0v−ρ1ρ⁢f⁢(s,τ+ρ)⁢𝑑τ). (3.11)
Remark 3.1

In particular, in the absence of drift, that is when H=0, one has ρ=ρ0:=1/m⁢e⁢s⁢(Ω) and the functional (3.11) reduces to

J0⁢(v) = 1ρ0⁢∫Ω(12⁢|∇v|2−∫0v−ρ0f⁢(s,τ+ρ0)⁢𝑑τ)
= 1ρ0⁢∫Ω(12⁢|∇v|2−∫ρ0vf⁢(x,ξ)⁢𝑑ξ)
= 1ρ0⁢∫Ω(12⁢|∇v|2−∫0vf⁢(x,ξ)⁢𝑑ξ)+1ρ0⁢∫Ω(∫0ρ0f⁢(x,ξ)⁢𝑑ξ).

Since the last term is a constant, we may say that in this case, the solution w of problem (1.1) minimizes the functional

∫Ω(12⁢|∇v|2−∫0vf⁢(x,ξ)⁢𝑑ξ),

which is the energy functional of the Dirichlet problem for the equation −Δ⁢w=f⁢(x,w).

3.3 Existence via Schauder’s fixed point theorem

If instead of Lipschitz continuity, we only assume a linear growth condition, and we add compactness, then the existence of solutions still holds based on Schauder’s fixed point theorem.

As in the case of Section 3.1, we give results first in terms of Ψ and next in terms of Φ.

Theorem 3.5

The problems (1.2) and (1.1) have at least one weak solution u∈ Eℒ and w=(u+1)⁢ρ∈L2⁢(Ω), respectively, if one of the following conditions holds:

(a)

Ψ:Eℒ→Eℒ′ is completely continuous and there are constants 0≤a0<1 and b0≥0 such that

|Ψ⁢(u)|Eℒ′≤a0⁢[u]+b0for ⁢u∈Eℒ. (3.12)
(b)

Ψ:Eℒ→Lρ2 is continuous, and there are constants 0≤a1<λ1 and b1≥0 such that

|Ψ⁢(u)|ρ≤a1⁢[u]+b1for ⁢u∈Eℒ. (3.13)
(c)

Ψ:Lρ2→Lρ2 is continuous, and there are constants 0≤a2<λ1 and b2≥0 such that

|Ψ⁢(u)|ρ≤a2⁢|u|ρ+b2for ⁢u∈Lρ2. (3.14)

Proof. (a) The operator ℒ−1⁢Ψ is completely continuous. In addition if R≥b0/(1−a0), then for every u∈Eℒ with [u]≤R, one has

[ℒ−1⁢Ψ⁢u]=|Ψ⁢u|Eℒ′≤a0⁢[u]+b0≤a0⁢R+b0≤R.

Hence the operator ℒ−1⁢Ψ maps the closed ball of Eℒ with center at the origin and radius R into itself. The conclusion follows now from Schauder’s fixed point theorem.

(b) The condition (C2) implies that the embedding Lρ2⊂Eℒ′ is compact, i.e. the injection j:Lρ2→Eℒ′ is completely continuous. Also (3.13) shows that Ψ is bounded (maps bounded sets into bounded sets). Hence the operator j⁢Ψ is completely continuous as a composition of two bounded and continuous operators where one of them, namely j, is completely continuous. In addition, from Poincaré’s inequality and (3.13),

|Ψ⁢(u)|Eℒ′≤1λ1⁢|Ψ⁢(u)|ρ≤1λ1⁢(a1⁢[u]+b1).

Hence (3.12) holds with a0=a1/λ1 and b0=b1/λ1. Thus we are in case (a).

(c) In this case, the operator Ψ⁢j0 is continuous from Eℒ to Lρ2. In addition

|Ψ⁢(u)|ρ≤a2⁢|u|ρ+b2≤a2λ1⁢[u]+b2for ⁢u∈Eℒ.

Hence we are in case (b) with a1=a2/λ1 and b1=b2. □

Theorem 3.6

The problems (1.2) and (1.1) have at least one weak solution u∈ Eℒ and w=(u+1)⁢ρ∈L2⁢(Ω), respectively, if one of the following conditions holds:

(d)

Φ:L2⁢(Ω)→L2⁢(Ω) is continuous from Lρ−12 to Lρ−12, and there are constants 0≤a<λ1 and h∈Lρ−12 such that

|Φ⁢(u)⁢(x)|≤a⁢|u⁢(x)|+h⁢(x)for ⁢u∈L2⁢(Ω)⁢ and a.a. ⁢x∈Ω. (3.15)
(e)

Φ⁢(u)⁢(x)=f⁢(x,u⁢(x)), where f:Ω×ℝ→ℝ satisfies the Carathéodory conditions, and there exist 0≤a<λ1 and h∈Lρ−12 such that

|f⁢(x,s)|≤a⁢|s|+h⁢(x)for every ⁢s∈ℝ⁢ and a.a. ⁢x∈Ω. (3.16)

Proof. (d) The case reduces to case (c) in Theorem 3.5. As in the proof of Theorem 3.2, we can show that Ψ maps Lρ2 into itself and that (3.14) holds. Next we prove that Ψ is continuous from Lρ2 to itself. Indeed, if uk→u in Lρ2, the ρ⁢uk→ρ⁢u in L2⁢(Ω), so ρ⁢uk→ρ⁢u in Lρ−12. Consequently, Φ⁢((uk+1)⁢ρ)→Φ⁢((u+1)⁢ρ) in Lρ−12, that is Ψ⁢(uk)→Ψ⁢(u) in Lρ2, as wished.

(e) We reduce this case to (d). Under the above conditions on f, the Nemytskii operator associated to f maps continuously L2⁢(Ω) into L2⁢(Ω) and |f⁢(x,v)|L2⁢(Ω)≤a⁢|v|L2⁢(Ω)+|h|L2⁢(Ω). The same is true for the function

g⁢(x,s)=1ρ⁢(x)⁢f⁢(x,s⁢ρ⁢(x))

which also satisfies the Carathéodory conditions and the growth condition

|g⁢(x,s)|≤a⁢|s|+h⁢(x)for every ⁢s∈ℝ⁢ and a.a. ⁢x∈Ω.

Thus the Nemytskii operator associated to g is continuous from L2⁢(Ω) to L2⁢(Ω). This implies that Φ is continuous from Lρ−12 to Lρ−12. Indeed, if uk→u in Lρ−12, then

uk/ρ→u/ρin ⁢L2⁢(Ω),whence ⁢g⁢(⋅,uk/ρ)→g⁢(⋅,u/ρ)in ⁢L2⁢(Ω).

Hence f⁢(⋅,uk)/ρ→f⁢(⋅,u)/ρ in L2⁢(Ω), that is f⁢(⋅,uk)=Φ⁢(uk)→f⁢(⋅,u)=Φ⁢(u) in Lρ−12. Also (3.16) clearly gives (3.15). Hence we are in case (d). □

3.4 Existence via Schaefer’s fixed point theorem

In the next result, the linear growth of the nonlinear reaction term is relaxed in case that in compensation, a sign condition holds for a part of the reaction term.

Let Φ=Φ0+Φ1 and correspondingly Ψ=Ψ0+Ψ1 and f=f0+f1.

We first state a general existence principle.

Theorem 3.7

Let Ψ0 be as in Theorem 3.5 (a). If in addition Ψ1 :Eℒ→Eℒ′ is completely continuous and

(Ψ1⁢(u),u)≤0for ⁢u∈Eℒ,

then the problems (1.2) and (1.1) have at least one weak solution u∈ Eℒ and w=(u+1)⁢ρ∈L2⁢(Ω), respectively. In addition [u]≤b0/(1−a0).

Proof. The operator ℒ−1⁢Ψ is completely continuous and for every solution u∈Eℒ∖{0} of the equation u=λ⁢ℒ−1⁢Ψ⁢(u) and any λ∈(0,1), one has

[u] = λ⁢(Ψ⁢(u),u)[u]=λ⁢(Ψ0⁢(u),u)[u]+λ⁢(Ψ1⁢(u),u)[u]
≤ λ⁢(Ψ0⁢(u),u)[u]≤λ⁢|Ψ0⁢(u)|Eℒ′<a0⁢[u]+b0.

Here we have assumed without loss of generality that b0>0. Hence [u]<b0/(1−a0), that is, the set of all solutions of the equations u=λ⁢ℒ−1⁢Ψ⁢(u) for λ∈(0,1), is bounded in Eℒ. Now Schaefer’s fixed point theorem guarantees the existence of a fixed point u∈ Eℒ of ℒ−1⁢Ψ with [u]≤b0/(1−a0). □

The next theorem gives us some sufficient conditions for the complete continuity of Ψ1.

Theorem 3.8

The operator Ψ1 is completely continuous from Eℒ to Eℒ′ if the compactness condition (Cq) holds for some q≥2, and one of the following conditions is satisfied:

(i)

Ψ1:Eℒ→Lρp is continuous and bounded for some p≥q/(q−1);

(ii)

Ψ1:Lρq→Lρp is continuous and bounded for some p≥q/(q−1);

(iii)

Φ1:Lρ1−qq→Lρ1−pp is continuous and bounded for some p≥q/(q−1);

(iv)

Φ1⁢(u)⁢(x)=f1⁢(x,u⁢(x)), where f1:Ω×ℝ→ℝ satisfies the Carathéodory conditions, and there exist a∈ℝ+ and h∈Lρ1−pp such that

|f1⁢(x,s)|≤a⁢ρp−qp⁢|s|qp+h⁢(x)for every ⁢s∈ℝ⁢ and a.a. ⁢x∈Ω. (3.17)

Proof. (i) From (Cq), the embedding Eℒ⊆Lρq is compact, and so is the embedding Lρq′⊆Eℒ′. For p≥q/(q−1)=q′, since ρ∈L∞⁢(Ω), one has Lρp⊆Lρq′, whence the compact inclusion Lρp⊆Eℒ′.

(ii) Use the embedding Eℒ⊆Lρq to reduce the case to (i).

(iii) Let uk→u in Lρq. Then ρ⁢uk→ρ⁢u in Lρ1−qq. Then, ρ⁢(uk+1)→ρ⁢(u+1) in Lρ1−qq, and by the assumption, Φ1⁢(ρ⁢(uk+1))→Φ1⁢(ρ⁢(u+1)) in Lρ1−pp. This yields Ψ1⁢(uk)=ρ−1⁢Φ1⁢(ρ⁢(uk+1))→Ψ1⁢(u) in Lρp. Thus we are in case (ii).

(iv) We are in case (iii). First observe that Φ1⁢(u)⁢(x)=ρ⁢(x)1−1/p⁢Ng⁢(v)⁢(x), where for each u∈Lρ1−qq, v:=ρ1/q−1⁢u∈Lq⁢(Ω), and

g⁢(x,s)=ρ⁢(x)1p−1⁢f1⁢(x,ρ⁢(x)1−1q⁢s).

The problem reduces to show that the Nemytskii operator Ng is well-defined from Lq⁢(Ω) to Lp⁢(Ω). Indeed, using (3.17) we have

|g⁢(x,s)|=ρ⁢(x)1p−1⁢|f1⁢(x,ρ⁢(x)1−1q⁢s)|≤a⁢|s|qp+h0⁢(x),

where h0=ρ1−pp⁢h∈Lp⁢(Ω).

Remark 3.2 (Positive solutions)

As it is well-known, the existence of nonnegative solutions of boundary values problems is closely connected with maximum principles (see [4, Chapter 2]). For our elliptic operator ℒ⁢u=−Δ⁢u−ε⁢∇u⋅∇H, the maximum principle holds, more exactly, if u∈C2⁢(Ω),ℒ⁢u≥0 and there is x0∈Ω with u⁢(x0)=infΩu, then u is constant on the connected component of Ω that contains x0. Consequently, assuming that Φ⁢(v)≥0 for every function v, and that a solution u of problem (1.2) is regular belonging to D⁢(ℒ), then u≥−1.Indeed, otherwise, we would have u=0 around ∂Ω, u≠0 on a compact subset of Ω and u⁢(x0)=minΩ⁡u<−1 for some x0∈Ω. This implies that u is constant u⁢(x0) on the connected component of Ω that contains x0. But this is impossible since u is zero on ∂Ω. The case of generalized solutions can be discussed similarly using the maximum principle for weak solutions.

4 Fokker-Planck reaction-diffusion systems

Under the assumptions of Theorem 3.3, the unique stationary solution of a single equation is a global minimum of the associated energy function. We now prove that in case of a system of equations and under suitable conditions, the stationary solution is a Nash type equilibrium with respect to the couple of energy functionals associated to the component equations.

For simplicity, we shall consider only systems of two equations, that is

{−Δ⁢w1+ε1⁢d⁢i⁢v⁢(w1⁢∇H1)=Φ1⁢(w1,w2)−Δ⁢w2+ε2⁢d⁢i⁢v⁢(w2⁢∇H2)=Φ2⁢(w1,w2)∫Ωw1=1,∫Ωw2=1. (4.1)

In this case, denoting by

ρi=eεi⁢Hi/|eεi⁢Hi|L1⁢(Ω)

and making the substitution

wi=(ui+1)⁢ρi,

for i=1,2, we arrive to the system

{−Δ⁢u1−ε1⁢∇u1⋅∇H1=Ψ1⁢(u1,u2)−Δ⁢u2−ε2⁢∇u2⋅∇H2=Ψ2⁢(u1,u2)∫Ωρ1⁢u1=0,∫Ωρ2⁢u2=0 (4.2)

where

Ψi⁢(u1,u2)=1ρi⁢Φi⁢((w1+1)⁢ρ1,(w2+1)⁢ρ2),i=1,2.

All the elements defined in Section 2 for one equation, ℒ, Eℒ, J and the scalar products and norms given by (2.2), are now duplicated for the two equations of the system, and we show it by an index i, i=1,2. We point out that we do not assume a variational structure on the whole system, but only for each component equation; thus, more exactly, we assume that there exist functions Θi:Eℒ1×Eℒ2→ℝ, i=1,2, bounded from above on bounded sets and such that for each i, Ψi is the Fréchet derivative of Θi⁢(u1,u2) with respect to the variable ui. Hence, the energy functionals are

Ji⁢(u1,u2)=12⁢[ui]i2−Θi⁢(u1,u2),i=1,2.

The analogue for systems of Theorem 3.1 case (a) is the following result. The reader can easily obtain the analogues for the cases (b)-(e).

Theorem 4.1

Assume that Ψi:Eℒ1×Eℒ2→Eℒ1′×Eℒ2′ and there are nonnegative constants ai⁢j, i,j=1,2, such that

|Ψi⁢(u1,u2)−Ψi⁢(v1,v2)|Eℒi′≤ai⁢1⁢[u1−v1]1+ai⁢2⁢[u2−v2]2,i=1,2, (4.3)

for all u1,v1∈Eℒ1 and u2,v2∈Eℒ2, and the spectral radius of the matrix M=(ai⁢j)i,j=1,2 is strictly less than 1. Then, problems (4.2) and (4.1) have unique weak solutions (u1,u2)∈ Eℒ1×Eℒ2 and (w1,w2)∈L2⁢(Ω)×L2⁢(Ω), wi=(ui+1)⁢ρi (i=1,2), respectively, and (u1,u2) is a Nash equilibrium of the pair of functionals (J1,J2), namely

J1⁢(u1,u2)=minEℒ1⁡J1⁢(⋅,u2),J2⁢(u1,u2)=minEℒ2⁡J2⁢(u1,⋅). (4.4)

Proof. First we prove the existence and the uniqueness of the solution to problem (4.2). Using (4.3) and the same arguments as in the proof of Theorem 3.1(a), we arrive to

[ℒi−1⁢Ψi⁢(u1,u2)−ℒi−1⁢Ψi⁢(v1,v2)]i≤ai⁢1⁢[u1−v1]1+ai⁢2⁢[u2−v2]2,i=1,2, (4.5)

which, using the matrix M, can be written in the matrix form

([ℒ1−1⁢Ψ1⁢(u1,u2)−ℒ1−1⁢Ψ1⁢(v1,v2)]1[ℒ2−1⁢Ψ2⁢(u1,u2)−ℒ2−1⁢Ψ2⁢(v1,v2)]2)≤M⁢([u1−v1]1[u2−v2]2).

Now, since the spectral radius of matrix M is strictly less than 1, the existence and uniqueness of the solution (u1,u2) follow from Perov’s fixed point theorem (see [18, Theorem 1]).

In order to use the weak form of Ekeland’s principle, we look for two balls Bi⊂Eℒi of positive radius Ri, i=1,2, with the property that

ℒi−1⁢Ψi⁢(B1×B2)⊂Bifor ⁢i=1,2. (4.6)

Taking v1=v2=0 in (4.5) and assuming that [ui]i≤Ri, i=1,2, we obtain

[ℒi−1⁢Ψi⁢(u1,u2)]i≤γi+ai⁢1⁢R1+ai⁢2⁢R2,i=1,2,

where γi=[ℒi−1⁢Ψi⁢(0,0)]i. Hence, in order to obtain the desired inclusions (4.6), it is enough to have γi+ai⁢1⁢R1+ai⁢2⁢R2≤Ri, i=1,2, or in the matrix form

(I−M)⁢(R1R2)≥(γ1γ2).

Multiplying on the left by (I−M)−1 (which is a positive matrix since the spectral radius of M is less than 1, see [18]) yields

(R1R2)≥(I−M)−1⁢(γ1γ2),

which shows that the desired numbers R1,R2 exist.

Next, we prove that the solution (u1,u2) is the Nash equilibrium of the pair of functionals (J1,J2). To this aim, we use an iterative procedure. Denote by Ji⁢i the Fréchet derivative of the functional Ji⁢(u1,u2) with respect to ui. To begin the iterative procedure, we fix an arbitrary element u2,0∈B2. At each step k≥1, u2,k−1∈B2 been found at the previous step k−1, first as in the proof of Theorem 3.3, we apply Ekeland’s principle in B1 to the functional J1⁢(⋅,u2,k−1) and find an element u1,k∈B1 such that

J1⁢(u1,k,u2,k−1)≤infB1J1⁢(⋅,u2,k−1)+1k,[ℒ1−1⁢J11⁢(u1,k,u2,k−1)]1≤1k. (4.7)

Next, we apply Ekeland’s principle in B2 to the functional J2⁢(u1,k,⋅) and obtain an element u2,k∈B2 with

J2⁢(u1,k,u2,k)≤infB2J2⁢(u1,k,⋅)+1k,[ℒ2−1⁢J22⁢(u1,k,u2,k)]2≤1k. (4.8)

Our aim now is to prove that the two sequences (ui,k)k, i=1,2 are Cauchy and so convergent. Let

αk=ℒ1−1⁢J11⁢(u1,k,u2,k−1)and ⁢βk=ℒ2−1⁢J22⁢(u1,k,u2,k).

Clearly αk→0 and βk→0 in Eℒ1 and Eℒ2, respectively. As in the case of only one equation, we have

u1,k−ℒ1−1⁢Ψ1⁢(u1,k,u2,k−1) = αk, (4.9)
u2,k−ℒ2−1⁢Ψ2⁢(u1,k,u2,k) = βk. (4.10)

By equation (4.9), we deduce

[u1,k+p−u1,k]1 ≤ [ℒ1−1⁢Ψ1⁢(u1,k+p,u2,k+p−1)−ℒ1−1⁢Ψ1⁢(u1,k,u2,k−1)]1+[αk+p−αk]1
≤ a11⁢[u1,k+p−u1,k]1+a12⁢[u2,k+p−1−u2,k−1]2+[αk+p−αk]1
≤ a11⁢[u1,k+p−u1,k]1+a12⁢[u2,k+p−u2,k]2
+a12⁢([u2,k+p−1−u2,k−1]2−[u2,k+p−u2,k]2)+[αk+p−αk]1.

By equation (4.10) we have

[u2,k+p−u2,k]2≤a21⁢[u1,k+p−u1,k]1+a22⁢[u2,k+p−u2,k]2+[βk+p−βk]2. (4.12)

Denote

δk,p=[u1,k+p−u1,k]1,ηk,p=[u2,k+p−u2,k]2,ξk,p=[αk+p−αk]1,χk,p=[βk+p−βk]2.

Obviously, ξk,p→0 and χk,p→0 as k→∞, uniformly with respect to p. Using the above notations, the inequalities (4) and (4.12) become

δk,p ≤ a11⁢δk,p+a12⁢ηk,p+ξk,p+a12⁢(ηk−1,p−ηk,p),
ηk,p ≤ a21⁢δk,p+a22⁢ηk,p+χk,p.

These can be put under the following matrix form

(δk,pηk,p)≤M⁢(δk,pηk,p)+(ξk,p+a12⁢(ηk−1,p−ηk,p)χk,p).

Hence

(δk,pηk,p)≤(I−M)−1⁢(ξk,p+a12⁢(ηk−1,p−ηk,p)χk,p).

Let (I−M)−1=(μi⁢j). Then

δk,p ≤ μ11⁢(ξk,p+a12⁢(ηk−1,p−ηk,p))+μ12⁢χk,p,
ηk,p ≤ μ21⁢(ξk,p+a12⁢(ηk−1,p−ηk,p))+μ22⁢χk,p.

The second inequality yields

ηk,p≤μ21⁢a121+μ21⁢a12⁢ηk−1,p+μ21⁢ξk,p+μ22⁢χk,p1+μ21⁢a12. (4.13)

Note that the sequence (uk,p)k is bounded uniformly with respect to p as follows from (4.10) (recall that u1,k∈B1 and u2,k∈B2). Consequently, by its definition, the sequence (ηk,p)k is also bounded uniformly with respect to p. Thus we can apply to it the following lemma proved in [16, Lemma 3.2].

Lemma 4.1

Let (xk,p)k and (yk,p)k be two sequences of real numbers depending on a parameter p such that the sequence (xk,p)k is bounded uniformly with respect to p, and

0≤xk,p≤λ⁢xk−1,p+yk,p (4.14)

for all k,p and some λ∈[0,1[. If yk,p→0 as k→∞ uniformly with respect to p, then xk,p→0 uniformly with respect to p.

Indeed, (4.13) reads as (4.14). Therefore ηk,p→0 uniformly with respect to p, which proves that the sequence (u2,k)k is Cauchy. Next, inequality (4) together with a11<1 (which is a consequence of the fact that the spectral radius of matrix M is less that 1) implies that (u1,k)k is Cauchy too. Let v1,v2 be the limits of the sequences (u1,k)k and (u2,k)k, respectively. Passing to the limit in (4.7) and (4.8), we obtain that (v1,v2) solves system (4.2). The uniqueness of the solution and the arbitrariness of R1,R2 imply that (v1,v2)=(u1,u2) and that (u1,u2) satisfies (4.4). □

Acknowledgements

The research is carried out within the national group GNAMPA of INdAM and supported by the INdAM-GNAMPA Project 2019 Metodi topologici per problemi differenziali nonlineari ed applicazioni.

References

  • [1] V. Barbu, Generalized solutions to nonlinear Fokker-Planck equations, J. Differential Equations 261 (2016), 2446–2471.
  • [2] M. Bengfort, H. Malchow and F.M. Hilker, The Fokker–Planck law of diffusion and pattern formation in heterogeneous environments, J. Math. Biol. 73 (2016), 683–704.
  • [3] H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer, New York, 2011.
  • [4] V.I. Bogachev, N.V. Krylov, M. Röckner and S.V. Shaposhnikov, Fokker–Planck–Kolmogorov Equations, Amer. Math. Soc., Providence, Rhode Island, 2015.
  • [5] T. Cardinali, R. Precup and P. Rubbioni, A unified existence theory for evolution equations and systems under nonlocal conditions, J. Math. Anal. Appl. 432 (2015), 1039–1057.
  • [6] T. Cardinali, R. Precup and P. Rubbioni,  Heterogeneous vectorial fixed point theorems, Mediterr. J. Math. 14 (2017), no. 2, Art. 83, 12 pp.
  • [7] L. Chupin, Fokker-Planck equation in bounded domain, Ann. Inst. Fourier, Grenoble 60 (2010), 217–255.
  • [8] D.G. de Figueiredo, Lectures on the Ekeland variational principle with applications and detours, Tata Institute of Fundamental Research, Bombay, Lectures on Mathematics and Physics-Mathematics 81, Springer, Berlin, 1989.
  • [9] J. Droniou, Non-coercive elliptic problems, Potential Anal. 17 (2002), 181–203.
  • [10] I. Ekeland, On the variational principle, J. Math. Anal. Appl. 47 (1974), 324–353.
  • [11] T.D. Franck, Nonlinear Fokker-Planck Equations. Fundamentals and Applications, Springer, Berlin, 2005.
  • [12] P. Jebelean and R. Precup, Poincaré inequalities in reflexive cones, Appl. Math. Lett. 24 (2011), 359–363.
  • [13] P.A. Markowich and C. Villani, On the trend to equilibrium for the Fokker-Planck equation: an interplay between physics and functional analysis, Matematica Contemporanea (SBM) 19 (2000), 1–31.
  • [14] S.G. Mikhlin, Linear Partial Differential Equations, Vysshaya Shkola, Moscow, 1977 [in Russian].
  • [15] A. Okubo and S.A. Levin, Diffusion and Ecological Problems: Modern Perspectives, Interdisciplinary Applied Mathematics, vol 14, 2nd edn., Springer, New York, 2001.
  • [16] R. Precup, Nash-type equilibria and periodic solutions to nonvariational systems, Adv. Nonlinear Anal. 3 (2014), 197–207.
  • [17] R. Precup, Linear and Semilinear Partial Differential Equations, De Gruyter, Berlin, 2013.
  • [18] R. Precup, The role of matrices that are convergent to zero in the study of semilinear operator systems, Math. Comput. Modelling 49 (2009), no. 3-4, 703–708.
  • [19] H. Risken, The Fokker-Planck Equation: Methods of Solution and Applications, Springer, Berlin, 1996.
  • [20] F. Sattin, Fick’s law and Fokker-Planck equation in inhomogeneous environments, Physics Letters A 372 (2008), 22, 3941–3945.
  • [21] E. Zeidler, Applied Functional Analysis: Applications to Mathematical Physics, Springer, Berlin, 1995.
2022

Related Posts