The small Deborah number limit for the compressible fluid-particle flows

Many natural phenomena are described by the fluid-particle two-phase flow models, including the droplets of spray, diesel engines, sedimentation analysis, biotechnology, medicine, and mineral processes [3, 4, 40], where one phase is considered as a suspension of particles within the other phase thought as a fluid. One of the models that has been widely considered is the Vlasov-Fokker-Planck equation coupled with the compressible Navier-Stokes equations (VFP-CNS):

where $f(t,x,v)$ is the density distribution function of the particles, $\rho(t,x)$, $u(t,x)$ denote the density and velocity of the fluid at time $t\geq 0$, position $x=(x_{1},x_{2},x_{3})\in\mathbb{R}^{3}$, velocity $v=(v_{1},v_{2},v_{3})\in\mathbb{R}^{3}$ respectively. $\tilde{P}(\rho)=\tilde{A}\rho^{\gamma}$ denotes the pressure function with $\tilde{A}>0,\,\gamma>1$, $\tilde{L}u=-\tilde{\mu}\Delta_{x}u-(\tilde{\mu}+\tilde{\lambda})\nabla_{x}\textup{div}_{x}u$ is the so-called the Lam$\acute{\text{e}}$ operator with constants $\tilde{\mu}+\tilde{\lambda}>0$ and $\tilde{\mu}>0$ being the dynamic viscosity of the fluid. $R$, $\mathcal{M}$ and $\rho_{P}$ are the radius, mass and associated mass density of one single spherical particle, respectively. $\rho_{F}$ is the mass density of the fluid, $k$ is the Boltzmann constant, and $\mathcal{T}>0$ is the temperature of the suspension that is assumed to be constant throughout this paper as in [22].

After re-normalization as in [10], we can deduce the scaled VFP-CNS system in the dimensionless form by introducing $\chi=\frac{\mathscr{P}}{\rho_{F}U^{2}}$, the Mach number Ma, Reynolds number Re, and Deborah number De:

where $\mathscr{P}$ is a pressure unit, $U$ is a macroscopic velocity unit, $P(\rho)=A\rho^{\gamma}$ with $A>0$ and $Lu=\frac{2}{9\tilde{\mu}}\tilde{L}u:=-\mu\Delta_{x}u-(\mu+\lambda)\nabla_{x}\textup{div}_{x}u$ with dimensionless constant $\mu>0,\mu+\lambda>0$. We refer to [27, 36] for more physical background about the dimensionless analysis as well as the Mach number Ma, Reynolds number Re, and Deborah number De.

In this paper, the small Deborah number limit in the light particles regime is considered for the scaled VFP-CNS system (1.2) in the sense that

Then the scaled VFP-CNS system (1.2) can be rewritten as

where $M:=M(v)$ is a global normalized Maxwellian equilibrium given by:

Inspired by [22], taking $\varepsilon\to 0$, the Navier-Stokes-Smoluchowski (NSS) equation is formally deduced (see Section 2) for $\big(n_{0}(t,x),u_{0}(t,x),\rho_{0}(t,x)\big)$:

Due to its physical significance, the fluid-particle two-phase model has attracted people’s attention for a long time, which can be traced back to the work of Caflish-Papanicolaou in [7]. In the seminal work [22, 23], Goudon-Jabin-Vasseur first established a systematic framework to study the particle-fluid two-phase flow, i.e., the VFP equation coupled with incompressible Navier-Stokes equation, where their study involved two distinct scalings for the light particles (density of the particle is much less than the fluid density, i.e., $\rho_{P}/\rho_{F}=O(\varepsilon^{2})$) and fine particles (density of the particle is of the same order as gas density, i.e., $\rho_{P}\approx\rho_{F}$), as informed by dimensional analysis; our model (1.4) follows the scaling proposed in [22], though the compressible Navier-Stokes equation is considered to be coupled with the VFP equation.

To clearly state our contributions in this paper, we begin with a comprehensive overview of the fluid-particle model, encompassing the previous results of both its well-posedness and hydrodynamic limit; furthermore, the mathematical challenges and contributions of our work are illustrated for comparison.

Previous results for “well-posedness” of the fluid-particle model: There are many results concerning the well-posedness of the VFP equation coupled with the compressible Navier-Stokes system. Mellet-Vasseur studied the existence of global weak solutions to the VFP-CNS system with Dirichlet or reflection boundary conditions in [32]. Li-Mu-Wang in [28] obtained the global well-posedness of a strong solution when the initial data is a small perturbation of some given equilibrium, along with the algebraic convergence rate of a solution toward the equilibrium, and similar results can be extended to the VFP equation coupled with the non-isentropic CNS equation by Mu-Wang in [34]. Recently, the existence of a global-in-time strong solution to the VFP-CNS system with specular reflection boundary conditions was proved by Li-Liu-Yang in [29]. Besides, in the presence of various boundary conditions, the global existence of weak solutions to the incompressible case has been shown in [8, 29, 39, 41], and we refer the readers to our previous work [21] and the references therein for more results regarding the VFP-INS model. As for VFP equation coupled with the Euler equations, Carrillo-Duan-Moussa [9] and Duan-Liu [16] proved the existence and large time behaviors of the classical solutions to the Cauchy problem, while the stability and asymptotic analysis of such a model was studied by Carrillo-Goudon in [10].

Previous results for “hydrodynamic limit” and “limiting two-phase model”: Another important research field of the fluid-particle model is its hydrodynamic limit. Carrillo-Goudon in [10] first studied the VFP equation coupled with the compressible Euler equations, where they formally derived the hydrodynamic model in the so-called “bubbling” and “flowing” regime with different scalings. The rigorous justifications for hydrodynamic limit in these two regimes were completed by Mellet-Vasseur [33] and Ballew [2], respectively. It is worth mentioning that both results relied on the relative entropy method (also called the “modulated energy” method for different asymptotic problems), which is reminiscent of the weak-strong uniqueness principle by Dafermos [5] and Lions [30]. Recently, Choi-Jung also applied a similar strategy to study the case of VFP equation coupled with compressible Navier-Stokes equations with a density-dependent viscosity [12]. The hydrodynamic limit results of the VFP equation coupled with the incompressible Navier-Stokes equations were established by Goudon-Jabin-Vasseur in [22, 23] via the weak compactness and relative entropy argument. In our previous work [21], we further rigorously justified the hydrodynamic limit with an explicit convergent rate by designing a new expansion form and applying the refined energy estimate.

Furthermore, to rigorously justify the hydrodynamic limit in a stronger topology (e.g., pointwise), it is essential to rely on the well-posedness and regularity theory of the limiting system. In fact, the limiting two-phase system in the so-called “bubbling” regime [10, Section 4.3] corresponds to the NSS system (1.6). Within this framework, Huang-Ding-Wen established the local existence and uniqueness of the strong solutions to the compressible NSS system in [24], which was later extended to global-in-time solutions for initial data with small energy near equilibrium in [14]. In the presence of a “small” external potential, the global existence of classical solutions was further proved by Ding-Huang-Li in [13]. For additional results on the NSS system, we refer the reader to [1, 11, 20]. Although this paper primarily focuses on the “bubbling” regime, another class of two-phase limiting models, arising from the “flowing” regime [10, Section 4.2], has attracted significant attention in recent years. In [37], Vasseur-Wen-Yu established the global existence of weak solutions for a bi-fluid two-phase model with a pressure law. This result was further generalized by Novotny-Pokorny in [35] and Wen in [38]. In [6], Bresch-Mucha-Zatorska proved the global existence of weak solutions to a related two-fluid compressible Stokes system for more general adiabatic index. For further developments on the limiting two-phase system in the “flowing” regime, we refer the reader to [18, 19, 26]. These results mentioned above motivate our future goal of rigorously justifying the hydrodynamic limit in the “flowing” regime.

Mathematical challenges and our contributions: As we discussed above, in the previous results concerning the hydrodynamic limit of the VFP-CNS system [2, 33], the proof essentially relied on the relative entropy argument. Albeit successful, the pointwise convergence seems to be tough to obtain through the “modulated energy” method. As a consequence, our main purpose in this paper is to rigorously justify the hydrodynamic limit from the VFP-CNS system (1.4) to the NSS system (1.6) in a pointwise sense with an explicit convergence rate. Specifically speaking, we start with looking for a special class of solutions to the scaled kinetic-fluid coupled system (1.4) in the form of Hilbert expansion, where, in contrast with the previous scalings corresponding to the flow regime [33] and bubble regime [2], we take the scaling in (1.4) from [22] to model the two-phase flow of light particles. The key point here is that the Hilbert expansion is taken around the classical solution to our limiting macroscopic system (1.6), hence, the essential part of the proof lies in the establishment of the uniform energy estimate for the remainder system (3.2), which is manageable thanks to its less singular and nonlinear property than the original coupled system (1.4). Note that this sort of strategy and expansion has been applied in justifying the stability of the Boltzmann equation near Couette flows in [17], the hydrodynamic limit of the self-organized kinetic equation coupled with a fluid equation in [25], as well as the VFP equation coupled with the incompressible Navier-Stokes equations in our preceding work [21]. However, even though the decomposition framework has proven highly effective for incompressible flows, extending this methodology to the compressible regime is highly non-trivial. The critical challenges, along with the novelties of our strategies, are illustrated below:

• •

  In contrast with the incompressible case, where the density is a constant, the appearances of the extra continuity equation and the density-dependent pressure function in the coupled compressible Navier-Stokes equation bring more difficulties in our analysis. The notable challenge arises from the strong couplings between the momentum equations and the continuity equation, leading to severe singularity in the remainder system such that the usual expansion form of $(f^{\varepsilon},u^{\varepsilon},\rho^{\varepsilon})$, directly derived from the classical Hilbert expansion, would no longer be sufficient to directly handle the singular terms via dissipation. Hence, we construct the expansion (1.14) based on our observation and formal analysis (see Section 2). Note that the primary distinction of the expansion in the compressible case lies in the extra density expansion of $\rho^{\varepsilon}=\rho_{0}+\varepsilon\rho=(1+h_{0})+\varepsilon\rho$. Still, balancing the nonlinear pressure terms gives rise to a distinctive singular term $\frac{1}{\varepsilon}R_{2}$ in (3.3) in the compressible case,

  To address this, we introduce a new auxiliary function (e.g., $B^{\gamma}_{\alpha}$ in Lemma 3.1)

  $$B^{\alpha}_{\gamma}=\partial_{x}^{\alpha}\left[(1+h_{0}+\varepsilon\rho)^{\gamma}\right]-\partial_{x}^{\alpha}\left[(1+h_{0})^{\gamma}\right]-\gamma\varepsilon(1+h_{0}+\varepsilon\rho)^{{\gamma}-1}\partial_{x}^{\alpha}\rho,$$

  and explicitly calculate its spatial derivatives before applying Taylor’s expansion to each resulting component, where an $\mathcal{O}(\varepsilon)$ term naturally emerges and exactly cancels out the singularity of $\frac{1}{\varepsilon}$ inherent in the pressure term. This crucial cancellation enables us to uniformly control the nonlinearity from the pressure term and successfully close the energy estimates. This newly-designed expansion profile provides the “optimal” order of convergence in the sense of aligning with the order of singularity in the remainder system (3.2). It is worth mentioning that the applications of such a novel strategy, i.e., re-designing the expansion form to match the singular behavior and close the energy estimate, are expected to be extensively fruitful, especially when the high-order singularity exists.
  

• •

  In addition, our study extends the previous convergence results by the relative entropy method [2, 33] to a pointwise convergence (see Corollary 1.1), thanks to the high-regularity estimates. To this end, we apply dedicated energy estimates to the reminder system (3.2) by proposing the refined energy/dissipation structures that incorporate higher regularity in both spatial and velocity variables (see Section 3.2). More specifically, our refined energy/dissipation functionals (3.9) are designed based on the well-established macro-micro decomposition [15] such that they consist of contributions from both macroscopic and microscopic parts, which also perfectly match our new expansion form (see (1.14)). The refined total energy estimate of the reminder system (3.2) can be found in Proposition 3.2.

The notations that will be used throughout this paper are introduced as follows:

(i) $A_{1}\lesssim B_{1}$ stands for $A_{1}\leq CB_{1}$ with generic constant $C>0$. $A_{1}\sim B_{1}$ stands for $C_{1}A_{1}\leq B_{1}\leq C_{2}A_{1}$ with some generic constants $C_{1},C_{2}>0$.

(ii) For multi-indices $\alpha=(\alpha_{1},\alpha_{2},\alpha_{3})$ and $\beta=(\beta_{1},\beta_{2},\beta_{3})$, we denote

(iii) For $d,e\in\mathbb{N}$, we denote the following inner-product and the spaces:

Let $\nu(v)=1+|v|^{2}$ and denote $\|\cdot\|_{\nu}$ by

(iv) By following [15], the velocity orthogonal projection $\mathbf{P}\textup{:}\,L^{2}_{v}\to\textup{Span}\{\sqrt{M},v_{1}\sqrt{M},v_{2}\sqrt{M},v_{3}\sqrt{M}\}$ is denoted as

with

The linearized Fokker-Planck operator $\mathcal{L}$ is denoted as

with $\textup{Ker}\mathcal{L}=\textup{Span}\{\sqrt{M}\}$ in $L^{2}_{v}$ norm, and $\mathcal{L}$ can be decomposed by

Notice that $\mathbf{I}-\mathbf{P}_{0}$, $\mathbf{I}-\mathbf{P}$ is self-adjoint in $H^{d}_{x,v}$, i.e., for any $f,g\in H^{d}_{x,v}$,

and it can be verified that

According to [9], the dissipative property holds for $\mathcal{L}$, i.e., there exists a constant $c_{0}>0$ such that

Let $h_{0}(t,x):=\rho_{0}(t,x)-1,m_{0}(t,x):=n_{0}(t,x)-1$, the NSS system (1.6) becomes

For the NSS system (1.11) above, we have the following well-posedness result:

We seek a solution to the VFP-CNS system (1.4) in the following form by employing $(m_{0},u_{0},h_{0})$ of the NSS system (1.11):

where $g$, $u$, $\rho$ also depend on $\varepsilon$, but, for simplicity, we will omit the dependence in the notation throughout the paper.

Then, the main theorem is presented as follows:

By applying the Sobolev embedding $H^{2}\hookrightarrow L^{\infty}$, we obtain the following corollary:

The rest of the paper is organized as follows: in Section 2, we present the formal analysis. In Section 3, we derive global-in-time a priori estimates for the remainder system. Finally, the proofs of the main Theorem and Corollary are illustrated in Section 4.

Based on the formal analysis in Section 2, we introduce the remainder terms of the VFP-CNS system (1.4):

where

and $(n_{0},u_{0},\rho_{0})$ is the solution to the NSS system (1.6). In other words, we seek a solution to the VFP-CNS system (1.4) in our newly-designed expansion form (1.14).