计算传质学数值计算的实现.pptVIP

计算传质学——数值计算的实现 离散化的基础 例：一维常微分方程的离散和求解 二维微分方程离散和求解 在第(i,j)网格节点离散化后形式为： 离散格式 其他差分格式： 二阶导数的差分格式 流场计算的SIMPLE算法（Semi-Implicit Method for Pressure-Linked Equations） 考虑二维流场的计算，连续性方程和N—S方程： 连续性方程离散化 构造速度修正方程和压力修正方程： 假设压力场p*，求得的速度为U*，V*； U*，V*与真实的U，V之间的差距为U/，V/； SIMPLE算法中的重要近似： 速度修正方程为： SIMPLE算法的计算步骤： 假设初始速度分布和压力分布； 用当前速度和压力计算各个方程系数，求解动量方程，得到U*和V*； 用U*和V*求解压力修正方程，得到p·； 计算修正后的p,U,V：，，如果不收敛，回到步骤2。 蒸馏塔板的模拟 模型： 湍流：k-ε模型 传质：c2- εc模型 UDF编程 #include udf.h #define ct 0.11 #define c1 1.8 #define c2 2.2 #define c3 0.72 #define c4 0.8 #define us 0.5 #define rhog 1.27 #define rhol 998.2 #define g 9.81 #define h