[物理学与PDEs]第2章第1节 理想流体力学方程组 1.3 理想流体力学方程组的数学结构
1. 局部音速 $c$: $c^2=\cfrac{\p p}{\p \rho}>0$.
2. 将理想流体力学方程组 $$\beex \bea \rho\cfrac{\p {\bf u}}{\p t} +(\rho {\bf u}\cdot\n){\bf u}+\n p&=\rho{\bf F},\\ \cfrac{1}{\rho c^2}\cfrac{\p p}{\p t} +\n\cdot{\bf u}+\cfrac{1}{\rho c^2}({\bf u}\cdot\n)p&=0,\\ \cfrac{\p S}{\p t}+({\bf u}\cdot\n)S&=0 \eea \eeex$$ 写成 $$\bee\label{2_1_2_dc} A_0\cfrac{\p U}{\p t} +\sum_{i=1}^3 A_i\cfrac{\p U}{\p x_i}=C, \eee$$ 其中 $U=(u_1,u_2,u_3,p,S)^T$, 则有 $$\beex \bea A_0=\sex{\ba{ccccc} \rho &&&&\\ &\rho&&&\\ &&\rho&&\\ &&&\cfrac{1}{\rho c^2}&\\ &&&&1 \ea},&\quad A_1=\sex{\ba{ccccc} \rho u_1&&&1&\\ &\rho u_1&&&\\ &&\rho u_1&&\\ 1&&&\cfrac{u_1}{\rho c^2}&\\ &&&&u_1 \ea},\\ A_2=\sex{\ba{ccccc} \rho u_2&&&&\\ &\rho u_2&&1&\\ &&\rho u_2&&\\ &1&&\cfrac{u_2}{\rho c^2}&\\ &&&&u_2 \ea},&\quad A_3=\sex{\ba{ccccc} \rho u_3&&&&\\ 0&\rho u_3&&&\\ &&\rho u_3&1&\\ &&1&\cfrac{u_3}{\rho c^2}&\\ &&&&u_3 \ea},\\ C=(\rho F_1,\rho F_2,\rho F_3,0,0)^T.& \eea \eeex$$
3. 当 $\rho>0$ 时, \eqref{2_1_2_dc} 为一阶拟线性对称双曲型偏微分方程组. 而可考虑 Cauchy 问题、初-边值问题.
4. 理想流体力学方程组可化为一阶拟线性对称双曲组 $$\bex \cfrac{\p L^0_{v_i}}{\p t} +\sum_{k=1}^3 \cfrac{\p }{\p x_k}L^k_{v_i}=0,\quad i=0,1,\cdot,4. \eex$$ 其中 $$\bex L_0=-\cfrac{p}{T},\quad L^k=-\cfrac{p}{T}u_k\ (k=1,2,3). \eex$$ 这里,
(1) $L_{v_iv_j}$ 为对称正定阵.
(2) $v_i\ (i=0,1,\cdots,4)$ 及 $L$ 为 $\rho,\rho u_1,\rho u_2,\rho u_3, \rho e+\cfrac{1}{2}\rho u^2$ 及 $\rho S$ 的 Legendre 变换.
5. 一般的守恒律方程组可化为一阶对称双曲组的一个充要条件
设有守恒律方程组 $$\bee\label{2_1_shl} \cfrac{\p U}{\p t}+\sum_{k=1}^3 \cfrac{\p}{\p x_k} B^k(U)=0, \eee$$ 其中 $$\bex U=(u_1,\cdots,u_n)^T,\quad B^k=(b^k_1,\cdots,b^k_n)^T. \eex$$ 则 \eqref{2_1_shl} 可通过未知函数变换 $$\bex U=U(V),\quad (u_i=u_i(v_1,\cdots,v_n),\ i=1,\cdots,n) \eex$$ 化为守恒律形式的一阶对称双曲组的充要条件为: 存在严格凸的标量 $W(U)$ 与向量函数 $H=(h_1(U),h_2(U),h_3(U))^T$, 使成立如下附加守恒律 $$\bex \cfrac{\p }{\p t}W(U)+\sum_{k=1}^3 \cfrac{\p}{\p x_k}h_k(U)=0. \eex$$ 这里, $W(U)$ 称为 \eqref{2_1_shl} 的熵函数, $H(U)$ 称为熵流函数.
[物理学与PDEs]第2章第1节 理想流体力学方程组 1.3 理想流体力学方程组的数学结构的更多相关文章
- [物理学与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 ...
- [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.1 变形梯度张量
$$\bex \rd{\bf y}={\bf F}\rd {\bf x}, \eex$$ 其中 ${\bf F}=\n_x{\bf y}=\sex{\cfrac{\p y_i}{\p x_j}}$ 为 ...
- [物理学与PDEs]第4章第3节 一维反应流体力学方程组 3.3 一维反应流体力学方程组的数学结构
一维理想反应流体力学方程组是一阶拟线性双曲组.
随机推荐
- C++中 #if 和 #ifdef 区别
以#开头的都是预编译指令,就是在正式编译之前,编译器做一些预处理的工作#if 条件语句程序段1 //如果条件语句成立,那么就编译程序段1#endif程序段2//如果条件不语句成立,那么就编译程序段2# ...
- Springboot配置文件解析器
@EnableScheduling @MapperScan(value = "com.****.dao") @EnableTransactionManagement @Enable ...
- canvas如何自适应屏幕大小
可以用JS监控屏幕大小,然后调整Canvas的大小.在代码中加入JS $(window).resize(resizeCanvas); function resizeCanvas() { ...
- 好程序员分享居中一个float元素
好程序员分享居中一个float元素,我们布局的时候,用margin来设置float元素的外边距来达到效果.对于,在文档流中的元素,我们很容易让它水平居中,只要给元素设置一个固定的宽度,用margin: ...
- 012_call和apply区别
一. function fn(a,b) { console.log(this); } fn.call(null,1,2); //call为参数方式 fn.apply(null,[1,2]); //ap ...
- ESP8266天线问题
http://www.icxbk.com/ask/detail/28832.html 目前市面上常见的外接天线包括 1.FPC天线,就是一小块柔性PCB,上面走一个铜线,下方不覆铜,然后一般带一个贴纸 ...
- i春秋 百度杯”CTF比赛 十月场 login
出现敏感的信息,然后进行登录 登录成功发现奇怪的show 然后把show放到发包里面试一下 出现了源码,审计代码开始 出flag的条件要user 等于春秋 然后进行login来源于反序列化后的logi ...
- 检测web界面不能访问后重启
检测并重启脚本:checkAndRestart.sh #!/bin/bash nowpath=$(cd ")";pwd) source $nowpath/omcparam.prop ...
- private,protected,public和default的区别
private,protected,public和default的区别 除了default以外,其他都是Java语言的关键字.default代表的是对类成员没有进行修饰的情况.它本身也代表了一种访问控 ...
- UIGestureRecognizer - BNR
继续上篇UITouch - BNR.该篇将实现线条选择.移动和删除操作. UIGestureRecognizer有一系列子类,每一个子类都用于识别特定的手势.当识别出一个手势时,手势识别器会拦截视图的 ...