Abstract
We study the initial boundary value problem of a compressible non-Newtonian fluid. The system describes the motion of the compressible viscous isentropic gas flow driven by the non-Newtonian self-gravitational force. The existence of strong solutions are derived in one dimensional bounded intervals by constructing a semi-discrete Galerkin scheme. Moreover, the uniqueness of solutions are also investigated. The main point of the study is that the viscosity term and potential term are fully nonlinear, and the initial vacuum is allowed.
1 Introduction
In mathematical physics, the Navier-Stokes equation is known as one of the most fundamental equations in fluid mechanics. The compressible isentropic Navier-Stokes equations, which are the basic models describing the evolution of a viscous compressible fluid in a domain x ∈ Ω, read as follows:
where the unknowns ρ, u represent the density and the velocity of the fluid, respectively. Here D(u) = (∇u + (∇u)T)/2 is the strain tensor and P(ρ) = aργ (a > 0, γ > 1) is the pressure, μ and λ are the viscosity constants which satisfy the physical requirements μ ≥ 0 and 2μ + 3λ ≥ 0. The Navier-Stokes equations are the equations governing the motion of usual fluids like water, air, oil etc., under quite general conditions, and they appear in the study of many important phenomena, either alone or coupled with other equations. For instance, they are used in theoretical studies in aeronautical sciences, in meteorology, in thermo-hydraulics, in the petroleum industry, in plasma physics, etc. From the point of view of continuum mechanics the Navier-Stokes equations are essentially the simplest equations describing the motion of a fluid, and they are derived under a quite simple physical assumption, namely, the existence of a linear local relation between stresses and strain rates. The compressible isentropic Navier-Stokes equations (1.1) are derived from the conservation laws of mass and the balance of momentum, (for details see [1, 2]). While the physical model leading to the Navier-Stokes equations is simple, the situation is quite different from the mathematical point of view. In particular, because of the nonlinearity, the mathematical study of these equations is difficult and requires the full power of modern functional analysis. A major question is whether the solution remains smooth all the time. These and other related questions are interesting not only for mathematical understanding of the equations but also for understanding the phenomenon of turbulence.
There are huge literatures on the study of the existence and behavior of solutions to Navier-Stokes equations. Some of the previous relevant works in this direction can be summarized as follows. For instance, the 1D version of (1) were addressed by Kazhikhov et al. in [3] for sufficiently smooth data, and by Hoff [4] for discontinuous initial data, where the data is uniformly away from the vacuum; the existence of global weak solutions for isentropic flow were investigated by Lions in [5] by using the weak convergence method. In [6], the authors employed a new method to prove the existence and uniqueness of local strong solutions in the case where the initial data satisfies some compatibility conditions. The dynamics of weak solutions and vacuum states were investigated in [7] for the 1D compressible Navier-Stokes equations with density-dependent viscosity in bounded spatial domains or periodic domains. For other results we refer the reader to [8-13] and the references cited therein.
The above references mainly concerned the fluid which the relation between the stress and rate of strain is linear, namely, the Newtonian fluid. The study of non-Newtonian fluids (the relation between the stess and rate of strain is not linear) mechanics is of great significance because such fluids describe more realistic phenomenon. These flows are frequently encountered in many physical and industrial processes [14], such as porous flows of oils and gases [15], biological fluid flows of blood [16], saliva and mucus, penetration grouting of cement mortar and mixing of massive particles and fluids in drug production [17]. Many studies are based on the field of non-Newtonian flows, both theoretically and experimentally. In [18], Ladyzhenskaya first proposed a special form for the incompressible model, namely that the viscous stress tensor Γij = (μ0 + μ1|E(∇u)p−2|)Eij(∇u). For μ0 = 0, if p < 2, it is a pseudo-plastic fluid, and when p > 2, it is a dilatant fluid (see also [19]). From a Physics point of view, the model captures the shear thinning fluid for the case of 1 < p < 2, and captures the shear thickening fluid for the case of p > 2. In [20], the trajectory attractor and global attractor for an autonomous non-Newtonian fluid in dimension two was studied. The existence and uniqueness of solutions for non-Newtonian fluids were established in [21] by applying the Ladyzhenskaya’s viscous stress tensor model. Then the global existence and exponential stability of solutions to the one-dimensional full non-Newtonian fluids were investigated in [22]. Recently, in [23], the authors present a decoupling multiple-relaxation-time lattice Boltzmann flux solver for simulating non-Newtonian power-law fluid flows.
On the other hand, another basiclly important case is when the motion of compressible viscous isentropic flow is driven by a self-gravitational force. However, at the present, little is known yet on the strong solutions to system (1) with non-Newtonian self-gravitational potential on bounded domain, even in one dimensional case. In this paper we focus on the following 1D system of compressible non-Newtonian equations
in ΩT = Ω × (0, T) with the initial and boundary conditions
Here, ρ, u, Φ denote the density, velocity and the non-Newtonian gravitational potential, respectively. P = aργ (a > 0, γ > 1) is the pressure, the initial density
We state the main results as follows:
Theorem 1.1
Assume that (ρ0, u0) satisfies the following conditions
If there is a functiong ∈ L2(Ω), such that the initial data satisfy the following compatibility condition:
then there exist a time T* ∈ (0, +∞) and a unique strong solution (ρ, u, Φ) to (2)-(4) such that
The rest of the paper is organized as follows. After stating the notations, in Section 2, we first present some useful lemmas, then the analysis of a priori estimates for smooth solutions are derived. In Section 3, we give the proof of existence, and finally complete the proof of uniqueness of the main theorem in Section 4.
In what follows, we use the following abbreviations for simplicity of notation:
Throughout this paper, we will omit the variables t, x of functions if it does not cause any confusion. We use C to denote a generic constant that may vary in different estimates.
2 A priori estimates for smooth solutions
In this section, we provide some known facts that will be used in the proof of the main result.
Lemma 2.2
([21]). Letρ0 ∈ H1(Ω), u0 ∈
Then there is a subsequences
To prove the existence of strong solutions, we require some more regular estimates. Next, we derive some a priori estimates for smooth solutions which are crucial to prove the local existence of strong solutions.
Let (ρ, u, Φ) be a smooth solution of (2)-(4) and ρ0 ≥ δ, where 0 < δ ≪ 1 is a positive number, m0 := ∫Ωρ0(x)dx be initial mass and m0 > 0. Throughout the paper, we will denote
As stated above, we need to estimate the uniform bound of the approximate solutions. Now, we consider the following linearized problem
with the initial and boundary conditions
where
and u0 ∈
We will prove the existence of solutions for (5)-(9) by virtue of the uniform estimates which do not depend on ε, depending only on M0, and prove the limit of the approximate solutions is the solution of problem (2)-(4) with vacuum.
It follows from (10) and Young’s inequality that there exists a constant C depending only on M0, such that
We construct an auxiliary function
Our derivation will be based on the local boundedness of Ψ(t). Before we estimate each term of Ψ, we need to do the estimate of |uxx|L2. Firstly, multiplying (7) by Φ and integrating over Ω and using Young’s inequality,
Then, it follows from equation (6) and (5) that
Taking it by L2-norm, using Young’s inequality, we obtain
which along with (11), implies that
Estimate for |ρ|H1
We are going to estimate the first term of Ψ(t). Multiplying (5) by ρ and integrating over Ω with respect to x, we obtain from Sobolev inequality
Differentiating (5) with respect to x, multiplying it by ρx and integrating over Ω on x, and using Sobolev inequality, we have
Together with (13), (14) and Gronwall’s inequality, it follows that
Using (5) we obtain
Besides, differentiating (7) with respect to time t, multiplying it by Φt and integrating over Ω with respect to x and using Young’s inequality, we have
where C is a positive constant depending only on M0.
Estimate for | u | W 0 1 , p
We turn to the second term of Ψ(t). Multiplying (6) by ut, integrating over ΩT, together with (5), (10), Sobolev inequality and Young’s inequality, we obtain
In the second inequality we have used
where C is a positive constant, depending only on M0.
Estimate for | ρ u t | L 2
We estimate the last term of Ψ(t). Differentiating (6) with respect to t yields
Multiplying it by ut and integrating over Ω with respect to x, we derive
Note that
Let
From (12), it follows that
Then, from (20) and (5), (19) can be rewritten as
Using Sobolev inequality, Young’s inequality, (6) and (12), we obtain
Substituting Ij(j = 1, 2, …, 9) into (21) and integrating over (τ, t) ⊂ (0, T), we have
Next, we estimate
Multiplying (6) by ut and integrating over Ω, we obtain
Since (ρ, u, Φ) is a smooth solution, we have
Then taking a limit on τ for (22) as τ → 0, we get
which combined with (12), (15)-(17), (18) and the definition of Ψ(t) leads to
where C, c̃ are positive constants, depending only on M0.
In view of this inequality, we can find a time T* > 0 and a constant C, such that
where C is a positive constant depending only on M0.
3 Proof of the existence
Since a priori estimates for higher regularity have been derived, the existence of strong solutions can be established by a standard argument, in the case of bounded domains, we construct approximate solutions via a semi-discrete Galerkin scheme, derive uniform bounds and thus obtain solutions by passing to the limit. Our method that constructed approximate systems is similar to that in [6]. To implement a semi-discrete Galerkin scheme, we take our basic function space as X =
Let ρ0, u0, Φ0, and f be functions satisfying the hypothes of theorem, assume for the moment that
The initial and boundary conditions are
Under the hypotheses of the theorem, similarly, for any fixed δ > 0, we may get the similar estimate
Since C does not depend on ε, δ and m (for any m > M, M is dependent on the approximate velocity of the initial condition). We can deduce from the uniform estimate(24) that (ρm, um, Φm) converges, up to an extraction of subsequences. Let m → ∞. We obtain the following estimates in the obvious weak sense
For each small δ > 0,
where
We deduce that (ρδ, uδ, Φδ) is a solution of the following initial and boundary value problem
where
Together with Lemma 2.2, there is s subsequence
where C is a positive constant depending only on M0.
4 Proof of the uniqueness
We now prove the uniqueness results. Let (ρ, u, Φ), (ρ̄, ū, Φ̄) be two solutions of the problem (2)-(4). Combining (2)1 and (2)2,
we show that
Multiplying the above equation by (u – ū) and integrating over Ω, we obtain
Since
Thus
Moreover, from (2)3, by a direct calculation, we have
Multiplying it by (Φ – Φ̄) and integrating over Ω, we obtain
Note that
Together with (27) and (28), we arrive at
Consequently, (26) can be rewritten as
On the other hand, from the conservation of mass equation (2)1, using the identity
Multiplying the above equation by (ρ – ρ̄) and integrating it over Ω, we obtain
Thus
Furthermore, (5)1 implies
Multiplying it by (P – P̄) and integrating over Ω, we get
From this and the Grownwall’s inequality, yields
which means
This complete the proof of the main theorem.
Acknowledgement
The authors would like to thank the anonymous referees for their valuable suggestions and comments. This work was supported by the National Natural Science Foundation of China (Nos. 11526105; 11572146), the Education Department Foundation of Liaoning Province (No. JQL201715411) and the Natural Science Foundation of Liaoning Province (No. 20180550585).
References
[1] Chorin A. J., Marsden J. E., A Mathematical Introduction to Fluid Mechanics, Springer-Verlag, New York, (1993).10.1007/978-1-4612-0883-9Search in Google Scholar
[2] Teman R., Navier-Stokes Equations, Theory and Numerical Analysis, North-Holland, Amsterdam, (1984).Search in Google Scholar
[3] Kazhikhov A. V., Shelukhin V V: Unique global solution with respect to time of initial-boundary value problems for one-dimensional equations of a viscous gas. J. Appl. Math. Mech., 1977, 41(2), 273-282.10.1016/0021-8928(77)90011-9Search in Google Scholar
[4] Hoff D., Global existence for 1D, compressible, isentropic Navier-Stokes equations with large initial data. Trans. Amer. Math. Soc., 1987, 303(1), 169-181.10.1090/S0002-9947-1987-0896014-6Search in Google Scholar
[5] Lions P. L., Mathematical topics in fluids mechanics, Vol.2, Oxford Lecture Series in Mathematics and Its Applications, Vol. 10, Clarendon Press, Oxford, (1998).Search in Google Scholar
[6] Choe H., Kim H., Strong solutions of the Navier-Stokes equations for isentropic compressible fluids, J.Differential Equations, 2003, 190, 504-523.10.1016/S0022-0396(03)00015-9Search in Google Scholar
[7] Li H. L., Li J., Xin Z. P., Vanishing of vacuum states and blow-up phenomena of the compressible Navier-Stokes equations, Commun. Math. Phys, 2008, 281, 401-444.10.1007/s00220-008-0495-4Search in Google Scholar
[8] Xin Z. P., Blow up of smooth solutions to the compressible Navier-Stokes equation with compact density, Comm. Pure Appl. Math., 1998, 51, 229-240.10.1002/(SICI)1097-0312(199803)51:3<229::AID-CPA1>3.0.CO;2-CSearch in Google Scholar
[9] Jiang S., Large-time behavior of solutions to the equations of a one-dimensional viscous polytropic ideal gas in unbounded domains. Commun. Math. Phys., 1999, 200, 181-193.10.1007/s002200050526Search in Google Scholar
[10] Feireisl E., Novotný A., Petzeltová H., On the existence of globally defined weak solution to the Navier-Stokes equations, J. Math. Fluid Mech., 2001, 3, 358-392.10.1007/PL00000976Search in Google Scholar
[11] Jiang S., Zhang P., On spherically symmetric solutions of the compressible isentropic Navier-Stokes equations, Commun. Math. Phys., 2001, 215, 559-581.10.1007/PL00005543Search in Google Scholar
[12] Yin J. P., Tan Z., Local existence of the strong solutions for the full Navier-Stokes-Poisson equations, Nonlinear Anal., 2009, 71, 2397-2415.10.1016/j.na.2009.01.074Search in Google Scholar
[13] Li J., Liang Z. L., Some uniform estimates and large-time behavior of solutions to one-dimensional compressible Navier-Stokes systems in unbounded domains with large data, Arch. Rational Mech. Anal., 2016, 220, 1195-1208.10.1007/s00205-015-0952-0Search in Google Scholar
[14] Chhabra R. P., Richardson J. F., Non-Newtonian Flow and Applied Rheology (Second Edition), Oxford, (2008).Search in Google Scholar
[15] Chevalier T., Rodts S., Chateau X., Chevalier C. and Coussot P., Breaking of non-Newtonian character in flows through a porous medium, Physical Review E, 2014, 89, 023002.10.1103/PhysRevE.89.023002Search in Google Scholar PubMed
[16] Yun B. M., Dasi L. P., Aidun C. K. and Yoganathan A. P., Computational modelling of flow through prosthetic heart valves using the entropic lattice-Boltzmann method, Journal of Fluid Mechanics, 2014, 743, 170-201.10.1017/jfm.2014.54Search in Google Scholar
[17] Gachelin J., Mino G., Berthet H., Lindner A., Rousselet A. and Clement E., Non-Newtonian viscosity of Escherichia coli suspensions, Physical Review Letters, 2013, 110, 268103.10.1103/PhysRevLett.110.268103Search in Google Scholar PubMed
[18] Ladyzhenskaya O. A., New equations for the description of viscous incompressible fluids and solvability in the large of the boundary value problems for them. In Boundary Value Problems of Mathematical Physics, vol. V, Amer. Math. Soc., Providence, RI. (1970) 95-118.Search in Google Scholar
[19] Böhme G., Non-Newtonian fluid Mechanics, Appl.Math.Mech., North-Holland, Amsterdam (1987).Search in Google Scholar
[20] Zhao C., Zhou S., Li Y., Trajectory attractor and global attractor for a two-dimensional incompressible non-Newtonian fluid. J. Math. Anal. Appl., 2007, 325, 1350-1362.10.1016/j.jmaa.2006.02.069Search in Google Scholar
[21] Yuan H. J., Xu X. J., Existence and uniqueness of solutions for a class of non-Newtonian fluids with singularity and vacuum, J. Differential Equations, 2008, 245, 2871-2916.10.1016/j.jde.2008.04.013Search in Google Scholar
[22] Qin Y., Liu X., Yang X., Global existence and exponential stability of solutions to the one-dimensional full non-Newtonian fluids. Nonlinear Anal., Real World Appl., 2012, 13, 607-633.10.1016/j.nonrwa.2011.07.053Search in Google Scholar
[23] Wang Y., Shu C., Yang L. M., Yuan H. Z., A Decoupling Multiple-Relaxation-Time Lattice Boltzmann Flux Solver for Non-Newtonian Power-Law Fluid Flows, Journal of Non-Newtonian Fluid Mechanics, 2016, 235, 20-28.10.1016/j.jnnfm.2016.03.010Search in Google Scholar
[24] Málek J., Nečas J., Rokyta M., Ru̇žička M., Weak and Measure-Valued Solutions to Evolutionary PDEs. Chapman and Hall, New York, (1996).10.1007/978-1-4899-6824-1Search in Google Scholar
© 2018 Song et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.