计算流体力学有限元方法及其编程详解

序
前言
第1章 有限元方法的基本思路
1.1 一维一次常微分方程的有限元数值求解
1.1.1 方程及精确解
1.1.2 有限元方法求解
1.1.3 程序编写
1.2 一维二次常微分方程的有限元数值求解
1.2.1 方程及精确解
1.2.2 有限元方法求解
1.2.3 程序编写
第2章 理想流体势流的有限元求解
2.1 求解实例和数学方程
2.1.1 求解实例
2.1.2 理想流体的特征及流动方程
2.1.3 边界条件
2.2 有限元方法求解二维拉普拉斯方程
2.2.1 计算区域的离散
2.2.2 插值函数及相关计算
2.2.3 加权余量方程
2.2.4 单元方程的建立
2.2.5 总体方程的组合
2.2.6 编程计算流程
2.3 程序编写
2.3.1 主程序
2.3.2 网格划分程序
2.3.3 网格图形显示程序
2.4 计算结果
第3章 牛顿流体流动的有限元求解
3.1 求解实例和数学方程
3.1.1 求解实例
3.1.2 数学方程
3.1.3 边界条件
3.2 速度?压力有限元求解
3.2.1 计算区域的离散
3.2.2 插值函数及其相关计算
3.2.3 加权余量方程
3.2.4 单元方程的建立
3.2.5 总体方程的组合
3.2.6 求解流程
3.3 速度?压力有限元程序
3.3.1 网格离散程序
3.3.2 主程序
3.3.3 单元Bei子块计算程序
3.3.4 单元Cei子块计算程序
3.3.5 单元Deij子块计算程序
3.3.6 单元Fei子块计算程序
3.3.7 网格细化程序
3.3.8 压力插值程序
3.3.9 矩形网格绘制程序
3.3.1 0计算结果
3.4 罚函数有限元求解
3.4.1 计算区域的离散
3.4.2 插值函数及其相关计算
3.4.3 加权余量方程
3.4.4 单元方程的建立
3.4.5 总体方程的组合
3.4.6 压力的计算
3.4.7 求解流程
3.5 罚函数有限元程序
3.5.1 主程序
3.5.2 单元DPeij子块计算程序
3.5.3 单元内结点压力计算程序
3.5.4 其他程序
3.5.5 计算结果
第4章 非牛顿流体流动的有限元求解
4.1 计算实例及数学方程
4.1.1 计算实例
4.1.2 数学方程
4.2 有限元方法求解方程
4.2.1 计算区域的离散
4.2.2 插值函数及其相关计算
4.2.3 加权余量方程
4.2.4 单元方程的建立
4.2.5 总体方程的组合
4.2.6 迭代求解流程
4.3 程序编写
4.3.1 网格生成程序
4.3.2 主程序
4.3.3 单元结点黏度计算程序
4.3.4 单元Deij子块计算程序
4.3.5 其他程序
4.4 计算结果分析
4.4.1 网格数量对计算精度的影响
4.4.2 求解问题的速度、压力、剪切速率和黏度分布
4.4.3 物性参数对出口速度分布的影响
4.4.4 入口压力对出口流量的影响
第5章 考虑惯性项影响的牛顿流体流动有限元求解
5.1 求解实例和数学方程
5.1.1 求解实例
5.1.2 数学方程
5.1.3 边界条件
5.2 有限元求解
5.2.1 计算区域的离散
5.2.2 插值函数及其相关计算
5.2.3 加权余量方程
5.2.4 单元方程的建立
5.2.5 总体方程的组合
5.2.6 非线性方程组的求解方法
5.3 相关程序编写
5.3.1 “速度项提出法”+“Newton Raphson迭代法”相关程序
5.3.2 “直接推导法”+“线性化交替迭代法”相关程序
5.4 结果分析
5.4.1 两组程序计算结果对比
5.4.2 惯性项影响分析
第6章 非牛顿流体非定常流动的有限元求解
6.1 求解实例和数学方程
6.1.1 求解实例
6.1.2 数学方程
6.1.3 边界条件
6.2 有限元求解
6.2.1 计算区域的离散
6.2.2 插值函数及其相关计算
6.2.3 加权余量方程
6.2.4 单元方程的建立
6.2.5 总体方程的组合
6.2.6 非定常问题非线性方程组的求解方法
6.3 相关程序编写
6.3.1 网格生成程序
6.3.2 主程序
6.3.3 单元方程子块计算程序
6.3.4 Bird?Carreau本构模型的单元内结点黏度计算程序
6.3.5 单元内速度积分程序及单元面积计算程序
6.3.6 其他程序
6.4 结果分析
第7章 与时间有关的热传导问题的有限元求解
7.1 求解实例和数学方程
7.1.1 求解实例
7.1.2 数学方程和边界条件
7.2 热传导方程的有限元求解
7.2.1 计算区域的离散
7.2.2 插值函数及其相关计算
7.2.3 加权余量方程
7.2.4 单元方程的建立
7.2.5 总体方程的组合
7.2.6 代入边界条件及迭代求解
7.3 相关程序编写
7.3.1 网格生成程序
7.3.2 主程序
7.3.3 单元温度积分计算程序
7.3.4 单元面积计算程序
7.3.5 热传导项CDe子块计算程序
7.3.6 时间项CDe子块计算程序
7.3.7 热传导边界项CDBe子块计算程序
7.3.8 其他程序
7.4 计算结果分析
7.4.1 区域温度变化
7.4.2 加热热流密度对升温过程的影响
7.4.3 空气温度对升温过程的影响
第8章 速度与温度耦合问题的有限元求解
8.1 求解实例和数学方程
8.1.1 求解实例
8.1.2 数学方程
8.1.3 边界条件
8.1.4 与剪切速率和温度有关的本构方程
8.2 能量方程的有限元求解
8.2.1 计算区域的离散
8.2.2 插值函数及其相关计算
8.2.3 加权余量方程
8.2.4 单元方程的建立
8.2.5 总体方程的组合
8.2.6 能量方程与N?S方程组耦合时的求解流程
8.3 相关程序
8.3.1 网格生成程序
8.3.2 主程序
8.3.3 单元结点黏度计算程序
8.3.4 单元Deij子块计算程序
8.3.5 单元Cei子块计算程序
8.3.6 单元Bei子块计算程序
8.3.7 单元Fei子块计算程序
8.3.8 单元CDe子块计算程序
8.3.9 单元DLe子块计算程序
8.3.10 单元NHe子块计算程序
8.3.11 单元CDBe子块计算程序
8.4 计算结果分析
8.4.1 计算结果
8.4.2 入口压力对温度分布的影响
8.4.3 流动区域收敛比对温度分布的影响
8.4.4 黏性耗散对温度分布的影响
参考文献