应用有限差分法计算2维欧拉方程.docVIP

基于非结构网格 二维Euler方程的Jameson求解方法 姓名：王司文 学号：sx1301102 摘要 本文介绍了基于CFD理论的求解二维可压缩流Euler方程的Jameson中心格式方法。在空间离散上采用的是有限体积法，时间上采用的是四步显式Runge－Kutta迭代求解。人工耗散项为守恒变量的二阶和四阶差分项。边界条件采用的是无反射边界条件，并采用当地时间步长进行加速收敛。最后对NACA0012翼型划分了三角形，并应用本文程序进行数值模拟，结果较为理想。 关键字：CFD，Jameson中心格式，Euler方程，有限体积法 Abstract A method for the numerical solution of the two-dimensional Euler equations has been developed. The cell-centred symmetric finite-volume spatial discretisation is applied in a general formulation. The integration in time, to a steady-state solution, is performed using an explicit, four-stage Runge－Kutta procedure. The artificial dissipation is constructed as a blending of second and fourth differences of the conserved variables. And in the boundary, there is none of the outgoing waves are reflected back into the computational domain. An acceleration technique called local time stepping is used. At last, standard test cases for both subsonic and supersonic flows have been used to validate the method. Key words：CFD, Jameson method, Euler equations, finite-volume 第一章 引言 在工程应用的推动下，计算流体力学随着计算机技术的发展和计算格式的不断更新而迅猛发展。在航空航天领域，CFD已经与地面试验和飞行试验共同构成了飞行器设计，飞行器性能分析和飞行器空气动力学设计的三大工具。由此可见，CFD数值模拟对于实际问题的求解是十分重要的。 在目前，实际流动问题的解多是通过求解Navier-Stokes方程来获得的，因为Navier-Stokes方程能够反映流体的质量守恒、动量守恒和能量守恒的规律。而由于本文针对的是二维无粘可压缩流，所以可以将Navier-Stokes方程简化为二维Euler方程，因此本文研究的是Euler方程的求解。 计算流体力学常用的求解方法有有限体积法、有限单元法和差分方法。本文所采用的求解方法是Jameson求解方法，该求解方法是中心格式的有限体积法。这种方法有两个特点：1、在空间的离散上，内场边的通量是通过将左右单元中心处的通量值取平均获得的；2、为了更好的捕捉间断点和提高格式的稳定性，在计算通量时加入了人工耗散项，这人工耗散项是由守恒变量的二阶和四阶差分项组成的。 Jameson格式可以方便地求解以非结构存储方式存储的网格。实际求解过程中，对初始化的全流场进行直接的时间积分直至收敛以获得稳定的数值解。本文的计算是基于非结构网格进行的，在时间上则通过四步Runge-Kutta迭代来获得定常解，边界条件则是采用无反射边界条件，并且采用了加速收敛技术。最后，分别对NACA0012翼型的跨音速和超音速绕流进行了数值模拟。 本文第二章是对Jameson求解方法的理论推导，第三章是提供算例和结果进行验证，第四章是总结与展望。 第二章 方法论述 本章具体介绍了二维Euler方程的Jameson求解方法，Jameson格式是一种有限体积的求解方法。这种方法就是对每个网格单元进行积分，然后通过计算得到流场结果的方法。为此，本章首先从二维Euler方程的控制方程着手，对其进行积分，得到积分方程。然后分别介绍了时间和空间上的离散，人工耗散项的计算，以及边界的处理。 2.1 控制方程 由第一章的介绍可知控制方程是Euler方程，对于二维的控制体的积分域为Ω，边界为S。则Euler方程可以改写为：