渗流分析的非正交网格有限体积法

将有限体积法的体积积分用Ｇｒｅｅｎ公式转化为单元边界积分，使之能应用于非正交网格的渗流分析计算；保留有限体积法计算简单实用、计算精度高、省机时省内存的特点，又有网格划分灵活方便，对复杂计算域适应性强的优点。     

NS方程的非结构化网格方法及其差分格式

采用同位非结构化网格上的有限体积法,对ＮａｖｉｅｒＳｔｏｋｅｓ 方程的ＳＩＭＰＬＥ 算法及差分格式进行了研究．提出了适用于非结构化网格的不用求解单元顶点变量值的二阶混合差分格式．该格式的优点在于：１） 减少了计算工作量;２） 避免了普通混合差分格式因使用简单的一阶迎风差分格式所引起的网格界面方向相关性的问题．最后,采用三角形网格,利用提出的方法及差分格式,对方腔内的驱动层流及绕翼型湍流进行了数值计算,计算结果与基准解或实验值的符合程度优于普通混合差分格式     

一种基于非结构化网格的流体体积法界面重构方案

基于三角形非结构化网格,发展了一种类似于结构化网格简单直线界面计算方法的界面重构方案.该方法为纯欧拉方法,根据位于两相界面网格的相邻网格流体信息在该网格内部构造平行于某边的界面线,通过界面线与网格边界围成的几何形状计算流体输运,采用修正系数保证全局的质量守恒.对均匀平直速度场中的同心双运动界面、旋转速度场中的马蹄形界面和剪切流场中的圆形运动界面模型等3个经典算例,使用速度反向还原法检测了所提方法,7600个非结构化网格上的界面构造效果达到了6400个结构化网格的精度.     

有限单元法的非结构化网格及其自动生成

本文介绍了有限单元法非结构化网格的基本原理及其自动形成方法。由于有限元法是一种离散的数值求解方法，其近似求解方法的精确度，很大程度上取决于所形成网格的质量。另外，对于工程中一些形状复杂的问题，一般的网格生成方法很难对其进行离散。非结构化网格及其自动生成，使复杂形状的工程问题容易地进行离散，改善所形成网格的质量，提高近似计算的精确度，并且在有限元网格修正自适应分析中具有重要作用。     

非规则网格有限体积法处理不规则边界的研究

提出用非规则网格有限体积法计算不规则边界区域的传热和流体流动问题,分析并解决了采用此方法计算时所遇到的困难;混合差分方案的有效实施,系数矩阵对角占优特性的保持以及压力修正方程的合理构造等,最后,通过算例证明了该方法能够准确的计算不规则边界区域的传热和流体流动问题。     

四边形网格的离散有限体积方法

讨论了椭圆方程在四边形网格上的离散有限体积方法．给出了离散解与精确解的H^1误差估计．最后我们给出了数值算例．     

基于非结构化网格的2.5维时间域激发极化反演

针对我国西部山区起伏地形,提出通过Delaunay三角剖分生成非结构化三角网格,可以根据实际地形情况,较好地控制节点分布,从而降低初始模型的几何离散误差。在反演中对模型参数施加最小光滑和背景约束信息,从而提高反演的分辨率。通过在正演求解稀疏矩阵时采用稳定型双共轭梯度法和计算系数矩阵时利用电位互换等手段,可有效地提高反演计算速度。最后,为验证反演方法的有效性,编制了跨平台反演软件;并对2例地电模型和实测数据进行反演试算。计算结果表明:反演耗费时间少,仅迭代3次拟合差便趋于稳定,反演结果能较好地刻画异常体形态。     

三维Euler方程的非结构网格非振荡有限体积法

考虑三维Euler方程,对非结构四面体网格给出了一类时空分离的非震荡有限体积方法。空间方面采用Roe通量进行离散,同时基于物理量守恒运用最小二乘和极值原理在每个四面体单元构造非振荡线性重构多项式,时间离散采用两步Runge-Kutta TVD方法。文末对三维球形爆炸问题进行数值模拟,得到了满意的结果。     

基于气泡堆积的非结构化网格生成技术

为了有效地提高网格的质量,基于物理思想发展了一种基于气泡堆积的非结构化网格生成算法.该算法不同于其他典型的非结构化网格生成算法,具体步骤为:首先按一定的规则向给定区域内添加具有虚拟质量的气泡,使其紧密堆积以充满整个计算区域;然后引入气泡之间的相互作用力得到各气泡的运动控制方程,通过求解该方程来调整各气泡的位置,反复迭代使所有气泡受力最小,最终达到整个系统的受力平衡,从而有效避免相邻气泡过于重叠或分离,实现调整节点位置优化网格的目的;最后通过Delaunay三角形化方法有效地联接这些气泡的中心,生成非结构化三角形网格.基于上述算法开发了相应的程序,将网格生成结果与Delaunay划分并结合拉普拉斯光顺化方法生成的网格进行了对比,发现网格的几何不规则度大大减小,质量得到了明显的提高.     

基于弹簧近似的非结构化网格自适应处理方法

首先对弹簧近似法进行了说明,并对节点弹簧法和边弹簧法进行了对比.在边弹簧法的基础上,从调节弹簧平衡长度的角度出发,提出了一种非结构化网格的自适应处理方法,并结合有限单元法,给出了网格自适应处理的流程.最后,对一维对流扩散方程、二维线性抛物型方程以及欧拉方程进行了求解.计算结果表明:文中提出的网格自适应处理方法可以使几个问题的网格节点布置更为合理,有效地降低数值误差,提高计算精度,特别对于算例中所涉及的二维非结构化网格,自适应化过程中并没有出现网格重叠现象.     

二维非结构网格的一个TVD型有限体积方法

文章将有关文献提出的有限差分算法应用到有限体积方法中,在无结构任意三角形网格上建立了一种守恒型有限体积的高分辨率算法,该方法简单易行,应用改进后的算法到二维浅水波方程的求解上,进行了二维部分溃坝、圆坝两个算例的数值模拟,模拟结果表明:改进后的算法既有较高的激波捕捉能力,又能处理界面任意复杂的计算问题。     

无结构三角形网格下运动界面追踪的WENO有限体积法

在处理运动界面追踪问题的流体体积函数(VOF)法的基础上,给出了一种无结构三角形网格下的高分辨率的运动界面捕捉方法.该方法采用高精度的加权本质无振荡(WENO)有限体积格式离散VOF函数的空间导数,采用三阶TVD Runge-kutta方法离散时间导数,采用Lax-Friedrichs通量作为数值流通量.数值试验结果表明,用该方法来进行旋转速度场和剪切速度场的运动界面追踪,可以得到与理论解非常一致的追踪结果.     

基于四边形剖分的多尺度有限体积方法

在一类多尺度方法（heterogeneous multi-scale method,HMM）的基础上,结合有限体积方法（finitevolume method,FVM）,提出多尺度有限体积方法（HMM-FVM）,用于数值计算一类多尺度抛物型偏微分方程,同时给出相关误差估计.计算结果表明,多尺度有限体积方法比传统的有限元法、有限体积法有效,既节省计算量,又有较高的精度.     

基于无结构网格的有限体积法和颗粒方法相结合求解湍流的联合PDF方程

本研究着力于发展一种新的有限体积法和颗粒方法的混合算法 ,用于湍流 (冷态流以及反应流 )的速度 频率 标量联合PDF(probabilitydensityfunction)输运模型方程的数值求解 .采用一种基于无结构网格的有限体积法 ,求解质量、动量、能量平均守恒方程和平均状态方程 ,为联合PDF的输运模型方程的求解提供其所需的平均速度场和平均密度场 ;颗粒方法中 ,用MonteCarlo算法来求解PDF方程 ,为平均场的计算提供能源源项、雷诺应力、热通量、比热比等信息 .两种方法之间的耦合是一个关键 ,为了验证该算法的可行性 ,论文还对二维V型槽后的湍流冷态流场进行了数值模拟 ,并与实验、大涡模拟、k 模型的结果进行了比较 ,结果表明 ,该算法的结果与实验结果吻合得较好 .     

弹性力学中三种有限体积法方案的对比

为了研究弹性力学中三种典型的有限体积法方案(FVM1,FVM2和FVM3)的精度差异问题,采用理论推导和算例验证的方法,建立了基于有限体积法原理的力平衡积分方程,利用三角形网格进行离散化,在计算出内部单元和边界单元对力平衡贡献的基础上,对每一种有限体积法方案确立了与有限单元法方案(FEM)刚度方程组相类似的线性方程组,对比所建立的各个方程组的待定系数.结果表明:方案FVM2与FEM精度相同,即为最优方案;在不考虑集中力的情况下,方案FVM3的精度也与FEM相同.最后以一个数值算例,对这三种有限体积法方案进行了对比分析,验证了上述分析.该结果对工程计算数值模拟具有一定的参考价值和指导意义.     

基于八叉树网格的直流电阻率法有限元数值模拟

在VC环境下,使用C语言,实现了基于八叉树网格的三维直流电阻率法有限元数值模拟程序。三个典型模型的算例被用来证明八叉树网格方案适合于复杂地电模型的数值模拟：数值模拟了球体模型,随着网格对球面边界的不断逼近,数值模拟的精度也逐步提高;计算了覆盖层下倾斜板状高阻体模型,在偶极-偶极视电阻率断面上,能清晰地看见高阻体模型的形态、埋深;仿真了联合剖面法测量三维地形下埋藏低阻椭球体的情况,由于地形的影响,在视电阻率曲线上不能发现低阻异常的存在,利用比较法进行地形改正后,低阻异常得到凸显。