Cartesian vector solutions for $ N $-dimensional non-isentropic Euler equations with Coriolis force and linear damping

Xitong Liu; Xiao Yong Wen; Manwai Yuen; Xitong Liu; Xiao Yong Wen; Manwai Yuen

doi:10.3934/math.2023877

AIMS Mathematics

2023, Volume 8, Issue 7: 17171-17196. doi: 10.3934/math.2023877

Previous Article Next Article

Theory article Special Issues

Cartesian vector solutions for $N$ -dimensional non-isentropic Euler equations with Coriolis force and linear damping

1.
Department of Mathematics and Information Technology, The Education University of Hong Kong, 10 Lo Ping Road, Tai Po, New Territories, Hong Kong, China
2.
School of Applied Science, Beijing Information Science and Technology University, Beijing 100192, China

Received: 03 March 2023 Revised: 27 March 2023 Accepted: 27 March 2023 Published: 17 May 2023
MSC : 35C05, 35Q31, 76B03, 76M60

In this paper, we construct and prove the existence of theoretical solutions to non-isentropic Euler equations with a time-dependent linear damping and Coriolis force in Cartesian form. New exact solutions can be acquired based on this form with examples presented in this paper. By constructing appropriate matrices $A(t)$ , and vectors ${\mathbf{b} }(t)$ , special cases of exact solutions, where entropy $s = \ln\rho$ , are obtained. This is the first matrix form solution of non-isentropic Euler equations to the best of the authors' knowledge.

Keywords:

non-isentropic fluids,
Euler equations,
Coriolis force,
linear-damping symmetric and anti-symmetric matrices,
curve integration,
Cartesian vector form solutions

Citation: Xitong Liu, Xiao Yong Wen, Manwai Yuen. Cartesian vector solutions for $N$ -dimensional non-isentropic Euler equations with Coriolis force and linear damping[J]. AIMS Mathematics, 2023, 8(7): 17171-17196. doi: 10.3934/math.2023877

Related Papers:

[1]	Jingye Zhao, Zonghua Wei, Jiahui Liu, Yongqiang Fan . Vanishing magnetic field limits of solutions to the non-isentropic Chaplygin gas magnetogasdynamics equations. AIMS Mathematics, 2025, 10(1): 1675-1703. doi: 10.3934/math.2025077
[2]	Noufe H. Aljahdaly . Study tsunamis through approximate solution of damped geophysical Korteweg-de Vries equation. AIMS Mathematics, 2024, 9(5): 10926-10934. doi: 10.3934/math.2024534
[3]	Shang Mengmeng . Large time behavior framework for the time-increasing weak solutions of bipolar hydrodynamic model of semiconductors. AIMS Mathematics, 2017, 2(1): 102-110. doi: 10.3934/Math.2017.1.102
[4]	Aidi Yao . Two-dimensional pseudo-steady supersonic flow around a sharp corner for the generalized Chaplygin gas. AIMS Mathematics, 2022, 7(7): 11732-11758. doi: 10.3934/math.2022654
[5]	Changdev P. Jadhav, Tanisha B. Dale, Vaijanath L. Chinchane, Asha B. Nale, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On solutions of fractional differential equations for the mechanical oscillations by using the Laplace transform. AIMS Mathematics, 2024, 9(11): 32629-32645. doi: 10.3934/math.20241562
[6]	Lining Tong, Li Chen, Simone Göttlich, Shu Wang . The global classical solution to compressible Euler system with velocity alignment. AIMS Mathematics, 2020, 5(6): 6673-6692. doi: 10.3934/math.2020429
[7]	Shaomin Wang, Cunji Yang, Guozhi Cha . On the variational principle and applications for a class of damped vibration systems with a small forcing term. AIMS Mathematics, 2023, 8(9): 22162-22177. doi: 10.3934/math.20231129
[8]	Ahmed E. Abouelregal, Khalil M. Khalil, Wael W. Mohammed, Doaa Atta . Thermal vibration in rotating nanobeams with temperature-dependent due to exposure to laser irradiation. AIMS Mathematics, 2022, 7(4): 6128-6152. doi: 10.3934/math.2022341
[9]	Waleed Hamali, Abdulah A. Alghamdi . Exact solutions to the fractional nonlinear phenomena in fluid dynamics via the Riccati-Bernoulli sub-ODE method. AIMS Mathematics, 2024, 9(11): 31142-31162. doi: 10.3934/math.20241501
[10]	Kunquan Li . Analytical solutions and asymptotic behaviors to the vacuum free boundary problem for 2D Navier-Stokes equations with degenerate viscosity. AIMS Mathematics, 2024, 9(5): 12412-12432. doi: 10.3934/math.2024607

Abstract

1. Introduction

The non-isentropic Euler equations in $\mathbf{R}^{N}$ in fluid dynamics with a time-dependent linear damping and Coriolis force can be expressed as follows:

$\begin{align} & \rho_{t}+\mathrm{div}(\rho{\mathbf{u}}) = 0, \end{align}$

(1.1)

$\begin{align} & (\rho{\mathbf{u}})_{t}+\mathrm{div}(\rho{\mathbf{u}} \otimes{\mathbf{u}})+\rho J{\mathbf{u}}+\alpha(t)\rho {\mathbf{u}}+\triangledown p = \mathbf{0}, \end{align}$

(1.2)

$\begin{align} & S_t+{\mathbf{u}}\cdot\triangledown S = 0, \end{align}$

(1.3)

where ${\mathbf{u}} = (u_{1}, u_{2}, \cdots, u_{N})^{T}$ is an $\mathit{N}$ -dimensional velocity field, $\rho({\mathbf{x}}, t)$ and $p({\mathbf{x}}, t) = e^S\rho^\gamma$ represent density and the pressure function respectively, $J^{T} = -J$ representing Corilis force is an anti-symmetric matrix. The damping term $\alpha(t)\rho{\mathbf{u}}$ with $\alpha(t) \geq 0$ as a coefficient of friction is proportional to the momentum.

For the special case when $\alpha(t) = 0$ , the equations are reduced to Euler equations extended and governed by Coriolis rotational force ^[1,2,3,4]. The theoretical global existence of the Euler equations with rotational forces can be referred to ^[5,6,7]. Further studies on stability and tropical cyclones driven by this model can be referred to ^{[8,9,10,11,12,13]}.

If $J = 0$ , (1.1)–(1.3) are reduced to non-isentropic linear-damped Euler equations, which provide an important model regarding to its physical behaviours. The system can also be used to describe compressible gas dynamics through a porous material driven by a friction force ^[14,15,16]. Weak solutions of the damped Euler equations are shown with asymptotic and large-time behavious in ^{[16,17,18,19]}. Chow, Fan, and Yuen, in 2017, constructed the solutions of Cartesian form with $J = 0$ in ^[20], which can be regarded as a special case in this article, while taking the parameter $\gamma$ and $2\alpha$ in ^[20] to be $\gamma+1$ and $\alpha$ respectively. For time-dependent damping, Dong and Li studied a class of analytical solutions with free-boundary ^[21] in 2022.

For the case with $J = 0$ and $\alpha(t) = 0$ , the system (1.1)–(1.3) is reduced to the Euler equations

$\begin{align} & \rho_{t}+\mathrm{div}(\rho{\mathbf{u}}) = 0, \end{align}$

(1.4)

$\begin{align} & (\rho{\mathbf{u}})_{t}+\mathrm{div}(\rho{\mathbf{u}} \otimes{\mathbf{u}})+\triangledown p = \mathbf{0}, \end{align}$

(1.5)

$\begin{align} & S_t+{\mathbf{u}}\cdot\triangledown S = 0. \end{align}$

(1.6)

There are lots of researches on Euler equations, for example, see ^{[22,23,24,25,26]}. Among all the topics, constructing analytical and exact solutions are crucial ^{[27,28,29,30,31,32,33,34]} with a common pattern of the velocity function ${\mathbf{u}}$ in linear form in many previous studies. For non-isentropic Euler equations, Barna and Mátyás presented the analytic solutions for one-dimensional Euler equations and three-dimensional Navier-Stokes equations with polytropic equation of state ^[34,35], which can be referred to by taking $n\neq\gamma$ and the viscosities to be zero respectively. Based on the linear form of velocity, An, Fan, and Yuen contributed with Cartesian rotational solutions to the $\mathit{N}$ -dimension isentropic compressible Euler equations (1.4)–(1.6) ^[36] in 2015:

$\begin{equation} {\mathbf{u}} = {\mathbf{b}}(t)+A(t){\mathbf{x}} , \end{equation}$

(1.7)

where ${\mathbf{b}}(t)$ and $A(t)$ are vector and matrix respectively. Further studies have shown the existence of general solutions in Cartesian form to isentropic Euler equations with damping and rotational forces in ^[20] and ^[37], respectively.

Referring to the many blowup pheonomena studies ^[38,39,40], the global solution is still complicated to look for.

2. Materials and methods

In this article, the existence of a form of Cartesian solutions to non-isentropic Euler equations with rational force and linear damping (1.1)–(1.3) is proven by adopting mainly techniques on matrices, vectors, and curve integration. Enforcing $e^S = \rho$ and regarding velocity field ${\mathbf{u}}$ as an linear transformation of ${\mathbf{x}}\in\mathbf{R}^{N}$ , the problem is equivalent to finding the pressure function $p$ , which leads us to a quadratic form and requirments on the matrix $A$ and vector ${\mathbf{b}}$ . With this finding, we can construct some special exact solutions, which could be utilized in benchmarks for testings, simulations of computing flows.

In the following sections, we will prove the existence of the non-isentropic damped Euler equations with Coriolis forces, which admit Cartesian solutions by using appropriate requirements on matrix $A$ and vector ${\mathbf{b}}$ . We will give examples on this first cartesian form solutions to non-isentropic Euler equations based on our finding.

3. Results

In this section, we consider the non-isentropic Euler equations. Suppose that the density $\rho$ and pressure $p$ satisfy the relation

$\begin{equation} p(\rho) = e^S\rho^{\gamma}\text{,} \end{equation}$

(3.1)

where the constant $\gamma = c_{p}/c_{u}\geq1$ , and $c_{p}$ and $c_{u}$ are the specific heats per unit mass under constant pressure and constant volume, respectively. Then we have the following theorem.

Theorem 3.1. If matrices $A$ with $tr(A) = 0$ and $B = (A_{t}+A^{2}+JA+\alpha(t) A)/2$ satisfy the matrix differential equations

$\begin{equation} B^{T} = B\mathit{\text{,}} \end{equation}$

(3.2)

$\begin{equation} B_{t}+BA+A^{T}B = 0, \end{equation}$

(3.3)

then the compressible Euler equations with a time-dependent linear damping and Coriolis force (1.1)–(1.3) have explicit solutions in the form

$\begin{equation} \mathbf{u} = \mathbf{b}(t)+A\mathbf{x}, \end{equation}$

(3.4)

$\begin{equation} \rho = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)]^{\frac{1}{\gamma}}, \end{equation}$

(3.5)

$\begin{equation} S = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)], \end{equation}$

(3.6)

where $\mu = (\frac{\gamma}{\gamma+1})^{\frac{1}{\gamma}}$ ; the vector function ${\mathbf{b}}(t)$ and scalar function $c(t)$ satisfy the ordinary differential equations:

$\begin{equation} (\mathbf{b}_{t}+A\mathbf{b}+J\mathbf{b}+\alpha(t)\mathbf{b})_t+A^{T}(\mathbf{b}_{t}+A\mathbf{b}+J\mathbf{b} +\alpha(t)\mathbf{b})+2B\mathbf{b} = \mathbf{0}, \end{equation}$

(3.7)

$\begin{equation} c_{t}-\mathbf{b}^{T}(\mathbf{b}_{t}+A\mathbf{b}+J\mathbf{b}+\alpha(t)\mathbf{b}) = 0. \end{equation}$

(3.8)

Proof. By (3.65), (3.5) and (3.6), $\rho > 0$ , $S = \ln\rho$ . Let

$\begin{equation} \bar{p} = \frac{\gamma+1}{\gamma}\rho^\gamma, \end{equation}$

(3.9)

$\begin{equation} \frac{\triangledown p}{\rho} = \frac{1}{\rho}\triangledown(e^s\rho^\gamma) = \frac{1}{\rho}\triangledown(\rho^{\gamma+1}) = (\gamma+1)\rho^{\gamma-1}\triangledown\rho = \triangledown(\frac{\gamma+1}{\gamma}\rho^\gamma) = \triangledown\bar{p}. \end{equation}$

(3.10)

With (3.9), the compressible Euler equations (1.2) and (1.3) can then be written as

$\begin{equation} \rho_{t}+\mathrm{div}(\rho{\mathbf{u}}) = 0, \end{equation}$

(3.11)

$\begin{equation} {\mathbf{u}}_{t}+({\mathbf{u}}\cdot\triangledown ){\mathbf{u}}+J{\mathbf{u}}+\alpha(t){\mathbf{u}}+\triangledown\bar{p} = \mathbf{0}, \end{equation}$

(3.12)

$\begin{equation} S_t+{\mathbf{u}}\cdot\triangledown S = 0. \end{equation}$

(3.13)

Owing to the equivalent relation (3.9) between $\bar{p}$ and $\rho$ , we mainly deal with $\bar{p}$ when solving Eqs (3.11) and (3.12). Substituting Eq (3.4) into Eq (3.12), we have

$\begin{align} & {\mathbf{u}}_{t}+({\mathbf{u}}\cdot\triangledown ){\mathbf{u}}+J{\mathbf{u}}+\alpha(t) {\mathbf{u}}+\triangledown\bar{p} \end{align}$

(3.14)

$\begin{align} & = {\mathbf{b}}_{t}+A_{t}{\mathbf{x}}+[({\mathbf{b} }+A{\mathbf{x}})\cdot\triangledown]({\mathbf{b}} +A{\mathbf{x}})+JA{\mathbf{x}}+\alpha(t) A{\mathbf{x}}+J{\mathbf{b} }+\alpha(t) {\mathbf{b}}+\triangledown\bar{p} \end{align}$

(3.15)

$\begin{align} & = {\mathbf{b}}_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A_{t}{\mathbf{x} }+({\mathbf{b}}\cdot\triangledown)A{\mathbf{x}} +(A{\mathbf{x}}\cdot\triangledown)A{\mathbf{x}} +JA{\mathbf{x}}+\alpha(t) A{\mathbf{x}}+\triangledown\bar{p} \end{align}$

(3.16)

$\begin{align} & = {\mathbf{b}}_{t}+(A+J+\alpha(t)){\mathbf{b}}+(A_{t}+A^{2} +JA+\alpha(t) A){\mathbf{x}}+\triangledown\bar{p} = \mathbf{0}. \end{align}$

(3.17)

Let

$\begin{equation} B = (b_{ij})_{N\times N} = \frac{1}{2}(A_{t}+A^{2}+JA+\alpha(t) A),{ \ \ \ \ \ } J = (g_{ij})_{N\times N}. \end{equation}$

(3.18)

Then the above equation can be written into a component form

$\begin{equation} Q_{i}(x_{1},\cdots,x_{N})\equiv-b_{it}-\alpha(t) b_i-\sum\limits_{k = 1}^{N}(a_{ik}b_{k}+g_{ik} b_{k}+2b_{ik}x_{k}) = \frac{\partial\bar{p}}{\partial x_{i}}\text{,} \ \ \ \ \ \ i = 1,2,\cdots,N. \end{equation}$

(3.19)

Then, the following sufficient and necessary compatible conditions of these $N$ equations,

$\begin{equation} \frac{\partial Q_{j}(x_{1},\cdots,x_{N})}{\partial x_{i}} = \frac{\partial Q_{i}(x_{1},\cdots,x_{N})}{\partial x_{j}},\ \ \ \ \ \ i,j = 1,2,\cdots ,N, \end{equation}$

(3.20)

lead to

$\begin{equation} b_{ji} = b_{ij},\ \ \ \ \ \ i,j = 1,2,\cdots,N, \end{equation}$

(3.21)

which implies that $B = \frac{1}{2}(A_{t}+A^{2}+JA+\alpha(t) A)$ is a symmetric matrix. Under the condition (3.20), $\bar{p}({\mathbf{x}})$ is a complete differential function,

$\begin{equation} d\bar{p}({\mathbf{x}}) = \sum\limits_{i = 1}^{N}\frac{\partial\bar{p} ({\mathbf{x}})}{\partial x_{i}}dx_{i} = \sum\limits_{i = 1}^{N}Q_{i}(x_{1} ,\cdots,x_{N})dx_{i}. \end{equation}$

(3.22)

Therefore we can choose a special integration route to obtain

$\begin{align} & \bar{p}({\mathbf{x}},t) = \sum\limits_{i = 1}^{N}\int_{(0,0,\cdots,0)} ^{(x_{1},x_{2},\cdots,x_{N})}Q_{i}(x_{1},x_{2},\cdots,x_{N})dx_{i} \end{align}$

(3.23)

$\begin{align} & = \int_{0}^{x_{1}}Q_{1}(x_{1},0,\cdots,0)dx_{1},+\int_{0}^{x_{2}}Q_{2} (x_{1},x_{2},0,\cdots,0)dx_{2}\\ & +\cdots+\int_{0}^{x_{N}}Q_{N}(x_{1},x_{2},\cdots,x_{N})dx_{N} \end{align}$

(3.24)

$\begin{align} & = -\sum\limits_{i = 1}^{N}\left[ b_{i,t}+\sum\limits_{k = 1}^{N}(a_{ik}b_{k}+g_{ik} b_{k})+\alpha(t) b_i\right] x_{i}-\sum\limits_{i = 1}^{N}b_{ii}x_{i}^{2}-2\sum\limits_{i,k = 1,\ i < k} ^{N}b_{ik}x_{i}x_{k}+c(t) \end{align}$

(3.25)

$\begin{align} & = -{\mathbf{x}}^{T}({\mathbf{b}}_{t}+J{\mathbf{b} }+A{\mathbf{b}}+\alpha(t) \mathbf{b})-{\mathbf{x}}^{T}B{\mathbf{x} }+c(t). \end{align}$

(3.26)

Next, we show that functions (3.4)–(3.6) satisfy (3.11). By (3.9), we have

$\begin{align} &\rho_t = (\mu\bar{p}^{\frac{1}{\gamma}})_t = \frac{\mu}{\gamma}\bar{p}^{\frac{1}{\gamma}-1}\bar{p}_t \text{,} \end{align}$

(3.27)

$\begin{align} &\rho \mathrm{tr}(A) = \mu\bar{p}^\frac{1}{\gamma}\mathrm{tr}(A) = \frac{\mu}{\gamma}\bar{p}^{\frac{1}{\gamma}-1}\gamma \mathrm{tr}(A)\bar{p}\text{,} \end{align}$

(3.28)

$\begin{align} &\triangledown\rho = \triangledown(\mu\bar{p}^{\frac{1}{\gamma}}) = \frac{\mu}{\gamma}\bar{p}^{\frac{1}{\gamma}-1}\triangledown\bar{p}\text{,} \end{align}$

(3.29)

$\begin{align} &{\mathbf{u}}\cdot\triangledown\rho = \frac{\mu}{\gamma}\bar{p}^{\frac{1}{\gamma}-1}{\mathbf{u}}^T\triangledown\bar{p}\text{.} \end{align}$

(3.30)

From Eqs (3.27)–(3.30), we have

$\begin{align} & \rho_{t}+\mathrm{div}(\rho{\mathbf{u}}) = \rho_{t}+\rho\mathrm{tr}(A)+{\mathbf{u}}\cdot\triangledown\rho\\ = &-\frac{\mu}{\gamma} \bar{p}^{\frac{1}{\gamma}-1}\{{\mathbf{x}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})_t +{\mathbf{x}}^TB_t{\mathbf{x}}-c_t(t)\\ &+\gamma \mathrm{tr}(A)[{\mathbf{x}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})+ {\mathbf{x}}^TB{\mathbf{x}}-c(t)]\\ &+({\mathbf{b}}+A{\mathbf{x}})^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}} +2B{\mathbf{x}})\} \end{align}$

(3.31)

$\begin{align} = &-\frac{\mu}{\gamma} \bar{p}^{\frac{1}{\gamma}-1}\{{\mathbf{x}}^{T}(B_t+\gamma \mathrm{tr}(A)B+2A^TB){\mathbf{x}}\\ &+{\mathbf{x}}^{T}[({\mathbf{b}}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})_t +(\gamma \mathrm{tr}(A)I+A^T)({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})+2B{\mathbf{b}}]\\ &-[c_t+\gamma \mathrm{tr}(A)c-{\mathbf{b}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})]\} \end{align}$

(3.32)

$\begin{align} = &-\frac{\mu}{\gamma} \bar{p}^{\frac{1}{\gamma}-1}\{{\mathbf{x}}^{T}[B_t+2A^TB]{\mathbf{x}}\\ &+{\mathbf{x}}^{T}[({\mathbf{b}}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})_t+A^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})+2B{\mathbf{b}}]\\ &-[c_t-{\mathbf{b}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})]\} = 0, \end{align}$

(3.33)

where we use the condition of the first term

$\begin{equation} \mathbf{x}^{T}(B_{t}+2A^{T}B)\mathbf{x} = 0\text{,} \end{equation}$

(3.34)

which is equivalent to

$\begin{equation} (B_{t}+2A^{T}B)^{T} = -(B_{t}+2A^{T}B)\text{,} \end{equation}$

(3.35)

that is,

$\begin{equation} B_{t}+BA+A^{T}B = 0\text{,} \end{equation}$

(3.36)

which is (3.3). The second and third terms are controlled to be $0$ with (3.7) and (3.8). By (3.6), we have

$\begin{equation} S = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x}+c(t)] = \ln\rho . \end{equation}$

(3.37)

From (3.9), (3.37) is equivalent to

$\begin{align} &S = \ln(\mu\bar{p}^\gamma) = \ln\mu+\frac{1}{\gamma}\ln\bar{p} \text{,} \end{align}$

(3.38)

$\begin{align} &S_t = \frac{(\ln\bar{p})_t}{\gamma} = \frac{1}{\gamma}\bar{p}^{-1}\bar{p}_t\text{,} \end{align}$

(3.39)

$\begin{align} &\triangledown S = \frac{1}{\gamma}\triangledown\ln\bar{p} = \frac{1}{\gamma}\bar{p}^{-1}\triangledown\bar{p}\text{.} \end{align}$

(3.40)

Substituting (3.4)–(3.6) and (3.38)–(3.40) to (3.13) and using (3.3), (3.7), and (3.8), we obtain by a similar argument used in obtaining Eq (3.33) that

$\begin{align} &S_t+{\mathbf{u}}\cdot\triangledown S = \frac{1}{\gamma}\bar{p}^{-1}(\bar{p}_t+{\mathbf{u}}^T\triangledown\bar{p}) \end{align}$

(3.41)

$\begin{align} = & \frac{1}{\gamma}\bar{p}^{-1}[-{\mathbf{x}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}} +\alpha(t){\mathbf{b}})_t-{\mathbf{x}}^TB_t{\mathbf{x}}+c_t(t)\\ &-({\mathbf{x}}^TA^T+{\mathbf{b}}^T)({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}} +2B{\mathbf{x}})] \end{align}$

(3.42)

$\begin{align} = &\frac{1}{\gamma\bar{p}}\{-{\mathbf{x}}^T[({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}} +\alpha(t){\mathbf{b}})_t+A^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}})+2B{\mathbf{b}}]\\ &-{\mathbf{x}}^T[B_t+2A^TB]{\mathbf{x}}+c_t(t)-{\mathbf{b}}^T({\mathbf{b}}_t+A{\mathbf{b}}+J{\mathbf{b}} +\alpha(t){\mathbf{b}})\} = 0. \end{align}$

(3.43)

□

We observe that Eq (3.3) is a $N^{2}$ matrix differential equation, which demands us to apply special reduction conditions to acquire solutions.

Corollary 3.1. If $A$ is an anti-symmetric matrix, that is

$\begin{equation} A^T = -A, \end{equation}$

(3.44)

and the following conditions are satisfied:

$\begin{equation} A_t+\alpha(t) A = 0\mathit{\text{,}} \end{equation}$

(3.45)

$\begin{equation} AJ = JA\mathit{\text{,}} \end{equation}$

(3.46)

$\begin{equation} B_t = 0\mathit{\text{,}} \end{equation}$

(3.47)

$\begin{equation} {\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} = \mathbf{0}\mathit{\text{,}} \end{equation}$

(3.48)

$\begin{equation} c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) = 0\mathit{\text{,}} \end{equation}$

(3.49)

then the compressible Euler equations (3.11)–(3.13) admit a general solution

$\begin{align} & \mathbf{u} = \mathbf{b}(t)+A\mathbf{x}, \end{align}$

(3.50)

$\begin{align} & \rho = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)]^{\frac{1}{\gamma}}, \end{align}$

(3.51)

$\begin{align} & S = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x}+c(t)]. \end{align}$

(3.52)

Proof. By (3.45) and (3.46),

$\begin{align} B^T & = \frac{1}{2}(A_t+A^2+JA+\alpha(t) A)^T \end{align}$

(3.53)

$\begin{align} & = \frac{1}{2}[(-A)(-A)+(-A)(-J)] \end{align}$

(3.54)

$\begin{align} & = \frac{1}{2}(A^2+JA) = B. \end{align}$

(3.55)

We can then simplify (3.3), (3.7), and (3.8) into (3.47), (3.48), and (3.49). Since matrix $A$ is anti-symmetric, we have

$\begin{equation} BA+A^{T}B = 0. \end{equation}$

(3.56)

By (3.47), we have

$\begin{equation} B_{t} = 0, \end{equation}$

(3.57)

$\begin{equation} B_{t}+BA+A^{T}B = 0. \end{equation}$

(3.58)

Thus, Eq (3.3) is ensured.

Since

$\begin{align} &B^T = B, \end{align}$

(3.59)

$\begin{align} &AJ = JA, \end{align}$

(3.60)

$\begin{align} &A^T+A = 0, \end{align}$

(3.61)

we have

$\begin{align} &(\mathbf{b}_{t}+A\mathbf{b}+J\mathbf{b}+\alpha(t)\mathbf{b})_t+A^{T}(\mathbf{b}_{t}+A\mathbf{b}+J\mathbf{b} +\alpha(t)\mathbf{b})+2B\mathbf{b} \end{align}$

(3.62)

$\begin{align} & = {\mathbf{b}}_{tt}+A_t{\mathbf{b}}+A{\mathbf{b}}_t-A{\mathbf{b}}_t-A(A\mathbf{b}+J\mathbf{b}+\alpha(t)\mathbf{b})\\ &+(A_t+A^2+JA+\alpha(t)A){\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.63)

$\begin{align} & = {\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} = {\mathbf{0}}. \end{align}$

(3.64)

Thus, Eq (3.64) is simplified to (3.48). □

Next, we give the following examples in $2$ to $N$ -dimension to demonstrate special cases of this corollary.

Remark 3.1. As (3.5) and (3.6) demand

$\begin{equation} -{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t) > 0 \end{equation}$

(3.65)

for the positivity of the argument of the logarithm and density, the solutions exist locally.

Example 3.1. When $\alpha = 0$ , we have the following examples:

2-dimensional Case: We take constant matrix

$\begin{equation} A = J = k_1\left[ \begin{array} [c]{cc} 0 & 1 \\ -1 & 0 \end{array} \right] {, \ } \mathbf{b} = k_2\left[ \begin{array} [c]{c} \cos(k_1t)\\ \sin(k_1t) \end{array}\right]{, \ }c(t) = 0\text{, } \end{equation}$

(3.66)

where $k_1$ and $k_2$ are arbitrary constants.

By (3.18),

$\begin{align} B& = \frac{1}{2}(A_{t}+A^{2}+JA+\alpha(t) A) = \frac{1}{2}(2A^2) = A^2 \end{align}$

(3.67)

$\begin{align} & = {k_1}^2\left[ \begin{array} [c]{cc} -1 & 0 \\ 0 & -1 \end{array} \right]. \end{align}$

(3.68)

Since $A$ is a constant matrix, $A_t = 0$ , taking $\alpha(t) = 0$ ,

$\begin{equation} B_t = \frac{d(A_t+A^2+JA+\alpha(t)A)}{2dt} = 0, \end{equation}$

(3.69)

Eqs (3.45) and (3.47) are satisfied. As $J\! = \!A$ , Eq (3.46) is guaranteed. Equation (3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.70)

$\begin{align} & = -{k_1}^2{\mathbf{b}}+{\mathbf{0}}+J{\mathbf{b}}_{t}+{\mathbf{0}} \end{align}$

(3.71)

$\begin{align} & = -{k_1}^2k_2\left[ \begin{array} [c]{c} \cos(k_1t)\\ \sin(k_1t) \end{array}\right]+{k_1}^2k_2\left[ \begin{array} [c]{cc} 0 & 1 \\ -1 & 0 \end{array} \right]\left[ \begin{array} [c]{c} -\sin(k_1t)\\ \cos(k_1t) \end{array}\right] = {\mathbf{0}}\text{, } \end{align}$

(3.72)

Eq (3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.73)

$\begin{align} & = 0-k_2\left[ \begin{array} [c]{c} \cos(k_1t)\\ \sin(k_1t) \end{array}\right]^T\left(k_1k_2\left[ \begin{array} [c]{c} -\sin(k_1t)\\ \cos(k_1t) \end{array}\right]+2k_1k_2\left[ \begin{array} [c]{cc} 0 & 1 \\ -1 & 0 \end{array} \right]\left[ \begin{array} [c]{c} \cos(k_1t)\\ \sin(k_1t) \end{array}\right]\right) \end{align}$

(3.74)

$\begin{align} & = -k_1{k_2}^2\left[ \begin{array} [c]{c} \cos(k_1t)\\ \sin(k_1t) \end{array}\right]^T\left[ \begin{array} [c]{c} \sin(k_1t)\\ -\cos(k_1t) \end{array}\right] = 0 \text{.} \end{align}$

(3.75)

we obtain the following solution:

$\begin{equation} \mathbf{u}(t) = \left[ \begin{array} [c]{c} k_2\cos(k_1t)+ k_1x_2 \\ k_2\sin(k_1t)-k_1x_1 \end{array} \right] , \end{equation}$

(3.76)

$\begin{align} \rho& = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)]^{\frac{1}{\gamma}} \end{align}$

(3.77)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+2A{\mathbf{b}})-\mathbf{x}^{T}A^2\mathbf{x}]^{\frac{1}{\gamma}} \end{align}$

(3.78)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}(k_1k_2\left[ \begin{array} [c]{c} -\sin(k_1t) \\ \cos(k_1t) \end{array} \right]+2k_1k_2\left[ \begin{array} [c]{c} -\sin(k_1t) \\ \cos(k_1t) \end{array} \right])-\mathbf{x}^{T}{k_1}^2\left[ \begin{array} [c]{cc} -1 & 0 \\ 0 & -1 \end{array} \right]\mathbf{x}]^{\frac{1}{\gamma}} \end{align}$

(3.79)

$\begin{align} & = \mu [{k_1}^2({x_1}^2+{x_2}^2)+k_1k_2(-\sin(k_1t)x_1+\cos(k_1t)x_2)]^\frac{1}{\gamma} , \end{align}$

(3.80)

$\begin{align} S& = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x}+c(t)] \end{align}$

(3.81)

$\begin{align} & = \ln\mu+\frac{1}{\gamma}\ln [{k_1}^2({x_1}^2+{x_2}^2)+k_1k_2(-\sin(k_1t)x_1+\cos(k_1t)x_2)] . \end{align}$

(3.82)

3-dimensional Case: We take constant matrix

$\begin{equation} A = J = k_1\left[ \begin{array} [c]{ccc} 0 & 1 & -1\\ -1 & 0 & 1\\ 1 & -1 & 0 \end{array} \right] {, \ } \mathbf{b} = k_2t\left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right]{, \ }c(t) = \frac{3{k_2}^2}{2}t^2\text{, } \end{equation}$

(3.83)

where $k_1$ and $k_2$ are arbitrary constants.

Since matrix $A$ is a constant matrix, (3.45)–(3.47) are satisfied. By using of (3.83), (3.48) and (3.49) are ensured. By (3.18),

$\begin{align} B& = \frac{1}{2}(A_{t}+A^{2}+JA+\alpha(t) A) = \frac{1}{2}(2A^2) = A^2 \end{align}$

(3.84)

$\begin{align} & = {k_1}^2\left[ \begin{array} [c]{ccc} -2 & 1 & 1\\ 1 & -2 & 1\\ 1 & 1 & -2 \end{array} \right]. \end{align}$

(3.85)

Since $A$ is a constant matrix, $A_t = 0$ , taking $\alpha(t) = 0$ ,

$\begin{equation} B_t = \frac{d(A_t+A^2+JA+\alpha(t)A)}{2dt} = 0, \end{equation}$

(3.86)

Eqs (3.45) and (3.47) are satisfied. As $J = A$ , Eq (3.46) is guaranteed. Equation (3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.87)

$\begin{align} & = {\mathbf{0}}+{\mathbf{0}}+J{\mathbf{b}}_{t}+{\mathbf{0}} \end{align}$

(3.88)

$\begin{align} & = k_1k_2\left[ \begin{array} [c]{ccc} 0 & 1&-1 \\ -1 & 0&1\\ 1&-1&0 \end{array} \right]\left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right] = {\mathbf{0}}\text{, } \end{align}$

(3.89)

Eq (3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.90)

$\begin{align} & = 3{k_2}^2t-k_2\left[ \begin{array} [c]{c} t\\ t\\ t \end{array}\right]^T\left(k_2\left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right]+2k_1k_2\left[ \begin{array} [c]{ccc} 0 & 1&-1 \\ -1 & 0&1\\ 1&-1&0 \end{array} \right]\left[ \begin{array} [c]{c} t\\ t\\ t \end{array}\right]\right) \end{align}$

(3.91)

$\begin{align} & = 3{k_2}^2t-3{k_2}^2t+0 = 0 \text{.} \end{align}$

(3.92)

Therefore we obtain the solution:

$\begin{equation} \mathbf{u}(t) = \left[ \begin{array} [c]{c} k_2 t+ k_1(x_2-x_3) \\ k_2 t+ k_1(x_3-x_1)\\ k_2 t+ k_1(x_1-x_2) \end{array} \right] , \end{equation}$

(3.93)

$\begin{align} \rho& = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)]^{\frac{1}{\gamma}} \end{align}$

(3.94)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+2A{\mathbf{b}})-\mathbf{x}^{T}A^2\mathbf{x}+\frac{3{k_2}^2}{2}t^2]^{\frac{1}{\gamma}} \end{align}$

(3.95)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}(k_2\left[ \begin{array} [c]{c} 1 \\ 1\\ 1 \end{array} \right]+0)-\mathbf{x}^{T}{k_1}^2\left[ \begin{array} [c]{ccc} -2 & 1 & 1\\ 1 & -2 & 1\\ 1 & 1 & -2 \end{array} \right]\mathbf{x}+\frac{3{k_2}^2}{2}t^2]^{\frac{1}{\gamma}} \end{align}$

(3.96)

$\begin{align} & = \mu [2{k_1}^2({x_1}^2+{x_2}^2+{x_3}^2-x_1x_2-x_1x_3-x_2x_3)-k_2(x_1+x_2+x_3)+\frac{3{k_2}^2}{2}t^2]^\frac{1}{\gamma} , \end{align}$

(3.97)

$\begin{align} S& = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x}+c(t)] \end{align}$

(3.98)

$\begin{align} & = \ln\mu+\frac{1}{\gamma}\ln [2{k_1}^2({x_1}^2+{x_2}^2+{x_3}^2-x_1x_2-x_1x_3-x_2x_3)-k_2(x_1+x_2+x_3)+\frac{3{k_2}^2}{2}t^2] . \end{align}$

(3.99)

Remark 3.2. The 3-dimensional example has the same setting with Example 5 in ^[37], which admits the same ${\mathbf{u}}$ solution but has different entropy and density.

4-dimensional Case: We take

$\begin{equation} A = J = k_1\left[ \begin{array} [c]{cccc} 0 & -2 & 1& 1 \\ 2 & 0 & 1& -3\\ -1 & -1 & 0& 2 \\ -1& 3 & -2&0 \end{array} \right], \end{equation}$

(3.100)

$\begin{equation} \mathbf{b} = k_2t\left[ \begin{array} [c]{c} 1\\ 1\\ 1\\ 1 \end{array}\right]{, \ }c(t) = 2{k_2}^2t^2\text{, } \end{equation}$

(3.101)

where $k_1$ and $k_2$ are arbitrary constants. By (3.18),

$\begin{align} B& = \frac{1}{2}(A_{t}+A^{2}+JA+\alpha(t) A) = \frac{1}{2}(2A^2) = A^2 \end{align}$

(3.102)

$\begin{align} & = {k_1}^2\left[ \begin{array} [c]{cccc} -6& 2& -4& 8 \\ 2& -14& 8& 4\\ -4& 8& -6& 2 \\ 8& 4& 2&- 14 \end{array} \right]. \end{align}$

(3.103)

Since $A$ is a constant matrix, $A_t = 0$ , taking $\alpha(t) = 0$ ,

$\begin{equation} B_t = \frac{d(A_t+A^2+JA+\alpha(t)A)}{2dt} = 0, \end{equation}$

(3.104)

Eqs (3.45) and (3.47) are satisfied. As $J = A$ , Eq (3.46) is guaranteed. Equation (3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.105)

$\begin{align} & = {\mathbf{0}}+{\mathbf{0}}+J{\mathbf{b}}_{t}+{\mathbf{0}} \end{align}$

(3.106)

$\begin{align} & = k_1k_2\left[ \begin{array} [c]{cccc} 0 & -2 & 1& 1 \\ 2 & 0 & 1& -3\\ -1 & -1 & 0& 2 \\ -1& 3 & -2&0 \end{array} \right]\left[ \begin{array} [c]{c} 1\\ 1\\ 1\\ 1 \end{array}\right] = {\mathbf{0}}\text{, } \end{align}$

(3.107)

Eq (3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.108)

$\begin{align} & = 4{k_2}^2t-k_2\left[ \begin{array} [c]{c} t\\ t\\ t\\ t \end{array}\right]^T\left(k_2\left[ \begin{array} [c]{c} 1\\ 1\\ 1\\ 1 \end{array}\right]+2k_1k_2\left[ \begin{array} [c]{cccc} 0 & -2 & 1& 1 \\ 2 & 0 & 1& -3\\ -1 & -1 & 0& 2 \\ -1& 3 & -2&0 \end{array} \right]\left[ \begin{array} [c]{c} t\\ t\\ t\\ t \end{array}\right]\right) \end{align}$

(3.109)

$\begin{align} & = 4{k_2}^2t-4{k_2}^2t+0 = 0 \text{.} \end{align}$

(3.110)

We have the following solutions:

$\begin{equation} \mathbf{u} = k_2t\left[ \begin{array} [c]{c} 1\\ 1\\ 1\\ 1 \end{array}\right]+k_1\left[ \begin{array} [c]{c} -2x_2+x_3+x_4\\ 2x_1+x_3-3x_4\\ -x_1-x_2+2x_4\\ -x_1+3x_2-2x_3 \end{array}\right] , \end{equation}$

(3.111)

$\begin{align} \rho& = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x} +c(t)]^{\frac{1}{\gamma}} \end{align}$

(3.112)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+2A{\mathbf{b}})-\mathbf{x}^{T}A^2\mathbf{x}+2{k_2}^2t^2]^{\frac{1}{\gamma}} \end{align}$

(3.113)

$\begin{align} & = \mu[-{\mathbf{x}}^{T}(k_2\left[ \begin{array} [c]{c} 1 \\ 1\\ 1\\ 1 \end{array} \right]+0)-\mathbf{x}^{T}{k_1}^2\left[ \begin{array} [c]{cccc} -6& 2& -4& 8 \\ 2& -14& 8& 4\\ -4& 8& -6& 2 \\ 8& 4& 2&- 14 \end{array} \right]\mathbf{x}+2{k_2}^2t^2]^{\frac{1}{\gamma}} \end{align}$

(3.114)

$\begin{align} & = \mu[-k_2(x_1+x_2+x_3+x_4)+k_1^2(6{x_1}^2+14{x_2}^2+6{x_3}^2+14{x_4}^2\\ &-4x_1x_2+8x_1x_3-16x_1x_4-16x_2x_3-8x_2x_4-4x_3x_4)+2{k_2}^2t^2]^{\frac{1}{\gamma}} , \end{align}$

(3.115)

$\begin{align} S& = \ln\mu+\frac{1}{\gamma}\ln[-{\mathbf{x}}^{T}({\mathbf{b} }_{t}+J{\mathbf{b}}+\alpha(t) {\mathbf{b}}+A{\mathbf{b}})-\mathbf{x}^{T}B\mathbf{x}+c(t)] \end{align}$

(3.116)

$\begin{align} & = \ln\mu+\frac{1}{\gamma}\ln[-k_2(x_1+x_2+x_3+x_4)+k_1^2(6{x_1}^2+10{x_2}^2+6{x_3}^2+14{x_4}^2\\ &-4x_1x_2+8x_1x_3-16x_1x_4-16x_2x_3-8x_2x_4-4x_3x_4)+2{k_2}^2t^2] . \end{align}$

(3.117)

Example 3.2. When $\alpha$ is a constant, we have the following examples.

2-dimensional Case: We take

$\begin{equation} A = -J = k_1e^{-\alpha t}\left[ \begin{array} [c]{cc} 0 & 1\\ -1 & 0 \end{array} \right] {, \ } \mathbf{b} = k_2e^{-\alpha t}\left[ \begin{array} [c]{c} 1\\ 1 \end{array}\right] {, \ } \end{equation}$

(3.118)

$c(t) = m > 0$ , where $k_1$ , $k_2$ , and $m$ are arbitrary constants. Then we get a solution

$\begin{align} &\mathbf{u}(t) = e^{-\alpha t}\left[ \begin{array} [c]{c} k_1 x_2+k_2 \\ -k_1 x_1+k_2 \end{array} \right], \end{align}$

(3.119)

$\begin{align} & \rho = \mu m^{\frac{1}{\gamma}}, \end{align}$

(3.120)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln m. \end{align}$

(3.121)

3-dimensional Case: We take

$\begin{equation} A = -J = k_1e^{-\alpha t}\left[ \begin{array} [c]{ccc} 0 & 1 & 1\\ -1 & 0 & 1\\ -1 & -1 & 0 \end{array} \right] {, \ } \mathbf{b} = k_2e^{-\alpha t}\left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right], \end{equation}$

(3.122)

$c(t) = m > 0$ , where $k_1$ , $k_2$ , and $m$ are arbitrary constants. Then we get a solution

$\begin{align} &\mathbf{u}(t) = e^{-\alpha t}\left[ \begin{array} [c]{c} k_1(x_2+x_3)+k_2 \\ k_1(x_3-x_1)+k_2 \\ -k_1(x_1+x_2)+k_2 \end{array} \right], \end{align}$

(3.123)

$\begin{align} & \rho = \mu m^{\frac{1}{\gamma}}, \end{align}$

(3.124)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln m. \end{align}$

(3.125)

4-dimensional Case: We take

$\begin{equation} A = -J = k_1 e^{-\alpha t}\left[ \begin{array} [c]{cccc} 0 & 1 & 1&1\\ -1 & 0 & 1&1\\ -1 & -1 & 0 & 1\\ -1&-1&-1&0 \end{array} \right] {, \ } \mathbf{b} = k_2e^{-\alpha t}\left[ \begin{array} [c]{c} 1\\ 1\\ 1\\ 1 \end{array}\right], \end{equation}$

(3.126)

$c(t) = m > 0$ , where $k_1$ , $k_2$ , and $m$ are arbitrary constants. Then we get a solution

$\begin{align} &\mathbf{u}(t) = e^{-\alpha t}\left[ \begin{array} [c]{c} k_1(x_2+x_3+x_4)+k_2 \\ k_1(x_3+x_4-x_1)+k_2 \\ k_1(x_4-x_1-x_2)+k_2 \\ -k_1(x_1+x_2+x_3)+k_2 \end{array} \right], \end{align}$

(3.127)

$\begin{align} & \rho = \mu m^{\frac{1}{\gamma}}, \end{align}$

(3.128)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln m. \end{align}$

(3.129)

$N$ -dimensional Case: We take

$\begin{equation} A = k_1e^{-\alpha t}\left[ \begin{array} [c]{ccccc} 0 & 1 & 1& \cdots&1\\ -1 & 0 &1&\cdots&1\\ -1 & -1 & \ddots& &1\\ \vdots & \vdots & &\ddots&1\\ -1 & -1&-1 &\cdots &0 \end{array}\right] \text{, }{\mathbf{b}} = k_2 e^{-\alpha t}\left[\begin{array} [c]{c} 1 \\ 1\\ \vdots\\ 1 \end{array}\right]\text{, } \end{equation}$

(3.130)

$\begin{equation} J = -A, \ \ \ \ c(t) = m > 0 \text{, } \end{equation}$

(3.131)

where $k_1$ , $k_2$ , and $m$ are arbitrary constants. Then we get a solution

$\begin{align} & u_i = e^{-\alpha t}[k_1(\sum\limits_{k = i+1}^{N}{x_k}-\sum\limits_{k = 1}^{i-1}{x_k})+k_2], \end{align}$

(3.132)

$\begin{align} & \rho = \mu m^{\frac{1}{\gamma}}, \end{align}$

(3.133)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln m. \end{align}$

(3.134)

Proof. Since N-dimensional case covers 2 to 4-dimensional cases, here gives the verification of N-dimensional case. (3.46) is guaranteed by $J = -A$ , with

$\begin{equation} A_t = -\alpha k_1e^{-\alpha t} = -\alpha A, \end{equation}$

(3.135)

(3.45) is satisfied. Therefore,

$\begin{equation} B = \frac{1}{2}(A_t+A^2+JA+\alpha A) = 0, \ \ B_t = 0, \end{equation}$

(3.136)

(3.47) is ensured. Substituting (3.130) and (3.131) into (3.48) and (3.49) produces

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.137)

$\begin{align} & = \alpha^2 {\mathbf{b}}+2\alpha A{\mathbf{b}}+(\alpha{\mathbf{b}}-A{\mathbf{b}})_t \end{align}$

(3.138)

$\begin{align} & = \alpha^2{\mathbf{b}}+2\alpha A{\mathbf{b}}-\alpha^2{\mathbf{b}}-2\alpha A{\mathbf{b}} = {\mathbf{0}}, \end{align}$

(3.139)

and

$\begin{align} &c_t-{\mathbf{b}}^T({\mathbf{b}}_t+\alpha{\mathbf{b}}+A{\mathbf{b}}+J{\mathbf{b}}) \end{align}$

(3.140)

$\begin{align} = &0-{\mathbf{b}}^T(-\alpha{\mathbf{b}}+\alpha{\mathbf{b}}-J{\mathbf{b}}+J{\mathbf{b}}) = 0. \end{align}$

(3.141)

□ When $\alpha(t)$ is not a constant, we have the following examples.

Example 3.3. ( $2$ -dimensional case) We take

$\begin{equation} A = t^{k_1}\left[ \begin{array} [c]{cc} 0 & 1\\ -1 & 0 \end{array} \right] {, \ } J = (t^{k_1}-\frac{k_2}{t^{k_1}})\left[ \begin{array} [c]{cc} 0 & -1\\ 1 & 0 \end{array} \right] {, \ } \mathbf{b} = t^{k_1}\left[ \begin{array} [c]{c} 1\\ 1 \end{array}\right] {, \ }c(t) = \beta\text{, }\alpha(t) = -\frac{k_1}{t}\text{, } \end{equation}$

(3.142)

where $k_1 < 0$ , $k_2$ , and $\beta$ are arbitrary constants. As

$\begin{equation} AJ = JA = (t^{2k_1}-k_2)\left[ \begin{array} [c]{cc} 1 & 0\\ 0 & 1 \end{array} \right]\text{, } \end{equation}$

(3.143)

(3.46) is satisfied.

Denoting

$\begin{equation} Q = (q_{ij})_{N\times N} = \left[ \begin{array} [c]{cc} 0 & 1\\ -1 & 0 \end{array} \right]\text{, }{\mathbf{w}} = \left[ \begin{array} [c]{c} 1\\ 1 \end{array}\right]\text{, } \end{equation}$

(3.144)

it is easy to see

$\begin{equation} A_t = k_1 t^{k_1-1}Q = -\alpha(t)A\text{, } \end{equation}$

(3.145)

and,

$\begin{align} &B = \frac{A^2+JA}{2} = \frac{(A+J)A}{2} = \frac{k_2}{2t^{k_1}}QA = -\frac{k_2}{2} I = B^T\text{, } \end{align}$

(3.146)

$\begin{align} &B_t = (-\frac{k_2}{2} I)_t = 0, \end{align}$

(3.147)

therefore, (3.47) is satisfied. Since

$\begin{equation} {\mathbf{b}}_t = -\alpha(t) {\mathbf{b}}, J = \frac{k_2}{t^{k_1}}Q-A\text{, } \end{equation}$

(3.148)

(3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.149)

$\begin{align} = &-(\alpha(t) {\mathbf{b}})_t-2\alpha(t) A{\mathbf{b}}+[(\frac{k_2}{t^{k_1}}Q-A)t^{k_1}{\mathbf{w}}]_t+(\alpha(t) {\mathbf{b}})_t \end{align}$

(3.150)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-(A{\mathbf{b}})_t \end{align}$

(3.151)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-A{\mathbf{b}}_t-A_t{\mathbf{b}} \end{align}$

(3.152)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}} = {\mathbf{0}}\text{, } \end{align}$

(3.153)

(3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.154)

$\begin{align} = &0-t^{k_1}{\mathbf{w}}^T[k_1 t^{k1-1}{\mathbf{w}}+t^{k_1}Qt^{k_1}{\mathbf{w}}+(\frac{k_2}{t^{k_1}}-t^{k_1})Qt^{k_1}{\mathbf{w}}-k_1 t^{k1-1}{\mathbf{w}}] \end{align}$

(3.155)

$\begin{align} = &0-t^{k_1}{\mathbf{w}}^T(\frac{k_2}{t^{k_1}}Q t^{k_1}{\mathbf{w}}) \end{align}$

(3.156)

$\begin{align} = &-k_2t^{k_1}{\mathbf{w}}^TQ {\mathbf{w}} \end{align}$

(3.157)

$\begin{align} = &-k_2t^{k_1}\sum\limits_{i = 1,j = 1}^{N}q_{ij} = 0 \text{.} \end{align}$

(3.158)

Then we get a solution

$\begin{align} &\mathbf{u}(t) = t^{k_1}\left[ \begin{array} [c]{c} 1+x_2 \\ 1-x_1 \end{array} \right] , \end{align}$

(3.159)

$\begin{align} & \rho = \mu [k_2(\frac{{x_1}^2+{x_2}^2}{2}-x_1-x_2)+\beta]^{\frac{1}{\gamma}} , \end{align}$

(3.160)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln [k_2(\frac{{x_1}^2+{x_2}^2}{2}-x_1-x_2)+\beta] . \end{align}$

(3.161)

Example 3.4 ( $3$ -dimensional case). We take

$\begin{equation} A = t^{k_1}\left[ \begin{array} [c]{ccc} 0 & 1 & 1\\ -1 & 0 & 1\\ -1 & -1 & 0 \end{array} \right] {, \ } J = (t^{k_1}-\frac{k_2}{t^{k_1}})\left[ \begin{array} [c]{ccc} 0 & -1& -1\\ 1 & 0 & -1\\ 1 & 1 & 0 \end{array}\right] {, \ } \mathbf{b} = t^{k_1}\left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right]\text{, }c(t) = \beta\text{, } \alpha(t) = -\frac{k_1}{t}\text{, } \end{equation}$

(3.162)

where $k_1 < 0$ , $k_2$ , and $\beta$ are arbitrary constants. As

$\begin{equation} AJ = JA = (t^{2k_1}-k_2)\left[ \begin{array} [c]{ccc} 2 & 1 & -1\\ 1 & 2 &1\\ -1 & 1 & 2 \end{array} \right]\text{, } \end{equation}$

(3.163)

(3.46) is satisfied.

Denoting

$\begin{equation} Q = (q_{ij})_{N\times N} = \left[ \begin{array} [c]{ccc} 0 & 1 & 1\\ -1 & 0 &1\\ -1 & -1 & 0 \end{array} \right]\text{, }{\mathbf{w}} = \left[ \begin{array} [c]{c} 1\\ 1\\ 1 \end{array}\right]\text{, } \end{equation}$

(3.164)

it is easy to see

$\begin{equation} A_t = -\alpha(t)A\text{, } \end{equation}$

(3.165)

$\begin{equation} B = \frac{A^2+JA}{2} = \frac{(A+J)A}{2} = \frac{k_2}{2t^{k_1}}QA = \frac{k_2}{2} Q^2 = B^T\text{, } \end{equation}$

(3.166)

$\begin{equation} B = \frac{A^2+JA}{2} = \frac{k_2}{2}\left[ \begin{array} [c]{ccc} -2 & -1 & 1\\ -1 & -2 &-1\\ 1 & -1 & -2 \end{array}\right], \end{equation}$

(3.167)

therefore

$\begin{equation} B^T = B\text{, }B_t = 0\text{, } \end{equation}$

(3.168)

(3.47) is satisfied.

Since

$\begin{equation} {\mathbf{b}}_t = -\alpha(t) {\mathbf{b}}\text{, } J = \frac{k_2}{t^{k_1}}Q-A\text{, } \end{equation}$

(3.169)

(3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.170)

$\begin{align} = &-(\alpha(t) {\mathbf{b}})_t-2\alpha(t) A{\mathbf{b}}+[(\frac{k_2}{t^{k_1}}Q-A)t^{k_1}{\mathbf{w}}]_t+(\alpha(t) {\mathbf{b}})_t \end{align}$

(3.171)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-(A{\mathbf{b}})_t \end{align}$

(3.172)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-A{\mathbf{b}}_t-A_t{\mathbf{b}} \end{align}$

(3.173)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}} = {\mathbf{0}}\text{, } \end{align}$

(3.174)

(3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.175)

$\begin{align} = &0-t^{k_1}{\mathbf{w}}^T[k_1 t^{k1-1}{\mathbf{w}}+t^{k_1}Qt^{k_1}{\mathbf{w}}+(\frac{k_2}{t^{k_1}}-t^{k_1})Qt^{k_1}{\mathbf{w}}-k_1 t^{k1-1}{\mathbf{w}}] \end{align}$

(3.176)

$\begin{align} = &0-t^{k_1}{\mathbf{w}}^T(\frac{k_2}{t^{k_1}}Q t^{k_1}{\mathbf{w}}) \end{align}$

(3.177)

$\begin{align} = &-k_2t^{k_1}{\mathbf{w}}^TQ {\mathbf{w}} \end{align}$

(3.178)

$\begin{align} = &-k_2t^{k_1}\sum\limits_{i = 1,j = 1}^{N}q_{ij} = 0 \text{.} \end{align}$

(3.179)

We then get a solution

$\begin{align} &\mathbf{u}(t) = t^{k_1}\left[ \begin{array} [c]{c} x_2+x_3+1 \\ x_3-x_1+1 \\ -x_1-x_2+1 \end{array} \right] , \end{align}$

(3.180)

$\begin{align} & \rho = \mu [k_2({x_1}^2+{x_2}^2+{x_3}^2+x_1x_2+x_2x_3-x_1 x_3-2x_1+2x_3)+\beta]^{\frac{1}{\gamma}} , \end{align}$

(3.181)

$\begin{align} & S = \ln\mu +\frac{1}{\gamma}\ln [k_2({x_1}^2+{x_2}^2+{x_3}^2+x_1x_2+x_2x_3-x_1 x_3-2x_1+2x_3)+\beta] . \end{align}$

(3.182)

Remark 3.3 ( $N$ -dimensional case). We can abtain N-dimensional solutions denoting

$\begin{equation} Q = (q_{ij})_{N\times N} = \left[ \begin{array} [c]{ccccc} 0 & 1 & 1& \cdots&1\\ -1 & 0 &1&\cdots&1\\ -1 & -1 & \ddots& &1\\ \vdots & \vdots & &\ddots&1\\ -1 & -1&-1 &\cdots &0 \end{array}\right] \text{, }{\mathbf{w}} = \left[\begin{array} [c]{c} 1 \\ 1\\ \vdots\\ 1 \end{array}\right]\text{, } \end{equation}$

(3.183)

and taking

$\begin{equation} A = f(t)Q {, \ }J = (\frac{k_1}{f(t)}-f(t))Q {, \ } {\mathbf{b}} = f(t){\mathbf{w}}{, \ } \alpha(t) = -\frac{\dot{f}(t)}{f(t)}\text{, } c(t) = \beta \text{,} \end{equation}$

(3.184)

where $\dot{f}(t)f(t)\leq0$ , $k_1$ and $\beta$ are arbitrary constants. As

$\begin{equation} AJ = JA = (k_1-f(t)^2)Q^2, \end{equation}$

(3.185)

(3.46) is satisfied. It is easy to see

$\begin{equation} A_t = -\alpha(t)A\text{, } \end{equation}$

(3.186)

$\begin{equation} B = \frac{A^2+JA}{2} = \frac{(A+J)A}{2} = \frac{k_1}{2f(t)}Qf(t)Q = \frac{k_1}{2} Q^2\text{, } \end{equation}$

(3.187)

therefore

$\begin{equation} B^T = B\text{, } B_t = 0\text{, } \end{equation}$

(3.188)

(3.47) are satisfied. Since

$\begin{equation} {\mathbf{b}}_t = -\alpha(t) {\mathbf{b}}\text{, }J = \frac{k_1}{f(t)}Q-A\text{, } \end{equation}$

(3.189)

(3.48) is satisfied by

$\begin{align} &{\mathbf{b}}_{tt}+2A_t{\mathbf{b}}+(J{\mathbf{b}}+\alpha(t){\mathbf{b}})_{t} \end{align}$

(3.190)

$\begin{align} = &-(\alpha(t) {\mathbf{b}})_t-2\alpha(t) A{\mathbf{b}}+[(\frac{k_1}{f(t)}-f(t)) Q f(t){\mathbf{w}}]_t+(\alpha(t) {\mathbf{b}})_t \end{align}$

(3.191)

$\begin{align} = &-2\alpha(t)A{\mathbf{b}}+(k_1Q{\mathbf{w}}-A{\mathbf{b}})_t \end{align}$

(3.192)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-(A{\mathbf{b}})_t \end{align}$

(3.193)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}-A{\mathbf{b}}_t-A_t{\mathbf{b}} \end{align}$

(3.194)

$\begin{align} = &-2\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}}+\alpha(t) A{\mathbf{b}} = {\mathbf{0}}\text{, } \end{align}$

(3.195)

(3.49) is satisfied by

$\begin{align} &c_{t}-{\mathbf{b}}^{T}(\mathbf{b}_{t}+A{\mathbf{b}}+J{\mathbf{b}}+\alpha(t){\mathbf{b}}) \end{align}$

(3.196)

$\begin{align} = &0-f(t){\mathbf{w}}^T[\dot{f}(t){\mathbf{w}}+f(t)Q{f}(t){\mathbf{w}}+(\frac{k_1}{f(t)}-f(t))Q{f}(t){\mathbf{w}}-\dot{f}(t){\mathbf{w}}] \end{align}$

(3.197)

$\begin{align} = &0-f(t){\mathbf{w}}^T(\frac{k_1}{f(t)}Q f(t){\mathbf{w}}) \end{align}$

(3.198)

$\begin{align} = &-k_1f(t){\mathbf{w}}^TQ{\mathbf{w}} \end{align}$

(3.199)

$\begin{align} = &-k_1f(t)\sum\limits_{i = 1,j = 1}^{N}q_{ij} = 0 \text{.} \end{align}$

(3.200)

4. Conclusions and discussion

In this paper, we construct the Cartesian solutions

${\mathbf{u}} = {\mathbf{b}}(t)+A(t){\mathbf{x}}$

for the non-isentropic Euler equations with a time-dependent linear damping and a rotational force. By constructing appropriate matices $A(t)$ and vectors ${\mathbf{b}}(t)$ , we obtain new theoretical new exact solutions, which are obtained under the requirement of entropy $S = \ln\rho$ . We then invite the scientific community to provide solutions with other forms of or more general form of entropy. The global existence of the solutions remains open, while the blowup phenomena are complicated to higher dimensional cases due to the existence of many temporal variables and the multiple requirements imposed on them.

Conflict of interest

The author declares there is no interest in relation to this article.

Supplementary

Verification of examples on Euler equations

For simplicity, we use the same $\bar{p}$ defined in (3.9), solutions of $\rho$ and $S$ in all dimensions are equivalent to

$\begin{align} &\rho = \mu \bar{p}^{\frac{1}{\gamma}} , \end{align}$

(4.1)

$\begin{align} &S = \ln\mu+\frac{1}{\gamma}\ln\bar{p}. \end{align}$

(4.2)

It is clear that from the theorem (3.5) and (3.6) and can be easily verified from substitution that all solutions satisfy $S = \ln\rho$ . Dividing $\rho$ from both sides of (1.2), we rewrite the Euler equations (1.1)–(1.3) as

$\begin{align} &\rho_{t}+\sum\limits_{k = 1}^{N}\frac{\partial}{\partial x_k}\rho u_k = 0 , \end{align}$

(4.3)

$\begin{align} &\frac{\partial u_{i}}{\partial t}+\sum\limits_{k = 1}^{N}u_k(\frac{\partial u_i}{\partial x_k}+j_{ik})+\alpha u_i+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_i}\rho^\gamma = 0 , \end{align}$

(4.4)

$\begin{align} &S_t+\sum\limits_{k = 1}^{N}u_k\frac{\partial}{\partial x_k}S = 0 . \end{align}$

(4.5)

Example 1

For 2-dimension case: Substituting (3.76)–(3.82) and (4.1) into (4.3) produces

$\begin{align} &\rho_{t}+\frac{\partial}{\partial x_1}(\rho u_1)+\frac{\partial}{\partial x_2}(\rho u_2) \end{align}$

(4.6)

$\begin{align} = &\mu(\bar{p}^{\frac{1}{\gamma}})_t+\frac{\partial}{\partial x_1}(\mu\bar{p}^{\frac{1}{\gamma}} u_1)+\frac{\partial}{\partial x_2}(\mu\bar{p}^{\frac{1}{\gamma}} u_2) \end{align}$

(4.7)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\bar{p}_t+(k_2\cos(k_1t)+k_1x_2)\frac{\partial}{\partial x_1}\mu\bar{p}^\frac{1}{\gamma}+(k_2\cos(k_1t)-k_1x_1)\frac{\partial}{\partial x_2}\mu\bar{p}^\frac{1}{\gamma} \end{align}$

(4.8)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[-k_1^2k_2\cos(k_1t)x_1-k_1^2k_2\sin(k_1t)x_2]\\ &+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\frac{\partial}{\partial x_1}\left[\bar{p}(k_2\cos(k_1t)+k_1x_2)\right]+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\frac{\partial}{\partial x_2}\left[\bar{p}(k_2\cos(k_1t)-k_1x_1)\right] \end{align}$

(4.9)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\{-k_1^2k_2\cos(k_1t)x_1-k_1^2k_2\sin(k_1t)x_2\\ &+\frac{\partial}{\partial x_1}[{k_1}^2({x_1}^2+{x_2}^2)+k_1k_2(-\sin(k_1t)x_1+\cos(k_1t)x_2)](k_2\cos(k_1t)+k_1x_2)\\ &+\frac{\partial}{\partial x_2}[{k_1}^2({x_1}^2+{x_2}^2)+k_1k_2(-\sin(k_1t)x_1+\cos(k_1t)x_2)](k_2\sin(k_1t)-k_1x_1)\} \end{align}$

(4.10)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[-k_1^2k_2\cos(k_1t)x_1-k_1^2k_2\sin(k_1t)x_2\\ &+(2k_1^2x_1-k_1k_2\sin(k_1t))(k_2\cos(k_1t)+k_1x_2)\\ &+(2k_1^2x_2+k_1k_2\cos(k_1t))(k_2\sin(k_1t)-k_1x_1)] = 0. \end{align}$

(4.11)

Substituting (3.76)–(3.82) into (4.4), the first momentum gives

$\begin{align} &\frac{\partial u_{1}}{\partial t}+u_1(\frac{\partial u_1}{\partial x_1}+j_{11})+u_2(\frac{\partial u_1}{\partial x_2}+j_{12})+\alpha u_1+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_1}\rho^\gamma \end{align}$

(4.12)

$\begin{align} = &\!-k_1k_2\sin(k_1t)+u_1(0+0)+(k_2\sin(k_1t)-k_1x_1)(k_1+k_1)+2k_1^2x_1-k_1k_2\sin(k_1t) \end{align}$

(4.13)

$\begin{align} = &0, \end{align}$

(4.14)

the second momentum gives

$\begin{align} &\frac{\partial u_{2}}{\partial t}+u_1(\frac{\partial u_2}{\partial x_1}+j_{11})+u_2(\frac{\partial u_2}{\partial x_2}+j_{12})+\alpha u_2+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_2}\rho^\gamma \end{align}$

(4.15)

$\begin{align} = &k_1k_2\sin(k_1t)+(k_2\cos(k_1t)-k_1x_2)(-k_1-k_1)+u_2(0+0)+2k_1^2x_2-k_1k_2\cos(k_1t) \end{align}$

(4.16)

$\begin{align} = &0. \end{align}$

(4.17)

Substituting (3.76)–(3.82) into (4.5) gives

$\begin{align} &S_t+u_1\frac{\partial}{\partial x_1}S+u_2\frac{\partial}{\partial x_2}S \end{align}$

(4.18)

$\begin{align} = &\frac{1}{\gamma\bar{p}}(-k_1^2k_2\cos(k_1t)x_1-k_1^2k_2\sin(k_1t)x_2)\\ &+\frac{1}{\gamma\bar{p}}(2k_1^2x_1-k_1k_2\sin(k_1t))(k_2\cos(k_1t)+k_1x_2)\\ &+\frac{1}{\gamma\bar{p}}(2k_1^2x_2+k_1k_2\cos(k_1t))(k_2\sin(k_1t)-k_1x_1) \end{align}$

(4.19)

$\begin{align} = &\frac{1}{\gamma\bar{p}}(-k_1^2k_2\cos(k_1t)x_1-k_1^2k_2\sin(k_1t)x_2)\nonumber\\ &+\frac{1}{\gamma\bar{p}}(k_1^2k_2\sin(k_1t)x_2+k_1^2k_2\cos(k_1t)x_1) = 0. \end{align}$

(4.20)

For 3-dimensional case: Substituting (3.93)–(3.99) and (4.1) into (4.3) produces

$\begin{align} &\rho_{t}+\frac{\partial}{\partial x_1}(\rho u_1)+\frac{\partial}{\partial x_2}(\rho u_2)+\frac{\partial}{\partial x_3}(\rho u_3) \end{align}$

(4.21)

$\begin{align} = &\mu(\bar{p}^{\frac{1}{\gamma}})_t+\frac{\partial}{\partial x_1}(\mu\bar{p}^{\frac{1}{\gamma}} u_1)+\frac{\partial}{\partial x_2}(\mu\bar{p}^{\frac{1}{\gamma}} u_2)+\frac{\partial}{\partial x_3}(\mu\bar{p}^{\frac{1}{\gamma}} u_3) \end{align}$

(4.22)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\bar{p}_t+[k_2t+k_1(x_2-x_3)]\frac{\partial}{\partial x_1}\mu\bar{p}^\frac{1}{\gamma}\\ &+[k_2t+k_1(x_3-x_1)]\frac{\partial}{\partial x_2}\mu\bar{p}^\frac{1}{\gamma}+[k_2t+k_1(x_1-x_2)]\frac{\partial}{\partial x_3}\mu\bar{p}^\frac{1}{\gamma} \end{align}$

(4.23)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}3k_2^2t +\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(x_2-x_3)]\frac{\partial\bar{p}}{\partial x_1}\\ &+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(x_3-x_1)]\frac{\partial\bar{p}}{\partial x_2}+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(x_1-x_2)]\frac{\partial\bar{p}}{\partial x_3} \end{align}$

(4.24)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\{3k_2^2t+[k_2t+k_1(x_2-x_3)]\frac{\partial\bar{p}}{\partial x_1}\\ &+[k_2t+k_1(x_3-x_1)]\frac{\partial\bar{p}}{\partial x_2}+[k_2t+k_1(x_1-x_2)]\frac{\partial\bar{p}}{\partial x_3}\} \end{align}$

(4.25)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\{3k_2^2t+[k_2t+k_1(x_2-x_3)][2k_2^2(2x_1-x_2-x_3)-k_2]\\ &+[k_2t+k_1(x_3-x_1)][2k_2^2(2x_2-x_1-x_3)-k_2]\\ &+[k_2t+k_1(x_1-x_2)][2k_2^2(2x_3-x_1-x_2)-k_2]\} = 0. \end{align}$

(4.26)

Substituting (3.93)–(3.99) into (4.4), the first momentum gives

$\begin{align} &\frac{\partial u_{1}}{\partial t}+u_1(\frac{\partial u_1}{\partial x_1}+j_{11})+u_2(\frac{\partial u_1}{\partial x_2}+j_{12})+u_3(\frac{\partial u_1}{\partial x_3}+j_{13})+\alpha u_1+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_1}\rho^\gamma \end{align}$

(4.27)

$\begin{align} = &k_2+u_1(0+0)+[k_2t+k_1(x_3-x_1)](k_1+k_1)+[k_2t+k_1(x_1-x_2)](-k_1-k_1)\\ &+2k_1^2(2x_1-x_2-x_3)-k_2 = 0, \end{align}$

(4.28)

the second momentum gives

$\begin{align} &\frac{\partial u_{2}}{\partial t}+u_1(\frac{\partial u_2}{\partial x_1}+j_{21})+u_2(\frac{\partial u_2}{\partial x_2}+j_{22})+u_3(\frac{\partial u_2}{\partial x_3}+j_{23})+\alpha u_2+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_2}\rho^\gamma \end{align}$

(4.29)

$\begin{align} = &k_2+[k_2t+k_1(x_2-x_3)](-k_1-k_1)+u_2(0+0)+[k_2t+k_1(x_1-x_2)](k_1+k_1)\\ &+2k_1^2(2x_2-x_1-x_3)-k_2 = 0, \end{align}$

(4.30)

the third momentum gives

$\begin{align} &\frac{\partial u_{3}}{\partial t}+u_1(\frac{\partial u_3}{\partial x_1}+j_{31})+u_2(\frac{\partial u_3}{\partial x_2}+j_{32})+u_3(\frac{\partial u_3}{\partial x_3}+j_{33})+\alpha u_3+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_3}\rho^\gamma \end{align}$

(4.31)

$\begin{align} = &k_2+[k_2t+k_1(x_2-x_3)](k_1+k_1)+[k_2t+k_1(x_3-x_1)](-k_1-k_1)+u_3(0+0)\\ &+2k_1^2(2x_3-x_1-x_2)-k_2 = 0. \end{align}$

(4.32)

Substituting (3.93)–(3.99) into (4.5) gives

$\begin{align} &S_t+u_1\frac{\partial}{\partial x_1}S+u_2\frac{\partial}{\partial x_2}S+u_3\frac{\partial}{\partial x_3}S \end{align}$

(4.33)

$\begin{align} = &\frac{1}{\gamma\bar{p}}3k_2^2t+[k_2t+k_1(x_2-x_3)]\frac{\partial}{\partial x_1}\ln\bar{p}\\ &+[k_2t+k_1(x_3-x_1)]\frac{\partial}{\partial x_2}\ln\bar{p}+[k_2t+k_1(x_1-x_2)]\frac{\partial}{\partial x_3}\ln\bar{p} \end{align}$

(4.34)

$\begin{align} = &\frac{1}{\gamma\bar{p}}\{3k_2^2t+[k_2t+k_1(x_2-x_3)][2k_1^2(2x_1-x_2-x_3)-k_2]\\ &+[k_2t+k_1(x_3-x_1)][2k_1^2(2x_2-x_1-x_3)-k_2]\\ &+[k_2t+k_1(x_1-x_2)][2k_1^2(2x_3-x_1-x_2)-k_2]\} = 0. \end{align}$

(4.35)

For 4-dimensional case: Substituting (3.111)–(3.117) and (4.1) into (4.3) produces

$\begin{align} &\rho_{t}+\frac{\partial}{\partial x_1}(\rho u_1)+\frac{\partial}{\partial x_2}(\rho u_2)+\frac{\partial}{\partial x_3}(\rho u_3)+\frac{\partial}{\partial x_4}(\rho u_4) \end{align}$

(4.36)

(4.37)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\bar{p}_t+[k_2t+k_1(-2x_2+x_3+x_4)]\frac{\partial}{\partial x_1}\mu\bar{p}^\frac{1}{\gamma}+[k_2t+k_1(2x_1+x_3-3x_4)]\frac{\partial}{\partial x_2}\mu\bar{p}^\frac{1}{\gamma}\\ &+[k_2t+k_1(-x_1-x_2+2x_4)]\frac{\partial}{\partial x_3}\mu\bar{p}^\frac{1}{\gamma}+[k_2t+k_1(-x_1+3x_2-2x_3)]\frac{\partial}{\partial x_4}\mu\bar{p}^\frac{1}{\gamma} \end{align}$

(4.38)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}4k_2^2t +\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(-2x_2+x_3+x_4)]\frac{\partial\bar{p}}{\partial x_1}+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(2x_1+x_3-3x_4)]\frac{\partial\bar{p}}{\partial x_2}\\ &+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(-x_1-x_2+2x_4)]\frac{\partial\bar{p}}{\partial x_3}+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[k_2t+k_1(-x_1+3x_2-2x_3)]\frac{\partial\bar{p}}{\partial x_4} \end{align}$

(4.39)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\{4k_2^2t+[k_2t+k_1(-2x_2+x_3+x_4)]\frac{\partial\bar{p}}{\partial x_1}+[k_2t+k_1(2x_1+x_3-3x_4)]\frac{\partial\bar{p}}{\partial x_2}\\ &+[k_2t+k_1(-x_1-x_2+2x_4)]\frac{\partial\bar{p}}{\partial x_3}+[k_2t+k_1(-x_1+3x_2-2x_3)]\frac{\partial\bar{p}}{\partial x_4}\} \end{align}$

(4.40)

$\begin{align} = &\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}\{4k_2^2t+[k_2t+k_1(-2x_2+x_3+x_4)][-k_2+k_1^2(12x_1-4x_2+8x_3-16x_4)]\\ &+[k_2t+k_1(2x_1+x_3-3x_4)][-k_2+k_1^2(28x_2-4x_1-16x_3-8x_4)]\\ &+[k_2t+k_1(-x_1-x_2+2x_4)][-k_2+k_1^2(12x_3+8x_1-16x_2-4x_4)]\\ &+[k_2t+k_1(-x_1+3x_2-2x_3)][-k_2+k_1^2(28x_4-16x_1-8x_2-4x_3)]\} = 0. \end{align}$

(4.41)

Substituting (3.111)–(3.117) into (4.4), the first momentum gives

$\begin{align} &\frac{\partial u_{1}}{\partial t}+u_1(\frac{\partial u_1}{\partial x_1}+j_{11})+u_2(\frac{\partial u_1}{\partial x_2}+j_{12})+u_3(\frac{\partial u_1}{\partial x_3}+j_{13})+u_4(\frac{\partial u_1}{\partial x_4}+j_{14})+\alpha u_1+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_1}\rho^\gamma \end{align}$

(4.42)

$\begin{align} = &k_2+u_1(0+0)+[k_2t+k_1(2x_1+x_3-3x_4)](-2k_1-2k_1)\\ &+[k_2t+k_1(-x_1-x_2+2x_4)](k_1+k_1)+[k_2t+k_1(-x_1+3x_2-2x_3)](k_1+k_1)\\ &-k_2+k_1^2(12x_1-4x_2+8x_3-16x_4) = 0, \end{align}$

(4.43)

the second momentum gives

$\begin{align} &\frac{\partial u_{2}}{\partial t}+u_1(\frac{\partial u_2}{\partial x_1}+j_{21})+u_2(\frac{\partial u_2}{\partial x_2}+j_{22})+u_3(\frac{\partial u_2}{\partial x_3}+j_{23})+u_4(\frac{\partial u_2}{\partial x_4}+j_{24})+\alpha u_2+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_2}\rho^\gamma \end{align}$

(4.44)

$\begin{align} = &k_2+[k_2t+k_1(-2x_2+x_3+x_4)](2k_1+2k_1)+u_2(0+0)\\ &+[k_2t+k_1(-x_1-x_2+2x_4)](k_1+k_1)+[k_2t+k_1(-x_1+3x_2-2x_3)](-3k_1-3k_1)\\ &-k_2+k_1^2(-4x_1+28x_2-16x_3-8x_4) = 0, \end{align}$

(4.45)

the third momentum gives

$\begin{align} &\frac{\partial u_{3}}{\partial t}+u_1(\frac{\partial u_3}{\partial x_1}+j_{31})+u_2(\frac{\partial u_3}{\partial x_2}+j_{32})+u_3(\frac{\partial u_3}{\partial x_3}+j_{33})+u_4(\frac{\partial u_3}{\partial x_4}+j_{34})+\alpha u_3+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_3}\rho^\gamma \end{align}$

(4.46)

$\begin{align} = &k_2+[k_2t+k_1(-2x_2+x_3+x_4)](-k_1-k_1)+[k_2t+k_1(2x_1+x_3-3x_4)](-k_1-k_1)\\ &+u_3(0+0)+[k_2t+k_1(-x_1+3x_2-2x_3)](2k_1+2k_1)\\ &-k_2+k_1^2(8x_1-16x_2+12x_3-4x_4) = 0, \end{align}$

(4.47)

the fourth momentum gives

$\begin{align} &\frac{\partial u_{4}}{\partial t}+u_1(\frac{\partial u_4}{\partial x_1}+j_{41})+u_2(\frac{\partial u_4}{\partial x_2}+j_{42})+u_3(\frac{\partial u_4}{\partial x_3}+j_{43})+u_4(\frac{\partial u_4}{\partial x_4}+j_{44})+\alpha u_4+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_4}\rho^\gamma \end{align}$

(4.48)

$\begin{align} = &k_2+[k_2t+k_1(-2x_2+x_3+x_4)](-k_1-k_1)\\ &+[k_2t+k_1(2x_1+x_3-3x_4)](3k_1+3k_1)+[k_2t+k_1(-x_1-x_2+2x_4)](-2k_1-2k_1)\\ &+u_4(0+0)-k_2+k_1^2(-16x_1-8x_2-4x_3+28x_4) = 0. \end{align}$

(4.49)

Substituting (3.111)–(3.117) into (4.5) gives

$\begin{align} &S_t+u_1\frac{\partial}{\partial x_1}S+u_2\frac{\partial}{\partial x_2}S+u_3\frac{\partial}{\partial x_3}S+u_4\frac{\partial}{\partial x_4}S \end{align}$

(4.50)

$\begin{align} = &\frac{\bar{p}_t}{\gamma\bar{p}}+u_1\frac{\partial}{\partial x_1}\frac{\ln\bar{p}}{\gamma}+u_2\frac{\partial}{\partial x_2}\frac{\ln\bar{p}}{\gamma}+u_3\frac{\partial}{\partial x_3}\frac{\ln\bar{p}}{\gamma}+u_4\frac{\partial}{\partial x_4}\frac{\ln\bar{p}}{\gamma} \end{align}$

(4.51)

$\begin{align} = &\frac{1}{\gamma\bar{p}}\{4k_2^2t+[k_2t+k_1(-2x_2+x_3+x_4)][-k_2+k_1^2(12x_1-4x_2+8x_3-16x_4)]\\ &+[k_2t+k_1(2x_1+x_3-3x_4)][-k_2+k_1^2(28x_2-4x_1-16x_3-8x_4)]\\ &+[k_2t+k_1(-x_1-x_2+2x_4)][-k_2+k_1^2(12x_3+8x_1-16x_2-4x_4)]\\ &+[k_2t+k_1(-x_1+3x_2-2x_3)][-k_2+k_1^2(28x_4-16x_1-8x_2-4x_3)]\} = 0. \end{align}$

(4.52)

Example 2

Since N-dimensional case covers 2 to 4-dimensional cases, here gives the verification of N-dimensional case. Substituting solutions into Euler equations, as $S$ is a constant, (4.5) is guaranteed. Since $\rho$ is also a constant, by

$\begin{align} &\rho_{t}+\sum\limits_{k = 1}^{N}\frac{\partial}{\partial x_k}\rho u_k \end{align}$

(4.53)

$\begin{align} = &0+\rho\sum\limits_{k = 1}^{N}\frac{\partial}{\partial x_k} u_k \end{align}$

(4.54)

$\begin{align} = &\rho e^{-\alpha t}\sum\limits_{k = 1}^{N}\frac{\partial}{\partial x_k}[k_1(\sum\limits_{g = k+1}^{N}{x_g}-\sum\limits_{g = 1}^{k-1}{x_g})+k_2] = 0, \end{align}$

(4.55)

Eq (4.3) is verified.

$\begin{align} \frac{\partial u_i}{\partial x_k}& = \frac{\partial }{\partial x_k}e^{-\alpha t}[k_1(\sum\limits_{k = i+1}^{N}{x_k}-\sum\limits_{k = 1}^{i-1}{x_k})+k_2] \end{align}$

(4.56)

$\begin{align} & = \left.\begin{cases} -k_1 e^{-\alpha t}, & \text{for } k < i \\ 0, & \text{for } k = i \\ k_1 e^{-\alpha t}, & \text{for } k > i\\ \end{cases}\right\} = -j_{ik}, \end{align}$

(4.57)

therefore,

(4.58)

$\begin{align} = &-\alpha u_i+0+\alpha u_i+0 = 0, \end{align}$

(4.59)

the $n$ -th momentum Eq (4.4) is satisfied.

Example 3

Substituting (3.159)–(3.161) into (4.3) produces

$\begin{align} &\rho_{t}+\frac{\partial}{\partial x_1}\rho u_1+\frac{\partial}{\partial x_2}\rho u_2 \end{align}$

(4.60)

$\begin{align} = &0+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}t^{k_1}(1+x_2)k_2(x_1-1)+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}t^{k_1}(1-x_1)k_2(x_2+1) = 0. \end{align}$

(4.61)

Substituting (3.159)–(3.161) into (4.4) gives

$\begin{align} &\frac{\partial u_{i}}{\partial t}+u_1(\frac{\partial u_i}{\partial x_k}+j_{i1})+u_2(\frac{\partial u_i}{\partial x_k}+j_{i2})+\alpha u_i+\frac{\gamma+1}{\gamma}\frac{\partial}{\partial x_i}\rho^\gamma \end{align}$

(4.62)

$\begin{align} = &u_1(\frac{\partial u_i}{\partial x_k}+j_{i1})+u_2(\frac{\partial u_i}{\partial x_k}+j_{i2})+0 \end{align}$

(4.63)

$\begin{align} = &k_1t^{k_1-1}\left[ \begin{array}[c]{c} 1+x_1\\ 1-x_2 \end{array} \right] +t^{k_1}(1+x_2)(t^{k_1}\left[ \begin{array}[c]{c} 0\\ -1 \end{array} \right]+\left[ \begin{array}[c]{c} 0\\ t^{k_1}-\frac{k_2}{t^{k_1}} \end{array} \right])+\\ &t^{k_1}(1-x_1)(t^{k_1}\left[ \begin{array}[c]{c} 1\\ 0 \end{array} \right]+\left[ \begin{array}[c]{c} -t^{k_1}+\frac{k_2}{t^{k_1}}\\ 0 \end{array} \right])-\frac{k_1}{t}t^{k_1}\left[ \begin{array}[c]{c} 1+x_1\\ 1-x_2 \end{array} \right]+\left[ \begin{array}[c]{c} k_2(x_1-1)\\ k_2(x_2+1) \end{array} \right] = 0. \end{align}$

(4.64)

Substituting (3.159)–(3.161) into (4.5) gives

$\begin{align} &S_t+u_1\frac{\partial}{\partial x_1}S+u_2\frac{\partial}{\partial x_2}S \end{align}$

(4.65)

$\begin{align} = &0+t_{k_1}(1+x_2)\frac{k_2(2x_1-1)}{\gamma\bar{p}}+t_{k_1}(1-x_1)\frac{k_2(2x_2+1)}{\gamma\bar{p}} = 0. \end{align}$

(4.66)

Example 4

Substituting (3.180)–(3.182) into (4.3) produces

$\begin{align} &\rho_{t}+\frac{\partial}{\partial x_1}\rho u_1+\frac{\partial}{\partial x_2}\rho u_2+\frac{\partial}{\partial x_3}\rho u_3 \end{align}$

(4.67)

$\begin{align} = &0+\frac{\mu}{\gamma}\bar{p}^{\frac{1-\gamma}{\gamma}}[t^{k_1}(x_2+x_3+1)k_2(2x_1+x_2-x_3-2)+t^{k_1}(x_3-x_1+1)k_2(2x_2+x_1+x_3)\\ &+t^{k_1}(-x_1-x_2+1)k_2(2x_3-x_1+x_2+2)] = 0. \end{align}$

(4.68)

Substituting (3.180)–(3.182) into (4.4), since

$\begin{equation} u_t = k_1t^{k_1-1} = \frac{k_1}{t}t^{k_1} = -\alpha(t)u, \end{equation}$

(4.69)

we have

(4.70)

$\begin{align} = &u_1(\frac{\partial u_i}{\partial x_k}+j_{i1})+u_2(\frac{\partial u_i}{\partial x_k}+j_{i2})+u_3(\frac{\partial u_i}{\partial x_k}+j_{i3})+0 \end{align}$

(4.71)

$\begin{align} = &t^{k_1}(x_2+x_3+1)(t^{k_1}\left[ \begin{array}[c]{c} 0\\ -1\\ -1 \end{array} \right]+\left[ \begin{array}[c]{c} 0\\ t^{k_1}-\frac{k_2}{t^{k_1}}\\ t^{k_1}-\frac{k_2}{t^{k_1}} \end{array} \right])+t^{k_1}(x_3-x_1+1)(t^{k_1}\left[ \begin{array}[c]{c} 1\\ 0\\ -1 \end{array} \right]+\left[ \begin{array}[c]{c} -t^{k_1}+\frac{k_2}{t^{k_1}}\\ 0\\ t^{k_1}-\frac{k_2}{t^{k_1}} \end{array} \right])+\\ &t^{k_1}(-x_1-x_2+1)(t^{k_1}\left[ \begin{array}[c]{c} 1\\ 1\\ 0 \end{array} \right]+\left[ \begin{array}[c]{c} -t^{k_1}+\frac{k_2}{t^{k_1}}\\ -t^{k_1}+\frac{k_2}{t^{k_1}}\\ 0 \end{array} \right]+\left[ \begin{array}[c]{c} k_2(2x_1+x_2-x_3-2)\\ k_2(2x_2+x_1+x_3)\\ k_2(2x_3-x_1+x_2-2) \end{array} \right]) = 0. \end{align}$

(4.72)

Substituting (3.180)–(3.182) into (4.5) gives

$\begin{align} &S_t+u_1\frac{\partial}{\partial x_1}S+u_2\frac{\partial}{\partial x_2}S+u_3\frac{\partial}{\partial x_3}S \end{align}$

(4.73)

$\begin{align} = &0+\frac{1}{\gamma\bar{p}}[t^{k_1}(x_2+x_3+1)k_2(2x_1+x_2-x_3-2)+t^{k_1}(x_3-x_1+1)k_2(2x_2+x_1+x_3)\\ &+t^{k_1}(-x_1-x_2+1)k_2(2x_3-x_1+x_2+2)] = 0. \end{align}$

(4.74)

References

[1]	M. Ghil, S. Childress, Topics in Geophysical Fluid Dynamics: Atmospheric Dynamics, Dynamo Theory, and Climate Dynamics, New York: Springer-Verlag, 1987. https://doi.org/10.1007/978-1-4612-1052-8
[2]	J. Pedlosky, Geophysical Fluid Dynamics, New York: Springer-Verlag, 1987. https://doi.org/10.1007/978-1-4612-4650-3
[3]	J. Marshall, R. A. Plumb, Atmosphere, Ocean, and Climate Dynamics: An Introductory Text, San Diego, CA: Academic Press, 2008.
[4]	F. V. Dolzhansky, Fundamentals of Geophysical Hydrodynamics, Berlin, Heidelberg: Springer, 2013. https://doi.org/10.1007/978-3-642-31034-8
[5]	H. Liu, E. Tadmor, Rotation prevents finite-time breakdown, Phys. D, 188 (2004), 262–276. https://doi.org/10.1016/j.physd.2003.07.006 doi: 10.1016/j.physd.2003.07.006
[6]	B. Cheng, E. Tadmor, Long-time existence of smooth solutions for the rapidly rotating shallow-water and Euler equations, SIAM J. Math. Anal., 39 (2008), 1668–1685. https://doi.org/10.1137/070693643 doi: 10.1137/070693643
[7]	B. Cheng, C. J. Xie, On the classical solutions of two dimensional inviscid rotating shallow water system, J. Differ. Equations, 250 (2011), 690–709. https://doi.org/10.1016/j.jde.2010.09.017 doi: 10.1016/j.jde.2010.09.017
[8]	O. S. Rozanova, M. K. Turzynsky, Nonlinear stability of localized and non-localized vortices in rotating compressible media, In: Theory, Numerics and Applications of Hyperbolic Problems II, Cham: Springer, 2018. https://doi.org/10.1007/978-3-319-91548-7_41
[9]	O. S. Rozanova, J. L. Yu, C. K. Hu, Typhoon eye trajectory based on a mathematical model: Comparing with observational data, Nonlinear Anal. Real World Appl., 11 (2010), 1847–1861. https://doi.org/10.1016/j.nonrwa.2009.04.011 doi: 10.1016/j.nonrwa.2009.04.011
[10]	O. S. Rozanova, J. L. Yu, C. K. Hu, On the position of vortex in a two-dimensional model of atmosphere, Nonlinear Anal. Real World Appl., 13 (2012), 1941–1954. https://doi.org/10.1016/j.nonrwa.2011.12.023 doi: 10.1016/j.nonrwa.2011.12.023
[11]	O. S. Rozanova, M. K. Turzynsky, On systems of nonlinear ODE arising in gas dynamics: Application to vortical motion, In: Differential and Difference Equations with Applications, Cham: Springer, 2018. https://doi.org/10.1007/978-3-319-75647-9_32
[12]	O. S. Rozanova, M. K. Turzynsky, Full classification of motions with uniform deformation on a rotating plane, AIP Conf. Proc. 2164 (2019), 090005. https://doi.org/10.1063/1.5130835
[13]	O. S. Rozanova, M. K. Turzynsky, The stability of vortices in gas on the l-plane: The influence of centrifugal force, In: Nonlinear Analysis and Boundary Value Problems, Cham: Springer, 2019. https://doi.org/10.1007/978-3-030-26987-6_9
[14]	L. Hsiao, T. P. Liu, Convergence to nonlinear diffusion waves for solutions of a system of hyperbolic conservation laws with damping, Comm. Math. Phys., 143 (1992), 599–605. https://doi.org/10.1007/BF02099268 doi: 10.1007/BF02099268
[15]	D. Y. Fang, J. Xu, Existence and asymptotic behavior of C1 solutions to the multi-dimensional compressible Euler equations with damping, Nonlinear Anal. Theory Methods Appl., 70 (2009), 244–261. https://doi.org/10.1016/j.na.2007.11.049 doi: 10.1016/j.na.2007.11.049
[16]	F. M. Huang, R. H. Pan, Convergence rate for compressible Euler equations with damping and vacuum, Arch. Ration. Mech. Anal., 166 (2003), 359–376. https://doi.org/10.1007/s00205-002-0234-5 doi: 10.1007/s00205-002-0234-5
[17]	W. K. Wang, T. Yang, The pointwise estimates of solutions for Euler equations with damping in multi-dimensions, J. Differ. Equations, 173 (2001), 410–450. https://doi.org/10.1006/jdeq.2000.3937 doi: 10.1006/jdeq.2000.3937
[18]	F. M. Huang, R. H. Pan, Asymptotic behavior of the solutions to the damped compressible Euler equations with vacuum, J. Differ. Equations, 220 (2006), 207–233. https://doi.org/10.1016/j.jde.2005.03.012 doi: 10.1016/j.jde.2005.03.012
[19]	F. M. Huang, P. Marcati, R. H. Pan, Convergence to the Barenblatt solution for the compressible Euler equations with damping and vacuum, Arch. Ration. Mech. Anal., 176 (2005), 1–24. https://doi.org/10.1007/s00205-004-0349-y doi: 10.1007/s00205-004-0349-y
[20]	K. W. Chow, E. G. Fan, M. W. Yuen, The analytical solutions for the $N$ -dimensional damped compressible Euler equations, Stud. Appl. Math., 138 (2017), 294–316. https://doi.org/10.1111/sapm.12154 doi: 10.1111/sapm.12154
[21]	J. W. Dong, J. J. Li, Analytical solutions to the compressible Euler equations with time-dependent damping and free boundaries, J. Math. Phys., 63 (2022), 101502. https://doi.org/10.1063/5.0089142 doi: 10.1063/5.0089142
[22]	S. Friedlander, M. M. Vishik, Lax pair formulation for the Euler equation, Phys. Lett. A, 148 (1990), 313–319. https://doi.org/10.1016/0375-9601(90)90809-3 doi: 10.1016/0375-9601(90)90809-3
[23]	Y. Li, A Lax pair for the two dimensional Euler equation, J. Math. Phys., 42 (2001), 3552–3553. https://doi.org/10.1063/1.1378305 doi: 10.1063/1.1378305
[24]	Y. Li, A. V. Yurov, Lax pairs and Darboux transformations for Euler equations, Stud. Appl. Math., 111 (2003), 101–113. https://doi.org/10.1111/1467-9590.t01-1-00229 doi: 10.1111/1467-9590.t01-1-00229
[25]	S. Y. Lou, M. Jia, X. Y. Tang, F. Huang, Vortices, circumfluence, symmetry groups, and Darboux transformations of the (2+1)-dimensional Euler equation, Phys. Rev. E, 75 (2007), 056318. https://doi.org/10.1103/PhysRevE.75.056318 doi: 10.1103/PhysRevE.75.056318
[26]	S. Y. Lou, M. Jia, F. Huang, X. Y. Tang, Bäcklund transformations, solitary waves, conoid waves and Bessel waves of the (2+1)-dimensional Euler equation, Internat. J. Theoret. Phys., 46 (2007), 2082–2095. https://doi.org/10.1007/s10773-006-9327-5 doi: 10.1007/s10773-006-9327-5
[27]	L. I. Sedov, Similarity and Dimensional Methods in Mechanics, Boca Raton: CRC Press, 1993. https://doi.org/10.1201/9780203739730
[28]	A. J. Majda, A. L. Bertozzi, Vorticity and Incompressible Flow, Cambridge: Cambridge University Press, 2001. https://doi.org/10.1017/CBO9780511613203
[29]	A. G. Ramn, Inverse Problems: Mathematical and Analytical Techniques with Applications to Engineering, Cham: Springer, 2005. https://doi.org/10.1007/b100958
[30]	P. G. Drazin, N. Riley, The Navier-Stokes Equations: A Classification of Flows and Exact Solutions, Cambridge: Cambridge University Press, 2006. https://doi.org/10.1017/CBO9780511526459
[31]	M. W. Yuen, Self-similar solutions with elliptic symmetry for the compressible Euler and Navier-Stokes equations in $R^N$ , Commun. Nonlinear Sci. Numer. Simul., 17 (2012), 4524–4528. https://doi.org/10.1016/j.cnsns.2012.05.022 doi: 10.1016/j.cnsns.2012.05.022
[32]	M. W. Yuen, Vortical and self-similar flows of 2D compressible Euler equations, Commun. Nonlinear Sci. Numer. Simul., 19 (2014), 2172–2180. https://doi.org/10.1016/j.cnsns.2013.11.008 doi: 10.1016/j.cnsns.2013.11.008
[33]	M. W. Yuen, Rotational and self-similar solutions for the compressible Euler equations in $R^3$ , Commun. Nonlinear Sci. Numer. Simul., 20 (2015), 634–640. https://doi.org/10.1016/j.cnsns.2014.06.027 doi: 10.1016/j.cnsns.2014.06.027
[34]	I. F. Barna, L. Mátyásb, Analytic solutions for the one-dimensional compressible Euler equation with heat conduction and with different kind of equations of state, Miskolc Math. Notes, 14 (2013), 785–799. https://doi.org/10.18514/MMN.2013.694 doi: 10.18514/MMN.2013.694
[35]	I. F. Barna, L. Mátyásb, Analytic solutions for the three-dimensional compressible Navier-Stokes equation, Fluid Dyn. Res., 46 (2014), 055508. https://doi.org/10.1088/0169-5983/46/5/055508 doi: 10.1088/0169-5983/46/5/055508
[36]	H. L. An, E. G. Fan, M. W. Yuen, The Cartesian vector solutions for the $N$ -dimensional compressible Euler equations, Stud. Appl. Math., 134 (2015), 101–119. https://doi.org/10.1111/sapm.12056 doi: 10.1111/sapm.12056
[37]	E. G. Fan, M. W. Yuen, A method for constructing special solutions for multidimensional generalization of Euler equations with Coriolis force, Chinese J. Phys., 72 (2021), 136–144. https://doi.org/10.1016/j.cjph.2021.03.013 doi: 10.1016/j.cjph.2021.03.013
[38]	T. T. Li, Y. Zhou, D. X. Kong, Weak linear degeneracy and global classical solutions for general quasilinear hyperbolic systems, Commun. Partial Differ. Equ., 19 (1994), 1263–1317. https://doi.org/10.1080/03605309408821055 doi: 10.1080/03605309408821055
[39]	T. T. Li, Y. Zhou, D. X. Kong, Global classical solutions for general quasilinear hyperbolic systems with decay initial data, Nonlinear Anal., 28 (1997), 1299–1332. https://doi.org/10.1016/0362-546X(95)00228-N doi: 10.1016/0362-546X(95)00228-N
[40]	G. Chen, R. H. Pan, S. G. Zhu, Singularity formation for the compressible Euler equations, SIAM J. Math. Anal., 49 (2017), 2591–2614. https://doi.org/10.1137/16M1062818 doi: 10.1137/16M1062818

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(1688) PDF downloads(59) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

Cartesian vector solutions for $N$ -dimensional non-isentropic Euler equations with Coriolis force and linear damping

Related Papers:

Abstract

1. Introduction

2. Materials and methods

3. Results

4. Conclusions and discussion

Conflict of interest

Supplementary

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Catalog

AIMS Mathematics

Cartesian vector solutions for N N -dimensional non-isentropic Euler equations with Coriolis force and linear damping

Related Papers:

Abstract

1. Introduction

2. Materials and methods

3. Results

4. Conclusions and discussion

Conflict of interest

Supplementary

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Cartesian vector solutions for $N$ -dimensional non-isentropic Euler equations with Coriolis force and linear damping