Home Localization and multiplicity in the homogenization of nonlinear problems
Article Open Access

Localization and multiplicity in the homogenization of nonlinear problems

  • Renata Bunoiu and Radu Precup EMAIL logo
Published/Copyright: May 22, 2019

Abstract

We propose a method for the localization of solutions for a class of nonlinear problems arising in the homogenization theory. The method combines concepts and results from the linear theory of PDEs, linear periodic homogenization theory, and nonlinear functional analysis. Particularly, we use the Moser-Harnack inequality, arguments of fixed point theory and Ekeland's variational principle. A significant gain in the homogenization theory of nonlinear problems is that our method makes possible the emergence of finitely or infinitely many solutions.

MSC 2010: 35B27; 35J25

1 Introduction

The aim of homogenization theory is that, starting from an initial problem depending on a small positive parameter ε related to the heterogeneous character of the problem, to find an approximation of its solution which solves a problem independent of ε, called the homogenized problem. From a mathematical point of view this is obtained by passing to the limit with ε ⟶ 0 in the initial problem.

One of the main features in classical homogenization problems is the uniqueness of the solution of the problem under study. This is crucial in order to justify the convergence of the whole sequence of solutions of the microscopic problem to the solution of the homogenized problem. In the nonlinear case, there is a wide class of problems for which the uniqueness of the solution fails. When several solutions of a problem exist, it is important to be able to carry on the analysis on one particular solution and this needs its localization. A lot of results for nonlinear problems in homogenization theory have been obtained in the last years (see for instance [1], [3], [7], [8] [10], [12], [13], [18]). In most cases, in order to pass to the limit, it was necessary to obtain upper estimates of the solutions. However, up to our knowledge, none of these works had to deal with estimates from below of the solutions.

The aim of this paper is to develop a method of localization in bounded sets for the solutions of nonlinear homogenization problems, which allows the emergence of multiple solutions for a given problem. The localization which will be given on the initial problem, for each ε small enough, will be preserved after passage to the limit with ε ⟶ 0. We will illustrate this method on a simple model, namely the homogenization of an elliptic nonlinear Dirichlet problem stated in a fixed domain, with periodic fast oscillating coefficients depending on the small parameter ε, which expresses the heterogeneous properties of the medium. We emphasize that our method applies to more involved problems in nonlinear homogenization theory.

The results are then extended to systems, for which we are able to localize each component of the solution, independently. Moreover, the solution appears as a Nash equilibrium with respect to the energy functionals associated to the equations of the system. In this case again, for a given nonlinear system, by the localization technique, we are able to prove the existence of finitely or infinitely many local Nash equilibrium solutions.

The starting point of our analysis is a linear microscopic pivot problem with a right-hand side, gϵL(Ω) which allows to estimate the solutions uε both from below and from above, independently of ε, via two constants mg and M g respectively. At this point, we underline the role of the Moser-Harnack inequality in establishing the lower bound mg (see Lemma 3.1), which, up to our knowledge, is new in the frame of homogenization theory.

The question we answer in this paper is how the source term g can be modified by a nonlinear source factor f(uε) with a feedback action, such that the solution uε of the new problem remains L-bounded independently of ε, between two given bounds r, R. The same question is also addressed for nonlinear systems.

Our study is motivated by real-world applications in physics, engineering and biology, from where at least two physical requirements emerge. In these cases the solution corresponds to a density rate.

First requirement: Find a suitable nonlinear state-dependent source factor f(uε), in order to guarantee that the density rate uε stays bounded between two a priori given bounds r, R, 0 < r < R. Such a requirement is natural in problems arising from medicine, where for instance the concentration rate of a drug in a tissue with a nonhomogeneous structure must remain between some prescribed limits. The same question can be addressed in problems of thermodynamics, when the heat diffusion in strongly heterogeneous media has to be limited.

Second requirement: The nonlinear state-dependent source factor f being given, find the bounds of the corresponding density rate uε.

Our approach makes use of methods and results from the linear theory of PDEs (weak and strong maximum principle, Poincaré and Moser-Harnack inequalities, compact embedding theorems), linear periodic homogenization theory, and concepts and techniques of nonlinear functional analysis (Banach's contraction principle and its vector version owed to Perov, Ekeland variational principle, Nash equilibrium).

The paper is organized as follows: In Section 2 we give the statement of the problem and we recall some basic results of periodic linear homogenization theory. In Section 3 we state and prove the main existence and localization result giving a variational characterization to the solution and, as a consequence, we obtain a multiplicity result for nonlinearities with a repeated suitable behavior. These solutions are local minima of the associated energy functionals. Section 4 is devoted to the study of a nonlinear system, for which we are able to localize componentwise the solution.

2 Statement of the problem and preliminaries results

We are concerned with the following nonlinear problem with homogeneous Dirichlet boundary conditions, and strongly periodic oscillating coefficients:

(2.1) d i v A ε ( x ) u ε = g ( x ) f u ε in  Ω u ε = 0 on  Ω .

Here Ωn is an open domain with sufficiently smooth boundary Ω. We make the following assumptions on the data:

Let Y = (0, 1) n. For α , β ϵ R , with 0 α β , let M α , β , Y be the set of all the matrices Aϵ(L(Y))n×n such that for any ξϵn,

(2.2) α| ξ |2(A(y)ξ,ξ)β| ξ |2,

almost everywhere in Y. For a Y-periodic symmetric matrix A ∈ M(α, β, Y), we set

Aε(χ)=A({ χε })a.e.inΩ.

The function g belongs to the space L(Ω) and the function f: is continuous.

Throughout this paper Ω is a fixed subdomain of Ω,ΩΩ. We first look for sufficient conditions on f and g such that, for given numbers 0 < r < R, problem (2.1) has a solution uε such that the following boundedness conditions

(2.3) 0uε(χ)Rfora.a.xϵΩ,
(2.4) ruε(χ)fora.a.xϵΩ,

are satisfied for each ε>0 small enough.

Then, we prove that, after passage to the limit in (2.1) with ε0, the limit u of the sequence uε solves uniquely the problem

(2.5) d i v A 0 u = g ( x ) f u in  Ω u = 0 on  Ω ,

which is therefore the homogenized problem of (2.1). The solution u of the homogenized problem, satisfies the same boundedness conditions. Here the entries of the constant matrix A0 are defined in (2.8)-(2.9) below. This localization is useful for the numerical computation of the solutions, especially when an initial approximation is required.

In addition, this localization result, applied to distinct pairs of numbers (rk, R k) yields finitely or infinitely many solutions for problem (2.1), provided that either

0<rk<Rk<rk+1<Rk+1

or

0<rk+1<Rk+1<rk<Rk.

Analogously, in Section 4, we discuss the localization of solutions for nonlinear systems. In this case the localization can be obtained componentwise, which is a significant gain for the treatment of systems in comparison with the classical approach.

2.1 The linear homogenization problem

We start by recalling a classical result from the homogenization of linear elliptic problems. With the domain Ω and the coefficients Aε defined as before, we consider the linear problem (see [2, Chapter 1]):

(2.6) d i v A ε ( x ) u ε = h ( x ) in  Ω u ε = 0 on  Ω .

If hϵL2(Ω) is given, then one can prove (see for instance [5, Chapter 6] and references therein) that problem (2.6) admits a unique solution uεϵH01(Ω) which converges to some u weakly in H01(Ω), and strongly in L2 (Ω). The limit u is the unique solution of the homogenized problem

(2.7) { div(A0u)=h(x)inΩu=0onΩ

with the constant homogenized coefficients given by

(2.8) A i j 0 = Y a i j k = 1 n a i k χ j y k d y ,

where, for i, j = 1,..., n, aij are the entries of the matrix A. The j-th components of the vectorial function , χjϵHper1(Y)(j=1,,n) are the unique weak solutions of the following cell problems:

(2.9) d i v y ( A ( y ) ( y χ j e j ) ) = 0 in  Y , M Y ( χ j ) = 0.

Here Hper1(Y) denotes the space of H1 functions which are Y-periodic and My(v) represents the mean value on the set Y of a function v.

When we make additional assumptions on the right-hand side h, one has more information on the corresponding solution of problem (2.6), as illustrated by the two remarks below.

Remark 2.1

If h(x)0fora.e.xϵΩ, then the weak maximum principle (see [9, Theorem 8.11) guarantees that for every ε > 0, one has uε(x) 0 a.e. in Ω.

Remark 2.2

hϵL(Ω) then, according to [9, Theorem 8.15], the solution uε belongs to L (Ω) and there is a constant Γh, independent of ε, such that

uε L(Ω)Γh.

In view of Remark 2.2, we define for each ε ∈ (0, 1), the solution operator

Sε:L(Ω)V:=H01(Ω)L(Ω),hSε(h):=uhε,

where uhε is the corresponding solution of problem (2.6).

The operator Sε is linear and due to Remark 2.1, it preserves the positivity and the order, more exactly

if0h1h2,then0Sε(h1)Sε(h2).

Searching now a weak solution uε ∈ V of problem (2.1) is thus equivalent to solving the nonlinear equation

(2.10) uε=Sε(gf(uε)).

For each ε, denoting

(2.11) Nε(v):=Sε(gf(v)),vϵV,

any solution of problem (2.10) is a fixed point of the operator Nε. The existence and the uniqueness of such a fixed point will be obtained by using Banach's contraction principle. Alternatively, the existence and the uniqueness of the solution of (2.1) can be obtained using the variational approach in connection with the energy functional associated to (2.1), namely

Jε:V,Jε(v)=12ΩAε(χ)vvΩg(χ)0v(χ)f(τ)dτ.

Similarly, we associate to problem (2.5) the energy functional

J:V,J(v)=12ΩA0vvΩg(χ)0v(χ)f(τ)dτ.

3 Main results for equations

3.1 A lower bound lemma

Our aim in this paper being the localization of solutions of problem (2.1), as shown in (2.3)–(2.4), we need bounds from above and from below. For the upper bounds, we use classical results from the literature, as mentioned in Remark 2.2. Getting lower bounds is a more difficult task and, up to our knowledge, it is new in the frame of the homogenization theory. To obtain such lower bound estimates, we make use of an auxiliary result, Lemma 3.1, whose proof is based on Moser-Harnack's inequality. The key auxiliary result is the following.

Lemma 3.1

For each h ∈ L (Ω) with h ≥ 0 in Ω and h > 0 on a subset of nonzero measure of Ω, there exist εh, > 0 and yh > 0 such that uε (x) ≥ yh for a.a. x in the fixed subdomain ΩΩ and every ε < εh.

Proof. The proof will be done by contradiction. Assume the contrary. Then for every integer ≥ ≥ 1, there exists εk<1/k with

(3.1) infΩuεk<1k.

Let pϵ be fixed such that 1p<n/(n2)ifn3andp=2ifn<3. Using Moser-Harnack's inequality (see [11, Theorem 12.1.2]) and covering arguments as in [15, Theorem 1.3], we have

(3.2) infΩuεkC uεk LP(Ω),

where the constant C > 0 only depends on p,Ω,Ω, and the ellipticity constants α and β in (2.2). From (3.1) and (3.2) we deduce that

(3.3) uεk LP(Ω)0ask.

Since p<n/(n2)<2n/(n2)=2, the injection of H01(Ω)inLp(Ω) is compact, and consequently up to a subsequence we nave uεkuinLp(Ω), Where u is the solution of the homogenized problem (2.7). Then from (3.3), u Lp(Ω)=0, that is u = 0 in Ω which is impossible in view of the strong maximum principle applied to the homogenized problem (2.7) in which h ≥ 0 and h is not the null function.

Coming back to our initial question, namely the localization of the solution for problem (2.1), let us define, for the given numbers 0 < r < R, the set

DrR={ vV:0v(x)Ra.e. inΩ,rv(x)a.e.inΩ }.

Note that the set DrR is a closed subset of H01(Ω). Indeed, if vkDrRandvkvinH01(Ω), then vkv in L2(Ω) and there is a subsequence (vkj) of the sequence (vk) such that vki(x)v(x) for a.a. x ∈ Ω (see [4, Theorem 4.9]). Then, from the definition of DrR, since

0vkj(x)Ra.einΩ,andrvkj(x)a.e.inΩ,

passing to the limit we obtain

0v(x)Ra.einΩ,andrv(x)a.e.inΩ,

that is vDrR.

3.2 Invariance of DrR through the operator Nε

The first step for applying Banach's contraction principle to the nonlinear operator (2.11) is to guarantee the invariance condition. To this aim we state our first hypothesis.

(h1) g(x)0 a.e in Ω,g(x)=0forxΩ/Ω, and g > 0 on a subset of nonzero measure; f ( τ ) ≥0 for every τ0.

Let vDrR. Since rv(x)Ra.e.inΩ and f is continuous, we have

mf:=minτ[ r,R ]f(τ)f(v(x))Mf:=maxτ[ r,R ]f(τ),a.e.inΩ.

Clearly 0mfMf. This, together with the property of g of being nonnegative in Ω and to be zero outside Ω, gives

g(x)mfg(x)f(v(x))g(x)Mf,a.e.inΩ.

Then using the positivity and monotonicity properties of Sε, Remark 2.2, by applying Nε to the previous inequalities we obtain

0Nε(v)(x)=Sε(gf(v))(x)Sε(gMf)(x)=MfSε(g)(x)MfΓg,fora.a.xΩ.

For the lower estimation, we make use of Lemma 3.1, applied for g, and we obtain:

(3.4) Nε(v)(x)Sε(gmf)(x)=mfSε(g)(x)mfyg,fora.a.xΩ.

We are now in position to derive a sufficient condition for Nε (v) ∈ DrR to hold, and thus for the invariance condition Nε(DrR)DrR to be verified, namely

(h2) MfΓgRandrmfyg.

Remark 3.1

If in addition the function f is nondecreasing on the interval [r, R] , then mf = f (r), M f = f (R) and condition (h2) becomes

(h2)' f(R)ΓgRandrf(r)yg.

This gives an answer to the second requirement stated in the Introduction: f being given, such bounds r, R always exist if the following asymptotic conditions are satisfied:

limτ+f(τ)τ=0,limτ0f(τ)τ=+.

3.3 The contraction condition

Assume that f is Lipschitz continuous on the interval [r, R] , more precisely that the following condition holds:

(h3) There exists a constant 0l<αλ1/ g L(Ω) such that

(3.5) | f(τ)fτ¯ || ττ¯ |forallτ,τ¯[ r,R ].

Here λ1 is the first eigenvalue of the Dirichlet problem for Δ.

We recall that imposing on f the Lipschitz continuity condition is natural for getting uniqueness in nonlinear problems. Nevertheless, in our case, this condition is assumed only locally, as is often the case in real-world applications, and mathematically it makes the possibility to study by localization problems having multiple solutions.

For any fixed ε > 0, we prove that, under hypothesis (h3), the operator Nε satisfies the contraction property on DrR. To this aim, let v1, v2 ∈ DrR and denote wi = Nε (vi), i = 1, 2. Then

d i v A ε w i = g f v i in  Ω w i = 0 on  Ω

in the weak sense, for i = 1, 2. These give

d i v A ε w 1 w 2 = g f v 1 f v 2 in  Ω w 1 w 2 = 0 on  Ω

in the weak sense. Thus, multiplying by w1 – w2 and integrating over Ω yield

(3.6) Ω A ε w 1 w 2 w 1 w 2 = Ω g f v 1 f v 2 w 1 w 2 .

For the left-hand side integral in (3.6), by the ellipticity condition (2.2) on Aε, one has

(3.7) α w 1 w 2 H 0 1 Ω 2 Ω A ε w 1 w 2 w 1 w 2 .

For the right-hand side integral in (3.6), since g is zero outside and by using (h3), we have

Ωg[ f(v1)f(v2) ](w1w2)=Ωg[ f(v1)f(v2) ](w1w2)lΩg| v1v2 || w1w2 |l g L(Ω) v1v2 L2(Ω) w1w2 L2(Ω).

Two times application of Poincaré's inequality with the sharp constant 1/λ1 yields

(3.8) Ωg[ f(v1)f(v2) ](w1w2)l g L(Ω)λ1 v1v2 H01(Ω) w1w2 H01(Ω).

Now (3.6)-(3.8) imply

(3.9) w1w2 H01(Ω)l g L(Ω)αλ1 v1v2 H01(Ω),

which in view of (h3) shows that N ε is a contraction on DrR, with the contraction constant

L:=l g L(Ω)αλ1<1.

3.4 Existence and localization result

We are now in position to state the existence and localization result for the nonlinear problem (2.1).

Theorem 3.2

Assume that conditions (h1)-(h3) hold. Then there exists ε0 > 0 such that

  1. for any ε < ε0, problem (2.1) has a unique solution uε in DrR, which is in DrR the unique minimum point of the energy functional JΕ.

  2. uεuasε0, weakly in H01(Ω) and strongly in L2 (Ω) , where u is the unique function in DrR which solves the homogenized problem (2.5). Also u is in DrR the unique minimum point of the energy functional J, and

Jε(uε)J(u)asε0.

  1. Proof. (a) The existence and uniqueness in DrR of uε follows from Banach's contraction principle, in view of the previous considerations. In order to prove the convergence of uε to u, note that, due to the boundedness in L(Ω)oftheset{ gf(uε) },theset{ uε } is bounded in H01(Ω), and so relatively compact in L2 (Ω). Thus there is a subsequence (uεk) with εk0 which is weakly convergent in H01(Ω) to some uΗ01(Ω), and strongly convergent in L2 (Ω) to u. Obviously, the limit u belongs to DrR. Now we are in position to pass to the limit in problem (2.1). The limit of the left-hand side of (2.1) is found as in the linear case (see [5, Chapter 9]), while for the right-hand side of (2.1) we use the uniform boundedness of (uεk) in L (Ω), its strong convergence in L2 (Ω), and the continuity of Nemytskii's operator associated to f from L2 (Ω) to itself. Thus, the limit u satisfies problem (2.5). The uniqueness in DrR of the solution of the homogenized problem, can be proved with similar arguments as above, by means of Banach's contraction principle, and implies that the entire sequence (uε) converges to u, weakly in H01(Ω), and strongly in L2 (Ω).

  2. (b) To show that for each fixed ε, uε minimizes Jε in DrR, let us first note that Jε is a C1 functional on V and under the standard identification of H1(Ω)toH01(Ω), one has

(3.10) (Jε)(v)=vNε(v).

Hence the fixed points of Nε are the critical points of Jε. It is easy to see that Jε is bounded from below on DrR. Then, using Ekeland's variational principle (see, e.g., [6], [17]), there is a minimizing sequence (vk) of elements from DrR such that

(3.11) Jε(vk)infDrRJε+1k

and

(3.12) Jε(vk)Jε(v)+1k vvk H01(Ω)

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

vt=vk(Jε)(vk),0<t<1.

Using (3.10), one has

vt=(1t)vk+tNε(vk).

Here, one has v k ∈ DrR and by the invariance property in Section 3.2, N ε (vk) also belongs to DrR. Since DrR is convex, it follows that vtDrR for every t ∈ (0, 1). Replacing v by vt into (3.12) and then dividing by t, yields

t1(Jε(vk)Jε(vkt(Jε)(vk)))1k (Jε)(vk) H01(Ω),

whence letting t go to zero, one finds

J ε v k , J ε v k 1 k J ε v k H 0 1 Ω ,

that is

J ε v k H 0 1 Ω 1 k .

Hence

(3.13) w k := J ε v k 0 in  H 0 1 Ω .

Using (3.10), we have wk = vk – Nε (vk) and by the contraction property of Nε, we find

vk+pvk H01(Ω) Nε(vk+p)Nε(vk) H01(Ω)+ wk+pwk H01(Ω)L vk+pvk H01(Ω)+ wk+pwk H01(Ω),

whence

vk+pvk H01(Ω)11L wk+pwk H01(Ω).

This ensures that (vk) is a Cauchy sequence, so it converges in to some vε. Clearly vεDrR and from (3.11) and (3.13), passing to the limit with k we obtain that vε minimizes Jε in DrR and (Jε)' (vε) = 0. Hence vε is a fixed point in DrR of Nε. The uniqueness of the fixed point, guaranteed by Banach's contraction principle, implies that vε = uε as desired. The proof of the fact that the limit u, solution of the homogenized problem, minimizes J in DrR is analogous.

Finally the convergence Jε(uε)J(u) follows from a standard result in homogenization theory (see [5, Section 8.2]) and the properties of Nemytskii's operator.

Remark 3.2

(Multiplicity) If the conditions (h1)-(h3) are satisfied for two pairs of numbers (r, R), let them be (r1,R1) and (r2, R2), with 0 < r1 < R1 < r2 < R2, then the sets Dr1R1 and Dr2R2 are disjoint, and according to Theorem 3.2, problem (2.1), as well as the corresponding homogenized problem (2.5), have two distinct solutions, one in Dr1R1 and the other in Dr2R2. This argument, extended to a finite or infinite sequence of pairs (r, R) , can be used in order to obtain multiple (finitely or infinitely many) solutions of problem (2.1) and of the corresponding homogenized problem (2.5).

Remark 3.3

Under the assumption that g vanishes in Ω/Ω, we were able to localize the solutions assuming the Lipschitz continuity condition for f only locally, in the interval [r, R]. It is precisely this local Lipschitz condition which allows us to obtain, via Theorem 3.2, multiple solutions. If the assumption that g vanishes in Ω\Ω' does not hold, one can however localize the solutions in the set DrR, but under the stronger condition that f is Lipschitz continous in [0, R] instead of [r, R]. Correspondingly, in this case M f is

Mf=maxτ[ 0,R ]f(τ)

and the estimation from below (3.4) should be driven as follows:

Nε(v)(x)Sε(1Ωgf(v))(x)Sε(1Ωgmf)(x)=mfSε(1Ωg)(x)mfy1Ωg.

Here 11Ω' is the characteristic function of This implies that yg in (3.4) should be replaced by y1Ωg.

Nevertheless, under this stronger Lipschitz condition on f, Theorem 3.2 cannot produce multiple solutions. Indeed, if Theorem 3.2 would apply to two pairs of numbers, let them be (r1, R1) and (r2, R2), with 0 < r1 < R1 < r2 < R2, then the sets Dr1R1 and Dr2R2 being disjoints, Nε would have two distinct fixed points, one in Dr1R1 and other in Dr2R2. On the other hand, f should be Lipschitz continuous on [0, R2] with the Lipschitz constant l<αλ1/ g L(Ω). Then it is easy to see that the contraction condition on Nε, namely inequality (3.9), holds true for every v1,v2BR2:={ vV:0v(x)R2 and this guarantees that Nε has at most one fixed point in B R2. Since both disjoint sets Dr1R1 and Dr2R2 are included in BR2 , we get a contradiction. This justifies our claim about multiplicity.

4 Componentwise localization for systems

The technique used in Section 3 for one equation can be extended to systems, for the localization of each component of the solution, individually. For simplicity we discuss only a system of two equations in Ω, namely

(4.1) { div(A1ε(x)u1ε)=g1(x)f1(u1ε,u2ε)div(A2ε(x)u2ε)=g2(x)f2(u1ε,u2ε)

under the Dirichlet boundary condition uiε=0onΩ,i=1,2.

We assume on the matrices Aiε, and on the functions gi, fi, similar hypotheses as in Section 2, allowing the ellipticity constants for the two equations to be different, namely αi and βi, i = 1,2. Then, solving the Dirichlet problem for system (4.1) in the weak sense is equivalent to finding in V x V the fixed point of the problem

u1ε=N1ε(u1ε,u2ε)u2ε=N2ε(u1ε,u2ε)

where Niε=Siε(gifi),i=1,2. Also, in terms of the energy functionals associated to the equations of the system, solving problem (4.1) is equivalent to solving the critical point problem

J11ε(u1ε,u2ε)=0J22ε(u1ε,u2ε)=0,

where by Jiiε(v1,v2) we have denoted the Fréchet derivative of Jiε(v1,v2) with respect to the variable vi, for i = 1,2, and the expressions of the functionals Jiε(v1,v2) are

J1ε(v1,v2)=12ΩA1ε(x)v1v1Ωg1(x)0v1(x)f1(τ,v2(x))dτ,J2ε(v1,v2)=12ΩA2ε(x)v2v2Ωg2(x)0v2(x)f2(v1,(x),τ)dτ.

We are interested to guarantee the existence and uniqueness of a solution (u1ε,u2ε) in a bounded subset of V x V, more exactly (u1ε,u1ε)D:=Dr1R1×Dr2R2, where 0 < r1 < R1 and 0 < r2 < R2 are given numbers.

Here are the hypotheses that duplicate the assumptions made in Section 3 for a single equation:

(H1) gi(x)0a.e in Ω,gi(x)=0for xΩ/Ω, and gi>0 on a subset of nonzero measure; fi(τ1,τ2)0 for every τ1,τ20.

(H2) MiΓgiRiandrimiygi, where

Mi=maxτ1[ r1,R1 ]τ2[ r2,R2 ]fi(τ1,τ2),mi=minτ1[ r1,R1 ]τ2[ r2,R2 ]fi(τ1,τ2).

(H3) There exist constants 0 ≤ lij; (i,j = 1, 2) such that

| fi(τ1,τ2)fi(τ¯1,τ¯2) |li1| τ1τ¯1 |+li2| τ2τ¯2 |forallτ1,τ¯1[ r1,R1 ]andτ2,τ¯2[ r2,R2 ],

and the spectral radius of the matrix L=[ Lij ]i,j=1,2, where

Lij:=lij gi L(Ω)αiλ1(i,j=1,2)

is less than one.

Theorem 4.1

Assume that conditions (H1)-(H3) hold. Then there exists ε0 > 0 such that for any ε < ε0 :

  1. (i) system (4.1) has a unique solution (u1ε,u2ε)D, which is in D the unique Nash equilibrium point with respect to the energy functionals J1ε,J2ε,i.e.,

J1ε(u1ε,u2ε)=minDr1R1J1ε(,u2ε)andJ2ε(u1ε,u2ε)=minDr2R2J2ε(u1ε,).

  1. uiεuiasε0, weakly in H01(Ω) and strongly in L2 (Ω), where (u1, u2) is the unique couple of functions in D which solves the homogenized system

(4.2) { div(A10u1)=g1(x)f1(u1,u2)div(A20u2)=g2(x)f2(u1,u2)

Also (u1,u2 ) is in D the unique Nash equilibrium with respect to the energy functionals J1, J2 associated to the equations of the homogenized system (4.2), and

Jiε(u1ε,u2ε)Ji(u1,u2)asε0,fori=1,2.

Proof. (a) Existence, uniqueness and localization: Using similar arguments as in the case of a single equation, we obtain for any (v1,v2),(v¯,v¯2)Dandi{ 1,2 }:

0Niε(v1,v2)(x)MiSiε(gi)(x)MiΓgi,fora.a.xΩ,

and

(4.3) Niε(v1,v2)(x)miSiε(gi)(x)miygi,fora.a.xΩ,

which in view of (H2) show that Nε(D)D,whereNε=(N1ε,N2ε).

Also

Niε(v1,v2)Niε(v¯1,v¯2) H01(Ω)Li1 v1v¯1 H01(Ω)+Li2 v2v¯2 H01(Ω),i=1,2

which, using the matrix L defined in (H3), can be written in the matrix form

N 1 ε v 1 , v 2 N 1 ε v ¯ 1 , v ¯ 2 H 0 1 Ω N 2 ε v 1 , v 2 N 2 ε v ¯ 1 , v ¯ 2 H 0 1 Ω L v 1 v ¯ 1 H 0 1 Ω v 2 v ¯ 2 H 0 1 Ω .

Then, from Perov's fixed point theorem (see [14, Theorem 1]), the operator N ε has in D a unique fixed point (u1ε,u2ε), solution of system (4.1).

(b) Nash equilibrium: In order to prove that the solution (u1ε,u2ε) is a Nash equilibrium, we use an iterative approximation scheme. We start with a fixed element u2,0εinDr2R2. At each step k(k1),u2,k1ε been determined at step k – 1, first we apply Ekeland' s principle to J1ε(,u2,k1ε) and find an u1,kεDr1R1 such that

(4.4) J1ε(u1,kε,u2,k1ε)infDr1R1J1ε(,u2,k1ε)+1k, J11ε(u1,kε,u2,k1ε) H01(Ω)1k.

Next we apply Ekeland’s principle to J2ε(u1,kε,) and obtain an u2,kεDr2R2 with

(4.5) J2ε(u1,kε,u2,kε)infDr2R2J2ε(u1,k,ε)+1k, J22ε(u1,kε,u2,kε) H01(Ω)1k.

Let

vk:=J11ε(u1,kε,u2,k1ε)andwk:=J22ε(u1,kε,u2,kε).

From (4.4) and (4.5), one has vk,wk0inH01(Ω). Also, as in the case of one equation, we have

(4.6) u 1 , k ε N 1 ε u 1 , k ε , u 2 , k 1 ε = v k u 2 , k ε N 2 ε u 1 , k ε , u 2 , k ε = w k .

The first equality in (4.6) written for k and k + p, with p\{ 0 }, yields

(4.7) u 1 , k + p ε u 1 , k ε H 0 1 Ω N 1 ε u 1 , k + p ε , u 2 , k + p 1 ε N 1 ε u 1 , k ε , u 2 , k 1 ε H 0 1 Ω + v k + p v k H 0 1 Ω L 11 u 1 , k + p ε u 1 , k ε H 0 1 Ω + L 12 u 2 , k + p 1 ε u 2 , k 1 ε H 0 1 Ω + v k + p v k H 0 1 Ω L 11 u 1 , k + p ε u 1 , k ε H 0 1 Ω + L 12 u 2 , k + p ε u 2 , k ε H 0 1 Ω + L 12 u 2 , k + p 1 ε u 2 , k 1 ε H 0 1 Ω u 2 , k + p ε u 2 , k ε H 0 1 Ω + v k + p v k H 0 1 Ω .

The second equality in (4.6) written for k and k + p gives

(4.8) u2,k+pεu2,kε H01(Ω)L21 u1,k+pεu1,kε H01(Ω)+L22 u2,k+pεu2,kε H01(Ω)+ wk+pwk H01(Ω)

Denote

ak,p= u1,k+pεu1,kε H01(Ω),bk,p= u2,k+pεu2,kε H01(Ω),ck,p= vk+pvk H01(Ω),dk,p= wk+pwk H01(Ω).

Clearly,

(4.9) ck,p0anddk,p0ask,uniformlywithrespecttop.

With these notations, (4.7) and (4.8) become

ak,pL11ak,p+L12bk,p+L12(bk1,pbk,p),bk,pL21ak,p+L22bk,p+dk,p,

which can be expressed in the matrix form

(4.10) a k , p b k , p L a k , p b k , p + c k , p + L 12 b k 1 , p b k , p d k , p .

Under the assumption that the spectral radius of the matrix L is less than one, the matrix I – L (where I is the unit matrix), is invertible and its inverse contains only nonnegative elements (see [14, Lemma 2]). Thus (4.10) gives

a k , p b k , p I L 1 c k , p + L 12 b k 1 , p b k , p d k , p .

Let (IL)1 = [ρi j]. Then

(4.11) a k , p ρ 11 c k , p + L 12 b k 1 , p b k , p + ρ 12 d k , p b k , p ρ 21 c k , p + L 12 b k 1 , p b k , p + ρ 22 d k , p .

From the second inequality, one has

(4.12) bk,pρ21L121+ρ21L12bk1,p+ρ21ck,p+ρ22dk,p1+ρ21L12.

Note that the sequence (bk,p)k1 is bounded uniformly with respect to p. Indeed, from (4.6),

u2,kεwk=N2ε(u2,kε,u2,kε),

whence

{ div(Aε(x)(u2,kεwk))=g2(x)f2(u1,kε,u2,kε)inΩu2,kεwk=0onΩ.

Since the right-hand side is bounded in L (Ω) independently of k, this implies

u2,kεwk H01(Ω)c0

for every k. As a result, the sequence (u2,kε)k1 is bounded in the norm of H01(Ω) by some constant c. Then, clearly bk,p2c, as claimed. Now we recall a lemma proved in 16, Lemma 3.2 .

Lemma 4.2

Let (xk , p ) , (yk,P) be two sequences of real numbers depending on a parameter p, such that

(xk , p) is bounded uniformly with respect to p,

and

(4.13) 0xk,pλxk1,p+yk,p

for all k, p and some λ[ 0,1 ).Ifyk,p0ask uniformly with respect to p, then xk,p0 uniformly with respect to p.

We apply the previous lemma for

χk,p=bk,p,yk,p=ρ21ck,p+ρ22dk,p1+ρ21L12andλ=ρ21L121+ρ21L12bk1,p.

According to (4.12), one has (4.13), while due to (4.9), one has yk,p0ask, uniformly with respect to p. Also, obviously, 0λ<1. It follows that bk,p0ask, uniformly with respect to p, that is (u2,kε) is a Cauchy sequence. Next, the first inequality in (4.11) implies that (u1,kε) is also a Cauchy sequence. Let v1ε,v2ε be the limits of the sequences (u1,kε),(u2,kε)ask, respectively. Clearly (v1ε,v2ε)D, and passing to the limit in (4.4), (4.5) we obtain that (v1ε,v2ε) solves system (4.1). The uniqueness of the solution in D implies that (v1ε,v2ε)=(u1ε,u2ε).

References

[1] G. Allaire and H. Hutridurga, Upscating nonlinear adsorption in periodic porous media-homogenization approach Appl. Anal. 96 (10) (2016), 2126–2161.10.1080/00036811.2015.1038254Search in Google Scholar

[2] A. Bensoussan, J.L. Lions and G. Papanicolaou, Asymptotic Analysis for Periodic Structures North Holland, Amsterdam, 1978.Search in Google Scholar

[3] A. Braides and A. Defranceschi, Homogenization of Multiple Integrals Oxford University Press, Oxford, 1998.10.1093/oso/9780198502463.001.0001Search in Google Scholar

[4] H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations Springer, New York, 2011.10.1007/978-0-387-70914-7Search in Google Scholar

[5] D. Cioranescu and P. Donato, An Introduction in Homogenization Oxford University Press, 1999.10.1093/oso/9780198565543.001.0001Search in Google Scholar

[6] I. Ekeland, On the variational principle, J. Math. Anal. Appl. 47(1974), 324–353.10.1016/0022-247X(74)90025-0Search in Google Scholar

[7] M. Gahn, P. Knabner and M. Neuss-Radu, Homogenization of reaction-diffusion processes in a two-component porous medium with a nonlinear flux condition at the interface, and application to metabolic processes in cells SIAM J. Appl. Math. 76 (2016), 1819–1843.10.1137/15M1018484Search in Google Scholar

[8] A. Gaudiello and T. A. Mel'nyk, Homogenization of a nonlinear monotone problem with nonlinear Signorini boundary conditions in a domain with highly rough boundary,]. Differential Equations 265 (2018), 5419–5454.10.1016/j.jde.2018.07.002Search in Google Scholar

[9] D. Gilbarg and N.S. Trudinger, Elliptic Partial Differential Equations of Second Order Springer, Berlin, 1998.Search in Google Scholar

[10] I. Graf, M. A. Peter and J. Sneyd, Homogenization of a nonlinear multiscale model of calcium dynamics in biological cells J. Math. Anal. Appl. 419 (2014), 28–47.10.1016/j.jmaa.2014.04.037Search in Google Scholar

[11] J. Jost, Partial Differential Equations Springer, New York, 2007.10.1007/978-0-387-49319-0Search in Google Scholar

[12] V.A. Khoa and A. Muntean, Asymptotic analysis of a semi-linear elliptic system in perforated domains: Well–posedness and correctors for the homogenization limit, J. Math. Anal. Appl. 439 (2016) 271–295.10.1016/j.jmaa.2016.02.068Search in Google Scholar

[13] A. Piatnitski and M. Ptashnyk, Homogenization of biomechanical models for plan tissues Multiscale Model. Simul. 15(2017), 339–387.10.1137/15M1046198Search in Google Scholar

[14] R. Precup, The role of matrices that are convergent to zero in the study of semilinear operator systems Math. Comp. Modelling 49(2009), 703–708.10.1016/j.mcm.2008.04.006Search in Google Scholar

[15] R. Precup, Moser–Harnack inequality, Krasnosel'skii type fixed point theorems in cones and elliptic problems To pol. Methods Nonlinear Anal. 40 (2012), 301–313.Search in Google Scholar

[16] R. Precup, Nash-type equilibria and periodic solutions to nonvariational systems Adv. Nonlinear Anal. 3 (2014), 197–207.10.1515/anona-2014-0006Search in Google Scholar

[17] V.D. Rădulescu, Qualitative Analysis of Nonlinear Elliptic Partial Differential Equations: Monotonicity, Analytic, and Variational Methods Hindawi, New York, 2008.10.1155/9789774540394Search in Google Scholar

[18] C. Timofte, Homogenization results for the calcium dynamics in living cells Math. Comput. Simulat. 133 (2017), 165–174.10.1016/j.matcom.2015.06.011Search in Google Scholar

Received: 2018-10-26
Accepted: 2018-11-24
Published Online: 2019-05-22

© 2020 R. Bunoiu and R. Precup, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 Public License.

Articles in the same Issue

  1. Frontmatter
  2. On the moving plane method for boundary blow-up solutions to semilinear elliptic equations
  3. Regularity of solutions of the parabolic normalized p-Laplace equation
  4. Cahn–Hilliard equation on the boundary with bulk condition of Allen–Cahn type
  5. Blow-up solutions for fully nonlinear equations: Existence, asymptotic estimates and uniqueness
  6. Radon measure-valued solutions of first order scalar conservation laws
  7. Ground state solutions for a semilinear elliptic problem with critical-subcritical growth
  8. Generalized solutions of variational problems and applications
  9. Existence and non-existence results for Kirchhoff-type problems with convolution nonlinearity
  10. Nonlinear Sherman-type inequalities
  11. Global regularity for systems with p-structure depending on the symmetric gradient
  12. Homogenization of a net of periodic critically scaled boundary obstacles related to reverse osmosis “nano-composite” membranes
  13. Noncoercive resonant (p,2)-equations with concave terms
  14. Evolutionary quasi-variational and variational inequalities with constraints on the derivatives
  15. Sharp estimates on the first Dirichlet eigenvalue of nonlinear elliptic operators via maximum principle
  16. Localization and multiplicity in the homogenization of nonlinear problems
  17. Remarks on a nonlinear nonlocal operator in Orlicz spaces
  18. A Picone identity for variable exponent operators and applications
  19. On the weakly degenerate Allen-Cahn equation
  20. Continuity results for parametric nonlinear singular Dirichlet problems
  21. Construction of type I blowup solutions for a higher order semilinear parabolic equation
  22. Singularly perturbed Choquard equations with nonlinearity satisfying Berestycki-Lions assumptions
  23. Comparison results for nonlinear divergence structure elliptic PDE’s
  24. Constant sign and nodal solutions for parametric (p, 2)-equations
  25. Monotonicity formulas for coupled elliptic gradient systems with applications
  26. Berestycki-Lions conditions on ground state solutions for a Nonlinear Schrödinger equation with variable potentials
  27. A class of semipositone p-Laplacian problems with a critical growth reaction term
  28. The role of superlinear damping in the construction of solutions to drift-diffusion problems with initial data in L1
  29. Reconstruction of Tesla micro-valve using topological sensitivity analysis
  30. Lewy-Stampacchia’s inequality for a pseudomonotone parabolic problem
  31. Global well-posedness of nonlinear wave equation with weak and strong damping terms and logarithmic source term
  32. Regularity Criteria for Navier-Stokes Equations with Slip Boundary Conditions on Non-flat Boundaries via Two Velocity Components
  33. Homoclinics for singular strong force Lagrangian systems
  34. A constructive method for convex solutions of a class of nonlinear Black-Scholes equations
  35. On a class of nonlocal nonlinear Schrödinger equations with potential well
  36. Superlinear Schrödinger–Kirchhoff type problems involving the fractional p–Laplacian and critical exponent
  37. Regularity for minimizers for functionals of double phase with variable exponents
  38. Boundary blow-up solutions to the Monge-Ampère equation: Sharp conditions and asymptotic behavior
  39. Homogenisation with error estimates of attractors for damped semi-linear anisotropic wave equations
  40. A-priori bounds for quasilinear problems in critical dimension
  41. Critical growth elliptic problems involving Hardy-Littlewood-Sobolev critical exponent in non-contractible domains
  42. On the Sobolev space of functions with derivative of logarithmic order
  43. On a logarithmic Hartree equation
  44. Critical elliptic systems involving multiple strongly–coupled Hardy–type terms
  45. Sharp conditions of global existence for nonlinear Schrödinger equation with a harmonic potential
  46. Existence for (p, q) critical systems in the Heisenberg group
  47. Periodic traveling fronts for partially degenerate reaction-diffusion systems with bistable and time-periodic nonlinearity
  48. Some hemivariational inequalities in the Euclidean space
  49. Existence of standing waves for quasi-linear Schrödinger equations on Tn
  50. Periodic solutions for second order differential equations with indefinite singularities
  51. On the Hölder continuity for a class of vectorial problems
  52. Bifurcations of nontrivial solutions of a cubic Helmholtz system
  53. On the exact multiplicity of stable ground states of non-Lipschitz semilinear elliptic equations for some classes of starshaped sets
  54. Sign-changing multi-bump solutions for the Chern-Simons-Schrödinger equations in ℝ2
  55. Positive solutions for diffusive Logistic equation with refuge
  56. Null controllability for a degenerate population model in divergence form via Carleman estimates
  57. Eigenvalues for a class of singular problems involving p(x)-Biharmonic operator and q(x)-Hardy potential
  58. On the convergence analysis of a time dependent elliptic equation with discontinuous coefficients
  59. Multiplicity and concentration results for magnetic relativistic Schrödinger equations
  60. Solvability of an infinite system of nonlinear integral equations of Volterra-Hammerstein type
  61. The superposition operator in the space of functions continuous and converging at infinity on the real half-axis
  62. Estimates by gap potentials of free homotopy decompositions of critical Sobolev maps
  63. Pseudo almost periodic solutions for a class of differential equation with delays depending on state
  64. Normalized multi-bump solutions for saturable Schrödinger equations
  65. Some inequalities and superposition operator in the space of regulated functions
  66. Area Integral Characterization of Hardy space H1L related to Degenerate Schrödinger Operators
  67. Bifurcation of time-periodic solutions for the incompressible flow of nematic liquid crystals in three dimension
  68. Morrey estimates for a class of elliptic equations with drift term
  69. A singularity as a break point for the multiplicity of solutions to quasilinear elliptic problems
  70. Global and non global solutions for a class of coupled parabolic systems
  71. On the analysis of a geometrically selective turbulence model
  72. Multiplicity of positive solutions for quasilinear elliptic equations involving critical nonlinearity
  73. Lack of smoothing for bounded solutions of a semilinear parabolic equation
  74. Gradient estimates for the fundamental solution of Lévy type operator
  75. π/4-tangentiality of solutions for one-dimensional Minkowski-curvature problems
  76. On the existence and multiplicity of solutions to fractional Lane-Emden elliptic systems involving measures
  77. Anisotropic problems with unbalanced growth
  78. On a fractional thin film equation
  79. Minimum action solutions of nonhomogeneous Schrödinger equations
  80. Global existence and blow-up of weak solutions for a class of fractional p-Laplacian evolution equations
  81. Optimal rearrangement problem and normalized obstacle problem in the fractional setting
  82. A few problems connected with invariant measures of Markov maps - verification of some claims and opinions that circulate in the literature
Downloaded on 27.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/anona-2020-0001/html
Scroll to top button