Processing math: 63%
Research article Special Issues

Analysis of the impact of offshore wind power on the Japanese energy grid

  • Received: 08 December 2022 Revised: 28 January 2023 Accepted: 03 February 2023 Published: 15 February 2023
  • As part of its economy-wide decarbonization target towards 2050, Japan plans to increase renewable generation, especially offshore wind, for which the country has a high potential. However, this resource is currently under-developed as available turbines are prone to shut-downs and can even suffer damage during the passage of typhoons. With new typhoon proof (T-class) turbines being currently developed by various companies, Japan now aims to develop 10 GW of offshore wind between 2021 and 2030, and 91 GW in the long-term. This research estimates the impact of integrating offshore wind into the Japanese main power grid using T-class turbines by considering three scenarios. First, a business-as-usual (BAU) case with 10 GW offshore wind capacity (following the 6th Strategic Energy Plan of Japan). Second, an offshore wind capacity of 91 GW. Third, the 91 GW offshore capacity being redistributed amongst regions to maximize its integration opportunities (Scenario 2). The simulations were carried out using the Energy System simulation model (EnSym). The results show that the BAU and Scenario 1 resulted in offshore wind achieving 1.7% and 7.28% of generation share, respectively, increasing to 9.77% for Scenario 2. Increasing the share of offshore wind in the energy mix mainly replaced liquefied natural gas (LNG).

    Citation: YingTung Chen, Kristina Knüpfer, Miguel Esteban, Tomoya Shibayama. Analysis of the impact of offshore wind power on the Japanese energy grid[J]. AIMS Energy, 2023, 11(1): 110-134. doi: 10.3934/energy.2023006

    Related Papers:

    [1] Xiao Qing Huang, Jia Feng Liao . Existence and asymptotic behavior for ground state sign-changing solutions of fractional Schrödinger-Poisson system with steep potential well. Communications in Analysis and Mechanics, 2024, 16(2): 307-333. doi: 10.3934/cam.2024015
    [2] Yan-Fei Yang, Chun-Lei Tang . Positive and sign-changing solutions for Kirchhoff equations with indefinite potential. Communications in Analysis and Mechanics, 2025, 17(1): 159-187. doi: 10.3934/cam.2025008
    [3] Fangyuan Dong . Multiple positive solutions for the logarithmic Schrödinger equation with a Coulomb potential. Communications in Analysis and Mechanics, 2024, 16(3): 487-508. doi: 10.3934/cam.2024023
    [4] Maomao Wu, Haidong Liu . Multiple solutions for quasi-linear elliptic equations with Berestycki-Lions type nonlinearity. Communications in Analysis and Mechanics, 2024, 16(2): 334-344. doi: 10.3934/cam.2024016
    [5] Yonghang Chang, Menglan Liao . Nonexistence of asymptotically free solutions for nonlinear Schrödinger system. Communications in Analysis and Mechanics, 2024, 16(2): 293-306. doi: 10.3934/cam.2024014
    [6] Shengbing Deng, Qiaoran Wu . Existence of normalized solutions for the Schrödinger equation. Communications in Analysis and Mechanics, 2023, 15(3): 575-585. doi: 10.3934/cam.2023028
    [7] Rui Sun, Weihua Deng . A generalized time fractional Schrödinger equation with signed potential. Communications in Analysis and Mechanics, 2024, 16(2): 262-277. doi: 10.3934/cam.2024012
    [8] Floyd L. Williams . From a magnetoacoustic system to a J-T black hole: A little trip down memory lane. Communications in Analysis and Mechanics, 2023, 15(3): 342-361. doi: 10.3934/cam.2023017
    [9] Sergey A. Rashkovskiy . Nonlinear Schrödinger equation with a short-range compensating field. Communications in Analysis and Mechanics, 2025, 17(2): 520-549. doi: 10.3934/cam.2025021
    [10] Xueqi Sun, Yongqiang Fu, Sihua Liang . Normalized solutions for pseudo-relativistic Schrödinger equations. Communications in Analysis and Mechanics, 2024, 16(1): 217-236. doi: 10.3934/cam.2024010
  • As part of its economy-wide decarbonization target towards 2050, Japan plans to increase renewable generation, especially offshore wind, for which the country has a high potential. However, this resource is currently under-developed as available turbines are prone to shut-downs and can even suffer damage during the passage of typhoons. With new typhoon proof (T-class) turbines being currently developed by various companies, Japan now aims to develop 10 GW of offshore wind between 2021 and 2030, and 91 GW in the long-term. This research estimates the impact of integrating offshore wind into the Japanese main power grid using T-class turbines by considering three scenarios. First, a business-as-usual (BAU) case with 10 GW offshore wind capacity (following the 6th Strategic Energy Plan of Japan). Second, an offshore wind capacity of 91 GW. Third, the 91 GW offshore capacity being redistributed amongst regions to maximize its integration opportunities (Scenario 2). The simulations were carried out using the Energy System simulation model (EnSym). The results show that the BAU and Scenario 1 resulted in offshore wind achieving 1.7% and 7.28% of generation share, respectively, increasing to 9.77% for Scenario 2. Increasing the share of offshore wind in the energy mix mainly replaced liquefied natural gas (LNG).



    In the past decades, the following Schrödinger-Poisson system

    {Δu+V(x)u+ϕu=f(x,u)inR3,Δϕ=u2inR3 (1.1)

    has been studied extensively by many authors, where V:R3R and fC(R3×R,R). This system can be used to describe the interaction of a charged particle with the electrostatic field in quantum mechanics. In this context, the unknown u and ϕ represent the wave functions related to the particle and electric potentials, respectively. Moreover, the local nonlinearity f(x,u) models the interaction among particles. We refer the reader to [6,20] for more details on its physical background.

    It is worth noting that system (1.1) is a nonlocal problem due to the appearance of the term ϕu, where ϕ=ϕu is presented in (1.4) below. This fact states that problem (1.1) is no longer a pointwise identity and brings some essential difficulties. For example, the term R3ϕuu2dx in the corresponding energy functional is homogeneous of degree four, then, compared with the local Schrödinger equation, it seems difficult to obtain the boundedness and compactness for any Palais-Smale sequence. In light of the previous observations, the existence of solutions for problem (1.1) have been widely studied and some open problems have been proposed [3,11,15,16,19,25,28,30,35].

    In what follows, we are particularly interested in the existence of sign-changing solutions (also known as nodal solutions) for problem (1.1). From this perspective, Wang and Zhou [29] were concerned with the existence and energy property of sign-changing solutions for problem (1.1) with f(x,u)=|u|p2u. By introducing appropriate compactness conditions on V, they used methods different from [5] to prove that the so-called sign-changing Nehari manifold is nonempty provided that 4<p<6. Then, combining some analytical techniques and the Brouwer degree theory, the existence of least energy sign-changing solutions was established. After that, the authors in [21] investigated sign-changing solutions of problem (1.1) when fC1(R,R) satisfied super-cubic and subcritical growth at infinity, superlinear growth at origin, and a well-known Nehari-type monotonicity condition. In particular, they established the energy doubling [31]. Moreover, the authors in [10,38] obtained the similar existence results if the nonlinearity f satisfied asymptotically cubic and three-linear growth, respectively. On the other hand, when f satisfies three-sublinear growth, the existence and multiplicity of sign-changing solutions can be obtained by invariant sets of descending flow [13,18]. For more interesting results, such as the Sobolev critical exponent or bounded domains, we refer to [1,24,27,34,36,37] and the references therein.

    According to the previous statements, we observe that the nonlinearities always satisfy superlinear growth or convexity (i.e. f(x,u)=|u|p2u, 2<p<6) provided that the sign-changing solution of Schrödinger-Poisson systems in the whole space R3 is considered. Once the nonlinearity is not constrained by the above forms, the methods mentioned previously cannot be directly used. Therefore, in present paper, we focus on a special type of nonlinearities; that is, the concave-convex type, such as f(x,u)=|u|p2u+|u|q2u with 4<p<6 and 1<q<2. The concave-convex nonlinearities were introduced in [2], where the authors proved the existence of infinitely many solutions with negative energy for local elliptic problems in bounded domains. After this work, a great attention has been paid to the existence of solutions to elliptic problems with concave-convex nonlinearities. For example, see [7,8,17,32] for local Schrodinger equations, and [9,14,22,23,26,33] for Schrodinger-Poisson systems.

    Note that only [7,8,17,33] involve the sign-changing solutions. More precisely, Bobkov [7] considered the following Schrödinger equation

    {Δu=λ|u|q2u+|u|γ2uin Ω,u=0inΩ,

    where ΩRN is a bounded connected domain with a smooth boundary, N1, 1<q<2<γ<2 and 2 is the well-known Sobolev critical exponent. They proved the existence of a sign-changing solution on the nonlocal interval λ(,λ0), where λ0 is determined by the variational principle of nonlinear spectral analysis through the fibering method. Moreover, the author in [8] obtained similar existence results and some interesting properties for the nodal solutions of the elliptic equation

    {Δu=λk(x)|u|q2u+h(x)|u|γ2uinΩ,u=0inΩ,

    where 1<q<2<γ<2, λR and the weight functions k, hL(Ω) satisfy the conditions essinfxΩk(x)>0 and essinfxΩh(x)>0. Note that the methods in [7,8] cannot be applied to the nonlocal elliptic problem (1.1). To this end, based on the setting of bounded domains, Yang and Ou [33] studied the following Schrödinger-Poisson system

    {Δu+ϕu=λ|u|p2u+|u|q2uinΩ,Δϕ=u2inΩ,u=0inΩ, (1.2)

    where Ω is a bounded domain with smooth boundary Ω in R3 and 1<p<2, 4<q<6, λ is a constant. By constrained variational method and quantitative deformation lemma, they obtained that the problem (1.2) has a nodal solution uλ with positive energy when λ<λ, λ is a constant. Here, we point out that if the bounded domain is involved, the embedding H10(Ω)Lp(Ω) is compact for 1p<2, which not only avoids the verification of compactness but also ensures the boundedness of the concave term. However, once the whole space is considered, these points cannot be directly determined. Therefore, motivated by the works described above, in this paper we focus on the following Schrödinger-Poisson system in the whole space R3 with concave-convex nonlinearities

    {Δu+V(x)u+ϕu=|u|p2u+λK(x)|u|q2uinR3,Δϕ=u2inR3, (1.3)

    where 1<q<2, 4<p<6, λ>0 and V, K satisfy the assumptions:

    (V) VC(R3,R) satisfies infxR3V(x)a>0 for each A>0, meas{xR3:V(x)A}<, where a is a constant and meas denotes the Lebesgue measure in R3;

    (K) K is positive and KL66q(R3).

    Here the condition (V) is similar to [17]. This condition also ensures the compactness of embedding HLp(R3), 2p<2, where H is the Hilbert space

    H={uH1(R3):R3V(x)u2dx<+}

    endowed with the norm

    u=(R3(|u|2+V(x)u2)dx)12[4,39,  Lemma 3.4].

    Meanwhile, we point out that the authors in [17] considered the local Schrödinger type equation in RN

    Δu+V(x)u=λ|u|q2u+μu+ν|u|p2u,

    where 1<q<2<p<2, N2 and λ,μ,ν are parameters. The above equation involves a combination of concave and convex terms. They obtained infinitely many nodal solutions by using the method of invariant sets. However, it seems that this method cannot be applied to problem (1.3). In order to overcome the previous difficulties, we introduce the condition (K), which guarantees a weak continuity result (see Lemma 2.2 below). Moreover, conditions (V) and (K) allow us to construct a suitable nonempty closed subset of sign-changing Nehari manifold similar to [33], and then a least energy sign-changing solution can be obtained.

    Before proceeding, we discuss the basic framework for dealing with our problem. The usual norm in the Lebesgue space Lr(R3) is denoted by |u|r=(R3|u|rdx)1r, r[1,+). It is well known that, by the Lax-Milgram theorem, when uH, there exists a unique ϕuD1,2(R3) such that Δϕu=u2, where

    ϕu(x)=14πR3u2(y)|xy|dy. (1.4)

    Substituting (1.4) into (1.3), we can rewrite system (1.3) as the following equivalent form

    Δu+V(x)u+ϕuu=|u|p2u+λK(x)|u|q2u   inR3. (1.5)

    Therefore, the energy functional associated with system (1.3) is defined by

    Iλ(u)=12R3(|u|2+V(x)u2)dx+14R3ϕuu2dx1pR3|u|pdxλqR3K(x)|u|qdx,uH.

    The functional Iλ(u) is well-defined for every uH and belongs to C1(H,R). Furthermore, for any vH,

    Iλ(u),v=R3(uv+V(x)uv)dx+R3ϕuuvdxR3|u|p2uvdxλR3K(x)|u|q2uvdx.

    As is well known, the solution of problem (1.5) is the critical point of the functional Iλ(u). Moreover, if uH is a solution of problem (1.5) and u±0, then u is a sign-changing solution of system (1.3), where

    u+(x)=max{u(x),0}andu(x)=min{u(x),0}.

    Naturally, we introduce the Nehari manifold of Iλ as

    Nλ={uH{0}:Iλ(u),u=0},

    which is related to the behavior of the map φu:rIλ(ru) (r>0) (see [12] for the introduction of this map). For uH, we have

    φu(r)=12r2||u||2+14r4R3ϕuu2dx1prp|u|ppλqrqR3K(x)|u|qdx.

    It is well known that, for any uH{0}, φu(r)=0 if and only if ruNλ, which also implies that φu(1)=0 if and only if uNλ. This manifold is always used to find the positive ground state solution. In order to obtain sign-changing solutions of problem (1.3), it is necessary to consider the sign-changing Nehari manifold

    Mλ={uH:u±0,Iλ(u),u±=0}.

    Hoeever, this manifold cannot be directly applied due to appearance of concave term λk(x)|u|q2u. As we will see, inspired by [33], we can construct the set MλMλ and prove this set is a nonempty closed set, where Mλ is defined by (2.5) below. We then show that the minimization problem mλ:=infuMλIλ(u) is attained by some uλMλ with positive energy. Finally, the classical deformation lemma [32, Lemma 2.3] states that the uλ is a weak solution of problem (1.3). Up to now, the main results can be stated as follows.

    Theorem 1.1. Assume that (V) and (K) hold. Then there exists a constant λ>0 (determined in (2.13)) such that for any λ(,λ), problem (1.3) possesses a least energy sign-changing solution uλ with positive energy.

    Remark 1.2. As mentioned previously, our Theorem 1.1 extends the result of [33] to the whole space R3. Moreover, for nonlinearities that do not involve concave terms, such as f(x,u)=|u|p2u (4<p<6), with the aid of classic techniques in [29], one can obtain that the sign-changing Nehari manifold is nonempty. However, in our paper, 1<q<2 and 4<p<6 mean that λK(x)|u|q2u is concave and |u|p2u is convex, which is different from [29]. At this point, it is difficult to directly prove that the set Mλ is nonempty. To this end, we carefully analyze the behavior of fu(s,t)=Iλ(su++tu) and then introduce the set Mλ. In particular, by determining some important lower bound estimates, we verify that Mλ. On the other hand, we point out that the range of parameter λ can be negative, which is similar to previous results.

    The remainder of this paper is organized as follows. In section 2, we present some preliminary lemmas that are crucial for proving our main results. Section 3 is devoted to proving Theorem 1.1.

    In this section, we present some preliminary lemmas that are crucial for proving our main results. First, we recall some well-known properties of ϕu that are a collection of results in [11,19].

    Lemma 2.1. For any uH, we get

    (i) there exists C>0 such that R3ϕuu2dxCu4;

    (ii) ϕu0, for any uH;

    (iii) ϕtu=t2ϕu, for any t>0 and uH;

    (iv) if unu in H, then ϕunϕu in D1,2(R3) and

    limnR3ϕunu2ndx=R3ϕuu2dx.

    Next, we verify a weak continuity of the concave term. The proof is similar to [32, Lemma 2.13], but we state the proof here for the readers convenience.

    Lemma 2.2. Assume that 1<q<2 and (K) hold, then the functional

    G:HR:uR3K(x)|u|qdx

    is weakly continuous.

    Proof. Undoubtedly, it is sufficient to prove that if unu in H, then R3K(x)|un|qdxR3K(x)|u|qdx as n. In fact, if unu in H, going if necessary to a subsequence, we can assume that unu a.e. on R3. Since unu in H, we get that {un} is bounded in L6(R3) and {uqn} is bounded in L6q(R3). Therefore, uqnuq in L6q(R3). Combining with (K) and the definition of weak convergence, we obtain

    R3K(x)|un|qdxR3K(x)|u|qdxas  n.

    Now, for any uH with u±0, we introduce the map fu:[0,+)×[0,+)R defined by fu(s,t)=Iλ(su++tu), i.e.,

    fu(s,t)=12s2||u+||2+14s4R3ϕu+(u+)2dx+12s2t2R3ϕu+(u)2dx1psp|u+|ppλqsqR3K(x)|u+|qdx+12t2||u||2+14t4R3ϕu(u)2dx1ptp|u|ppλqtqR3K(x)|u|qdx. (2.1)

    Here we have used the fact that

    R3ϕu+|u|2dx=R3ϕu|u+|2dx.

    Moreover, we have that

    fu(s,t)=(Iλ(su++tu),u+,Iλ(su++tu),u)=(1sIλ(su++tu),su+,1tIλ(su++tu),tu),

    which implies that for any uH with u±0, su++tuMλ if and only if the pair (s,t) is a critical point of fu.

    Lemma 2.3. Assume that 1<q<2, 4<p<6 and the assumption (K) hold, then there exists a constant λ1>0 such that for any uH with u±0, there holds that

    (i) if λ(0,λ1), then for any fixed t0, fu(s,t) has exactly two critical points, 0<s1(t)<s2(t); s1(t) is the minimum point and s2(t) is the maximum point; moreover, if λ0, then for any fixed t0, fu(s,t) has exactly one critical point, s3(t)>0, and it is the maximum point;

    (ii) if λ(0,λ1), then for any fixed s0, fu(s,t) has exactly two critical points, 0<t1(s)<t2(s); t1(s) is the minimum point and t2(s) is the maximum point; moreover, if λ0, then for any fixed s0, fu(s,t) has exactly one critical point, t3(s)>0, and it is the maximum point.

    Proof. We define fμ(s,t):[0,+)×[0,+)R by

    fμ(s,t)=12||su++tu||2+14μR3ϕsu++tu(su++tu)2dx1p|su++tu|ppλqR3K(x)|su++tu|qdx,

    where μ is a nonnegative parameter.

    (i) For any fixed t0, a direct calculation gives

    fμs(s,t)=s||u+||2+μs3R3ϕu+(u+)2dx+μst2R3ϕu+(u)2dxsp1|u+|ppλsq1R3K(x)|u+|qdx=sq1(s2q||u+||2+μs4qR3ϕu+(u+)2dx+μs2qt2R3ϕu+(u)2dxspq|u+|ppλR3K(x)|u+|qdx).

    Then, if s>0, fμs(s,t)=0 is equivalent to

    βμ(s)=s2q||u+||2+μs4qR3ϕu+(u+)2dx+μs2qt2R3ϕu+(u)2dxspq|u+|ppλR3K(x)|u+|qdx=0.

    For βμ(s), we can obtain that

    βμ(s)=s1q((2q)||u+||2+μ(4q)s2R3ϕu+(u+)2dx+μ(2q)t2R3ϕu+(u)2dx(pq)sp2|u+|pp).

    Clearly, for any fixed t0, 1<q<2 and 4<p<6, we can infer that βμ has exactly one critical point sμ>0, where sμ is related to t for any μ0. Moreover, βμ is strictly increasing in (0,sμ) and strictly decreasing in (sμ,+).

    Noting that if μ=0, we have

    β0(s)=s2q||u+||2spq|u+|ppλR3K(x)|u+|qdx

    and

    β0(s)=s1q((2q)||u+||2(pq)sp2|u+|pp).

    Hence,

    s0=((2q)||u+||2(pq)|u+|pp)1p2

    and

    β0(s0)=p2pq(2qpq)2qp2||u+||2(pq)p2|u+|p(2q)p2pλR3K(x)|u+|qdx.

    Let

    αμ(s)=s2q||u+||2+μs4qR3ϕu+(u+)2dx+μs2qt2R3ϕu+(u)2dxspq|u+|pp

    and

    λ+1=infuH,u±0α0(s0)R3K(x)|u+|qdx,

    then it follows from Sobolev embedding that

    λ+1p2pq(2qpq)2qp21Sq6Sp(2q)p2p|K|66q>0. (2.2)

    Therefore, if λ(0,λ+1), we can deduce βμ(s0)>β0(s0)>0 for any uH with u±0. For any fixed t0, if λ(0,λ+1), there exist unique s1(t) and s2(t) with 0<s1(t)<s2(t), such that βμ(s)=0 and βμ(s1(t))>0, βμ(s2(t))<0. On the other hand, when λ0, there is a unique s3(t)>0 such that βμ(s)=0 and βμ(s3(t))<0.

    Finally, considering that

    2fμs2(s,t)=||u+||2+3μs2R3ϕu+(u+)2dx+μt2R3ϕu+(u)2dx(p1)sp2|u+|ppλ(q1)sq2R3K(x)|u+|qdx=sq2(s2q||u+||2+3μs4qR3ϕu+(u+)2dx+μs2qt2R3ϕu+(u)2dx(p1)spq|u+|ppλ(q1)R3K(x)|u+|qdx),

    we can find that

    2fμs2(s,t)=(q1)sq2βμ(s)+sq1βμ(s).

    Note that fμs(s,t)=0 is equivalent to βμ(s)=0, then we have

    2fμs2(s,t)=sq1βμ(s)  if  fμs(s,t)=0.

    Thus, the facts that βμ(s1(t))>0,βμ(s2(t))<0 and βμ(s3(t))<0 signify that

    2fμs2(s1(t),t)>0,2fμs2(s2(t),t)<0and2fμs2(s3(t),t)<0.

    In particular, when μ=1, the results still hold.

    (ii) For any fixed s0, let

    λ1=infuH,u±0α0(s0)R3K(x)|u|qdx=infuH,u±0p2pq(2qpq)2qp2||u||2(pq)p2|u|p(2q)p2pR3K(x)|u|qdx. (2.3)

    Analogously with the proof (i), we can derive the conclusion.

    At last, it is easy to see that λ+1=λ1. Indeed, u±0 indicates that (u)±0, which yields that λ+1=λ1. Let λ1=λ+1=λ1, from (i) and (ii), then the proof is completed.

    Let

    λ2=infuH{0}p2pq(2qpq)2qp2||u||2(pq)p2|u|p(2q)p2pR3K(x)|u|qdx.

    Undoubtedly,

    λ1λ2p2pq(2qpq)2qp21Sq6Sp(2q)p2p|K|66q>0. (2.4)

    Here, the notation Sp represents the embedding constant of HLp(R3), which has a value depending on p[2,6]. According to Lemma 2.3, the following corollary is a direct result.

    Corollary 2.4. Assume that 1<q<2, 4<p<6, the assumption (K) and 0<λ<λ2 hold, then for any uH{0}, φu(r) has exactly two critical points, 0<r1(u)<r2(u) and φu(r1(u))>0, φu(r2(u))<0. On the other hand, when λ0, φu(r) has exactly one critical point, r3(u)>0, and φu(r3(u))<0.

    Lemma 2.5. Assume that 1<q<2, 4<p<6, the assumption (K) and λ<λ1 hold, then for any uMλ, (2fu/s2)(1,1)0 and (2fu/t2)(1,1)0. Moreover, Iλ(u)+ as ||u||+, i.e., the functional Iλ is coercive and bounded from below on Nλ.

    Proof. If 0<λ<λ1, from Lemma 2.3, it follows that fu(s,1) has exactly two critical points s1(1), s2(1) and 2fus2(s1(1),1)>0, 2fus2(s2(1),1)<0. Since uMλ, we have fus(1,1)=0, which means that s1(1)=1 or s2(1)=1. Hence, 2fus2(1,1)0. Analogously, we can conclude that 2fut2(1,1)0.

    If λ0, it follows from Lemma 2.3 that fu(s,1) has exactly one critical point s3(1) and 2fus2(s3(1),1)<0. Combining with uMλ, we have fus(1,1)=0, which shows that s3(1)=1. Therefore, we get 2fus2(1,1)0. Similarly, we can deduce that t3(1)=1, 2fut2(1,t3(1))<0 and the claim is clearly true.

    Note that uMλNλ, then the Sobolev embedding indicates that

    Iλ(u)=Iλ(u)14Iλ(u),u=14||u||2+(141p)R3|u|pdx+λ(141q)R3K(x)|u|qdx14||u||2+λ(141q)|K|66q|u|q614||u||2+λ1(141q)|K|66qSq6||u||q.

    Hence, combining 1<q<2 and 4<p<6, we derive Iλ(u)+ as ||u||+. That is, the functional Iλ(u) is coercive and bounded from below on Nλ. The proof is completed.

    Similarly, we obtain the following result.

    Corollary 2.6. Assume that 1<q<2, 4<p<6, the assumption (K) and λ<λ2 hold, then for any uNλ, φu(1)0.

    In what follows, we construct the following sets

    Mλ={uMλ:2fus2(1,1)<0,2fut2(1,1)<0}

    and

    Mλ={uMλ:2fus2(1,1)<0,2fut2(1,1)<0,φu(1)<0}. (2.5)

    Unquestionably, MλMλ. According to the properties of fu mentioned above, we can verify that the set Mλ is nonempty and Mλ=Mλ (see Lemma 2.9). To this end, we first get the following fact.

    Lemma 2.7. Assume that 1<q<2, 4<p<6 and the assumption (K) hold, there exists σ>0, which is independent of u and λ, such that

    ||u±||>σ>0

    for any uMλ.

    Proof. For any uMλ, from 2fus2(1,1)<0, 2fut2(1,1)<0 and Sobolev embedding, it follows that

    (2q)||u±||2<(2q)||u±||2+(4q)R3ϕu±(u±)2dx+(2q)R3ϕu+(u)2dx<(pq)|u±|pp(pq)Spp||u±||p,

    which implies

    ||u±||>(2q(pq)Spp)1p2:=σ>0. (2.6)

    Hence, the proof is finished.

    In order to prove that Mλ, let us define

    λ3=infuMλ{(p2)||u+||2+(p4)R3ϕu+(u+)2dx(pq)R3K(x)|u+|qdx,(p2)||u||2+(p4)R3ϕu(u)2dx(pq)R3K(x)|u|qdx}.

    From Sobolev embedding and Lemma 2.7, it follows that

    (p2)||u±||2+(p4)R3ϕu±(u±)2dx(pq)R3K(x)|u±|qdx(p2)||u±||2q(pq)Sq6|K|66q>(p2)σ2q(pq)Sq6|K|66q>0,

    where σ is given by (2.6). Therefore, we have

    λ3(p2)σ2q(pq)Sq6|K|66q>0. (2.7)

    Furthermore, we compute that

    2fus2(s,t)=||u+||2+3s2R3ϕu+(u+)2dx+t2R3ϕu+(u)2dx(p1)sp2|u+|ppλ(q1)sq2R3K(x)|u+|qdx,
    2fut2(s,t)=||u||2+3t2R3ϕu(u)2dx+s2R3ϕu+(u)2dx(p1)tp2|u|ppλ(q1)tq2R3K(x)|u|qdx

    and

    2fust(s,t)=2stR3ϕu+(u)2dx,2futs(s,t)=2stR3ϕu+(u)2dx.

    For any uMλ, we obtain

    2fus2(1,1)=||u+||2+3R3ϕu+(u+)2dx+R3ϕu+(u)2dx(p1)|u+|ppλ(q1)R3K(x)|u+|qdx=(2q)||u+||2+(4q)R3ϕu+(u+)2dx+(2q)R3ϕu+(u)2dx(pq)|u+|pp=(2p)||u+||2+(4p)R3ϕu+(u+)2dx+(2p)R3ϕu+(u)2dxλ(qp)R3K(x)|u+|qdx

    and

    2fut2(1,1)=||u||2+3R3ϕu(u)2dx+R3ϕu+(u)2dx(p1)|u|ppλ(q1)R3K(x)|u|qdx=(2q)||u||2+(4q)R3ϕu(u)2dx+(2q)R3ϕu+(u)2dx(pq)|u|pp=(2p)||u||2+(4p)R3ϕu(u)2dx+(2p)R3ϕu+(u)2dxλ(qp)R3K(x)|u|qdx.

    Lemma 2.8. Assume that 1<q<2, 4<p<6, the assumption (K) and λ<λ1 hold, then for any uH with u±0, there exists a unique pair (su,tu)R+×R+ such that suu++tuuMλ. Moreover, if λ<λ3, then Iλ(suu++tuu)=maxs,t>0Iλ(su++tu).

    Proof. First, we only prove the case of 0<λ<λ1 since the proof of λ0 is very similar. Let uH with u±0, from the proof of Lemma 2.3, then we have that fus(s,t) satisfies the conditions:

    (i) fut(s,t2(s))=0 for all s0;

    (ii) fut(s,t) is continuous and has continuous partial derivatives in [0,+)×[0,+);

    (iii) 2fut2(s,t2(s))<0 for all s0.

    Hence, we can obtain that if 0<λ<λ1, fut(s,t)=0 determines an implicit function t2(s) with continuous derivative on [0,+) by using the implicit function theorem. Analogously, if 0<λ<λ1, fus(s,t)=0 determines an implicit function s2(t) with continuous derivative on [0,+).

    On the other hand, for every s0, from fut(s,t2(s))=0 and fut(s,t)0 for sufficiently large t>0, we can show that

    t2(s)<s  for large enough s. (2.8)

    Otherwise, if t2(s)s, where s is large enough, it follows from the definition of fut(s,t) that fut(s,t2(s))<0, which contradicts fut(s,t2(s))=0. Similarly, we get

    s2(t)<t for sufficiently large t. (2.9)

    Therefore, by (2.8), (2.9), t2(0)>0,s2(0)>0, the continuity of t_2(s) and s_2(t) , we conclude that the curves of t_2(s) and s_2(t) must intersect at some point (s_u, t_u)\in \mathbb{R}^+\times\mathbb{R}^+ . That is, \frac{\partial f_u}{\partial t}(s_u, t_u) = \frac{\partial f_u}{\partial s}(s_u, t_u) = 0 . Additionally, noting that

    t_2^\prime(s) = -\left(\frac{\partial ^2f_u}{\partial t\partial s}\bigg/\frac{\partial ^2f_u}{\partial t^2}\right)(s, t_2(s)) > 0

    for any s > 0 , we obtain that the function t_2(s) is strictly increasing in (0, +\infty) . Similarly, the function s_2(t) is strictly increasing in (0, +\infty) . Consequently, there is a unique pair (s_u, t_u)\in \mathbb{R}^+\times\mathbb{R}^+ such that

    \frac{\partial f_u}{\partial s}(s_u, t_u) = \frac{\partial f_u}{\partial t}(s_u, t_u) = 0

    and

    \frac{\partial ^2f_u}{\partial s^2}(s_u, t_u) < 0, \; \frac{\partial ^2f_u}{\partial t^2}(s_u, t_u) < 0;

    that is, s_uu^++t_uu^-\in \mathcal{M}_\lambda^- .

    Next, we prove that (s_u, t_u) is the unique maximum point of f_u(s, t) on [0, +\infty)\times[0, +\infty) . In fact, if u\in\mathcal{M}_\lambda^- , we only show that (s_u, t_u) = (1, 1) is the pair of numbers such that I_\lambda(s_uu^++t_uu^-) = \max\limits_{s, t > 0}I_\lambda(su^++tu^-) . Define

    H(u) = \left(\frac{\partial ^2f_u}{\partial s^2}\frac{\partial ^2f_u}{\partial t^2}-\frac{\partial ^2f_u}{\partial t\partial s}\frac{\partial ^2f_u}{\partial s\partial t}\right)\bigg|_{(1, 1)}.

    If we verify that H(u) > 0 , then (1, 1) is a local maximum point of f_u(s, t) . Combining uniqueness of (s_u, t_u) , we have (1, 1) as a global maximum point of f_u(s, t) . Let u\in\mathcal{M}_\lambda^- , then

    \begin{equation*} \begin{split} &\quad \; H(u)\\&\! = \left(\frac{\partial ^2f_u}{\partial s^2}\frac{\partial ^2f_u}{\partial t^2}\!-\frac{\partial ^2f_u}{\partial t\partial s}\frac{\partial ^2f_u}{\partial s\partial t}\right)\bigg|_{(1, 1)}\\&\! = \!\Big((2-p)||u^+||^2\!\!+(4-p)\!\int_{\mathbb{R}^{3}}\!\phi_{u^+} {(u^+)}^2\mathrm{d}x\!+(2-p)\!\int_{\mathbb{R}^{3}}\!\phi_{u^+} {(u^-)}^2\mathrm{d}x\!-\lambda(q-p)\!\int_{\mathbb{R}^{3}}\!K(x)|u^+|^{q}\mathrm{d}x\!\Big)\\&\!\times\!\Big((2-p)||u^-||^2\!+(4-p)\!\int_{\mathbb{R}^{3}}\!\phi_{u^-} {(u^-)}^2\mathrm{d}x\!+(2-p)\!\int_{\mathbb{R}^{3}}\!\phi_{u^+} {(u^-)}^2\mathrm{d}x\!-\lambda(q-p)\!\int_{\mathbb{R}^{3}}\!K(x)|u^-|^{q}\mathrm{d}x\!\Big)\\&\!-4\!\left(\int_{\mathbb{R}^{3}}\phi_{u^+} {(u^-)}^2\mathrm{d}x\right)^2. \end{split} \end{equation*}

    From \frac{\partial ^2f_u}{\partial s^2}(1, 1) < 0 , \frac{\partial ^2f_u}{\partial t^2}(1, 1) < 0 , if \lambda < \lambda_3 , we derive

    \begin{equation*} \begin{split} &\quad-\frac{\partial ^2f_u}{\partial s^2}(1, 1)-2\int_{\mathbb{R}^{3}}\phi_{u^+}{(u^-)}^2\mathrm{d}x\\& = (p-2)||u^+||^2+(p-4)\int_{\mathbb{R}^{3}}\phi_{u^+} {(u^+)}^2\mathrm{d}x+(p-4)\int_{\mathbb{R}^{3}}\phi_{u^+} {(u^-)}^2\mathrm{d}x\\&\quad-\lambda(p-q)\int_{\mathbb{R}^{3}}K(x)|u^+|^{q}\mathrm{d}x\\& > (p-2)||u^+||^2+(p-4)\int_{\mathbb{R}^{3}}\phi_{u^+} {(u^+)}^2\mathrm{d}x-\lambda_3(p-q)\int_{\mathbb{R}^{3}}K(x)|u^+|^{q}\mathrm{d}x\\ & > 0 \end{split} \end{equation*}

    and

    \begin{equation*} \begin{split} &\quad-\frac{\partial ^2f_u}{\partial t^2}(1, 1)-2\int_{\mathbb{R}^{3}}\phi_{u^+}{(u^-)}^2\mathrm{d}x\\& > (p-2)||u^-||^2+(p-4)\int_{\mathbb{R}^{3}}\phi_{u^-} {(u^-)}^2\mathrm{d}x-\lambda_3(p-q)\int_{\mathbb{R}^{3}}K(x)|u^-|^q\mathrm{d}x\\& > 0, \end{split} \end{equation*}

    which show that H(u) > 0 .

    If u\not\in\mathcal{M}_\lambda^- , then there exists a unique pair (s_u^\prime, t_u^\prime) of positive numbers such that s_u^\prime u^++t_u^\prime u^-\in\mathcal{M}_\lambda^- . Let v = s_u^\prime u^++t_u^\prime u^- , i.e., v\in\mathcal{M}_\lambda^- . Repeat the above steps and we will get H(v) > 0 . Hence, the proof is completed.

    Lemma 2.9. If 1 < q < 2 , 4 < p < 6 , the assumption (K) and \lambda < \min{\{\lambda_1, \lambda_2, \lambda_3\}} hold, then \mathcal{M}_\lambda^*\neq \emptyset . Moreover, we get \mathcal{M}_\lambda^* = \mathcal{M}_\lambda^- .

    Proof. By the definitions of \mathcal{M}_\lambda^- and \mathcal{M}_\lambda^* , \mathcal{M}_\lambda^*\subset\mathcal{M}_\lambda^- is obvious. Hence, we only need to prove that if \lambda < \min{\{\lambda_1, \lambda_2, \lambda_3\}} , then \mathcal{M}_\lambda^-\subset\mathcal{M}_\lambda^* . That is, for any u\in \mathcal{M}_\lambda^- , \varphi_u reaches its maximum at point r = 1 . It follows from \lambda < \lambda_1 and Lemma 2.8 that \mathcal{M}_\lambda^-\neq \emptyset , and from Lemma 2.8, for any u\in\mathcal{M}_\lambda^- , we obtain H(u) > 0 when \lambda < \lambda_3 . Combining f_u(r, r) = \varphi_u(r) , it follows that r = 1 is a maximum of \varphi_u . Therefore, \mathcal{M}_\lambda^-\subset\mathcal{M}_\lambda^* . This completes the proof of Lemma 2.9.

    Corollary 2.10. If 1 < q < 2 , 4 < p < 6 , the assumption (K) and \lambda < \min{\{\lambda_1, \lambda_2, \lambda_3\}} hold, for u\in H and u^\pm\neq 0 , then there exists a unique pair (s_u, t_u)\in\mathbb{R}^+\times\mathbb{R}^+ such that s_uu^++t_uu^-\in\mathcal{M}_\lambda^* and I_\lambda(s_uu^++t_uu^-) = \max\limits_{s, t > 0}I_\lambda(su^++tu^-).

    Lemma 2.11. If 1 < q < 2 , 4 < p < 6 , the assumption (K) and \lambda < \lambda_4 hold, for all u\in H\backslash\{0\} , then there exists r_u > 0 such that \varphi_u(r_u) > 0 , where \lambda_4 > 0 .

    Proof. Fixed u\in H\backslash\{0\} , let

    \begin{equation*} E_u(r) = \frac{1}{2}r^2||u||^2-\frac{1}{p}r^p|u|_p^p \end{equation*}

    for any r\ge0 , then we have

    \begin{equation} \begin{split} \varphi_u(r)& = \frac{1}{2}r^2||u||^2+\frac{1}{4}r^4\int_{\mathbb{R}^{3}}\phi_{u} {u}^2\mathrm{d}x-\frac{1}{p}r^p|u|^p_p-\frac{\lambda}{q}r^q\int_{\mathbb{R}^{3}}K(x)|u|^{q}\mathrm{d}x\\&\ge \frac{1}{2}r^2||u||^2-\frac{1}{p}r^p|u|^p_p-\frac{\lambda}{q}r^q\int_{\mathbb{R}^{3}}K(x)|u|^{q}\mathrm{d}x\\& = E_u(r)-\frac{\lambda}{q}r^q\int_{\mathbb{R}^{3}}K(x)|u|^{q}\mathrm{d}x. \end{split} \end{equation} (2.10)

    Considering E_u(r) , we obtain that there is a unique r_1(u) = \Big(\frac{||u||^2}{|u|_p^p}\Big)^{\frac{1}{p-2}} > 0 such that E_u(r) achieves its maximum at r_1(u) and the maximum value is E_u(r_1(u)) = \frac{p-2}{2p}\Big(\frac{||u||}{\; |u|_p}\Big)^{\frac{2p}{p-2}} . Moreover, from Sobolev embedding and (2.10), it is clear to calculate that

    \begin{equation} \begin{split} \varphi_u(r_1(u))&\ge E_u(r_1(u))-\frac{\lambda}{q}(r_1(u))^q\int_{\mathbb{R}^{3}}K(x)|u|^{q}\mathrm{d}x\\&\ge E_u(r_1(u))-\frac{\lambda}{q}(r_1(u))^q|K|_{\frac{6}{6-q}}S_6^q||u||^q\\& = E_u(r_1(u))-\frac{\lambda}{q}S_6^q|K|_{\frac{6}{6-q}} \Big(\frac{2p}{p-2}\Big)^{\frac{q}{2}}(E_u(r_1(u)))^{\frac{q}{2}}\\& = (E_u(r_1(u)))^{\frac{q}{2}}\Big((E_u(r_1(u)))^{\frac{2-q}{2}} -\frac{\lambda}{q}S_6^q|K|_{\frac{6}{6-q}}\Big(\frac{2p}{p-2}\Big)^{\frac{q}{2}}\Big). \end{split} \end{equation} (2.11)

    Consequently, by taking

    \begin{equation} \begin{split} \lambda_4& = \frac{(p-2)q}{2p|K|_{\frac{6}{6-q}}S_6^q}\inf\limits_{u\in H\backslash\{0\}}\left(\frac{||u||}{\; |u|_p}\right)^{\frac{p(2-q)}{p-2}}\\&\ge\frac{(p-2)q}{2p|K|_{\frac{6}{6-q}}S_6^qS_p^{\frac{p(2-q)}{p-2}}} > 0, \end{split} \end{equation} (2.12)

    we conclude that if \lambda < \lambda_4 , it holds

    \begin{equation*} \begin{split} \frac{\lambda}{q}S_6^q|K|_{\frac{6}{6-q}}\left(\frac{2p}{p-2}\right)^ {\frac{q}{2}}& < \frac{\lambda_4}{q}S_6^q|K|_{\frac{6}{6-q}} \left(\frac{2p}{p-2}\right)^{\frac{q}{2}}\\&\le \frac{1}{q}S_6^q|K|_{\frac{6}{6-q}}\left(\frac{2p}{p-2}\right)^ {\frac{q}{2}}\frac{(p-2)q}{2pS_6^q|K|_{\frac{6}{6-q}}} \left(\frac{||u||}{\; |u|_p}\right)^{\frac{p(2-q)}{p-2}}\\& = (E_u(r_1(u)))^{\frac{2-q}{2}} \end{split} \end{equation*}

    for any u\in H\backslash\{0\} . This together with (2.11) yields that \varphi_u(r_1(u)) > 0 for any \lambda < \lambda_4 .

    Let

    \begin{equation} \lambda^* = \min{\{\lambda_1, \lambda_2, \lambda_3, \lambda_4\}}, \end{equation} (2.13)

    then it follows from (2.2), (2.4), (2.7) and (2.12) that \lambda^* > 0 . Now, we consider the properties of the set \mathcal{M}_\lambda^* .

    Lemma 2.12. If 1 < q < 2 , 4 < p < 6 , \lambda < \lambda^* and the assumptions (V) and (K) hold, then \mathcal{M}_\lambda^* is a closed set.

    Proof. Letting \{u_n\}\subset\mathcal{M}_\lambda^* satisfy u_n\rightarrow u_0 as n\rightarrow \infty in H , we now prove that u_0\in\mathcal{M}_\lambda^* . From \{u_n\}\subset\mathcal{M}_\lambda^* , we obtain

    \begin{equation} \langle I_\lambda^\prime(u_0), u_0^\pm\rangle = \lim\limits_{n\rightarrow \infty}\langle I_\lambda^\prime(u_n), u_n^\pm\rangle = 0, \end{equation} (2.14)
    \begin{equation} \frac{\partial ^2f_{u_0}}{\partial s^2}(1, 1) = \lim\limits_{n\rightarrow \infty}\frac{\partial ^2f_{u_n}}{\partial s^2}(1, 1)\le 0, \end{equation} (2.15)
    \begin{equation} \frac{\partial ^2f_{u_0}}{\partial t^2}(1, 1) = \lim\limits_{n\rightarrow \infty}\frac{\partial ^2f_{u_n}}{\partial t^2}(1, 1)\le 0, \end{equation} (2.16)
    \begin{equation} \varphi_{u_0}^{\prime\prime}(1) = \lim\limits_{n\rightarrow \infty}\varphi_{u_n}^{\prime\prime}(1)\le 0. \end{equation} (2.17)

    From Lemma 2.7, it follows that ||u_n^\pm|| > \sigma > 0 for any u_n\in\mathcal{M}_\lambda^- and hence ||u_0^\pm|| = \lim\limits_{n\rightarrow \infty}||u_n^\pm|| > \sigma > 0 , which indicates u_0^\pm\neq 0 . Using this and (2.14), we obtain u_0\in\mathcal{M}_\lambda and r = 1 is a critical point of \varphi_{u_0} . Consequently, by (2.15)-(2.17), Lemma 2.5 and Corollary 2.6, we derive that

    \frac{\partial ^2f_{u_0}}{\partial s^2}(1, 1) < 0, \; \frac{\partial ^2f_{u_0}}{\partial t^2}(1, 1) < 0, \; \varphi_{u_0}^{\prime\prime}(1) < 0.

    Hence, u_0\in\mathcal{M}_\lambda^* and \mathcal{M}_\lambda^* is a closed set.

    Lemma 2.13. If 1 < q < 2 , 4 < p < 6 , \lambda < \lambda^* and assumptions (V) and (K) hold, then the infimum m_\lambda: = \inf\limits_{u\in\mathcal{M}_\lambda^*}I_\lambda(u) can be achieved by some u_\lambda\in\mathcal{M}_\lambda^* and m_\lambda > 0 .

    Proof. According to Lemma 2.5, m_\lambda > -\infty when \lambda < \lambda^* . Let \{u_n\}\subset\mathcal{M}_\lambda^* be a minimizing sequence for the functional I_\lambda , namely I_\lambda(u_n)\rightarrow m_\lambda as n\rightarrow \infty . Since the functional I_\lambda is coercive on \mathcal{M}_\lambda^* , then \{u_n\} is bounded in H . Going if necessary to a subsequence, we may assume that

    u_n\rightharpoonup u_\lambda\; \rm in\; \textit H, \; \textit u_\textit n\rightarrow \textit u_\lambda\; in\; \textit L^{\textit p}(\mathbb{R}^{3}).

    Now, we first claim that u_\lambda^\pm\neq 0 . In fact, from Lemma 2.2, Lemma 2.7 and the convergence of \{u_n\} in L^p(\mathbb{R}^{3}) , for any \lambda < \lambda^* , we conclude that

    \begin{equation*} \begin{split} |u_\lambda^\pm|_p^p+\lambda\int_{\mathbb{R}^{3}}K(x)|u^\pm_\lambda|^q\mathrm{d}x& = \lim\limits_{n\rightarrow \infty}\Big(|u_n^\pm|_p^p+\lambda\int_{\mathbb{R}^{3}}K(x)|u^\pm_n|^q\mathrm{d}x\Big)\\& = \lim\limits_{n\rightarrow \infty}\Big(||u_n^\pm||^2+\int_{\mathbb{R}^{3}}\phi_{u_n^\pm}{(u_n^\pm)^2}\mathrm{d}x+\int_{\mathbb{R}^{3}}\phi_{u_n^+}{(u_n^-)^2}\mathrm{d}x\Big)\\&\ge \lim\limits_{n\rightarrow \infty}{||u_n^\pm||^2} > \sigma^2 > 0. \end{split} \end{equation*}

    This means that u_\lambda^\pm\neq 0 for any \lambda < \lambda^* .

    Next, we proof that u_n\rightarrow u_\lambda in H . Arguing by contradiction, suppose that

    ||u_\lambda^+|| < \lim\limits_{n\rightarrow \infty}\inf||u_n^+||\; \rm or\; ||\textit u_ \lambda^-|| < \lim\limits_{\textit n\rightarrow \infty}\inf||\textit u_\textit n^-||.

    From Corollary 2.10, there exists a unique pair (s_{u_\lambda}, t_{u_\lambda}) such that \tilde{u}_\lambda = s_{u_\lambda}u_\lambda^++t_{u_\lambda}u_\lambda^-\in\mathcal{M}_\lambda^* and I_\lambda(u_n^++u_n^-) = \max\limits_{s, t > 0}I_\lambda(s_{u_\lambda}u_n^++t_{u_\lambda}u_n^-) . Consequently,

    m_\lambda\leq I_\lambda(\tilde{u}_\lambda) < \lim\limits_{n\rightarrow \infty}\inf I_\lambda(s_{u_\lambda}u_n^++t_{u_\lambda}u_n^-) \leq\lim\limits_{n\rightarrow \infty}\inf I_\lambda(u_n^++u_n^-) = m_\lambda.

    That is, we get a contradiction. Therefore, u_n\rightarrow u_\lambda in H and m_\lambda is achieved by u_\lambda . Combining the fact that \mathcal{M}_\lambda^* is closed, so u_\lambda\in\mathcal{M}_\lambda^* .

    Finally, it follows from u_\lambda\in\mathcal{M}_\lambda^* that \varphi_{u_\lambda}(r) reached its global maximum at r = 1 . By this and Lemma 2.11, we can deduce that \varphi_{u_\lambda}(1) > 0 , i.e., m_\lambda > 0 . This finishes the proof of Lemma 2.13.

    The main aim of this section is to prove our results. Thanks to Lemma 2.13, it suffices to check that the minimizer u_\lambda for m_\lambda is a sign-changing of problem (1.3).

    Proof of Theorem 1.1. Since u_\lambda\in\mathcal{M}_\lambda^* , according to Corollary 2.10, we obtain that

    \begin{equation*} I_\lambda(su_\lambda^++tu_\lambda^-) < I_\lambda(u_\lambda^++u_\lambda^-) = m_\lambda, \; \rm for\; (\textit s, \textit t)\in (\mathbb{R}^+\times\mathbb{R}^+)\backslash {(1, 1)}. \end{equation*}

    Moreover, we get I_\lambda(u_\lambda) > 0 , \varphi_{u_\lambda}^{\prime\prime}(1, 1) < 0 , (\partial ^2f_{u_\lambda}/\partial s^2)(1, 1) < 0 and (\partial ^2f_{u_\lambda}/\partial t^2)(1, 1) < 0 .

    Let D = (1-\delta, 1+\delta)\times(1-\delta, 1+\delta) and h: D\rightarrow H by h(s, t) = su_\lambda^++tu_\lambda^- for any (s, t)\in D . Then there is a constant 0 < \delta < 1 such that

    \begin{equation} 0 < m: = \max\limits_{\partial D}I_\lambda(h(s, t)) < m_\lambda, \; \; \; \max\limits_{(s, t)\in D}\frac{\partial ^2f_{h(s, t)}}{\partial s^2}(1, 1) < 0, \end{equation} (3.1)
    \begin{equation} \max\limits_{(s, t)\in D}\frac{\partial ^2f_{h(s, t)}}{\partial t^2}(1, 1) < 0, \; \; \; \max\limits_{(s, t)\in D}\varphi_{h(s, t)}^{\prime\prime}(1) < 0. \end{equation} (3.2)

    By the quantitative deformation lemma, we prove that I_\lambda^\prime(u_\lambda) = 0 . Suppose by contradiction that I_\lambda^\prime(u_\lambda)\neq 0 , then there exist \lambda_1 > 0 and \xi > 0 such that

    ||I_\lambda^\prime(v)||\ge \lambda_1\; \rm for\; all\; \textit v\in\textit H, \; ||\textit v-\textit u_\lambda||\le 3\xi.

    Let \varepsilon = \min\{\frac{m_\lambda-m}{3}, \frac{\lambda_1\xi}{8}\} and s_\xi = \{u\in H:\; ||u-u_\lambda||\le \xi\} , then the deformation lemma (see[32], Lemma 2.3) shows that there is a deformation \eta\in C([0, 1]\times H, H) such that

    (i) \eta(d, u) = u if u\not \in I_\lambda^-([m_\lambda-2\varepsilon, m_\lambda+2\varepsilon])\cap s_{2\xi}, \; d\in[0, 1] ;

    (ii) I_\lambda(\eta(d, u))\le I_\lambda(u) for all u\in H , d\in[0, 1] ;

    (iii) I_\lambda(\eta(d, u)) < m_\lambda , \forall u\in I_\lambda^{m_\lambda}\cap s_\xi , \forall d\in(0, 1] .

    First, we need to prove that

    \begin{equation} \max\limits_{(s, t)\in D}I_\lambda(\eta(d, h(s, t))) < m_\lambda \; \rm for\; all \; \textit d\in (0, 1]. \end{equation} (3.3)

    In fact, for any d\in(0, 1] , it follows from Corollary 2.10 and (ii) that

    \max\limits_{\{(s, t)\in D:\; h(s, t)\not\in s_{\xi}\}}I_\lambda(\eta(d, h(s, t)))\le \max\limits_{\{(s, t)\in D:\; h(s, t)\not\in s_{\xi}\}}I_\lambda(h(s, t)) < m_\lambda.

    Moreover, Corollary 2.10 and (iii) imply that

    \max\limits_{\{(s, t)\in D:\; h(s, t)\in s_{\xi}\}}I_\lambda(\eta(d, h(s, t))) < m_\lambda\; \rm for\; all \; \textit d\in (0, 1].

    Hence, (3.3) holds. From the continuity of \eta and (3.1)-(3.2), there exists a constant d_0\in(0, 1] such that

    \begin{equation} \begin{split} &\max\limits_{(s, t)\in D}\frac{\partial ^2f_{\eta(d_0, h(s, t))}}{\partial s^2}(1, 1) < 0, \\& \max\limits_{(s, t)\in D}\frac{\partial ^2f_{\eta(d_0, h(s, t))}}{\partial t^2}(1, 1) < 0, \\&\max\limits_{(s, t)\in D}\varphi^{\prime\prime}_{\eta_{(d_0, h(s, t))}}(1) < 0. \end{split} \end{equation} (3.4)

    In the following, we prove that \eta(d_0, h(D))\cap\mathcal{M}_\lambda^*\neq\emptyset , which contradicts the definition of m_\lambda . In fact, let g(s, t) = \eta(d_0, h(s, t)) and

    \psi_1(s, t) = \left(\langle I_\lambda^\prime(h(s, t)), u_\lambda^+\rangle, \langle I_\lambda^\prime(h(s, t)), u_\lambda^-\rangle\right),
    \psi_2(s, t) = \left(\frac{1}{s}\langle I_\lambda^\prime(g(s, t)), g^+(s, t)\rangle, \frac{1}{t}\langle I_\lambda^\prime(g(s, t)), g^-(s, t)\rangle\right).

    Then Corollary 2.10 and the degree theory yield deg (\psi_1, D, 0) = 1 . On the other hand, we know \varepsilon < \frac{m_\lambda-m}{3} , m < m_\lambda-2\varepsilon . Hence, from (i) we have \eta(d, h(s, t)) = h(s, t) for d\in(0, 1] , (s, t)\in\partial D , and it follows that

    \psi_1(s, t) = \psi_2(s, t)\; {\rm for\; any}\; ( s, t)\in\partial \textit D.

    Combining the homotopy invariance property of the degree, we get deg (\psi_2, D, 0) = deg (\psi_1, D, 0) = 1 . That is, there exists (s_0, t_0)\in D such that \psi_2(s_0, t_0) = 0 . Therefore, using (3.4) and \psi_2(s_0, t_0) = 0 , we have \eta(d_0, h(s_0, t_0))\in\mathcal{M}_\lambda^* , i.e., \eta(d_0, h(D))\cap\mathcal{M}_\lambda^*\neq\emptyset . From this, u_\lambda is a critical point of I_\lambda , i.e., I_\lambda^\prime(u_\lambda) = 0 .

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This research is supported by National Natural Science Foundation of China [No.11971393].

    The authors declare there is no conflict of interest.



    [1] United Nations Framework Convention on Climate Change (UNFCCC) (2016) Paris Agreement. Available from: https://unfccc.int/sites/default/files/english_paris_agreement.pdf.
    [2] Intergovernmental Panel on Climate Change (IPCC) (2015) Climate change 2014 mitigation of climate change: Summary for policymakers and technical summary. Available from: https://www.ipcc.ch/site/assets/uploads/2018/03/WGIIIAR5_SPM_TS_Volume-3.pdf.
    [3] Kawakami A (2020) Suga's 2050 zero-carbon goal thrusts Japan into green tech race. NIKKEI Asia. Available from: https://asia.nikkei.com/Spotlight/Environment/Suga-s-2050-zero-carbon-goal-thrusts-Japan-into-green-tech-race.
    [4] Ministry of Economy, Trade and Industry (METI) (2021) Outline of 6th Strategic Energy Plan. Available from: https://www.enecho.meti.go.jp/en/category/others/basic_plan/pdf/6th_outline.pdf.
    [5] Agency for Natural Resources and Energy (ANRE) (2018) Japan's strategic energy plan. Available from: https://www.numo.or.jp/topics/1-1Nakanishi.pdf.
    [6] Ministry of Economy, Trade and Industry (METI) (2021) Japan's energy 2020—10 questions for understanding the current energy situation. Available from: https://www.enecho.meti.go.jp/en/category/brochures/pdf/japan_energy_2020.pdf.
    [7] Kato J (2020) Initiatives to make offshore wind power the primary power source. JWPA. Available from: https://reglobal.co/wp-content/uploads/2021/01/Offshore-Wind-in-Japan.pdf.
    [8] Ministry of Economy, Trade and Industry (METI)(2020) Methods of the METI for the realization of a hydrogen society. Available from: https://www.env.go.jp/seisaku/list/ondanka_saisei/lowcarbon-h2-sc/events/PDF/shiryou06.pdf.
    [9] Hokkaido Government (2020) Hokkaido hydrogen society realization strategy vision. Available from: https://www.pref.hokkaido.lg.jp/fs/5/6/0/1/2/5/1/_/suisovisionkaiteiban_soan.pdf.
    [10] METI Kyushu (2022) Kyushu hydrogen guidebook. Available from: https://www.kyushu.meti.go.jp/seisaku/kankyo/jirei/2022_hydrogen/pdf/1_all.pdf.
    [11] Kimura K (2018) Analysis of wind power costs in Japan. Renewable Energy Institute (Renewable EI). Available from: https://www.renewable-ei.org/en/activities/reports/img/pdf/20180125/JapanWindPowerCostReport_EN_20180124.pdf.
    [12] Ishihara T, Yamaguchi A, Takahara K, et al. (2005) An analysis of damaged wind turbines by typhoon Maemi in 2003. The Sixth Asia-Pacific Conference on Wind Engineering. Available from: file: ///C: /Users/User/Downloads/Turbine-Damage-Maemi-Typhone-Analysis. pdf.
    [13] Chen X, Xu JZ (2016) Structural failure analysis of wind turbines impacted by super typhoon Usagi. Engineering Failure Analysis. 60: 391–404.https://doi.org/10.1016/j.engfailanal.2015.11.028
    [14] Siemens Gamesa offshore. Available from: https://www.siemensgamesa.com/products-and-services/offshore.
    [15] Vestas V174-9.5 MW turbine. Available from: https://www.vestas.com/en/products/offshore/v174-9-5-mw-.
    [16] Vestas V236-15 MW turbine. Available from: https://www.vestas.com/en/products/offshore/V236-15MW.
    [17] Durakovic A (2021) MingYang launches 16 MW offshore wind turbine. Available from: https://www.offshorewind.biz/2021/08/20/mingyang-launches-16-mw-offshore-wind-turbine/.
    [18] IECRE-Renewable Energy (2022) Type Certificate RNA—Rotor Nacelle Assembly. Adapted to offshore wind conditions, issued to General Electric Renovables España. Available from: https://certificates.iecre.org/#/deliverables/CERT/1587187/view.
    [19] Energy Northern Perspective (2021) Siemens Gamesa to install first of its SG 11.0-200 DD machines at Ørsted's Gode Wind 3. Available from: https://energynorthern.com/2021/12/20/siemens-gamesa-to-install-first-of-its-sg-11-0-200-dd-machines-at-orsteds-gode-wind-3/.
    [20] Largue P (2020) MHI Vestas installs flagship V174-9.5 MW prototype. Power Eng Int. Available from: https://www.powerengineeringint.com/renewables/mhi-vestas-installs-flagship-v174-9-5-mw-prototype/.
    [21] Durakovic A (2021) SG 11.0-200 DD turbine model secures full type certificate early. offshoreWIND. biz. Available from: https://www.offshorewind.biz/2021/05/31/sg-11-0-200-dd-turbine-model-secures-full-type-certificate-early/.
    [22] Durakovic A (2021) Vestas 15 MW prototype offshore wind turbine to spin in Denmark. offshoreWIND. biz. Available from: https://www.offshorewind.biz/2021/10/15/vestas-15-mw-prototype-offshore-wind-turbine-to-spin-in-denmark/.
    [23] Nehls G (2021) MingYang Smart Energy launches MySE 16.0-242 offshore hybrid drive wind turbine. CompositesWorld. Available from: https://www.compositesworld.com/news/mingyang-smart-energy-launches-myse-160-242-offshore-hybrid-drive-wind-turbine.
    [24] Energy Rich Japan Report (2003) Available from: www.energyrichjapan.info.
    [25] Esteban M, Zhang Q, Utama A, et al. (2010) Methodology to estimate the output of a dual solar-wind renewable energy system in Japan. Energy Policy 38: 7793–7802.https://doi.org/10.1016/j.enpol.2010.08.039 doi: 10.1016/j.enpol.2010.08.039
    [26] Esteban M, Zhang Q, Utama A (2012) Estimation of the energy storage requirement of a future 100% renewable energy system in Japan. Energy Policy 47: 22–31.https://doi.org/10.1016/j.enpol.2012.03.078 doi: 10.1016/j.enpol.2012.03.078
    [27] Esteban M, Portugal Pereira J (2014) Post-disaster resilience of a 100% renewable energy system in Japan. Energy 68: 756–764.https://doi.org/10.1016/j.energy.2014.02.045 doi: 10.1016/j.energy.2014.02.045
    [28] Esteban M, Portugal Pereira J, Mclellan BC, et al. (2018) 100% renewable energy system in Japan: Smoothening and ancillary services. Appl Energy 224: 698–707.https://doi.org/10.1016/j.apenergy.2018.04.067 doi: 10.1016/j.apenergy.2018.04.067
    [29] Wakiyama T, Kuriyama A (2018) Assessment of renewable energy expansion potential and its implications on reforming Japan's electricity system. Energy Policy 115: 302–316.https://doi.org/10.1016/j.enpol.2018.01.024. doi: 10.1016/j.enpol.2018.01.024
    [30] Yoshihara K, Ohashi H (2017) Assessing the impact of renewable energy sources: Simulation analysis of the Japanese electricity market. Research Institute of Economy, Trade and Industry. Discussion Paper Series 17-E-063. Available from: https://www.rieti.go.jp/en/publications/summary/17040014.html.
    [31] Knüpfer K, Rogalski N, Knüpfer A, et al. (2022) A reliable energy system for Japan with merit order dispatch, high variable renewable share and no nuclear power. Appl Energy 328.https://doi.org/10.1016/j.apenergy.2022.119840
    [32] Tsuchiya H (2012) Electricity supply largely from solar and wind resources in Japan. Renewable Energy 48: 318–325.https://doi.org/10.1016/j.renene.2012.05.011 doi: 10.1016/j.renene.2012.05.011
    [33] Heard BP, Brook BW, Wigley TML, et al. (2017) Burden of proof: A comprehensive review of feasibility of 100% renewable-electricity systems. Renewable Sustainable Energy Rev 76: 1122–1133.https://doi.org/10.1016/j.rser.2017.03.114 doi: 10.1016/j.rser.2017.03.114
    [34] Ministry of Economy, Trade and Industry (METI) (2021) About renewable energy in 2030. Available from: https://www.meti.go.jp/shingikai/enecho/denryoku_gas/saisei_kano/pdf/031_02_00.pdf.
    [35] Organization for Cross-regional Coordination of Transmission Operators (OCCTO) (2019) Aggregation of electricity supply plans, Fiscal Year 2019, Figure 3–4: Composition of installed power generation capacity for each regional service area. Available from: http://www.occto.or.jp/en/information_disclosure/supply_plan/files/supplyplan_2019.pdf.
    [36] Agency for Natural Resources and Energy (ANRE) (2022) Website of the introduction status of renewable energy electricity. Available from: https://www.fit-portal.go.jp/PublicInfoSummary.
    [37] Japan Atomic Industrial Forum (JAIF) (2019) Current status of nuclear power plants in Japan. Available from: https://www.jaif.or.jp/cms_admin/wp-content/uploads/2019/12/jp-npps-operation20191205_en.pdf.
    [38] Knüpfer K, Dumlao SMG, Esteban M, et al. (2021) Analysis of PV subsidy schemes, installed capacity and actual generation in Japan from 1994 to 2019. Energies 8: 2128.https://doi.org/10.3390/en14082128 doi: 10.3390/en14082128
    [39] Ministry of Economy, Trade and Industry (METI) (2020) Overview of the vision for offshore wind power industry. Available from: https://www.meti.go.jp/shingikai/energy_environment/yojo_furyoku/pdf/002_02_e01_01.pdf.
    [40] Department of Energy (DOE) (2020) Global Energy storage database. Available from: https://sandia.gov/ess-ssl/gesdb/public/index.html.
    [41] Japan Meteorological Agency (JMA) (2020) Past meteorological data download landing page. Available from: https://www.data.jma.go.jp/gmd/risk/obsdl/index.php.
    [42] Japan Meteorological Agency (JMA) (2021) List of regional meteorological measurement station. Available from: https://www.jma-net.go.jp/common/catalogue/format/ObdObs_Amedas_ObsList_format.pdf.
    [43] Dhar A, Naeth MA, Jennings PD, et al. (2020) Geothermal energy resources: potential environmental impact and land reclamation. Environ Rev 28: 415–427.https://doi.org/10.1139/er-2019-0069 doi: 10.1139/er-2019-0069
    [44] USACE, USFWS (2008) Hydropower rehabilitations, dissolved oxygen and minimum flow at Wolf Creek Dam, Kentucky and Center hill and dale hollow dams, Tennessee. Draft Environmental Impact Assessment, 13. Available from: https://books.google.co.jp/books?id = nQMzAQAAMAAJ & printsec = frontcover & hl = ja#v = onepage & q & f = false.
    [45] Agency for Natural Resources and Energy (ANRE) (2020) Power survey statistics of 2020. Available from: https://www.enecho.meti.go.jp/statistics/electric_power/ep002/pdf/2020/0-2020.pdf.
    [46] TEPCO (2020) Hourly demand progress. Available from: https://www.tepco.co.jp/en/forecast/html/download-e.html.
    [47] Kainou K (2016) Development of policy impact assessment model for the regulatory reform policy of Japanese electricity market. Research Institute of Economy, Trade and Industry. Discussion Paper Series 16-J-012. Available from: https://www.rieti.go.jp/jp/publications/dp/16j012.pdf.
    [48] Japan Atomic Industrial Forum (JAIF) (2022) Current status of nuclear power plants in Japan. Available from: https://www.jaif.or.jp/cms_admin/wp-content/uploads/2022/04/2022-02.pdf.
    [49] JPOWER (2022) Geothermal Power Generation Business. Available from: https://www.jpower.co.jp/bs/renewable_energy/geothermal/.
    [50] Japan Organization for Metals and Energy Security (JOGMEC) (2019) Geothermalin Japan. Available from: https://geothermal.jogmec.go.jp/information/plant_japan/.
  • This article has been cited by:

    1. Xueqi Sun, Yongqiang Fu, Sihua Liang, Normalized solutions for pseudo-relativistic Schrödinger equations, 2024, 16, 2836-3310, 217, 10.3934/cam.2024010
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2700) PDF downloads(150) Cited by(5)

Figures and Tables

Figures(7)  /  Tables(17)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog