论文摘要
非结构网格有限体积法适合应用于复杂几何结构的流场计算,而SIMPLE系列算法适合应用于工业上广泛存在的低速流、不可压流计算,欧拉-拉格朗日方法适合应用于含稀疏颗粒相的流场计算。本文正是以这些理论为基础展开工作,包括CGNS网格前后处理,压力修正算法的实现和拓展,化学组分源项模型,颗粒相与气相之间的相互作用等。分述如下:本文提出了两种将CGNS非结构网格点数据结构转化为非结构网格面数据结构的快速算法,以作为本文操作非结构网格的基础:内部网格面的整数hash表算法和内部双面壁边界网格面的实数hash表算法。由于CGNS非结构网格存储的是网格结点数据,例如网格结点坐标、每一个网格由哪些结点构成等数据信息;而有限体积法离散过程中,需要用到大量的面积分,所以更多的需要网格面数据结构。内部网格面的构造根据网格面两侧的网格在这个公共网格面上的结点序号相同的特点设计算法,双面壁网格面则根据重合在一起的两个边界网格面中心距离接近于零,或者说小于某一个特征小量,比如该网格面面积的算术平方根来设计算法。两种hash表算法具有快速,内存消耗量小于后续有限体积求解器内存峰值的特点,使前处理算法不会成为整个计算的瓶颈。在压力修正方程的构造过程中,考虑温度修正量对密度修正量的影响,通过静焓方程构造温度修正量与压力修正量的关系,而这种做法的实际效果相当于增强了压力修正量方程线性代数方程组的主对角线。对于可压缩流计算,在流量计算时,在线性插值和迎风格式的基础上,网格面上的密度值采用低阶和高阶的混合格式,动量方程中的压力项处理为滞后的压力梯度项,即线性插值计算方法,和网格面上的压力积分之和,网格面上的压力值则可以采用低阶和高阶混合格式。这些做法兼顾了分辨率和计算的稳定性。在喷雾燃烧数值模拟中考虑了湍流反应组分源项模型如EBU模型、ED模型、Arrinius/EBU/ED模型和Presumed PDF模型。对于颗粒相和气相之间的相互作用,本文一是将颗粒相对气相的反作用阻力中的部分阻力系数剥离,由拉格朗日系统转入欧拉系统,加入到气相动量方程线性代数方程组的主对角线上,这一系数为正实数,增强了主对角线的系数,也增强了含颗粒流的气相迭代计算稳定性;二是在采用时间推进方法计算稳态流场过程中,对耦合的颗粒相与连续相之间计算采用显式和隐式的计算方法进行了研究,这样做的目的在于相对传统的显式时间计算方式,可以采用一个较大的时间步长值进行迭代,从而减少时间推进时的计算量,使得喷雾燃烧数值模拟的计算稳定性得以增强同时减少计算量。在课题组研制的自主知识产权流体力学软件开发过程中,本文作者承担了部分功能模块的研发工作,取得了上述研究成果。并就一些算例与商用软件和经典计算结果予以对比,验证了本文算法的正确性。
论文目录
相关论文文献
- [1].三维空间非结构网格生成方法研究[J]. 计算机工程与应用 2017(22)
- [2].基于非结构网格间断有限元方法在可压缩流体中的应用[J]. 工程热物理学报 2012(07)
- [3].并行燃烧数值模拟计算优化——面向自适应非结构网格的动态负载平衡方法[J]. 计算机工程与应用 2013(21)
- [4].基于神威·太湖之光的非结构网格众核优化技术[J]. 计算机工程 2019(06)
- [5].基于非结构网格的三维水流数值模拟研究[J]. 西北农林科技大学学报(自然科学版) 2017(11)
- [6].非结构网格布点方法研究进展[J]. 计算机工程与应用 2008(32)
- [7].非结构网格催化壁面的隐式处理方法[J]. 推进技术 2018(02)
- [8].非结构网格二阶有限体积法中黏性通量离散格式精度分析与改进[J]. 力学学报 2018(03)
- [9].基于非结构网格的气动弹性数值方法研究[J]. 航空动力学报 2009(09)
- [10].非结构网格边粘性公式[J]. 计算物理 2008(01)
- [11].非结构网格下稀疏下三角方程求解器众核优化技术研究[J]. 计算机科学 2019(S1)
- [12].面向非结构网格应用并行程序的编程工具[J]. 计算机科学与探索 2018(02)
- [13].基于非结构网格的高效求解方法研究[J]. 计算力学学报 2012(02)
- [14].稀疏非结构网格上的亚声速流高精度数值模拟[J]. 航空学报 2009(02)
- [15].二维自适应非结构网格DSMC并行算法研究[J]. 计算力学学报 2009(02)
- [16].基于二维非结构网格的全隐式,迎风格式欧拉方程研究[J]. 科学技术与工程 2008(03)
- [17].一种TVD方法在三维非结构网格中的应用[J]. 航空计算技术 2008(05)
- [18].基于伴随方程方法的非结构网格自适应技术及应用[J]. 空气动力学学报 2011(03)
- [19].基于非结构网格的溃坝水流干湿变化过程数值模拟[J]. 水力发电学报 2008(05)
- [20].混合非结构网格下壁面最短距离的快速计算方法[J]. 西北工业大学学报 2014(04)
- [21].基于非结构网格流场超大规模并行计算[J]. 空气动力学学报 2019(01)
- [22].基于非结构网格的三维中子输运方程求解方法[J]. 原子能科学技术 2014(10)
- [23].二维潮流计算中非结构有限体积法研究和应用[J]. 海洋工程 2009(04)
- [24].基于非结构网格的近海海域污染物迁移三维数值模拟并行算法[J]. 计算机应用 2019(S1)
- [25].基于非结构网格隐式算法的GPU加速研究[J]. 计算机系统应用 2018(05)
- [26].基于非结构网格的风压分配方法[J]. 石家庄铁道大学学报(自然科学版) 2014(04)
- [27].非结构网格上p型多重网格法流场数值模拟[J]. 计算力学学报 2011(03)
- [28].基于非结构网格求解器的CHN-T1标模气动特性计算研究[J]. 空气动力学学报 2019(02)
- [29].自然对流换热的高精度非结构网格有限体积法[J]. 计算物理 2013(03)
- [30].Lagrange非结构网格高阶交错型守恒气体动力学格式[J]. 应用数学和力学 2014(01)