Pure Mathematics
Vol. 11  No. 09 ( 2021 ), Article ID: 45348 , 8 pages
10.12677/PM.2021.119184

简单波和二维非线性波动系统的特征分解

冯恬

长安大学理学院,陕西 西安

收稿日期:2021年8月13日;录用日期:2021年9月15日;发布日期:2021年9月22日

摘要

本文考虑具有广义Chaplygin气体的二维非线性波动系统,对可压缩欧拉系统进行简化,得到了自相似平面上压力变量和倾斜角变量沿着特征线的简单波和一些特征分解,这些结果是Courant和Friedrichs专著中关于可约方程的著名定理的推广。

关键词

非线性波动系统,特征分解,简单波

Simple Waves and Characteristic Decompositions of the Two-Dimensional Nonlinear Wave System

Tian Feng

School of Science, Chang’an University, Xi’an Shaanxi

Received: Aug. 13th, 2021; accepted: Sep. 15th, 2021; published: Sep. 22nd, 2021

ABSTRACT

In this paper, we consider a two-dimensional nonlinear wave system with generalized Chaplygin gas which is simplified model of the compressible Euler system and get simple waves and some characteristic decompositions for the pressure variable and inclination angle variables along the characteristic curves in the self-similar plane. These results are generalization of the well-known theorem on reducible equations in Courant and Friedrichs’s monograph.

Keywords:Nonlinear Wave Systems, Characteristic Decomposition, Simple Waves

Copyright © 2021 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

黎曼不变量的存在在双曲守恒律系统解的构造中起着基础性的作用,例如波动方程中D’Alembert公式的构造,以及奇点的发展证明 [1]。对于二维定常等熵无旋欧拉方程组,由于黎曼不变量的存在,任何与定常状态相邻的解都是一个简单波 [2]。黎曼不变量帮助我们得出任何与常状态相邻的双曲状态的解都是简单波的结论,然而,由于黎曼不变量对于n > 2可能一般不存在,这种方法有局限性。

Dai和Zhang [3] 首次为压力梯度系统揭示了特征分解是作为构建全局平滑解斑块的强大工具。Li、Zhang和Zheng [4] 使用特征分解技术研究了准线性双曲系统,证明了一个非常重要的结果,绝热欧拉系统在自相似平面上相邻的一个恒定状态是一个简单波,其中物理变量 ( u , v , c , p , ρ ) 沿着波的直线特征族是恒定的。Hu和Sheng [5] 建立了广义2阶拟线性严格双曲方程组特征分解存在的一个充分条件。通过使用特征分解技术,我们不仅发现了黎曼变量,而且还发现了在一些双曲系统中所谓的Riemann变量。因此特征分解方法是处理拟线性双曲系统问题的一种有效方法 [4] [6] [7] [8] [9]。

本文考虑二维非线性波动系统

{ ρ t + ( ρ u ) x + ( ρ v ) y = 0 , ( ρ u ) t + p x = 0 , ( ρ v ) t + p y = 0 , (1.1)

其中 ρ 为密度, ( u , v ) 为速度,状态方程 p = p ( ρ ) ρ 的赋值函数,系统(1.1)是通过忽略气体动力学中二维可压缩Euler方程速度 ( u , v ) 的二次项而得到的,Čanić、Keyfitz、Kim [10] 分别推导了非线性波动系统和可压缩气体动力学绝热方程的混合系统,并证明了所产生的混合系统具有复杂的非线性相关性。Hu和Wang [11] 提出了压力的特征分解,并利用分解建立先验估计,构造了带有Chaplygin气体的系统(1.1)的半双曲片的全局解。本文我们推导了压力变量和角变量的一些有趣的特征分解。

在本文中,我们感兴趣的是修正的Chaplygin气体

p = A ρ B ρ α , (1.2)

其中参数 α > 0 , A > 0 , B > 0 ,它是用来描述当前宇宙的加速膨胀 [12]。注意到当 A = 0 , α = 1 时,(1.2)由Chaplygin [13] 和Tsien [14] 引入,作为空气动力学中计算飞机机翼升力的合适数学近似。另一方面,对于 A = 0 , 0 < α < 1 ,该模型给出了从初始尘埃类物质到渐进宇宙常数的宇宙学演化。这一广义模型以前已经被研究 [15]。

为了进一步研究系统(1.1)的内在特征,我们的主要系统是(1.1)在自相似变量 ( ξ , η ) = ( x / t , y / t ) 下的系统

{ ξ ρ ξ + η ρ η ( ρ u ) ξ ( ρ v ) η = 0 , ξ ( ρ u ) ξ + η ( ρ u ) η p ξ = 0 , ξ ( ρ v ) ξ + η ( ρ v ) η p η = 0. (1.3)

我们写(1.3)为矩阵形式

L W ξ + Q W η = 0 , (1.4)

其中

a ( p ) = 1 ρ ( p ) = A + B α ρ α 1 , W = ( p , u , v ) , L = ( u ξ a ρ 0 ξ u a a ρ ξ 0 ξ v 0 a ρ ξ ) Q = ( v η 0 a ρ η u a ρ η 0 η v a 0 a ρ η ) .

通过简单计算,我们得到(1.4)的特征方程

(1.5)

并且特征值为

λ 0 = η ξ , λ ± = ξ η ± a ( ξ 2 + η 2 a ) ξ 2 a . (1.6)

系统(1.4)在 ξ 2 + η 2 = a 处改变了类型,且当 ξ 2 + η 2 > a 时为双曲型。我们注意到系统(1.1)的特征值与 ( u , v ) 无关,这意味着(1.1)可以转化为一个仅关于p的偏微分方程。

2. 预备知识

我们首先将p从 ( ρ u ) ( ρ v ) 中解耦得到一个二阶拟线性方程

( a ( p ) ξ 2 ) p ξ ξ 2 ξ η p ξ η + ( a ( p ) η 2 ) p η η + b ( p ) ( ξ p ξ + η p η ) 2 2 ( ξ p ξ + η p η ) = 0 , (2.1)

其中 b ( p ) = ρ ( p ) ρ ( p ) = a ( p ) a ( p ) ,(2.1)的特征值是

Λ ± = λ ± = ξ η ± a ( ξ 2 + η 2 a ) ξ 2 a . (2.2)

2.1. 四个角变量

在这里,我们使用了倾角变量 α , β [8] ,它首次提出是在Courant和Friedrichs的专著中 [2] ,同时流角可以被定义为 σ : = ( α + β ) / 2 ,并且马赫角 ω : = ( α β ) / 2 ,我们可以用这些变量表示 ( ξ , η )

ξ = a cos σ sin ω , η = a sin σ sin ω . (2.3)

注意到这与 [4] 的结果一致,现在我们再次推导细节。因为 tan α = Λ + , tan β = Λ

tan α = ξ η + a ( ξ 2 + η 2 a ) ξ 2 a , tan β = ξ η a ( ξ 2 + η 2 a ) ξ 2 a . (2.4)

对上述公式进行相减和相乘,我们有

tan α tan β = 2 a ( ξ 2 + η 2 a ) ξ 2 a , tan α × tan β = η 2 a ξ 2 a . (2.5)

为了消除 η 2 ,我们注意到 η 2 a = tan α × tan β × ( ξ 2 a ) ,并且得到

( tan α tan β ) ( ξ 2 a ) = 2 a [ ξ 2 + tan α × tan β × ( ξ 2 a ) ] .

然后我们有

ξ a = ± cos α ± cos β sin ( α β ) . (2.6)

虽然每种情况都是有意义的,但是为了方便,我们选择

ξ a = cos α + cos β sin ( α β ) . (2.7)

最后将 α β 的表达式代入(2.7),我们得到(2.3)的第一个方程,应用 ( tan α + tan β ) ( ξ 2 a ) = 2 ξ η ,我们获得(2.3)的第二个方程。

2.2. 一阶特征形式

首先,我们引入一些基本符号

+ = ξ + Λ + η , = ξ + Λ η , 0 = ξ ξ η η , ¯ + = cos α ξ + sin α η , ¯ = cos β ξ + sin β η , ¯ 0 = cos σ ξ + sin σ η , ¯ = sin σ ξ cos σ η . (2.8)

因此,我们有

¯ + = cos α + , ¯ = cos β , ¯ 0 = 0 ξ 2 + η 2 , ¯ + + ¯ = 2 cos ω ¯ 0 , ¯ + ¯ = 2 sin ω ¯ . (2.9)

性质2.1 对于非线性波动系统(1.1),有下式成立

¯ ± ω = sin 2 ω a + b tan ω 2 ¯ ± p , ¯ ± σ = sin 2 ω a , ¯ α = 2 sin 2 ω a + b tan ω 2 ¯ p , ¯ + α = b tan ω 2 ¯ + p , ¯ + β = 2 sin 2 ω a b tan ω 2 ¯ + p , ¯ β = b tan ω 2 ¯ p , ¯ σ = sin ω a , ¯ 0 σ = 0. (2.10)

证明:证明的方法以前已经介绍 [7] ,因此对我们来说很容易证明这些公式,我们首先证明 ¯ ± σ ¯ ± ω ,微分公式(2.3),我们有

¯ ± ξ = a ¯ ± p 2 a cos σ sin ω + a sin σ sin ω ¯ ± σ + cos σ cos ω ¯ ± ω sin 2 ω , (2.11)

¯ ± η = a ¯ ± p 2 a sin σ sin ω a cos σ sin ω ¯ ± σ sin σ cos ω ¯ ± ω sin 2 ω . (2.12)

结合(2.11)和(2.12),我们得到

cos ( α σ ) = a ¯ + p 2 a 1 sin ω + a cos ω sin 2 ω ¯ + ω , cos ( β σ ) = a ¯ p 2 a 1 sin ω + a cos ω sin 2 ω ¯ ω .

通过计算发现

¯ ± ω = sin 2 ω a + b tan ω 2 ¯ ± p . (2.13)

同理,我们有

¯ ± σ = sin 2 ω a . (2.14)

由于 α = σ + ω , β = σ ω ,所以 ¯ ± α , ¯ ± β 的结果也可以获得,根据(2.9),我们可以直接得到余下公式的证明。

3. 特征分解

这一节,我们将给出变量 p , α , β 的二阶特征分解,其中p的分解是最重要的。

定理3.1 压力p在非线性波动系统中满足

{ ¯ ¯ + p = + p ( b 4 cos 2 ω + p + ( b 2 + b 4 cos 2 ω ) p + sin 2 ω a ) , ¯ + ¯ p = p ( b 4 cos 2 ω p + ( b 2 + b 4 cos 2 ω ) + p + sin 2 ω a ) . (3.1)

证明:回顾(2.8)中 ¯ ± 的定义,我们获得

p ξ = sin α ¯ p sin β ¯ + p sin 2 ω , p η = cos β ¯ + p cos α ¯ p sin 2 ω . (3.2)

¯ ¯ + p = ( c o s β ξ + s i n β η ) ( c o s α ξ + s i n α η ) p = ( s i n α p ξ ¯ α + cos α p η ¯ α ) + ( c o s α ¯ p ξ + sin α ¯ p η ) . (3.3)

将(3.2)代入(3.3),我们有

s i n α p ξ ¯ α + cos α p η ¯ α = b tan ω 2 sin 2 ω ( ¯ p ) 2 + b tan ω 2 tan 2 ω ¯ p ¯ + p + 2 sin 2 ω a tan 2 ω ¯ + p 2 sin 2 ω a sin 2 ω ¯ p . (3.4)

c o s α ¯ p ξ + sin α ¯ p η = cos α cos β ( p ξ ξ + ( t a n α + t a n β ) p ξ η + t a n α t a n β p η η ) = cos α cos β a ξ 2 ( ( a ξ 2 ) p ξ ξ 2 ξ η p ξ η + ( a η 2 ) p η η ) = b 4 cos 2 ω ( ¯ + p + ¯ p ) 2 + 2 sin 2 ω a sin 2 ω ( ¯ + p + ¯ p ) . (3.5)

然后我们将(3.4)和(3.5)代入(3.3)得到

¯ ¯ + p = b 4 cos 2 ω ( ¯ + p ) 2 + ( b tan ω 2 tan 2 ω + b 2 cos 2 ω ) ¯ p ¯ + p + 2 sin 2 ω ( cos 2 ω + 1 ) a sin 2 ω ¯ + p ¯ + p + ( 2 sin 2 ω a sin 2 ω + 2 sin 2 ω a sin 2 ω ) ¯ p + ( b tan ω 2 sin 2 ω + b 4 cos 2 ω ) ( ¯ p ) 2 = b 4 cos 2 ω ( ¯ + p ) 2 + ( b 2 + b 4 cos 2 ω ) ¯ p ¯ + p + sin 2 ω a ¯ + p = ¯ + p ( b 4 cos 2 ω ¯ + p + ( b 2 + b 4 cos 2 ω ) ¯ p + sin 2 ω a ) . (3.6)

¯ + ¯ p 类似可得。

我们现在返回来考虑系统(1.3)的第二个和第三个方程,它们可以写成

{ 0 ( ρ u ) + p ξ = 0 , 0 ( ρ v ) + p η = 0. (3.7)

通过一系列推导,我们可以计算出 ( u , v )

进一步,我们使用定理3.1可以获得 α , β 的二阶特征分解。

定理3.2 对于变量 α β ,我们有以下二阶特征分解

{ ¯ ¯ + α + M 1 ¯ + α = 0 , ¯ + ¯ β + M 2 ¯ β = 0 , ¯ + ¯ α + M 3 ¯ α = ( 2 b sin 2 ω a b 2 3 tan ω a ) ¯ + α , ¯ ¯ + β + M 4 ¯ + β = ( 2 b sin 2 ω a b 2 3 tan ω a ) ¯ β . (3.8)

M 1 = 1 sin 2 ω ¯ + α ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ α + 2 b sin 2 ω a b 2 + 2 tan ω a , M 2 = 1 sin 2 ω ¯ β + ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ + β 2 b sin 2 ω a b 2 2 tan ω a , M 3 = 1 sin 2 ω ¯ α ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ + α + tan ω a sin 2 ω a , M 4 = 1 sin 2 ω ¯ + β + ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ β + tan ω a sin 2 ω a .

证明:由性质2.1,我们有

{ ¯ + p = 2 b tan ω ¯ + α = 2 b tan ω ¯ + β 2 sin 2 ω a b , ¯ p = 2 b tan ω ¯ α 2 sin 2 ω a b = 2 b tan ω ¯ β . (3.9)

由定理3.1,性质2.1和(3.9),我们可以计算 ¯ ¯ + α

¯ ¯ + α = ¯ ( b tan ω 2 ¯ + p ) = 1 2 [ ( b ¯ p tan ω + b cos 2 ω ¯ ω ) ¯ + p + b tan ω ¯ ¯ + p ] = ¯ + α [ 1 sin 2 ω ¯ + α + ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ α 2 b sin 2 ω a b 2 2 tan ω a ] . (3.10)

类似地我们有

¯ + ¯ β = ¯ β [ 1 sin 2 ω ¯ β ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ + β + 2 b sin 2 ω a b 2 + 2 tan ω a ] . (3.11)

以相同的方式可以计算 ¯ + ¯ α

¯ + ¯ α = ¯ + ( 2 sin 2 ω a + b tan ω 2 ¯ p ) = 2 sin 2 ω a ¯ + ω a sin 2 ω a a ¯ + p + b tan ω 2 ¯ + p ¯ p + b 2 cos 2 ω ¯ + ω ¯ p + b tan ω 2 ¯ + ¯ p = ¯ α [ 1 sin 2 ω ¯ α + ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ + α tan ω a + sin 2 ω a ] + ¯ + α ( 2 b sin 2 ω a b 2 3 tan ω a ) .

同理可得

¯ ¯ + β = ¯ + β [ 1 sin 2 ω ¯ + β ( 2 b b 2 tan ω + cos 2 ω + 4 sin 2 ω ) ¯ β tan ω a + sin 2 ω a ] + ¯ β ( 2 b sin 2 ω a b 2 3 tan ω a ) .

定理3.2证明完成。

同时,特征分解使我们得出以下结论。

定理3.3 在二维非线性波动系统的自相似平面上,任何与常状态相邻的双曲状态都是一个简单波区域。

文章引用

冯 恬. 简单波和二维非线性波动系统的特征分解
Simple Waves and Characteristic Decompositions of the Two-Dimensional Nonlinear Wave System[J]. 理论数学, 2021, 11(09): 1657-1664. https://doi.org/10.12677/PM.2021.119184

参考文献

  1. 1. Lax, P. (1964) Development of Singularities of Solutions of Nonlinear Hyperbolic Partial Differential Equations. Journal of Mathematical Physics, 5, 611-613. https://doi.org/10.1063/1.1704154

  2. 2. Courant, R. and Friedrichs, K.O. (1948) Supersonic Flow and Shock Waves. Interscience, New York.

  3. 3. Dai, Z.H. and Zhang, T. (2000) Existence of Global Smooth Solution for a Degenerate Goursat Problem of Gas Dynamics. Archive for Rational Mechanics and Analysis, 155, 277-298. https://doi.org/10.1007/s002050000113

  4. 4. Li, J.Q., Zhang, T. and Zheng, Y.X. (2006) Simple Waves and a Characteristic Decomposition of the Two Dimensional Compressible Euler Equations. Communications in Mathematical Physics, 267, 1-12. https://doi.org/10.1007/s00220-006-0033-1

  5. 5. Hu, Y.B. and Sheng, W.C. (2012) Characteristic Decomposition of the 2× 2 Quasilinear Strictly Hyperbolic Systems. Applied Mathematics Letters, 25, 262-267.

  6. 6. Chen, X. and Zheng, Y.X. (2010) The Interaction of Rare Faction Waves of the Two-Dimensional Euler Equations. Indiana University Mathematics Journal, 59, 231-256. https://doi.org/10.1512/iumj.2010.59.3752

  7. 7. Li, J.Q., Yang, Z.C. and Zheng, Y.X. (2011) Characteristic Decompositions and Interactions of Rare Faction Waves of 2-D Euler Equations. Journal of Differential Equations, 250, 782-798. https://doi.org/10.1016/j.jde.2010.07.009

  8. 8. Li, J.Q. and Zheng, Y.X. (2009) Interaction of Rare Faction Waves of the Two-Dimensional Self-Similar Euler Equations. Archive for Rational Mechanics and Analysis, 193, 623-657. https://doi.org/10.1007/s00205-008-0140-6

  9. 9. Zafar, M. and Sharma, V.D. (2014) Characteristic Decomposition of Compressible Euler Equations for a Non-Ideal Gas in Two-Dimensions. Journal of Mathematical Physics, 55, 1-12. https://doi.org/10.1063/1.4896080

  10. 10. Čanić, S., Keyfitz, B.L. and Kim, E.H. (2001) Mixed Hyperbolic-Elliptic Systems in Self-Similar Flows. Boletim da Sociedade Brasileira de Matemática, 32, 377-399. https://doi.org/10.1007/BF01233673

  11. 11. Hu, Y.B. and Wang, G.D. (2014) Semi-Hyperbolic Patches of Solutions to the Two-Dimensional Nonlinear Wave System for Chaplygin Gases. Journal of Differential Equations, 257, 1567-1590. https://doi.org/10.1016/j.jde.2014.05.020

  12. 12. Benaoum, H.B. (2002) Accelerated Universe from Modified Chaplygin Gas and Tachyonic Fluid.

  13. 13. Chaplygin, S. (1904) On Gas Jets. Scientific Memoirs, Moscow University Mathematic Physics, 21, 1-121.

  14. 14. Tsien, H.S. (1939) Two Dimensional Subsonic Flow of Compressible Fluids. Journal of the Aeronautical Sciences, 6, 399-407. https://doi.org/10.2514/8.916

  15. 15. Bento, M.C., Bertolami, O. and Sen, A.A. (2002) Generalized Chaplygin Gas, Accelerated Expansion, and Dark-Energy- Matter Unification. Physical Review D, 66, Article ID: 043507. https://doi.org/10.1103/PhysRevD.66.043507

期刊菜单