Recruitment strategies for cervical cancer screening in three Mediterranean low and middle-income countries: Albania, Montenegro, and Morocco

Elisa Camussi; Lina Jaramillo; Roberta Castagno; Marta Dotti; Gianluigi Ferrante; Latifa Belakhel; Youssef Chami Khazraji; Alban Ylli; Kozeta Filipi; Đjurđjica Ostojić; Milica Stanisic; Luigi Bisanti; Livia Giordano; Elisa Camussi; Lina Jaramillo; Roberta Castagno; Marta Dotti; Gianluigi Ferrante; Latifa Belakhel; Youssef Chami Khazraji; Alban Ylli; Kozeta Filipi; Đjurđjica Ostojić; Milica Stanisic; Luigi Bisanti; Livia Giordano

doi:10.3934/medsci.2024009

AIMS Medical Science

2024, Volume 11, Issue 2: 99-112. doi: 10.3934/medsci.2024009

Previous Article Next Article

Research article Topical Sections

Recruitment strategies for cervical cancer screening in three Mediterranean low and middle-income countries: Albania, Montenegro, and Morocco

1.
Epidemiology and Screening Department – CPO Piedmont, AOU Città della Salute e della Scienza di Torino - Via Cavour 31, 10123 Turin, Italy
2.
Non-communicable diseases unit, Epidemiology and Disease Control Department, Ministry of Health, 335 Ave Mohammed V, Rabat 10020, Morocco
3.
Lalla Salma Foundation, Prevention and Treatment of Cancer, Villa No. 1 Touarga Fouaka Mechouar Said 10700 Rabat, Morocco
4.
Institute of Public Health of Albania, Rruga Aleksander Moisiu, nr 80, Tirane, Albania
5.
Institute for Public Health of Montenegro, bb Džona Džeksona, Podgorica, Montenegro
6.
Epidemiologia & Prevenzione Scientific Board, via Giusti 4, 21053 Castellanza (VA), Italy

Received: 11 March 2024 Revised: 15 May 2024 Accepted: 23 May 2024 Published: 31 May 2024

Introduction

Cervical cancer (CC) poses a substantial burden in low-and middle-income countries (LMICs), where challenges in implementing effective screening programs and achieving high participation rates persist.

Aims

This study sought to compare different strategies for recruiting women for CC screening in Albania, Montenegro, and Morocco, and compared usual care (ongoing invitation method) with an alternative approach (intervention strategy).

Methods

Within each country, the following comparisons were made: face-to-face (FF) invitations versus phone calls (PCs) in Albania, PCs versus letter invitations in Montenegro, and FF invitations to women attending healthcare centers versus a combined approach termed “Invitation made in Morocco” (utilizing PC and FF for hard-to-reach women) in Morocco. Questionnaires that assessed facilitators and barriers to participation were administered to women who either attended or refused screening.

Results

In Albania, significant differences in the examination coverage were observed between the invitation methods (PC: 46.1% vs. FF: 87.1%, p < 0.01) and between the rural and urban settings (rural: 89.1% vs. urban: 76.3%, p < 0.01). In Montenegro, the coverage varied based on the recruitment method (PC: 17.7% vs. letter invitation: 7.6%; p < 0.01), the setting (urban: 28.3% vs. rural: 13.2%; p < 0.01), and age (<34 years: 10.9% vs. 34+: 9.6%, p < 0.01). In Morocco, no significant differences were observed. Common screening facilitators included awareness of CC prevention and understanding the benefits of early diagnosis, while key barriers included a limited perception of personal CC risk and the fear of testing positive.

Discussion

FF appeared to be effective in promoting participation, but its broader implementation raised sustainability concerns. PC invitations proved feasible, albeit necessitating updates to population registries. Restricting FF contacts for hard-to-reach communities may enhance the affordability and equity.

Keywords:

Citation: Elisa Camussi, Lina Jaramillo, Roberta Castagno, Marta Dotti, Gianluigi Ferrante, Latifa Belakhel, Youssef Chami Khazraji, Alban Ylli, Kozeta Filipi, Đjurđjica Ostojić, Milica Stanisic, Luigi Bisanti, Livia Giordano. Recruitment strategies for cervical cancer screening in three Mediterranean low and middle-income countries: Albania, Montenegro, and Morocco[J]. AIMS Medical Science, 2024, 11(2): 99-112. doi: 10.3934/medsci.2024009

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

Introduction

Aims

Methods

Results

Discussion

Abbreviations

CC:

Cervical cancer;

EuMedCN:

European Mediterranean Cancer Network;

FF:

Face-to-face;

HICs:

High-income countries;

HPV:

Human Papilloma Virus;

IMM:

Invitation made in Morocco;

LMICs:

Low-and middle-income countries;

PC:

Phone call;

WHO:

World Health Organization

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)

Author contributions

Camussi Elisa: statistical analysis, writing, review & editing; Jaramillo Lina: study design, statistical analysis, writing; Castagno Roberta: writing, review & editing; Dotti Marta: writing; Ferrante Gianluigi: writing, review & editing; Belakhel Latifa: study design, writing, review & editing; Chami Khazraji Youssef: study design, writing, review & editing; Ylli Alban: study design, writing, review & editing; Filipi Kozeta: study design, writing, review & editing; Ostojić Đjurđjica: study design, writing, review & editing; Stanisic Milica: study design, writing, review & editing; Bisanti Luigi: writing, review & editing; Giordano Livia: study design, validation, writing, review & editing, supervision. All authors have read and approved the final version of the manuscript for publication.

Conflict of interest

The authors report there are no conflicting interests to declare.

References

[1]	Bray F, Laversanne M, Sung H, et al. (2024) Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 74: 229-263. https://doi.org/10.3322/caac.21834
[2]	Hull R, Mbele M, Makhafola T, et al. (2020) Cervical cancer in low and middle-income countries. Oncol Lett 20: 2058-2074. https://doi.org/10.3892/ol.2020.11754
[3]	Ebrahimi N, Yousefi Z, Khosravi G, et al. (2023) Human papillomavirus vaccination in low- and middle-income countries: progression, barriers, and future prospective. Front Immunol 14: 1150238. https://doi.org/10.3389/fimmu.2023.1150238
[4]	Duncan J, Harris M, Skyers N, et al. (2021) A call for low-and middle-income countries to commit to the elimination of cervical cancer. Lancet Reg Health Am 2. https://doi.org/10.1016/j.lana.2021.100036
[5]	Bruni L, Serrano B, Roura E, et al. (2022) Cervical cancer screening programmes and age-specific coverage estimates for 202 countries and territories worldwide: a review and synthetic analysis. Lancet Glob Health 10: e1115-e1127. https://doi.org/10.1016/S2214-109X(22)00241-8
[6]	Srinath A, van Merode F, Rao SV, et al. (2023) Barriers to cervical cancer and breast cancer screening uptake in low- and middle-income countries: a systematic review. Health Policy Plan 38: 509-527. https://doi.org/10.1093/heapol/czac104
[7]	Devarapalli P, Labani S, Nagarjuna N, et al. (2018) Barriers affecting uptake of cervical cancer screening in low and middle income countries: a systematic review. Indian J Cancer 55: 318-326. https://doi.org/10.4103/ijc.IJC_253_18
[8]	Petersen Z, Jaca A, Ginindza TG, et al. (2022) Barriers to uptake of cervical cancer screening services in low-and-middle-income countries: a systematic review. BMC Womens Health 22: 486. https://doi.org/10.1186/s12905-022-02043-y
[9]	World Health OrganizationGlobal strategy to accelerate the elimination of cervical cancer as a public health problem (2020).
[10]	Giordano L, Bisanti L, Salamina G, et al. (2016) The EUROMED CANCER network: state-of-art of cancer screening programmes in non-EU Mediterranean countries. Eur J Public Health 26: 83-89. https://doi.org/10.1093/eurpub/ckv107
[11]	Union for the MediterraneanThe WoRTH project: women's right to health. Available from: https://ufmsecretariat.org/project/womens-right-to-health-the-worth-project
[12]	Tavasoli SM, Pefoyo AJ, Hader J, et al. (2016) Impact of invitation and reminder letters on cervical cancer screening participation rates in an organized screening program. Prev Med 88: 230-236. https://doi.org/10.1016/j.ypmed.2016.04.019
[13]	Camilloni L, Ferroni E, Cendales BJ, et al. (2013) Methods to increase participation in organised screening programs: a systematic review. BMC Public Health 13: 1-16. https://doi.org/10.1186/1471-2458-13-464
[14]	Staley H, Shiraz A, Shreeve N, et al. (2021) Interventions targeted at women to encourage the uptake of cervical screening. Cochrane Database Syst Rev 9: CD002834. https://doi.org/10.1002/14651858.CD002834.pub3
[15]	Antinyan A, Bertoni M, Corazzini L (2021) Cervical cancer screening invitations in low and middle income countries: evidence from Armenia. Soc Sci Med 273: 113739. https://doi.org/10.1016/j.socscimed.2021.113739
[16]	Devi S, Joshi S (2023) The effect of multimodal interventions regarding early cervical cancer diagnosis on the women's knowledge, attitude and participation in cervical screening program. Asian Pac J Cancer Prev 24: 3949-3956. https://doi.org/10.31557/APJCP.2023.24.11.3949
[17]	Decker KM, Turner D, Demers AA, et al. (2013) Evaluating the effectiveness of cervical cancer screening invitation letters. J Womens Health 22: 687-693. https://doi.org/10.1089/jwh.2012.4203
[18]	Rees I, Jones D, Chen H, et al. (2018) Interventions to improve the uptake of cervical cancer screening among lower socioeconomic groups: a systematic review. Prev Med J 111: 323-335. https://doi.org/10.1016/j.ypmed.2017.11.019
[19]	Vale DB, Teixeira JC, Bragança JF, et al. (2021) Elimination of cervical cancer in low- and middle-income countries: inequality of access and fragile healthcare systems. Int J Gynaecol Obstet 152: 7-11. https://doi.org/10.1002/ijgo.13458
[20]	Orwat J, Caputo N, Key W, et al. (2017) Comparing rural and urban cervical and breast cancer screening rates in a privately insured population. Soc Work Public Health 32: 311-323. https://doi.org/10.1080/19371918.2017.1289872
[21]	Diendéré J, Kiemtoré S, Coulibaly A, et al. (2023) Low attendance in cervical cancer screening, geographical disparities and sociodemographic determinants of screening uptake among adult women in Burkina Faso: results from the first nationwide population-based survey. Rev Epidemiol Sante 71: 101845. https://doi.org/10.1016/j.respe.2023.101845
[22]	Donovan J, O'Donovan C, Nagraj S (2019) The role of community health workers in cervical cancer screening in low-income and middle-income countries: a systematic scoping review of the literature. BMJ Glob Health 4: e001452. https://doi.org/10.1136/bmjgh-2019-001452
[23]	Darj E, Chalise P, Shakya S (2019) Barriers and facilitators to cervical cancer screening in Nepal: a qualitative study. Sex Reprod Healthc 20: 20-26. https://doi.org/10.1016/j.srhc.2019.02.001
[24]	Black E, Hyslop F, Richmond R (2019) Barriers and facilitators to uptake of cervical cancer screening among women in Uganda: a systematic review. BMC Womens Health 19: 108. https://doi.org/10.1186/s12905-019-0809-z
[25]	Liebermann EJ, VanDevanter N, Shirazian T, et al. (2020) Barriers to cervical cancer screening and treatment in the Dominican Republic: perspectives of focus group participants in the Santo Domingo area. J Transcult Nurs 31: 121-127. https://doi.org/10.1177/1043659619846247
[26]	Zhang M, Sit JWH, Chan DNS, et al. (2022) Educational interventions to promote cervical cancer screening among rural populations: a systematic review. Int J Environ Res Public Health 19: 6874. https://doi.org/10.3390/ijerph19116874
[27]	Zhang D, Advani S, Waller J, et al. (2020) Mobile technologies and cervical cancer screening in low- and middle-income countries: a systematic review. JCO Glob Oncol 6: 617-627. https://doi.org/10.1200/JGO.19.00201
[28]	Uy C, Lopez J, Trinh-Shevrin C, et al. (2017) Text messaging interventions on cancer screening rates: a systematic review. J Med Internet Res 19: e296. https://doi.org/10.2196/jmir.7893

Reader Comments

Your name:*

Email:*
© 2024 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 Medical Science

0.4

Metrics

Article views(1320) PDF downloads(45) Cited by(1)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Tables(3)

AIMS Medical Science

Recruitment strategies for cervical cancer screening in three Mediterranean low and middle-income countries: Albania, Montenegro, and Morocco

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

Figures and Tables

Other Articles By Authors

Catalog

AIMS Medical Science

Recruitment strategies for cervical cancer screening in three Mediterranean low and middle-income countries: Albania, Montenegro, and Morocco

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

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog