[物理学与PDEs]第2章第4节 激波 4.1 间断连接条件
1. 守恒律方程 $$\bex \cfrac{\p f}{\p t}+\cfrac{\p q}{\p x}=0 \eex$$ 在间断线上应满足 ``间断连接条件'': $$\bex [f]\cfrac{\rd x}{\rd t}=[q]. \eex$$
2. 对一维理想流体力学方程组 $$\beex \bea \cfrac{\p\rho}{\p t}+\cfrac{\p}{\p x}(\rho u)&=0,\\ \cfrac{\p}{\p t}(\rho u) +\cfrac{\p}{\p x}(\rho u^2+p)&=\rho F,\\ \cfrac{\p}{\p t}\sex{\rho e+\cfrac{1}{2}\rho u^2} +\cfrac{\p}{\p x}\sez{\sex{ \rho e+\cfrac{1}{2}\rho u^2+p }u}&=\rho Fu, \eea \eeex$$ 其间断连接条件为 $$\beex \bea [\rho]\cfrac{\rd x}{\rd t}&=[\rho u],\\ [\rho u]\cfrac{\rd x}{\rd t}&=[\rho u^2+p],\\ \sez{\rho e+\cfrac{1}{2}\rho u^2}\cfrac{\rd x}{\rd t} &=\sez{\sex{\rho e+\cfrac{1}{2}\rho u^2}u}. \eea \eeex$$ 此称为 Rankine-Hugoniot (R.H.) 条件.
3. 设 $U=\cfrac{\rd x}{\rd t}$ 为间断的传播速度, 记 $$\bex v_\pm=u_\pm-U, \eex$$ 则 R.H. 条件可化为 $$\beex \bea \rho_-v_-&=\rho_+v_+,\\ \rho_-v_-^2+p_-&= \rho_+v_+^2+p_+,\\ \sex{\rho_-e_-+\cfrac{1}{2}\rho_-v_-^2+p_-}v_-&=\sex{\rho_+e_++\cfrac{1}{2}\rho_+v_+^2+p_+}v_+. \eea \eeex$$
4. 记 $m=\rho_-v_-=\rho_+v_+$, 则
(1) 若 $m=0$, 则 $x=x(t)$ 为接触间断 (contact discontinuity), 此时, $$\bex v_-=v_+=0\ra u_+=u_-=U, \eex$$ 该间断线随流体一以同一速度运动, 无流体越过间断线.
(2) 若 $m\neq 0$, 则 $x=x(t)$ 为激波, 在越过激波时, 由 R.H. 条件可导出各热力学量应满足的方程. 比如 $$\bex H(\tau,p;\tau_0,p_0)\equiv e(\tau,p)-e_0(\tau,p) +\cfrac{1}{2}(p_0+p)(\tau-\tau_0)=0. \eex$$ 此称为 Hugoniot 方程或热力学激波条件 (只依赖于热力学量 $\tau$, $p$). 另外, $v_-$, $v_+$ 同号. 若同为负号, 则 $u_-,u_+<U$, 而流体自右向左越过激波, 而激波相对于流体来说向右运动, 称为右传播激波; 若同为负号, 则称为左传播激波 (书 P 134).
[物理学与PDEs]第2章第4节 激波 4.1 间断连接条件的更多相关文章
- [物理学与PDEs]第2章第4节 激波 4.2 熵条件
1. R.H. 条件仅仅给出了越过激波时的能量守恒定律, 即热力学第一定律; 但客观的流体运动过程还需满足热力学第二定律, 即越过激波是个熵增过程: $$\bex S_1>S_0\quad(0 ...
- [物理学与PDEs]第5章第1节 引言
1. 弹性力学是研究弹性体在荷载的作用下, 其内力 (应力) 和变形所满足的规律的学科. 2. 荷载主要有两种, 一是作用在弹性体上的机械力 (本章讨论); 二是由温度等各种能导致弹性体变形的物理 ...
- [物理学与PDEs]第4章第1节 引言
1. 本章讨论可燃流体在流动过程中同时伴随着燃烧现象的情况. 2. 燃烧有两种, 一种是爆燃 (deflagration): 火焰低速向前传播, 此时流体微元通常是未燃气体.已燃气体的混合物; 一 ...
- [物理学与PDEs]第5章第6节 弹性静力学方程组的定解问题
5. 6 弹性静力学方程组的定解问题 5. 6. 1 线性弹性静力学方程组 1. 线性弹性静力学方程组 $$\bee\label{5_6_1_le} -\sum_{j,k,l}a_{ijkl}\cf ...
- [物理学与PDEs]第5章第5节 弹性动力学方程组及其数学结构
5.5.1 线性弹性动力学方程组 1. 线性弹性动力学方程组 $$\beex \bea 0&=\rho_0\cfrac{\p{\bf v}}{\p t}-\Div_x{\bf P}-\r ...
- [物理学与PDEs]第5章第4节 本构方程 - 应力与变形之间的关系
5. 4 本构方程 - 应力与变形之间的关系 5.4.1. 本构关系的一般形式 1. 若 Cauchy 应力张量 ${\bf T}$ 满足 $$\bex {\bf T}({\bf y})=\hat{\ ...
- [物理学与PDEs]第5章第3节 守恒定律, 应力张量
5. 3 守恒定律, 应力张量 5. 3. 1 质量守恒定律 $$\bex \cfrac{\p \rho}{\p t}+\Div_y(\rho{\bf v})=0. \eex$$ 5. 3. 2 应 ...
- [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.3 位移梯度张量与无穷小应变张量
1. 位移向量 $$\bex {\bf u}={\bf y}-{\bf x}. \eex$$ 2. 位移梯度张量 $$\bex \n_x{\bf u}={\bf F}-{\bf I}. \eex$ ...
- [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.2 Cauchy - Green 应变张量
1. 引理 (极分解): 设 $|{\bf F}|\neq 0$, 则存在正交阵 ${\bf R}$ 及对称正定阵 ${\bf U},{\bf V}$ 使得 $$\bex {\bf F}={\bf ...
随机推荐
- 使用Visual Studio Code进行ABAP开发
长期以来,我们都使用SAP GUI进行ABAP编码工作,事务代码SE38甚至成了ABAP的代名词. SAP GUI的代码编辑能力和一些专业的IDE比较起来难免相形见绌,为了给开发者们更好的体验,SAP ...
- 合并两个有序链表的golang实现
将两个有序链表合并为一个新的有序链表并返回.新链表是通过拼接给定的两个链表的所有节点组成的. 输入:->->, ->-> 输出:->->->->-> ...
- 打印进度条>>>>
i+1: 当前的数量 300: 总数量 import sys print("下载中...") def process(curr, count): cursor_count = c ...
- tian
上次后来没继续在微信上聊,是因为快过年了,想趁那段时间结合年假做点东西.接下来阳历三四月份就受美国制裁.结果接下来制裁.fang. 16年的那次主要是生气,在一块儿 另外经济上也有问题. 我也想过不再 ...
- 研究好vif 和vshow
另外从源头上处理的???,怎么自己排查出错误??必须 ??https://www.jb51.net/article/124116.htm
- web框架开发-快速认识Django中间件
中间件 中间件的概念 中间件顾名思义,是介于request与response处理之间的一道处理过程,相对比较轻量级,并且在全局上改变django的输入与输出. 因为改变的是全局,所以需要谨慎实用,用不 ...
- 利用XShell上传、下载文件(使用sz与rz命令)
XSHELL工具上传文件到Linux以及下载文件到本地(Windows) Xshell很好用,然后有时候想在windows和linux上传或下载某个文件,其实有个很简单的方法就是rz,sz.首先你 ...
- linux下 启动node 和关闭node
1.用forever 进行管理 npm install -g forever forever start app.js //启动 forever stop app.js //关闭 2.用自带的服务n ...
- Linux内存管理 (23)一个内存Oops解析
专题:Linux内存管理专题 关键词:DataAbort.fsr.pte.backtrace.stack. 在内存相关实际应用中,内存异常访问是一种常见的问题. 本文结合异常T32栈回溯.Oops ...
- java获取真实的IP地址工具类
在实际项目中,有调用微信支付完成支付功能,在微信支付的请求参数中需要传递一个本机的ip地址,java代码运行环境目前为windows10以及centos7. 以下为获取ip地址工具类: package ...