SIMPLE算法求解方腔内粘性不可压流动学习资料.doc

《SIMPLE算法求解方腔内粘性不可压流动学习资料.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、SIMPLE算法求解方腔内粘性不可压流动目录一、问题描述1二、离散格式2交错网格3方程离散4三、SIMPLE算法基本思想7边界条件处理8虚拟网格处理9方程求解10输出变量处理12SIMPLE算法流程图14四、程序中主要变量的意义14五、计算结果与讨论16函数最大值16变量等值线图17主要结论21六、源程序21一、问题描述假设的方腔内充满粘性不可压缩流体，左、右、下壁固定，上壁以运动，试求时的定常解，方腔如图1所示。图1方腔内流动示意图二、离散格式本算例采用求解不可压缩流动的经典算法，即SIMPLE算法，求解方腔内粘性不可压缩流体运动的定常解。SIMPLE算法的全称为S

2、emi-ImplicitMethodforPressure-LinkedEquations，即求解压力关联方程的半隐式算法。采用SIMPLE算法时，为了避免中心差分格式将“棋盘”型参量分布误认为是均匀分布，需要用交错网格对计算域进行离散。交错网格交错网格如图2所示，压力、密度等物理量存储在控制体的中心，这个控制体称为主控制体。速度分量分别存储在主控制体的和位置处，标记为位置，再分别以此为中心，划分速度分量u、v的控制体。采用空间均匀网格，等间距离散整个求解域，如图3所示。图2交错网格示意图图3求解域离散示意图图3中阴影部分代表方腔内的流动区域，阴影区域的边界代表方腔的

3、上、下、左、右壁面，阴影区域外面的网格节点是为边界处理需要而设定的虚拟网格节点，后面介绍边界处理方法时详细论述。方程离散无量纲化的守恒型不可压缩方程为其积分形式为图4主控制体图5速度u控制体图6速度v控制体采用有限体积法离散方程，连续性方程在主控制体上离散X方向动量方程在速度u控制体上离散，时间采用前差Y方向动量方程在速度v控制体上离散，时间采用前差其中，数值通量通量分别定义在主控制体的中心和角点，如图所示，并按照如下方法离散通量分别定义在主控制体的中心和角点，如图所示，并按照如下方式离散通量和的某些项冻结于M时间层，使离散化之后的方程对是线性的。将离散化之后的和代入

4、离散后的x方向和y方向的动量方程，整理之后得离散后的动量方程如下其中以上是SIMPLE算法中离散化的动量方程三、SIMPLE算法基本思想SIMPLE算法是一种解决压力-速度耦合问题的“半隐式”算法。首先给定M时刻猜测的速度场，用于计算离散动量方程中的系数和常数项。给定M+1时刻猜测的压力场估计值，迭代求解离散动量方程，得到M+1时刻速度场的估计值，速度场的估计值满足如下离散方程。一般地，速度场不满足离散的连续性方程，因而需要对速度场和压力场进行修正。M+1时刻的修正值和估计值有如下关系其中，和分别速度和压力的修正量，修正量亦满足离散的动量方程编号为（i,j）的速度修正

5、量不仅与压力修正量有关，还与邻近点的速度修正量有关。SIMPLE算法的重要假定：速度的改变只与压力的改变有关，忽略邻近点对速度修正的影响。因而得到如下速度修正量修正后的速度分量将修正后的速度分量代入离散后的连续性方程，得到压力修正方程其中采用迭代法求解压力修正方程，得到压力修正量，代入修正公式得到M+1时刻的速度场和压力场。将M+1时刻的速度场和压力场作为新的猜测的速度场和猜测的压力场估计值，采用上述方法计算下一个时刻的速度场和压力场，直到满足收敛条件。收敛判据为很小的正实数，视计算的精度要求而定。本算例中取。若，则，此时，从而来自于离散动量的满足离散的连续性方程。因

6、此可以作为收敛判据。边界条件处理首先对计算区域离散，并流动边界之外扩充一个虚拟网格，将真实流动的离散域包围，如图7所示。图7虚拟网格扩充示意图压力p的存储位置由图中符号代表，速度分量u的存储位置由图中符号代表，速度分量v的存储位置由图中符号代表。阴影区域表示真实的流动区域，阴影区域外的网格为虚拟网格。速度分量u在x方向虚拟网格上的存储位置分别记为，在y方向虚拟网格上的存储位置分别记为。速度分量v在x方向虚拟网格上的存储位置分别记为，在y方向虚拟网格上的存储位置分别记为。压力和压力修正量在x方向虚拟网格上的存储位置分别记为，在y方向虚拟网格上的存储位置分别记为。虚拟网格

7、处理静止壁面的虚拟网格按图8所示方式处理，可以满足粘性边界条件，以下边界为例说明。图8静止壁面虚拟网格处理示意图边界外的虚拟网格存储位置上的速度分量，由边界内的真实网格对应存储位置上的速度分量沿边界对称得到。这样可以保证虚拟网格与真实网格在对应存储位置上的速度分量大小相等，方向相反，二者在边界上的线性平均值为零，满足粘性边界条件。上壁面有速度，为运动壁面，速度分量。速度分量v采用前述对称方法处理。速度分量u的处理方法略微不同。上壁面速度不随时间变化，虚拟网格上速度分量由线性插值得到，可以保证上壁面满足边界条件，即。存储位置位于恰好位于静止壁面上的速度