Processing math: 100%
Research article Special Issues

Global existence and boundedness of chemotaxis-fluid equations to the coupled Solow-Swan model

  • In this paper, we consider the following Keller-Segel-(Navier)-Stokes system to the coupled Solow-Swan model

    {nt+un=Δnχ(nc)+μ1nμ2nk,xΩ,t>0,ct+uc=Δcc+μ3cαw1α,xΩ,t>0,wt+uw=Δww+n,xΩ,t>0,ut+κ(uu)=ΔuP+nΦ,u=0,xΩ,t>0,

    in a smooth bounded domain ΩRN(N=2,3) with no-flux boundary for n,c,w and no-slip boundary for u, where the parameters χ>0,α(0,1),μ1R,μ20,μ3>0 and κ{0,1},kN. Due to the interference of the fractional nonlinear term of the Solow-Swan model, we use the Moser-Trudinger inequality to obtain the global existence of the solution for two-dimensional case without logistic source. For three-dimensional case, we control the required estimation with the help of the negative term of logistic source to obtain the boundedness and asymptotic behavior. In the process of estimating the corresponding term, we find the order of the negative term of the logistic source is related to the spatial dimension, and we give the decay estimate of the corresponding solutions when μ1<0 or μ1=0,μ2>0.

    Citation: Jie Wu, Zheng Yang. Global existence and boundedness of chemotaxis-fluid equations to the coupled Solow-Swan model[J]. AIMS Mathematics, 2023, 8(8): 17914-17942. doi: 10.3934/math.2023912

    Related Papers:

    [1] Yanlin Li, Nasser Bin Turki, Sharief Deshmukh, Olga Belova . Euclidean hypersurfaces isometric to spheres. AIMS Mathematics, 2024, 9(10): 28306-28319. doi: 10.3934/math.20241373
    [2] Nasser Bin Turki, Sharief Deshmukh, Olga Belova . A note on closed vector fields. AIMS Mathematics, 2024, 9(1): 1509-1522. doi: 10.3934/math.2024074
    [3] Meraj Ali Khan, Ali H. Alkhaldi, Mohd. Aquib . Estimation of eigenvalues for the α-Laplace operator on pseudo-slant submanifolds of generalized Sasakian space forms. AIMS Mathematics, 2022, 7(9): 16054-16066. doi: 10.3934/math.2022879
    [4] Sharief Deshmukh, Mohammed Guediri . Some new characterizations of spheres and Euclidean spaces using conformal vector fields. AIMS Mathematics, 2024, 9(10): 28765-28777. doi: 10.3934/math.20241395
    [5] Songting Yin . Some rigidity theorems on Finsler manifolds. AIMS Mathematics, 2021, 6(3): 3025-3036. doi: 10.3934/math.2021184
    [6] Ibrahim Aldayel . Value of first eigenvalue of some minimal hypersurfaces embedded in the unit sphere. AIMS Mathematics, 2023, 8(11): 26532-26542. doi: 10.3934/math.20231355
    [7] Amira Ishan . On concurrent vector fields on Riemannian manifolds. AIMS Mathematics, 2023, 8(10): 25097-25103. doi: 10.3934/math.20231281
    [8] Ibrahim Al-Dayel, Meraj Ali Khan . Ricci curvature of contact CR-warped product submanifolds in generalized Sasakian space forms admitting nearly Sasakian structure. AIMS Mathematics, 2021, 6(3): 2132-2151. doi: 10.3934/math.2021130
    [9] Hanan Alohali, Sharief Deshmukh . Some generic hypersurfaces in a Euclidean space. AIMS Mathematics, 2024, 9(6): 15008-15023. doi: 10.3934/math.2024727
    [10] Ibrahim Al-Dayel, Sharief Deshmukh, Olga Belova . Characterizing non-totally geodesic spheres in a unit sphere. AIMS Mathematics, 2023, 8(9): 21359-21370. doi: 10.3934/math.20231088
  • In this paper, we consider the following Keller-Segel-(Navier)-Stokes system to the coupled Solow-Swan model

    {nt+un=Δnχ(nc)+μ1nμ2nk,xΩ,t>0,ct+uc=Δcc+μ3cαw1α,xΩ,t>0,wt+uw=Δww+n,xΩ,t>0,ut+κ(uu)=ΔuP+nΦ,u=0,xΩ,t>0,

    in a smooth bounded domain ΩRN(N=2,3) with no-flux boundary for n,c,w and no-slip boundary for u, where the parameters χ>0,α(0,1),μ1R,μ20,μ3>0 and κ{0,1},kN. Due to the interference of the fractional nonlinear term of the Solow-Swan model, we use the Moser-Trudinger inequality to obtain the global existence of the solution for two-dimensional case without logistic source. For three-dimensional case, we control the required estimation with the help of the negative term of logistic source to obtain the boundedness and asymptotic behavior. In the process of estimating the corresponding term, we find the order of the negative term of the logistic source is related to the spatial dimension, and we give the decay estimate of the corresponding solutions when μ1<0 or μ1=0,μ2>0.



    The Laplace operator Δ acting on smooth functions of a Riemannian manifold (Mn,g), is defined by Δf=div(f), where f is the gradient of f. It is known that Δ is a self adjoint elliptic operator with respect to the inner product (,) defined on the algebra of smooth functions C(Mn) with compact support by

    (f,h)=MnfhdVg,

    where dVg is the volume form on Mn with respect to the metric g. If Δf=λf, for a constant λ, then λ is said to be an eigenvalue of the Laplace operator Δ, the negative sign in the definition chosen so that for non-zero eigenvalue λ, λ>0. The set of eigenvalues λi of the Laplace operator Δ on a Riemannian manifold (Mn,g) is called the spectrum of (Mn,g). Spectra of known Riemannian manifolds such as the sphere Sn(c) and the real projective space RPn are known, and a nice description could be found in [2] (Chapter 2, Section 5). Two Riemannian manifolds of same dimension (Mn,g) and (¯Mn,¯g) having the same spectra are said to be isospectral manifolds, whereas they are said to be isometric if there exists a diffeomorphism ϕ:Mn¯Mn that preserves the metric ϕ(¯g)=gϕ, that is, ϕ is an isometry. In the mid Nineteenth century, it was an open question whether isospectral Riemannian manifolds are isometric. We see that physicists followed this question under the topic "On hearing the shape of a drum" or "Can one hear the shape of a drum?" (cf. [5,11]), while Milnor constructed two flat tori in dimension 16, which are isospectral but not isometric (cf. [12]). In [12], the author constructed two 2-dimensional compact manifolds of constant negative curvature which are isospectral and not isometric. Further, Ejiri constructed two non-flat compact Riemannian manifolds which are isospectral but not isometric (cf. [6]). This initiated an interest in comparing spectra of two Riemannian manifolds, for further results in this direction, we refer to (cf. [1,7,8,9,10,15]). One of the natural questions is under what conditions two isospectral Riemannian manifolds are isometric? We know that the n-sphere Sn(c) has first non-zero eigenvalue λ1=nc, and, supposing an n-dimensional Riemannian manifold (Mn,g) also has non-zero eigenvalue λ1=nc, we reduce the above general question to the following specific question: Under what condition is an n-dimensional Riemannian manifold (Mn,g) that has non-zero eigenvalue λ1=nc for a positive constant c isometric to Sn(c)? In order to answer this question, naturally, we need a smooth function ρ on an n-dimensional Riemannian manifold (Mn,g) that satisfies Δρ=ncρ. In order to address this issue of finding a smooth function, we consider two situations:

    (ⅰ) Considering (Mn,g) as an isometrically immersed hypersurface in the Euclidean space (Rn+1,,) with isometric immersion f:(Mn,g)(Rn+1,,), though there are several known functions on the hypersurface (Mn,g), namely, the mean curvature function, the scalar curvature function, as well as the support function, but they cannot be eigenfunctions of the Laplace operator corresponding to a non-zero eigenvalue, while (Mn,g) is isometric to the sphere Sn(c). Then, we go for the alternative, namely, for a constant unit vector field a on the Euclidean space Rn+1, we get a smooth function ρ=f,a on the hypersurface (Mn,g), which we require to satisfy Δρ=ncρ, c>0.

    (ⅱ) We seek the Riemannian manifold (Mn,g), and there exists a non-trivial concircular vector field ξ with potential function ρ and we require that the potential function ρ satisfies Δρ=ncρ, c>0.

    In this article, we explore the above two situations that an n-dimensional Riemannian manifold (Mn,g) has an eigenvalue of the Laplace operator same as the first non-zero eigenvalue of the sphere Sn(c), and find an additional condition so that these two are isometric (See Theorems 1, 2, 3). In the language of physics, "yes, we could hear a node through an apparatus to predict the shape of a drum".

    On an n-dimensional Riemannian manifold (Mn,g), denote by the Riemannian connection. The curvature tensor R of (Mn,g) is given by

    R(X,Y)Z=[X,Y]Z[X,Y]ZX,Y,ZΓ(TMn), (2.1)

    where Γ(TMn) is the space of smooth sections of the tangent bundle TMn. The Ricci tensor Ric of (Mn,g) is a symmetric tensor given by

    Ric(X,Y)=nl=1g(R(el,X)Y,el)X,YΓ(TMn), (2.2)

    where {e1,..,en} is a local orthonormal frame (or lof) on (Mn,g). The Ricci operator Q of (Mn,g) is related to the Ricci tensor by

    Ric(X,Y)=g(QX,Y)X,YΓ(TMn)

    and therefore Q is a symmetric (1,1) tensor field. The scalar curvature τ of the Riemannian manifold (Mn,g) is given by

    τ=nl=1Ric(el,el).

    The following formula is well known (cf. [2]),

    12τ=nl=1(elQ)(el),

    where τ is the gradient of τ, and the covariant derivative is given by

    (XQ)(Y)=XQYQ(XY).

    Given a smooth function f:MnR on a Riemannian manifold (Mn,g), the Laplace operator Δ acts on f, given by

    Δf=div(f),

    where f is the gradient of f and

    divX=nl=1g(elX,el).

    If (Mn,g) is closed, then Stokes's theorem implies

    Mn(divX)dVg=0,

    where dVg is the volume element of (Mn,g).

    Given a symmetric (1,1) tensor field T on an n-dimensional Riemannian manifold (Mn,g) with trace t, that is,

    t=nl=1g(Tel,el),

    then the Cauchy–Schwartz inequality is

    T21nt2, (2.3)

    where

    T2=nl=1g(Tel,Tel).

    Moreover, the equality in (2.3) holds if and only if

    T=tnI,

    where I is the identity (1,1) tensor.

    Suppose on an n-dimensional closed Riemannian manifold (Mn,g) there is a smooth function f that satisfies Δf=ncf for a constant c. Then we have fΔf=ncf2. Then, integrating by parts, the last equation leads to

    Mnf2dVg=ncMnf2dVg. (2.4)

    In this section, we consider an n-dimensional Riemannian manifold (Mn,g) that admits an isometric immersion f:(Mn,g)(Rn+1,,) into the Euclidean space (Rn+1,,), where , is the Euclidean inner product. We denote by ς the unit normal to the hypersurface (Mn,g), and by B the shape operator with respect to the isometric immersion f. Then, denoting the Euclidean connection on (Rn+1,,) by ¯ and the Riemannian connection on the hypersurface (Mn,g) by , we have the following fundamental equations for the hypersurface (cf. [2])

    ¯XY=XY+gBX,YςX,YΓ(TMn) (3.1)

    and

    ¯Xς=BXXΓ(TMn). (3.2)

    The mean curvature β of the hypersurface (Mn,g) is given by

    β=1nnl=1g(Bel,el), (3.3)

    where {e1,..,en} is a local orthonormal frame on (Mn,g).

    The curvature tensor, the Ricci tensor, and the scalar curvature of the hypersurface (Mn,g) are given by (cf. [2])

    R(X,Y)Z=g(BY,Z)BXg(BX,Z)BXX,Y,ZΓ(TMn),
    Ric(X,Y)=nβg(BX,Y)g(BX.BY)X,YΓ(TMn)

    and

    τ=n2β2B2.

    Also, the shape operator B has the following wonderful property (Codazzi equation of the hypersurfaces in flat spaces)

    (XB)(Y)=(YB)(X)X,YΓ(TMn). (3.4)

    Using Eqs (3.3) and (3.4) as well as the symmetry of the shape operator, we have for any XΓ(TMn)

    nX(β)=nl=1g(XBel,el)+nl=1g(Bel,Xel)=nl=1g((XB)(el)+B(Xel),el)+nl=1g(Bel,Xel)=nl=1g((elB)(el),X)+2nl=1g(Bel,Xel).

    Note that Bel=jg(Bel,ej)ej and Xel=kωkl(X)ek, where g(Bel,ej) is symmetric while the connection forms ωkl are skew symmetric. Therefore,

    nl=1g(Bel,Xel)=jklg(Bel,ej)ωkl(X)g(ej,ek)=nl=1g(Bel,ej)ωjl(X)=0.

    Thus, above equation becomes

    nX(β)=nl=1g((elB)(el),X),

    that is, we have

    nl=1(elB)(el)=nβ. (3.5)

    Treating the isometric immersion f:(Mn,g)(Rn+1,,) as a position vector of points of Mn in Rn+1, and defining σ=f,ς, called the support function of the hypersurface (Mn,g), we express f as

    f=ξ+σς,

    where ξΓ(TMn) is tangential to (Mn,g). Differentiating equation (2.4), while using Eqs (3.1) and (3.2), we have upon equating the tangential and normal parts

    Xξ=X+σBX (3.6)

    and

    σ=Bξ.

    Taking a constant unit vector field a on the Euclidean space Rn+1 (for instance a coordinate vector field), we define a smooth function h on the hypersurface (Mn,g), by h=a,ς. Denoting the tangential part of a to the hypersurface (Mn,g) by ζ, we have

    a=ζ+hς. (3.7)

    Differentiating the above equation with respect to XΓ(TMn), while using Eqs (3.1) and (3.2), we have upon equating the tangential and normal parts

    Xζ=hBXXΓ(TMn) (3.8)

    and

    h=Bζ. (3.9)

    Now, we prove the main result of this section.

    Theorem 1. An n-dimensional compact and connected isometrically immersed hypersurface f:(Mn,g)(Rn+1,,) in the Euclidean space (Rn+1,,) with mean curvature β and a constant unit vector a=ζ+hς on Rn+1, where the function ρ=f,a satisfies Δρ=ncρ for a positive constant c, is isometric to the sphere Sn(c) if and only if the Ricci curvature Ric(ζ,ζ) satisfies

    MnRic(ζ,ζ)dVgn(n1)Mnh2β2dVg.

    Proof. Consider an n-dimensional compact and connected Riemannian manifold (Mn,g) that admits an isometric immersion f:(Mn,g)(Rn+1,,) in the Euclidean space (Rn+1,,) with shape operator B, mean curvature β, and a constant unit vector a=ζ+hς on Rn+1, where the function ρ=f,a satisfies

    Δρ=ncρ (3.10)

    for a positive constant c. Also, the Ricci curvature Ric(ζ,ζ) satisfies

    MnRic(ζ,ζ)dVgn(n1)Mnh2β2dVg. (3.11)

    Now, differentiating ρ=f,ain the direction of XΓ(TMn), we get X(ρ)=X,a=f,ζ. This gives us the gradient of ρ as

    ρ=ζ. (3.12)

    The Hessian operator Hρ of the function ρ is given by HρX=Xρ, XΓ(TMn), and using Eqs (3.8) and (3.12), we arrive at

    HρX=hBXXΓ(TMn). (3.13)

    Taking the trace in the above equation and taking account of Eq (3.10), we get

    cρ=hβ (3.14)

    and therefore, through Eq (3.13), we conclude

    HρX+cρX=hBXhβXXΓ(TMn).

    From the above equation, we reach

    Hρ+cρI2=h2BβI2. (3.15)

    Next, using Eq (3.8), we have

    h(BXβX)=XζhβX,

    which yields

    h2BβI2=ζ2+nh2β22hβdivζ.

    Inserting divζ=nhβ (an outcome of Eq (3.8)), in the above equation, we arrive at

    h2BβI2=ζ2nh2β2. (3.16)

    Recalling the following well known integral formula (cf. [16])

    Mn(Ric(ζ,ζ)+12|£ζg|2ζ2(divζ)2)dVg=0,

    and integrating Eq (3.16) while using the above integral formula, we conclude

    Mnh2BβI2dVg=Mn(Ric(ζ,ζ)+12|£ζg|2(divζ)2nh2β2)dVg. (3.17)

    Using Eq (3.8), we compute

    (£ζg)(X,Y)=2hg(BX,Y)X,YΓ(TMn)

    and consequently, we have

    12|£ζg|2=2h2B2.

    Thus, inserting above equation and divζ=nhβ in Eq (3.17) confirms

    Mnh2BβI2dVg=Mn(Ric(ζ,ζ)+2h2B2n2h2β2nh2β2)dVg,

    that is,

    Mnh2BβI2dVg=Mn(Ric(ζ,ζ)+2h2(B2nβ2)n(n1)h2β2)dVg. (3.18)

    For a local orthonormal frame {e1,..,en}, we have

    BβI2=kg(Bekβek,Bekβek)=B2+nβ22βkg(Bek,ek)=B2nβ2.

    Utilizing the above equation in (3.18), we arrive at

    Mnh2BβI2dVg=Mn(n(n1)h2β2Ric(ζ,ζ))dVg.

    Inserting from Eq (3.15) in the above equation, we have

    MnHρ+cρI2dVg=Mn(n(n1)h2β2Ric(ζ,ζ))dVg

    and treating it with inequality (3.11) allows us to reach the conclusion

    Hρ=cρI.

    Note that ρ satisfies Eq (3.10), that is, Δρ=ncρ for a non-zero constant c. We claim that ρ can not be a constant, for if it were, Equation (3.10) will imply ρ=0, and then Eq (3.14) will imply hβ=0. Note that by Eq (3.6) we have divξ=n(1+σβ), and therefore on the compact hypersurface (Mn,g), we have (cf. [4])

    Mn(1+σβ)=0,

    which does not allow β=0. Hence, in the situation where ρ is a constant, we have h=0, and, also, by Eq (3.12), ζ=0, and by virtue of Eq (3.7), we will reach the conclusion a=0, contrary to the assumption that a is a unit vector. Thus, ρ is a non-constant function which satisfies Obata's equation (3.8) proving that (Mn,g) is isometric to the sphere Sn(c) (cf. [13,14]).

    Conversely, consider the isometric immersion f:Sn(c)(Rn+1,,) of the sphere Sn(c) in the Euclidean space (Rn+1,,) given by f(x)=x. Then, the unit normal ς=cf, the shape operator B=cI, and the mean curvature β=c. Consider the unit constant vector a given by the first Euclidean coordinate vector field, that is,

    a=x1=ζ+hς,

    where ζ is tangent to the sphere Sn(c) and h=a,ς=a,cf=cf,a. Thus, defining ρ=f,a, we have

    h=cρ. (3.19)

    Now, differentiating (3.19) in the direction of XΓ(TSn(c)) and equating the tangential and normal parts, we confirm

    Xζ=chXh=cζ (3.20)

    Using Eqs (3.19) and (3.20), we see cζ=cρ, that is,

    ζ=ρ, (3.21)

    which, in view of the first equation in (3.20), provides

    Δρ=divζ=nch=ncρ. (3.22)

    Finally, using Eq (3.21), the Ricci curvature Ric(ζ,ζ) of the sphere Sn(c) is given by

    Ric(ζ,ζ)=(n1)cζ2=(n1)cρ2. (3.23)

    However, Equations (3.19) and (3.22) confirm

    Sn(c)ρ2dVg=ncSn(c)ρ2dVg=nSn(c)h2dVg=ncSn(c)h2β2dVg.

    Now, integrating Eq (3.23) while using above equation yields

    Sn(c)Ric(ζ,ζ)dVg=n(n1)Sn(c)h2β2dVg.

    Hence, the converse also holds.

    Next, we prove the following result for the complete hypersurface (Mn,g) of the Euclidean space (Rn+1,,).

    Theorem 2. An n-dimensional complete and simply connected isometrically immersed hypersurface f:(Mn,g)(Rn+1,,), n>1, in the Euclidean space (Rn+1,,) with mean curvature β and a constant unit vector a=ζ+hς on Rn+1, where the function h=a,ς0 satisfies Δh=nch for a positive constant c, is isometric to the sphere Sn(c) if and only if the mean curvature β is a constant along the integral curves of ζ and β2c holds.

    Proof. Consider an n-dimensional complete and simply connected Riemannian manifold (Mn,g) that admits an isometric immersion f:(Mn,g)(Rn+1,,) in the Euclidean space (Rn+1,,) such that the function h=a,ς0 satisfies

    Δh=nch, (3.24)

    with mean curvature β satisfying

    ζ(β)=0

    and

    β2c. (3.25)

    We use Eqs (3.5) and (3.8), the symmetry of the shape operator B, and a local orthonormal frame {e1,..,en} in order to find div(Bζ),

    div(Bζ)=nl=1g(elBζ,el)=nl=1g((elB)(ζ)+B(hBel),el)=nl=1g(ζ,(elB)(el))+hB2=nζ(β)+hB2.

    Using Eq (3.15), we get

    div(Bζ)=hB2.

    Now, taking the divergence in Eq (3.9) and using the above equation with Eq (3.24) yields

    hB2=nch,

    that is,

    h2(B2nβ2)=nh2(cβ2). (3.26)

    Combining above equation with inequality (3.25), while keeping in view Cauchy–Schwartz's inequality B2nβ2, we get

    h2(B2nβ2)=0.

    Since h0, we get

    B2=nβ2,

    which, being an inequality in Cauchy–Schwartz's inequality B2nβ2, we must have

    B=βI. (3.27)

    The above equation implies

    (XB)(Y)=X(β)YX,YΓ(TMn),

    which gives

    nl=1(elB)(el)=β.

    Combining the above equation with Eq (3.5) yields

    (n1)β=0

    and, as n>1, we get that β is a constant, and by virtue of Eqs (3.26) and (3.27), we have

    β2=c.

    Now, using Eq (3.27) in the expression of the curvature tensor of the hypersurface with the above equation gives

    R(X,Y)Z=c{g(Y,Z)Xg(X,Z)Y}X,,Y,ZΓ(TMn),

    that is, (Mn,g) is a complete and simply connected space of constant positive curvature c. Hence, (Mn,g) is isometric to Sn(c). The converse is trivial.

    Consider an n-dimensional Riemannian manifold (Mn,g) that possesses a concircular vector field ξ (cf. [3]), that is, the vector field satisfies

    Xξ=σXXΓ(TMn), (4.1)

    where σ is a smooth function, called the potential function of the concircular vector field. A concircular vector field is said to be non-trivial if the potential function σ0. Using Eq (4.1), we immediately have

    divξ=nσ. (4.2)

    In this section, we are interested in an n-dimensional compact Riemannian manifold (Mn,g) that possesses a non-trivial concircular vector field ξ with potential function σ satisfying

    Δσ=ncσ,

    where c>0 is a constant, that is, σ is an eigenfunction of the Laplace operator with eigenvalue the same as the first non-zero eigenvalue of the sphere Sn(c), and we find a condition under which (Mn,g) is isometric to the sphere Sn(c).

    Before we approach this issue, we first prepare an auxiliary result to prove the main result. First, for Eq (4.1), using Eq (2.1) immediately gives the following expression of the curvature tensor, namely

    R(X,Y)ξ=X(σ)YY(σ)XX,YΓ(TMn).

    Taking the trace in the above equation and using Eq (2.2), we reach

    Ric(X,ξ)=(n1)X(σ) (4.3)

    and this equation gives the following expression for the Ricci operator Q operating on ξ, namely

    Q(ξ)=(n1)σ,

    where σ is the gradient of the potential function σ.

    In the following paragraph, we show that for each concircular vector field ξ on a connected Riemannian manifold (Mn,g) there corresponds a smooth function f, which we call a concircular function of the concircular vector field ξ. Note that Eq (4.2) implies

    R(X,ξ)ξ=X(σ)ξξ(σ)XXΓ(TMn)

    and the operator R(X,ξ)ξ, XΓ(TMn) is symmetric in X, and, therefore, the above equation implies

    X(σ)g(ξ,Y)=Y(σ)g(ξ,X)XΓ(TMn).

    The above equation implies

    X(σ)ξ=g(ξ,X)σ

    and taking the inner product in the above equation with σ and replacing X by ξ, we conclude

    (ξ(σ))2=ξ2σ2,

    that is

    ξ2σ2=g(ξ,σ)2.

    This proves that vector fields σ and ξ are parallel, and, consequently, there exists a smooth function f such that

    σ=fξ.

    We call this function f the concircular function of the concircular vector field ξ.

    First, we prove the following proposition.

    Proposition 1. Let ξ be a non-trivial concircular vector field on an n-dimensional compact Riemannian manifold (Mn,g) with potential function σ and concircular function f. If the potential function σ satisfies

    Δσ=ncσ

    for a positive constant c, then

    Mn(Hσ21n(Δσ)2)dVg=n1nMn(ξ(f))2dVg.

    Proof. Let ξ be a non-trivial concircular vector field on an n-dimensional compact Riemannian manifold (Mn,g) with potential function σ and concircular function f, and the potential function σ satisfies

    Δσ=ncσ (4.4)

    for a positive constant c. Using Eqs (4.2) and (4.3), we have

    Ric(ξ,ξ)=(n1)ξ(σ)=(n1)[div(σξ)nσ2]

    and integrating the above equation, we confirm

    MnRic(ξ,ξ)dVg=n(n1)Mnσ2dVg.

    Using the integral formula in [16], we have for a vector field ζ on (Mn,g)

    Mn(Ric(ζ,ζ)+12|£ζg|2ζ2(divζ)2)dVg=0.

    Replacing ζ in the above equation by σ and noting that

    12|£ζg|2=2Hσ2ζ2=Hσ2,

    we conclude

    Mn(Ric(σ,σ)+Hσ2(Δσ)2)dVg=0.

    Thus, we have

    Mn(Hσ21n(Δσ)2)dVg=Mn(n1n(Δσ)2Ric(σ,σ))dVg

    and, inserting Eq (4.4), we reach

    Mn(Hσ21n(Δσ)2)dVg=Mn(n(n1)c2σ2Ric(σ,σ))dVg.

    Now, inserting from Eq (4.4) in the above equation takes us to

    Mn(Hσ21n(Δσ)2)dVg=Mn(n(n1)c2σ2f2Ric(ξ,ξ))dVg. (4.5)

    Using Eqs (4.3) and (4.4), we have

    Ric(ξ,ξ)=(n1)ξ(σ)=(n1)g(ξ,σ)=(n1)fξ2

    and using this in Eq (4.5) leads to

    Mn(Hσ21n(Δσ)2)dVg=(n1)Mn(nc2σ2+f3ξ2)dVg. (4.6)

    Note that taking the divergence on both sides of Eq (4.4) and using Eq (4.2) gives

    Δσ=ξ(f)+nfσ

    and combining this with Eq (4.4) allows us to conclude

    ξ(f)=nσ(c+f). (4.7)

    Using Eqs (4.2), (4.4), and (4.7), we compute

    div(f2σξ)=ξ(f2σ)+nf2σ2=f2g(ξ,σ)+2σfξ(f)+nf2σ2=f3ξ22nσ2f(c+f)+nf2σ2

    and, integrating the above equation, gives

    Mnf3ξ2dVg=Mn(nf2σ2+2ncfσ2)dVg.

    Inserting the above equation into Eq (4.6) leads to

    Mn(Hσ21n(Δσ)2)dVg=n(n1)Mnσ2(c+f)2dVg

    and combining it with Eq (4.7) yields

    Mn(Hσ21n(Δσ)2)dVg=(n1)nMn(ξ(f))2dVg.

    As a straightforward application of the above result, we have the following theorem.

    Theorem 3. An n-dimensional compact and connected Riemannian manifold (Mn,g) admits a non-trivial concircular vector field ξ with potential function σ and concircular function f such that the potential function σ satisfies

    Δσ=ncσ

    for a positive constant c, and the concircular function f is a constant along the integral curves of ξ if and only if (Mn,g) is isometric to Sn(c).

    Proof. Suppose an n-dimensional compact and connected Riemannian manifold (Mn,g) admits a non-trivial concircular vector field ξ with potential function σ and concircular function f such that the potential function σ satisfies

    Δσ=ncσ (4.8)

    for a positive constant c, and the concircular function f is a constant along the integral curves of ξ. Then, by Proposition 1 we have

    Mn(Hσ21n(Δσ)2)dVg=0. (4.9)

    The Cauchy–Schwartz inequality implies

    Hσ21n(Δσ)2 (4.10)

    and equality holds if and only if

    Hσ=ΔσnI. (4.11)

    In view of inequality (4.10) and Eq (4.9), we are ready to conclude the equality

    Hσ2=1n(Δσ)2

    and, therefore, Eq (4.11) holds. Combining Eqs (4.8) and (4.11), we arrive at

    Hσ=cσI. (4.12)

    Note that the potential function σ can not be a constant, for if it were a constant, the above equation would imply σ=0, which is contrary to the assumption that ξ is a non-trivial concircular vector field. Thus, Equation (4.12) is Obata's equation, and therefore (Mn,g) is isometric to the sphere Sn(c) (cf. [13,14]).

    Conversely, take a constant unit vector a on the Euclidean space (Rn+1,,) while treating Sn(c) as a hypersurface of (Rn+1,,) with unit normal ς, shape operator B=cI, and expressing a as

    a=ξ+hς, (4.13)

    where ξ is tangent to Sn(c) and h=a,ς. Differentiating equation (4.13) with respect to XΓ(TSn(c)) and equating the tangential and normal parts, we arrive at

    Xξ=chXh=cξ. (4.14)

    This confirms that ξ is a concircular vector field on Sn(c) with potential function σ=ch, and the second equation gives σ=cξ. This proves that the concircular function f=c. Moreover, if the potential function σ=0, we get h=0, and by the second equation in (4.14), we get ξ=0. In this case, Eq (4.13) confirms a=0, a contradiction to the fact that a is a unit vector. Thus, the potential function σ0, that is, the concircular vector field ξ on Sn(c) is non-trivial. Note that divξ=nch=nσ, and, combining it with the equation σ=cξ, we get Δσ=ncσ with c a positive constant. Hence, the converse holds.

    We have initiated the study of an n-dimensional compact Riemannian manifold (Mn,g) that has an eigenvalue nc for a positive constant c of the Laplace operator the same as the first non-zero eigenvalue of the n-sphere Sn(c) of constant curvature c, and searched for an additional condition under which (Mn,g) is isometric to the sphere Sn(c). The main aim was to find an appropriate smooth function on (Mn,g) that will become the eigenfunction of the Laplace operator with eigenvalue nc as seen in Theorems 1 and 3. Naturally, the scope of this study is quite modest, for instance one can consider an n-dimensional compact Riemannian manifold (Mn,g) that admits a torse forming vector field ξ (cf. [17]). Recall that a torse forming vector field ξ on (Mn,g) satisfies

    Xξ=σX+ω(X)ξXΓ(TMn),

    where σ is a smooth function defined on Mn called the conformal scalar and ω is a smooth 1-form on Mn called the generating form of the a torse forming vector field ξ. It will be an interesting question to consider torse forming vector field ξ on an n-dimensional compact Riemannian manifold (Mn,g) such that its conformal scalar σ satisfies Δσ=ncσ for a positive constant c, and find conditions under which (Mn,g) is isometric to Sn(c).

    We know that the second non-zero eigenvalue of the sphere Sn(c) is given by λ2=2(n+1)c, and another aspect of our work could be, if there is a smooth function f on an n-dimensional compact Riemannian manifold (Mn,g) such that Δf=2(n+1)cf, that is, (Mn,g) has an eigenvalue same as second non-zero eigenvalue of the sphere Sn(c), to find additional condition on (Mn,g) so that (Mn,g) is isometric to Sn(c).

    Sharief Deshmukh: Conceptualization, Methodology, Writing-original draft, Writing-review and editing, Supervision; Amira Ishan: Conceptualization, Methodology, Writing-review and editing; Olga Belova: Formal analysis, Writing-original draft, Writing-review and editing. All authors have read and agreed to the published version of the manuscript.

    The authors would like to acknowledge the Deanship of Graduate Studies and Scientific Research, Taif University for funding this work.

    The authors declare no conflicts of interest in this paper.



    [1] N. Bellomo, A. Bellouquid, Y. Tao, M. Winkler, Toward a mathematical theorey of Keller-Segel models of pattern formation in biological tissues, Math. Mod. Meth. Appl. Sci., 25 (2015), 1663–1763. https://doi.org/10.1142/S021820251550044X doi: 10.1142/S021820251550044X
    [2] P. Biler, Global solutions to some parabolic-elliptic systems of chemotaxis, Adv. Math. Sci. Appl., 9 (1999), 347–359. https://doi.org/10.2307/3857479 doi: 10.2307/3857479
    [3] T. Black, Global generalized solutions to a parabolic-elliptic Keller-Segel system with singular sensitivity, Discrete Cont. Dyn.-S, 13 (2020), 119–137. https://doi.org/10.1016/j.urology.2019.11.020 doi: 10.1016/j.urology.2019.11.020
    [4] S. Chang, P. Yang, Conformal deformation of metrics on S2, J. Differ. Geom., 27 (1988), 259–296. https://doi.org/10.1038/scientificamerican1088-27 doi: 10.1038/scientificamerican1088-27
    [5] X. Cao, Global classical solutions in chemotaxis-Navier-Stokes system with rotational flux term, J. Differ. Eq., 261 (2016), 6883–6914. https://doi.org/10.1016/j.jde.2016.09.007 doi: 10.1016/j.jde.2016.09.007
    [6] X. Cao, J. Lankeit, Global classical small-data solutions for a three-dimensional chemotaxis Navier-Stokes system involving matrix-valued sensitivities, Calc. Var. Partial Diff. Eq., 55 (2016), 55–107. https://doi.org/10.2216/0031-8884-55.1.107 doi: 10.2216/0031-8884-55.1.107
    [7] Y. Chiyo, M. Marras, Y. Tanaka, T. Yokota, Blow-up phenomena in a parabolic-elliptic-elliptic attraction-repulsion chemotaxis system with superlinear logistic degradation, Nonlinear Anal., 212 (2021), 112550. https://doi.org/10.1016/j.na.2021.112550 doi: 10.1016/j.na.2021.112550
    [8] M. DiFrancesco, A. Lorz, P. A. Markowich, Chemotaxis-fluid coupled model for swimming bacteria with nonlinear diffusion: Global existence and asymptotic behavior, Discrete Cont. Dyn.-A, 28 (2010), 1437–1453. https://doi.org/10.1055/s-0029-1218690 doi: 10.1055/s-0029-1218690
    [9] R. Duan, X. Li, Z. Xiang, Global existence and large time behavior for a two-dimensional chemotaxis-Navier-Stokes system, J. Differ. Equations, 263 (2017), 6284–6316. https://doi.org/10.1016/j.jde.2017.07.015 doi: 10.1016/j.jde.2017.07.015
    [10] R. Duan, A. Lorz, P. A. Markowich, Global solutions to the coupled chemotaxis-fluid equations, Commun. Part. Diff. Eq., 35 (2010), 1635–1673. https://doi.org/10.1080/03605302.2010.497199 doi: 10.1080/03605302.2010.497199
    [11] R. Duan, Z. Xiang, A note on global existence for the chemotaxis-Stokes model with nonlinear diffusion, Int. Math. Res. Not., 2014, 1833–1852. http://doi.org/10.1093/imrn/rns270
    [12] M. Fuest, Finite-time blow-up in a two-dimensional Keller-Segel system with an environmental dependent logistic source, Nonlinear Anal. Real World Appl., 52 (2020), 103022. https://doi.org/10.1016/j.nonrwa.2019.103022 doi: 10.1016/j.nonrwa.2019.103022
    [13] K. Fujie, Boundedness in a fully parabolic chemotaxis system with singular sensitivity, J. Math. Anal. Appl., 424 (2015), 675–684. https://doi.org/10.1016/j.jmaa.2014.11.045 doi: 10.1016/j.jmaa.2014.11.045
    [14] K. Fujie, T. Senba, Global existence and boundedness in a parabolic-elliptic Keller-Segel system with general sensitivity, Discrete Cont. Dyn.-B, 21 (2016), 81–102. https://doi.org/10.21714/2179-8834/2016v21n4p81-102 doi: 10.21714/2179-8834/2016v21n4p81-102
    [15] K. Fujie, M. Winkler, T. Yokota, Boundedness of solutions to parabolic-elliptic Keller-Segel systems with signal-dependent sensitivity, Math. Method. Appl. Sci., 38 (2015), 1212–1224. https://doi.org/10.1111/ecog.01398 doi: 10.1111/ecog.01398
    [16] J. J. Neto, J. Claeyssen, Capital-induced labor migration in a spatial solow model, J. Econ., 115 (2015), 25–47. https://doi.org/10.1007/s00712-014-0404-6 doi: 10.1007/s00712-014-0404-6
    [17] J. Juchem Neto, J. Claeyssen, S. Pôrto Júnior, Economic agglomerations and spatio-temporal cycles in a spatial growth model with capital transport cost, Physica A, 494 (2018), 76–86. https://doi.org/10.1016/j.physa.2017.12.036 doi: 10.1016/j.physa.2017.12.036
    [18] J. Juchem Neto, J. Claeyssen, S. Pôrto Júnior, Returns to scale in a spatial Solow-Swan economic model, Physica A, 533 (2019), 122055. https://doi.org/10.1016/j.physa.2019.122055 doi: 10.1016/j.physa.2019.122055
    [19] D. Henry, Geometric theory of semilinear parabolic equations, Springer-Verlag, New York, 1981.
    [20] D. Horstmann, M. Winkler, Boundedness vs. blow-up in a chemotaxis system, J. Differ. Equations, 215 (2005), 52–107.
    [21] Y. Ke, J. Zheng, An optimal result for global existence in a three-dimensional Keller-Segel-Navier-Stokes system involving tensor-valued sensitivity with saturation, Calc. Var. Partial Dif., 58 (2019), 1–27.
    [22] E. F. Keller, L. A. Segel, Initiation of slime model aggregation viewed as an instability, J. Theor. Biol., 26 (1970), 399–415. https://doi.org/10.1016/0022-5193(70)90092-5 doi: 10.1016/0022-5193(70)90092-5
    [23] O. A. Ladyzenskaya, V. A. Solonnikov, N. N. Ural'ceva, Linear and quasi-linear equations of parabolic type, Amer. Math. Soc. Trans., Providence, 1968.
    [24] J. Lankeit, A new approach toward boundedness in a two-dimensional parabolic chemotaxis system with singular sensitivity, Math. Method. Appl. Sci., 39 (2016), 394–404. https://doi.org/10.1016/j.tins.2016.04.002 doi: 10.1016/j.tins.2016.04.002
    [25] X. Li, Global classical solutions in a Keller-Segel(-Navier)-Stokes system modeling coral fertilization, J. Differ. Equations, 267 (2019), 6290–6315. https://doi.org/10.1016/j.jde.2019.06.021 doi: 10.1016/j.jde.2019.06.021
    [26] B. Li, Y. Li, On a chemotaxis-type Solow-Swan model for economic growth with capital-induced labor migration, J. Math. Anal. Appl., 511 (2022), 126080. https://doi.org/10.1016/j.jmaa.2022.126080 doi: 10.1016/j.jmaa.2022.126080
    [27] M. Li, Z. Xiang, G. Zhou, The stability analysis of a 2D Keller-Segel-Navier-Stokes system in fast signal diffusion, Eur. J. Appl. Math., 34 (2022), 160–209. http://doi.org/10.1017/S0956792522000067 doi: 10.1017/S0956792522000067
    [28] K. Lin, C. Mu, L. Wang, Large-time behavior of an attraction-repulsion chemotaxis system, J. Math. Anal. Appl., 426 (2015), 105–124. https://doi.org/10.1016/j.jmaa.2014.12.052 doi: 10.1016/j.jmaa.2014.12.052
    [29] F. Dai, B. Liu, Boundedness and asymptotic behavior in a Keller-Segel(-Navier) system with indirect signal production, J. Differ. Equations, 314 (2022), 201–250. https://doi.org/10.1016/j.jde.2022.01.015 doi: 10.1016/j.jde.2022.01.015
    [30] F. Dai, B. Liu, Global weak solutions in a three-dimensional Keller-Segel-Navier-Stokes system with indirect signal production, J. Differ. Equations, 333 (2022), 436–488. https://doi.org/10.1016/j.jde.2022.06.015 doi: 10.1016/j.jde.2022.06.015
    [31] J. Liu, Y. Wang, Global weak solutions in a three-dimensional Keller-Segel-Navier-Stokes system involving a tensor-valued sensitivity with saturation, J. Differ. Equations, 262 (2017), 5271–5305. https://doi.org/10.1016/j.jde.2017.01.024 doi: 10.1016/j.jde.2017.01.024
    [32] S. Liu, L. Wang, Global boundedness of a chemotaxis model with logistic growth and general indirect signal production, J. Math. Anal. Appl., 505 (2022), 125613. https://doi.org/10.1016/j.jmaa.2021.125613 doi: 10.1016/j.jmaa.2021.125613
    [33] X. Liu, Y. Zhang, Y. Han, Small-data solutions of chemotaixs-fluid system with indirect signal production, J. Math. Anal. Appl., 508 (2022), 125908. https://doi.org/10.1016/j.jmaa.2021.125613 doi: 10.1016/j.jmaa.2021.125613
    [34] X. Liu, J. Zheng, Convergence rates of solutions in apredator-preysystem with indirect pursuit-evasion interaction in domains of arbitrary dimension, Discrete Cont. Dyn.-B, 28 (2023), 2269–2293. https://doi.org/10.3934/dcdsb.2022168 doi: 10.3934/dcdsb.2022168
    [35] N. Mizoguchi, P. Souplet, Nondegeneracy of blow-up points for the parabolic Keller-Segel system, Ann. I. H. Poincaré-An., 31 (2014), 851–875. https://doi.org/10.1111/1911-3846.12048 doi: 10.1111/1911-3846.12048
    [36] M. Mizukami, T. Yokota, A unified method for boundedness in fully parabolic chemotaxis systems with signal-dependent sensitivity, Math. Nachr., 290 (2017), 2648–2660. https://doi.org/10.1002/mana.201600399 doi: 10.1002/mana.201600399
    [37] T. Nagai, T. Senba, Global existence and blow-up of radial solutions to a parabolic-elliptic system of chemotaxis, Adv. Math. Sci. Appl., 8 (1998), 145–156. https://doi.org/10.1016/S0030-4018(98)00425-8 doi: 10.1016/S0030-4018(98)00425-8
    [38] T. Nagai, T. Senba, K. Yoshida, Application of the Trudinger-Moser inequality to a parabolic system of chemotaxis, Funkc. Ekvac., 40 (1997), 411–433. https://doi.org/10.1016/S0304-3932(97)00048-2 doi: 10.1016/S0304-3932(97)00048-2
    [39] K. Osaki, T. Tsujikawa, A. Yagi, M. Mimura, Exponential attractor for a chemotaxis-growth system of equations, Nonlinear Anal.-Theor., 51 (2002), 119–144. https://doi.org/10.1016/S0362-546X(01)00815-X doi: 10.1016/S0362-546X(01)00815-X
    [40] M. M. Porzio, V. Vespri, Hölder estimates for local solutions of some doubly nonlinear degenerate parabolic equation, J. Differ. Equations, 103 (1993), 146–178. https://doi.org/10.1006/jdeq.1993.1045 doi: 10.1006/jdeq.1993.1045
    [41] Y. Peng, Z. Xiang, Global existence and boundedness in a 3D Keller-Segel-Stokes system with nonlinear diffusion and rotational flux, Z. Angew. Math. Phys., 68 (2017), 68.
    [42] Y. Peng, Z. Xiang, Global existence and convergence rates to a chemotaxis-fluids system with mixed boundary conditions, J. Differ. Equations, 267 (2019), 1277–1321. https://doi.org/10.1016/j.jde.2019.02.007 doi: 10.1016/j.jde.2019.02.007
    [43] Y. Peng, Z. Xiang, Global solution to the coupled Chemotaxis-Fluids system in a 3D unbounded domain with boundary, Math. Mod. Meth. Appl. Sci., 28 (2018), 869–920. https://doi.org/10.1142/S0218202518500239 doi: 10.1142/S0218202518500239
    [44] Y. Shen, Preliminary of global differential geometry, 3 Eds., Higher Eduction Press, Bei Jing, 2009.
    [45] V. A. Solonnikov, Schauder estimate for the evolutionary generalized Stokes problem, In: Nonlinear Equations and Spectral Theory, Providence, Rhode Island, 2007,165–200.
    [46] Y. Tao, Z. Wang, Competing effects of attraction vs. repulsion in chemotaxis, Math. Mod. Meth. Appl. Sci., 23 (2013), 1–36.
    [47] Y. Tao, M. Winkler, Blow-up prevension by quadratic degradation in a two-dimensional Keller-Segel-Navier-Stokes system, Z. Angew. Math. Phys., 67 (2016), 138.
    [48] Y. Tao, M. Winkler, Global existence and boundedness in a Keller-Segel-Stokes model with arbitrary porous medium diffusion, Discrete Cont. Dyn.-A, 32 (2012), 1901–1914.
    [49] I. Tuval, L. Cisneros, C. Dombrowski, C. W. Wolgemuth, J. O. Kessler, R. E. Goldstein, Bacterial swimming and oxygen transport near contact line, Proc. Natl. Acad. Sci. USA, 102 (2005), 2277–2282. https://doi.org/10.1073/pnas.0406724102 doi: 10.1073/pnas.0406724102
    [50] N. Trudinger, On embeddings into Orlicz spaces and some applications, J. Math. Mech., 17 (1967), 473–483. https://doi.org/10.1512/iumj.1968.17.17028 doi: 10.1512/iumj.1968.17.17028
    [51] Y. Wang, X. Cao, Global classical solutions of a 3D chemotaxis-Stokes system with rotation, Discrete Cont. Dyn.-B, 20 (2015), 3235–3254. https://doi.org/10.3934/dcdsb.2015.20.3235 doi: 10.3934/dcdsb.2015.20.3235
    [52] Y. Wang, M. Winkler, Z. Xiang, Global classical solutions in a two-dimensional chemotaxis Navier-Stokes system with subcritical sensitivity, Ann. Sci. Norm.-Sci., 18 (2018), 421–466.
    [53] Y. Wang, M. Winkler, Z. Xiang, Global mass-preserving solutions to a chemotaxis-fluid model involving Dirichlet boundary conditions for the signal, Anal. Appl., 20 (2022), 141–170.
    [54] Y. Wang, M. Winkler, Z. Xiang, Global solvability in a three-dimensional Keller-Segel-Stokes system involving arbitrary superlinear logistic degradation, Adv. Nonlinear Anal., 10 (2021), 707–731. https://doi.org/10.1002/pchj.457 doi: 10.1002/pchj.457
    [55] Y. Wang, M. Winkler, Z. Xiang, Immediate regularization of measure-type population densities in a two-dimensional chemotaxis system with signal consumption, Sci. China Math., 64 (2021), 725–746. https://doi.org/10.1007/s11425-020-1708-0 doi: 10.1007/s11425-020-1708-0
    [56] Y. Wang, M. Winkler, Z. Xiang, Local energy estimates and global solvability in a three-dimensional chemotaxis-fluid system with prescribed signal on the boundary, Commun. Part. Diff. Eq., 46 (2021), 1058–1091. https://doi.org/10.1080/03605302.2020.1870236 doi: 10.1080/03605302.2020.1870236
    [57] Y. Wang, M. Winkler, Z. Xiang, The fast signal diffusion limit in Keller-Segel(-fluid) systems, Calc. Var. Partial Dif., 58 (2019), 196. https://doi.org/10.1007/s00526-019-1656-3 doi: 10.1007/s00526-019-1656-3
    [58] Y. Wang, Z. Xiang, Global existence and boundedness in a Keller-Segel-Stokes system involving a tensor-valued sensitivity with saturation, J. Differ. Equations, 259 (2015), 7578–7609. https://doi.org/10.1016/j.jde.2015.08.027 doi: 10.1016/j.jde.2015.08.027
    [59] Y. Wang, Z. Xiang, Global existence and boundedness in a Keller-Segel-Stokes system involving a tensor-valued sensitivity with saturation: The 3D case, J. Differ. Equations, 261 (2016), 4944–4973. https://doi.org/10.1016/j.jde.2016.07.010 doi: 10.1016/j.jde.2016.07.010
    [60] Y. Wang, L. Yang, Boundedness in a chemotaxis-fluid system involving a saturated sensitivity and indirect signal production mechanism, J. Differ. Equations, 287 (2021), 460–490. https://doi.org/10.1016/j.jde.2021.04.001 doi: 10.1016/j.jde.2021.04.001
    [61] M. Winkler, Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model, J. Differ. Equations, 248 (2010), 2889–2905. https://doi.org/10.1016/j.jde.2010.02.008 doi: 10.1016/j.jde.2010.02.008
    [62] M. Winkler, A three-dimensional Keller-Segel-Navier-Stokes system with logistic source: Global weak solutions and asymptotic stabilization, J. Funct. Anal., 276 (2019), 1339–1401. https://doi.org/10.1016/j.jfa.2018.12.009 doi: 10.1016/j.jfa.2018.12.009
    [63] M. Winkler, Global large-data solutions in a chemotaxis-(Navier-) Stokes system modeling cellular swimming in fluid drops, Commun. Part. Diff. Eq., 37 (2012), 319–351. https://doi.org/10.1080/03605302.2011.591865 doi: 10.1080/03605302.2011.591865
    [64] M. Winkler, Global mass-preserving solutions in a two-dimensional chemotaxis-Stokes system with rotation flux components, J. Evol. Equ., 18 (2018), 1267–1289. https://doi.org/10.1007/s00028-018-0440-8 doi: 10.1007/s00028-018-0440-8
    [65] M. Winkler, Global solutions in a fully parabolic chemotaxis system with singular sensitivity, Math. Method. Appl. Sci., 34 (2011), 176–190. https://doi.org/10.1002/mma.1346 doi: 10.1002/mma.1346
    [66] M. Winkler, Global weak solutions in a three-dimensional chemotaxis-Navier-Stokes syste, Ann. I. H. Poincaré-An., 33 (2016), 1329–1352. https://doi.org/10.1016/j.anihpc.2015.05.002 doi: 10.1016/j.anihpc.2015.05.002
    [67] M. Winkler, How far do chemotaxis-driven forces influence regularity in the Navier-Stokes system? T. Am. Math. Soc., 369 (2017), 3067–3125. https://doi.org/10.1090/tran/6733
    [68] M. Winkler, Small-mass solutions in the two-dimensional Keller-Segel system coupled to the Navier-Stokes equations, SIAM J. Math. Anal., 52 (2020), 2041–2080. https://doi.org/10.1137/19M1264199 doi: 10.1137/19M1264199
    [69] M. Winkler, Stabilization in a two-dimensional chemotaxis-Navier-Stokes system, Arch. Ration. Mech. Anal., 211 (2014), 455–487. https://doi.org/10.1007/s00205-013-0678-9 doi: 10.1007/s00205-013-0678-9
    [70] J. Wu, H. Natal, Boundedness and asymptotic behavior to a chemotaxis-fluid system with singular sensitivity and logistic source, J. Math. Anal. Appl., 484 (2020), 123748. https://doi.org/10.1016/j.jmaa.2019.123748 doi: 10.1016/j.jmaa.2019.123748
    [71] J. Wu, C. Wu, A note on the global existence of a two-dimensional chemotaxis-Navier-Stokes system, Appl. Anal., 98 (2019), 1224–1235. https://doi.org/10.1080/00036811.2017.1419199 doi: 10.1080/00036811.2017.1419199
    [72] P. Yu, Blow up prevention by saturated chemotaxis sensitivity in a 2D Keller-Segel-Stokes system, Acta Appl. Math., 169 (2020), 475–497. https://doi.org/10.1007/s10440-019-00307-8 doi: 10.1007/s10440-019-00307-8
    [73] Q. Zhang, Y. Li, Global weak solutions for the three-dimensional chemotaxis-Navier-Stokes system with nonlinear diffusion, J. Differ. Equations, 259 (2015), 3730–3754. https://doi.org/10.1016/j.jde.2015.05.012 doi: 10.1016/j.jde.2015.05.012
    [74] W. Zhang, P. Niu, S. Liu, Large time behavior in a chemotaxis model with logistic growth and indirect signal production, Nonlinear Anal., Real Word Appl., 50 (2019), 484–497. https://doi.org/10.1016/j.nonrwa.2019.05.002 doi: 10.1016/j.nonrwa.2019.05.002
    [75] X. Zhao, S. Zheng, Global boundedness of solutions in a parabolic-parabolic chemotaxis system with singular sensitivity, J. Math. Anal. Appl., 443 (2016), 445–452. https://doi.org/10.1016/j.jmaa.2016.05.036 doi: 10.1016/j.jmaa.2016.05.036
    [76] J. Zheng, A new result for the global existence (and boundedness) and regularity of a three-dimensional Keller-Segel-Navier-Stokes system modeling coral fertilization, J. Differ. Equations, 272 (2021), 164–202. https://doi.org/10.1016/j.jde.2020.09.029 doi: 10.1016/j.jde.2020.09.029
    [77] J. Zheng, An optimal result for global existence and boundedness in a three-dimensional Keller-Segel-Stokes system with nonlinear diffusion, J. Differ. Equations, 267 (2019), 2385–2415. https://doi.org/10.1016/j.jde.2019.03.013 doi: 10.1016/j.jde.2019.03.013
    [78] J. Zheng, Boundedness of solutions to a quasilinear parabolic-elliptic Keller-Segel system with logistic source, J. Differ. Equations, 259 (2015), 120–140. https://doi.org/10.1016/j.jde.2015.02.003 doi: 10.1016/j.jde.2015.02.003
    [79] J. Zheng, Eventual smoothness and stabilization in a three-dimensional Keller-Segel-Navier-Stokes system with rotational flux, Calc. Var. Partial Dif., 61 (2022), 52. https://doi.org/10.1007/s00526-021-02164-6 doi: 10.1007/s00526-021-02164-6
  • 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(2001) PDF downloads(126) Cited by(31)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog