第六章算法应用.ppt

TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 对流方程式 差分方程式 对流出项采用隐式法 TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 对流方程式 差分方程式 对流出项采用隐式法 von Neumann 稳定性分析 结论 c=u△t/△z，对于所有的Θ都是|C|1，因此方程无条件稳定 TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 预测步 动量守恒方程的对流项 质量、能量守恒方程的对流项 TRAC-BF1的方法 SETS法的简化法 对流项使用施主元差分、流出项隐式处理 方法稳定，但是不满足守恒性 在修正阶段为满足质量守恒定律就要修正空泡份额 方法特征 修正步 预测步 为了满足质量守恒，修正空泡份额。对汽相、液相的质量守恒方程进行了隐式法差分 动量守恒方程的对流项 质量、能量守恒方程的对流项 隐式法 3 半隐式法 2 显式法 1 半隐解法 把相变化项与压力传播项这种时间尺度短的项进行隐式差 把对流移动的项进行显式差分的方法 时间步长受到与对流移动相关的Courant条件的限制 缺点 半隐解法 把相变化项与压力传播项这种时间尺度短的项进行隐式差 把对流移动的项进行显式差分的方法 时间步长受到与对流移动相关的Courant条件的限制 缺点 缺点 全隐解法 把所有的项都进行隐式差分的一种方法 放宽了对时间步长的限制，能够取相对较大的时间步长 使得编程变得复杂，增大了一个时间步长的计算量 MINCS 程序：日本核能研究所开发 隐式解法 基于Newton法的解法 基于逐次代入法的解法 CATHARE 程序：法国CEA、EDF以及FRAMATOME共同开发 PHOENICE 程序：英国Spalding开发 1.基于Newton法的解法 解法中差分方程式一律采用Newton法进行求解 采用两流体模型，空间网格数是N时，X和F是6N维的向量。J是Jacobi矩阵 通常，迭代的初值用前一个时间步的解Xn，求收敛时的Xn+1。 1.基于Newton法的解法 解法中差分方程式一律采用Newton法进行求解 采用两流体模型，空间网格数是N时，X和F是6N维的向量。J是Jacobi矩阵 优点：单管的时候的Jacobi矩阵是一个以6行6列的子矩阵为单位的组合三对角矩阵 缺点：Jacobi矩阵的形式复杂（网格数多的时候），求解需要很多的计算时间；解法用在多维问题的求解中是有困难的 通常，迭代的初值用前一个时间步的解Xn，求收敛时的Xn+1。 2.逐次代入法 逐次代入法 PHOENICE 程序 典型例子 两相流分析领域 SIMPLEST-ANL法 IPSA法 美国阿贡国家实验室开发的COMMIX-1A程序中采用的 2.逐次代入法 求解思路 MINCS 程序：求解一维两相流分析模型基本方程和结构关系式的程序 两流体模型 有限差分 全隐式处理 Newton法求解 1.差分式 MINCS程序的空间离散化的方法以交错网格法的有限差分法为基础 对流项使用施主元法进行差分化 运动量的对流项近似于下面所示的迎风差分和中心差分的一次结合 1.差分式 守恒定律 质量守恒 动量守恒 能量守恒 2.差分式的解法 守恒方程 2.差分式的解法 状态方程 守恒方程 引入 引入 2.差分式的解法 状态方程 守恒方程 引入 引入 2.差分式的解法 状态方程 守恒方程 引入 引入 Newton法的迭代求解 谢谢 期待您的建议! 求解速度 算法求解思路 求解内能 迭代求解压力与含气率 获得新值 半隐式法 2 显式法 1 隐式法 3 半隐式法：除了对流项使用显式求解以外，其他的相间的相互作用以及压力传播项全部采用隐式法进行求解 半隐式法：除了对流项使用显式求解以外，其他的相间的相互作用以及压力传播项全部采用隐式法进行求解 时间步长理论上只受伴随着流体流动的Courant条件限制 典型例子 SMAC法 SOLA法 SOLA-ICE ICE法 半隐式法：除了对流项使用显式求解以外，其他的相间的相互作用以及压力传播项