无网格！这是无网格CFD

计算流体动力学从业者最头疼的就是高质量分网问题，这通常占据整个仿真任务70%的时间。尽管在过去三十年算法的准确性和速度有了巨大进步，对复杂、真实世界模拟仿真来说分网还是最耗时和CFD仿真过程可靠度最低的部分。无网格方法因此提供了一个可行的替代基于网格的流体计算方法，并且不需要传统的网格结构这样就解决了很多分网相关的问题。下面就来介绍当前CFD中主流的无网格方法。

那么，网格是什么呢？

网格或格子定义为分析域或模型的离散单元格或单元。所有的流动变量和其他变量都在这些离散单元格中心求解。整个过程将物理域分解为更小的子域（单元/单元格）称为分网。这些单元格分组形成边界区域并且在这些区域施加边界条件。不仅仅产生高质量网格并保持它是一项繁重任务，还可能有其他影响诸如：

• 收敛速率

• 结果的准确性

• 所需CPU时间

尽管近年来有很多软件有自动网格划分功能但大多数CFD从业者还是手动进行网格划分。使用自动网格划分方法用户还是需要提供基本的输入，如单元尺寸，需要划分网格的区域，求解器来进行网格生成。可是这并非对所有案例都是可行的，对复杂几何模型难于实现自动。因此出现了CFD分析方法中的“无网格CFD”。

什么是无网格方法?

无网格方法用于建立整个问题域的代数方程系统而不使用域离散的预定义网格。问题域内节点是分散的，节点在边界上也是分散的节点组来代表（不离散）问题域及其边界。无网格意味着无需节点间关系的信息，与传统的有限体积或有限差分方法需要这种关系不同。

图形显示翼型划分了网格（左侧），右侧相同的翼型则被点包围。无网格方法在这些点处求解偏微分方程

为什么用无网格?

很多现有数值方法诸如有限体积方法，有限差分方法需要网格。在这样的网格中，每个点有一个固定数量的预定义相邻的点，与相邻点的连接关系用于定义数学运算诸如导数，并使用这个信息，求解整个域方程。

但对仿真材料可以移动或经历大变形（例如移动网格问题），网格的连接不引入误差则难以保持。尽管这种情况下仿真中可以重新划分网格，最终还会导致进一步增加误差。无网格方法可以避免这种问题。无网格方法的其他优势：

• 节省分网的时间

• 复杂几何模型的仿真非常容易，而复杂几何模型难于分网，可能花费几个星期来分网。

• 无需额外人员帮助划分网格。

不同的无网格方法:

近年来开发了很多无网格方法，我们来看一下发展的状况。

光滑粒子流体动力学方法(SmoothedParticle Hydrodynamics(SPH))

SPH，是用于天文学的最老的无网格方法之一，随后越来越多用于流体流动研究。这种方法，节点视作物理颗粒，带有质量和密度等信息可以随时间移动。这种方法中相邻颗粒任何属性值或其导数都是独立的。颗粒可以按任何顺序使用而没有影响，颗粒可以来回移动或甚至交换位置。该方法的域离散、场函数近似和数值求解基本步骤总结如下：

• 连续体分解为任意分布的一组没有连接（无网格）的颗粒；

• 采用场函数近似的积分表示法；

• 引入颗粒近似转换积分表示法为有限求和

径向基函数(Radial BasisFunctions(RBF)) 

RBF开发为一种求解偏微分方程无网格方法，源于认识到径向基函数插值对任何组节点任意维度都是光滑和准确的。RBF,其值依赖于到原点或任何其他指定点的距离的函数，来通过基于单一变量函数（径向基函数）线性组合近似多变量函数。通常用来仅知道有限数量点（或太难于评估的情况）近似函数或数据。RBF一些常用类型有：

• 高斯函数

• 多二次函数

• 逆二次函数

• 逆多二次函数

有限点集方法(Finite PointsetMethod(FPM))

FPM，是使用拉格朗日法的颗粒方法，流体被有限数量的颗粒（点）代替，是非静止颗粒。这些颗粒随流体速度移动承载流体量，诸如密度、速度、压力等。类似的边界可以通过有限数量的边界颗粒近似，边界条件施加在上面。如同SPH方法，FPM也不使用刚性相邻节点 /点（如有限体积方法FVM）列表。因此允许所有的点/颗粒移动，每个时间步长相邻点列表重新计算得到。

这种方法适于复杂几何模型，带有自由表面，多相流的流动问题。

FPM方法比广泛使用的无网格方法SPH有一些优势。SPH主要难点在于施加边界条件。FPM方法通过使用移动最小二乘或最小二乘法非常自然地把颗粒置于边界并施加边界条件而解决这个困难。

一些商业无网格CFD软件

让我们了解一些商业无网格CFD软件吧。

XFlow:

这是商业无网格软件之一，开发者是Next Limit Dynamics，现在属于达索系统。它使用格子玻尔兹曼方程及基于无网格颗粒的运动求解器。 XFlow能够求解下列问题：

• 移动边界问题

• 多相流

• 流固耦合

• 瞬态分析

• 大涡模拟

• 声学

• 非牛顿流体

风力涡轮机的CFD仿真(来源: XFlow)

XFlow仿真能够评估涡轮效率，预测叶片载荷，尾迹湍流强度，或风电场相互干扰。

NOGRID:

NOGRID, 始于2006年用于CFD分析。无网格软件NOGRID使用有限点集方法和N-S方程来求解CFD问题。这个软件的特色功能包括：

• 多相流

• 非牛顿流动

• 流固耦合

NOGRID软件在基于网格方法由于必须进行网格重划分而受限的情况非常有优势。使用快速稳健的NOGRID求解器通常的建模和计算时间可以显著缩短。

加速的船舶CFD仿真显示速度动画，使用NOGRID仿真的 (来源：NOGRID)

DualSPHysics:

这个软件基于光滑粒子流体动力学方法，开发用于研究欧拉方法遇到困难的自由表面流动现象，诸如水波或溃坝对海洋结构的影响。这是开源软件产品由约翰霍普金斯大学（美国），维戈大学（西班牙）和曼彻斯特大学（英国）研究者联合开发的。这个软件还可以使用GPU来进行求解。下面案例显示漂浮刚体相互作用问题。可以解决.

Algodoo:

Algodoo是由Algoryx Simulation AB开发的。这是一个二维仿真框架使用SPH方法主要用于教学目的的软件。易于使用可以创建良好的教育目的和学习物理原理的可视化图形。对学科教师和学生是一个好工具。

这个工具不像高端CFD求解器能够给出精确结果。但以非常少的时间给出可视化原理和进行很多不同仿真来说还是一个很好的软件。

除了流体仿真Algodoo还支持结构、光学的仿真。

无网格技术节省了大量的网格划分时间和精力。无网格CFD技术蓬勃发展未来前景光明，因为它解决了每个CFD工程师面对的最大困难之一—“分网”