非结构网格Euler方程解的有限体积法

摘要

非结构网格一般采用三角形（2-D）或四面体（3-D）单元。当前常用的方法主要有叉树法，推进面法，Delaunay法。该文中采用的是Delaunay方法。Euler方程的求解，本文采用有限体积法。对于由有限体积法得到的半离散化方程，可采用显式格式和隐式格式进行求解。显式格式常采用的是Runge-Kutta显示时间多步推进格式，该文采用五步显式格式。隐式格式求解常采用标准迭代方法，如Richardsons迭代，Jacobi迭代Causs一Seidel迭代等，GMRES---广义最小余量怯。为了加速收敛性采用预条件GMRES方法。