本文用python3.6对一维激波管(sod、lax激波管)问题的几个经典算法进行了实现,包括Lax-Friedrichs、Lax-Wendoff、Mac-Cormack、Godunov、Roe等格式,激波管问题的初值采用了sod和lax激波管的初值(具体见代码),推进时间为0.2s,网格数视不同算法而定,参考书目见参考文献,其中水鸿寿、Toro、Laney的书帮助最大,具体代码参...
激波管是一种用于研究气体动力学的实验装置,其工作原理是通过活塞驱动驶入管道中的气体产生激波。此时,激波管内气体的参数会发生变化,需要通过一维Euler方程来描述。 解决一维Euler方程激波管问题的关键在于求解激波管内部气体的状态参数,包括气体密度、速度、压力和内能等。还需要考虑激波管两侧介质之间的压力、密度、速度...
让我们先来了解一下什么是一维激波管问题。一维激波管是由两个不同状态的流体组成的管道,其中一个是初始状态,另一个是激波通过后的状态。该问题的目标是确定激波传播的速度和管道中流体的各个参数的变化情况。 为了解决这个问题,我们可以使用Matlab编程来模拟激波管中的流体流动。首先,我们需要定义一些初始参数,如初始...
一维问题(激波管问题)数值解与其计算程序流程.doc,PAGE -A. 一维问题数值解与计算程序 一维问题,即激波管问题,是一个典型的一维可压缩无黏气体动力学问题,并有 解析解。对它采用二阶精度两步差分格式进行数值求解。同时,为了初学者入门和练习方便,这里给出了用语言和编
空间上 分别使用 0次 1次 2次正交多项式函数空间离散数值解时间上 分别使用 Euler向前离散, 二阶SSPRK 三阶SSPRK 限制器 TVB 限制器 数值结果代码DG code
一维激波管(Lax shock tube)问题的数值求解 文章目录 一维激波管(Lax shock tube)问题的数值求解 @[toc] 问题描述 一般的守恒算法格式 精确解 常用守恒格式(先介绍标量) Lax-Friedrichs格式 Roe格式 Lax-Wendroff格式 ENO格式简介(先介绍标量) 5阶WENO格式 空间离散 Lax-Friedrichs通量分裂 原函数方法重构左值右值....
0/0 收藏人数: 1 评论次数: 0 文档热度: 文档分类: 待分类 系统标签: 激波管sph模拟wij质点方程 [文章编号]10012246X(2003)0220153204一维激波管问题的SPH模拟徐立,孙锦山(中国工程物理研究院,北京研究生部,北京100088)[摘要][关键词]叙述SPH的基本原理,讨论了实现一维激波管问题SPH模拟的过程,并给出算例和说明...
常用激波管问题来检验 SPH = 1 [,] ()式对 求导得到场量的梯度(假设 可微) 方法和程序4 5 . 本文以模拟一维激波管问题为例, 6 ! W N ! ! A 详细介绍SPH 实现的过程. () ( ,) () A ! = m W ! ! I . 8 i $ = 1 1 SPH 基本原理 实践中有时为求更加精确,可将上式重新写成 ! !
处理好;2)能自然处理分裂成碎片的自由表面;3)容易加入如多相流、真实状态方程、可压缩性、辐射和固化等复杂物理现象;4)处理二维和三维复杂几何结构方便.另外,程序代码逻辑和结构简单.一维激波管问题因为存在数值间断,所以在数值模拟中具有特殊的意义.常用激波管问题来检验SPH方法和程序[4,5].本文以模拟一维激波管...
ENO格式在捕捉激波间断问题中需要计算3个模板多项式,而这些模板多项式仅在间断附近有较大差别,在其它较为平滑的位置差别很小.针对这一特点,对ENO格式进行优化,达到既节约计算时间又不影响计算效果的目的.之后,以一维激波管问题为例,在压强梯度大于预设临界值的位置用标准ENO格式计算,而在平滑的位置上,直接用预设模板...