首页> 中文学位 >基于FEniCS的二维非结构网格Navier-Stokes方程求解器研究
【6h】

基于FEniCS的二维非结构网格Navier-Stokes方程求解器研究

代理获取

目录

声明

第1章 绪论

1.1 研究背景与意义

1.2 国内外研究进展

1.2.1 CFD求解器研究进展

1.2.2 CFD有限元方法的研究进展

1.2.3 FEniCS简介

1.2.4 稳定化方法研究进展

1.2.5 自适应有限元研究进展

1.3 开发环境

1.4 论文的主要创新点

1.5 论文的结构组织

第2章 数值计算方法与验证

2.1 有限单元法

2.1.1 弱形式

2.1.2 有限单元形式

2.2 速度-压力耦合问题中的分裂格式

2.2.1 Chorin投影法

2.2.2 增量压力修正法

2.3 Stokes问题与混合有限元法

2.3.1 Stokes问题数值稳定条件

2.4 求解器数值验证

2.4.1 案例1:方腔流动

2.4.2 案例2:圆柱绕流

2.5 小结

第3章 流体求解中的稳定化方法

3.1 抑制对流占优的稳定化方法

3.1.1 标准Galerkin方法

3.1.2 人工粘性方法

3.1.3 SUPG方法

3.2 对流方程数值计算中的稳定化方法

3.2.1 问题描述

3.2.2 连续初始条件数值模拟

3.2.3 不连续初始条件数值模拟

3.3 不可压缩流体模拟中的稳定性方法

3.3.1 圆柱绕流数值模拟中的稳定性方法

3.3.2 亚临界区域绕流计算中的稳定性方法

3.4 小结

第4章 自适应有限元方法

4.1 自适应有限元数值方法

4.1.1 后验误差估计

4.1.2 一致分裂格式

4.1.3 伴随方程的求解

4.1.4 自适应方法技术路线

4.2 圆柱绕流自适应计算

4.3 方型体绕流自适应计算

4.4 小结

第5章 结论与展望

5.1 工作总结

5.2 展望

致谢

参考文献

攻读学位期间发表的学术论文

展开▼

摘要

在随着现代计算集群的计算速度和储存容量大幅度提高,借助于数值水池为基础的分析评估、预报及优化设计的虚拟综合集成现代计算流体力学体系,为复杂海洋环境下船舶航行性能与海洋平台结构安全性的研究提供了所必需的虚拟试验环境和技术手段。复杂几何外形下,非结构化网格能较快速的完成前处理工作,并能达到一定的数值精度,是现代CFD技术发展应用的主流方向。有限元法在复杂边界条件,非结构化网格上适应性良好,是求解Navier-Stokes方程的主流数值方法之一,其具有良好的科研价值和工程应用背景。 本文使用Python语言,基于FEniCS等科学计算库下开发了算子分裂格式,并对比了不同稳定化格式在对流占优问题的计算精度、效率与稳定性;基于伴随理论,开发基于误差指示器的非结构网格自动加密模块,发展了一种适用于有限元法的非结构网格自适应技术。对目标函数较为敏感的网格区域加密并优化网格,简化了非结构化网格前处理流程,提高了绕流问题的计算精度。提出了将一致分裂格式应用到自适应有限元中,简化了误差项,相较于经典的投影法和增量压力修正格式提高了网格自适应中的计算精度与效率。 本文的主要内容如下: 1.给出了有限元在泛函空间上的数学定义与符号规定,以及Navier-Stokes方程的弱形式推导。选择合适的算子分裂格式与压力场与速度场的分裂形式,对压力场-速度场选择合适的混合单元,验证了有限元方法在方腔流,圆柱绕流算例上的数值精度; 2.讨论了混合有限单元中的对流占优问题,以及两种常用的抑制对流占优的稳定性方法,分别为人工粘性法与SUPG方法。在对流方程,圆柱绕流问题中将标准Galerkin方法与上述稳定性方法在计算精度与计算效率进行了对比与讨论。最后,讨论了以上三种方法在临界区域内流场中的计算情况; 3.基于一致分裂格式,针对物体边界网格密度不足的问题,开发自适应模块,应用在有限元流场求解程序上.讨论了网格加密比率、最大加密迭代次数以及初始网格密度在圆柱绕流问题上的误差分析.依照不同的受力构建了误差指示器,本文发现圆柱体/方型体绕流问题在不同的误差指标上的计算误差存在显著差异.最后讨论了方形体绕流计算中不同分裂格式的计算效率与计算精度情况。 本文针对使用有限元法求解不可压NS方程问题所存在的几大困难,在前人现有工作的基础上,编程实现了高效稳定的稳定化方法与网格自动加密算法,可用于解决船舶与海洋工程中常见的流体问题,为结构物设计优化提供帮助。

著录项

相似文献

  • 中文文献
  • 外文文献
  • 专利
代理获取

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号