Perturbational blowup solutions to the compressible Euler equations with damping

Background The N-dimensional isentropic compressible Euler system with a damping term is one of the most fundamental equations in fluid dynamics. Since it does not have a general solution in a closed form for arbitrary well-posed initial value problems. Constructing exact solutions to the system is a useful way to obtain important information on the properties of its solutions. Method In this article, we construct two families of exact solutions for the one-dimensional isentropic compressible Euler equations with damping by the perturbational method. The two families of exact solutions found include the cases \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma >1$$\end{document}γ>1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma =1$$\end{document}γ=1, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma$$\end{document}γ is the adiabatic constant. Results With analysis of the key ordinary differential equation, we show that the classes of solutions include both blowup type and global existence type when the parameters are suitably chosen. Moreover, in the blowup cases, we show that the singularities are of essential type in the sense that they cannot be smoothed by redefining values at the odd points. Conclusion The two families of exact solutions obtained in this paper can be useful to study of related numerical methods and algorithms such as the finite difference method, the finite element method and the finite volume method that are applied by scientists to simulate the fluids for applications.

by the adiabatic γ-law. The constant α ≥ 0 is the damping coefficient.
System (1) is one of the most fundamental equations in fluid dynamics. Many interesting fluid dynamic phenomena can be described by system (1) (Lions , 1998a;Lions , 1998b). The Euler equations (α = 0) are also the special case of the noted Navier-Stokes equations, whose problem of whether there is a formation of singularity is still open and long-standing. Thus, the singularity formation in fluid mechanics has been attracting the attention of a number of researchers (Sideris 1985;Xin 1998;Suzuki 2013;Lei et al. 2013;Li and Wang 2006;Li et al. 2013).
Among others, we mention that in 2003, Sideris-Thomases-Wang (Sideris et al. 2003) obtained results for the three dimensional compressible Euler equations with a linear damping term with assumption γ > 1, that is, system (1) with N = 3 and γ > 1. They discovered that damping prevents the formation of singularities in small amplitude flows, but large solutions may still break down. They formulated the Euler system as a symmetric hyperbolic system, established the finite speed of propagation of the solution, and some energy estimates to obtain local existence as well as global existence of the solution. For larger solution, they showed that the solution will blow up in a finite time by establishing certain differential inequalities.
In this article, we consider the one dimensional case of system (1): More precisely, we apply the perturbational method to obtain the following main results. Theorem 1 For system (3) with γ > 1 and α > 0, one has the following family of exact solutions with parameters ξ , ρ(0, 0) > 0, a 0 > 0 and a 1 : where ρ γ −1 (t, 0) is given by a(t) and b(t) satisfy the following ordinary differential equations: Remark 2 The ordinary differential equation (O.D.E.) (6) will be analyzed in section 2 and it is well-known by the theory of ordinary differential equations that the solutions of system (7) exist and is C 2 as long as f and g, which are functions of ä, ȧ and a, are continuous.
Theorem 3 For the family of exact solutions in Theorem 1, we have the following five cases.
(iii) If ξ < 0, then the solution (4) blows up on a finite time.
Moreover, we show that the singularity formations in the cases iii), iv) and v) above are of essential type in the sense that the singularities cannot be smoothed by redefining values at the odd points. This is an improvement of the corresponding results in Yuen (2011).

Theorem 4 For system
where a(t) and b(t) satisfy the following ordinary differential equations: Theorem 5 For the family of exact solutions in Theorem 4, we have the following five cases.

Analysis of an O.D.E.
Consider the following initial value problem.
where γ ≥ 1, α > 0 and ξ ∈ R are constants. We set Lemma 6 For system (12), if T * is finite, then the one-sided limit Proof Note that we always have Suppose lim t→T * a(t) > 0. Then we can extend the solution of (12) to [0, T * + ε) by solving the following system.
This contradicts the definition of T * . Thus, the lemma is established. (12), we have the following three cases.

Thus,
As A/α < 0, we have a(t) < 0 for all sufficiently large t. This is impossible as T * = +∞. Thus, Case 3. is established.

Remark 8
The case for ξ = 0 will be analyzed in the proof of Theorem 3.

Proofs of the Theorems
Proof of Theorem 1 We divide the proof into steps.
Step 1. In the first step, we show a lemma.

Thus, we have
Taking integration with respect to x, we obtain On the other hand, multiplying ρ γ −2 on both sides of (3) 1 , we get From (33), we have and Substituting (36), (32) and (35) into (34), one obtains the relation claimed in the lemma.
Step 2. We set where c := c(t) and b := b(t) are functions of t. Then, (29) is transformed to where we arrange the terms according to the coefficients of x.
Step 3. We use the Hubble transformation: and set the coefficient of (38) to be zero. Thus,

Note that we have the novel identity
Multiplying the both sides of (40) by a γ +1 , it becomes for some constant ξ.
Step 4. With (39), we set the coefficient of x in (38) to be zero. Thus, b satisfies where Last Step. With (39) and setting the coefficient of 1 in (38) to be zero, we are required to solve where (38) ... a a + αä a + γȧä a 2 + γ αȧ 2 a 2 = 0. Next, we prove Theorem 3 as follows.
Proof of Theorem 3 For ξ > 0, case i) and case ii) of Theorem 3 follow from Case 1. and Case 2. of Lemma 7.
For ξ < 0, by Case 3. of Lemmas 6 and7, there exists a finite T * > 0 such that the onesided limit of a(t) is zero as t approaches to T * . It remains to show T * is not a removable singularity of ȧ/a. To this end, suppose one has Then, Thus, the singularity is of essential type and case iii) of Theorem 3 is proved.
For ξ = 0, (6) 1 becomes which can be solved by using integral factor. The solution is Thus, a(T ) = 0 if a 1 > 0. Also, a(T ) = 0 if a 1 < 0 and a 0 < −a 1 /α, where T := 1 α ln a 1 a 1 +a 0 α > 0. As (T, x) is an essential singularity of u(t, x) for any x. Thus, cases iv) and v) of Theorem 3 are established. The proof is complete.

Proof of Theorems 4 and 5
The corresponding relation of Lemma 9 for γ = 1 is With similar steps, one can obtain the family of exact solutions in Theorem 4.
Note that (10) is a special case of (12) and the arguments in the proof of Theorem 3 hold for γ = 1. Thus, the results for Theorem 5 follows.

Conclusion
The complicated Euler equations with a damping term (1) do not have a general solution in a closed form for arbitrary well-posed initial value problems. Thus, numerical methods and algorithms such as the finite difference method, the finite element method and the finite volume method are applied by scientists to simulate the fluids for applications in real world. Thus, our exact solutions in this article provide concrete examples for researchers to test their numerical methods and algorithms.