该博客详细介绍了如何利用Matlab编程实现MacCormack方法来求解一维激波管问题。通过设置不同的初始条件，模拟了在激波管中高压和低压气体碰撞的过程，探讨了气体动态行为。文章包含关键的控制方程、边界条件、初始条件，以及采用的人工黏性项。并提供了主函数shockwave_1D及多个辅助函数的代码实现，用于数值求解和结果输出。最后展示了t=0.4s时的密度、速度和压强的计算结果。
激波管两端封闭 内部充满气体 直管中有一薄膜将激波管隔开 在薄膜两侧充有均匀理想气体
薄膜两侧气体的压力密度不同当t<0时 气体处于静止状态
当t>=0时薄膜瞬时突然破裂气体从高压端冲向低压端
控制方程有一维欧拉方程和理想气体方程：
边界条件：
初始条件：
MacCormack两步差分迭代格式：
人工黏性项：
----------------------------------------------------------------------------
用到七个函数
主函数 shockwave_1D
子函数
timestep、MacCormack、U2F、boundarycondition、output
用来后处理的函数 post
---------------------------------------------------------------------------
最低0.47元/天 解锁文章
扫一扫
1592
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
