Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
Research article

Simulations of vehicle-induced mixing and near-road aerosol microphysics using computational fluid dynamics

  • Understanding the fate of ultrafine particles (UFP), especially from combustion sources, is essential to assess their impact on health and climate. Here, we present simulations of the behavior of UFP in the near-roadway environment (up to 300 m downwind) based on a model with coupled computational fluid dynamics (CFD) and aerosol microphysics. It is found that vehicle-induced mixing (VIM) caused by the combined effect of vehicle wake formation and production of turbulent kinetic energy plays an important role in downwind dilution of pollutants. Various methodologies for simulating VIM are explored, and a computationally efficient approach based on an effective roughness length for vehicle-induced mixing is proposed. Whereas the coagulation behavior of ultrafine particles is relatively well understood, condensation and/or evaporation can have equally large or larger impacts on the number and sizes of particles downwind of a roadway. Through a set of sensitivity simulations, we show that the particle losses are potentially significant via evaporation but depend strongly on several parameters or processes that are poorly understood and difficult to fully constrain for on-road traffic using measurements: the volatility distribution of organic species and gas-phase concentrations of semi-volatile organics.

    Citation: Satbir Singh, Peter J. Adams, Albert A. Presto. Simulations of vehicle-induced mixing and near-road aerosol microphysics using computational fluid dynamics[J]. AIMS Environmental Science, 2018, 5(5): 315-339. doi: 10.3934/environsci.2018.5.315

    Related Papers:

    [1] Yitian Wang, Xiaoping Liu, Yuxuan Chen . Semilinear pseudo-parabolic equations on manifolds with conical singularities. Electronic Research Archive, 2021, 29(6): 3687-3720. doi: 10.3934/era.2021057
    [2] Yaning Li, Yuting Yang . The critical exponents for a semilinear fractional pseudo-parabolic equation with nonlinear memory in a bounded domain. Electronic Research Archive, 2023, 31(5): 2555-2567. doi: 10.3934/era.2023129
    [3] Mingyou Zhang, Qingsong Zhao, Yu Liu, Wenke Li . Finite time blow-up and global existence of solutions for semilinear parabolic equations with nonlinear dynamical boundary condition. Electronic Research Archive, 2020, 28(1): 369-381. doi: 10.3934/era.2020021
    [4] Jun Zhou . Initial boundary value problem for a inhomogeneous pseudo-parabolic equation. Electronic Research Archive, 2020, 28(1): 67-90. doi: 10.3934/era.2020005
    [5] Qianqian Zhu, Yaojun Ye, Shuting Chang . Blow-up upper and lower bounds for solutions of a class of higher order nonlinear pseudo-parabolic equations. Electronic Research Archive, 2024, 32(2): 945-961. doi: 10.3934/era.2024046
    [6] Fanqi Zeng, Wenli Geng, Ke An Liu, Boya Wang . Differential Harnack estimates for the semilinear parabolic equation with three exponents on Rn. Electronic Research Archive, 2025, 33(1): 142-157. doi: 10.3934/era.2025008
    [7] Yang Cao, Qiuting Zhao . Initial boundary value problem of a class of mixed pseudo-parabolic Kirchhoff equations. Electronic Research Archive, 2021, 29(6): 3833-3851. doi: 10.3934/era.2021064
    [8] Milena Dimova, Natalia Kolkovska, Nikolai Kutev . Global behavior of the solutions to nonlinear Klein-Gordon equation with critical initial energy. Electronic Research Archive, 2020, 28(2): 671-689. doi: 10.3934/era.2020035
    [9] Lianbing She, Nan Liu, Xin Li, Renhai Wang . Three types of weak pullback attractors for lattice pseudo-parabolic equations driven by locally Lipschitz noise. Electronic Research Archive, 2021, 29(5): 3097-3119. doi: 10.3934/era.2021028
    [10] Hui Yang, Futao Ma, Wenjie Gao, Yuzhu Han . Blow-up properties of solutions to a class of p-Kirchhoff evolution equations. Electronic Research Archive, 2022, 30(7): 2663-2680. doi: 10.3934/era.2022136
  • Understanding the fate of ultrafine particles (UFP), especially from combustion sources, is essential to assess their impact on health and climate. Here, we present simulations of the behavior of UFP in the near-roadway environment (up to 300 m downwind) based on a model with coupled computational fluid dynamics (CFD) and aerosol microphysics. It is found that vehicle-induced mixing (VIM) caused by the combined effect of vehicle wake formation and production of turbulent kinetic energy plays an important role in downwind dilution of pollutants. Various methodologies for simulating VIM are explored, and a computationally efficient approach based on an effective roughness length for vehicle-induced mixing is proposed. Whereas the coagulation behavior of ultrafine particles is relatively well understood, condensation and/or evaporation can have equally large or larger impacts on the number and sizes of particles downwind of a roadway. Through a set of sensitivity simulations, we show that the particle losses are potentially significant via evaporation but depend strongly on several parameters or processes that are poorly understood and difficult to fully constrain for on-road traffic using measurements: the volatility distribution of organic species and gas-phase concentrations of semi-volatile organics.


    In this paper, we consider the following initial boundary value problem of semilinear pseudo-parabolic equations with conical degeneration

    utΔButΔBu=|u|p1u, (xb,˜x)intB,t>0, (1.1)
    u(xb,˜x,0)=u0,                 (xb,˜x)intB, (1.2)
    u(0,˜x,t)=0,                    (0,˜x)B,t0, (1.3)

    where u0H1,n22,0(B), 1<p<n+2n2 and n=l+12 is the dimension of B, lN. Here the domain B=[0,1)×X is regarded as the local model near the conical singularity on conical singular manifolds, where XRl is a closed compact C manifold. Denoting the interior of B by intB and the boundary of B by B:={0}×X. We use the coordinates (xb,˜x):=(xb,x1,x2,...,xl)B for 0xb<1, ˜xX near B. The conical Laplacian operator is defined as

    ΔB=2B=(xbxb)2+2x1+...+2xl,

    which is the totally characteristic degeneracy operators on a stretched conical manifold, and B=(xbxb,x1,...,xl) denotes the corresponding gradient operator with conical degeneracy on the boundary B. In particular, we intend to investigate problem (1.1)-(1.3) in the weighted Mellin-Sobolev spaces H1,n22,0(B), and the definition of such distribution spaces will be introduced in Section 2.

    The classical pseudo-parabolic equation

    utΔutΔu=|u|p1u, xΩ,t>0, (1.4)

    defined on a bounded domain ΩRn with smooth boundary appeared in various physical and biological phenomena. For example, taking u as the flow velocity, the homogeneous form of the model equation (1.4) was introduced to study the incompressible simple fluids with fading memory and the non-steady flow with the Rivlin-Ericksen tensors [47,16]. One can also know more about other applications by referring to [23,24,5].

    It is well known that equation (1.4) in the domains contained in classical Euclidean space with regular boundary has been well investigated. Cao et al. [9] considered the Cauchy problem of following model

    tukΔut=Δu+up, xRn,t>0,

    and obtained the critical global existence exponent and the critical Fujita exponent by integral representation and contraction mapping principle. Subsequently, its uniqueness was proved by Khomrutai in [27] for the case 0<p<1. Furthermore, Li and Du [32] considered the Cauchy problem of

    utkΔut=Δu+|x|σup, xRn,t>0,

    and achieved the global existence and blowup in finite time of solutions with the critical Fujita exponent and the second critical exponent respectively. They also showed that the inhomogeneous term |x|σ affects the decay asymptotic behavior of solutions and accelerates the blowup of solutions. Khomrutai [28] studied the Cauchy problem of sublinear pseudo-parabolic equation

    tuΔtu=Δu+V(x,t)up, xRn,t>0,

    where V(x,t)λ(t)|x|σ is a non-autonomous and unbounded potential function with 0<p<1, and established the global existence of solutions by approximation and monotonicity argument. They also derived the precise grow-up rate of solutions and critical growth exponent. In order to figure out the effects of small perturbation on the dynamical of diffusion and reaction, Cao and Yin [8] considered the following Cauchy problem

    tukΔut=Δu+up+f(x), xRn,t>0,

    and revealed that small perturbation may develop large variation of solutions as time evolves. We also recommend that the reader refer to [43] to learn more about the effects of the power index of nonlinearity on the dynamical behavior of the solution. Different from above studies that focus on the influence of the nonlinearities especially the power index on the global well-posedness of the solution, [50,49,36] comprehensively studied equation (1.4) by considering the influences of the initial data on the global well-posedness and corresponding properties of solution. Depending on the potential well depth, they classified the initial data to subcritical initial energy level J(u0)<d, critical initial energy level J(u0)=d and supcritical initial energy level J(u0)>d, and proved the global existence, finite time blow up and asymptotic behavior of solutions with J(u0)d. Moreover, thanks to the comparison principle, the global existence and nonexistence of solutions were also obtained at supcritical initial energy level J(u0)>d. When the nonlinear effects are dominated by the logarithmic term, Chen [14] investigated the following nonlinear pseudo-parabolic equation, i.e.,

    tuΔuΔtu=uln|u|, xΩ,t>0

    and proved the global existence and the finite time blow up of solutions under the subcritical and critical initial energy case, respectively. Focusing on the high initial energy level, Xu and Wang et al [51] studied the problem proposed in [50] and gave a sufficient condition on initial data leading to blow up in finite time by the potential well method, at the same time, they also estimated the upper bound of the blowup time. As an important method to reveal the influence of initial data on the dynamical behavior of solutions, the potential well theory can be applied not only to the study of the problem of parabolic equations, but also to the study of the problem for various types of nonlinear evolution equations or systems. Xu and Lian et al [48] investigated the global well-posedness of solutions for coupled parabolic systems in the variational framework, and the initial data leading to the global existence or finite time blow up of the solution are divided. Chen and Xu [15] considered a class of damped fourth-order nonlinear wave equations with logarithmic sources. By examining the effect of weak nonlinear sources on the blow up of the solution, they revealed the confrontation mechanism between the damping structure and the nonlinear source and found the initial data that caused the solution to blow up in infinite time. For related results of polynomial nonlinear sources, we refer to [52]. Furthermore, we suggest the readers refer to [53] for the study of high order nonlinear wave equations, [54,34] for the study of damped nonlinear wave equation problems using improved potential well methods at high initial energy levels, and so on, which are representative recent results, of course we can not list all of the results obtained by the potential well theory here due to the huge amount.

    Actually, geometric singularities have attracted considerable interest and have become the focus of extensive physical and mathematical research in recent years. To find static solutions of Einstein's equations coupled to brane sources, Michele [39] studied the generalizations of the so-called "football" shaped extra dimensions scenario to include two codimension branes, which can be transformed into the mathematical problem of solving the Liouville equation with singularities, where the function space he constructed can be described as a sphere with conical singularities at the brane locations. After that some cone solitons (in the case of compact surfaces) were found by Hamilton [25], where cone singularities also arise naturally on the study of such kind of solitons. Not only being widely applied in cosmology and physics, the cone singular manifold itself also brings a lot of interesting topics to pure mathematics, such as the analytic proof of the cobordism theorem [31]. Conical singularities become a hotspot mainly for reasons of two aspects. Firstly, a manifold with conical singularities is one of the most fundamental stratified spaces and the investigation on it is motivated by the desire of understanding the dynamic behavior of the solution of nonlinear evolution equations on such stratified space. Topologically these spaces are of iterated cone type, in which, due to the conical singularity, the classical differential operator cannot be applied to such manifolds. Secondly, the methods developed for the domains with smooth boundaries cannot be directly applied to domains with singularities. It is a challenge and also an interesting problem in the community to restitute the conclusions established on the smooth domain for the problems defined on the conical space.

    Inspired by above, it is natural to bring some ideas and develop techniques to establish a comprehensive understanding of operator theory on the manifolds with conical singularities, which was first explored by Kondrat'ev in [29] by introducing the celebrated Mellin-Sobolev spaces Hs,γp(B) as the work space for the partial differential equations, then the theory on the conically degenerate pseudo-differential operator and the weighted Sobolev space on the conical manifolds were summarized by Schulze and Egorov in [45] and [18]. With the development of the study related to the singularity problem, including the study on the partial differential equation on manifolds with conical singularities, there occurs a large number of related results about various kinds of initial boundary value problems for evolution equations. For example, when using the porous medium equation to describe the flow of a substance in a porous medium material, the medium usually shows various irregular shapes in different regions, among which the most essential case is that the boundary of region includes the conical singularities. In order to discuss such situation for the porous medium equation, Lian and Liu [33] studied the initial-boundary value problem of the porous medium equation

    ut=Δum+V(x)up, xD,t>0

    in a cone D=(0,)×Sn1, and they proved that if the nonlinear power index p belongs to a suitable interval then the problem has no global non-negative solutions for any non-negative initial datum u0 unless u00. Beside that, they also showed that this problem has global solutions for some u00 when the nonlinear power index p is out of that interval. Considering the following porous medium equation

    utΔum=f(u,t), xD,t>0 (1.5)

    Roidos and Schrohe [44] obtained some results about the existence, uniqueness and maximal Lp-regularity of a short time solution and showed the short time asymptotic behavior of the solution near the conical point. In addition, the behavior for large times of non-negative solutions to the linear Dirichlet problem of equation (1.5) in cone-like domains was obtained by Andreucci [2]. For more related work, we refer to Laptev [19,30] for the high-order evolution inequalities in cone-like domains and Mazzeo et al [38] for the Ricci flow on asymptotically conical surfaces.

    As a differential operator reflecting the diffusion form on the conical singular manifold, the emergence of the cone operator ΔB brings about the first problem that needs to be solved urgently is the existence of solution of the differential equation. Unlike the usual smooth domain, the appearance of cone singularities makes the classical embedding theorem fail, which baffles the proof of the existence theorem. Therefore, in order to overcome the difficulty and obtain the existence of the weak solution by variational method in H1,n22,0(B), Chen et al [12] considered the nonlinear Dirichlet boundary value problems on manifolds with conical singularities

    ΔBu=|u|p1u, xintB,

    and obtained the existence of non-trivial weak solution. Moreover, they also established the well-known cone Sobolev inequality and Poincarˊe inequality in the weighted Sobolev spaces. Subsequently, Chen et al [13] extended this result to the nonlinear elliptic equations with a general nonlinear source and the critical Sobolev exponents respectively. These works describe the mechanism by which cone differential operators act on the regularity of solutions of differential equations. Of course, differential operators do more than affect the regularity of solutions. The influence of degenerate differential operators on the solutions of nonlinear elliptic equations is also reflected in many other aspects such as the eigenvalue problems and the existence of multiple solutions. Far from being complete, we refer the readers to [22,42,40,41,3] and references therein. It is worth to mention here that the change of the variational structure of equations caused by the non-classical forms of differential operators can bring big challenges to the application of variational techniques [21,7]. Further more, when utilizing the variational techniques the geometrical feature needs to be taken into consideration [20]. In fact we need to overcome the difficulties mentioned above in the application of variational techniques to the parabolic version.

    In order to understand the effect of different initial data belonging to H1,n22,0(B) on the well-posedness of the solution, Chen and Liu [11] investigated the following conical degenerate parabolic equation on the conical manifold

    utΔBu=|u|p1u, xintB,t>0,

    and obtained not only the existence of global solutions with exponential decay, but also the blow up in finite time under low initial energy level and critical initial energy level. Recently, Mohsen and Morteza [1] studied the semilinear conical-degenerate parabolic equation

    tuΔBu+V(x)u=g(x)|u|p1u, xintB,t>0,

    where V(x)L(intB)C(intB) is the positive potential function and gL(intB)C(intB) is a non-negative weighted function. Then they got the results of global solutions with exponential decay and showed the finite time blow up of solutions on manifolds with conical singularities under subcritical initial energy level J(u0)<d.

    Our goal is to obtian local and global well-posedness of solutions to problem (1.1)-(1.3). In details, by the potential well method, we classify the initial datum and give a threshold condition, which tells us that as long as the initial datum falls into the specified invariant set and the initial energy satisfies J(u0)d, the solution exists as a global one or blow up in a finite time. Moreover, for global solutions, we give an estimation of the asymptotic behavior of them. For the solution that blows up in finite time, we estimate the lower bound of the blowup time. Different from previous proofs in the potential well framework, we integrate the proofs of the sub-critical and critical initial energy cases as a whole part, which makes the results much more concise. Last but not least, we also investigate the finite time blow up of solutions to problem (1.1)-(1.3) at high initial energy levels. By giving a sufficient condition tied to the initial data in weighted Sobolev space H1,n22,0(B), the theorem can not only explain what kind of initial data cause the solution to blow up in finite time, but also get the corresponding upper bound estimate of the blowup time.

    The content of this paper is arranged as follows. In Section 2, we give the geometric description of conical singularities, the definitions of the weighted Sobolev spaces and several propositions of the manifold with conical singularities. Then we introduce the potential well structure for problem (1.1)-(1.3) and prove a series of corresponding properties in Section 3. Section 4 is concerned with the local existence and uniqueness theory. In Section 5, we not only prove the invariant manifolds, global existence and decay of solutions to describe the corresponding asymptotic behavior, but also prove the finite time blow up of solutions and estimate the lower bound of blowup time in Theorem 5.2. In Section 6, we give a sufficient condition to obtain the finite time blow up of the solution in Theorem 6.4. In particular, we also estimate the upper bound of the blowup time of the solution. Finally, some remarks and acknowledgements about this paper are given.

    In this section, main definitions of the manifold with conical singularities together with a brief description of its properties are given, for more details we refer to [45,18] and the references therein. Furthermore, we introduce some functional inequalities on the manifold with conical singularity, for more applications of these inequalities one can refer to [12,13].

    For l2, let XSl be a bounded open set in a unit sphere of R1+lˆx, and define the straight cone XΔ by

    XΔ={ˆxR1+l | ˆx=0 or ˆx|ˆx|X}.

    The polar coordinates (ρ,θ) gives us a description of XΔ{0} in the form X=R+×X, which is called the open stretched cone with the base X, and {0}×X is the boundary of X.

    Now we extend it to a more general situation by describing the singular space associated with a manifold with conical singularities. A finite dimensional manifold B with finite conical singularities B0={b1,b2,...,bN} has the following two properties:

    a) BB0 is a C manifold.

    b) Any biB0 (i=1,2,...,N) has an open neighborhood G in B, such that there is a homeomorphism χ:GXΔ for some closed compact C manifold X=X(bi), and φ restricts a diffeomorphism φ:G{0}X.

    By above assumptions we can define the stretched manifold associated with B. Let B be a C manifold with compact C boundary BbiB0X=X(bi) for which there exists a diffeomorphism

    BB0BB:=intB.

    Furthermore, the restriction of this diffeomorphism to Gibi is also a diffeomorphism GibiUiX(bi), where GiB is an open neighborhood near bi, and UiB is a collar neiborhood with Ui[0,1)×X(bi).

    The typical differential operators on a manifold with conical singularities, i.e., the so-called Fuchsian type operators in a neighborhood of xb=0, have the following form

    A:=xμbμk=0ak(xb)(xbxb)k

    with (xb,˜x)X and ak(xb)C(¯R+,Diffμk(X))[45,18]. For such singular operators, we introduce the following cone weighted Sobolev space.

    Definition 2.1 (The space Hm,γp(Rn+)). For (xb,˜x)R+×Rl:=Rn+, mN, γR and 1<p<, assume u(xb,˜x)D(Rn+), where the dual (C0(Rn+))=D(Rn+) is the space of all distributions in Rn+. We denote the spaces

    Hm,γp(Rn+):={uD(Rn+) | xNpγb(xbxb)kα˜xuLp(Rn+)}

    for any kN, multi-index αNn with k+|α|m.

    Therefore, Hm,γp(Rn+) is a Banach space with the following norm

    uHm,γp(Rn+)=k+|α|m(R+RNxNb|xγb(xbxb)kα˜xu(xb,˜x)|pdxbxbd˜x)1p.

    Definition 2.2 ([45] The space Hm,γp,0(B)). We give the definition of Hm,γp,0(B) as follows

    (ⅰ) Let X be a closed compact C manifold covered by open neighborhoods O={O1,...,ON} of the coordinate. Let the subordinate partition of unity {ψ1,...,ψN} be fixed and charts χj:OjRn, j=1,...,N. Then we say that uHm,γp(X) if and only if uD(X), whose norm is defined as follows

    uHm,γp(X)=(Nj=1(1×χj)1ψjupHm,γp(Rn+))1p,

    where 1×χj:C0(R+×Rn)C0(R+×Oj) is the pull-back function with respect to 1×χj:R+×OjR+×Rn. And the closure of C0(X) in Hm,γp(X) is denoted as Hm,γp,0(X).

    (ⅱ) Let B be the stretched cone manifolds. Then Hm,γp(B) denotes the subspace of all uWm,ploc(intB) such that

    Hm,γp(B)={uWm,ploc(intB) | ωuHm,γp,0(X)}

    for any cut-off function ω(xb) supported by a collar neighborhood of xb(0,1). Moreover, the subspace Hm,γp,0(B) of Hm,γp(B) is defined as follows

    Hm,γp,0(B)=[ω]Hm,γp,0(X)+[1ω]Wm,p0(intB),

    where the classical Sobolev space Wm,p0(intB) denotes the closure of C0(intB) in Wm,p(˜B) for ˜B as a closed compact C manifold of dimension n containing B as a submanifold with boundary.

    Proposition 1 (Cone Sobolev inequality [12]). Assuming 1p<n,  1p=1p1n then for all u(xb,˜x)C0(Rn+) the following estimate

    uLγp(Rn+)c1(xbxb)uLγp(Rn+)+(c1+c2)ni=1xiuLγp(Rn+)+c2c3uLγp(Rn+) (2.1)

    holds, where γ,γR are constants with γ=γ1, and

    c1=(n1)pn(np),
    c2=(n1)p|(n1)(γ1)(n1)pnp|1nn(np)

    and

    c3=(n1)pnp.

    Moreover, for u(xb,˜x)H1,γp,0(Rn+), there holds

    uLγp(Rn+)CuH1,γp,0(Rn+),

    where C=c1+c2.

    Proposition 2 (Cone Poincarˊe inequality [12]). Let B=[0,1)×XRn+ is bounded and 1<p<, γ is a constant. If u(xb,˜x)H1,γp,0(B) then

    u(xb,˜x)Lγp(B)cBu(xb,˜x)Lγp(B),

    where the optimal constant c depends only on B and p.

    Proposition 3 (Cone Hölder inequality [13]). If uLnpp(B),vLnpp(B) with p,p(1,) and 1p+1p=1, then we have the following cone type Hölder inequality

    B|uv|dxbxbd˜x(B|u|pdxbxbd˜x)1p(B|v|pdxbxbd˜x)1p.

    Proposition 4 (Eigenvalue problem [13]). There exist 0<λ1<λ2λ3...λk..., such that for each k1, the following Dirichlet problem

    {ΔBψk=λkψk,  (xb,˜x)intB,ψk=0,  (xb,˜x)B,

    admits a non-trivial solution in H1,n22,0(B). Moreover, {ψk}k1 constitute an orthonormal basis of the Hilbert space H1,n22,0(B).

    In order to state our main results, we shall introduce some definitions and notations as follows. In the sequel, for convenience we denote

    (u,v)B=Buvdxbxbd˜x   and   uLnpp(B)=(B|u|pdxbxbd˜x)1p.

    Furthermore, we denote H1,n22,0(B):=((B,B)B+(,)B)12. Of course, the norm H1,n22,0(B) is equivalent to the norm Ln22(B), which can be stemmed from Proposition 2.

    Throughout the paper, C will be used to denote various positive constants, whose value may change from line to line, and its dependence on other variables will be emphasised only if needed.

    Definition 3.1. (Weak solution). A function u is called a weak solution to problem (1.1)-(1.3) on [0,T]×B, if it satisfies

    (ⅰ) uL(0,T;H1,n22,0(B)) and utL2(0,T;H1,n22,0(B));

    (ⅱ) u(0)=u0;

    (ⅲ) for any ηH1,n22,0(B), the identity

    But(t)ηdxbxbd˜x+BBut(t)Bηdxbxbd˜x+BBu(t)Bηdxbxbd˜x
    =B|u(t)|p1u(t)ηdxbxbd˜x

    holds for a.e. t[0,T].

    Let us introduce the following functionals on the cone Sobolev space H1,n22,0(B) as the potential energy functional

    J(u)=12B|Bu|2dxbxbd˜x1p+1B|u|p+1dxbxbd˜x (3.1)

    and the so-called Nehari functional

    I(u)=B|Bu|2dxbxbd˜xB|u|p+1dxbxbd˜x. (3.2)

    Then J(u) and I(u) are well-defined and belong to C1(H1,N22,0(B),R).

    The weak solution u(t) in Definition 3.1 satisfies the conservation of energy, i.e.,

    t0ut(τ)2H1,n22,0(B)dτ+J(u(t))=J(u0),   0t<T. (3.3)

    By making use of the functionals above, we define the potential well depth d as follows

    d=inf (3.4)

    where the Nehari manifold

    \begin{align} \mathcal{N} = \left\{u\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})\ \Big|\ I(u(t)) = 0,\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u(t)|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\neq0\right\} \end{align} (3.5)

    separates the whole space \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) into the following two unbounded manifolds

    \begin{align} \mathcal{N}_{+} = \left\{u\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})\ \Big|\ I(u(t)) > 0\right\} \end{align} (3.6)

    and

    \begin{align} \mathcal{N}_{-} = \left\{u\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})\ \Big|\ I(u(t)) < 0\right\}. \end{align} (3.7)

    Furthermore, we introduce the following potential well

    \begin{align} \mathcal{W}: = \mathcal{N}_{+} \cup \{0\} \end{align} (3.8)

    and the outside of the corresponding potential well

    \begin{align} \mathcal{V}: = \mathcal{N}_{-}. \end{align} (3.9)

    Now, we give some corresponding properties of the potential well as follows.

    Lemma 3.2 (The properties of the energy functional J(u) ). Assume that u\in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) and \|\nabla_{\mathbb{B}}u\|_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\neq0, we have

    (i) \lim_{\lambda\rightarrow0}j(\lambda) = 0 and \lim_{\lambda\rightarrow +\infty}j(\lambda) = -\infty, where j(\lambda): = J(\lambda u);

    (ii) there exist a unique \lambda^{*} = \lambda^{*}(u)>0, such that j'(\lambda^{*}) = 0;

    (iii) j(\lambda) is strictly increasing on 0\leq\lambda<\lambda^{*}, strictly decreasing on \lambda>\lambda^{*} and takes the maximum at \lambda = \lambda^{*};

    (iv) i(\lambda)>0 for 0\leq\lambda<\lambda^{*}, i(\lambda)<0 for \lambda>\lambda^{*} and i(\lambda^{*}) = 0, where i(\lambda): = I(\lambda u).

    Proof. (ⅰ) From the definition of J(u), we know that

    \begin{align*} j(\lambda) = J(\lambda u) = \frac{\lambda^{2}}{2}\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}-\frac{\lambda^{p+1}}{p+1} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}, \end{align*}

    which gives \lim_{\lambda\rightarrow0}j(\lambda) = 0 and \lim_{\lambda\rightarrow +\infty}j(\lambda) = -\infty.

    (ⅱ) An easy calculation shows that

    \begin{align} j'(\lambda) = \lambda\left(\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}-\lambda^{p-1} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\right). \end{align} (3.10)

    Then taking j'(\lambda) = 0 we obtain that

    \begin{align*} \lambda^{*} = \left(\frac{\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}}{\int_{\mathbb{B}}|u|^{p+1} \frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}}\right)^{\frac{1}{p-1}} > 0. \end{align*}

    (ⅲ) By a direct calculation, (3.10) gives j'(\lambda)>0, for 0<\lambda<\lambda^{*},\ j'(\lambda)<0 for \lambda^{*}<\lambda<\infty. Hence, the conclusion of (ⅲ) holds.

    (ⅳ) The conclusion follows from

    \begin{align*} i(\lambda) = I(\lambda u) = \lambda^{2}\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}-\lambda^{p+1} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} = \lambda j'(\lambda). \end{align*}

    Next we give the relationship between \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} and I(u) in the following lemma.

    Lemma 3.3 (The properties of the Nehari functional I(u) ). Suppose that u\in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) and r = \left(\frac{1}{(c^2+1)C^{p+1}_{*}}\right)^{\frac{1}{p-1}} , where c and C_{*} are the optimal constants of the cone Poincar \acute{e} inequality and cone Sobolev inequality respectively.

    (i) If 0<\|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}<r, then I(u)>0.

    (ii) If I(u)<0, then \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}>r.

    (iii) If I(u) = 0, then \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\geq r or \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = 0.

    Proof. (ⅰ) From 0<\|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}<r, we have

    \begin{align*} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}&\leq C^{p+1}_{*}\|u\|^{p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} < C^{p+1}_{*}r^{p-1}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\& = \frac{1}{c^2+1}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} \le \|\nabla_{\mathbb{B}}u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}. \end{align*}

    Then by the definitions of r and I(u) we obtain that I(u)>0.

    (ⅱ) It is easy to see that \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\neq 0 by I(u)<0 . Combining cone Poincar \acute{\rm e} inequality and cone Sobolev inequality, it follows that

    \begin{align*} \|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} & = \int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}+\int_{\mathbb{B}}|u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ &\leq (c^2+1)\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ & < (c^2+1)\int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \notag\\ &\leq (c^2+1)C^{p+1}_{*}\|u\|^{p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, \notag \end{align*}

    then we get

    \begin{align*} \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > \left(\frac{1}{ (c^2+1)C^{p+1}_{*}}\right)^{\frac{1}{p-1}} = r. \end{align*}

    (ⅲ) If I(u) = 0 and \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\neq 0 , then by

    \begin{align*} \|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}& = \int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}+\int_{\mathbb{B}}|u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ &\leq (c^2+1)\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ & = (c^2+1)\int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \notag\\ &\leq (c^2+1)C^{p+1}_{*}\|u\|^{p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, \notag \end{align*}

    we get \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\geq \left(\frac{1}{ (c^2+1)C^{p+1}_{*}}\right)^{\frac{1}{p-1}} = r.

    The depth d can be estimated as follows.

    Lemma 3.4 (The potential well depth). Suppose that 1<p<\frac{n+2}{n-2} and d is defined as (3.4). Then we have

    \begin{align} d = \inf\limits_{u\in\mathcal{N}}J(u) = \frac{p-1}{2(p+1)}\left(\frac{1}{ (c^2+1)^{\frac{p+1}{2}}C^{p+1}_{*}}\right)^{\frac{2}{p-1}}, \end{align} (3.11)

    where c is the best coefficient of the cone Poincar \acute{ e} inequality.

    Proof. Suppose that u\in \mathcal{N}, then \|u\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\ge r from Lemma 3.3. Thus according to the definition of J(u) , I(u) and cone Poincar \acute{\rm e} inequality, we arrive at

    \begin{align*} J(u)& = \left(\frac{1}{2}-\frac{1}{p+1}\right)\int_{\mathbb{B}}|\nabla_{\mathbb{B}}u|^{2}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}+\frac{1}{p+1}I(u) \\ &\geq\left(\frac{1}{2}-\frac{1}{p+1}\right)\frac{1}{c^2+1}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} \notag\\ &\geq \frac{p-1}{2(p+1)(c^2+1)}r^{2},\notag \end{align*}

    which gives (3.11).

    Lemma 3.5 (Osgood Lemma, [10]). Let \rho:[t_0,T]\rightarrow [0,\alpha] be a measurable function, \gamma is a locally integrable, positive function defined on [t_0,T] , \mu:[0,\alpha]\rightarrow [0,+\infty) is a nondecreasing, continuous function and \mu(0) = 0 , a\ge0 is a constant. If

    \rho(t)\le a+\int^{t}_{t_0}\gamma(s)\mu(\rho(s)){\rm d}s,

    holds almost everywhere for t\in[t_0,T] , then

    -M(\rho(t))+M(a)\le \int^{t}_{t_0}\gamma(s){\rm d}s,

    is true almost everywhere for t\in[t_0,T] when a>0 , where M(x) = \int^{\alpha}_{x}\frac{{\rm d}s}{\mu(s)} . In addition, when a = 0 and M(0) = \infty , \rho(t) = 0 holds almost everywhere for t\in[t_0,T] .

    In this section, we prove the local existence of the solution of problem (1.1)-(1.3). The local existence theorem is given as follows,

    Theorem 4.1 (Local existence). Suppose that u_{0}\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) . Then there exist T>0 and a unique weak solution u\in C([0,T],\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) of problem (1.1)-(1.3) on [0,T]\times\mathbb{B} . Moreover, if

    T_{\max} = \sup\{T > 0: u = u(t)\ \mathit{\text{exists on}}\ [0,T]\} < \infty,

    then

    \lim\limits_{t\rightarrow T_{\max}}\| u(t)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = \infty.

    Proof. We divide the proof into 2 steps.

    Step 1. Local exisence. We prove the local existence of the solution to problem (1.1)-(1.3) by virtue of the Galerkin method and the compactness property [35]. For the initial data u_0\in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) , let R^2: = 2\|\nabla_\mathbb{B}u_0\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} . For every m\ge1 , let \Psi_{m} = \text{Span}\{\psi_{1}, \psi_{2},...,\psi_{m}\} , where \{\psi_{j}\} is the orthonormal complete system of eigenfuctions of -\Delta_{\mathbb{B}} in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) such that \|\psi_{j}\|_{L^{\frac{n}{2}}_{2}(\mathbb{B})} = 1 for all j . Then, \{\psi_{j}\} is orthonormal and complete in L^{\frac{n}{2}}_{2}(\mathbb{B}) and \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) by Proposition 4, and we denote by \{\lambda_{j}\} the corresponding eigenvalues. Let

    \begin{align} u_{m0} = \sum\limits^{m}_{j = 1}(\nabla_{\mathbb{B}}u_0,\nabla_{\mathbb{B}}\psi_{j})_{\mathbb{B}}\psi_{j},\ m\in \mathbb{N}^+ \end{align} (4.1)

    such that u_{m0}\in\Psi_{m} and

    \begin{align} u_{m0}\rightarrow u_0\ \text{in}\ \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})\ \text{as}\ m\rightarrow \infty. \end{align} (4.2)

    For all m\ge1 , we seek C^1 -continuous functions g_{m1}(t),...,g_{mm}(t) to form an approximate solution to problem (1.1)-(1.3) of the form

    \begin{align} u_{m}(t, x_b, \widetilde{x}) = \sum^{m}_{j = 1}g_{mj}(t)\psi_{j}(x_b, \widetilde{x}), \end{align} (4.3)

    which solves the problem

    \begin{align} \begin{cases} (u_{mt}(t),\psi_{j})_{\mathbb{B}} +(\nabla_{\mathbb{B}}u_{mt}(t),\nabla_{\mathbb{B}}\psi_{j})_{\mathbb{B}} +(\nabla_{\mathbb{B}}u_m(t),\nabla_{\mathbb{B}}\psi_{j})_{\mathbb{B}} = (|u_m(t)|^{p-1}u_m(t),\psi_{j})_{\mathbb{B}},\\ u_{m}(0,x_b, \widetilde{x}) = u_{m0}, \end{cases} \end{align} (4.4)

    for j = 1,2,...,m and t\ge0 . Problem (4.4) is equivalent to the following systems of ODEs

    \begin{align} \begin{cases} (1+\lambda_j)g_{mjt}(t)+\lambda_jg_{mj}(t) = (|u_m(t)|^{p-1}u_m(t),\psi_j)_{\mathbb{B}},\\ g_{mj}(0) = (\nabla_{\mathbb{B}}u_0,\nabla_{\mathbb{B}}\psi_j)_{\mathbb{B}},\ j = 1,2,...,m. \end{cases} \end{align} (4.5)

    One can deduce that for any fixed m there exists a t_m>0 and a unique solution g_{mj}\in C^1[0,t_m] of the Cauchy problem (4.5) by the Cauchy-Peano theorem since f(t): = (|u_m(t)|^{p-1}u_m(t),\psi_j)_{\mathbb{B}} is continuous with respect to t . Multiplying (4.4) by g_{mjt}(t) and summing for j , we have

    \begin{align*} &(u_{mt}(t), u_{mt}(t))_{\mathbb{B}}+(\nabla_{\mathbb{B}}u_{mt}(t),\nabla_{\mathbb{B}}u_{mt}(t))_{\mathbb{B}}+(\nabla_{\mathbb{B}}u_m(t),\nabla_{\mathbb{B}}u_{mt}(t))_{\mathbb{B}} \\& = (|u_m(t)|^{p-1}u_m(t), u_{mt}(t))_{\mathbb{B}}, \end{align*}

    which tells us that for all t\in[0,t_m] ,

    \begin{align} &\|u_{mt}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\|\nabla_{\mathbb{B}}u_{mt}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\frac{1}{2}\frac{{\rm d}}{{\rm d}t}\|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ = &\int_{\mathbb{B}}|u_m(t)|^{p-1}u_m(t)u_{mt}(t)\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}. \end{align} (4.6)

    For the last term in equation (4.6), by using the cone Hölder inequality (Proposition 1), cone Sobolev inequality (Proposition 3) and Young's inequality, we deduce

    \begin{align} \int_{\mathbb{B}}|u_m(t)|^{p-1}u_m(t)u_{mt}(t)\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \le&\int_{\mathbb{B}}|u_m(t)|^{p}|u_{mt}(t)|\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ \le&\|u_{mt}(t)\|_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\|u_m(t)\|^p_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ \le&C^{p+1}_*\|\nabla_\mathbb{B}u_{mt}(t)\|_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\|\nabla_\mathbb{B}u_m(t)\|^p_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ \le&\frac12\|\nabla_\mathbb{B}u_{mt}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\frac{C^{2(p+1)}_*}{2}\|\nabla_\mathbb{B}u_m(t)\|^{2p}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}, \end{align} (4.7)

    which together with (4.6) gives,

    \begin{align} \|u_{mt}(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\frac{{\rm d}}{{\rm d}t}\|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\le C^{2(p+1)}_*\|\nabla_\mathbb{B}u_m(t)\|^{2p}_{L^{\frac{n}{2}}_{2}(\mathbb{B})},t\in[0,t_m]. \end{align} (4.8)

    Integrating (4.8) over [0,t] , we obtain

    \begin{align*} &\|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\int^{t}_{0}\|u_{mt}(s)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}s \\&\le \|\nabla_{\mathbb{B}}u_{m}(0)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} +C_*^{2(p+1)}\int^{t}_{0}\|\nabla_\mathbb{B}u_m(s)\|^{2p}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}{\rm d}s, \end{align*}

    which combining the formula (4.2) and the fact R^2: = 2\|\nabla_\mathbb{B}u_0\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} shows that there exists some m_0\in \mathbb N^+ such that

    \begin{align} \|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\int^{t}_{0}\|u_{mt}(s)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}s \le R^2 +C_*^{2(p+1)}\int^{t}_{0}\|\nabla_\mathbb{B}u_m(s)\|^{2p}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}{\rm d}s, \end{align} (4.9)

    for t\in[0,t_m] and m\ge m_0 . In order to estimate the first term in (4.9), one can apply the Osgood lemma [10,Lemma 5.2.1] as follows.

    For any fixed m\ge m_0 , let \rho(t): = \|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} , t\in[0,t_m] , since g_{mj}\in C^1[0,t_m] , we deduce that \max\limits_{t\in[0,t_m]}|g_{mj}(t)| is bounded, j = 1,...,m , then

    \rho(t) = \sum\limits^m_{j = 1}\lambda_j|g_{mj}(t)|^2\le m\max\limits_{1\le j\le m}\left(\lambda_j\max\limits_{t\in[0,t_m]} |g_{mj}(t)|^2\right): = \alpha_m < +\infty,

    which implies that \rho:[0,t_m]\rightarrow [0,\alpha_m] is a measurable function. Now we pick

    \begin{align*} \gamma(s)&\equiv C^{2(p+1)}_*:[0,t_m]\rightarrow \mathbb{R}\ \text{is a locally integrable, positive function},\\ \mu(s)& = s^p:[0,\alpha_{m}]\rightarrow [0,+\infty)\ \text{is a continuous, non-decreasing function}, \end{align*}

    which satisfies \mu(0) = 0 , and a = R^2 , then,

    \begin{align*} M(\nu) = \int^{\alpha_m}_\nu\frac{{\rm d}s}{\mu(s)} = \frac{1}{p-1}\left(\nu^{-(p-1)}-\alpha_m^{-(p-1)}\right). \end{align*}

    Therefore, we have

    \begin{align*} -M(\rho(t))+M(a)\le \int^t_0\gamma(s){\rm d}s = C^{2(p+1)}_* t, \end{align*}

    that is,

    \begin{align*} \frac{1}{p-1}\left(\alpha_m^{-(p-1)}-\|\nabla_{\mathbb{B}}u_{m}(t)\|^{-2(p-1)}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\right) +\frac{1}{p-1}\left(R^{-2(p-1)}-\alpha_m^{-(p-1)}\right)\le C^{2(p+1)}_* t. \end{align*}

    By a simple calculation, we obtain

    \begin{align} \|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\leq \left[R^{2(1-p)}-(p-1)C_*^{2(p+1)} t \right]^{-\frac{1}{p-1}}, \ \ t\in[0,t_m]. \end{align} (4.10)

    Taking

    \begin{align} T = T(R): = \frac{R^{2(1-p)}}{2 (p-1)C_*^{2(p+1)}}, \end{align} (4.11)

    it follows from estimates (4.9)-(4.11) that

    \begin{equation} \|\nabla_{\mathbb{B}}u_{m}(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\int^{t}_{0}\|u_{mt}(s)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}s\le C(R),\ \ t\in[0, T],\ m\ge m_0, \end{equation} (4.12)

    where

    \begin{equation*} C(R) = R^2\left(1+\frac{2^{\frac{1}{p-1}}}{p-1}\right). \end{equation*}

    Then combining with (4.12) and Proposition 1 we obtain

    \begin{equation} \||u_m(t)|^{p-1}u_m(t)\|_{L^{\frac{np}{p+1}}_{\frac{p+1}{p}}(\mathbb{B})} = \|u_m(t)\|^p_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\leq C_*^{p}\|\nabla_{\mathbb{B}}u_{m}(t)\|^p_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\leq \left(C_*\sqrt{C(R)}\right)^{p} \end{equation} (4.13)

    for all t\in [0, T] and m\geq m_0 , which means

    \begin{align*} &\{u_m\}\ \ \hbox{is bounded in}\ \ L^\infty(0, T; \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})),\\ &\{u_{mt}\}\ \ \hbox{is bounded in}\ \ L^2(0, T; \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})),\\ &\{|u_m|^{p-1}u_m\}\ \ \hbox{is bounded in}\ \ L^\infty(0, T; L^{\frac{np}{p+1}}_{\frac{p+1}{p}}(\mathbb{B})). \end{align*}

    Hence, by the Aubin-Lions-Simon Lemma [46] and the weak compactness there exist a u\in C([0, T]; \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) and a subsequence of \{u_{m}\} , which is still denoted by \{u_{m}\} , such that

    \begin{equation} \begin{split} &u_{m}\rightarrow u\ \hbox{in}\ L^{\infty}(0,T;\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}))\ \hbox{weakly star}, \\ &u_{m}\rightarrow u\ \hbox{in}\ C([0, T]; L^{\frac{n}{2}}_{2}(\mathbb{B}))\ \hbox{and a.e. in}\ [0, T]\times \text{int}\mathbb{B},\\ & u_{mt}\rightarrow u_t\ \hbox{in}\ L^2(0, T;\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}))\ \hbox{weakly},\\ & |u_{m}|^{p-1}u_m\rightarrow \chi \ \hbox{in}\ L^\infty(0,T; L^{\frac{np}{p+1}}_{\frac{p+1}{p}}(\mathbb{B}))\ \hbox{weakly star}, \end{split} \end{equation} (4.14)

    which together with the Lions Lemma [35,Chap. 1,p12] deduce that

    \begin{equation*} \chi = |u|^{p-1}u. \end{equation*}

    Then for each j fixed, let m\rightarrow \infty in (4.4), we have

    \begin{align*} &\int_{\mathbb{B}}u_{t}(t) \psi_{j}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} + \int_{\mathbb{B}}\nabla_{\mathbb{B}}u(t)\nabla_{\mathbb{B}} \psi_{j}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} +\int_{\mathbb{B}}\nabla_{\mathbb{B}}u_{t}(t)\nabla_{\mathbb{B}}\psi_{j}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \\& = \int_{\mathbb{B}}|u(t)|^{p-1}u(t)\psi_{j}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \end{align*}

    for a.e. t\in [0, T] and every j = 1, 2, \cdots . By the fact that \{\psi_{j}\} is the complete orthonormal basis in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) , we have

    \begin{align*} (u_{t}(t), \eta)_{\mathbb{B}} +(\nabla_{\mathbb{B}}u(t), \nabla_{\mathbb{B}}\eta)_{\mathbb{B}} +(\nabla_{\mathbb{B}}u_{t}(t), \nabla_{\mathbb{B}}\eta)_{\mathbb{B}} = \int_{\mathbb{B}}|u(t)|^{p-1}u(t) \eta\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \end{align*}

    for a.e. t\in [0, T] and every \eta\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) . By the convergence (4.14) and formula (4.12), we have

    \begin{equation} \|\nabla_{\mathbb{B}}u(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\int^{t}_{0}\|u_{t}(s)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}s\le C(R),\ \ t\in[0, T]. \end{equation} (4.15)

    Moreover, it follows from the fact u_m\rightarrow u in C([0, T]; L^{\frac{n}{2}}_{2}(\mathbb{B})) that

    \begin{equation*} u_{m0} = u_m(0)\rightarrow u(0)\ \ \hbox{in}\ \ L^{\frac{n}{2}}_{2}(\mathbb{B}), \end{equation*}

    which combining with (4.2) implies that

    \begin{equation*} u(0) = u_0\ \ \hbox{in}\ \ \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}). \end{equation*}

    Thus, u\in C([0, T]; \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) is a weak solution of problem (1.1)-(1.3) on [0,T]\times\mathbb{B} .

    Step 2. Uniqueness. Suppose that there are two solutions u and v to problem (1.1)-(1.3) corresponding to the initial data u_0\in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) . Then, w = v-v solves the following problem

    \begin{align} &w_{t}-\Delta_{\mathbb{B}} w_{t}-\Delta_{\mathbb{B}}w = |u|^{p-1}u-|v|^{p-1}v, && (x_b, \widetilde{x})\in \text{int}\mathbb{B},t > 0, \end{align} (4.16)
    \begin{align} &w(x_b, \widetilde{x},0) = 0,&& (x_b, \widetilde{x})\in \text{int}\mathbb{B}, \end{align} (4.17)
    \begin{align} &w(0, \widetilde{x},t) = 0, && (0, \widetilde{x})\in\partial \mathbb{B},t\ge0. \end{align} (4.18)

    Since w\in C([0, T]; \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) , one can multiply both side of equation (4.16) by w(t) and show that

    \begin{equation} \begin{split} & \frac{\mathrm d}{\mathrm d t}\left(\|w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\|\nabla_{\mathbb{B}}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\right)+ 2\|\nabla_{\mathbb{B}}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ = \ &2\int_{\mathbb{B}}\left(|u|^{p-1}u-|v|^{p-1}v\right)w \frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}, \ \ t\in [0, T], \end{split} \end{equation} (4.19)

    For the last term in equation (4.19), by using the cone Hölder inequality (Proposition 1), cone Sobolev inequality (Proposition 3) and estimate (4.15), we deduce

    \begin{equation*} \begin{split} &2\int_{\mathbb{B}}\left(|u|^{p-1}u-|v|^{p-1}v\right)w \frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ = \ & 2p\int_{\mathbb{B}}\int_0^1|\theta u+(1-\theta)v|^{p-1}|w|^2 {\rm d} \theta \frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} \\ \leq\ & 2p\int_{\mathbb{B}}\left(|u|^{p-1}+|v|^{p-1}\right)|w|^2\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\\ \leq\ & 2p \left( \|u\|^{p-1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}+ \|v\|^{p-1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\right) \|w\|^2_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ \leq\ & 2pC_*^{p+1} \left( \|\nabla_\mathbb{B}u\|^{p-1}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+ \|\nabla_\mathbb{B}v\|^{p-1}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\right) \|\nabla_\mathbb{B}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ \leq\ & C_1(R)\|\nabla_\mathbb{B}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}, \ \ t\in [0, T], \end{split} \end{equation*}

    with C_1(R) = 4p C_*^{p+1}(C(R))^{\frac{p-1}{2}} , which together with (4.19) gives,

    \begin{equation} \frac{\mathrm d}{\mathrm d t}\left(\|w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\|\nabla_{\mathbb{B}}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\right)\leq C_1(R)\left(\|w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} +\|\nabla_{\mathbb{B}}w\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\right),\ \ t\in[0,T]. \end{equation} (4.20)

    Applying the Gronwall inequality to (4.20) and making use of (4.17), we have

    \begin{equation*} \|w(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\leq e^{C_1(R) t} \|w(0)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = 0,\ \ \forall t\in [0, T], \end{equation*}

    which leads to the uniqueness of weak solution.

    Concerning formula (4.11) we observe that the local existence time T merely depends on the norms of the initial data. Therefore, by using the similar idea as shown in [4], the solution can be continued as long as \|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} remains bounded. In details, taking u(\cdot,T) as the initial data and repeating above argument, we know that the problem (1.1)-(1.3) has a unique weak solution on the interval [0,T_1] ( T_1>T ). After an iteration process, we get a single increasing sequence \{T_k\}^{\infty}_{k = 1} such that the problem (1.1)-(1.3) has a unique weak solution on [0,T_k] , where T_k has two possibilities for sequence \{T_k\}^{\infty}_{k = 1} , that either T_{\max} = \lim\limits_{k\rightarrow \infty}T_k is finite, or T_{\max} = \infty . Moreover, if T_{\max} = \infty , then the problem (1.1)-(1.3) possesses a unique global solution. If T_{\max}<+\infty , then

    \begin{align} \limsup\limits_{t\rightarrow T_{\max}}\|\nabla_{\mathbb{B}} u(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} = +\infty. \end{align} (4.21)

    In fact, if

    \sup\limits_{t\in [0,T_{\max})}\|\nabla_{\mathbb{B}} u(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\le M(T_{\max}),

    then for any t_0\in(0,T_{\max}) , by taking u(\cdot,t_0) as the initial data, we know from the above argument that there exists a T^*>0 dependent on M and independent on t_0 such that the problem (1.1)-(1.3) has a unique weak solution on the interval [t_0,t_0+T^*] . Due to the arbitrariness of t_0\in(0,T_{\max}) , the weak solution of problem (1.1)-(1.3) can be extended to the interval [0,T_{\max}+\varepsilon] with arbitrary positive \varepsilon>0 , which contradicts the maximal interval [0,T_{\max}) . This contradiction shows that (4.21) holds.

    Moreover, we give the following corollary.

    Corollary 1 (Blow-up of the weak solution). Suppose that u is a weak solution of problem (1.1)-(1.3), which can be ensured by Theorem 4.1, if T_{\max}<\infty , then we have

    \begin{align*} \lim\limits_{t\rightarrow T_{\max}}\|u(t)\|^2_{L^{\frac{n}{s}}_{s}(\mathbb{B})} = \infty\ \ \ \mathit{\text{for}}\ \ \ s\geq \max \left\{1,\frac{n(p-1)}{2}\right\}. \end{align*}

    Proof. Recalling the definition of energy functional J(u) in (3.1) and the energy relation (3.3), we get

    \begin{align*} \frac{1}{2}\|\nabla_\mathbb{B}u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} \leq \frac{1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}+J(u_0)\ \text{for all}\ t \in [0,T_{\max}). \end{align*}

    Combining with the following Gagliardo-Nirenberg interpolation inequality

    \begin{align*} \|u\|_{L^{\frac{n}{B_1}}_{B_1}(\mathbb{B})} \leq C\|\nabla_\mathbb{B}u\|^{a}_{L^{\frac{n}{B_2}}_{B_2}(\mathbb{B})} \|u\|^{1-a}_{L^{\frac{n}{B_3}}_{B_3}(\mathbb{B})} \end{align*}

    for \frac{1}{B_1} = \left(\frac{1}{B_2}-\frac{1}{n}\right)a+\frac{1-a}{B_3}\ \text{and}\ 0<a<1 , we have

    \begin{align*} \frac{1}{2}\|\nabla_\mathbb{B}u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}-J(u_0) \leq \frac{1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})} \leq C\|\nabla_\mathbb{B}u\|^{(p+1)a}_{L^{\frac{n}{2}}_{2}(\mathbb{B})} \|u\|^{(p+1)(1-a)}_{L^{\frac{n}{s}}_{s}(\mathbb{B})}, \end{align*}

    where B_1 = p+1 , B_2 = 2 and B_3 = s with \frac{n(p-1)}{2}<s<p+1 such that 0<a = \frac{2n(p-1-s)}{(p+1)(sn+2s-ns)}<\frac{2}{p+1} . Thus we have

    \begin{align*} \frac{1}{2}\|\nabla_\mathbb{B}u\|^{2-(p+1)a}_{L^{\frac{n}{2}}_{2}(\mathbb{B})} -J(u_0)\|\nabla_\mathbb{B}u\|^{-(p+1)a}_{L^{\frac{n}{2}}_{2}(\mathbb{B})} \leq C\|u\|^{(p+1)(1-a)}_{L^{\frac{n}{s}}_{s}(\mathbb{B})}, \end{align*}

    which implies that \lim\limits_{t\rightarrow T_{\max}}\|u(t)\|^2_{L^{\frac{n}{s}}_{s}(\mathbb{B})} = \infty.

    In this section, we study the well-posedness of solutions of problem (1.1)-(1.3) in the case of sub-critical and critical initial energy levels. The succeeding result is given to show the invariant sets of the solution for problem (1.1)-(1.3).

    Lemma 5.1. Suppose that u_{0}\in \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) and J(u_0)<d , then

    (i) all weak solutions of problem (1.1)-(1.3) belong to \mathcal{W} provided that u_0\in \mathcal{W};

    (ii) all weak solutions of problem (1.1)-(1.3) belong to \mathcal{V} provided that u_0\in \mathcal{V}.

    Proof. (ⅰ) Suppose that J(u_0)<d , I(u_0)>0 and u(t) is the corresponding solution of problem (1.1)-(1.3). If u_0 = 0 , then u(t) = 0 , i.e., u(t)\in \mathcal{W} . If u_0\in \mathcal{N}_+ , i.e., I(u_{0})>0 , we claim that u(t)\in \mathcal{N}_+ for 0\le t<T . By reduction to absurdity, provided that there exists a t_{0}\in (0,T) such that u(t_{0})\in \mathcal{N} for the first time and u(t)\in \mathcal{N}_+ for 0\le t<t_0 , i.e., I(u(t))>0 for t\in [0,t_0) and I(u_0) = 0 . Since u(t_{0}) is a solution of problem (1.1)-(1.3), by the conservation of energy, it follows that

    \begin{align} \int^{t_0}_{0}\|u_t(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+J(u(t_0)) = J(u_{0}) < d \end{align} (5.1)

    for any t_{0}\in (0,T) . However, by the definition of d and I(u(t_{0})) = 0 , it implies that J(u(t_{0}))\ge d , which contradicts (5.1). Thus, u(t)\in \mathcal{N}_+\subset\mathcal{W} for any t\in [0,T) .

    (ⅱ) Similar to the proof of (i), we can obtain that u(t)\in \mathcal{V} for t\in [0,T) provided u_0\in \mathcal{V} .

    Theorem 5.2 (Global existence and asymptotic behavior). Suppose that u_{0}\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) and J(u_0)\le d . If u_0\in \mathcal{W}, then problem (1.1)-(1.3) admits a unique global weak solution u\in L^{\infty}(0,\infty;\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) with u_{t}\in L^{2}(0,\infty;\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) and satisfies

    \begin{equation} \|u(t)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\int^{t}_{0}\|u_t(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau < \left(1+\frac{2(p+1)(c^2+1)}{p-1}\right)d,\ \ t\ge0. \end{equation} (5.2)

    Moreover, the solution satisfies the estimate

    \begin{align*} \|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\leq e^{-2\beta t}\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})},\ \ \ 0\leq t < \infty, \end{align*}

    where \beta = 1-C^{p+1}_{*}\left(\frac{2(p+1)(c^2+1)}{p-1}J(u_0)\right)^{\frac{p-1}{2}}.

    Proof. We divide the proof into two parts, which are the global existence and the asymptotic behavior.

    Part Ⅰ: Global existence.

    First, we give the global existence of the solution of problem (1.1)-(1.3) for J(u_0)<d .

    Since the conclusion is trivial when u_0 = 0 , we only consider the case u_{0}\in \mathcal{W}\backslash\{0\} . From Theorem 4.1, let u be the weak solution of problem (1.1)-(1.3) corresponding to the initial data u_0 . It follows from (3.1) and (3.2) that

    \begin{equation} J(u(t)) = \frac{p-1}{2(p+1)}\|\nabla_{\mathbb{B}}u(t)\|^{2}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\frac{1}{p+1}I(u(t)) \end{equation} (5.3)

    and

    \begin{align} \int^{t}_{0}\|u_t(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+J(u(t)) = J(u_{0}) < d. \end{align} (5.4)

    for all t\in [0, T_{\max}) . Then by Lemma 5.1 we deduce that

    \begin{align*} u(t)\in\mathcal{W},\ \text{for all}\ 0 < t < T_{\max}, \end{align*}

    which implies

    \begin{equation} I(u(t)) > 0\ \ \hbox{for all}\ \ 0 < t < T_{\max}. \end{equation} (5.5)

    Thus, the combination of (5.3)-(5.5) and Proposition 1 shows that

    \begin{equation} \|u(t)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\int^{t}_{0}\|u_t(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau < \left(1+\frac{2(p+1)(c^2+1)}{p-1}\right)d,\ \ 0 < t < T_{\max}. \end{equation} (5.6)

    Therefore, by virtue of the Continuation Principle, it follows T_{\max} = \infty , i.e., u \in L^{\infty}(0,\infty;\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})) with u_{t}\in L^{2}(0,\infty; \mathcal{H}^{\frac{n}{2}}_{2}(\mathbb{B})) is a global weak solution of problem (1.1)-(1.3) corresponding to initial data u_0 and satisfies the estimate (5.2).

    For the case J(u_0) = d , we apply the idea of scalar transformation to get the corresponding global existence result.

    Let \gamma_k = 1-\frac{1}{k} , k\in \mathbb{N}^+ . Then there exists a sequence \{u_{k}(0)\} such that u_{k}(0) = \gamma_{k}u_0 . Consider the corresponding initial boundary value problems

    \begin{align} &u_{t}-\Delta_{\mathbb{B}} u_{t}-\Delta_{\mathbb{B}}u = |u|^{p-1}u,&& (x_b, \widetilde{x})\in \text{int}\mathbb{B},t > 0, \end{align} (5.7)
    \begin{align} &u(x_b, \widetilde{x},0) = u_{k}(0),&& (x_b, \widetilde{x})\in \text{int}\mathbb{B}, \end{align} (5.8)
    \begin{align} &u(0, \widetilde{x},t) = 0,&& (0, \widetilde{x})\in\partial \mathbb{B},t\ge0, \end{align} (5.9)

    according to Theorem 4.1 and the estimate (5.6), we know that the weak solution u to problem (5.7)-(5.9) corresponding to the initial data u_k(0) satisfies

    \begin{align} \int^{t}_{0}\|u_{t}(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+J(u(t)) = J(u_k(0)), 0 < t < T_{\max}, \end{align} (5.10)

    where T_{\max} is the maximal existence time of the solution u .

    By I(u_0)>0 , it follows from (ⅲ)-(ⅳ) in Lemma 3.2 that \lambda = \lambda(u_0)\in(0,1) and \lambda^*>1 , which implies that I(u_0)>I(\gamma_ku_0)>0 and J(\gamma_ku_0)<J(u_0)\le d . By the continuity of I(u) , J(u) with respect to u , we can choose a sufficiently large k such that \|u_{k}(0)-\gamma_{k}u_0\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}<\frac{1}{k} with I(u_{k}(0))>0 and J(u_{k}(0))<d , which means u_{k}(t)\in \mathcal{W} for [0,T_{k}) by Lemma 5.1.

    Since \lim\limits_{k\rightarrow \infty}\gamma_{k} = 1 , then for u_{k}(0)\in \mathcal{W} we have

    \begin{align} u_{k}(0)\rightarrow u_0\ \text{strongly in}\ \mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}),\ \text{as}\ k\rightarrow \infty, \end{align} (5.11)

    and

    \begin{align} \int^{t}_{0}\|u_{kt}(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+J(u_{k}(t)) = J(u_{k}(0)) < J(u_0) = d,\ \ 0 < t < \tilde{T}, \end{align} (5.12)

    then we can obtain the boundness of u_{k} by the same way in the previous case J(u_0)<d , namely (5.6), and the global existence result follows.

    Part Ⅱ: Asymptotic behavior.

    In this part, we start with the claim that u(t)\in\mathcal{W} for t\in[0,\infty) when J(u)\le d .

    If J(u_0)<d , then by the assumption that I(u_0)>0 we can derive immediately that u(t)\in\mathcal{W} for t>0 from Lemma 5.1.

    If J(u_0) = d , then we conclude that there exists an enough small t_0 such that I(u(t))>0 for t\in[0,t_0) through the continuity of I(u(t)) with respect to t and I(u_0)>0 . Next, we assert that I(u(t))>0 for t\in[t_0,\infty) . Arguing by contradiction, assuming that t_1\in[t_0,\infty) is the first time such that I(u(t_1)) = 0 , then by the definition of the potential well depth, we have

    \begin{align*} J(u(t_1))\ge d. \end{align*}

    In addition, by the energy identity (3.3), we have

    \begin{align*} 0 < J(u(t_1)) = J(u_{0})-\int^{t_1}_{0}\|u_t(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\le d, \end{align*}

    which tells us that J(u(t_1)) = d and u_t\equiv 0 for 0<t<t_1 . On the other hand, as u is a global weak solution of problem (1.1)-(1.3), we directly obtain the following equation by multiplying (1.1) by u and integrating over \mathbb{B}

    (u_{t},u )_{\mathbb{B}} +(\nabla_{\mathbb{B}}u_{t}, \nabla_{\mathbb{B}} u )_{\mathbb{B}} +(\nabla_{\mathbb{B}}u, \nabla_{\mathbb{B}}u)_{\mathbb{B}} = \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}, t\ge0.

    Then from the definition of I(u) and I(u(t))>0 for t\in[0,t_1) , it follows that

    \begin{align} (u_{t},u )_{\mathbb{B}} +(\nabla_{\mathbb{B}}u_{t}, \nabla_{\mathbb{B}} u )_{\mathbb{B}} = -I(u) < 0, t\in [0,t_1), \end{align} (5.13)

    that is,

    \frac{{\rm d}}{{\rm d}t}\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = -2I(u) < 0, t\in[0,t_1),

    which contradicts with u_t\equiv 0 for 0<t<t_1 . Hence, we conclude I(u(t))>0 for t\ge0 .

    Since u(t)\in \mathcal{W} for t\ge0 of the case J(u_0)\le d , namely I(u(t))>0 for 0\leq t<\infty , and it is easy to derive that the energy functional J(u(t)) is non-increasing from the energy identity (3.3), then we have

    \begin{align} J(u_0)&\geq J(u(t)) \\ & = \left(\frac{1}{2}-\frac{1}{p+1}\right)\|\nabla_{\mathbb{B}}u\|^{2}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+ \frac{1}{p+1}I(u) \\ & > \frac{p-1}{2(p+1)(c^2+1)}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, \end{align} (5.14)

    where c is the optimal constant in Proposition 2. Moreover,

    \begin{align} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x} &\leq C^{p+1}_{*}\|u\|^{p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ & = C^{p+1}_{*}\left(\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)^{\frac{p-1}{2}}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (5.15)

    Then from (5.14) we define

    \begin{align*} \alpha: = &C^{p+1}_{*}\left(\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)^{\frac{p-1}{2}}\\ < &C^{p+1}_{*}\left(\frac{2(p+1))(c^2+1)}{p-1}J(u_0)\right)^{\frac{p-1}{2}}\\ \le&C^{p+1}_{*}\left(\frac{2(p+1))(c^2+1)}{p-1}d\right)^{\frac{p-1}{2}} = \frac{1}{c^2+1} < 1. \end{align*}

    Hence, taking \beta: = 1-\alpha>0 , we obtain from (5.15) that

    \begin{align*} \int_{\mathbb{B}}|u|^{p+1}\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\le\alpha\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = (1-\beta)\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, \end{align*}

    which gives

    \begin{align} \beta\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\le I(u(t)). \end{align} (5.16)

    On the other hand, by the definition of I(u) we have

    \begin{align} \frac{{\rm d}}{{\rm d}t}\left(\int_{\mathbb{B}}(|u(t)|^{2}+|\nabla_{\mathbb{B}}u(t)|^{2})\frac{{\rm d}x_b}{x_b}{\rm d}\widetilde{x}\right) = -2I(u(t))\leq -2\beta\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (5.17)

    Then by Gronwall's inequality we can obtain

    \begin{align*} \|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} \leq e^{-2\beta t}\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align*}

    Theorem 5.3 (Finite time blow up and lower bound estimate of blowup time). Suppose that u_{0}\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) and J(u_0)\le d . If u_0\in \mathcal{V} , then u(t) blows up in finite time, i.e., there exists a T >0 such that

    \begin{align*} \lim\limits_{t\rightarrow T^{-}}\int^{t}_{0}\|u(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau = +\infty,\ \ \ 0 < t < T. \end{align*}

    Moreover, T is bounded below, which can be estimated by

    \begin{align*} T\geq \frac{\|u_0\|^{-p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-r^{-p+1}}{(p-1)C_*^{p+1}}. \end{align*}

    Proof. Firstly, by Theorem 4.1, we already have the local existence for u_0\in\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B}) , then arguing by contradiction, we suppose that the solution u exists globally in time and here we just consider nontrivial u since the trivial ones do not agree the initial condition, thus are not solutions of problem (1.1)-(1.3). Then we define an auxiliary function as

    \begin{align} M(t): = \int^{t}_{0}\|u(\tau)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau, t\in[0,\infty). \end{align} (5.18)

    For t\in[0,\infty) we can compute its derivative as follows

    \begin{align} M'(t)& = \|u(t)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} = \|u(t)\|^{2}_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\|\nabla_{\mathbb{B}} u(t)\|^{2}_{L^{\frac{n}{2}}_{2}(\mathbb{B})} \end{align} (5.19)

    and further by (1.1) and (3.2) we have

    \begin{align} M''(t) = 2(u_{t}(t), u(t))_{\mathbb{B}}+2(\nabla_\mathbb{B} u_{t}(t),\nabla_\mathbb{B} u(t))_{\mathbb{B}} = -2I(u(t)). \end{align} (5.20)

    Next, we will reveal that the solution actually does not exist globally by showing that M(t) tends to infinity in finite time. However, it is easier to demonstrate that M^{-\gamma}(t) has a zero point, where the exponent -\gamma is a negative constant. Therefore, we take the latter as the proof scheme. Before that, we claim that M^{-\gamma}(t) is concave for sufficiently large t by constructing a differential inequality with M(t) . Recalling (3.3) and (5.3), for t\in[0,\infty) we can obtain by (5.19) and (5.20) that

    \begin{align} M''(t)& = -2(p+1)J(u)+(p-1) \| \nabla_{\mathbb{B}}u(t)\|^{2}_{L^{\frac{n}{2}}_{2}(\mathbb{B})} \\ &\geq2(p+1)\left(\int^{t}_{0}\| u_t\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau-J(u_0)\right)+\frac{(p-1) M'(t)}{c^2+1}, \end{align} (5.21)

    where c is the optimal constant of the Cone Poincar \acute{\rm e} inequality (see Proposition 2). Notice that

    \begin{align*} \left(\int^{t}_{0}(u,u_t)_{\mathbb{B}}{\rm d}\tau\right)^{2} = &\left(\frac{1}{2}\int^{t}_{0}\frac{{\rm d}}{{\rm d}\tau}\| u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right)^{2} \notag \\ = &\frac{1}{4}\left(\|u\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-2\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\| u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+ \|u_0\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right) \\ = &\frac{1}{4}\left((M'(t))^{2}-2M'(t)\| u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\| u_0\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right), \end{align*}

    hence

    \begin{align} (M'(t))^{2} = 4\left(\int^{t}_{0}(u,u_t)_{\mathbb{B}}{\rm d}\tau\right)^{2}+2\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}M'(t)-\| u_0\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (5.22)

    From (5.18), (5.21) and (5.22), making use of the Hölder inequality, we have for t\in[0,\infty)

    \begin{align} &M(t)M''(t)-\frac{p+1}{2}(M'(t))^{2} \\ \geq& \int^{t}_{0}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\left(2(p+1)\left( \int^{t}_{0}\|u_t\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau-J(u_0)\right)+\frac{(p-1)M'(t)}{c^2+1}\right) \\ &-\frac{p+1}{2}\left(4\left(\int^{t}_{0}(u,u_t)_{\mathbb{B}}{\rm d}\tau\right)^{2}+2M'(t)\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} -\|u_0\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right) \\ = &2(p+1)\left(\int^{t}_{0}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\int^{t}_{0}\|u_t\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau -\left(\int^{t}_{0}(u,u_t)_{\mathbb{B}}{\rm d}\tau\right)^{2}\right) \\ &-2(p+1)J(u_0)M(t)+\frac{(p-1)M(t)M'(t)}{c^2+1} \\ &-(p+1)M'(t)\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\frac{p+1}{2}\|u_0\|^{4}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ \geq&-2(p+1)J(u_{0})M(t)+\frac{(p-1)M(t)M'(t)}{c^2+1}-(p+1)M'(t)\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (5.23)

    To estimate the right hand side of (5.23), we claim first that u(t)\in \mathcal{V} for t\in [0,\infty) , and the proof of this claim is similar with the proof in the second step of Theorem 5.2, for which we omit it here. Therefore, we have I(u(t))<0 for t\ge0 , then by the density of the real number we can derive that there exists a \delta>0 such that

    \begin{align} &M''(t) = -2I(u(t))\ge\delta,t\ge0, \end{align} (5.24)

    Integrating over the both sides of (5.24) from 0 to t , we get

    \begin{align} &M'(t)\ge\delta t+M'(0)\ge\delta t,t\ge0. \end{align} (5.25)

    By applying the same operation that we took for (5.24), we have

    \begin{align} &M(t)\ge\delta t^2+M(0)\ge\delta t^2,t\ge0. \end{align} (5.26)

    Thus by (5.23), for sufficiently large t>0 we have

    \begin{align} &M(t)M''(t)-\frac{p+1}{2}(M'(t))^{2}\\ \ge&-2(p+1)J(u_{0})M(t)+\frac{(p-1)M(t)M'(t)}{c^2+1}-(p+1)M'(t)\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ = &M(t)\left(\frac{p-1}{2(c^2+1)}M'(t)-2(p+1)J(u_0)\right)\\ & \quad +M'(t)\left(\frac{p-1}{2(c^2+1)}M(t)-(p+1)\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right) > 0. \end{align} (5.27)

    By direct computation we can get that

    \begin{align} (M^{-\gamma}(t))'' = -\gamma M^{-\gamma-2}(t)\left(M(t)M''(t)-(\gamma+1)(M'(t))^{2}\right). \end{align} (5.28)

    Let \gamma = \frac{p-1}{2}>0 and N(t): = M^{-\gamma}(t) , then (5.28) implies that N''(t)<0 for sufficiently large t>0 according to (5.27) and the facts that \gamma>0 and M(t)\ge0 , which implies that N(t) is concave for sufficiently large t>0 .

    Note first that there exists a sufficiently small \tilde{t}>0 such that the non-trivial u exists locally for t\in(0,\tilde{t}] by Theorem 4.1, which claims that N(t) = M^{-\gamma}(t)>0 for t\in(0,\tilde{t}] . Then we assert that N(t) is decreasing for t\in(0,\infty) . Recalling that M'(t)>0 for t\in(0,\infty) , hence

    \begin{align*} N'(t) = -\gamma M^{-\gamma-1}M'(t) < 0,t\in(0,\infty), \end{align*}

    which makes N(t) keep falling for t>0 and guarantees that the zero point can be reached.

    In fact, the case of the asymptote won't occur since N(t) is concave for sufficiently large t>0 , which forces the curve of N(t) to hit the t -axis for sufficiently large t>0 and ensures the existence of the zero point. Hence there exists a 0< T<\infty such that

    \begin{align*} N(t) = M(t)^{-\gamma}\rightarrow0,t\rightarrow T^-, \end{align*}

    which means also M(t) = \int^{t}_{0}\|u\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\rightarrow \infty \ \text{as}\ \ t\rightarrow T^{-}, and contradicts the hypothesis of global existence, thus u blows up in finite time.

    Next, we seek the lower bound of the blow up time. By (5.20) we have

    \begin{align} M''(t) = -2I(u(t)) = -2\|\nabla_{\mathbb{B}}u(t)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+2\|u(t)\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}. \end{align} (5.29)

    Then Proposition 2, (5.19) and (5.29) imply

    \begin{align} M''(t)\leq2C^{p+1}_*(M'(t))^{\frac{p+1}{2}}. \end{align} (5.30)

    By (5.25) we have already known M'(t)>0 for t\in[0,T) , thus we can divide (5.30) by (M'(t))^{\frac{p+1}{2}} and get the following inequality,

    \begin{align} \frac{M''(t)}{(M'(t))^{\frac{p+1}{2}}}\leq 2C^{p+1}_*. \end{align} (5.31)

    Integrating the inequality (5.31) from 0 to t , we have

    \begin{align} (M'(0))^{-\frac{p-1}{2}}-(M'(t))^{-\frac{p-1}{2}}\leq (p-1)C_*^{p+1}t. \end{align} (5.32)

    Let t\rightarrow T in (5.32), since M'(t) = \|u(t)\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}>r^2 by (ii) of Lemma 3.3, we can conclude that

    \begin{align*} T\geq \frac{\|u_0\|^{-p+1}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-r^{-p+1}}{(p-1)C_*^{p+1}}. \end{align*}

    In this section we give the finite time blow up result for the sup-critical initial energy case, for which we introduce the following three lemmas first.

    Lemma 6.1 ([26]). Suppose that a positive, twice-differentiable function \psi(t) satisfies

    \begin{align*} \psi''(t)\psi(t)-(1+\theta)(\psi'(t))^2\ge0,\ t > 0,\psi(t)\in C^2,\psi(t) > 0 \end{align*}

    where \theta>0 is a constant. If \psi(0)>0 and \psi'(0)>0 , then there exists a t_1 with 0<t_1\le \frac{\psi(0)}{\theta\psi'(0)} such that \psi(t) tends to infinity as t\rightarrow t_1 .

    Lemma 6.2. Suppose that J(u_0)>0 and u_0\in\mathcal{V} , then the map

    \begin{align*} t\mapsto \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} \end{align*}

    increases strictly while u(t)\in \mathcal{V} for t\in[0,T_0] , where T_0 is a positive constant.

    Proof. Firstly, an auxiliary function is defined as follows

    \begin{align} F(t): = \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (6.1)

    Then it follows from (1.1) that

    \begin{align} \begin{split} F'(t) = & 2(u_{t}(t), u(t))_{\mathbb{B}}+2(\nabla_{\mathbb{B}}u_{t}(t),\nabla_{\mathbb{B}}u(t))_{\mathbb{B}} = -2I(u). \end{split} \end{align} (6.2)

    Hence by u(t)\in \mathcal{V} for t\in[0,T_0] we get

    \begin{align} F'(t) > 0,t\in[0, T_0], \end{align} (6.3)

    which implies that the map

    t\mapsto \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}

    is strictly increasing for t\in[0,T_0] .

    Lemma 6.3 (Invariant set \mathcal{V} ). Assume that u_0 satisfies

    \begin{equation} \frac{(p-1)}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > J(u_0) > 0, \end{equation} (6.4)

    where c is the optimal constant in Proposition 2, then u(t)\in \mathcal{V} for t\in[0,T) , where T\le +\infty is the maximal existence time of the solution.

    Proof. First we claim that u_0\in \mathcal{V} by (6.4), i.e. u(t)\in \mathcal{V} for t = 0 . By the definition of J(u(t)) (3.1) we have

    \begin{align} J(u_0) = &\frac{1}{2}\|\nabla_{\mathbb{B}} u_0\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} - \frac{1}{p+1}\|u_0\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ = &\left(\frac{1}{2}- \frac{1}{p+1}\right)\|\nabla_{\mathbb{B}} u_0\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\frac{1}{p+1}I(u_0)\\ \ge&\frac{(p-1)}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\frac{1}{p+1}I(u_0), \end{align} (6.5)

    thus by (6.4) and (6.5) we can drive that I(u_0)<0 , namely u_0\in \mathcal{V} .

    Then we prove that u(t)\in \mathcal{V} for t\in(0,T) . Since I(u_0)<0, we have I(u(t))<0 on [0,t_0) for sufficiently small t_0>0 according to the continuity of I(u(t)) with respect to t . Arguing by contradiction, if I(u(t)) doesn't remain negative on [0,T) , then there must exist a t_1\in [t_0,\infty) such that I(u(t_1)) = 0 for the first time. It follows from Lemma 6.2 that the map t\mapsto \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} increases strictly for t\in[0,t_1) , by which and (6.4) we get

    \begin{align} \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > \|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > \frac{2(c^2+1)(p+1)}{(p-1)} J(u_0),\ t\in(0,t_1), \end{align} (6.6)

    then from the continuity of \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} with respect to t , we can obtain

    \begin{align} \|u(t_1)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > \|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} > \frac{2(c^2+1)(p+1)}{p-1} J(u_0). \end{align} (6.7)

    On the other hand, recalling the definition of J(u(t)) and (3.3), we have

    \begin{align} J(u_0) = &J(u(t_1)) +\int_0^{t_{1}} \| u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} {\rm d}\tau\\ \ge &\frac{1}{2}\|\nabla_{\mathbb{B}} u(t_1)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})} - \frac{1}{p+1}\|u(t_1)\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ = &\left(\frac{1}{2}- \frac{1}{p+1}\right)\|\nabla_{\mathbb{B}} u(t_1)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+\frac{1}{p+1}I(u(t_1)), \end{align} (6.8)

    then applying the Cone Poincar \acute{\rm e} inequality (Proposition 2) to (6.8) with I(u(t_1)) = 0 , we have

    \begin{align} J(u_0)&\ge\frac{p-1}{2(p+1)}\|\nabla_{\mathbb{B}} u(t_1)\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ &\ge\frac{(p-1)}{2(c^2+1)(p+1)}\|u(t_1)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})},\\ &\ge\frac{(p-1)}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, \end{align} (6.9)

    which contradicts (6.7), thus we can drive that u(t)\in\mathcal{V} for t\in [0,T) .

    Next we give the blow up results for arbitrary positive initial energy as follows:

    Theorem 6.4 (Finite time Blow up with J(u_0)>0 ). Let u_{0}\in\mathcal{H}^{1,\frac{n+1}{2}}_{2,0}(\mathbb{B}) . If u_0 satisfies (6.4), then u(t) blows up in finite time. Furthermore, the upper bound of blowup time can be estimated by

    \begin{align*} T\le\frac{4\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\sigma}, \end{align*}

    where \sigma: = \frac{p-1}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-J(u_0) .

    Proof. Similar as Theorem 5.2, we also divide the proof into two parts, which are the finite time blow up and the estimating of the upper bound of the blow up time.

    Part Ⅰ: Finite time blow up.

    Firstly, Theorem 4.1 has asserted the local existence of the solution, then arguing by contradiction, we suppose that u(t) exists globally in time, i.e. the maximal existence time T = \infty .

    Now we take a sufficiently small \varepsilon>0 and a positive constant c_0>0 such that

    \begin{align} c_0 > \frac{1}{4}\varepsilon^{-2}\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} \end{align} (6.10)

    and define a new auxiliary function for t\in[0,\infty) as follows

    \begin{align} P(t): = (M(t))^2+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} M(t)+c_0, \end{align} (6.11)

    where M(t) is the auxiliary function defined before in (5.18). Hence for t\in[0,\infty) ,

    \begin{align} P'(t) = \left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M'(t) \end{align} (6.12)

    and

    \begin{align*} P''(t) = \left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)+2(M'(t))^2. \end{align*}

    Set \xi: = 4c_0-\varepsilon^{-2}\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}, then (6.10) indicates \xi>0. Thus for t\in[0,\infty) we have

    \begin{align*} (P'(t))^2& = \left(4(M(t))^2+4\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}M(t)+\varepsilon^{-2}\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)(M'(t))^2\\ & = \left(4(M(t))^2+4\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}M(t)+4c_0-\xi\right)(M'(t))^2\\ & = \left(4P(t)-\xi\right)(M'(t))^2, \end{align*}

    which tells

    \begin{align} 4P(t)(M'(t))^2 = (P'(t))^2+\xi(M'(t))^2,t\in[0,\infty). \end{align} (6.13)

    By (6.13), for t\in[0,\infty) we get

    \begin{align} &2P(t)P''(t)\\ = &2\left(\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)+2(M'(t))^2\right)P(t)\\ = &2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)+4P(t)(M'(t))^2\\ = &2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)+(P'(t))^2+\xi(M'(t))^2. \end{align} (6.14)

    Then from (6.13) and (6.14), for t\in[0,\infty) we can derive that

    \begin{align} &2P(t)P''(t)-(1+\beta)(P'(t))^2\\ = & 2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)+\xi(M'(t))^2-\beta(P'(t))^2\\ = & 2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)+\xi(M'(t))^2-\beta(4P(t)-\xi)(M'(t))^2\\ = & 2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)-4\beta P(t)(M'(t))^2+\xi(1+\beta)(M'(t))^2\\ > & 2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)-4\beta P(t)(M'(t))^2, \end{align} (6.15)

    where \beta>0 is a positive constant that will be determined in the sequel.

    Next, we estimate the term M''(t) . Testing the both sides of (1.1) by u , for t\in[0,\infty) we have

    \begin{align} \frac{1}{2}\frac{{\rm d}}{{\rm d}t}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}& = -I(u)\\ & = -\|\nabla_{\mathbb{B}} u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}+ \|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ & = -2J(u)+\frac{p-1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}. \end{align} (6.16)

    By (6.4), we can take \beta such that

    \begin{align} 1 < \beta < \frac{(p-1)\|u_0\|^{2}_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{2(c^2+1)(p+1)J(u_0)}. \end{align} (6.17)

    Notice by Lemma 6.3 that I(u(t))<0 for t\in[0,+\infty) , then combining (6.16) with (3.3), for t\in[0,\infty) we see

    \begin{align} &\frac{1}{2}\frac{{\rm d}}{{\rm d}t}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ = &-2J(u)+\frac{p-1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ = &2(\beta-1)J(u)-2\beta J(u)+\frac{p-1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ \geq&-2\beta J(u_0)+2\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\frac{p-1}{p+1}\|u\|^{p+1}_{L^{\frac{n}{p+1}}_{p+1}(\mathbb{B})}\\ = &-2\beta J(u_0)+2\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau-\frac{p-1}{p+1}I(u)+\frac{p-1}{p+1}\|\nabla_{\mathbb{B}} u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\\ > &-2\beta J(u_0)+2\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\frac{p-1}{p+1}\|\nabla_{\mathbb{B}} u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}. \end{align} (6.18)

    The application of Cone Poincar \acute{\rm e} inequality asserts

    \begin{align} \frac{p-1}{p+1}\|\nabla_{\mathbb{B}} u\|^2_{L^{\frac{n}{2}}_{2}(\mathbb{B})}\geq \frac{(p-1)}{(c^2+1)(p+1)}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (6.19)

    Putting (6.19) into (6.18), for t\in[0,\infty) we have

    \begin{align} M''(t) = &\frac{{\rm d}}{{\rm d}t}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ > & -4\beta J(u_0)+4\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\frac{2(p-1)}{(c^2+1)(p+1)}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}. \end{align} (6.20)

    By (6.16), H \ddot{\rm o} lder and Young's inequalities, for t\in[0,\infty) we estimate the term (M'(t))^2 as follows

    \begin{align} (M'(t))^2 = &\|u\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\\ = &\left(\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+2\int^t_0(u,u_t)_{\mathbb{B}}{\rm d}\tau\right)^2\\ \leq& \left(\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+2\left(\int^t_0\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right)^{\frac{1}{2}}\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right)^{\frac{1}{2}}\right)^2\\ = &\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+4\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\left(\int^t_0\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right)^{\frac{1}{2}}\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right)^{\frac{1}{2}}\\ &+4M(t)\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\\ \leq& \|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+2\varepsilon\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}M(t)+2\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\\ &+4M(t)\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau. \end{align} (6.21)

    Then from (6.20), (6.21) and (6.15), for t\in[0,\infty) we have

    \begin{align*} &2P''(t)P(t)-(1+\beta)(P'(t))^2\\ > & 2\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)M''(t)P(t)-4\beta P(t)(M'(t))^2\\ > & I_1I_2-I_3I_4, \end{align*}

    where

    \begin{align*} I_1: = &2P(t)\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right),\\ I_2: = &-4\beta J(u_0)+4\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\frac{2(p-1)}{(c^2+1)(p+1)}\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})},\\ I_3: = &4\beta P(t),\\ I_4: = &\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+2\varepsilon\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}M(t)+2\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\\ &+4M(t)\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau. \end{align*}

    Taking \gamma: = \frac{2(p-1)}{(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-4\beta J(u_0), then (6.17) ensures \gamma>0. Choosing \varepsilon such that

    \begin{align*} \varepsilon < \frac{\gamma}{2\beta \|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}, \end{align*}

    then recalling that \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})} is decreasing with respect to t by Lemma 6.2, for t\in[0,\infty) we obtain

    \begin{align*} &2P''(t)P(t)-(1+\beta)(P'(t))^2\\ > &I_1I_2-I_3I_4\\ = &I_1\left(4\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\gamma \right)-I_3I_4\\ > &I_1\left(4\beta \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+2\beta\varepsilon\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)-I_3I_4\\ = &4\beta P(t)\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)\left(2 \int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\varepsilon\|u_0\|^2\right)-I_3I_4\\ = &I_3\left(\left(2M(t)+\varepsilon^{-1}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)\left(2\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\varepsilon\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)-I_4\right)\\ = &0. \end{align*}

    Thus

    \begin{align*} P''(t)P(t)-\frac{1+\beta}{2}(P'(t))^2 > 0,t\in[0,\infty). \end{align*}

    Since P(0) = c_0>\frac{1}{4}\varepsilon^{-2}\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}>0 and P'(0) = \varepsilon^{-1}\|u_0\|^4_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}>0, by Lemma 6.1 we can conclude that there exists a 0<T<\infty such that

    \begin{align*} \lim\limits_{t\rightarrow T}P(t) = +\infty. \end{align*}

    According to the definition of P(t) , i.e. (6.11), we can conclude that

    \begin{align*} \lim\limits_{t\rightarrow T}M(t) = +\infty, \end{align*}

    which claims the blow up of the solution.

    Part Ⅱ: Upper bound of the blow up time

    In Part I we have draw the conclusion that the maximal existence time T is finite, next, we estimate the upper bound of blowup time. For t\in[0,T) , we define

    \begin{align} \psi(t): = \int^t_0\|u(\tau)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+(T-t)\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\mu(t+\nu)^2, \end{align} (6.22)

    where \mu>0 and \nu>0 are constants, which will be determined later in the process of argumentation. Thus we have

    \begin{align} \psi'(t) = \|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+2\mu(t+\nu), \end{align} (6.23)

    then it follows from (3.3), Proposition 2 and (6.4) that

    \begin{align} \psi''(t) = &-2I(u(t))+2\mu\\ \ge&(p-1)\|\nabla_{\mathbb{B}} u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-2(p+1)J(u)\\ \ge&\frac{p-1}{c^2+1}\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-2(p+1)J(u)\\ = &\frac{p-1}{c^2+1}\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-2(p+1)J(u_0)+2(p+1)\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\\ \ge&2(p+1)\left(\frac{p-1}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-J(u_0)+\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\right) > 0. \end{align} (6.24)

    Furthermore, we obtain \psi'(t)\ge\psi'(0) = 2\mu(t+\nu)>0 via (3.3), which implies \psi(t)\ge\psi(0) = T\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\mu\nu^2>0 for all t\in[0,T) .

    On the other hand, we can derive that

    \begin{align} -\frac{1}{4}(\psi'(t))^2 = &-\left(\frac{1}{2}\left(\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}- \|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)+\mu(t+\nu)\right)^2\\ = &{\left(\int^t_0\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu(t+\nu)^2\right)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right)}\\ &-{\left(\frac{1}{2}\left(\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)+\mu(t+\nu)\right)^2}\\ &-\left(\psi(t)-(T-t)\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right)\\ = &I_5-I_6-\left(\psi(t)-(T-t)\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right), \end{align} (6.25)

    where

    \begin{align*} &I_5: = {\left(\int^t_0\|u\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu(t+\nu)^2\right)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right)},\\ &I_6: = {\left(\frac{1}{2}\left(\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)+\mu(t+\nu)\right)^2}. \end{align*}

    To estimate (6.25) clearly, we will show that I_5-I_6>0 ,

    \begin{align*} I_5-I_6 = &I_5-\left(\frac{1}{2}(\|u(t)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})})+\mu(t+\nu)\right)^2\\ = &I_5-\left(\frac{1}{2}\int^t_0\frac{d}{d\tau}\|u(\tau)\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu(t+\nu)\right)^2\\ = &I_5-\left(\int^t_0(u,u_t)_{\mathbb{B}}{\rm d}\tau+\mu(t+\nu)\right)^2\\ \ge& I_5-\left(\int^t_0\|u(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\|u_t(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu(t+\nu)\right)^2, \end{align*}

    then by the Cauchy-Schwartz inequality we get

    \begin{align*} I_5-I_6 \ge& I_5-\left(\int^t_0\|u(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau\int^t_0\|u_t(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu(t+\nu)\right)^2\\ = & I_5-\left(k_1(t)k_2(t)+\mu(t+\nu)\right)^2\\ = &\left((k_1(t))^2+\mu(t+\nu)^2\right)\left((k_2(t))^2+\mu\right) -\left(k_1(t)k_2(t)+\mu(t+\nu)\right)^2\\ = &\left(\sqrt{\mu}k_1(t)\right)^2-2\sqrt{\mu}k_1(t)\sqrt{\mu}(t+\nu)k_2(t) +\left(\sqrt{\mu}(t+\nu)k_2(t)\right)^2\\ = &\left(\sqrt{\mu}k_1(t)-\sqrt{\mu}(t+\nu)k_2(t)\right)^2\\ \ge&0, \end{align*}

    where k_1(t): = \int^t_0\|u(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau and k_2(t): = \int^t_0\|u_t(\tau)\|_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau . Hence,

    \begin{align} -(\psi'(t))^2\ge&-4\left(\psi(t)-(T-t)\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right)\\ \ge&-4\psi(t)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right). \end{align} (6.26)

    Then by (6.22), (6.24) and (6.26), we achieve

    \begin{align} &\psi(t)\psi''(t)-\frac{p+1}{2}(\psi'(t))^2\\ \ge&\psi(t)\left(\psi''(t)-2(p+1)\left(\int^t_0\|u_t\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}{\rm d}\tau+\mu\right)\right)\\ \ge&2(p+1)\psi(t)\left(\frac{p-1}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-J(u_0)-\mu\right). \end{align} (6.27)

    Let \sigma: = \frac{p-1}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-J(u_0) , then we can take a sufficiently small \mu\in\left(0,\sigma\right] such that

    \begin{align} \frac{p-1}{2(c^2+1)(p+1)}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}-J(u_0)-\mu\ge0, \end{align} (6.28)

    hence for t\in[0,\infty) , from (6.27) and (6.28) we have

    \begin{align*} \psi(t)\psi''(t)-\frac{p+1}{2}(\psi'(t))^2\ge0. \end{align*}

    It is easy to verify that \psi(0) = T\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}+\mu\nu^2>0 , \psi'(0)\ge2\mu\nu>0 , which implies that the conditions of Lemma 6.1 are satisfied. Then applying Lemma 6.1, we can derive that

    \begin{align} T\le\frac{2\psi(0)}{(p-1)\psi'(0)} \le\frac{\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\mu\nu}T+\frac{\nu}{p-1}. \end{align} (6.29)

    Let \nu be large enough such that

    \begin{align} \nu\in\left(\frac{\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\mu},+\infty\right), \end{align} (6.30)

    then it follows from (6.29) that

    \begin{align} T\le\frac{\mu\nu^2}{(p-1)\mu\nu-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}. \end{align} (6.31)

    Note that \mu\in (0,\sigma] by (6.28). To estimate the upper bound of the blow up time, we can revise the range of \mu according to (6.30) and define the following set to describe the pair (\mu,\nu)

    \begin{align*} \mathfrak{M}: = \left\{(\nu,\mu)\ \Bigg|\ \nu\in\left(\frac{\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\sigma},+\infty\right), \mu\in \left(\frac{\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\nu},\sigma\right]\right\}. \end{align*}

    From the above discussions, we get

    \begin{align} T\le\inf\limits_{(\mu,\nu)\in\mathfrak{M}}\frac{\mu\nu^2}{(p-1)\mu\nu-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}. \end{align} (6.32)

    Next we define

    \begin{align*} f(\mu,\nu): = \frac{\mu\nu^2}{(p-1)\mu\nu-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}, \end{align*}

    where (\mu,\nu)\in\mathfrak{M} . Thus the estimation of the blow-up time turns into seeking the minimal value of f(\mu,\nu) . First differentating f(\mu,\nu) with respect to \mu , we have

    \begin{align*} \frac{\partial}{\partial\mu}f(\mu,\nu) = \frac{-\nu^{2}\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{\left((p-1) \mu \nu-\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}\right)^{2}} < 0, \end{align*}

    which means that f(\mu,\nu) is decreasing with respect to \mu , for which we have

    \begin{align*} \inf\limits_{(\mu,\nu)\in\mathfrak{M}}f(\mu,\nu) = \inf\limits_{\nu}f(\sigma,\nu), \end{align*}

    where \nu\in\left(\frac{\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\sigma},+\infty\right). Then we compute the partial derivative of f(\sigma,\nu) with respect to \nu and let it equal zero, we can obtain the minimum

    \begin{align*} \nu_{\min} = \frac{2\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1) \sigma}, \end{align*}

    thus we have

    \begin{align*} \inf\limits_{(\mu,\nu)\in\mathfrak{M}}f(\mu,\nu) = f(\sigma,\nu_{\min}) = \frac{4\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\sigma}, \end{align*}

    which means also

    \begin{align*} T\le\frac{4\|u_0\|^2_{\mathcal{H}^{1,\frac{n}{2}}_{2,0}(\mathbb{B})}}{(p-1)\sigma}. \end{align*}

    In the first version [37] finished in 2016, XP Liu obtained the threshold conditions of the global existence and nonexistence for solutions together with the asymptotic behavior at both subcritical and critical initial energy levels. But this version [37] did not consider the local existence nor the dynamic behavior at the sup-critical initial energy level. Then, the authors realized that the framework of a family of potential wells brought some problems by unnecessary technical complexities without receiving the benefits expected. After XP Liu left HEU in 2016, YT Wang and YX Chen continued to reconsider problem (1.1)-(1.3) and reconstructed the conclusions again in the framework of a single potential well to replace the family version, and finally completed this second version in 2020. This new version of the present paper not only gives the same results with those in the first version [37] by a single potential well, but also proves the local existence and the finite time blow up at sup-critical initial energy level with the estimation of the blowup time at all initial energy levels.

    The authors are grateful to Professor Hua Chen in Wuhan University, who strongly encouraged this work when he visited Harbin Engineering University in 2017 and 2018.

    [1] Pope CA, Ezzati M, Dockery DW (2009) Fine-particulate air pollution and life expectancy in the United States. N Engl J Med 360: 376–386. doi: 10.1056/NEJMsa0805646
    [2] Adams PJ, Seinfeld JH (2003) Disproportionate impact of particulate emissions on global cloud condensation nuclei concentrations. Geophys Res Lett 30: 335–343.
    [3] Kanakidou M, Seinfeld JH, Pandis SN, et al. (2005) Organic aerosol and global climate modeling: a review. Atmos Chem Phys 5: 1053–1123. doi: 10.5194/acp-5-1053-2005
    [4] Spracklen DV, Carslaw KS, Poschl U, et al. (2011) Global cloud condensation nuclei influenced by carbonaceous combustion aerosol. Atmos Chem Phys 11: 9067–9087. doi: 10.5194/acp-11-9067-2011
    [5] Karner AA, Eisinger DS, Niemeier DA (2010) Near-roadway air quality: synthesizing the findings from real-world data. Environ Sci Technol 44: 5334–5344. doi: 10.1021/es100008x
    [6] Padro-Martinez LT, Patton AP, Trull JB, et al. (2012) Mobile monitoring of particle number concentration and other traffic-related air pollutants in a near-highway neighborhood over the course of a year. Atmos Environ 61: 253–264. doi: 10.1016/j.atmosenv.2012.06.088
    [7] Perkins JL, Padro-Martinez LT, Durant JL (2013) Particle number emission factors for an urban highway tunnel. Atmos Environ 74: 326–337. doi: 10.1016/j.atmosenv.2013.03.046
    [8] Ostro B, Goldberg D, Reynolds P, et al. (2015) Associations of mortality with long-term exposures to fine and ultrafine particles, species and sources: results from the California Teachers Study cohort. Environ Health Perspect 123: 549–556. doi: 10.1289/ehp.1408565
    [9] Smith SJ, Bond TC (2014) Two hundred fifty years of aerosols and climate: the end of the age of aerosols. Atmos Chem Phys 14: 537–549. doi: 10.5194/acp-14-537-2014
    [10] Bond TC (2001) Spectral dependence of visible light absorption by carbonaceous particles emitted from coal combustion. Geophys Res Lett 28: 4075–4078. doi: 10.1029/2001GL013652
    [11] Lack DA, Cappa CD (2010) Impact of brown and clear carbon on light absorption enhancement, single scatter albedo and adsorption wavelength dependance of black carbon. Atmos Chem Phys 10: 4207–4220. doi: 10.5194/acp-10-4207-2010
    [12] Smith SJ, Mizrahi A (2013) Near-term climate mitigation by short-lived forcers. Proc Natl Acad Sci 110: 14202–14206. doi: 10.1073/pnas.1308470110
    [13] Dallmann TR, Harley RA, Kirchstetter TW (2011) Effects of diesel particle filter retrofits and accelerated fleet turnover on drayage truck emissions at the Port of Oakland. Environ Sci Technol 45: 10773–10779. doi: 10.1021/es202609q
    [14] Dallmann TR, Kirchstetter TW, DeMartini SJ, et al. (2013) Quantifying on-road emissions from gasoline-powered motor vehicles: accounting for the presence of medium- and heavy-duty diesel trucks. Environ Sci Technol 47: 13873–13881. doi: 10.1021/es402875u
    [15] May AA, Nguyen NT, Presto AA, et al. (2014) Gas- and particle-phase primary emissions from in-use, on-road gasoline and diesel vehicles. Atmos Environ 88: 247–260. doi: 10.1016/j.atmosenv.2014.01.046
    [16] Robert MA, Bergen SV, Kleeman MJ, et al. (2007) Size and composition distributions of particulate matter emissions: Part 1 - light-duty gasoline vehicles. J Air Waste Manage 57: 1414– 1428. doi: 10.3155/1047-3289.57.12.1414
    [17] Robert MA, Kleeman MJ, Jakober CA (2007) Size and composition distributions of particulate matter emissions: Part 2 - heavy-duty diesel vehicles. J Air Waste Manage 57: 1429–1438. doi: 10.3155/1047-3289.57.12.1429
    [18] Schauer JJ, Kleeman MJ, Cass GR, et al. (1999) Measurement of emissions from air pollution sources. 2. C2 through C30 organic compounds from medium duty diesel trucks. Environ Sci Technol 33: 1578–1587.
    [19] Schauer JJ, Kleeman MJ, Cass GR, et al. (2002). Measurement of emissions from air pollution sources. 5. C1 - C32 organic compounds from gasoline-powered motor vehicles. Environ Sci Technol 36: 1169–1180.
    [20] Adams PJ, Seinfeld JH (2002) Predicting global aerosol size distributions in general circulation models. J Geophys Res 107: 4370 – 4392. doi: 10.1029/2001JD001010
    [21] Pierce JR, Theodoristi G, Adams PJ, et al. (2009) Parameterization of the effect of sub-grid scale aerosol dynamics on aerosol number emission rates. J Aerosol Sci 40: 285–293. doi: 10.1016/j.jaerosci.2008.11.005
    [22] Sahlodin AM, Sotudeh-Gharebagh R, Zhu Y (2007) Modeling of dispersion near roadways based on the vehicle-induced turbulence concept. Atmos Environ 41: 92–102. doi: 10.1016/j.atmosenv.2006.08.004
    [23] Wang YJ, Zhang KM (2009) Modeling near-road air quality using a computational fluid dynamics model, CFD-VIT-RIT. Environ Sci Technol 43: 7778–7783. doi: 10.1021/es9014844
    [24] Zhang KM, Wexler AS (2004) Evolution of particle number distribution near roadways. part 1: analysis of aerosol dynamics and its implications for engine emission measurements. Atmos Environ 38: 6643–6653.
    [25] Zhang KM, Wexler AS, Zhu YF, et al. (2004) Evolution of particle number distribution near roadways. part ii: the 'road-to-ambient' process. Atmos Environ 38: 6655–6665.
    [26] Zhu YF, Hinds WC, Kim S, et al. (2002) Study of ultrafine particles near a major highway with heavy-duty diesel traffic. Atmos Environ 36: 4323–4335. doi: 10.1016/S1352-2310(02)00354-0
    [27] Zhu YF, Hinds WC, Kim S, et al. (2002) Concentration and size distribution of ultrafine particles near a major highway. J Air Waste Manage Assoc 52: 1032–1042. doi: 10.1080/10473289.2002.10470842
    [28] Zhu YF, Kuhn T, Mayo P, et al. (2006) Comparison of daytime and nighttime concentration profiles and size distributions of ultrafine particles near a major highway. Environ Sci Technol 40: 2531–2536. doi: 10.1021/es0516514
    [29] Zhu YF, Hinds WC, Shen S, et al. (2004) Seasonal trends of concentration and size distribution of ultra fine particles near major highways in los angeles. Aerosol Sci Tech 38: 5–13. doi: 10.1080/02786820390229156
    [30] Huang L, Gong SL, Gordon M, et al. (2014) Aerosol-CFD modelling of ultrafine and black carbon particle emission, dilution, and growth near roadways. Atmos Chem Phys Discuss 14: 12235–12278. doi: 10.5194/acpd-14-12235-2014
    [31] Biswas S, Ntziachristos L, Moore KF, et al. (2007) Particle volatility in the vicinity of a freeway with heavy-duty diesel traffic. Atmos Environ 41: 3479–3493. doi: 10.1016/j.atmosenv.2006.11.059
    [32] Salmond JA, Williams DE, Laing G, et al. (2012) The influence of vegetation on the horizontal and vertical distribution of pollutants in a street canyon. Sci Total Environ 443: 287–298.
    [33] Neft I, Scungio M, Culver N, et al. (2016) Simulations of aerosol filtration by vegetation: Validation of existing models with available lab data and application to near-roadway scenario. Aerosol Sci Tech 50: 937–946. doi: 10.1080/02786826.2016.1206653
    [34] Salmond JA, Williams DE, Laing G, et al. (2013) The influence of vegetation on the horizontal and vertical distribution of pollutants in a street canyon. Sci Total Environ 443: 287–298. doi: 10.1016/j.scitotenv.2012.10.101
    [35] Scungio M, Arpino F, Stabile L, et al. (2013) Numerical simulation of ultrafine particle dispersion in urban street canyons with the spalart-allmaras turbulence model. Aerosol Air Qual Res 13: 1423–1437. doi: 10.4209/aaqr.2012.11.0306
    [36] Steffens JT, Wang YJ, Zhang KM (2012) Exploration of effects of a vegetation barrier on particle size distributions in near-road environment. Atmos Environ 50: 120–128. doi: 10.1016/j.atmosenv.2011.12.051
    [37] Steffens JT, Heist DK, Perry SG, et al. (2013) Modeling the effects of a solid barriers on pollutant dispersion under various atmospheric stability conditions. Atmos Environ 69: 76–85. doi: 10.1016/j.atmosenv.2012.11.051
    [38] Chang CH, Meroney RN (2003) Concentration and flow distributions in urban street canyons: wind tunnel and computational data. J Wind Eng Ind Aerod 91: 1141–1154. doi: 10.1016/S0167-6105(03)00056-4
    [39] Solazzo E, Cai X, Vardoulakis S (2008) Modelling wind flow and vehicle-induced turbulence in urban streets. Atmos Environ 42: 4918–4931. doi: 10.1016/j.atmosenv.2008.02.032
    [40] Hagler GSW, Tang W, Freeman MJ, et al. (2011) Model evaluation of roadside barrier impact on near-road air pollution. Atmos Environ 45: 2522–2530. doi: 10.1016/j.atmosenv.2011.02.030
    [41] Hanna SR, Tehranian S, Carissimo B, et al. (2002) Comparisons of model simulations with observations of mean flow and turbulence within simple obstacle arrays . Atmos Environ 36: 5067–5079. doi: 10.1016/S1352-2310(02)00566-6
    [42] Gowardhan AA, Pardyjak ER, Senocak I, et al. (2011) A CFD-based wind solver for an urban fast response transport and dispersion model. Environ Fluid Mech 11: 439–464. doi: 10.1007/s10652-011-9211-6
    [43] Kota SH, Ying Q, Zhang Y (2013) Simulating near-road reactive dispersion of gaseous air pollutants using a three-dimensional Eulerian model. Sci Total Environ 454–455: 348–357.
    [44] Wang YJ, DenBleyker A, McDonald-Buller E, et al. (2011) Modeling the chemical evolution of nitrogen oxides near roadways. Atmos Environ 45: 43–52. doi: 10.1016/j.atmosenv.2010.09.050
    [45] Wang YJ, Nguyen MT, Steffens JT, et al. (2013) Modeling multi-scale aerosol dynamics and micro-environment air quality near a large highway intersection using the ctag model. Sci Total Environ 443: 375–386. doi: 10.1016/j.scitotenv.2012.10.102
    [46] Wang YJ, Zhang KM (2012) Coupled turbulence and aerosol dynamics modeling of vehicle exhaust plumes using the CTAG model. Atmos Environ 59: 284–293. doi: 10.1016/j.atmosenv.2012.04.062
    [47] Richards K, Senecal P, Pomraning E (2013) A Three-Dimensional Computational Fluid Dynamics Program for Transient or Steady State Flow with Complex Geometries. Convergent Science, Inc. Convergent Science, Inc., Middleton,WI (version 2.1.0 ed.).
    [48] Issa R (1986) Solution of the implicitly discretised fluid flow equations by operator-splitting. J Comput Phys 62: 40–65. doi: 10.1016/0021-9991(86)90099-9
    [49] Tzivion S, Feingold G, Levin Z (1987) An efficient numerical solution to the stochastic collection equation. J Atmos Sci 44: 3139–3149. doi: 10.1175/1520-0469(1987)044<3139:AENSTT>2.0.CO;2
    [50] Tzivion S, Feingold G, Levin Z (1989) The evolution of raindrop spectra, part ii, collisional collection/breakup and evaporation in a rainshaft. J Atmos Sci 46: 3312–3328. doi: 10.1175/1520-0469(1989)046<3312:TEORSP>2.0.CO;2
    [51] Lee YH, Adams PJ (2012) A fast and efficient version of the two-moment aerosol sectional (TOMAS) global aerosol microphysics model. Aerosol Sci Tech 46: 678–689. doi: 10.1080/02786826.2011.643259
    [52] Shrivastava MK, Stanier CO, Robinson AL (2006) Modeling semi volatile organic aerosol mass emissions from combustion systems. Environ Sci Technol 40: 2671–2677. doi: 10.1021/es0522231
    [53] Singh S, Adams PJ, Misquitta A, et al. (2014) Computational analysis of particle nucleation in dilution tunnels: Effects of flow configuration and tunnel geometry. Aerosol Sci Tech 48: 638–648.  doi: 10.1080/02786826.2014.910291
    [54] Rader DJ, McMurry PH, Smith S (1987) Evaporation rates of monodispersed organic aerosols in the 0.02- to 0.2-μm-diameter range. Aerosol Sci Tech 6: 247–260.
    [55] May AA, Presto AA, Hennigan CJ, et al. (2013) Gas-particle partitioning of primary organic aerosol emissions: (1) gasoline vehicle exhaust. Atmos Environ 77: 128–139. doi: 10.1016/j.atmosenv.2013.04.060
    [56] May AA, Presto AA, Hennigan CJ, et al. (2013) Gas-particle partitioning of primary organic aerosol emissions: (2) diesel vehicles. Environ Sci Technol 47: 8288–8296. doi: 10.1021/es400782j
    [57] Kalthoff N, Baumer D, Corsmeier U, et al. (2005) Vehicle-induced turbulence near a motorway. Atmos Environ 39: 5737–5749. doi: 10.1016/j.atmosenv.2004.06.048
    [58] Norris SE, Richards PJ (2011) Appropriate boundary conditions for computational wind engineering models revisited. J Wind Eng Ind Aerod 99: 257–266. doi: 10.1016/j.jweia.2010.12.008
    [59] Zhang KM,Wexler AS, Niemeier DA, et al. (2005) Evolution of particle number distribution near roadways. part iii: traffic analysis and on-road size resolved particulate emission factors. Atmos Environ 39: 4155–4166.
    [60] Lipsky EM, Robinson AL (2005) Design and evaluation of a portable dilution sampling system for measuring fine particle emissions from combustion systems. Aerosol Sci Tech 39: 542–553. doi: 10.1080/027868291004850
    [61] Jung JG, Pandis SN, Adams PJ (2008) Evaluation of nucleation theories in sulfur-rich environment. Aerosol Sci Technol 42: 495–504. doi: 10.1080/02786820802187085
    [62] Zhao Y, Nguyen NT, Presto AA, et al. (2015) Intermediate volatility organic compound emissions from on-road diesel vehicles: chemical composition, emission factors, and estimated secondary organic aeros production. Environ Sci Technol 49: 11516–11526. doi: 10.1021/acs.est.5b02841
    [63] Zhao Y, Nguyen NT, Presto AA, et al. (2016) Intermediate volatility organic compound emissions from on-road gasoline vehicles and small o_-road gasoline engines. Environ Sci Technol 50: 4554–4563. doi: 10.1021/acs.est.5b06247
    [64] Gkatzelis GI, Papanastasiou DK, Florou K, et al. (2016) Measurement of nonvolatile particle number size distribution. Atmos Meas Tech 9: 103–114. doi: 10.5194/amt-9-103-2016
  • This article has been cited by:

    1. Ngo Tran Vu, Dao Bao Dung, Huynh Thi Hoang Dung, General decay and blow-up results for a class of nonlinear pseudo-parabolic equations with viscoelastic term, 2022, 516, 0022247X, 126557, 10.1016/j.jmaa.2022.126557
  • Reader Comments
  • © 2018 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(4873) PDF downloads(756) Cited by(3)

Figures and Tables

Figures(11)  /  Tables(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog