For the isentropic compressible fluids in one-space dimension, we prove that the Navier-Stokes equations with density-dependent viscosity have neither forward nor backward self-similar strong solutions with finite kinetic energy. Moreover, we obtain the same result for the nonisentropic compressible gas flow, that is, for the fluid dynamics of the Navier-Stokes equations coupled with a transport
equation of entropy. These results generalize those in Guo and Jiang's work (2006) where the one-dimensional compressible fluids with constant viscosity are considered.
1. Introduction
Self-similar solutions have attracted much attention in mathematical physics because understanding them is fundamental and important for investigating the well-posedness, regularity, and asymptotic behavior of differential equations in physics. Since the pioneering work of Leray [1], self-similar solutions of the Navier-Stokes equations for incompressible fluids have been widely studied in different settings (e.g., [2, page 207]; [3, page 120]; [4–10]; [11, Chapter 23]; [12–20]). On the contrary, studies on the self-similar solutions of the compressible Navier-Stokes equations have been limited partially due to the complicated nonlinearities in the equations (see [21–24]).
In one-space dimension, the isentropic compressible fluid flow is governed by the Navier-Stokes equations:
(1)ρt+(ρu)x=0,(ρu)t+(ρu2)x+P(ρ)x=(μ(ρ)ux)x,(x,t)∈ℝ×ℝ+,
where ρ=ρ(x,t) and u=u(x,t) are the density and velocity of the fluid, μ(ρ) and P(ρ) denote the density-dependent viscosity and pressure, respectively, and the subscripts mean partial derivations. Guo and Jiang [21] considered (1) with constant viscosity, μ(ρ)≡μ>0, and linear density-dependent pressure, P(ρ)=aρ, where a>0 is a constant, and proved that there exist neither forward nor backward self-similar solutions with finite total energy. Their investigation generalized the results for 3D incompressible fluids in Nečas et al.’s work [6] to the 1D compressible case with P(ρ)=ργ, where γ=1. The problem with γ>1, however, is open. From a physical point of view, one can derive the compressible Navier-Stokes equations from the Boltzmann equations by exploiting the Chapman-Enskog expansion up to the second order and then find that the viscosity depends on the temperature. If considering an isentropic process, this dependence can be translated into that on the density, such as μ(ρ)=ρθ, where θ>0 is a constant (see [25]). Okada et al. [26] pointed out that, because of the hard sphere interaction, the relation between indices θ and γ is θ=(γ-1)/2. In the first part of this paper, we are concerned with (1) where
(2)μ(ρ)=ρθ,P(ρ)=ργ,γ-12≤θ<γ+12,γ≥1.
When considering an ideal compressible gas flow, particularly in the thermodynamic analysis with exergy loss and entropy generation, both the viscosity and pressure rely on the entropy, so it is necessary to extend the nonisentropic fluid dynamics to include the transport of entropy (see [13, 27–35]). We consider the following coupled system of the Navier-Stokes equations with an entropy transport equation in a pure form:
(3)ρst+ρusx=0,(4)ρt+(ρu)x=0,(5)(ρu)t+(ρu2)x+P(ρ,s)x=(μ(ρ,s)ux)x,(x,t)∈ℝ×ℝ+,
where s=s(x,t) is the entropy of the fluid and μ(ρ,s) and P(ρ,s) denote the density-entropy-dependent viscosity and pressure, respectively. In this system, we assume that
(6)μ(ρ,s)=ρθes,P(ρ,s)=ργes,0≤θ<γ+12,γ≥1.
Navier-Stokes equations enjoy a scaling property: if (ρ,u) solves (1)-(2), then
(7)(ρ(λ),u(λ))=(λaρ(λcx,λdt),λbu(λcx,λdt))
does so for any λ>0, by setting a=1/(γ-θ), b=(γ-1)/2(γ-θ), c=(γ+1-2θ)/2(γ-θ), and d=1. Note that, from (2), a≥c>0 and b≥0. Solution (ρ,u) is called forward self-similar if
(8)(ρ,u)=(ρ(λ),u(λ)),for everyλ>0.
In that case, ρ(x,t) and u(x,t) are decided by their values at the instant of t=1/λ:
(9)ρ(x,t)=1taQ(xtc),u(x,t)=1tbU(xtc),t>0,
where Q(y)=ρ(y,1) and U(y)=u(y,1) are defined on ℝ. In the same manner, the backward self-similar solutions are of the form:
(10)ρ(x,t)=1(T-t)aQ(x(T-t)c),u(x,t)=1(T-t)bU(x(T-t)c),0<t<T,
where Q(y)=ρ(y,T-1) and U(y)=u(y,T-1) for T>1. Substitution of (9) or (10) into (1) gives
(11)(QU)′-c(yQ)′+(c-a)Q=0,(12)(QθU′)′-(QU2)′+c(yQU)′-(Qγ)′+(a+b-c)QU=0,
for forward self-similar solutions, or
(13)(QU)′+c(yQ)′-(c-a)Q=0,(14)(QθU′)′-(QU2)′-c(yQU)′-(Qγ)′-(a+b-c)QU=0,
for backward self-similar solutions, respectively. In comparison with those in Guo and Jiang [21], forward (backward) self-similar equations above process necessary modifications and additional difficulties. For instance, (11) and (13) have solutions with an additional integral term, and thus the modified blow-up analysis needs an L∞ estimate on the density and a new large-scale argument on the energy. In addition, conditions on θ and γ proposed in (2) are directly related to the energy estimate.
Mellet and Vasseur [25] obtained the global existence of strong solutions for the Cauchy problem of (1) with positive initial density having (possibly different) positive limits at x=±∞. Precisely, fix constant positive density ρ+>0 and ρ->0, and let ρ¯(x) be a smooth monotone function satisfying
(15)ρ¯(x)=ρ±when±x≥1,ρ¯(x)>0,∀x∈ℝ.
Assume that the initial data ρ(x,t)∣t=0=ρ0(x) and u(x,t)∣t=0=u0(x) satisfy
(16)0<κ_0≤ρ0(x)≤κ¯0<∞,ρ0-ρ¯∈H1(ℝ),u0∈H1(ℝ),
for some constants κ_0 and κ¯0. Assume also that μ(ρ) and P(ρ) verify (2). Mellet and Vasseur [25] proved that there exists a global strong solution (ρ,u) of (1) on ℝ×ℝ+ such that for every T>0:
(17)ρ-ρ¯∈L∞(0,T;H1(ℝ)),u∈L∞(0,T;H1(ℝ))∩L2(0,T;H2(ℝ)).
Moreover, for every T>0, there exist uniform bounds away from zero with respect to all strong solutions having the same initial data. Precisely, there exist some constants C(T), κ_(T), and κ¯(T) depending only on T, ρ0(x), and u0(x) such that the following bounds hold uniformly for any strong solution (ρ,u):
(18)∥ρ-ρ¯∥L∞(0,T;H1(ℝ))≤C(T),∥u∥L∞(0,T;H1(ℝ))≤C(T),(19)0<κ_(T)≤ρ(x,t)≤κ¯(T),∀(x,t)∈ℝ×[0,T].
Define
(20)p(ρ∣ρ¯)=1γ-1ργ-1γ-1ρ¯γ-γγ-1ρ¯γ-1(ρ-ρ¯)
as the relative potential energy density of (1), and
(21)E(t)=∫ℝ12ρ(x,t)u2(x,t)dx
as the kinetic energy. Note that, since p is strictly convex, p(ρ∣ρ¯) is nonnegative for every ρ, and p(ρ∣ρ¯)=0 if and only if ρ=ρ¯. Mellet and Vasseur [25] also showed that, if the initial total energy is finite, that is, the sum of the kinetic and potential energy at time 0 satisfies
(22)∫ℝ[12ρ0u02+p(ρ0∣ρ¯)]dx<+∞,
then the following global-energy estimate on (-∞,∞)×[0,T] holds uniformly with respect to all strong solutions; that is, for every T>0, there exists a positive constant C(T) depending only on T, ρ0(x), and u0(x) such that
(23)sup[0,T]∫ℝ[12ρu2+p(ρ∣ρ¯)]dx+∫0T∫ℝμ(ρ)|ux|2dxdt≤C(T)
holds for any strong solution (ρ,u). Correspondingly, for R>0, 0<t1<T, and some constant C(R,t1,T), we call
(24)supt1≤t≤T∫-RR[12ρu2+p(ρ∣ρ¯)]dx+∫t1T∫-RRμ(ρ)|ux|2dxdt≤C(R,t1,T)
the local-energy estimate on [-R,R]×[t1,T]. Note that the global-energy estimate implies the local-energy estimate.
The main result for the self-similar solutions of the isentropic compressible Navier-Stokes equations is as follows.
Theorem 1.
Assume that μ(ρ) and P(ρ) in (1) verify (2). Then the following statements are true.
There is no self-similar strong solution satisfying the global-energy estimate (23).
If there is a forward (backward) self-similar strong solution satisfying the local-energy estimate (24), then its kinetic energy (21) blows up as t↓0+(t↑T-).
For the self-similar solutions of the coupled system of the nonisentropic compressible Navier-Stokes equations with an entropy transport equation, the main result is as follows.
Theorem 2.
Assume that μ(ρ,s) and P(ρ,s) in (3)–(5) verify (6). Then the following statements are true.
There is no self-similar strong solution satisfying the global-energy estimate (23).
If there is a forward (backward) self-similar strong solution satisfying the local-energy estimate (24), then its kinetic energy (21) blows up as t↓0+(t↑T-).
Theorem 1 is proved in Section 2 and Theorem 2 in Section 3.
2. Proof of Theorem 1
Any self-similar solution of (1) is either forward or backward, so we first prove Theorem 1 for forward and then for backward self-similar solutions.
2.1. Forward Self-Similar SolutionsLemma 3.
If (Q,U) solves (11)-(12), then the corresponding strong solution (ρ,u) defined by (9) of (1) does not satisfy the global-energy estimate (23).
Proof.
From (11),
(25)Q(y)U(y)=cyQ(y)+(a-c)∫0yQ(z)dz+C0,
where C0 is an arbitrary constant. From (19),
(26)0<κ_(1)≤ρ(y,1)=Q(y)≤κ¯(1),∀y∈ℝ.
Since a≥c>0, (25) and (26) imply that, for y≥Y0>0 where Y0 is large enough,
(27)Q(y)U(y)=c2yQ(y)+cy2Q(y)+(a-c)×∫0yQ(z)dz+C0≥c2yQ(y)+c2yQ(y)+C0≥c2yQ(y)+c2Y0κ_(1)+C0≥c2yQ(y).
Thus, from (9) and (21), for any t>0,
(28)E(t)=∫ℝ12·1taQ(xtc)·1t2bU2(xtc)dx=12ta+2b-c∫ℝ1Q(y)[Q(y)U(y)]2dy≥12ta+2b-c∫Y0+∞1Q(y)[c2yQ(y)]2dy≥c2κ_(1)8ta+2b-c∫Y0+∞y2dy=+∞.
This proves the lemma.
Lemma 4.
If (Q,U) solves (11)-(12) and the corresponding strong solution (ρ,u) defined by (9) of (1) satisfies the local-energy estimate (24), then as t↓0, the kinetic energy (21) must blow up.
Proof.
Similar to the proof of Lemma 3, for any t>0 and R>Y0tc>0,
(29)E(t)≥c2κ_(1)8ta+2b-c∫Y0R/tcy2dy=c2κ_(1)24ta+2b-c(R3t3c-Y03)=c2κ_(1)24ta+2b+2c[R3-(Y0tc)3]⟶+∞,ast↓0.
This proves the lemma.
2.2. Backward Self-Similar SolutionsLemma 5.
If (Q,U) solves (13)-(14), then the corresponding strong solution (ρ,u) defined by (10) of (1) does not satisfy the global-energy estimate (23).
Proof.
Fix T>1. From (13),
(30)Q(y)U(y)=-cyQ(y)-(a-c)∫0yQ(z)dz+C0,
where C0 is an arbitrary constant. From (19),
(31)0<κ_(T)≤ρ(y,T-1)=Q(y)≤κ¯(T),∀y∈ℝ.
Since a≥c>0, (30) and (31) imply that, for y≥Y0>0 where Y0 is large enough,
(32)Q(y)U(y)=-c2yQ(y)-cy2Q(y)-(a-c)×∫0yQ(z)dz+C0≤-c2yQ(y).
Thus, from (10) and (21), for any t<T,
(33)E(t)=∫ℝ12·1(T-t)aQ(x(T-t)c)·1(T-t)2bU2×(x(T-t)c)dx≥c2κ_(T)8(T-t)a+2b-c∫Y0+∞y2dy=+∞.
This proves the lemma.
Lemma 6.
If (Q,U) solves (13)-(14) and the corresponding strong solution (ρ,u) defined by (10) of (1) satisfies the local-energy estimate (24), then as t↑T, the kinetic energy (21) must blow up.
Proof.
Recalling the proofs of Lemmas 4 and 5, for R>Y0(T-t)c>0,
(34)E(t)≥c2κ_(T)8(T-t)a+2b-c∫Y0R/(T-t)cy2dy=c2κ_(T)24(T-t)a+2b+2c×[R3-Y03(T-t)3c]⟶+∞,ast↑T.
This proves the lemma.
Now, Theorem 1 follows from the four lemmas above.
3. Proof of Theorem 2
If (ρ,u,s) solves (3)–(6), then
(35)(ρ(λ),u(λ))=(λls(λcx,λdt),λls(λcx,λdt))λaρ(λcx,λdt),λbu(λcx,λdt),λls(λcx,λdt))
does so for any λ>0, by setting a=1/(γ-θ), b=(γ-1)/2(γ-θ), c=(γ+1-2θ)/2(γ-θ), d=1, and l=0. Note that, from (6), a>0, b≥0, and c>0. The forward self-similar solutions have the following form:
(36)ρ(x,t)=1taQ(xtc),u(x,t)=1tbU(xtc),s(x,t)=S(xtc),t>0,
where Q(y)=ρ(y,1), U(y)=u(y,1), and S(y)=s(y,1). The backward self-similar solutions are
(37)ρ(x,t)=1(T-t)aQ(x(T-t)c),u(x,t)=1(T-t)bQ(x(T-t)c),s(x,t)=S(x(T-t)c),0<t<T,
where Q(y)=ρ(y,T-1), U(y)=u(y,T-1), and S(y)=s(y,T-1) for T>1.
Lions [29] investigated the coupled system of the Navier-Stokes equations with an entropy transport equation in a pure form and obtained the existence of weak solutions satisfying (19) and (23).
Proof of Theorem 2.
Suppose that (ρ,u,s) is a forward self-similar solution. Inserting (36) into (3), one gets U(y)=cy, and thus u(x,t)=cx/t. Therefore, for any t>0, (36), (19), and (21) yield
(38)E(t)=∫ℝ12·1taQ(xtc)·(cxt)2≥c2κ_(1)2ta+2∫ℝx2dx=+∞.
This means that the global-energy estimate (23) does not hold and that the kinetic energy (21) blows up as t↓0.
The case of backward self-similar solutions can be proved similarly, so Theorem 2 is proved.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is partially supported by Zhejiang Provincial Natural Science Foundation of China (no. LQ13G030018) and National Natural Science Foundation of China (nos. 11001049 and 11226184).
LerayJ.Sur le mouvement d'un liquide visqueux emplissant l'espace193463119324810.1007/BF02547354MR1555394BatchelorG. K.1999Cambridge, UKCambridge University PressMR1744638SedovL. I.199310thBoca Raton, Fla, USACRC PressMR0108122BarrazaO. A.Self-similar solutions in weak Lp-spaces of the Navier-Stokes equations199612241143910.4171/RMI/202MR1402672CannoneM.PlanchonF.Self-similar solutions for Navier-Stokes equations in ℝ31996211-217919310.1080/03605309608821179MR1373769NečasJ.RužičkaM.ŠverákV.On Leray's self-similar solutions of the Navier-Stokes equations1996176228329410.1007/BF02551584MR1397564PlanchonF.Asymptotic behavior of global solutions to the Navier-Stokes equations in ℝ31998141719310.4171/RMI/235MR1639283TsaiT.-P.On Leray's self-similar solutions of the Navier-Stokes equations satisfying local energy estimates19981431295110.1007/s002050050099MR1643650TasiT.-P.Forward discrete self-similar solutions of the Navier-Stokes equationssubmitted, http://arxiv.org/abs/1210.2783MillerJ. R.O'LearyM.SchonbekM.Nonexistence of singular pseudo-self-similar solutions of the Navier-Stokes system2001319480981510.1007/PL00004460MR1825409Lemarié-RieussetP. G.2002431Boca Raton, Fla, USAChapman & Hall/CRCChapman & Hall/CRC Research Notes in Mathematics10.1201/9781420035674MR1938147ChaeD.Nonexistence of asymptotically self-similar singularities in the Euler and the Navier-Stokes equations2007338243544910.1007/s00208-007-0082-6MR2302070ChaeD.Notes on the incompressible Euler and related equations on RN200930551352610.1007/s11401-009-0107-4MR2601484ChaeD.Nonexistence of self-similar singularities in the ideal magnetohydrodynamics200919431011102710.1007/s00205-008-0182-9MR2563631BrandoleseL.Fine properties of self-similar solutions of the Navier-Stokes equations2009192337540110.1007/s00205-008-0149-xMR2505358JiaoX.-Y.Some similarity reduction solutions to two-dimensional incompressible Navier-Stokes equation200952338939410.1088/0253-6102/52/3/02MR2641067BarnaI. F.Self-similar solutions of three-dimensional Navier—stokes equation20115647457502-s2.0-8005479543210.1088/0253-6102/56/4/25KarchG.PilarczykD.Asymptotic stability of Landau solutions to Navier-Stokes system2011202111513110.1007/s00205-011-0409-zMR2835864ŠverákV.On Landau's solutions of the Navier-Stokes equations2011179120822810.1007/s10958-011-0590-5MR3014106JiaH.ŠverákV.Local-in-space estimates near initial time for weak solutions of the Navier-Stokes equations and forward self-similar solutions201310.1007/s00222-013-0468-xGuoZ.JiangS.Self-similar solutions to the isothermal compressible Navier-Stokes equations200671565866910.1093/imamat/hxl013MR2268881GuoZ.XinZ.Analytical solutions to the compressible Navier-Stokes equations with density-dependent viscosity coefficients and free boundaries2012253111910.1016/j.jde.2012.03.023MR2917399YuenM.Self-similar solutions with elliptic symmetry for the compressible Euler and Navier-Stokes equations in ℝN201217124524452810.1016/j.cnsns.2012.05.022MR2960245AnH.YuenM.Supplement to ‘Self-similar solutions with elliptic symmetry for the compressible Euler and Navier-Stokes equations in ℝN’20131861558156110.1016/j.cnsns.2012.10.001MR3016906MelletA.VasseurA.Existence and uniqueness of global strong solutions for one-dimensional compressible Navier-Stokes equations2008394134413652-s2.0-4474908600110.1137/060658199OkadaM.Matušocircu-NečasováS.MakinoT.Free boundary problem for the equation of one-dimensional motion of compressible gas with density-dependent viscosity200248120MR1980822TengH.KinoshitaC. M.MasutaniS. M.ZhouJ.Entropy generation in multicomponent reacting flows199812032262322-s2.0-0006615556HirschfelderJ. O.CurtissC. F.BirdR. B.1954New York, NY, USAJohn Wiley & SonsLionsP.-L.19962New York, NY, USAOxford University PressOxford Lecture Series in Mathematics and Its ApplicationsMR1637634NishidaK.TakagiT.KinoshitaS.Analysis of entropy generation and exergy loss during combustion200229869874LiZ. W.ChouS. K.ShuC.YangW. M.Entropy generation during microcombustion20059782-s2.0-2144443214310.1063/1.1876573084914YapıcıH.KayataşN.AlbayrakB.BaştürkG.Numerical calculation of local entropy generation in a methaneair burner20054618851919LiorN.Sarmiento-DarkinW.Al-SharqawiH. S.The exergy fields in transport processes: their calculation and use20063155535782-s2.0-2914445543810.1016/j.energy.2005.05.009StanciuD.MarinescuM.DobrovicescuA.The influence of swirl angle on the irreversibilities in turbulent diffusion flames20071041431532-s2.0-38549155599SafariM.SheikhiM. R. H.JanbozorgiM.MetghalchiH.Entropy transport equation in large eddy simulation for exergy analysis of turbulent combustion systems20101234344442-s2.0-7795356593710.3390/e12030434