计算流体力学

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

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

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

${\displaystyle {\frac {\partial }{\partial t}}\iiint QdV+\iint Fd\mathbf {A} =0,}$

其中${\displaystyle Q}$是被守恒的变量的向量，而${\displaystyle F}$是通量的向量，${\displaystyle V}$是控制体积元体积，${\displaystyle A}$是控制元表面积。

• 有限元法（Finite Element Method, FEM）。这个方法在固体结构分析中很流行，但是也可以用于流体。但是，FEM表述需要特殊处理来保证保守的解。
• 有限差分 (Finite Difference Method, FDM)方法。这个方法有历史上的意义而且易于编程。现在只在特殊化的代码中使用。
• 边界元方法。流体占据的边界被分割成表面网格。

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

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

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

大涡流模拟（Large eddy simulation, 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）和一个内牛顿迭代一起使用。对于非线性问题不幸的是，正交基不能通过短循环来构造（相在一般的共轭梯度法中那样），而必须存储整个的序列的向量。

过去幾十年間，计算流体力学進展快速，除了應用在空氣動力學、渦輪、航海工程、航空載具、氣象學模擬、環境工程學等重要科學領域；在过去的几年中，計算流體力學在建筑设计中越来越重要。计算流体力学（和热傳學）是以電腦計算纳维尔-斯托克斯方程等能量、動量方程，以求得流动系统或热系统的控制微分方程 [1]，可透過有限元素法或有限體積法求解。这些技术广泛用于工程设计或者分析和流体交互作用的设备，例如车辆，飞行器，泵，化学装置和通风系统等等。该技术也用于计算机图形学，因为动画家不能像控制固态人物那样简单的处理流体，因而必须使用CFD技术结合交互工具，达成动画影片或游戏中的流体模拟。

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

计算流体力学可用於分析建築物廢氣對環境的影響，預測建築物中的煙霧和火災風險、量化室内環境品質以及设计自然通风系统等方面，這對於建築性能模擬有重要幫助[2]。建築性能模擬的目的是量化与建筑物的設計、建造、營運和控制等各個方面[2]，而计算流体力学正好提供具有科學根據的評估與計算。諸多建築能耗模擬軟體如ESP-r[3]、ENVI-met[4]、TRNSYS[5]都有求解建築計算流體力學問題之功能。

以前，大多数与建筑有关的流體力學问题，例如通風分析，風荷载，風环境等，都是使用风洞测试进行检查的，但如今，所有这些测试都可以使用计算流体力学有效地完成。计算流体力学可以在相對較短的时间内解决所有上述问题，并且與較舊的方法（實驗方法）相比，它更經濟，更有效。[6]计算流体力学可以用于预测建筑物内部和周围的空气流动，热传递和污染物输送。 计算流体力学可以輕鬆更改不同组件的类型和位置，送风条件和系统控制时间表，从而检查各种供暖通空調（HVAC）系统的有效性和效率。此外，计算流体力学通过建模和最佳化建筑场地平面图和室内布局，帮助制定被动的加热/冷却/通风策略（例如自然通风）[7]。 在全球范围内，建筑部门约占总能耗的40％[8]，加強建築節能設計對於總體節能有巨大效益。

由于建筑部门占总消耗量的巨大比例，因此有必要研究建筑物的最佳配置以减少建筑物的能源消耗。为了实现这一目标，计算流体力学可以发挥重要作用。建築性能模擬和CFD程序是重要的建筑设计工具，可用于评估建筑性能，包括热舒适度，室内空气质量机械系统效率和能耗[9]。

建筑物中的计算流体力学主要用於建築通风分析。

通过调节室内空气参数（空气温度，相对湿度，空气速度和空气中的化学物质浓度）来进行建筑物的通风研究，以找到室内温度可接受的热舒适环境。 计算流体力学可有效計算室内空气流動相關参数，用以预测建筑物的通风性能[10]。

这些空气参数对于设计舒适的室内以及室外环境中建筑物的良好整合至关重要。这是因为适当的通风系统的设计和控制策略的开发需要有关以下参数的详细信息；

• 空气流动
• 污染物扩散
• 温度分布

图1（a）：建筑物周围的流动（側視圖）

图1（b）：建筑物周围的流动（俯視圖）

上述資訊对于建筑师设计建筑結構與建築外殼也是有相當大的幫助。隨著電腦技術的進步，计算流体力学技术在建筑物中得到了广泛的成功应用[11]。

最近，通风及其相关领域已成为风能工程的重要组成部分。可以使用风洞调查（实验性地）或通过计算流体力学建模（理论上）来进行通风研究。在某些应用中，自然通风系统可能比强制通风系统更可取，因为它消除或减少了机械通风系统，从而既可提供风扇能量又可节省第一笔费用。在当今时代，由于开发了许多计算流体力学软件和其他建筑物性能模拟软件，因此评估建筑物中自然/强制通风系统的可能性变得更加容易。 计算流体力学分析比实验方法非常有用，因为可以在后处理中找到变量之间的其他相关关系。通过实验或数值获得的数据可通过两种方式使用：[12]

1. 更好的用户舒适度
2. 它提供了用作建筑物热量平衡计算输入的数据

在选择建筑场地和位置时，风荷载起着重要作用。如果两座建筑物在某个位置并排存在并且有一定的间隙，那么当通过该间隙吹动建筑物两端的风量为时，在第一瞬间分别绕过每座建筑物的流量总和必须大于或等于围绕单个建筑物的末端以压力损失为代价。

因此，将有压力的建筑物进入间隙，这将导致建筑物侧面的风荷载更高。当风吹过高层建筑的表面时，前面的向下流动会产生涡流（如图1所示）。靠近地面的反方向风速可能是参考风速的140％。因此，如果在该区域中存在任何建筑物，则该建筑物可能会受到损坏（尤其是建筑物的屋顶可能会受到严重损坏）。如果在建筑物的建造初期考虑到风荷载的影响，就可以成功地防止对建筑物的这种破坏。在建筑的早期，所有这些风荷载效应都是通过风洞试验确定的，但是今天所有这些试验都可以通过CFD分析成功地进行。为建筑物提供宜人的环境的重要性正在增加。经常要求建筑师和风力工程师在建筑物的形成阶段和施工的规划阶段仔细检查设计（方向，位置，位置以及周围建筑物之间的间隙）。 [12]因此，通过计算流体力学分析，可以找到适合建筑物方向，位置和局部最佳化問題（局部风速，对流系数等）。