非结构任意多边形网格辐射扩散方程有限体积格式

本文基于非结构任意多边形网格体系,给出了求解辐射扩散方程的中心型有限体积格式,格式中出现的网格节点未知量由相邻的网格中心未知量加权给出,综合考虑网格几何及扩散系数的影响,给出了节点未知量的一种加权方式,数值实验表明格式在各种非结构网格上具有较强的适应性.

任意多边形的非结构网格快速生成算法研究

对于包含大尺度运动边界的CFD数值模拟,网格重构是其中的关键,快速稳定的网格生成技术是其中的重要组成部分。建立了基于有向边的适用于任意多边形的快速三角初始化算法;证明了最长边剖分网格细化算法在一定条件下发散,并结合Delaunay边交换技术使细化算法封闭;建立了基于顶点弹簧理论的网格优化方法,以提高网格生成的质量。结果表明,算法具有较好的鲁棒性和高效性。

气泡堆积法生成曲线多边形区域非结构化网格及其应用

对气泡堆积法进行改进,发展了一种基于气泡堆积法生成曲线多边形区域非结构化网格的算法,将曲线边界映射为直线完成气泡添加和位置动态调整,再通过弧长参数化的方法将气泡位置映射回曲线边界,避免了移动调整过程中气泡偏离曲线的问题.提出了一种简便的判断点在复杂曲线多边形区域内的方法,简单且易于编程.通过添加不同大小的顶点气泡以及在内部设置人工点源,利用加权平均法实现了网格的局部加密.应用基于非结构化同位网格的SIMPLE算法对环扇形空腔顶盖驱动流进行了数值模拟,不同雷诺数下的计算结果与文献结果吻合较好.

任意多面体的快速非结构网格生成算法

对于包含复杂气动外形的CFD数值模拟,网格生成是其中的关键,快速稳定的网格生成技术是其中的重要组成部分。建立了基于可视面的适用于任意多面体的快速初始化算法;改进了Delaunay生成算法后使用Delaunay改进生成算法细化网格,讨论了网格质量判定依据对网格生成的影响,通过开发Delaunay面交换技术优化网格生成过程;建立了基于顶点弹簧理论的网格优化方法,以提高网格生成的质量。结果表明,本文建立的算法效率较高。

非结构多边形网格解扭和重分方法

在多维流体动力学计算中，流体运动和计算网格的关系可以分为两种情况。一是Lagrangian方法，即网格跟随流体运动；二是Eulerian方法，即流体流过固定；下动的网格。一般计算网格的运动是任意的。这就对应于任意Lagrangian—Eulerian（ALE）方法。ALE方法的核心是通过调整网格运动，使得数值模拟的精度、效率有所提高。它的主要步骤是：显式Lagrangian步；网格重分，即得到新的计算网格；物理量重映，即将Lagrangian步的计算结果变换到新网格上。在这3步中，较少研究网格重分。数值模拟和网格重分的一个基本前提是网格是合理的，或者说网格不能发生翻转，网格应当是凸的。而Lagrangian步数值模拟会造成网格扭曲，因此在网格重分前进行网格解扭是十分必要的。文中描述了通用的网格解扭、重分算法，使得解扭、重分后的网格有较好的几何品质，同时尽可能接近Lagrangian网格。

求解N-S方程的多边形网格与三角形网格对比

本文介绍了基于边的数据结果离散对流扩散方程的方法,编制了适用于任意多边形单元网格的有限容积法计算程序,利用方腔自然对流、外掠后台阶和移动顶盖驱等算例考察了三角形网格与多边形网格在数值计算方面的性能差别。结果表明,采用三角形网格收敛较慢且不容易获得网格无关的解,而采用与三角形网格对偶的多边形网格则收敛较快且易获得网格无关的解。

任意各向异性介质三维非结构谱元法直流电阻率正演模拟研究

各向异性是地电异常解释中不可忽视的因素,广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究,结合非结构谱元法建立模拟算法,充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力,提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值,实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度,进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征,并通过计算针对不同各向异性参数的视电阻率极性图,探究地下介质各向异性特征识别方法.最后,我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用.

多介质拉氏自适应流体动力学软件LAD2D研制及其应用

拉氏方法是内爆动力学过程数值模拟的主要方法.针对高温、高压、多介质和大变形等内爆问题,采用非结构任意多边形网格底层管理、计算过程中网格邻域可变技术,以及拉氏自适应网格加密方法和层次化、模块化程序设计思想,自主研发非结构拉氏自适应网格流体动力学软件LAD2D.从物理模型、计算方法、程序设计、程序验证与确认、大变形问题数值模拟等方面系统地介绍LAD2D.LAD2D对多介质爆轰弹塑性流体大变形问题有很强的适应能力.

数值模拟法分析压裂生产井的裂缝开度

非常规油气藏水力压裂后的数值模拟中,裂缝开度普遍采用的是均匀分布的方式。但是实验室的试验结果表明,有效应力发生变化会导致裂缝开度分布并不均匀。而在绝大多数数值模型中往往不考虑这一现象。本文调研了生产过程中压力对裂缝开度的影响作用。在数值模拟模型中,裂缝开度的作用比较明确,与此同时,考虑了闭合应力的影响。通过一个非结构化的泰森多边形网格完成离散型裂缝建模,该模型考虑了复杂的裂缝几何形状。