论文部分内容阅读
在磁流体力学数值模拟中有许多数值算法，其中迎风格式是根据方程组特征波的传播方向计算通量，能够正确的反映物理信息的传播方向，获得较好的间断捕捉效果。迎风格式可以分为FDS（Flux Difference Splitting）和FVS(Flux VectorSplitting)两类，本文所采用的AUSM(Advection Upstream Split Mach number三种格式属于结合FDS和FVS的混合通量分裂格式。AUSM算法的核心思想是把方程组的流通量分为流动项和压力项，然后根据网格界面的马赫数计算流通量，不用计算通量Jacobian矩阵的特征值和特征向量，方便易行。 本文通过Ryu and Jone激波管、rotor问题、Orzag-Tang vortex问题三个算例的对比，对AUSM、AUSM+、AUSMPW+三种数值格式的优劣进行评价。AUSM+比AUSM和AUSMPW+的稳定性好，AUSMPW+捕捉间断的能力强，但是振荡太大。同时，也对比了HDC（Hyperbolic Divergence Cleaning）和ACA(ArtificialCompressibility Analogy)两种磁场散度处理方法在AUSM三种算法下的作用，发现ACA和算例本身有关系、人为因素很大，HDC优于ACA。最后本文选择较为合适的AUSM+算法，运用六片网格技术，使用多级场和parker解作为输入初步应用到日冕结构的数值模拟上，模拟的结果能够较好的再现日冕结构，得到日冕的大致形态。同时HDC方法能够较好的控制磁场散度误差，具备在空间物理问题中进一步应用的价值。