QUICK格式在非结构化网格中的应用

在有限容积法的基础上发展了非结构化网格的对流项二次迎风插值（QUICK）格式．详细推导了扩散项采用格林函数法，对流项采用改进的QUICK格式离散方程，对项盖驱动流和圆柱绕流问题进行了计算，讨论了不同Re下计算的准确性和格式的收敛性，并与高精度结构化网格计算结果进行对比分析．结果表明，该格式的临界网格Peclet数为8／3左右，与中心差分相比较，该格式计算精度相当，对流稳定性好，收敛速度高．同等条件下较结构化网格对复杂区域的模拟更接近实际测量结果，是一种对复杂区域计算有应用前景的对流格式．     

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

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

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

结构化与非结构化网格中的导热计算

用基元有限容积方法求解热传导问题,同时采用结构化网格和非结构化网格,并和精确解相比较.虽然从两种网格形式,都能得到满意的结果,但结构化网格较三角形非结构化网格在规则区域更精确一些.为了提高非结构化网格的精度,二次扩散项是十分重要的.     

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

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

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

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

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

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

二维各向异性非结构化网格的自动生成与应用

根据生成适体坐标网络空间变换的思想，结合Delaunay非结构化网络生成技术，实现了各向异性非结构化网格的生成。部分算例表明，通过给定测度矩阵，可以控制生成预期的网格。 针对各向异性导热问题的控制容积有限元法（CVFE）离散数值计算表明，用各向异性网格比各向同性网格收敛快。     

非结构化网格中圆管流动的传热计算

圆管道内充分发展的层流和传热是热传导方程应用中一个有趣而重要的分支．本文给出几种在不同流道几何形状和热边界条件下,用非结构化网格方法进行数值计算的结果,并与结构化网格的计算结果进行比较分析．从而证明本方法和所给程序的正确性和实用性．     

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

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

通风分裂算法在有限元非结构网格中的推广

将在结构网格中发展起来的通风分裂算法推广到有限体积非结构网格中,并成功地进行了定常跨声速二维欧拉方程的数值计算．此格式具有流通量矢量分裂格式的简单性及通量差分分裂格式的高阶精度．最后给出的ＮＡＣＡ００１２翼型绕流的结果与实验结果符合的较好．     

采用非结构网格的IDEAL算法实现及定压边界条件处理

在非结构网格系统中采用有限容积法对不可压N-S方程进行了离散,其中对流项采用二阶迎风格式,扩散项采用超松弛校正格式.通过修正的动量插值格式构造了压力方程,利用IDEAL算法处理压力与速度的耦合.利用Taylor展式导出了邻近边界内节点的压力计算公式,进而将边界静压引入到压力方程的迭代计算中,确保了边界信息向内部区域的有效传递.根据质量守恒定律分别导出了压力入口和出口边界的速度计算公式,结合伯努利定律给出了总压边界的处理方法.通过压力驱动管内层流、三通管层流以及浮力驱动流等几个典型算例对文中方法进行了考核,并将计算结果与文献中的基准解和实验值进行了对比.结果表明,不管是静压边界还是总压边界,采用该处理方法均能获得令人满意的结果,尤其对于大空间自然对流问题,可以在保证精度的前提下节省大量的存储空间和运算时间.     

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

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

污染物大气扩散数值模拟中伪扩散的控制

用数值模拟的方法预测污染物在大气中的扩散已经被广泛采用,但数值模拟产生的伪扩散一直是影响模拟结果精度的重要因素.为此分析了数值模拟中伪扩散产生的原因,用一阶迎风差分和对流项二次迎风插值格式对一稳态污染物扩散模型进行了计算,并比较了计算结果,得出了有实际意义的结论.     

控制体有限元方法中质量加权迎风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散．为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式．通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法对回流问题数值模拟的假扩散．     

利用高精度迎风紧致群速度控制法模拟激波与柱形界面相互作用

 利用带有群速度控制的迎风紧致法和Level Set 法相结合,对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟,并与实验结果相比较。从模拟结果看,数值模拟的结果与实验结果基本相同,正确地模拟出界面的位置、激波的强度和速度,这表明文中给出的方法是合理可行的。     

非结构网格海洋模式在洞头人工渔礁区选址中的应用

应用一种非结构网格海洋模式(FVCOM)建立了洞头列岛及其临近海域的水动力数值模式,水平方向三角网格的采用使模型对大陆海洋线以及岛礁边界的拟合达到了很高精度.经流速、流向和潮位验证后,精度良好的模型被用于人工渔礁区选址的研究.以大潮底层最大流速和大潮底层平均流速为主要判别指标,结合海域地形和水深分布,选定了底层最大流速小于0.6 m/s、平均流速小于0.4 m/s的海域为适宜人工渔礁建设的礁区范围.     

隐式有限体积TVD格式在羽流场中的应用

火箭羽流会造成噪声、烟雾、热辐射、环境污染和信号衰减等多种效应,因此研究排气羽流问题至关重要。采用三维隐式有限体积TVD格式,通过求解定常可压缩Euler方程,数值模拟了来流马赫数为1．0、1．5、2．0和喷管出口压力比为1．0、2．0、3．0时的火箭羽流场,并对各流场参数和噪声进行了计算和比较。所采用的二阶格式具有分辨率高、收敛快速稳定、可较好模拟复杂计算区域和流场波系等优点。声压级在相交激波点附近出现剧烈变化,喷管出口处的激波和排气羽流造成强噪声。     

HLLC方法在二维非结构动网格上的应用

以自适应非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的Godunov一阶精度方法求解Euler方程，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，空间和时间都是一阶精度，对含有运动边界的多体干扰流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步进行多体分离数值模拟打下了基础。     

双线性等参元对称有限体格式在二维三温辐射热传导问题的应用

首先给出了二维三温辐射热传导问题的数学模型和有关数据,在四边形网格下,利用牛顿方法线性化,双线性等参元对称有限体格式离散,从而获得了与实际物理现象相符的数值结果,也验证了该格式的有效性.