计算流体力学-族谱新闻-族谱网

计算流体力学

2020-10-16

出处：族谱网

作者：阿族小谱

浏览:428次

转发:0次

评论:0

方法论所有这些方法都遵循同样的基本的程序。问题的几何（物理界限）被定义。流体占据的体积被分成离散胞腔（网格）。物理建模得到定义-例如，运动方程+焓+种类不灭边界条件被定义。这涉及到液体在问题的边界行为和性质。对于暂态问题，初始条件也要定义。方程作为静态或者暂态被重复求解。最后答案的分析和可视化。离散化方法给定离散化的稳定性通常在数值上建立，而不是像在简单的线性问题上那样可以解析的建立。必须特别小心才能保证离散化能够漂亮的处理不连续的解。欧拉方程和纳维-斯托克斯方程两者都可能有冲击波和接触表面。在使用中的一些离散化方法包括：有限体积法（FiniteVolumeMethod）。这是“经典”或者说标准的方法，在商用软件和研究用程序中最为常见。控制方程在离散的控制体积上求解。这个积分方法导致了一个本身就保守（也就是说，密度这样的量保持了物理上守恒）的方法。有限元法（FEM）。这个方法在固体结构分析...

方法论

所有这些方法都遵循同样的基本的程序。

问题的几何（物理界限）被定义。

流体占据的体积被分成离散胞腔（网格）。

物理建模得到定义 - 例如，运动方程+焓+种类不灭

边界条件被定义。这涉及到液体在问题的边界行为和性质。对于暂态问题，初始条件也要定义。

方程作为静态或者暂态被重复求解。

最后答案的分析和可视化。

离散化方法

给定离散化的稳定性通常在数值上建立，而不是像在简单的线性问题上那样可以解析的建立。必须特别小心才能保证离散化能够漂亮的处理不连续的解。欧拉方程和纳维-斯托克斯方程两者都可能有冲击波和接触表面。

在使用中的一些离散化方法包括：

有限体积法（Finite Volume Method）。这是“经典”或者说标准的方法，在商用软件和研究用程序中最为常见。控制方程在离散的控制体积上求解。这个积分方法导致了一个本身就保守（也就是说，密度这样的量保持了物理上守恒）的方法。

有限元法（FEM）。这个方法在固体结构分析中很流行，但是也可以用于流体。但是，FEM表述需要特殊处理来保证保守的解。

有限差分方法。这个方法有历史上的意义而且易于编程。现在只在特殊化的代码中使用。

边界元方法。流体占据的边界被分割成表面网格。

湍流模型

直接数值模拟（DNS）捕捉了所有相关尺度的湍流运动，所以不需要对最小的尺度建立另外的模型。该途径极其耗费资源，在现代的机器上对于复杂的问题甚至可能无法达成，因而对于流体运动最小尺度的模型还是有必要的。

雷诺平均纳维－斯托克斯方程

雷诺平均纳维－斯托克斯方程（RANS）是湍流最古老的方法。将流体运动的各物理量表达为一个统计平均值和一个脉动值的和，代入原始的N-S方程后再取统计平均，从而得到关于统计平均物理量的控制方程——RANS方程。在动量方程中，出现了一个新的应力项，它是脉动速度的二阶相关，称为雷诺应力。这个二阶张量是新的未知量，必须通过一定的模型使整个方程组得到封闭（即使得方程个数与未知量个数相等），才能求解。对于这个张量有很多不同的封闭模型，常见的如零阶模型，一阶模型和二阶模型。目前二阶模型，如k-epsilon模型和k-omiga模型，是工程中最为常见的湍流模型。

下面是常见的误解：RANS方程不适用于有一个时变平均流的流场，因为这些方程是‘时间-平均’的。事实上，统计不稳定（在非静态）流可以同等对待。这有时被称为URANS。没有任何雷诺平均方法中的先天缺陷使得这个情况不能被处理，但是用于闭合（封闭）方程的湍流模型必须要在均值中的变化发生的时间相对于包含大部分能量的湍流运动的时间尺度相对较大时才正确。

大涡流模拟

大涡流模拟（LES）是一种技术，其中更小的涡流被滤掉并用亚格点尺度模型来建模，而更大的能量的涡流则被模拟出来。该方法通常比RANS模型需要更细化的网格，但是比DNS解所需要的网格粗的多。

分离涡流模拟

分离涡流模拟（Detached eddy simulation，DES)是RANS模型的一个修改，其中模型在细到足够使用LES计算的区域切换到亚格点尺度表述。接近固体边界的区域和涡流长度尺度小于最大格点尺寸的区域被交给RANS模式的解。随着涡流长度尺度超出格点尺寸，区域使用LES模式。所以格点分辨率不像纯LES要求那样高，因此大大减少了计算开销。虽然DES最初为Spalart-Allmaras模型而建立（Spalart等人, 1997年），它可以用其它RANS模型实现（Strelets, 2001年），通过适当的修改RANS模型中隐含或显式地涉及到的长度尺度。所以虽然基于Spalart-Allamaras模型的DES像带墙模型的LES一样作用，基于其它模型的DES（譬如两方程模型）行为就像混合RANS-LES模型。格点的生成比简单的RANS或LES更为复杂因为要做RANS-LES切换。DES是一种非区域途径并会产生一个跨越解的RANS和LES区域的单一的光滑速度场。

求解算法

离散化完成后要进行的方程系统的基本求解由很多数值线性代数的为人熟知的许多算法完成。可以使用静态的迭代方法，譬如对称高斯-赛得尔（Gauss-Seidel）或者渐进松弛（successive overrelaxation），或者克雷洛夫子空间法（Krylov subspaces）。在后者中，解的余数在非线性算子的一个子空间的正交基上最小化。克雷洛夫子空间法通常和一个预条件算子（preconditioner）和一个内牛顿迭代一起使用。对于非线性问题不幸的是，正交基不能通过短循环来构造（相在一般的共轭梯度法中那样），而必须存储整个的序列的矢量。

应用

这些技术广泛用于工程设计或者分析和流体相互作用的设备，例如车辆，飞行器，泵，化学装置和通风系统等等。该技术也用于计算机图形学，因为动画家不能像控制固态人物那样简单的处理流体，因而必须使用CFD技术结合交互工具，达成动画影片或游戏中的流体模拟。

有许多求解纳维-斯托克斯方程的商用软件包。这些商用包的例子包括下面这些（按字母序排列）： AVL/FIRE、CFX、Fluent、KIVA、NUMECA、Phoenics、和STAR-CD。其它软件包作为附件或CFD工具的补充产品使用。这些包括用于后置处理的FieldView和求解细节化学动力学的KINetics。

计算流体力学知名专家

软件

CFD-ACE+多重物理耦合之计算流体力学商用分析软件，前身即为CFDRC。现在版权为法国ESI Group集团拥有。

OpenFOAM现在GPL版权的前商用软件。内部链接：OpenFOAM

Palabos现在GPL版权的前商用软件

NASA提供CFD软件给美国居民

Fire Dynamics Simulator消防工程应用软件

GERRISGPL版权不可压缩流体求解器。内部链接：Gerris.

管道系统液体流有一个功能完备的演示，它使用样品流体来做评估，可以下载使用

Atmos

Stoner流水线模拟器

流水线工作室包括一个评估版的下载（~ 132 M）