双曲型方程组
双曲型偏微分方程组的数值解法研究
双曲型偏微分方程组的数值解法研究双曲型偏微分方程组是描述波动、传播、传输等现象的常见数学模型之一，在各个科学领域中都有广泛的应用。
双曲型偏微分方程组通常具有复杂的特征，其解析解往往难以求得，因此需要用数值方法求解。
本文将介绍双曲型偏微分方程组的数值解法，并分析其优缺点，以及应用举例。
双曲型偏微分方程组的数值解法可以分为两类，即有限差分方法和有限元方法。
有限差分方法是将区域分割成网格，通过在网格上构建差分格式来近似微分方程，进而求解数值解。
有限元方法则是利用变分原理，将微分方程转化为弱形式，再通过有限元空间的数值逼近来求解数值解。
下面我们将分别介绍这两类方法。
有限差分方法是求解偏微分方程最常用的数值方法之一。
这类方法的基本思想是将区域划分成网格，通过差分逼近微分算子，将微分方程转化为代数方程组，进而求解数值解。
通常有限差分方法分为显式和隐式两种。
显式差分方法是根据精确度和稳定性的需求，选择合适的差分格式，将数值解的某一时刻的计算公式，仅由该时刻之前的数值解和已知的初值组成，计算简单，但存在较为严格的稳定性限制。
隐式差分方法则以更加严格的精确性和稳定性为代价，使用迭代法求解非线性代数方程组，计算复杂，但稳定性更加优良。
有限差分是求解双曲型偏微分方程最常见的数值方法之一。
虽然有限差分法计算公式简单，但是稳定性限制较高，当空间步长、时间步长不足以满足稳定性条件时，容易产生不稳定性及不合理的解，这是有限差分法的致命弱点之一。
此时有限元法常被作为替代方法。
有限元方法是求解双曲型偏微分方程另一种常用的数值方法。
有限元法基于变分原理，把求解微分方程转化为求最小值问题。
首先，将问题的定义域划分为若干子区域，然后在每个子区域内选取适当的试函数，通过构造一个弱变分解，就可以得到一个线性代数方程组。
有限元法具有更广泛的适用范围，解高维复杂结构问题时可以体现其独特性。
虽然有限元法可以处理不规则区域，但是计算量较大，常会出现稳定性的问题。
求双曲型偏微分方程组数值解的特征线法
特征线法是求解双曲型偏微分方程组数值解的一种常用方法，其基本思想是利用偏微分方程组的特征线方程来构造数值格式，从而求得偏微分方程组的数值解。
本文将介绍双曲型偏微分方程组的基本概念，特征线法的基本原理和数值求解过程，并结合实际问题进行案例分析，以便读者深入了解特征线法在实际工程和科学计算中的应用。
一、双曲型偏微分方程组的基本概念1. 双曲型偏微分方程组的定义和特点双曲型偏微分方程组是指具有双曲型特征的偏微分方程组，其在数学和物理上具有广泛的应用。
双曲型偏微分方程组的特点是方程中存在两个不同的特征方向，解的行为由这两个特征方向共同决定。
双曲型偏微分方程组的基本形式可以表示为：∂u/∂t + A(u)∂u/∂x = 0其中u是未知函数，t和x分别是时间和空间变量，A(u)是一个矩阵函数。
2. 双曲型偏微分方程组的物理意义和工程应用双曲型偏微分方程组描述了许多波动现象和守恒定律，因此在物理学、工程学和科学计算中有着重要的应用。
天气预报中的气象方程、弹性波动方程、流体力学方程等都可以用双曲型偏微分方程组描述，因此求解双曲型偏微分方程组的数值方法对于实际问题具有重要意义。
二、特征线法的基本原理和数值求解过程特征线法是一种求解双曲型偏微分方程组数值解的有效方法，其基本原理是利用偏微分方程组的特征线方程来构造数值格式，从而求得偏微分方程组的数值解。
特征线法的基本步骤包括确定特征线方程、构造数值格式、进行离散化和求解差分方程等。
1. 确定特征线方程双曲型偏微分方程组的特征线方程可以通过对方程进行特征分解得到，一般形式为：dx/∂t = λ1(u)du/∂t = λ2(u)其中λ1(u)和λ2(u)分别为特征线方程的两个特征方向，通过求解特征线方程可以确定数值方法的稳定性和收敛性。
2. 构造数值格式特征线法利用特征线方程构造数值格式，一般采用有限差分法或有限体积法进行离散化。
特征线法的数值格式应该满足守恒性、稳定性和收敛性等基本要求。
拟线性双曲型方程组的奇性
Ｖ ． Ｎ ｏ． ３６ ３ Ｍ ａｙ．２ ０８ ０
文 章 编 号 ： Ｏ Ｏ ３ ７ ２ ０ ） ３ ０ １ —０ １ ０ 一２ ６ （ ０ ８ ０ — ０ ０ ３
拟 线性 双 曲型方 程组 的奇 性
刘 法 贵 葛云 飞。 ，
（ ． 北 水 利 水 电学 院 数 学 系 ， 州 ４ ０ １ ；． １华 郑 ５０ １２ 河南 经 贸 职 业 学 院 ， 州 ４ ０ ５ ） 郑 ５ ０ ３
（ ２）
其 中 “一 （ ｌ … ， ， “ “ ， “ ） Ａ（ ）一 （ “ ） 口 （ ）∈ Ｃ ， （ 口 （） ， “ ， 声 ）∈ Ｃ （ ）且 Ｃ 模 有界． Ｒ
由严 格 双 曲 的定 义 ， “ Ａ（ ）存 在 个 实 而 互 异 的 特 征 （ ） （ ） … ， （ ） “ ， “ ， “ ．设 ｚ “ ）一 （ （ ） … ， （ ｚ “， ｚ “ ）（ （ ） 相应 地 ， （ ）一 （ “ ， ， “ ） 表 示 对 应 于 （ ） 一 １ ２ … ， “ ｒ （ ） … ｒ （ ） ） “（ ， ， ）的 左 （ 应 地 ， ） 征 相 右 特
］ ； ＋Ａ姑 （ ） ： Ａ２（ ） 一 ０ “ “ “ ＋ ｌ “ “＂ ，
收 稿 日期 ： ０ ７ ０ — ２ ２０— １ ９
基金项 目： 河南 省 基 础 科 学 研 究 课 题 ； 南 省 教 育 厅 自然科 学 基 金 （ ０ ８ １ ０ １ ） 河 ２ ０ Ａ １０ １ 作 者 简 介 ： 法 贵 （ ９ ， 事 偏 微 分 方 程 研 究 刘 １６ 一 ， 河 华 博 从
摘 要 ： 考虑一类 拟线性退 化特征双曲型方程组 Ｃ ｕｈ 问题 ， ａｃｙ 利用特征线法和局部延拓法证 明了其经典解 的
一阶半线性双曲型方程组的精确边界能控性与能观性
ａ ‘ ｔｘ ｃ， 亡＋Ａ（ ） ，
ｕ
Ｏ ｕ
＝
，
，
札， 一）
（） １
其 中 札：（１… ， 是 （ Ｘ 的未知 向量 函数 ， 札， 札 ） ｔ ） ，
Ａ ｔ ） ｉ ｛ １ ，）・ （ ） （ Ｘ ：ｄ ｇＡ （ ，一， ｔ ｝ ， ａ ｔ ，
Ｚ
、Ｊ
建立 了一阶半线性 双曲组 的整体精确边界能控性及能观性 理论 ，对于相应的控制及观测 时间给 出
了精 确 的 估 计 ， 并 揭示 了精 确 能 控 性 与 能 观 性 之 问 隐 含 的某 种 对 偶 关 系 。本 文 还 揭 示 了在 非 自治 系统 的情 形 下 ，精 确 能控 性 与 能 观 性 的 对 偶 关 系 有 可 能会 丧 失 。 关键 词 ： 阶 半 线性 双 曲型 方 程 组 ； 精 确 边 界 能控 性 ；精 确 边 界 能 观 性 ；对 偶性 一
ｆｔＸ 札 Ｊ Ｍ１ （， ， ）
其中 １ 为一个与 （，， ） ｔＸ 札 无关的正常数 。 设初始条件为
（） ４
ｔ ｏ： 札＝ （） ０ Ｘ Ｌ ＝ｔ ， ，
终端条件 为 ｔ ｔ ＋Ｔ： 札＝砂 ， ０ Ｘ Ｌ ＝ ｏ （）
收稿 日期： ０ ％０ —５ 作者简介：叶莎莎（９ １ １ ２ 向：偏微分方程
ｕ
，ｌ
ｌ ｌ
１
，Ｊ
，ｌ
分类号 ： ＡＭ Ｓ ２ ０ １ ５ ５ ； ７ ５ ； ９ ２ ； ３ ０ （０ ０ Ｌ ０ ３ Ｂ ５ ４ Ｊ ０ ９ Ｂ ５ ３
中图分类号： ７ ．７ Ｏ１５２
文献标识码： Ａ
、Ｊ
Ｕ
１ 引 言
考 虑 一 阶半 线 性 双 曲型万 程 组
一类带线性耗散项的2×2双曲型方程组解的衰减估计
？ １ ， Ｕ 一 （ ， ） ｒ 其 中 Ａ 一 （ Ｉ ＼ 口 ， ＇ ） Ｕ ： ， ） ， Ｂ ＝ （ ＼ Ｕ ：
ｂ ｌ
定理 １ 对 任 意 给定 的非 负整 数 ， 假 设 Ｕ０
∈Ｌ （ Ｒ） ｎＣ （ Ｒ） 且 Ｕ。有 紧支 集 ， 则 对 任 意 正整 数 Ｎ， 问题 （ １ ） 的解满 足下 列估 计
‘
Ｉ Ｉ （ Ｒ ） ， 其中Ｂ （ ｘ ｔ ） 一（ １ ＋
１ Ｔ
） ＿ 。 。 。 ， 常数 Ｃ仅 依
若记 Ｇ（ ， ￡ ）一 （ 则 上 式 可
以 简 写 为
Ｒ（ ） Ｇ
，
赖 于 Ｃ、 Ｎ 与 Ｕｏ （ ｚ ）的支集 。
２ 定 理 的 证 明
文［ １ ］ 中利用 Ｆ ｏ ｕ ｒ ｉ ｅ ｒ变 换 方 法 和 Ｇｒ ｅ ｅ ｎ函
数得 到 了带有 第 二声 速 的一维 线性 热 弹性方 程组 解 的衰 减估 计 。我 们 知 道 ， 衰 减 性 由特 征值 的 实 部 决定 ， 其 虚部 只影 响方 向 。从文 ［ １ ］ 的分 析 中我
Ｖｏ Ｉ ． ２ ６ Ｎｏ ． ２
２Ｏ１ ３
・
大 学教 学 ・
一
类带线性耗散项 的 ２ ×２ 双 曲型方程组解 的衰减估计
蒋 盼 盼
（ 南京 航 空 航 天 大 学 理 学 院 ， 江苏 南京 ２ １ １ １ ０ ６ ）
摘
要： 本文利用 Ｆ ｏ ｕ ｒ ｉ ｅ ｒ 分 析 方 法 研 究 了一 类 带 线 性 耗 散 项 的 ２ × ２双 曲 型 方 程 组 Ｃ ａ ｕ ｃ ｈ ｙ
Ｖｏ １ ． ２ ６ Ｎｏ ． ２
【计算流体力学】第3讲-差分方法1
a2u j
a3u j1+a4u j+2
扰动波传播方向
… j-2 j-1 j j+1 …
更多地使用上游信息
一般双曲守恒律方程
u f (u) 0 t x
f (u) f (u) f (u)
u f + f 0 t x x
df (u) 0 du
df (u) 0 du
例：
f 1 f u
u x j
时间积分，计算 出下一时刻的值
u lim u(x x) u(x) u j1 u j
x j x0
x
x
沿各自方向一维离散
➢多维方程的差分法： 维数分裂
u f1(u) f2 (u) 0 t x y
u
1. 构建差分格式
x j
已知均匀网格点上物理量的分布为uj ，
f1
x
f1
x
f2
y
f2
y
RAE2822翼型周 围的网格
问题： 原先需要计算2次导数，变换后需要计算4次，计算量增加 ✓利用坐标变换的性质，可以合并
14
坐标变换Jocabian系数的计算
已知 x x( ,)
y
y(
,)
需计算： x ,y ,x ,y
Step 1: 利用差分（或其他方法）计算出
网格间距变化要缓慢，否则会带 来较大误差
12
方法2） 在非等距网格上直接构造差分格式 （不易推广到高维）
原理： 直接进行Taylor展开，构造格式 格式系数是坐标（或网格间距）的函数
u x
j
a1u j2
a2u j1 a3u j
a4u j1 O(3 )
… j-2 j-1 j
双曲型方程求解方法及其应用
双曲型方程求解方法及其应用一、双曲型方程简介双曲型方程是一类二阶偏微分方程，其基本形式为:$$\dfrac{\partial^2 u}{\partial t^2}-\dfrac{\partial^2 u}{\partial x^2}=0$$双曲型方程的特点是存在两个独立的传播方向，解的形式通常是由两个波的叠加而成。
由于双曲型方程与空间和时间的关系有关，因此在物理、工程和科学领域中有着广泛的应用。
其中，双曲型方程的求解方法是求解偏微分方程的重要研究内容之一。
二、双曲型方程的求解方法对于双曲型方程，我们需要采取适当的数学工具来解决。
下面介绍几种常用的双曲型方程求解方法。
1. 分离变量法分离变量法是求解偏微分方程的常用方法之一，对于双曲型方程也可以采用分离变量法求解。
例如，我们可以假设$u(x,t)=X(x)T(t)$，将偏微分方程代入得到：$$\dfrac{T''}{T}=\dfrac{X''}{X}=-k^2$$这是两个常微分方程，可以通过求解得到$T(t)$和$X(x)$的通解，再合并得到$u(x,t)$的通解。
其中，使用的边界条件和初值条件对应具体问题的不同而有所不同。
2. 特征线法特征线法是一种求解双曲型偏微分方程的有效方法。
其基本思想是沿着方程组的特征线进行积分，将原方程转化为一维常微分方程。
例如，对于双曲型方程$\dfrac{\partial^2 u}{\partial t^2}-\dfrac{\partial^2 u}{\partial x^2}=0$，经过变换得到：$$\dfrac{\mathrm{d}u}{\mathrm{d}t}+\dfrac{\mathrm{d}u}{\mathrm{d}x}=0$$将$\dfrac{\mathrm{d}x}{\mathrm{d}t}=1$和$\dfrac{\mathrm{d}t}{\mathrm{d}u}=1$代入得到方程:$$\dfrac{\mathrm{d}x}{\mathrm{d}u}=\dfrac{1}{2},\dfrac{\mathrm{d}t}{\mathrm{ d}u}=-\dfrac{1}{2}$$由此可以得到$x=t+c_1,u=c_2$为特征线，设$u=f(x-t)$，则原方程变成$\dfrac{\mathrm{d}^2 f}{\mathrm{d} x^2}=0$，通解为$f(x-t)=k_1 x+k_2$，因此原方程的通解为$u(x,t)=k_1 x+k_2$。
线性严格双曲型方程组的初边值问题及其应用
第29卷第3期2008年6月华 北 水 利 水 电 学 院 学 报Journa l of Nort h China Institut e of W ate r Conservancy and Hydroe l ec tric Powe rVol 129No 13Jun .2008收稿日期3基金项目河南省教育厅自然科学研究资助项目();河南省基础研究资助项目作者简介刘法贵(65—),男,河南内乡人,教授,博士,主要从事偏微分方程方面的研究文章编号:1002-5634(2008)03-0099-03线性严格双曲型方程组的初边值问题及其应用刘法贵,赵 娟,郭吉强(华北水利水电学院,河南郑州450011)摘 要:考虑线性严格双曲型方程组初边值问题,利用特征线方法和局部解延拓方法,证明了其整体经典解的存在性和解对初始数据、边界数据的连续依赖性.关键词:线性双曲型方程组;初边值问题;整体经典解中图分类号:O175.27;O175.8 文献标识码:A 对于含2个自变量的一阶线性严格双曲型方程组U t +A (x,t)U x +B (x,t)U +C (x,t)=0(1)的Cauchy 问题整体经典解的存在惟一性和解的连续依赖性在许多教学参考书中都有论述[1-2],但其初边值问题则很少涉及[3-4].其中U =(u 1,u 2,…,u n )T,A (x,t),B (x,t)是具有光滑元的n 阶矩阵,C (x,t )=(c 1(x,t ),c 2(x,t ),…,c n (x,t ))T是光滑向量函数.利用特征线方法和局部解延拓方法,在区域Ω={(x,t)0≤x ≤l,t ≥0}上证明了方程组(1)整体经典解的存在惟一性和解对初始数据、边界数据的连续依赖性.1 主要结果由严格双曲型假设,不妨设方程组(1)的特征根和左特征向量分别为λ1(x,t)<λ2(x,t)<…<λm (x,t)<0<λm +1(x,t )<…<λn (x,t ) (2)L (x,t )=(l 1(x,t )l 2(x,t ),…,l n (x,t ))T式中l i (x,t)为行向量.由L ≠0,通过引进可逆变换V =LU ,方程组(1)可化为V t +Λ(x,t )V x =E (x,t )V +F (x,t )(3)Λ(x,t )=diag (λ1(x,t ),λ2(x,t ),…,λn (x,t ))E (x,t)=(e ij (x,t))n ×n ,F (x,t)=(f i (x,t))1×n由于方程组(1)与方程组(3)等价,因此,考虑方程组(3)具有初始条件(4)和边界条件(5)的初边值问题t =0∶v i =φi (x),0≤x ≤l,i =1,2,…,n(4) x =0∶v j =g j (t ),t ≥0,(j =m +1,m +2,…,n )x =l ∶v i =g i (t),t ≥0,(i =1,2,…,m )(5)式中:φk (x )∈C 1([0,l ]);g k (t )∈C 1(R +);k =1,2,…,n;g k (t )具有有界的C 1模.定理1 假设式(2)成立,初始条件和边界条件满足相容性条件,λi (x,t),e ij (x,t),f i (x,t)∈C 1且具有有界的C 1模,那么初边值问题(3)—(5)在区域Ω上存在惟一的整体经典解,且解连续依赖于初始数据和边界数据.2 定理的证明2.1 经典解整体存在的惟一性为证明经典解整体存在的惟一性,只须证明V (x,t )的C 1模在Ω上有界[5].设过点(t,x )的特征线为L i ∶x =x i (t), d xd t=λi (x,t)由假设式(2),过原点及(0,l )的特征线将区域Ω划分为4个子区域D 1={(x,t )x n (t,0)≤x ≤x 1(t ,l ),t ≥0}D 2={(x,t )0≤x ≤x 1(t,l ),t ≥0}-D 1D 3=D 1={(x,t )x n (t,0)≤x ≤x 1(t,0),t ≥0}:2008-0-20:2008A110011.:19.D 4={(x,t )x 1(t ,l )≤x ≤x n (t,0),t ≥0}在区域D 1上,利用有关Cauchy 问题的结论,易知初边值问题存在惟一的整体经典解,且其解连续依赖于初始条件.在区域D 2上,过点(t,x )向下引第1,2,…,m 条特征线交于区间[0,l]内一点(0,x i ),而第m +1,m +2,…,n 条特征线交于t 轴上点(t j ,0).这样,沿特征线积分方程组(3)得v i (x,t )=φi (x i )+∫0t∑nk =1e ik(x i (s ),s )v k (x i (s ),s )d s +∫tf i(x i (s),s)d s, (i =1,2,…,m )v j (x,t)=g j (t j )+∫t j t∑nk =1e jk(x i (s),s)v k (x i (s),s)d s +∫t jt f j(x j(s),s)d s,j =m +1,m +2,…,n由此,容易得到M (t )≤M 1+M 2∫tM (t )d t这里M (t)=max 1≤k ≤nsupxv k (x,t)M i (i =1,2,…)表示正常数(可以依赖于t ).因此,由B ell m an 不等式,即得M (t)≤M 3.在区域D 3上,过点(t,x )向下引第1,2,…,m 条特征线交于直线x =l 上点(t j ,0),第m +1,m +2,…,n 条特征线交于区间[0,l ]内一点(0,x i ).则沿特征线积分方程组(3),类似可得到式(6).在区域D 4上,过点(t ,x)向下引特征线分别交直线x =0和x =l 于(t i ,0)(i =m +1,m +2,…,n )和(t j ,l )(j =1,2,…,m ).则沿特征线积分方程组(3)得v i (x,t )=g i (t i )+∫t it∑nk =1e ik(x i (s ),s )v k (x i (s ),s )d s +∫t itf i(x i(s ),s )d sv j (x,t )=g j (t j )+∫t jt∑nk =1e jk (x j (s ),s )v k (x i (s ),s )d s +∫t jtf j(x j(s ),s )d s 由此,由假设条件,同样M (t)≤M 3成立.因此,在区域Ω上,V (x,t )的C 0模有界.为证明V x (x,t)的C 0模在区域Ω上有界,对方程组(3)两端关于x 求导,得(V x )t +Λ(x,t )(V x )x =(E -Λx )V x +E x V +F x (6)则类似估计V (x,)的模有界性的方法可得到V x ≤M ,因此,初边值问题(3)—(5)在区域Ω上存在惟一的整体经典解[5]2.2 解对初始数据和边界数据的连续依赖性设初边值问题(3)—(5)对应的经典解为V (x,t ),方程组(3)对应条件(4)和(5)的经典解记为V (x,t ).令Y (x,t )=V (x,t )-V (x,t )y i (x )=φi (x )-φi (x )p i (t )=g i (t )-g i (t )则Y t +Λ(x,t)Y x =E (x,t)Y (7)t =0∶Y i =y i (x ),0≤x ≤l,i =1,2,…,n(8)x =0∶Y j =p j (t ),t ≥0,(j =m +1,m +2,…,n )x =l ∶Y i =p i (t ),t ≥0,(i =1,2,…,m )(9)在区域D 1上,利用有关Cauchy 问题的结论,易知初边值问题的经典解连续依赖于初始条件.在区域D 2上,沿特征线积分方程组(7)得y i (x,t)=y i (x i )+∫0t∑nk =1e ik(x i (s),s)y k (x i (s),s)d s,i =1,2,…,my j (x,t)=p j (t j )+∫t jt∑nk =1e jk(x i (s),s)y k (x i (s),s)d s,j =m +1,m +2,…,n由此,容易得到‖y i (x,t)‖≤‖y i (x)‖+M 5∫t‖y ‖d s,i =1,2,…,m ‖y j (x,t)‖≤‖p j ‖+M 6∫t‖y ‖d s,j =m +1,m +2,…,n因此,在区域D 2上,经典解连续依赖于初始条件和边界条件.同理,在区域D 3和D 4上,初边值问题的经典解连续依赖于初始条件和边界条件.定理1证毕.注记1:方程组(3)的线性性质是必要的.对于非线性问题,其经典解可能会在有限时间内产生奇性.例如方程组 5u 5t +5u 5x -u +w -1=05v 5t +25u 5x -5w 5x-2u +w -v 2+2w v -w 2=05w 5t +25u 5x -5w 5x-2u +w =0(10)其特征值为 λ1=-1, λ2=0, λ3=1对应的左特征向量分别为=(,,),=(,,),3=(,,)则此方程组为半线性严格双曲型方程组考虑方程组()的初边值问题001 华 北 水 利 水 电 学 院 学 报 2008年6月t C 04.l 110-1l 201-1l 100.11 t =0∶u =u 0(x ),v =v 0(x ),w =w 0(x ),x ≥0x =0∶u =g (t ),t ≥0x =l ∶u =w(11)其中u 0(x ),v 0(x ),w 0(x )∈C 1([0,l ]),g (t )∈C 1(R +)具有有界的C 1模,且v 0(x )-w 0(x )=1(12)容易计算(v -w )t =(v -w )2(13)则得(v -w )(x,t )=(1-t )-1因此,当t →1-时,(v -w )(x,t )→∞.所以解在t =1处发生爆破.注记2:利用上述方法可以考虑方程组(3)具有如下边界条件的初边值问题x =0∶v j =f j (v 1,v 2,…,v m ,t ),(j =m +1,m +2,…,n )x =l ∶v i =f i (v m+1,v m +2,…,v n ,t ),(i =1,2,…,m )其中f j ,f i ∈C 1,且f j (0,0,…,0,t )=0,f i (0,0,…,0,t)=0能得到相仿的结论.3 应 用当管道中流量迅速改变(例如用阀门进行调节)时,水的加速或减速是伴有很大的压力变化,这种现象称为水锤(wate r ha mme r ).考虑到实用上的原因,探讨水流变化速率与所产生的水锤压力之间的关系是重要的.在研究水锤时,应该将水与管道看作是可以变形的,否则,管道中的水流突然停止流动意味着管道中水的全部质量的无限大的减速,这就会迫使阀门处产生无限大的压力.Jouko w sky -A llievi (出口处流量变化引起)水锤问题作为一维问题来研究,其方程组表达式为5v 5t +5<5x =-2kv 5<5t +c 25v5x=0(14)式中:v (x,t)为流速;<=p ρ+g h;p (x,t)为水的压力强度;ρ(x,t )为水的密度;h 为管道的高程;g 为重力加速度;c >0为常数.容易计算方程组(14)的特征值为λ1=-c, λ2=c考虑下面的初边值问题5v 5t +5<5x=-2kv 5<5t +c 25v5x=0t =0∶v =f (x ),φ=g (x ),x ∈[0,l ]x =0∶φ=cv +h 1(t ),t ≥0-x =l ∶φ=-cv +h 2(t ),t ≥0(15)式中:f (x ),g (x )∈C 1;h 1(t ),h 2(t )∈C 1,且具有有界的C 1模.定理2 假设f (x ),g (x ),h 1(t ),h 2(t )满足相容性条件,则初边值问题(15)在所考虑的区域上存在惟一的连续整体经典解,且经典解连续依赖于初始条件和边界条件.参 考 文 献[1]戴嘉尊.数学物理方程[M ].南京:东南大学出版社,2002.[2]陈祖墀.偏微分方程[M ].北京:中国科学技术出版社,1993.[3]谷超豪,李大潜.数学物理方程[M ].北京:高等教育出版社[M ],2002.[4]刘法贵.数学物理方程[M ].郑州:黄河水利出版社,2007.[5]L i Ta tian,YuWenci .BoundaryVa lue P roblem forQ uasili n 2ea r Hyperbolic Syste m s [M ].US A:Duke Unive rsity Pre ss,1985.I n it i a l Bounda r y V a lue Pr oblem a nd Its App lic a t i on s for Str i c t L i nea r Hyper bol ic Equa t i on sL IU Fa 2gui,ZHA O Juan,G UO J i 2qiang(North China I nstitute of W ater Conserv ancy and Hydr oelectric Po wer,Zhengzhou 450011,China )Ab stra ct:Considering the initia l boundary value p roblem for strict linea r hy pe rbolic equa tions,and usi ng the cha racteristicsm ethod and l oca l extansion m ethod,the gl obal existence of c l a ssi ca l soluti ons and conti nuous dependence of cla ssi ca l sol u ti ons are proved .K y y q ;y ;101第29卷第3期刘法贵等: 线性严格双曲型方程组的初边值问题及其应用 e w or d s:linea r h p erbo lic e ua ti o n s in itia l boun da r valu e p r o b le m gl o bal classical s o lu tion。
双曲型偏微分方程组的数学模型及算法研究
双曲型偏微分方程组的数学模型及算法研究双曲型偏微分方程组是数学中的一个重要分支，其广泛应用于物理学、工程学、生物学、经济学以及计算机科学等多个领域。
本文将介绍双曲型偏微分方程组（以下简称双曲型PDE）的基本概念，数学模型及其算法研究。
一、双曲型偏微分方程组的基本概念双曲型偏微分方程可以简单地表示为：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=f(x,t)$$其中，$u(x,t)$是待求函数，$f(x,t)$是已知函数，$a$是常数。
对于双曲型偏微分方程中的函数$u(x,t)$，其趋势和形状通常会随着时间或空间的变化而发生变化。
这种性质决定了双曲型PDE的求解方法与其它类型偏微分方程组不同。
二、双曲型偏微分方程的数学模型在实际问题中，双曲型偏微分方程可以用来描述声波、水波、热传导等现象。
以声波方程为例，我们可以得到：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=0$$此时$f(x,t)=0$。
该方程表示了声波在空气中的传播，其中$a$是声速，$u(x,t)$表示声波幅度。
可以看到，随着时间的推移，声波的幅度会发生变化，而空气中声波的传播速度$a$是固定不变的。
这种性质决定了声波传播方程是一个双曲型偏微分方程。
同样地，在热传导问题中，我们也可以得到一个双曲型偏微分方程模型：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=k\frac{\partial u}{\partial t}$$其中，$k$是热扩散系数。
这个方程描述了热传导的过程。
可以看到，随着时间的变化，温度分布图的形状和趋势也会随之改变。
双曲型偏微分方程模型的重要性在于其可以精确描述相应现象的物理过程，从而为实际应用提供基础和便利。
“非齐次拟线性双曲型方程组整体经典解”一文的注记
考 虑 如 下 常 微 分 方 程 的 Ｃ ｕ ｈ 题 ａ ｃｙ问
（） ｙ＿ ｌ ｆ ］ Ｙ－ ２ ， 刈 ｆ ［ ｙ（ ） ［ Ｙ ］ ｆ
（） ７
【 （ ）＝Ｙ Ｙ０ 。
引 理 ， 明 了 Ｃｕｈ 证 ａ ｃｙ问 题 （ ） 体 经 典 解 的存 在 １ 整 性 ， 考 虑 了解 的破 裂 问题 ， 富 了 “ 齐 次 拟线 性 并 丰 非 双 曲型方 程组整 体经典 解 ” 文 的主要结论 ． 一
维普资讯
第２ ８卷 第 ６期
２０ ０ ７年 ｌ ２月
华
北
水
利
水 电
学
院
学
报
Ｖ０ ． ． １２８ Ｎｏ ６
Ｄｅ ． ２ ０７ ｃ ０
Ｊ ｕ ｎ ｌｏ ａｈ Ｃｈ ｎ ｎ ｔ ｕｅ ｏ ａ ｅ ｎ ｅｖ ｎ ｙ ａ ｄ Ｈｙ ｒ ｅｅ ｔ ｃ Ｐｏ ｒ ｏ ｒ ａ ｆＮｏ ｉ ａ Ｉ ｓｉ ｔ ｆＷ ｔ ｒＣｏ ｓ ｒ ａ ｃ ｎ ｄ ｏ ｌｃｒ ｗｅ ｔ ｉ
要 ： 用 常 微 分 方 程 Ｃ ｕｈ 利 ａｃｙ问题 的 比较 引理 ， 明 了 一 类 拟 线 性 双 曲 型 方 程 Ｃ ｕｈ 证 ａｃｙ问 题 整 体 经 典 解 的
存 在 性 ， 富 了“ 齐 次 拟 线 性 双 曲 型 方 程 组整 体 经 典 解 ” 丰 非 一文 的 主要 结 论 ． 关 键 词： 比较 引 理 ； 拟线 性 双 曲型 方 程 Ｃ ｕｈ ａｃ ｙ问题 ； 体 经 典 解 整
Ｙ ￡一 ∞ （） 定 理 １的证明 （ ２ １） 在方 程 （ ） 边对 求 导 ， １两 得
则 Ｃｕｈ ａ ｃｙ问题 （ ） ｔ ０上存 在惟 一整体 经典解 ． １在 ＞ －
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1
1) 一阶常系数偏微方程组
U U A 0 t x U (u1 , u 2 ,...... u m )T
如果矩阵A 可以被对角化： A S 1 ΛS
U U S1 ΛS 0 t x S U U ΛS 0 t x
diag(1 , 2 ,...... m )
第2 讲 双曲型方程组及其间断解
§2.4 双曲型方程及其数学性质
1. 双曲方程边界条件的提法
考虑方程组： 令：
f(U) A U U f(U) 0 t x
U (u1 , u2 ,...... u m )T
U U A 0 t x
如果矩阵A 能通过相似变换对角化
双曲型
Copyright by Li Xinliang
考虑一维无粘流动（Euler方程），初始时 刻（t=0）流动状态如下：
可推广到一般的双曲型方程组
3
Copyright by Li Xinliang
2） 一维Euler方程
U F(U) 0 t x
U ( , u, E )T
A F(U) U
u 2 F(U) u p ( E p)u
A S 1 ΛS
令：
Ru
d
c
同理，沿特征线2：
dx / dt u c
du c d 0 d d u 1 c S d 2 2 dR / d 沿特征线1： R不变来自对于等熵完全气体
2c R u 1 2c S u 1
寻找积分因子：
R R( , u )
dR R u R d u
R R c 1; u
沿特征线1： dx / dt u c 有：
d (u c) d t x
第1个方程转化为 则有
du c d 0 d d
保持不变
Riemann 不变量 知识点，牢记！
7 Copyright by Li Xinliang
一维均熵流动沿特征线Riemann不变量保持不变
例2.1： 有限振幅波的传播问题
目的： 学会如何运用Riemann不变量解题
A sin x 0 x 2 u ( x,0) 0 others ( x,0) 1; p( x,0) 1
给定3个边界条件 给定2个边界条件 无需给定边界条件 给定1个边界条件
4
知识点
Copyright by Li Xinliang
2. 双曲型方程组的特征方程
U U A(U) 0 t x U (u1 , u 2 ,...... u m )T
变系数方程组的情况
令：
A S 1 ΛS
U U S1 ΛS 0 t x S U U ΛS 0 t x ωk ( S( U U Λ )0 t x
令
1 S ... m
（行向量）
x x( sk ) 满足： dx / dsk k t t ( sk ) t sk
令： 即：
V SU
有
v j x 0
V V Λ 0 t x
v j t
j
m个方程完全解耦， 可独立求解
x jt 0 有m 条特征线： m个特征相容关系式： v j const.
如果矩阵A能够（相似变换）对角化，则原方程是双曲型的
2 Copyright by Li Xinliang
diag(1 , 2 , 3 )
1 u, 2 u c, 3 u c
对于左边界：
条件
u 0 and u 0 and u 0 and u 0 and u c u c u c u c
描述
边界条件设定
超音速入口 亚音速入口 超音速出口 亚音速出口
则有：
R 0
S 0
R R( )
S S ( )
Riemann不变量
沿特征线： Riemann不变量保持不变
6
Copyright by Li Xinliang
常见情况讨论: 一维等（均）熵运动
U U B 0 t x
U u
u B c 2 /
u
例如，膨胀波
矩阵B的特征值
1, 2 u c
c S c
u u c[ t (u c) x ] [ t (u c) x ] 0 u u c[ (u c) ] [ (u c) ] 0 t x t x
对于两自变量情况,可化为：
11
21
u1 u 12 2 0
u1 u 22 2 0
2个常微方程
如果存在积分因子，使得
* du * du dR 111 1 1 12 2
* du * du dS 221 1 2 22 2
双曲方程边界条件提法
U U A 0 t x
v j t
j
v j x
0
j=1
变换成为了彼此独立的n个单波方程
j=2
方法： 独立给定j个方程的边界条件
如果 j>0, 则在左端给定vj的边界条件 如果 j<0, 则在右端给定vj的边界条件
A B
特点： 左、右边界总共给定n个边界条件，各自的个数视特征 值的符号确定
U U k )0 t x
U j t k U j x )0
在x-t空间引入曲线：
d dx dt dsk x dsk t dsk
j
kj (
j
kj
dU j dsk
0
（变系数情况）虽然不能解耦，但还能转换成常微方程
5 Copyright by Li Xinliang