试将一维理想磁流体力学方程组 (5. 10)-(5. 16) 化为一阶拟线性对称双曲组的形式.

解答: 由 (5. 12),(5. 16) 知 $$\beex \bea 0&=\cfrac{\p p}{\p \rho}\sex{\cfrac{\p \rho}{\p t}+u_1\cfrac{\p \rho}{\p x}+\rho \cfrac{\p u_1}{\p x}} +\cfrac{\p\rho}{\p S}\sex{\cfrac{\p S}{\p t}+u_1\cfrac{\p S}{\p x}}\\ &=\cfrac{\p p}{\p t}+u_1\cfrac{\p p}{\p x}+\tilde c^2\rho \cfrac{\p u_1}{\p x}, \eea \eeex$$ 而 $$\bex \cfrac{1}{\tilde c^2\rho }\cfrac{\p p}{\p t}+\cfrac{u_1}{\tilde c^2\rho}\cfrac{\p p}{\p x}+\cfrac{\p u_1}{\p x}=0.  \eex$$ 于是 (5. 10)-(5. 16) 为 $$\beex \bea \mu_0\cfrac{\p H_2}{\p t}+\mu_0u_1\cfrac{\p H_2}{\p x} +\mu_0H_2\cfrac{\p u_1}{\p x} -\mu_0H_1\cfrac{\p u_2}{\p x}&=0,\\ \mu_0\cfrac{\p H_3}{\p t} +\mu_0u_1\cfrac{\p H_3}{\p x} +\mu_0H_3\cfrac{\p u_1}{\p x} -\mu_0H_1\cfrac{\p u_3}{\p x}&=0,\\ \cfrac{1}{\tilde c^2\rho }\cfrac{\p p}{\p t}+\cfrac{u_1}{\tilde c^2\rho}\cfrac{\p p}{\p x}+\cfrac{\p u_1}{\p x}&=0,\\ \cfrac{1}{\rho}\cfrac{\p u_1}{\p t} +\cfrac{u_1}{\rho}\cfrac{\p u_1}{\p t} +\cfrac{\p p}{\p x} +\mu_0\sex{H_2\cfrac{\p H_2}{\p x}+H_3\cfrac{\p H_3}{\p x}}&=F_1,\\ \rho \cfrac{\p u_2}{\p t} +\rho u_1\cfrac{\p u_2}{\p x}-\mu_0H_1\cfrac{\p H_2}{\p x}&=F_2,\\ \rho \cfrac{\p u_3}{\p t}+\rho u_1\cfrac{\p u_3}{\p x} -\mu_0H_1\cfrac{\p H_3}{\p x}&=F_3,\\ \cfrac{\p S}{\p t}+u_1\cfrac{\p S}{\p x}&=0; \eea \eeex$$ 可写成 $$\bex A(U)\cfrac{\p U}{\p t}+B(U)\cfrac{\p U}{\p x}=C, \eex$$ 其中 $$\beex \bea U&=(H_2,H_3,p,u_1,u_2,u_3,S)^T,\\ A(U)&=\diag(\mu_0,\mu_0,\cfrac{1}{\tilde c^2\rho},\cfrac{1}{\rho},\rho,\rho,1),\\ B(U)&=\sex{\ba{ccccccc} \mu_0u_1&0&0&\mu_0H_2&-\mu_0H_1&0&0\\ 0&\mu_0u_1&0&\mu_0H_3&0&-\mu_0H_1&0\\ 0&0&\cfrac{u_1}{\tilde c^2 \rho}&1&0&0&0\\ \mu_0H_2&\mu_0H_3&1&\cfrac{u_1}{\rho}&0&0&0\\ -\mu_0H_1&0&0&0&\rho u_1&0&0\\ 0&-\mu_0H_1&0&0&0&\rho u_1&0\\ 0&0&0&0&0&0&u_1 \ea},\\ C&=(0,0,0,F_1,F_2,F_3,0)^T. \eea \eeex$$

[物理学与PDEs]第3章习题5 一维理想磁流体力学方程组的数学结构的更多相关文章

  1. [物理学与PDEs]第4章习题3 一维理想反应流体力学方程组的数学结构

    证明: Euler 坐标系下的一维反应流体力学方程组 (3. 10)-(3. 13) 也是一个一阶拟线性双曲型方程组. 证明: 由 (3. 10), (3. 12), (3. 13) 知 $$\bex ...

  2. [物理学与PDEs]第4章习题4 一维理想反应流体力学方程组的守恒律形式及其 R.H. 条件

    写出在忽略粘性与热传导性, 即设 $\mu=\mu'=\kappa=0$ 的情况, 在 Euler 坐标系下具守恒律形式的一维反应流动力学方程组. 由此求出在解的强间断线上应满足的 R.H. 条件 ( ...

  3. [物理学与PDEs]第3章习题6 Lagrange 坐标下的一维理想磁流体力学方程组的数学结构

    试讨论 Lagrange 形式下的一维理想磁流体力学方程组 (5. 33)-(5. 39) 的类型. 解答: 由 (5. 33), (5. 39) 知 $$\bex 0=\cfrac{\p p}{\p ...

  4. [物理学与PDEs]第2章习题10 一维理想流体力学方程组的 Lagrange 形式

    试证明: 一维理想流体力学方程组的 Lagrange 形式 (5. 22)-(5. 24) 也可写成如下形式 $$\beex \bea \cfrac{\p \tau}{\p t}-\cfrac{\p ...

  5. [物理学与PDEs]第2章习题8 一维定常粘性不可压缩流体的求解

    考察固定在 $y=0$ 与 $y=1$ 处两个平板之间的定常粘性不可压缩流体沿 $x$ 方向的流动. 设 $p=p(x)$, 且已知 $p(0) =p_1$, $p(L)=p_2$, $p_1> ...

  6. [物理学与PDEs]第2章习题7 一维不可压理想流体的求解

    设有以 $x$ 轴为轴向的等轴截面管道, 其中充满着沿 $x$ 方向流动的不可压缩的理想流体, 在每一横截面上流体的状态相同, 且 $p=p(x)$. 若已知 $p(0) =p_1$, $p(L)=p ...

  7. [物理学与PDEs]第2章习题6 有旋的 Navier-Stokes 方程组

    试证明: 由 Navier-Stokes 方程组描述的流体运动一般总是有旋的, 即若 $\rot{\bf u}={\bf 0}$, 则 Navier-Stokes 方程组 (3. 4)-(3. 5) ...

  8. [物理学与PDEs]第3章习题参考解答

    [物理学与PDEs]第3章习题1 只有一个非零分量的磁场 [物理学与PDEs]第3章习题2 仅受重力作用的定常不可压流理想流体沿沿流线的一个守恒量 [物理学与PDEs]第3章习题3电磁场的矢势在 Lo ...

  9. [物理学与PDEs]第4章习题参考解答

    [物理学与PDEs]第4章习题1 反应力学方程组形式的化约 - 动量方程与未燃流体质量平衡方程 [物理学与PDEs]第4章习题2 反应力学方程组形式的化约 - 能量守恒方程 [物理学与PDEs]第4章 ...

随机推荐

  1. idea软件破解汉化

    →http://idea.lanyus.com/上可以找到最新的破解补丁,下载并放到软件的bin目录下  →更改bin目录下的两个文件:Idea.exe.vmoptions和Idea64.exe.vm ...

  2. c++中 . 和 -> 的区别是什么?

    主要用于访问类的成员,->主要用于类类型的指针访问类的成员,而.运算符,主要用于类类型的对象访问类的成员. 例如: class A { public :int a } A ma; A *p=&a ...

  3. 超哥笔记--shell 基本命令(4)

    一 linux 命令行的组成结构 自定义命令行结构 PS1变量来控制 \u \W 最后一位工作目录 \w 绝对路径工作目录 \t 显示24h制的时间 \h PS1="[\u@\h \w \t ...

  4. F. Multicolored Markers(数学思维)

    思维:思维就是将大的矩形放在小矩形里面,让大矩形的宽和长尽量靠近. 很容易得到 (a+b)% i = 0 的话, 保证了大矩形的形成,同时里面表示了两种情况:1, a % i =0, b % i=0; ...

  5. 25 python 初学(socket,socketserver)

    参考blog :www.cnblogs.com/yuanchenqi/articles/5692716.html 1. sk = socket.socket() 里面有两个重要的参数,family 和 ...

  6. c++面经积累<2>

    4.类成员初始化方式:列表初始化和赋值初始化 赋值初始化通过在函数体内进行赋值,列表初始化,在构造函数后面加上冒号,使用初始化列表进行初始化.在函数体内进行初始化,是在所有的数据成员被分配内存空间后进 ...

  7. EF Core 多对多配置

    1.配置2个数据表 T_Authors ,T_Books 2.新建控制台项目,安装EF驱动 PM> Install-Package Pomelo.EntityFrameworkCore.Mysq ...

  8. C#调用Delphi DLL获取字符串(C# IntPtr 与 string互转 )

    前因后果 调用一门锁的dll实现读取酒店IC卡数据,直接用Readme里的方法出错. 函数声明: 一.读卡函数    ************************ Delphi 调用 ****** ...

  9. Asp.Net Core SignalR 用泛型Hub优雅的调用前端方法及传参

    继续学习 最近一直在使用Asp.Net Core SignalR(下面成SignalR Core)为小程序提供websocket支持,前端时间也发了一个学习笔记,在使用过程中稍微看了下它的源码,不得不 ...

  10. AI要被祭天了!删Bug,删着删着把全部代码都删了

    近日,美国版的“大众点评”,本想训练 AI 来消除 bug,结果它把所有内容删除了... Yelp 在其最新更新的 App 中写道: “我们为本周使用该app遇到问题的用户致歉.我们训练了一个神经网络 ...