[物理学与PDEs]第3章习题6 Lagrange 坐标下的一维理想磁流体力学方程组的数学结构
试讨论 Lagrange 形式下的一维理想磁流体力学方程组 (5. 33)-(5. 39) 的类型.
解答: 由 (5. 33), (5. 39) 知 $$\bex 0=\cfrac{\p p}{\p \tau}\sex{\cfrac{\p \tau}{\p t'}-\cfrac{\p u_1}{\p m}}+\cfrac{\p p}{\p S}\cfrac{\p S}{\p t'} =\cfrac{\p p}{\p t'}-p'(\tau)\cfrac{\p u_1}{\p m}, \eex$$ 而 $$\bex \cfrac{-1}{p'(\tau)}\cfrac{\p p}{\p t'}+\cfrac{\p u_1}{\p m}=0. \eex$$ 于是 (5. 33)-(5. 39) 为 $$\beex \bea \cfrac{-1}{p'(\tau)}\cfrac{\p p}{\p t'}+\cfrac{\p u_1}{\p m} &=0,\\ \cfrac{\mu_0}{\rho}\cfrac{\p H_2}{\p t'} +\mu_0H_2\cfrac{\p u_1}{\p m}-\mu_0H_1\cfrac{\p u_2}{\p m} &=0,\\ \cfrac{\mu_0}{\rho}\cfrac{\p H_3}{\p t'} +\mu_0H_3\cfrac{\p u_1}{\p m} -\mu_0H_1\cfrac{\ pu_3}{\p m} &=0,\\ \cfrac{\p u_1}{\p t'} +\cfrac{\p\rho}{\p m} +\mu_0H_2\cfrac{\p H_2}{\p m} +\mu_0H_3\cfrac{\p H_3}{\p m}&=F_1,\\ \cfrac{\p u_2}{\p t'}-\mu_0H_1\cfrac{\p H_2}{\p m}&=F_2,\\ \cfrac{\p u_3}{\p t'}-\mu_0H_1\cfrac{\p H_3}{\p m}&=F_3,\\ \cfrac{\p S}{\p t'}&=0; \eea \eeex$$ 其可化为 $$\bex A(U)\cfrac{\p U}{\p t'}+B(U)\cfrac{\p U}{\p m}=C, \eex$$ 其中 $$\beex \bea U&=(p,H_2,H_3,u_1,u_2,u_3,S)^T,\\ A(U)&=\diag\sex{\cfrac{-1}{p'(\tau)},\cfrac{\mu_0}{\rho},\cfrac{\mu_0}{\rho}, 1,1,1,1},\\ B(U)&=\sex{\ba{ccccccc} 0&0&0&1&0&0&0\\ 0&0&0&\mu_0H_2&-\mu_0H_1&0&0\\ 0&0&0&\mu_0H_3&0&-\mu_0H_1&0\\ 1&\mu_0H_2\mu_0H_3&0&0&0&0\\ 0&-\mu_0H_1&0&0&0&0&0\\ 0&0&-\mu_0H_1&0&0&0&0\\ 0&0&0&0&0&0&0 \ea},\\ C&=(0,0,0,F_1,F_2,F_3,0)^T. \eea \eeex$$ 故 Lagrange 形式下的一维理想磁流体力学方程组 (5. 33)-(5. 39) 是一阶对称双曲组.
[物理学与PDEs]第3章习题6 Lagrange 坐标下的一维理想磁流体力学方程组的数学结构的更多相关文章
- [物理学与PDEs]第3章习题5 一维理想磁流体力学方程组的数学结构
试将一维理想磁流体力学方程组 (5. 10)-(5. 16) 化为一阶拟线性对称双曲组的形式. 解答: 由 (5. 12),(5. 16) 知 $$\beex \bea 0&=\cfrac{\ ...
- [物理学与PDEs]第3章习题参考解答
[物理学与PDEs]第3章习题1 只有一个非零分量的磁场 [物理学与PDEs]第3章习题2 仅受重力作用的定常不可压流理想流体沿沿流线的一个守恒量 [物理学与PDEs]第3章习题3电磁场的矢势在 Lo ...
- [物理学与PDEs]第2章习题11 Lagrange 形式的一维理想流体力学方程组在强间断线上的间断连接条件
对由第 10 题给出的 Lagrange 形式的一维理想流体力学方程组, 给出解在强间断线上应满足的间断连接条件 (假设体积力 $F\equiv 0$). 解答: $$\beex \bea \sez{ ...
- [物理学与PDEs]第1章习题12 Coulomb 规范下电磁场的标势、矢势满足的方程
试给出在 Coulomb 规范下, 电磁场的标势 $\phi$ 与矢势 ${\bf A}$ 所满足的方程. 解答: 真空中的 Maxwell 方程组为 $$\bee\label{1_10_12:eq} ...
- [物理学与PDEs]第1章习题3 常场强下电势的定解问题
在一场强为 ${\bf E}_0$ (${\bf E}_0$ 为常向量) 的电场中, 置入一个半径为 $R$ 的导电球体, 试导出球外电势所满足的方程及相应的定解条件. 解答: 设导电球体为 $B_R ...
- [物理学与PDEs]第1章习题参考解答
[物理学与PDEs]第1章习题1 无限长直线的电场强度与电势 [物理学与PDEs]第1章习题2 均匀带电球面的电场强度与电势 [物理学与PDEs]第1章习题3 常场强下电势的定解问题 [物理学与PDE ...
- [物理学与PDEs]第2章习题参考解答
[物理学与PDEs]第2章习题1 无旋时的 Euler 方程 [物理学与PDEs]第2章习题2 质量力有势时的能量方程 [物理学与PDEs]第2章习题3 Laplace 方程的 Neumann 问题 ...
- [物理学与PDEs]第4章习题参考解答
[物理学与PDEs]第4章习题1 反应力学方程组形式的化约 - 动量方程与未燃流体质量平衡方程 [物理学与PDEs]第4章习题2 反应力学方程组形式的化约 - 能量守恒方程 [物理学与PDEs]第4章 ...
- [物理学与PDEs]第5章习题参考解答
[物理学与PDEs]第5章习题1 矩阵的极分解 [物理学与PDEs]第5章习题2 Jacobian 的物质导数 [物理学与PDEs]第5章习题3 第二 Piola 应力张量的对称性 [物理学与PDEs ...
随机推荐
- log4j控制指定包下的日志
最近观察日志发现如下两个问题: 1.项目用的是springboot项目,整合了rabbitmq,项目启动后,会自动监控rabbitmq谅解是否正常,导致控制台一直输出监控日志,此时就想阻止该类日志输出 ...
- GIL:全局解释器锁 VS 用户程序锁
既然有了GIL锁,CPython还要多线程干什么? ''' GIL:全局解释器锁的来历 四核:同一时刻真正有四个任务在运行,多核的意义在于此 单核:看上去是并发的,因为进行了上下文切换,单核永远是串行 ...
- day19-网络编程基础(二)
今天没有很多概念性的东西,主要是方法性的东西以及编程的一些方法吧 今日份目录 1.UDP传输的特点以及实验 2.UTP与UDP传输的区别 3.基于tcp的low版带验证功能的FTP小程序 4.基于so ...
- Django-CRM项目学习(三)-stark的页面展示
开始今日份整理,主要为stark组件的增删改以及model_form组件的使用以及form组件的回顾 1.list_display_link功能 功能:客制化字段进行跳转到指定的页面(编辑页面) 1. ...
- Jenkins + Ansible + Gitlab之ansible篇
Ansible介绍 什么是Ansible? Ansible是一个开源部署工具 开发语言:Python 特点:SSH协议通信,全平台,无需要编译,模块化部署管理 作用:推送Playbook进行远程节点快 ...
- linux目录1
1.linux常用快捷键 2.linux命令之vmstat 3.linux命令之df 4.linux命令之netstat 5.linux命令之 tar 6.The authenticity of ho ...
- mac 利用svn下载远程代码出现Agreeing to the Xcode/iOS license requires admin privileges, please re-run as root via sudo.
终端输出的信息:Agreeing to the Xcode/iOS license requires admin privileges, please re-run as root via sudo. ...
- warn_alloc():page allocation failure问题分析
关键词:warn_alloc().__GFP_XXX.order.CMA等等. 在内存申请的时候经常会遇到类似“ xxx: page allocation failure: order:10...”类 ...
- [Spark]如何设置使得spark程序不输出 INFO级别的内容
Spark程序在运行的时候,总是输出很多INFO级别内容 查看了网上的一些文章,进行了试验. 发现在 /etc/spark/conf 目录下,有一个 log4j.properties.template ...
- TypeError: 'Item' object has no attribute '__getitem__'
Error Msg: Traceback (most recent call last): File "start.py", line 8, in <module> E ...