Loading [MathJax]/jax/output/SVG/jax.js
Research article

Dynamic study of the pathogen-immune system interaction with natural delaying effects and protein therapy

  • Received: 19 September 2021 Revised: 30 November 2021 Accepted: 03 December 2021 Published: 14 February 2022
  • MSC : 92F05, 34K20, 34K18

  • This study aims to propose and analyze a mathematical model of the competitive interaction of the pathogen-immune system. Some effects of the existence of natural delays and the addition of therapeutic proteins are considered in the model. A delay arises from the indirect response of the host body when a pathogen invades. The other comes from the maturation of immune cells to produce immune memory cells since the immune system and antigenic substances responsible for provoking the production of immune memory cells. Analytical investigations suggest several sufficient conditions for the existence of a positive steady-state solution. There is a critical pair of delays at which oscillatory behavior appears around the positive steady-state solution. Numerical simulations were carried out to describe the results of the analysis and show that the proposed model can describe the speed of pathogen eradication due to the addition of therapeutic proteins as antigenic substances.

    Citation: Kasbawati, Yuliana Jao, Nur Erawaty. Dynamic study of the pathogen-immune system interaction with natural delaying effects and protein therapy[J]. AIMS Mathematics, 2022, 7(5): 7471-7488. doi: 10.3934/math.2022419

    Related Papers:

    [1] Lovelesh Sharma . Brezis Nirenberg type results for local non-local problems under mixed boundary conditions. Communications in Analysis and Mechanics, 2024, 16(4): 872-895. doi: 10.3934/cam.2024038
    [2] Shu Wang . Global well-posedness and viscosity vanishing limit of a new initial-boundary value problem on two/three-dimensional incompressible Navier-Stokes equations and/or Boussinesq equations. Communications in Analysis and Mechanics, 2025, 17(2): 582-605. doi: 10.3934/cam.2025023
    [3] Yuxuan Chen . Global dynamical behavior of solutions for finite degenerate fourth-order parabolic equations with mean curvature nonlinearity. Communications in Analysis and Mechanics, 2023, 15(4): 658-694. doi: 10.3934/cam.2023033
    [4] Eleonora Amoroso, Angela Sciammetta, Patrick Winkert . Anisotropic (p,q)-Laplacian problems with superlinear nonlinearities. Communications in Analysis and Mechanics, 2024, 16(1): 1-23. doi: 10.3934/cam.2024001
    [5] Leandro Tavares . Solutions for a class of problems driven by an anisotropic (p,q)-Laplacian type operator. Communications in Analysis and Mechanics, 2023, 15(3): 533-550. doi: 10.3934/cam.2023026
    [6] Efstratios Stratoglou, Alexandre Anahory Simoes, Leonardo J. Colombo . Reduction in optimal control with broken symmetry for collision and obstacle avoidance of multi-agent system on Lie groups. Communications in Analysis and Mechanics, 2023, 15(2): 1-23. doi: 10.3934/cam.2023001
    [7] Farrukh Dekhkonov . On one boundary control problem for a pseudo-parabolic equation in a two-dimensional domain. Communications in Analysis and Mechanics, 2025, 17(1): 1-14. doi: 10.3934/cam.2025001
    [8] Xiao Han, Hui Wei . Multiplicity of the large periodic solutions to a super-linear wave equation with general variable coefficient. Communications in Analysis and Mechanics, 2024, 16(2): 278-292. doi: 10.3934/cam.2024013
    [9] Huiyang Xu . Existence and blow-up of solutions for finitely degenerate semilinear parabolic equations with singular potentials. Communications in Analysis and Mechanics, 2023, 15(2): 132-161. doi: 10.3934/cam.2023008
    [10] Yang Liu . Global attractors for a nonlinear plate equation modeling the oscillations of suspension bridges. Communications in Analysis and Mechanics, 2023, 15(3): 436-456. doi: 10.3934/cam.2023021
  • This study aims to propose and analyze a mathematical model of the competitive interaction of the pathogen-immune system. Some effects of the existence of natural delays and the addition of therapeutic proteins are considered in the model. A delay arises from the indirect response of the host body when a pathogen invades. The other comes from the maturation of immune cells to produce immune memory cells since the immune system and antigenic substances responsible for provoking the production of immune memory cells. Analytical investigations suggest several sufficient conditions for the existence of a positive steady-state solution. There is a critical pair of delays at which oscillatory behavior appears around the positive steady-state solution. Numerical simulations were carried out to describe the results of the analysis and show that the proposed model can describe the speed of pathogen eradication due to the addition of therapeutic proteins as antigenic substances.



    The study of solutions to superlinear problems driven by the double phase operator is a new and important topic, since it sheds light on a range of applications in the field of mathematical physics such as elasticity theory, strongly anisotropic materials, Lavrentiev's phenomenon, etc. (see [1,2,3]).

    In the present paper, we study the existence and multiplicity of solutions for the following double–phase problems with mixed boundary conditions:

    {D(u)+|u|p(x)2u+b(x)|u|q(x)2u=f(x,u)   in Ω,u=0  on Λ1,(|u|p(x)2u+b(x)|u|q(x)2u)ν=g(x,u)   on Λ2, (1.1)

    where ΩRN(N2) is a bounded domain with Lipschitz boundary Ω, Λ1,Λ2 are disjoint open subsets of Ω such that Ω=¯Λ1¯Λ2 and Λ1, 1<p(x)<q(x)<N for all x¯Ω, b:¯Ω[0,+) is Lipschitz continuous, ν denotes the outer unit normal of Ω at the point xΛ2, f:Ω×RR and g:Λ2×RR are Carathéodory functions, and D is the double phase variable exponents operator given by

    D(u):=div(|u|p(x)2u+b(x)|u|q(x)2u),  for  uW1,D(Ω). (1.2)

    Note that the differential operator defined above in (1.2) is called the double phase operator with variable exponents, which is a natural generalization of the classical double phase operator when p and q are constant functions

    udiv(|u|p2u+b(x)|u|q2u).

    From the physical point of view, while studying the behavior of strongly anisotropic materials, V.V. Zhikov [2] in 1986 discovered that their hardening properties changed radically point by point, what is known as the Lavrentiev phenomenon [3]. To describe this phenomenon, he initially introduced the functional

    vΩ(|v|p+b(x)|v|q)dx, (1.3)

    where the integrand changes its ellipticity and growth properties according to the point in the domain. In the framework of mathematics, the functional (1.3) has been investigated by many authors with respect to regularity and nonstandard growth. For instance, we refer to the papers of P. Baroni et al. [4,5], P. Baroni et al. [6], G. Cupini [7], and the references therein.

    Multiple authors have recently concentrated on the study of double phase problems in the case when the exponents p and q are constants, and a plethora of results have been obtained; see, for example, W. Liu and G. Dai [8], M. El Ahmadi et al. [9], L. Gasiński and P. Winkert [10], N. Cui and H.R. Sun [11], Y. Yang et al. [12], and the references therein. For example, N. Cui and H.R. Sun [11] considered the following problem in the particular case: p(x)=p, q(x)=q, and λ=1

    {D(u)+|u|p2u+b(x)|u|q2u=f(x,u)     in Ω,(|u|p2u+b(x)|u|q2u)ν=g(x,u)  on Ω,

    where D(u):=div(|u|p2u+b(x)|u|q2u). The authors have proved the existence and multiplicity of nontrivial weak solutions for the above problem with superlinear nonlinearity. Their approach was based on critical point theory with Cerami condition.

    Very recently, Y. Yang et al. [12] considered the problem (1.1) in the particular case of p(x)=p and q(x)=q. Based on the maximum principle and homological local linking, they proved the existence of at least two bounded nontrivial weak solutions.

    The main novelty of the current paper is the combination of the double phase variable exponents operator with mixed boundary conditions, that is, the Dirichlet condition on Λ1 and the Steklov condition on Λ2, which is different from [13]. To the best of our knowledge, there are only a few results related to the study of such problems.

    To state our results, we make the subsequent hypotheses on f and g:

    (H0) There exist C1,C2>0, s1C+(Ω), and s2C+(Λ2) such that

    (ⅰ) |f(x,t)|C1(1+|t|s1(x)1)   for all (x,t)Ω×R,

    (ⅱ) |g(x,t)|C2(1+|t|s2(x)1)   for all (x,t)Λ2×R.

    (H1) (ⅰ) 1<p+q+<s1s+1<p(x) for all  xΩ, $

    (ⅱ) 1<p+q+<s2s+2<p(x) for all  xΛ2,

    where

    p(x):={Np(x)Np(x)       if p(x)<N,   if p(x)N,p(x):={(N1)p(x)Np(x)      if p(x)<N,     if p(x)N..

    (H2) (ⅰ) lim inf|t|F(x,t)|t|q+=+ uniformly a.e. xΩ,

    (ⅱ) lim inf|t|G(x,t)|t|q+=+ uniformly a.e. xΛ2,

    where F(x,t)=t0f(x,s)ds and G(x,t)=t0g(x,s)ds.

    (H3) (ⅰ) There exist c1,r10 and l1L(Ω) with l1(x)>Np such as

    |F(x,t)|l(x)c1|t|l(x)pF(x,t),

    for all (x,t)Ω×R, |t|r1 and F(x,t):=1q+f(x,t)tF(x,t)0.

    (ⅱ) There exist c2,r20 and l2L(Λ2) with l2(x)>N1p1 such as

    |G(x,t)|l(x)c2|t|l(x)pG(x,t),

    for all (x,t)Λ2×R, |t|r2 and G(x,t):=1q+g(x,t)tG(x,t)0.

    (H4) (ⅰ) f(x,t)=(|t|p+1) as t0 uniformly for a.e. xΩ,

    (ⅱ) g(x,t)=(|t|p+1) as t0 uniformly for a.e. xΛ2.

    (H5) (ⅰ) f(x,t)=f(x,t) for all (x,t)Ω×R,

    (ⅱ) g(x,t)=g(x,t) for all (x,t)Λ2×R.

    Let us consider ϕ:X0R the Euler functional corresponding to problem (1.1), which is defined as follows:

    ϕ(u)=I(u)φ(u),

    where

    I(u)=Ω(1p(x)|u|p(x)+b(x)q(x)|u|q(x))dx+Ω(1p(x)|u|p(x)+b(x)q(x)|u|q(x))dx,

    and

    φ(u)=ΩF(x,u)dx+Λ2G(x,u)dσ,

    with X0 will be defined in preliminaries and dσ is the measure on the boundary.

    Then, it follows from the hypothesis (H0) that the functional ϕC1(X0,R), and its Fréchet derivative is

    ϕ(u),v=Ω(|u|p(x)2+b(x)|u|q(x)2)u.v dx+Ω(|u|p(x)2+b(x)|u|q(x)2)u.v dxΩf(x,u)vdxΛ2g(x,u)vdσ,

    for any u,vX0. It is clear that any critical point of ϕ is a weak solution to the problem (1.1).

    Now, we present the main results of this paper.

    Theorem 1. Suppose that (H0),(H1),(H2),(H3) and (H4) hold. Then problem (1.1) has at least one nontrivial weak solution.

    Theorem 2. Suppose that (H0),(H1),(H2),(H3) and (H5) hold. Then problem (1.1) possesses a sequence of weak solutions (un) such that ϕ(un)+ as n+.

    To study double phase problems, we need some definitions and basic properties of W1,D(Ω), which are called Musielak–Orlicz–Sobolev spaces. For more details, see [14,15,16,17,18,19] and references therein.

    First, we recall the definition of variable exponent Lebesgue space. For pC+(¯Ω):={pC(¯Ω):p:=infx¯Ωp(x)>1}, we designate the variable exponent Lebesgue space by

    Lp(x)(Ω)={u:ΩR is measurable and Ω|u(x)|p(x)dx<+},

    equipped with the Luxemburg norm

    |u|p(x)=inf{μ>0:Ω|u(x)μ|p(x)dx1}.

    Proposition 1. [20]

    1. The Sobolev space (Lp(x)(Ω),|.|p(x)) is defined as the dual space Lq(x)(Ω), where q(x) is conjugate to p(x), i.e., 1p(x)+1q(x)=1. For any uLp(x)(Ω) and vLq(x)(Ω), we have

    |Ωuv dx|(1p+1q)|u|p(x)|v|q(x)2|u|p(x)|v|q(x).

    2. If p1,p2C+(¯Ω), p1(x)p2(x), for all x¯Ω, then Lp2(x)(Ω)Lp1(x)(Ω) and the embedding is continuous.

    Let pC+(Ω):={pC(Ω):p:=infxΩp(x)>1} and denote by dσ the Lebesgue measure on the boundary. We define

    Lp(x)(Ω)={u:ΩR is measurable and Ω|u|p(x)dσ<+},

    with the norm

    |u|p(x),Ω=inf{α>0:Ω|u(x)α|q(x)dσ1}.

    Now, we give the main properties of the Musielak–Orlicz–Sobolev functional space that we will use in the rest of this paper. Denote by N(Ω) the set of all generalized N-functions. Let us denote by

    D:Ω×[0,+)[0,+)

    the function defined as

    D(x,t)=tp(x)+b(x)tq(x), for all  (x,t)Ω×[0,+),

    where the weight function b(.) and the variable exponents p(.),q(.)C+(¯Ω) satisfies the following hypothesis:

    p(x)<q(x)<N,Nq(x)N+q(x)1<p(x)  for all x¯Ω  and  0b(.)L(Ω). (2.1)

    Note that the role of assuming the inequality Nq(x)N+q(x)1<p(x) is to ensure that q(x)<p(x) and q(x)<p(x) for all x¯Ω, where p(x)=Np(x)Np(x) and p(x)=(N1)p(x)Np(x).

    It is clear that D is a generalized N-function, locally integrable, and

    D(x,2t)2q+D(x,t), for all  (x,t)Ω×[0,+),

    which is called condition (Δ2).

    We designate the Musielak–Orlicz space by

    LD(Ω)={u:ΩR is measurable and ΩD(x,|u|)dx<+},

    equipped with the so-called Luxemburg norm

    |u|D=inf{μ>0:ΩD(x,|uμ|)dx1}.

    The Musielak–Orlicz–Sobolev space W1,D(Ω) is defined as

    W1,D(Ω)={uLD(Ω):|u|LD(Ω)},

    endowed with the norm

    u1,D=|u|D+|u|D.

    With such norms, LD(Ω) and W1,D(Ω) are separable, uniformly convex, and reflexive Banach spaces.

    On LD(Ω), we consider the function ρ:LD(Ω)R defined by

    ρ(u)=Ω(|u|p(x)+b(x)|u|q(x))dx.

    The relationship between ρ and |.|D is established by the next result.

    Proposition 2. (see [16]) For uLD(Ω), (un)LD(Ω), and μ>0, we have

    1. For u0, |u|D=μρ(uμ)=1;

    2. |u|D<1(=1,>1)ρ(u)<1(=1,>1);

    3. |u|D>1|u|pDρ(u)|u|q+D;

    4. |u|D<1|u|q+Dρ(u)|u|pD;

    5. limn+|un|D=0limn+ρ(un)=0 and limn+|un|D=+limn+ρ(un)=+.

    On W1,D(Ω), we introduce the equivalent norm by

    u:=inf{μ>0:Ω[|uλ|p(x)+b(x)|uλ|q(x)+|uλ|p(x)+b(x)|uμ|q(x)]dx1}. (2.2)

    Similar to Proposition (2), we have

    Proposition 3. (see [16]) Let

    ˆρ(u)=Ω(|u|p(x)+b(x)|u|q(x))dx+Ω(|u|p(x)+b(x)|u|q(x))dx.

    For uW1,D(Ω), (un)W1,D(Ω), and μ>0, we have

    1. For u0, u=μˆρ(uμ)=1;

    2. u<1(=1,>1)ˆρ(u)<1(=1,>1);

    3. u>1upˆρ(u)uq+;

    4. u<1uq+ˆρ(u)up;

    5. limn+un=0limn+ˆρ(un)=0 and limn+un=+limn+ˆρ(un)=+.

    We recall that problem (1.1) has a mixed boundary condition. For this, our Banach space workspace is given by

    X0:={uW1,D(Ω):u=0  on  Λ2},

    endowed with the equivalent norm (2.2). Obviously, since X0 is a closed subspace of W1,D(Ω), then (X0,.) is a reflexive Banach space.

    Proposition 4. (see [16]) Let hypothesis (2.1) be satisfied. Then the following embeddings hold:

    1. There is a continuous embedding LD(Ω)Lr(x)(Ω) for rC(¯Ω) with 1r(x)p(x) for all x¯Ω.

    2. There is a compact embedding W1,D(Ω)Lr(x)(Ω) for rC(¯Ω) with 1r(x)<p(x) for all x¯Ω.

    3. If pC+(¯Ω)W1,γ(Ω)for someγN. Then, there is a continuous embedding W1,D(Ω)Lr(x)(Ω) for rC(Ω) with 1r(x)p(x) for all xΩ.

    4. There is a compact embedding W1,D(Ω)Lr(x)(Ω) for rC(Ω) with 1r(x)<p(x) for all xΩ.

    It is important to note that when we replace W1,D(Ω) by X0 in Proposition 4, the embeddings 2 and 4 remain valid.

    Let A:W1,D(Ω)(W1,D(Ω)) be defined by

    A(u),v=Ω(|u|p(x)2+b(x)|u|q(x)2)u.vdx+Ω(|u|p(x)2+b(x)|u|q(x)2)u.vdx,

    for all u,vW1,D(Ω), where (W1,D(Ω)) denotes the dual space of W1,D(Ω) and .,. stands for the duality pairing between W1,D(Ω) and (W1,D(Ω)).

    Proposition 5. (see [16,Proposition 3.4]) Let hypothesis (2.1) be satisfied.

    1. A:W1,D(Ω)(W1,D(Ω)) is a continuous, bounded, and strictly monotone operator.

    2. A:W1,D(Ω)(W1,D(Ω)) satisfies the (S+)-property, i.e., if unu in W1,D(Ω) and ¯limn+A(un)A(u),unu0, then unu in W1,D(Ω).

    Definition 1. Let uX0. We say that u is a weak solution to the problem (1.1) if

    Ω(|u|p(x)2+b(x)|u|q(x)2)u.v dx+Ω(|u|p(x)2+b(x)|u|q(x)2)u.v dxΩf(x,u)vdxΛ2g(x,u)vdσ,=0,

    for all vX0.

    Now, we give the definition of the Cerami condition that was first introduced by G. Cerami in [21].

    Definition 2. Let (X,.) be a real Banach space and ϕC1(X,R). We say that ϕ satisfies the Cerami condition (we denote (C)condition) in X, if any sequence (un)X such that (ϕ(un)) is bounded and ϕ(un)(1+un)0 as n+ has a strong convergent subsequence in X.

    Remark 1. 1. It is clear from the above definition that if ϕ satisfies the (PS)-condition, then it satisfies the (C)-condition. However, there are functionals that satisfy the (C)-condition but do not satisfy the (PS)-condition (see [21]). Consequently, the (PS)-condition implies the (C)-condition.

    2. The (C)-condition and the (PS)-condition are equivalent if ϕ is bounded below (see [22]).

    Next, we present the following theorems, which will play a fundamental role in the proof of the main theorems.

    Theorem 3. (see [23]) Let (X,.) be a real Banach space; ϕC1(X,R) satisfies the (C)-condition; ϕ(0)=0, and the following conditions hold:

    1. There exist positive constants ρ and α such that ϕ(u)α for any uX with u=ρ.

    2. There exists a function eX such that e>ρ and ϕ(e)0.

    Then, the functional ϕ has a critical value cα, that is, there exists uX such that ϕ(u)=c and ϕ(u)=0 in X.

    Let X be a real, reflexive, and separable Banach space. Then there exist {ej}jNX and {ej}jNX such that

    X=¯span{ej:j=1,2,...},X=¯span{ej:j=1,2,...},

    and ei,ej=1 if i=j, ei,ej=0 if ij.

    We denote Xj=span{ej}, Yk=kj=1Xj, and Zk=¯+j=kXj.

    Theorem 4. (see [24]) Assume that X is a Banach space, and let ϕ:XR be an even functional of class C1(X,R) that satisfies the (C)condition. For every kN, there exists γk>ηk>0 such that

    (A1) bk:=inf{ϕ(u):uZk,u=ηk}+ as k+;

    (A2) ck:=max{ϕ(u):uYk,u=γk}0.

    Then, ϕ has a sequence of critical values tending to +.

    First of all, we are going to show that the functional ϕ fulfills the (C)-condition.

    Lemma 3.1. If assumptions (H0),(H1),(H2) and (H3) hold, then the functional ϕ satisfies the (C)-condition.

    Proof. Let (un)X0 be a Cerami sequence for ϕ, namely,

    (ϕ(un))  is bounded  and ϕ(un)X0(1+un)0, (3.1)

    which implies that

    sup|ϕ(un)|M  and ϕ(un),un=n(1), (3.2)

    where limn+n(1)=0 and M>0.

    We need to prove the boundedness of the sequence (un) in X0. To this end, assume to the contrary, that the sequence (un) is unbounded in X0. Without loss of generality, we can assume that un>1. By virtue of (H3), for n large enough, we have

    M+1ϕ(un)1q+ϕ(un),un=Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dx+Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dxΩF(x,un)dxΛ2G(x,un)dσ+1q+Ωf(x,un)undx+1q+Λ2g(x,un)undσ1q+Ω(|un|p(x)+b(x)|un|q(x))dx1q+Ω(|un|p(x)+b(x)|un|q(x))dx1q+Ω(|un|p(x)+b(x)|un|q(x))dx+1q+Ω(|un|p(x)+b(x)|un|q(x))dx1q+Ω(|un|p(x)+b(x)|un|q(x))dx1q+Ω(|un|p(x)+b(x)|un|q(x))dx+ΩF(x,un)dx+Λ2G(x,un)dσ,

    where F(x,un):=1q+f(x,un)unF(x,un)0 and G(x,un)=1q+g(x,un)unG(x,un)0.

    Then, we obtain

    M+1ΩF(x,un)dx+Λ2G(x,un)dσ,

    which implies

    M+1ΩF(x,un)dx, (3.3)

    and

    M+1Λ2G(x,un)dσ. (3.4)

    On the other hand, by Proposition 3, we have

    Mϕ(un)=Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dx+Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dxΩF(x,un)dxΛ2G(x,un)dσ1q+(Ω(|un|p(x)+b(x)|un|q(x))dx+Ω(|un|p(x)+b(x)|un|q(x))dx)ΩF(x,un)dxΛ2G(x,un)dσ1q+ˆρ(un)ΩF(x,un)dxΛ2G(x,un)dσ.

    Because un>1, we can obtain

    M1q+unpΩF(x,un)dxΛ2G(x,un)dσ. (3.5)

    Since un+ as n+, we deduce that

    ΩF(x,un)dx+Λ2G(x,un)dσ1q+unpM+    as n+. (3.6)

    Furthermore, using Proposition 3, we have

    ϕ(un)=Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dx+Ω(1p(x)|un|p(x)+b(x)q(x)|un|q(x))dxΩF(x,un)dxΛ2G(x,un)dσ1punq+ΩF(x,un)dxΛ2G(x,un)dσ.

    Then, we obtain

    ϕ(un)+ΩF(x,un)dx+Λ2G(x,un)dσ1punq+. (3.7)

    In view of condition (H2), there exist T1,T2>0 such that

    F(x,t)>|t|q+ for all  xΩ  and  |t|>T1,G(x,t)>|t|q+ for all  xΛ2  and  |t|>T2.

    Since F(x,.) and G(x,.) are continuous functions on [T1,T1] and [T2,T2], respectively, there exist C0,C0>0 such that

    |F(x,t)|C0 for all  (x,t)Ω×[T1,T1],|G(x,t)|C0for all  (x,t)Λ2×[T2,T2].

    Then, there exist two real numbers K and K, such that

    F(x,t)K  for all (x,t)Ω×R,G(x,t)K  for all (x,t)Λ2×R.

    Hence,

    F(x,un)K1punq+0,G(x,un)K1punq+0, (3.8)

    for all (x,n)¯Ω×N.

    Put βn=unun, so βn=1. Up to subsequences, for some βX0, we have

    βnβ  in X0,βnβ      in Ls(x)(Ω),βnβ      in Lr(x)(Λ2),βn(x)β(x) a.e., in Ω,βn(x)β(x) a.e., in Λ2, (3.9)

    for s(x)<p(x) and r(x)<p(x).

    Define the sets Ω0={xΩ:β(x)0} and Γ={xΛ2:β(x)0}.

    Obviously, since un+ as n+, we have

    |un(x)|=|βn(x)|un+,

    for any xΩ0Γ.

    Therefore, due to (H2), for all xΩ0Γ, we deduce

    F(x,un)1punq+=pF(x,un)|un(x)|q+|βn(x)|q++,G(x,un)1punq+=pG(x,un)|un(x)|q+|βn(x)|q++. (3.10)

    Thus, |Ω0|=0 and |Γ|=0. In fact, suppose by contradiction that |Ω0|0 or |Γ|0. Using (3.6), (3.7), (3.10), and Fatou's lemma, we get

    1=lim infn+ ΩF(x,un)dx+Λ2G(x,un)dσϕ(un)+ΩF(x,un)dx+Λ2G(x,un)dσlim infn+ΩF(x,un)dx+Λ2G(x,un)dσ1punq+lim infn+[Ω0F(x,un)1punq+dx+ΓG(x,un)1punq+dσ]lim supn+Ω0K1punq+dxlim supn+ΓK1punq+dσlim infn+Ω0F(x,un)1punq+dx+lim infn+ΓG(x,un)1punq+dσlim supn+Ω0K1punq+dxlim supn+ΓK1punq+dσ=lim infn+Ω0F(x,un)K1punq+dx+lim infn+Γg(x,un)K1punq+dσΩ0lim infn+F(x,un)K1punq+dx+Γlim infn+G(x,un)K1punq+dσΩ0lim infn+F(x,un)1punq+dx+Γlim infn+G(x,un)1punq+dσΩ0lim supn+K1punq+dxΓlim supn+K1punq+dx=+,

    which is a contradiction. Therefore, β(x)=0 for a.e. xΩ and for a.e. xΛ2.

    From (3.5) and (3.9), respectively, we can deduce that

    βn0    in Ls(x)(Ω),βn0    in Lr(x)(Λ2),βn(x)0    a.e. in Ω,βn(x)0    a.e. in Λ2, (3.11)

    for s(x)<p(x), r(x)<p(x), and

    0<1q+lim supn+[Ω|F(x,un)|unpdx+Λ2|G(x,un)|unpdσ]lim supn+Ω|F(x,un)|unpdx+lim supn+Λ2|G(x,un)|unpdσ. (3.12)

    Using (H0) and (H1), we obtain

    {0|un(x)|r1}|F(x,un)|unp dxC1{0|un(x)|r1}|un|+1s1(x)|un|s(x)unp dxC1|un|1unp+C1s1{0|un(x)|r1}|un|s1(x)p|βn|pdxC1C3ununp+Crsp1s|βn|ppCC3unp1+C1C4rsp1s1βnp0, as  n+, (3.13)

    where C3,C4>0, s is either s+1 or s1 and r1 comes from (H3).

    Put l1(x)=l(x)l(x)1. Since l1L(Ω) with l1(x)>Np, it follows that l1(x)p<p(x).

    On the other hand, by virtue of hypothesis (H3)(i), (3.3), and (3.11), we deduce

    {|un(x)|r1}|F(x,un)|unpdx2[{|un(x)|r1}(|F(x,un)||un|p)l1(x)dx]1l1(x)[{|un(x)|r1}|βn|l1(x)pdx]1l1(x)2c1l1(x)1[{|un(x)|r1}F(x,un)dx]1l1(x)[{|un(x)|r1}|βn|l1(x)pdx]1l1(x)2c1l1(x)1[ΩF(x,un)dx]1l1(x)[Ω|βn|l1(x)pdx]1l1(x)2c1l1(x)1(M+1)1l(x)[Ω|βn|l(x)p dx]1l(x)0 as  n+.

    Combining this with (3.13), we obtain

    Ω|F(x,un)|unpdx={0|un(x)|r1}|F(x,un)|unpdx+{|un(x)|r1}|F(x,un)|unpdx0, as  n+. (3.14)

    Similarly, let l2(x)=l2(x)l2(x)1. Since l2L(Ω) with l2(x)>N1p1, it follows that l2(x)p<p(x). Then, by (H3)(ii), (3.4), and (3.11), we can prove in a similar way that

    Λ2|G(x,un)|unpdσ0,    as  n+. (3.15)

    Consequently, combining (3.14) with (3.15), we obtain

    Ω|F(x,un)|unpdx+Λ2|G(x,un)|unpdσ0,    as  n+,

    which is a contradiction to (3.12). Thus, (un) is bounded in X0.

    Finally, we need to prove that any (C)-sequence has a convergent subsequence. Let (un)X0 be a (C)-sequence. Then, (un) is bounded in X0. Passing to the limit, if necessary, to a subsequence, from Proposition 4, we have

    unu in X0, unu  in Ls1(x)(Ω), unu  in Ls2(x)(Λ2),un(x)u(x)  a.e. xΩ, un(x)u(x)  a.e. xΛ2, (3.16)

    for 1s1(x)<p(x) and 1s2(x)<p(x). It is easy to check from (H0), (3.16) and Hölder's inequality that

    |Ωf(x,un)(unu)dx|C1|1+|un|s1(x)1|s1(x)|unu|s1(x)0 as  n+, (3.17)

    and

    |Λ2g(x,un)(unu)dσ|C2|1+|un|s2(x)1|s2(x)|unu|s2(x)0 as  n+, (3.18)

    where 1s1(x)+1s1(x)=1 and 1s2(x)+1s2(x)=1.

    Next, since unu, from (3.1), we have

    ϕ(un),unu0, as  n+. (3.19)

    Then

    ϕ(un),unu=A(un),unuΩf(x,un)(unu)dxΛ2g(x,un)(unu)dσ0 as  n+,

    where A is given in Proposition 5.

    Finally, the combination of (3.17), (3.18), and (3.19) implies

    A(un),unu0 as  n+.

    Since the operator A satisfies the (S+) property in view of Proposition 5, we can obtain that unu in X0. The proof is complete.

    Proof of Theorem 1

    Let us check that the functional ϕ satisfies the geometric conditions of the mountain pass in Theorem 3. By Lemma 3.1, ϕ satisfies the (C)condition. According to the definition of ϕ, we have ϕ(0)=0. Then, to apply Theorem 3, it remains to prove that

    (ⅰ) There exist positive constants ρ and α such that ϕ(u)α for any uX0 with u=ρ.

    (ⅱ) There exists a function eX0 such that e>ρ and ϕ(e)0.

    For (i), let u<1. Then, by Proposition 3, we have

    ϕ(u)1p+ˆρ(u)ΩF(x,u)dxΛ2G(x,u)dσ1p+uq+ΩF(x,u)dxΛ2G(x,u)dσ. (3.20)

    Using (H0) and (H4), for ε>0 be small enough, there exist C1(ε),C2(ε)>0 such that

    F(x,t)ε|t|p++C1(ε)|t|s1(x),(x,t)Ω×R,G(x,t)ε|t|p++C2(ε)|t|s2(x),(x,t)Λ2×R. (3.21)

    Since p+<s1(x)<p(x) and p+<s2(x)<p(x) for all x¯Ω and for all x¯Λ2 in view of condition (H1), we have from Proposition 4 that

    X0Lp+(Ω),X0Ls1(x)(Ω),X0Lp+(Λ2),X0Ls2(x)(Λ2).

    So, there exist ci>0(i=3,...6) such that

    |u|p+c3u,|u|s1(x)c4u,  uX0|u|p+,Λ2c6u,|u|s(x)c7u,  uX0.

    Therefore, by (3.20) and (3.21), for u<1 sufficiently small, we obtain

    ϕ(u)1q+uq+εΩ|u|p+dxC1(ε)Ω|u|s1(x)dxεΛ2|u|p+dσC2(ε)Λ2|u|s2(x)dσ1q+uq+εcp+3up+C1(ε)cs14us1εcp+5up+C2(ε)cs26us2.

    Since s1>p+ in view of condition (H1) and u<1, then us1<up+. Thus, we obtain

    ϕ(u)1q+uq+εcp+3up+C1(ε)cs14up+εcp+5up+C2(ε)cs26us21q+uq+(εcp+3+C1(ε)cs14+εcp+5)up+C2(ε)cs26us2.

    Since s2>q+p+, then by the standard argument, there exist positive constants ρ and α such that ϕ(u)α for any uX0 with u=ρ.

    Next, we affirm that there exists eX0 with u>ρ such that

    ϕ(e)<0. (3.22)

    In fact, from (H2), it follows that for every k>0, there exist constants Tk and Tk such that

    F(x,t)>k|t|q+ for all  xΩ  and  |t|>Tk,G(x,t)>k|t|q+ for all  xΛ2  and  |t|>Tk.

    Since F(x,.) and G(x,.) are continuous functions on [Tk,Tk] and [Tk,Tk], respectively, there exist constants C0,C0>0 such that

    |F(x,t)|C0 for all  (x,t)Ω×[Tk,Tk],|G(x,t)|C0 for all  (x,t)Λ2×[Tk,Tk].

    Thus,

    F(x,t)k|t|q+C0, for all  (x,t)Ω×R,G(x,t)k|t|q+C0, for all  (x,t)Λ2×R. (3.23)

    Let wX0{0} such that w=1 and l>1 be large enough. Using the above inequality, we obtain

    ϕ(lw)=Ω(1p(x)|lw|p(x)+b(x)q(x)|lw|q(x))dx+Ω(1p(x)|lw|p(x)+b(x)q(x)|lw|q(x))dxΩF(x,lw)dxΛ2G(x,lw)dσlq+pˆρ(w)klq+Ω|w|q+dxklq+Λ2|w|q+dσ+C0|Ω|+C0|Λ2|lq+pklq+Ω|w|q+dxklq+Λ2|w|q+dσ+C0|Ω|+C0|Λ2|=lq+(1pkΩ|w|q+dxkΛ2|w|q+dσ)+C0|Ω|+C0|Λ2|.

    As

    1pkΩ|w|q+dxkΛ2|w|q+dσ<0,

    for k large enough, we deduce

    ϕ(lw), as  l+.

    Thus, there exist t0>1 and e=t0wX0¯Bρ(0) such that ϕ(e)<0.

    Proof of Theorem 2

    To prove Theorem 2, we need the following auxiliary lemmas:

    Lemma 3.2. (see [25,26]) For sC+(¯Ω) and rC+(¯Λ2) such that s(x)<p(x) for all x¯Ω and r(x)<p(x) for all x¯Λ2. Let

    δk=sup{|u|s(x):u=1,uZk},δk=sup{|u|s(x),Λ2:u=1,uZk}.

    Then, limk+δk=limk+δk=0.

    Lemma 3.3. (see [27]) For all sC+(¯Ω) (rC+(¯Λ2)) and uLs(x)(Ω) (vLr(x)(Λ2)), there exists yΩ (zΛ2) such that

    Ω|u|s(x)dx=|u|s(y)s(x),Λ2|u|r(x)dσ=|u|r(z)r(x),Λ2. (3.24)

    Now, we return to the proof of Theorem 2. To this end, based on Fountain Theorem 4, we will show that the problem (1.1) possesses infinitely many weak solutions with unbounded energy. Evidently, according to (H5), ϕ is an even functional. By Lemma 3.1, we know that ϕ satisfies the (C)-condition. Then, to prove Theorem 2, it only remains to verify the following assertions:

    (A1) bk:=inf{ϕ(u):uZk,u=ηk}+ as k+,

    (A2) ck:=max{ϕ(u):uYk,u=γk}0.

    (A1) For any uZk such that u=ηk>1. It follows from (H0), Proposition 3, and Lemma 3.3 that

    ϕ(u)=Ω(1p(x)|u|p(x)+b(x)q(x)|u|q(x))dx+Ω(1p(x)|u|p(x)+b(x)q(x)|u|q(x))dxΩF(x,u)dxΛ2G(x,u)dσ1q+upC1Ω|u|dxC1Ω|u|s1(x)s1(x)dxC2Λ2|u|dσC2Λ2|u|s2(x)s2(x)dσ1q+upC1c1uC1s1|u|s1(y)s1(x)C2c2uC2s2|u|s2(z)s2(x),λ21q+upc3uC1s1|u|s1(y)s1(x)C2s2|u|s2(z)s2(x),λ2,

    where c3=max{C1c1,C2c2}.

    Then, it follows that

    ϕ(u){1q+upc3uC1s1C2s2   if |u|s1(x)1,|u|s2(x)11q+upc3uC1s1(δku)s+1C2s2(δku)s+2  if |u|s1(x)>1,|u|s2(x)>11q+upc3uC1s1C2s2(δku)s+2   if |u|s1(x)1,|u|s2(x)>11q+upc3uC1s1(δku)s+1C2s2   if |u|s1(x)>1,|u|s2(x)11q+upc3u2C1s1(δku)s+12C2s2(δku)s+22C1s12C2s212q+up+(14q+up2C1s1δs+1kus+1)+(14q+up2C1s2(δk)s+2us+2)c3u2C1s12C2s2. (3.25)

    Let us consider the following equations:

    14q+tp2C1s1δs+1kts+1=0, (3.26)

    and

    14q+tp2C1s2(δk)s+2ts+2=0. (3.27)

    Let ak and dk be the two non-zero solutions of (3.26) and (3.27), respectively. Then, we obtain

    ak=(4q+2C1s1δs+1k)1ps+1+   and   dk=(4q+2C2s2(δk)s+2)1ps+2+   as k+.

    We fix ηk as follows

    ηk=min{ak,dk}.

    Then, by Lemma 3.2, (3.25) and s+1,s+2>q+>p, we obtain

    ϕ(u)12q+ηpkc3ηkC6+   as  k+,

    where C6>0. Hence, (A1) holds.

    (A2) In view of Proposition 3 and (3.23), for uYk with u|>1, we have

    ϕ(u)1puq+ΩF(x,u)dxΛ2G(x,u)dσ1puq+kΩ|u|q+dxkΛ2|u|q+dσ+C0|Ω|+C0|Λ2|1puq+k(|u|q+q++|u|q+q+,Λ2)+C0|Ω|+C0|Λ2|.

    Since dimYk<, then all norms are equivalent in Yk. Therefore, as 1p<1, for k large enough, we obtain

    ϕ(u) as  u+.

    Finally, the assertion (A2) is also valid.

    This completes the proof.

    Mahmoud El Ahmadi: Writing-original draft, Writing-review & editing; Mohammed Barghouthe: Formal Analysis, Methodology; Anass Lamaizi: Formal Analysis; Mohammed Berrajaa: Supervision, Validation.

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

    The authors would like to thank the referees for their valuable comments and suggestions, which have improved the quality of this paper.

    The authors declare there is no conflict of interest.



    [1] D. D. Chaplin, Overview of the immune response, J. Allergy. Clin. Immun., 125 (2010), S3–S23. https://doi.org/10.1016/j.jaci.2009.12.980 doi: 10.1016/j.jaci.2009.12.980
    [2] C. A. Jr. Janeway, P. Travers, M. Walport, et al., Immunobiology: The Immune System in Health and Disease. 5th edition. New York: Garland Science, 2001.
    [3] C. R. Maldini, G. I. Ellis, J. L. Riley, CAR T cells for infection, autoimmunity and allotransplantation, Nat. Rev. Immunol., 18 (2018), 605–616. https://doi.org/10.1038/s41577-018-0042-2 doi: 10.1038/s41577-018-0042-2
    [4] L. B. Nicholson, The immune system, Essays Biochem., 60 (2016), 275–301. https://doi.org/10.1042/EBC20160017 doi: 10.1042/EBC20160017
    [5] J. M. Carton, W. R. Strohl, Protein therapeutics (introduction to biopharmaceuticals): Introduction to Biological and Small Molecule Drug Research and Development, Elsevier, 2013,127–159, ISBN 9780123971760, https://doi.org/10.1016/B978-0-12-397176-0.00004-2.10.1016/B978-0-12-397176-0.00004-2
    [6] M. Lever, T. D. C. Hirata, P. Russo, H. I. Nakaya, Systems immunology, Theoretical and Applied Aspects of Systems Biology., Springer International Publishing (2018), 159–173. https://doi.org/10.1007/978-3-319-74974-7.
    [7] N. Chirmule, V. Jawa, B. Mibohm, Immunogenicity to therapeutic proteins: impact on PK/PD and efficacy, AAPS J., 14 (2012), 296–302. https://doi.org/10.1208/s12248-012-9340-y doi: 10.1208/s12248-012-9340-y
    [8] K. Bloem, B. Hernández-Breijo, A. Martínez-Feito, T. Rispens, Immunogenicity of therapeutic antibodies: Monitoring antidrug antibodies in a clinical context, Ther Drug Monit., 4 (2017), 327–332. https://doi.org/10.1097/FTD.0000000000000404
    [9] J. J. P. Ruixo, P. Ma, A. T. Chow, The utility of modeling and simulation approaches to evaluate immunogenicity effect on the therapeutic protein pharmacokinetics, AAPS J., 15 (2013), 172–182. https://doi.org/10.1208/s12248-012-9424-8 doi: 10.1208/s12248-012-9424-8
    [10] R. Eftimie, J. J. Gillard, D. A. Cantrell, Mathematical models for immunology: Current state of the art and future research directions, B. Math. Biol., 78 (2016), 2091–2134. https://doi.org/10.1007/s11538-016-0214-9 doi: 10.1007/s11538-016-0214-9
    [11] S. Banerjee, S. Khajanchi, S. Chaudhuri, A mathematical model to elucidate brain tumor abrogation by immunotherapy with T11 target structure, PLoS ONE, 10 (2015), e0123611, https://doi.org/10.1371/journal.pone.0123611. doi: 10.1371/journal.pone.0123611
    [12] J. Reyes-Silveyra, A. R. Mikler, Modeling immune response and its effect on infectious disease outbreak dynamic, Theor. Biol. Med. Model., 13 (2016). https://doi.org/10.1186/s12976-016-0033-6 doi: 10.1186/s12976-016-0033-6
    [13] S. Kathman, T. M. Thway, L. Zhou, S. Lee, Utility of a bayesian mathematical model to predict the impact of immunogenicity on pharmacokinetics of therapeutic proteins, AAPS J., 18 (2016), 424–431. https://doi.org/10.1208/s12248-015-9853-2 doi: 10.1208/s12248-015-9853-2
    [14] W. L. Duan, H. Fang, C. Zeng, The stability analysis of tumor-immune responses to chemotherapy system with gaussian white noises. Chaos, Soliton. Fract., 127 (2019), 96–102. https://doi.org/10.1016/j.chaos.2019.06.030. doi: 10.1016/j.chaos.2019.06.030
    [15] G. A. Bocharov, V. Volpert, B. Ludewig, A. Meyerhans, Mathematical modeling of the immune system in homeostasis, infection and disease, Front. Immunol., 10 (2020), 2944. https://doi.org/10.3389/fimmu.2019.02944 doi: 10.3389/fimmu.2019.02944
    [16] G. A. Bocharov, D. S. Grebennikov, R. S. Savinkov, Mathematical immunology: From phenomenological to multiphysics modelling, Russ. J. Numer. Anal. Mathematical M., 35 (2020), 203–213. https://doi.org/10.1515/rnam-2020-0017 doi: 10.1515/rnam-2020-0017
    [17] S. A. Alharbi, A. Z. Rambely, Dynamic behaviour and stabilisation to boost the immune system by complex interaction between tumour cells and vitamins intervention, Adv. Differ. Equ-Ny., 1 (2020), 412. https://doi.org/10.1186/s13662-020-02869-6 doi: 10.1186/s13662-020-02869-6
    [18] A. Fenton, J. Lello, M. B. Bonsall, Pathogen responses to host immunity: The impact of time delays and memory on the evolution of virulence, Proc. R. Soc. B., 273 (2006). https://doi.org/10.1098/rspb.2006.3552 doi: 10.1098/rspb.2006.3552
    [19] F. A. Rihan, D. H. A. Rahman, Delay differential model for tumour-immune dynamics with HIV infection of CD4+ T-cells, Int. J. Comput. Math., 90 (2013), 594–614, http://dx.doi.org/10.1080/00207160.2012.726354. doi: 10.1080/00207160.2012.726354
    [20] F. A. Rihan, D. H. Abdelrahman, F. Al-Maskari, F. Ibrahim, M. A. Abdeen, Delay differential model for tumour-immune response with chemoimmunotherapy and optimal control, Comput. Math. Method. M., 2014 (2014), Article ID 982978. http://dx.doi.org/10.1155/2014/982978.10.1155/2014/982978.
    [21] S. Khajanchi, S. Banerjee, Stability and bifurcation analysis of delay induced tumor immune interaction model, Appl. Math. Comput., 248 (2014), 652–671. https://doi.org/10.1016/j.amc.2014.10.009. doi: 10.1016/j.amc.2014.10.009
    [22] S. Kayan, H. Merdan, R. Yafia, S. Goktepe, Bifurcation analysis of a modified tumor-immune system interaction model involving time delay, Math. Model. Nat. Pheno., 12 (2017), 120–145. https://doi.org/10.1051/mmnp/201712508 doi: 10.1051/mmnp/201712508
    [23] F. A. Rihan, S. Lakshmanan, H. Maurer, Optimal control of tumour-immune model with time-delay and immuno-chemotherapy. Appl. Math. Comput., 353 (2019), 147–165. https://doi.org/10.1016/j.amc.2019.02.002. doi: 10.1016/j.amc.2019.02.002
    [24] J. P. Mendonça, I. Gleria, M. L. Lyra, Delay-induced bifurcations and chaos in a two-dimensional model for the immune response, Physica A, 517 (2019), 484–490. https://doi.org/10.1016/j.physa.2018.11.039 doi: 10.1016/j.physa.2018.11.039
    [25] F. Fatehi, Y. N. Kyrychko, K. B. Blyuss, Time-delayed model of autoimmune dynamics, Math. Biosci. Eng., 16 (2019), 5613–5639. https://doi.org/10.3934/mbe.2019279 doi: 10.3934/mbe.2019279
    [26] Kasbawati, Mariani, N. Erawaty, N. Aris, A mathematical study of effects of delays arising from the interaction of anti-drug antibody and therapeutic protein in the immune response system, AIMS Math., 5 (2020), 7191–7213. https://doi.org/10.3934/math.2020460 doi: 10.3934/math.2020460
    [27] P. Das, P. Das, S. Das, Effects of delayed immune-activation in the dynamics of tumor-immune interactions, Math. Model. Nat. Pheno., 15 (2020), 45. https://doi.org/10.1051/mmnp/2020001 doi: 10.1051/mmnp/2020001
    [28] Q. Tang, G. Zhang, Stability and Hopf bifurcations in a competitive tumour-immune system with intrinsic recruitment delay and chemotherapy, Math. Biosci. Eng., 18 (2021), 1941–1965. https://doi.org/10.3934/mbe.2021101 doi: 10.3934/mbe.2021101
    [29] W. L. Duan, L. Lin, Noise and delay enhanced stability in tumor-immune responses to chemotherapy system, Chaos, Soliton. Fract., 148 (2021), 111019. https://doi.org/10.1016/j.chaos.2021.111019. doi: 10.1016/j.chaos.2021.111019
  • Reader Comments
  • © 2022 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(1903) PDF downloads(87) Cited by(2)

Figures and Tables

Figures(3)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog