CFD经验-流沙.docVIP

CFD为何物？CFD其实只是一个缩写，它的全称是Computational Fluid Dynamics，翻译成中文为“计算流体动力学”，国内习惯称之为计算流体力学。在这里可以分析一下它的名字。首先它是流体动力学。我们知道，流体力学分为静力学和动力学两种。静力学研究的是流体静止条件下的状态，而动力学则表示流体在流动过程中的运动状态。既然是流体动力学，当然研究的是流体的运动。流体运动起来了之后，由于我们的研究对象处于宏观状态，所以可以应用牛顿第二定律。其实NS方程都可以从牛顿第二定律推导。守恒定律是普遍通用的，通过从质量守恒、能量守恒、动量守恒的角度对流体微团进行分析，可以从流体流动过程中抽象出数学方程。 那么如何计算呢？本来CFD的控制方程是封闭的，5个方程（1个质量守恒方程、3个动量守恒方程、1个能量守恒方程），5个未知数（u、v、w、P、T），理论上是可以直接离散形成封闭的代数方程组进行求解，然而不幸的是，这5个方程是高度非线性的，再加上流体流动中形成的湍流涡尺度相差过大（最小涡和最大涡尺寸相差太大），因此如果直接对控制方程离散求解的话，则要求模型网格尺寸足够的小（小于最小涡尺寸，现实中最小涡尺寸非常小），这样无疑极大的增加计算量，以至于无法在工业应用中进行推广。于是人们想出了采用雷诺平均NS方程（RANS），将湍流中的速度脉动进行平均。由于引入了新的变量，结果导致控制方程不封闭。人们为了使控制方程重新封闭，引入了湍流模型，湍流模型虽然是对湍流状态的一种假设，但是的确解决了直接模拟过于浪费计算时间的问题，因此在工业上得到了广泛的应用。 最主要的优势莫过于CFD能大幅降低试验成本。CFD通常都是利用计算机进行模拟，因此减少了试验设备需求。更突出的优势在于，CFD中的模型可以根据需要随时修改。另一个不得不提的优势在于，很多不可重现、或重现代价非常大的情况，可以用CFD进行再现。如火灾、爆炸、一些自然灾害等。 CFD计算的准确性还存在很大的局限。首先它受制于当前计算机的计算能力，虽然如今计算机技术发展迅速，计算能力呈指数增长，但是我们对现实世界的模拟要求的细节也是越来越精密。要得到高精度的计算结果，当然对细节的考虑是越多越好，这无疑极大的增加了计算量。其次，数学理论的研究也在一定程度上制约了CFD的精度，高阶离散算法的收敛性、稳定性都不太容易解决。 CFD应用的致命伤：精度与可信度 CFD若要想在工程中得到广泛的应用，必须克服两大难点：准确性与可信性。在工程上，尤其是一些关键的工程中，谁也不敢轻易的应用一些精度与可信度得不到保证的数据。与固体应力计算使用有限单元法不同，目前主流的CFD软件几乎都是采用的有限体积法（除了CFX采用混合有限元法与有限体积法外，FLUENT、STAR-CD、Phonecis、Flow-3D等都是采用的有限体积法）。在计算量上来说，相同网格数量的模型，有限体积法消耗的内存要少于有限元法。在有限单元网格中存在的高次单元，其单元节点位于网格边的中点及网格体的中心，但是有限体积法中的高阶格式，其并非在网格单元中添加节点，而地更多的利用周围的节点。但是正因为如此，有限体积法计算精度要低于有限元法（在相同网格数量情况下）。影响CFD计算精度及可信度的原因自然不可能全怪罪于算法，更多的是问题存在于使用者及客观环境。CFD软件是一个黑盒子，利用CFD软件解决工程问题，软件使用者对于数据流向并不清楚，实际上对于非CFD专业的人事来说，也不必完全清楚CFD的内部运作方式，但是，如何有效的利用当前的软件，如何最大限度的发挥当前软件的计算性能，将计算结果精度及可信度提高，仍然是非常重要的，也是每一个从事CFD工程应用的人必须注意的。最需要注意的部分包括下面一些内容。一、精度1、算法导致的精度问题。一般来说，高阶算法的精度要高于低阶精度。但是收敛性却相反，采用高阶算法要比低阶算法收敛更困难一些。在一些高速流动情况中，采用迎风格式比中心差分格式能更好的收敛，在扩散占优的流动中则相反。以FLUENT为例，其具有一阶迎风格式与二阶迎风格式、幂律格式、QUICK格式以及三阶MUSCL格式。通常一阶迎风格式用于初步求解，较少用于最终计算结果的获得。QUICK格式在结构网格中具有三阶精度且收敛性较好，但是在非结构网格中只有二阶精度。二阶迎风格式在实际工程中用得非常多。三阶MUSCL格式用得较少，收敛性不是很好。2、边界条件。边界条件会对计算结果产生本质影响。也就是说，不确切的边界条件会导致不确切的计算结果，错误的边界条件一定得不到正确的结果。在实际工程中，能做为边界的位置的信息一定是确切的。换句话说，用户使用不确切的边界值，责任不在软件。当然有时候受条件限制得不到边界物理量，但是，软件的使用者应当对自己输入的边界值负责。其实相比较