This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Bouss...This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).展开更多
We study a spatiotemporal EIT problem with a dynamical boundary condition for the fractional Dirichlet-to-Neumann operator with a critical exponent.There are three major ingredients in this paper.The first is the fini...We study a spatiotemporal EIT problem with a dynamical boundary condition for the fractional Dirichlet-to-Neumann operator with a critical exponent.There are three major ingredients in this paper.The first is the finite time blowup and the decay estimate of the global solution with a lower-energy initial value.The second ingredient is the L^(q)(2 ≤q <∞) estimate of the global solution applying the Moser iteration,which allows us to show that any global solution is a classical solution.The third,which is the main ingredient of this paper,explores the long time asymptotic behavior of global solutions close to the stationary solution and the bubbling phenomenons by means of a concentration compactness principle.展开更多
In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-...In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.展开更多
Motivated by recent advances made in the study of dividend control and risk management problems involving the U.S.bankruptcy code,in this paper we follow[44]to revisit the De Finetti dividend control problem under the...Motivated by recent advances made in the study of dividend control and risk management problems involving the U.S.bankruptcy code,in this paper we follow[44]to revisit the De Finetti dividend control problem under the reorganization process and the regulator's intervention documented in U.S.Chapter 11 bankruptcy.We do this by further accommodating the fixed transaction costs on dividends to imitate the real-world procedure of dividend payments.Incorporating the fixed transaction costs transforms the targeting optimal dividend problem into an impulse control problem rather than a singular control problem,and hence computations and proofs that are distinct from[44]are needed.To account for the financial stress that is due to the more subtle concept of Chapter 11 bankruptcy,the surplus process after dividends is driven by a piece-wise spectrally negative Lévy process with endogenous regime switching.Some explicit expressions of the expected net present values under a double barrier dividend strategy,new to the literature,are established in terms of scale functions.With the help of these expressions,we are able to characterize the optimal strategy among the set of admissible double barrier dividend strategies.When the tail of the Lévy measure is log-convex,this optimal double barrier dividend strategy is then verified as the optimal dividend strategy,solving our optimal impulse control problem.展开更多
In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in...In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in the periodic domain.展开更多
In this article, certain Marcinkiewicz integral operators associated to surfaces of revolution on product domains were studied. The Lp boundedness for these operators are established under some rather weak conditions ...In this article, certain Marcinkiewicz integral operators associated to surfaces of revolution on product domains were studied. The Lp boundedness for these operators are established under some rather weak conditions on kernels. The main results essentially improve and extend some known results.展开更多
This paper is concerned with certain multilinear commutators of BMO functions and multilinear singular integral operators with non-smooth kernels. By the sharp maximal functions estimates, the weighted norm inequaliti...This paper is concerned with certain multilinear commutators of BMO functions and multilinear singular integral operators with non-smooth kernels. By the sharp maximal functions estimates, the weighted norm inequalities for this kind of commutators are established.展开更多
This article is concerned with the partial regularity for the weak solutions of stationary Navier-Stokes system under the controllable growth condition.By A-harmonic approximation technique,the optimal regularity is o...This article is concerned with the partial regularity for the weak solutions of stationary Navier-Stokes system under the controllable growth condition.By A-harmonic approximation technique,the optimal regularity is obtained.展开更多
Let （M1, F1） and （M2, F2） be two strongly pseudoconvex complex Finsler man- ifolds. The doubly wraped product complex Finsler manifold （f2 M1 x h M2, F） of （M1, F1） and （M2, F2） is the product manifold M1 x ...Let （M1, F1） and （M2, F2） be two strongly pseudoconvex complex Finsler man- ifolds. The doubly wraped product complex Finsler manifold （f2 M1 x h M2, F） of （M1, F1） and （M2, F2） is the product manifold M1 x M2 endowed with the warped product complex 2 2 Finsler metric F2 = f2F1 ＋ fl F2, where fl and f2 are positive smooth functions on M1 and M2, respectively. In this paper, the most often used complex Finsler connections, holomorphic curvature, Ricci scalar curvature, and real geodesics of the DWP-complex Finsler manifold are derived in terms of the corresponding objects of its components. Necessary and sufficient conditions for the DWP-complex Finsler manifold to be K/ihler Finsler （resp., weakly K/ihler Finsler, complex Berwald, weakly complex Berwald, complex Landsberg） manifold are ob- tained, respectively. It is proved that if （M1, F1） and （M2,F2） are projectively flat, then the DWP-complex Finsler manifold is projectively flat if and only if fl and f2 are positive constants.展开更多
In this paper, the authors establish the LV-mapping properties for a class of singular integrals along surfaces in Rn of the form {Ф（lul）u＇ ： u ε ]t^n} as well as the related maimal operators provided that the f...In this paper, the authors establish the LV-mapping properties for a class of singular integrals along surfaces in Rn of the form {Ф（lul）u＇ ： u ε ]t^n} as well as the related maimal operators provided that the function Ф satisfies certain oscillatory integral estimates of Van der Corput type, and the integral kernels are given by the radial function h E ε△γ（R＋） for γ 〉 1 and the sphere function ΩεFβ（S^n-1） for someβ 〉 0 which is distinct from HI（Sn-1）.展开更多
In this paper, firstly using different method and technique we derive the cor-responding integral representation formulas of （0, q）（q 〉 0） differential forms for the twotypes of the bounded domains in complex sub...In this paper, firstly using different method and technique we derive the cor-responding integral representation formulas of （0, q）（q 〉 0） differential forms for the twotypes of the bounded domains in complex submanifolds with codimension-m. Secondly weobtain the unified integral representation formulas of （0, q）（q 〉 0） differential forms for thegeneral bounded domain in complex submanifold with codimension-m, which include Hatzi-afratis formula, i.e. Koppelman type integral formula for the bounded domain with smoothboundary in analytic varieties. In particular, when m -- 0, we obtain the unified integralrepresentation formulas of （0, q）（q 〉 0） differential forms for general bounded domain in Cn,which are the generalization and the embodiment of Koppelman-Leray formula.展开更多
In this article, we mainly study the local equation of energy for weak solutions of 3D MHD equations. We define a dissipation term D（u, B） that steins from an eventual lack of smoothness in the solution, and then ob...In this article, we mainly study the local equation of energy for weak solutions of 3D MHD equations. We define a dissipation term D（u, B） that steins from an eventual lack of smoothness in the solution, and then obtain a local equation of energy for weak solutions of 3D MHD equations. Finally, we consider the 2D case at the end of this article.展开更多
In this paper, we consider the following integral system： u（x） = R n v q （y） ｜ x y ｜ nα dy, v（x） = R n u p （y） ｜ x y ｜ nμ dy, （0.1） where 0 〈 α, μ 〈 n; p, q ≥ 1. Using the method of moving planes...In this paper, we consider the following integral system： u（x） = R n v q （y） ｜ x y ｜ nα dy, v（x） = R n u p （y） ｜ x y ｜ nμ dy, （0.1） where 0 〈 α, μ 〈 n; p, q ≥ 1. Using the method of moving planes in an integral form which was recently introduced by Chen, Li, and Ou in [2, 4, 8], we show that all positive solutions of （0.1） are radially symmetric and decreasing with respect to some point under some general conditions of integrability. The results essentially improve and extend previously known results [4, 8].展开更多
This work consider boundary integrability of the weak solutions of a non-Newtonian compressible fluids in a bounded domain in dimension three, which has the constitutive equartions as ■The existence result of weak so...This work consider boundary integrability of the weak solutions of a non-Newtonian compressible fluids in a bounded domain in dimension three, which has the constitutive equartions as ■The existence result of weak solutions can be get based on Galerkin approximation. With the linear operator B constructed by BOGOVSKII, we show that the density ■is square integrable up to the boundary.展开更多
An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-...An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.展开更多
This paper is devoted to studying the commutators of the multilinear singular integral operators with the non-smooth kernels and the weighted Lipschitz functions. Some mapping properties for two types of commutators o...This paper is devoted to studying the commutators of the multilinear singular integral operators with the non-smooth kernels and the weighted Lipschitz functions. Some mapping properties for two types of commutators on the weighted Lebesgue spaces, which extend and generalize some previous results, are obtained.展开更多
The atom-bond connectivity（ABC） index provides a good model for the stability of linear and branched alkanes as well as the strain energy of cycloalkanes,which is defined as ABC（G） =∑ uv∈E（G） √d u＋dv-2 dudv,...The atom-bond connectivity（ABC） index provides a good model for the stability of linear and branched alkanes as well as the strain energy of cycloalkanes,which is defined as ABC（G） =∑ uv∈E（G） √d u＋dv-2 dudv,where du denotes the degree of a vertex u in G.A chemical graph is a graph in which no vertex has degree greater than 4.In this paper,we obtain the sharp upper and lower bounds on ABC index of chemical bicyclic graphs.展开更多
In this paper, we study the extension of isometries between the unit spheres of complex Banach spaces lp（Γ） and lp（△）（p 〉1）. We first derive the representation of isometries between the unit spheres of comple...In this paper, we study the extension of isometries between the unit spheres of complex Banach spaces lp（Γ） and lp（△）（p 〉1）. We first derive the representation of isometries between the unit spheres of complex Banach spaces lp（Γ） and lp（△）. Then we arrive at a conclusion that any surjective isometry between the unit spheres of complex Banach spaces lp（Γ）and lp（△） can be extended to be a linear isometry on the whole space.展开更多
In this paper, we consider the heat flow for the Hsystem with constant mean curvature in higher dimensions. We give sufficient conditions on the initial data such that the heat flow develops finite time singularity. W...In this paper, we consider the heat flow for the Hsystem with constant mean curvature in higher dimensions. We give sufficient conditions on the initial data such that the heat flow develops finite time singularity. We also provide a new set of initial data to guarantee the existence of global regular solution to the heat flow, that converges to zero in W 1,n with the decay rate t 2/（2-n） as time goes to infinity.展开更多
In this article, we are concerned with the strong solutions of the coupled Navier-Stokes-Poisson equations for isentropic compressible fluids in a domain Ω R^3. We prove the local existence of unique strong solution...In this article, we are concerned with the strong solutions of the coupled Navier-Stokes-Poisson equations for isentropic compressible fluids in a domain Ω R^3. We prove the local existence of unique strong solutions provided that the initial data u0 and u0 satisfy a nature compatibility condition. The important point in this article is that we allow the initial vacuum： the initial density may vanish in an open subset of Ω. This is achieved by getting some uniform estimates and using a Schauder fixed point theorem.展开更多
基金supported by National Natural Science Foundation of China(12071391,12231016)the Guangdong Basic and Applied Basic Research Foundation(2022A1515010860)。
文摘This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).
基金the NNSF of China(12071391)the Guangdong Basic and Applied Basic Research Foundation (2022A1515010069)。
文摘We study a spatiotemporal EIT problem with a dynamical boundary condition for the fractional Dirichlet-to-Neumann operator with a critical exponent.There are three major ingredients in this paper.The first is the finite time blowup and the decay estimate of the global solution with a lower-energy initial value.The second ingredient is the L^(q)(2 ≤q <∞) estimate of the global solution applying the Moser iteration,which allows us to show that any global solution is a classical solution.The third,which is the main ingredient of this paper,explores the long time asymptotic behavior of global solutions close to the stationary solution and the bubbling phenomenons by means of a concentration compactness principle.
基金supported by National Natural Science Foundation of China(12271277)the Open Research Fund of Key Laboratory of Nonlinear Analysis&Applications(Central China Normal University),Ministry of Education,China.
文摘In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.
基金the financial support from the National Natural Science Foundation of China(12171405 and 11661074)the Program for New Century Excellent Talents in Fujian Province University+2 种基金the financial support from the Characteristic&Preponderant Discipline of Key Construction Universities in Zhejiang Province(Zhejiang Gongshang University-Statistics)Collaborative Innovation Center of Statistical Data Engineering Technology&ApplicationDigital+Discipline Construction Project(SZJ2022B004)。
文摘Motivated by recent advances made in the study of dividend control and risk management problems involving the U.S.bankruptcy code,in this paper we follow[44]to revisit the De Finetti dividend control problem under the reorganization process and the regulator's intervention documented in U.S.Chapter 11 bankruptcy.We do this by further accommodating the fixed transaction costs on dividends to imitate the real-world procedure of dividend payments.Incorporating the fixed transaction costs transforms the targeting optimal dividend problem into an impulse control problem rather than a singular control problem,and hence computations and proofs that are distinct from[44]are needed.To account for the financial stress that is due to the more subtle concept of Chapter 11 bankruptcy,the surplus process after dividends is driven by a piece-wise spectrally negative Lévy process with endogenous regime switching.Some explicit expressions of the expected net present values under a double barrier dividend strategy,new to the literature,are established in terms of scale functions.With the help of these expressions,we are able to characterize the optimal strategy among the set of admissible double barrier dividend strategies.When the tail of the Lévy measure is log-convex,this optimal double barrier dividend strategy is then verified as the optimal dividend strategy,solving our optimal impulse control problem.
基金support by the NSFC(12071391,12231016)the Guangdong Basic and Applied Basic Research Foundation(2022A1515010860)support by the China Postdoctoral Science Foundation(2023M742401)。
文摘In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in the periodic domain.
基金Supported by the NSF of China (G10571122) the NFS of Fujian Province of China (Z0511004)
文摘In this article, certain Marcinkiewicz integral operators associated to surfaces of revolution on product domains were studied. The Lp boundedness for these operators are established under some rather weak conditions on kernels. The main results essentially improve and extend some known results.
基金Supported by the National Natural Science Foundation of China (10771054, 10771221, 11071200)the Youth Foundation of Wuyi University (No. xq0930)
文摘This paper is concerned with certain multilinear commutators of BMO functions and multilinear singular integral operators with non-smooth kernels. By the sharp maximal functions estimates, the weighted norm inequalities for this kind of commutators are established.
基金NNSF of China(10531020)the Program of 985 Innovation Engineering on Information in Xiamen University(2004-2007)
文摘This article is concerned with the partial regularity for the weak solutions of stationary Navier-Stokes system under the controllable growth condition.By A-harmonic approximation technique,the optimal regularity is obtained.
基金supported by Program for New Century Excellent Talents in University(NCET-13-0510)National Natural Science Foundation of China(11271304,11571288,11461064)+1 种基金the Fujian Province Natural Science Funds for Distinguished Young Scholar(2013J06001)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry
文摘Let （M1, F1） and （M2, F2） be two strongly pseudoconvex complex Finsler man- ifolds. The doubly wraped product complex Finsler manifold （f2 M1 x h M2, F） of （M1, F1） and （M2, F2） is the product manifold M1 x M2 endowed with the warped product complex 2 2 Finsler metric F2 = f2F1 ＋ fl F2, where fl and f2 are positive smooth functions on M1 and M2, respectively. In this paper, the most often used complex Finsler connections, holomorphic curvature, Ricci scalar curvature, and real geodesics of the DWP-complex Finsler manifold are derived in terms of the corresponding objects of its components. Necessary and sufficient conditions for the DWP-complex Finsler manifold to be K/ihler Finsler （resp., weakly K/ihler Finsler, complex Berwald, weakly complex Berwald, complex Landsberg） manifold are ob- tained, respectively. It is proved that if （M1, F1） and （M2,F2） are projectively flat, then the DWP-complex Finsler manifold is projectively flat if and only if fl and f2 are positive constants.
基金Supported by the National Natural Science Foundation of China(11071200,11371295)
文摘In this paper, the authors establish the LV-mapping properties for a class of singular integrals along surfaces in Rn of the form {Ф（lul）u＇ ： u ε ]t^n} as well as the related maimal operators provided that the function Ф satisfies certain oscillatory integral estimates of Van der Corput type, and the integral kernels are given by the radial function h E ε△γ（R＋） for γ 〉 1 and the sphere function ΩεFβ（S^n-1） for someβ 〉 0 which is distinct from HI（Sn-1）.
文摘In this paper, firstly using different method and technique we derive the cor-responding integral representation formulas of （0, q）（q 〉 0） differential forms for the twotypes of the bounded domains in complex submanifolds with codimension-m. Secondly weobtain the unified integral representation formulas of （0, q）（q 〉 0） differential forms for thegeneral bounded domain in complex submanifold with codimension-m, which include Hatzi-afratis formula, i.e. Koppelman type integral formula for the bounded domain with smoothboundary in analytic varieties. In particular, when m -- 0, we obtain the unified integralrepresentation formulas of （0, q）（q 〉 0） differential forms for general bounded domain in Cn,which are the generalization and the embodiment of Koppelman-Leray formula.
基金Supported by NSFC (10976026)supported by the Fundamental Research Funds for the Central Universities (11QZR18)the Research Funds for high-level talents of Huaqiao University (12BS232)
文摘In this article, we mainly study the local equation of energy for weak solutions of 3D MHD equations. We define a dissipation term D（u, B） that steins from an eventual lack of smoothness in the solution, and then obtain a local equation of energy for weak solutions of 3D MHD equations. Finally, we consider the 2D case at the end of this article.
基金Supported by National Natural Science Foundation of China-NSAF (10976026)
文摘In this paper, we consider the following integral system： u（x） = R n v q （y） ｜ x y ｜ nα dy, v（x） = R n u p （y） ｜ x y ｜ nμ dy, （0.1） where 0 〈 α, μ 〈 n; p, q ≥ 1. Using the method of moving planes in an integral form which was recently introduced by Chen, Li, and Ou in [2, 4, 8], we show that all positive solutions of （0.1） are radially symmetric and decreasing with respect to some point under some general conditions of integrability. The results essentially improve and extend previously known results [4, 8].
基金supported by the National Natural Science Foundation of China(11271305,11531010)
文摘This work consider boundary integrability of the weak solutions of a non-Newtonian compressible fluids in a bounded domain in dimension three, which has the constitutive equartions as ■The existence result of weak solutions can be get based on Galerkin approximation. With the linear operator B constructed by BOGOVSKII, we show that the density ■is square integrable up to the boundary.
文摘An FDTD system associated with uniaxial perfectly matched layer（UPML） for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic（TM） mode to Maxwell＇s equations is obtained by Yee＇s algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.
基金Supported by the National Natural Science Foundation of China (10771054,11071200)the NFS of Fujian Province of China (No. 2010J01013)
文摘This paper is devoted to studying the commutators of the multilinear singular integral operators with the non-smooth kernels and the weighted Lipschitz functions. Some mapping properties for two types of commutators on the weighted Lebesgue spaces, which extend and generalize some previous results, are obtained.
基金Supported by the National Natural Science Foundation of China(11071272,10831001,11171279,11101087)the Young Talent Foundation of Fuzhou University(XRC-1154)
文摘The atom-bond connectivity（ABC） index provides a good model for the stability of linear and branched alkanes as well as the strain energy of cycloalkanes,which is defined as ABC（G） =∑ uv∈E（G） √d u＋dv-2 dudv,where du denotes the degree of a vertex u in G.A chemical graph is a graph in which no vertex has degree greater than 4.In this paper,we obtain the sharp upper and lower bounds on ABC index of chemical bicyclic graphs.
基金supported by Higher Educational Science and Technology Program Foundation of Shandong Province(J11LA02)Young and Middle-Aged Scientists Research Foundation of Shandong Province(BS2010SF004)Higher Educational Science and Technology Program Foundation of Shandong Province(J10LA53)
文摘In this paper, we study the extension of isometries between the unit spheres of complex Banach spaces lp（Γ） and lp（△）（p 〉1）. We first derive the representation of isometries between the unit spheres of complex Banach spaces lp（Γ） and lp（△）. Then we arrive at a conclusion that any surjective isometry between the unit spheres of complex Banach spaces lp（Γ）and lp（△） can be extended to be a linear isometry on the whole space.
文摘In this paper, we consider the heat flow for the Hsystem with constant mean curvature in higher dimensions. We give sufficient conditions on the initial data such that the heat flow develops finite time singularity. We also provide a new set of initial data to guarantee the existence of global regular solution to the heat flow, that converges to zero in W 1,n with the decay rate t 2/（2-n） as time goes to infinity.
基金Supported by National Natural Science Foundation of China-NSAF (10976026)
文摘In this article, we are concerned with the strong solutions of the coupled Navier-Stokes-Poisson equations for isentropic compressible fluids in a domain Ω R^3. We prove the local existence of unique strong solutions provided that the initial data u0 and u0 satisfy a nature compatibility condition. The important point in this article is that we allow the initial vacuum： the initial density may vanish in an open subset of Ω. This is achieved by getting some uniform estimates and using a Schauder fixed point theorem.
