[发明专利]一种应用于高精度间断迦辽金流体仿真的粘性项处理方法

说明书

本发明公开了一种应用于高精度间断迦辽金流体仿真的粘性项处理方法，用于解决迦辽金方法在计算Navier‑Stokes方程时出现的计算精度低于理论精度的问题，从而精确捕捉流场中用于工程应用的升力、阻力、速度、密度和压力等信息。包括如下步骤：对空间进行剖分得到计算网格，采用高精度间断迦辽金框架对Navier‑Stokes方程进行处理，得到半离散形式的控制方程；定义新的提升算子，采用新的提升算子计算粘性通量，积分后获得粘性项的计算结果；采用迭代方式进行方程的求解计算，获得仿真结果。本发明在节省计算量的同时，有效保持了高阶格式的计算精度，计算精度高于理论精度。

技术领域

本发明涉及计算流体力学技术流域，具体涉及到一种应用于高精度间断迦辽金流体仿真的粘性项处理方法。

背景技术

高精度间断迦辽金仿真计算作为计算流体力学的重要分支，具备良好的色散和耗散特性，适合求解多尺度复杂问题，可以利用较少的计算资源捕捉精细的流场计算信息，因此得到了越来越多学者的高度重视。

高精度间断迦辽金方法最早只能求解含有一阶偏导数的双曲守恒型方程，如Euler方程，但是流体仿真面临的大多数问题都是基于Navier-Stokes方程，为了求解Navier-Stokes方程，并获取流场中的升力、阻力、速度、密度和压力等信息，可以通过引入Navier-Stokes方程中的控制变量的梯度和梯度的导数作为中间变量，将原始的二阶方程降为两个一阶系统，然后按照间断迦辽金框架进行离散处理。传统方法需要计算和存储变量梯度以及梯度的导数，在增加计算方程个数和实际计算量的同时，会降低高精度情况下的计算精度。

发明内容

本发明的目的是通过采用新的提升算子的方式，有效的去除了粘性项计算中梯度以及梯度导数的计算过程，不仅节省了计算量，在提高效率的同时，可以有效的保持高阶格式的计算精度。

为了实现上述目的，本发明采用如下技术方案：

一种应用于高精度间断迦辽金流体仿真的粘性项处理方法，采用高精度间断迦辽金框架对空间进行剖分生成计算网格，通过迭代积分计算粘性通量，生成获取升力、阻力、速度、密度和压力的网格，具体步骤如下：

步骤一：对空间进行剖分得到计算网格，采用高精度间断迦辽金框架对Navier-Stokes方程进行处理，得到半离散形式的控制方程。

步骤二：定义新的提升算子，采用新的提升算子计算粘性通量，积分后获得粘性项的计算结果。

步骤三：采用迭代方式进行方程的求解计算，获得仿真结果。

与现有技术相比，本发明的有益效果是：

1、相比于以往的方法，本方法只需计算提升算子，就可以获得粘性项的计算结果，省去了梯度和梯度导数的计算，计算过程简化，计算效率提高。

2、相比以往的方法，保持了高阶格式的计算精度，整个计算的精度相比以往方法有较大提升。

附图说明

本发明将通过例子并参照附图的方式说明，其中：

图1是本发明的流程示意图；

图2是采用本发明方法得到的密度等值线示意图；

图3是采用传统方法得到的密度等值线示意图；

图4是本发明方法与传统方法的误差和精度对比的对比表。

具体实施方式

本说明书中公开的所有特征，或公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合。

本说明书（包括任何附加权利要求、摘要和附图）中公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换。即，除非特别叙述，每个特征只是一系列等效或类似特征中的一个例子而已。