[物理学与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 一维反应流体力学方程组的数学结构
一维理想反应流体力学方程组是一阶拟线性双曲组.
随机推荐
- [原创]GDB调试指南-断点设置
前言 上篇<GDB调试指南-启动调试>我们讲到了GDB启动调试的多种方式,分别应用于多种场景.今天我们来介绍一下断点设置的多种方式. 为何要设置断点 在介绍之前,我们首先需要了解,为什么需 ...
- How-to: Do Real-Time Log Analytics with Apache Kafka, Cloudera Search, and Hue
Cloudera recently announced formal support for Apache Kafka. This simple use case illustrates how to ...
- SpringMVC项目读取不到外部CSS文件的解决办法及总结
昨天,在写一个新需求的时候,因为没有前端同事的对接,无奈只有自己写css,js放到ssm结构的后台项目中,因为之前开发都是把前端文件放在一个专门的服务器上,一直没有关注服务端项目读取静态资源的问题,运 ...
- Linux内存管理 (4)分配物理页面
专题:Linux内存管理专题 关键词:分配掩码.伙伴系统.水位(watermark).空闲伙伴块合并. 我们知道Linux内存管理是以页为单位进行的,对内存的管理是通过伙伴系统进行. 从Linux内存 ...
- .net core 的图片处理及二维码的生成及解析
写代码这事,掐指算来已经十有余年. 从html到css到javascript到vbscript到c#,从兴趣到职业,生活总是失落与惊喜并存. 绝大部分时候,出发并不是因为知道该到哪里去,只是知道不能再 ...
- Auto Layout: Programmatic Constraints - BNR
继续Auto Layout - BNR篇. 打开BNRDetailViewController.m文件,重载viewDidLoad方法来创建UIImageView对象.当你想要给通过加载NIB文件创建 ...
- python使用redis
版本: python 3.5 redis 3.0.1(redis的安装 pip install redis) 1.连接 import redis r = redis.Redis(host='192.1 ...
- 一篇博客带你入门Flask
一. Python 现阶段三大主流Web框架 Django Tornado Flask 对比 1.Django 主要特点是大而全,集成了很多组件,例如: Models Admin Form 等等, 不 ...
- React Native之配置URL Scheme(iOS Android)
React Native之配置URL Scheme(iOS Android) 一,需求分析 1.1,需要在网站中打开/唤起app,或其他app中打开app,则需要设置URL Scheme.比如微信的是 ...
- vue-cli的跨域配置(自己总结)