基于非结构网格二维欧拉方程的求解方法

基于非结构网格

二维Euler方程的Jameson

求解方法

姓名：

学号：

摘要

本文介绍了基于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 方程可以改写为：

⎰⎰Ω=-+Ω∂∂

s

Gdx Fdy Wd t 0)( （2-1） 其中：X,Y 是笛卡尔坐标系，W 为守恒变量矢量：

⎥⎥⎥⎥

⎦

⎤⎢⎢⎢⎢⎣⎡=E v u W ρρρρ （2-2） F,G 为流量矢量：

⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=UH UV P U U F ρρρρ2 ⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎣⎡+=VH P V UV V G ρρρρ2 （2-3） ρ，P ，H 和 E 分别为密度，压强，单位体积的总焓和单位体积的总能。U,V 分别为速度在X,Y 方向的分量。且对于理想气体，总焓和总能可以表示为：

2/)()1/(22V U P E ++-=ργρ （2-4）

P E H +=ρρ （2-5）

将计算域划分为有限个互不重叠的单元，并将积分守恒方程应用于每个单元，由于各个单元的面积不随时间变化，所以可将（2-1）式改写为：

⎰⎰Ω

Ω--=∂∂d Gdx Fdy t W S )( （2-6） 该方程对空间的离散是应用有限体积方法，然后对得到的半离散方程（时间的离散并未完成）按时间步长推进从而得到精确解。

2.2 空间离散

在进行时间离散之前，对方程（2-6）空间离散得到一个关于时间的常微分方程：

k k k Q dt

dW Ω-=/ （2-7） 其中：Ωk 是第k 个单元的面积，W k 是守恒变量矢量，k Q 是通量积分的离散近似值，可写成：

∑=∆-∆=

kedges

i i k x G y F Q 1)( （2-8）

其中： a b i x x x -=∆a b i y y y -=∆ （2-9）

这是第k 个单元所有边界的累计总和。采用中心格式的有限体积法，把守恒通量都控制在单元中心上，流体流经第i 条边的值由两个相邻单元中心点（k 和p ）的平均值决定：

2/)(p k i W W W += （2-10）

下标的意义如下图：