alphaFeedback

WorksMathematics

Uniform Regularity for the Compressible Navier-Stokes-Poisson System under Navier-Slip Boundary Conditions

AI contributions
Not recorded for this paper
View details

Abstract

In this paper, we establish the existence of a local-in-time strong solution to the compressible Navier–Stokes–Poisson system in a three-dimensional bounded smooth domain, subject to the general Navier–slip boundary condition. The solution is shown to be uniformly bounded in conormal Sobolev spaces, with uniform regularity estimates that are independent of the viscosity coefficients. As a consequence, we obtain the inviscid limit to the corresponding compressible Euler–Poisson system via a strong compactness argument.

Introduction and main results

The compressible Navier–Stokes–Poisson (NSP) system is a classical model that couples the Navier–Stokes equations with a self-consistent Poisson equation, and is commonly employed to describe the motion of charged particles such as electrons or holes. In this paper, we study the initial-boundary value problem for the compressible Navier–Stokes–Poisson system in a three-dimensional smooth bounded domain Ω⊂R3\Omega\subset\mathbb{R}^{3}:

{ρtε+div⁡(ρεuε)=0,ρεutε+ρεuε⋅∇uε+∇pε=μεΔuε+(λ+μ)ε∇div⁡uε+ρε∇ϕε,Δϕε=ρε−1,(1)\begin{cases} \rho^{\varepsilon}_{t}+\operatorname{div}(\rho^{\varepsilon}u^{\varepsilon})=0,\\ \rho^{\varepsilon}u^{\varepsilon}_{t}+\rho^{\varepsilon}u^{\varepsilon}\cdot\nabla u^{\varepsilon}+\nabla p^{\varepsilon}=\mu\varepsilon\Delta u^{\varepsilon}+(\lambda+\mu)\varepsilon\nabla\operatorname{div}u^{\varepsilon}+\rho^{\varepsilon}\nabla\phi^{\varepsilon},\\ \Delta\phi^{\varepsilon}=\rho^{\varepsilon}-1, \tag*{(1)} \end{cases}

where ρε\rho^{\varepsilon}, uεu^{\varepsilon}, and ϕε\phi^{\varepsilon} denote the density, velocity, and electric potential, respectively. The pressure is given by

pε=(ρε)γ,γ>1,p^{\varepsilon}=(\rho^{\varepsilon})^{\gamma},\qquad\gamma>1,

and the viscosity coefficients satisfy

μ>0,2μ+3λ>0.\mu>0,\qquad2\mu+3\lambda>0.

Here 0<ε≤10<\varepsilon\leq1 is the viscosity parameter, and the background charge density in the Poisson equation is normalized to be 1.

We impose a general Navier–slip boundary condition on the velocity. In its stress form, this condition is written as

uε⋅n=0,(Suε⋅n)τ=−(Auε)τon ∂Ω,u^\varepsilon\cdot n = 0,\qquad\left(Su^\varepsilon\cdot n\right)_\tau= -\left(Au^\varepsilon\right)_\tau\quad\text{on } \partial\Omega,

where

Suε=12(∇uε+(∇uε)T).Su^\varepsilon= \frac{1}{2}\left(\nabla u^\varepsilon+ \left(\nabla u^\varepsilon\right)^{T}\right).

nn is the unit outward normal vector to ∂Ω\partial\Omega, AA is a given smooth symmetric matrix on the boundary, and (⋅)τ(\cdot)_\tau denotes the tangential component. For smooth solutions, using the identity

(2S(v)n−(∇×v)×n)τ=−(2S(n)v)τ,\left(2S(v)n - (\nabla\times v)\times n\right)_\tau= -\left(2S(n)v\right)_\tau,

this boundary condition can be equivalently rewritten in the vorticity form

uε⋅n=0,n×ωε=[Buε]τon ∂Ω,(2)u^\varepsilon\cdot n = 0,\qquad n \times\omega^\varepsilon= \left[Bu^\varepsilon\right]_\tau\quad\text{on } \partial\Omega, \tag*{(2)}

where

ωε=∇×uε,B=2(A−S(n)).\omega^\varepsilon= \nabla\times u^\varepsilon,\qquad B = 2(A-S(n)).

This vorticity formulation is more convenient for the normal derivative and boundary vorticity estimates. For the electric potential, we impose the homogeneous Neumann boundary condition

∂nϕε=0on ∂Ω.(3)\partial_n\phi^\varepsilon= 0 \quad\text{on } \partial\Omega. \tag*{(3)}

As ε→0\varepsilon\to0, the viscous terms in (1) formally vanish, and the corresponding limiting system is the compressible Euler–Poisson system

{ρt+div⁡(ρu)=0,ρut+ρu⋅∇u+∇p=ρ∇ϕ,x∈Ω,t>0,Δϕ=ρ−1,(4)\begin{cases} \rho_t + \operatorname{div}(\rho u) = 0,\\ \rho u_t + \rho u\cdot\nabla u + \nabla p = \rho\nabla\phi,\qquad x\in\Omega,\quad t>0,\\ \Delta\phi= \rho- 1, \tag*{(4)} \end{cases}

where p=ργp=\rho^\gamma. The limiting boundary conditions are

u⋅n=0,∂nϕ=0on ∂Ω.(5)u\cdot n = 0,\qquad\partial_n\phi= 0 \quad\text{on } \partial\Omega. \tag*{(5)}

Thus, our aim is to establish uniform estimates for (1)–(3) in a conormal Sobolev framework and to justify the inviscid limit from (1) to (4).

The vanishing viscosity limit for the Navier–Stokes equations is a classical problem in mathematical fluid mechanics. In the absence of boundaries, the limit can be justified in the standard Sobolev framework as long as the corresponding Euler solution remains smooth [15]. In domains with physical boundaries, however, the problem becomes much more delicate because viscous and inviscid flows usually satisfy different boundary conditions. For the no-slip boundary condition, the inviscid limit is closely related to the Prandtl boundary layer theory and remains largely open, except in some special settings such as analytic data, vorticity initially separated from the boundary, or certain steady regimes [24, 25, 21, 23]. In contrast, Navier–slip type boundary conditions are more compatible with the inviscid boundary condition and therefore make it possible to establish uniform regularity estimates. For incompressible flows, uniform regularity and the inviscid limit under slip or Navier-type boundary conditions were studied in [35, 2, 10, 36]. A decisive progress was made by Masmoudi and Rousset [22], who developed conormal Sobolev estimates in general domains and justified the vanishing viscosity limit by compactness. This conormal approach was later adapted to the compressible Navier–Stokes equations with general Navier–slip boundary conditions by Wang, Xin and Yong [34], and to the full compressible Navier–Stokes system by Wang [33]. These works provide the main analytical background for the present paper, where conormal derivatives are used to capture tangential regularity while the missing normal regularity is recovered from the equations and the boundary structure.

The main object of this paper is the compressible Navier–Stokes–Poisson system. This model describes the motion of a compressible charged fluid coupled with a self-consistent electrostatic field and appears in plasma physics, semiconductor theory, ionized gas dynamics and related hydrodynamic models of charged particles. A systematic introduction to the compressible Navier–Stokes–Poisson equations can be found in Hsiao and Li [9]. Compared with the usual compressible Navier–Stokes equations, the Poisson equation introduces a nonlocal elliptic coupling into the momentum equation, and the electric force has to be estimated together with the fluid variables. This changes both the energy structure and the regularity analysis, especially in the presence of physical boundaries.

There is a large literature on the well-posedness and large-time behavior of the Navier–Stokes–Poisson system. Weak solutions were studied, for example, in [5, 16, 37], where the main difficulty is to combine the compactness of the fluid variables with the elliptic structure of the self-consistent potential. For strong solutions in the Cauchy problem, global existence and decay around equilibrium states have been obtained in various perturbative frameworks [8, 17, 26, 32, 28, 3]. These works show that the electric field may produce additional dissipation or dispersive effects, but at the same time creates new coupling terms in the high-order energy estimates. For problems with boundaries or nontrivial background profiles, stability of steady states, global well-posedness and decay have also been investigated in bounded or exterior domains [27, 19, 20, 4]. In particular, the treatment of boundary conditions for the electric potential is essential, since different choices of boundary conditions may lead to different boundary mechanisms. There are also related models with additional physical effects, such as damping term, heat conductivity, capillarity or quantum/Korteweg-type corrections [29, 30, 14, 1]. These variants further illustrate that the Poisson coupling interacts strongly with the dissipative and dispersive structures of the fluid equations.

Singular limits for the Navier–Stokes–Poisson system form another important direction. The quasi-neutral limit, the vanishing Debye length limit, the relaxation limit and the vanishing viscosity limit have been studied by compactness methods, relative entropy methods, dispersive estimates, boundary layer analysis and energy estimates. In the whole space or periodic domains, the convergence from Navier–Stokes–Poisson type systems to incompressible Euler, Euler–Poisson or related limiting systems has been considered in [31, 6, 11, 18]. For half-space or boundary problems, the interaction between the Poisson equation and the boundary conditions becomes more subtle. Ju and Xu [12] studied quasi-neutral and zero-viscosity limits for the Navier–Stokes–Poisson equations in the half-space. Ju, Luo and Xu [13] justified singular limits for a viscous plasma under Navier–slip boundary condition for the velocity and Dirichlet boundary condition for the electric potential; in that problem, the main difficulty is the strong density boundary layer generated by the prescribed boundary potential and its interaction with the weak velocity boundary layer. In contrast, when the electric potential satisfies a homogeneous Neumann boundary condition, the electrostatic energy is naturally compatible with integration by parts, and the strong boundary layer mechanism caused by the prescribed potential is absent. Nevertheless, uniform high-order estimates remain nontrivial because one still has to control the nonlocal electric field together with the compressible fluid variables.

Recently, uniform regularity and vanishing viscosity limits for Navier–Stokes–Poisson type systems have attracted further attention. For the compressible Navier–Stokes–Poisson system with damping term under the slip boundary condition, Gao, Peng, Wei and Yao [7] established global-in-time uniform bounds and obtained the vanishing viscosity limit with convergence rates. The damping term plays a stabilizing role in their analysis and allows global-in-time estimates. The present paper is concerned with a different setting: we consider the compressible Navier–Stokes– Stokes–Poisson system without damping term in a general three-dimensional bounded domain, equipped with the Navier–slip boundary condition for the velocity and the homogeneous Neumann boundary condition for the electrostatic potential. Therefore, the problem is local in time, but it retains the essential difficulties caused by the curved boundary, the compressibility and the nonlocal Poisson coupling.

The goal of this paper is to establish uniform conormal regularity estimates for the Navier–Stokes–Poisson system, independently of the viscosity, and then to justify the vanishing viscosity limit to the corresponding Euler–Poisson system. Compared with the compressible Navier–Stokes theory [34, 33], the new difficulty is the electric force ρ∇ϕ\rho\nabla\phi, which requires elliptic estimates for the Poisson equation at the same conormal level as the fluid estimates. In particular, one has to control ∇ϕ\nabla\phi in both conormal Sobolev norms and L∞L^\infty-type norms, and to show that the Poisson coupling can be incorporated into the energy method without losing derivatives. The homogeneous Neumann boundary condition for ϕ\phi is crucial in this respect. Our analysis shows that the conormal framework remains effective for the Navier–Stokes–Poisson coupling in general bounded domains and yields, by compactness, the inviscid Euler–Poisson limit.

We now state the main results. The system and boundary conditions have been given above. Throughout this paper, for each 0<ε≤10 < \varepsilon\le1, the initial potential ϕ0\phi_0 is determined by

Δϕ0=ρ0−1,∂nϕ0=0on ∂Ω.\Delta\phi_0=\rho_0-1,\qquad\partial_n\phi_0=0\quad\text{on }\partial\Omega.

We also assume the compatibility condition

∫Ω(ρ0−1) dx=0.\int_\Omega(\rho_0-1)\,dx=0.

We define the functional space adapted to the viscous Navier–Stokes–Poisson system. For 0<ε≤10 < \varepsilon\le1, set

∥(p,u,∇ϕ)(t)∥Xmε2:=∥(p,u,∇ϕ)(t)∥Hm2+∥∇u(t)∥Hm−12+∑k=0m−2∥∂tk∇p(t)∥m−1−k2+∥Δp(t)∥H12+∥(∇u,∇ϕ)(t)∥H1,∞2+ε∥∇∂tm−1p(t)∥L22+ε∥Δp(t)∥H22.(6)\begin{aligned} \|(p,u,\nabla\phi)(t)\|_{X_m^\varepsilon}^2 &:= \|(p,u,\nabla\phi)(t)\|_{\mathcal{H}^m}^2+\|\nabla u(t)\|_{\mathcal{H}^{m-1}}^2\\ &\quad+\sum_{k=0}^{m-2}\|\partial_t^k\nabla p(t)\|_{m-1-k}^2+\|\Delta p(t)\|_{\mathcal{H}^1}^2\\ &\quad+\|(\nabla u,\nabla\phi)(t)\|_{\mathcal{H}^{1,\infty}}^2+\varepsilon\|\nabla\partial_t^{m-1}p(t)\|_{L^2}^2+\varepsilon\|\Delta p(t)\|_{\mathcal{H}^2}^2. \tag*{(6)} \end{aligned}

For the limiting Euler–Poisson system (4), we use the notation

∥(p,u,∇ϕ)(t)∥Xm02:=∥(p,u,∇ϕ)(t)∥Hm2+∥∇u(t)∥Hm−12+∑k=0m−2∥∂tk∇p(t)∥m−1−k2+∥Δp(t)∥H12+∥(∇u,∇ϕ)(t)∥H1,∞2.(7)\begin{aligned} \|(p,u,\nabla\phi)(t)\|_{X_m^0}^2 &:= \|(p,u,\nabla\phi)(t)\|_{\mathcal{H}^m}^2+\|\nabla u(t)\|_{\mathcal{H}^{m-1}}^2\\ &\quad+\sum_{k=0}^{m-2}\|\partial_t^k\nabla p(t)\|_{m-1-k}^2+\|\Delta p(t)\|_{\mathcal{H}^1}^2\\ &\quad+\|(\nabla u,\nabla\phi)(t)\|_{\mathcal{H}^{1,\infty}}^2. \tag*{(7)} \end{aligned}

Here ∥⋅∥Hm\|\cdot\|_{\mathcal{H}^m} denotes the space-time conormal Sobolev norm, while ∥⋅∥k\|\cdot\|_k denotes the pure spatial conormal Sobolev norm. Moreover, ∥⋅∥H1,∞\|\cdot\|_{\mathcal{H}^{1,\infty}} denotes the corresponding conormal L∞L^\infty-type norm. Their precise definitions are given in § 2.

Compared with the space used for the compressible Navier–Stokes equations, the space in (1.6) contains the additional electric-field term ∥∇ϕ(t)∥Hm2\|\nabla\phi(t)\|_{\mathcal{H}^m}^2, which is required to control the force ρ∇ϕ\rho\nabla\phi coupled through the Poisson equation.

To complete the setup of our functional framework, we define the uniform solution space Xmε(T)X_m^\varepsilon(T) as

Xmε(T)={(p,u,∇ϕ)∈L∞([0,T];L2(Ω)) | ess sup0≤t≤T∥(p,u,∇ϕ)(t)∥Xmε<+∞}.(8)X_m^\varepsilon(T)=\left\{(p,u,\nabla\phi)\in L^\infty([0,T];L^2(\Omega))\ \middle|\ \mathop{\mathrm{ess\,sup}}_{0\le t\le T}\|(p,u,\nabla\phi)(t)\|_{X_m^\varepsilon}<+\infty\right\}. \tag*{(8)}

We assume that the initial data satisfy

0<1C^0≤ρ0(x)≤C^0<∞,(9)0 < \frac{1}{\widehat{C}_0} \leq\rho_0(x) \leq\widehat{C}_0 < \infty, \tag*{(9)}

and

sup⁡0≤ε≤1∥(p0,u0,∇ϕ0)∥Xεm≤C~0,(10)\sup_{0\leq\varepsilon\leq1}\left\|(p_0,u_0,\nabla\phi_0)\right\|_{X_\varepsilon^m}\leq\widetilde{C}_0, \tag*{(10)}

where

p0=(ρ0)γ.p_0=(\rho_0)^\gamma.

References

  1. [1]P. Antonelli and S. Spirito, Global existence of weak solutions to the Navier–Stokes–Korteweg equations, Annales de l’Institut Henri Poincaré C, Analyse Non Linéaire, 39 (2022), 171–200.arxiv.org/abs/1903.02441
  2. [2]H. Beirão da Veiga and F. Crispo, Concerning the W^{k,p}-inviscid limit for 3-D flows under a slip boundary condition, Journal of Mathematical Fluid Mechanics, 13 (2011), 117–135.
  3. [3]Q. Y. Bie, Q. Wang and Z.-A. Yao, Optimal decay rate for the compressible Navier–Stokes–Poisson system in the critical L^p framework, Journal of Differential Equations, 263 (2017), 8391–8417.DOI
  4. [4]Y. Z. Chen, B. Huang and X. D. Shi, Global well-posedness of classical solutions to the compressible Navier–Stokes–Poisson equations with slip boundary conditions in 3D bounded domains, Journal of Mathematical Fluid Mechanics, 26 (2024), Article 41.arxiv.org/abs/2102.07938
  5. [5]D. Donatelli, Local and global existence for the coupled Navier–Stokes–Poisson problem, Quarterly of Applied Mathematics, 61 (2003), 345–361.DOI
  6. [6]D. Donatelli and P. Marcati, A quasineutral type limit for the Navier–Stokes–Poisson system with large data, Nonlinearity, 21 (2008), 135–148.DOI
  7. [7]J. C. Gao, L. Y. Peng, Y. C. Wei and Z.-A. Yao, Global uniform regularity and vanishing viscosity limit for the compressible Navier–Stokes–Poisson equations with damping under slip boundary condition, Journal of Differential Equations, 460 (2026), 114067.DOI
  8. [8]C. C. Hao and H. L. Li, Global existence for compressible Navier–Stokes–Poisson equations in three and higher dimensions, Journal of Differential Equations, 246 (2009), 4791–4812.DOI
  9. [9]L. Hsiao and H. L. Li, Compressible Navier–Stokes–Poisson equations, Acta Mathematica Scientia, 30 (2010), 1937–1948.DOI
  10. [10]D. Iftimie and F. Sueur, Viscous boundary layers for the Navier–Stokes equations with the Navier slip conditions, Archive for Rational Mechanics and Analysis, 199 (2011), 145–175.
  11. [11]Q. C. Ju, F. C. Li and H. L. Li, The quasineutral limit of compressible Navier–Stokes–Poisson system with heat conductivity and general initial data, Journal of Differential Equations, 247 (2009), 203–224.arxiv.org/abs/0905.2872
  12. [12]Q. C. Ju and X. Xu, Quasi-neutral and zero-viscosity limits of Navier–Stokes–Poisson equations in the half-space, Journal of Differential Equations, 264 (2018), 867–896.DOI
  13. [13]Q. C. Ju, T. Luo and X. Xu, Singular limits for the Navier–Stokes–Poisson equations of the viscous plasma with the strong density boundary layer, Science China Mathematics, 66 (2023), 1495–1528.
  14. [14]A. Jingel, Global weak solutions to compressible Navier–Stokes equations for quantum fluids, SIAM Journal on Mathematical Analysis, 42 (2010), 1025–1045.
  15. [15]T. Kato, Nonstationary flows of viscous and ideal fluids in R^3, Journal of Functional Analysis, 9 (1972), 296–305.DOI
  16. [16]T. Kobayashi and T. Suzuki, Weak solutions to the Navier–Stokes–Poisson equation, Advances in Mathematical Sciences and Applications, 18 (2008), 141–168.
  17. [17]H. L. Li, A. Matsumura and G. J. Zhang, Optimal decay rate of the compressible Navier–Stokes–Poisson system in ℝ^3, Archive for Rational Mechanics and Analysis, 196 (2010), 681–713.DOI
  18. [18]Y. Li, Q. C. Ju and W. Q. Xu, Quasi-neutral limit of the full Navier–Stokes–Fourier–Poisson system, Journal of Differential Equations, 258 (2015), 3661–3687.DOI
  19. [19]H. R. Liu, T. Luo and H. Zhong, Global solutions to compressible Navier–Stokes–Poisson and Euler–Poisson equations of plasma on exterior domains, Journal of Differential Equations, 269 (2020), 9936–10001.DOI
  20. [20]H. R. Liu and H. Zhong, Global solutions to the initial boundary problem of 3-D compressible Navier–Stokes–Poisson equations on bounded domains, Zeitschrift für angewandte Mathematik und Physik, 72 (2021), Article 78.DOI
  21. [21]Y. Maekawa, On the inviscid limit problem of the vorticity equations for viscous incompressible flows in the half-plane, Communications on Pure and Applied Mathematics, 67 (2014), 1045–1128.DOI
  22. [22]N. Masmoudi and F. Rousset, Uniform regularity for the Navier–Stokes equation with Navier boundary condition, Archive for Rational Mechanics and Analysis, 203 (2012), 529–575.
  23. [23]T. T. Nguyen and T. T. Nguyen, The inviscid limit of Navier–Stokes equations for analytic data on the half-space, Archive for Rational Mechanics and Analysis, 230 (2018), 1103–1129.DOI
  24. [24]M. Sammartino and R. E. Caflisch, Zero viscosity limit for analytic solutions of the Navier–Stokes equation on a half-space. I. Existence for Euler and Prandtl equations, Communications in Mathematical Physics, 192 (1998), 433–461.DOI
  25. [25]M. Sammartino and R. E. Caflisch, Zero viscosity limit for analytic solutions of the Navier–Stokes equation on a half-space. II. Construction of the Navier–Stokes solution, Communications in Mathematical Physics, 192 (1998), 463–491.DOI
  26. [26]Z. Tan and G. C. Wu, Global existence for the non-isentropic compressible Navier–Stokes–Poisson system in three and higher dimensions, Nonlinear Analysis: Real World Applications, 13 (2012), 650–664.DOI
  27. [27]Z. Tan, Y. Wang and Y. Wang, Stability of steady states of the Navier–Stokes–Poisson equations with nonflat doping profile, SIAM Journal on Mathematical Analysis, 47 (2015), 179–209.arxiv.org/abs/1506.02199
  28. [28]Z. Tan, T. Yang, H. J. Zhao and Q. Y. Zou, Global solutions to the one-dimensional compressible Navier–Stokes–Poisson equations with large data, SIAM Journal on Mathematical Analysis, 45 (2013), 547–571.DOI
  29. [29]D. Wang and G.-Q. Chen, Formation of singularities in compressible Euler–Poisson fluids with heat diffusion and damping relaxation, Journal of Differential Equations, 144 (1998), 44–65.DOI
  30. [30]D. Wang, Global solutions and relaxation limits of Euler–Poisson equations, Zeitschrift für angewandte Mathematik und Physik, 52 (2001), 620–630.
  31. [31]S. Wang and S. Jiang, The convergence of the Navier–Stokes–Poisson system to the incompressible Euler equations, Communications in Partial Differential Equations, 31 (2006), 571–591.
  32. [32]Y. Wang, Decay of the Navier–Stokes–Poisson equations, Journal of Differential Equations, 253 (2012), 273–297.arxiv.org/abs/1112.4902
  33. [33]Y. Wang, Uniform regularity and vanishing dissipation limit for the full compressible Navier–Stokes system in three-dimensional bounded domain, Archive for Rational Mechanics and Analysis, 221 (2016), 1345–1415.DOI
  34. [34]Y. Wang, Z. P. Xin and Y. Yong, Uniform regularity and vanishing viscosity limit for the compressible Navier–Stokes with general Navier–slip boundary conditions in three-dimensional domains, SIAM Journal on Mathematical Analysis, 47 (2015), 4123–4191.DOI
  35. [35]Y. L. Xiao and Z. P. Xin, On the vanishing viscosity limit for the 3D Navier–Stokes equations with a slip boundary condition, Communications on Pure and Applied Mathematics, 60 (2007), 1027–1055.DOI
  36. [36]Y. L. Xiao and Z. P. Xin, On the inviscid limit of the 3D Navier–Stokes equations with generalized Navier–slip boundary conditions, Communications in Mathematics and Statistics, 1 (2013), 259–279.arxiv.org/abs/1301.0693
  37. [37]Y. H. Zhang and Z. Tan, On the existence of solutions to the Navier–Stokes–Poisson equations of a two-dimensional compressible flow, Mathematical Methods in the Applied Sciences, 30 (2007), 305–329.DOI

Paper details

vlnc:W_68dhcdnd
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents