流动问题场模拟数值计算主要包括以下步骤：建立物理与数学模型；选择坐标系；建立网格；选择建立离散方程的方法；选择对流项与扩散项的格式；边界条件的离散；选择流场计算方法；求解代数方程组；后处理与准确性判定。
1.建立物理与数学模型
首先，对所研究的实际问题作出一定的简化假设，以确立其物理模型。例如，当物理过程中流体的物性变化不大时可作常物性的假定，物理量的场在某一方向上变化相对于其他两个方向很小时可作2D假定等。考虑流动是多相流还是单相流动问题。大多数工程应用只需要考虑单相流动，而对于散热器晃动问题，需要考虑冷却液—空气的两相流动问题。在建立物理模型时，一般应考虑以下诸方面的因素。
（1）空间维数 空间维数是2D或3D。
（2）时间因素 时间是定常的还是非定常的。
（3）流动形态 流动形态是层流还是湍流。对于湍流，要选定相应的湍流模型。
（4）物性参数 物性为常物性还是变物性，可压缩的还是不可压缩的。
（5）过程类型 过程类型是抛物线型、双曲线型还是椭圆型。
（6）边界条件 边界条件是常规的一、二、三类边界条件，还是耦合的边界条件。
根据所确定的物理模型写出该过程的控制方程及相应的定解条件（初始条件及边界条件），建立数学模型。
2.选择坐标系
进行物理问题的数值计算时，最佳的坐标系是坐标轴与计算区域的边界相适应的坐标系。根据在空间任意一点上三个坐标面是否互相垂直，可区分为正交曲线坐标系与非正交曲线坐标系两大类。常用的正交曲线坐标系是：直角坐标系、圆柱坐标系、球坐标系。
3.建立网格
数值计算中用离散的网格代替原物理问题中的连续空间，网格中的节点则是所求解物理量的几何位置。从网格的构造来说，网络可以分为结构化（Structured）网格和非结构化（Un-structured）网格。
4.确定建立离散方程的方法
有限体积法保证了离散方程的守恒特性，物理意义明确，使用最广泛。有限体积法将守恒型的控制方程对区域离散后形成的控制体积作积分，对于节点间物理量的变化特性作出假设，从而得出节点间物理量间的代数方程式。
5.选择对流项与扩散项的离散格式(www.xing528.com)
在将控制方程对控制体积作积分的过程中，需要对所求解的变量在两个节点之间的变化特性作出假设，而不同的假设会导致不同的离散格式。在不同模型的流场计算中，代表具有扩散作用的二阶导数项采用具有二阶截差的中心差分形式，对流项采用不同的离散格式，非定常流动中的时间离散采用多步Runge-Kuta格式。
6.边界条件和初始条件的离散
初始条件和边界条件是控制方程有定解的前提，控制方程与相应的初始条件和边界条件的组合构成一个物理过程完整的数学描述，对初始条件和边界条件的处理，直接影响到计算结果的精度。
初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况，而边界条件是在求解区域的边界上所求解的变量或其导数随地点和时间的变化规律。对于一般性的开口计算区域，边界条件的类型有进口边界条件、固体边界条件、对称边界条件和出口边界条件四种。在CFD计算域内的流动是由边界条件驱动的，求解实际问题的过程就是将边界线或边界面上的数据，外推扩展到计算域内部的过程，因此提供符合物理实际且适定的边界条件是极其重要的，否则求解过程将很难进行。初始条件和边界条件的离散是将连续型的初始条件和边界条件转化为特定节点上的值。
7.选择流场计算算法
离散方程的数目与待求的未知量的数目相等，问题是封闭的，但是并非每一个因变量都有以之为函数的方程，如压力。在动量方程中，压力隐含在源项中，对速度场有很强的影响力。在具体计算中，压力场通过连续方程间接确定，当正确的压力场代入动量方程并解出速度场之后，该速度场应当满足连续方程；如果不满足，则需要对预先给出的压力场进行修正，直到解出的速度场满足连续方程为止，这就需要采用一些高级的算法来完成此项任务，即必须对离散方程进行某种调整，并且对各未知量（如速度、压力、温度等）的求解顺序及方式进行特殊处理。求解方法可分为耦合式解法和分离式解法两类。
基于原始变量的分离式解法的主要思路是顺序地、逐个地求解各变量代数方程组。目前使用最广泛的是1972年由Patanker和Splding提出的压力耦合方程组的半隐式算法SIMPLE。该算法已经很成熟，在应用上经过了很广泛的验证，其基本思想是：对于给定的压力场（该场可以是假定的值，或是上一次迭代计算所得到的结果），求解离散形式的动量方程，得出速度场。由于压力场是假定的或不准确的，这样，得到的速度场一般不满足连续方程，因此必须对给定的压力场进行修正。修正的原则是：与修正后的压力场相对应的速度场能满足本次迭代的连续方程。根据此原则，把由动量方程的离散形式所规定的压力与速度的关系代入连续方程的离散形式，从而得到压力修正方程，由压力修正方程得出压力修正值。根据修正后的压力场，求得新的速度场，然后检查速度场是否收敛。若不收敛，则用修正后的压力值作为给定的压力场，开始下一次迭代的计算，如此反复，直到获得收敛的解。
压力隐式算子分割法PISO是Issa于1986年提出的。PISO算法与SIMPLE算法的不同之处在于：SIMPLE算法是两步算法，即一步预测，另一步修正；而PISO算法增加了一步修正，即包含一步预测和两步修正，在完成了第一步修正得到速度、压力参数值后，寻求二次改进值，目的是使它们更好地同时满足动量方程和连续方程，PISO算法由于使用了一预测、一修正、再修正三步，从而可加快单个迭代步中的收敛速度。
耦合显式算法是由Fluent公司与NASA联合开发的，主要用来求解可压缩流动问题。该算法对N-S方程组进行联立求解，空间离散采用通量差分分裂格式，时间离散采用多步Runge-Kutta格式。该算法稳定性好，占用内存少，应用极为广泛。
耦合隐式算法是Fluent特有的算法。该算法也对N-S方程组进行联立求解，由于采用隐式格式，因此计算精度与收敛性要优于耦合显式算法，但此时占用内存较多。
使用哪种算法进行求解，要视具体问题而定。对于瞬态问题，PISO算法的优势明显；对于稳态问题，SIMPLE算法的优势明显。不可压缩的定常流动采用SIMPLE算法，非定常流动采用PISO算法；可压缩定常流动和非定常流动采用耦合隐式算法。例如，研究散热器晃动和列车远场空气动力噪声时，由于是瞬态流动，因此采用PISO算法；研究涡轮增压器叶轮和超声速拟似冲击波流动时，由于可压缩流动，因此采用耦合隐式算法。
8.求解代数方程组
将控制方程在所求解的网格上离散后，原问题将转化为如何求解规模巨大的代数方程组。方程组的阶数取决于所要求解问题的空间维度和离散网格的疏密程度，以及单个节点上所要求解的场变量个数。结构化网格上生成的代数方程组可以利用三对角矩阵算法（TD-MA）求解。当对流项采用高阶格式时，应采用五对角阵算法（PDMA）求解，从而提高收敛速度。对非结构化网格上生成的代数方程组求解时，一般采用点迭代法或共轭梯度法。
9.后处理与准确性评判
通过上述求解过程得出了各计算节点上的解后，可采用线值图、矢量图、等值线图、流线图、云图等方式将计算结果表示出来。
对数值计算的结果从物理过程的角度进行分析，并通过试验对数值仿真准确性进行评判。对非线性方程进行严密的数学求解是非常困难的，也是现代数学所致力于研究的。所以大量工程实际问题的数值计算结果仍属于“后验”的范畴，计算所得的结果是否可靠应得到实验的验证。
免责声明：以上内容源自网络，版权归原作者所有，如有侵犯您的原创版权请告知，我们将尽快删除相关内容。
