基于同位网格的一种新的定压边界处理方法

将SIMPLE-P算法中处理压力边界的思想推广到同位网格中,综合考虑压力和速度边界条件的数学特性,重新推导边界控制体的压力方程,借鉴IDEAL算法处理压力和速度耦合的思想,讨论定压边界流动传热问题的求解步骤。通过Kovasznay流、压力驱圆管层流和浮力驱动流等几个典型算例对新方法和SIMPLE-P算法进行考核并将计算结果与文献中的数值解和实验值进行对比。结果表明:新方法是准确有效的,对于定压边界问题的求解较之SIM-PLE-P算法具有更好的收敛性能;对于大空间自然对流问题使用新方法在较小的计算区域就可获得高精度的结果,从而节省了大量的存储空间和运算时间。     

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

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

非正交同位网格中的SIMPLE算法

利用有限体积法建立了非正交同位网格系统中的压力校正方程,详细给出了压力校正方程的离散形式,并对正方腔及倾斜角为45°的斜方腔内流场进行了数值模拟计算.结果表明,非正交同位网格系统中的SIMPLE算法具有编程简单,计算量、存储量小,计算效率高,收敛速度快的优点,可用于几何形状比较复杂流场的计算.     

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

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

非结构网格上浅水方程的LU-SGS隐式算法

针对浅水数值模拟的特点,建立了可以在非结构网格上求解浅水方程的LU-SGS隐式算法.对复杂地形的底坡项进行了修正,并利用淹没节点法处理移动边界问题.利用算例对此隐式算法进行了模拟验证,并与Roe显式算法进行了比较.结果表明,隐式算法可以有效地缩短计算时间,具有良好的计算和收敛效果.长江口实际潮流模拟结果也证明了此隐式算法处理实际浅水问题的能力.     

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

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

基于非结构网格的三维浅水模型

采用水位积分平衡法并引入边壁滑移系数,解决σ坐标系下有限体积离散时产生的静水压力项与底坡项的平衡问题和边壁阻力的模拟问题,建立非结构网格上三维浅水模型.风生流试验计算结果表明,模拟的纵向流速的垂向分布与理论分布吻合较好;Holtz丁坝试验计算结果表明,引入边壁滑移系数能够模拟坝体附近上游小回流区和下游小回流区;淹没丁坝...     

对流项二次迎风插值格式在非结构化网格中的应用

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

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

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

对流换热的高精度非结构网格算法

以层流流场的非结构网格高精度算法为基础,加入湍流模型和能量方程,实现了适合复杂区域中层流和湍流对流与换热问题的非结构网格高精度计算。运用发展的算法计算自然对流和后台阶突扩分离流,同时采用FLUENT软件在网格和参数设定相同的条件下对比计算,运用实验结果或经典数值结果比对考核。和FLUENT的计算结果相比,本文算法中高精度格式的计算结果与实验结果或经典数值结果吻合程度整体上更高,验证了本文算法的正确有效性和高精度。     

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

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

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

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

SIMPLE算法中压力修正方程边界条件确定新探

不可压缩粘性流体Ｎａｖｅｒ－Ｓｔｏｋｅｓ方程ＳＩＭＰＬＥ算法求解中压力修正方程起着很关键的作用。但其边界处理的特殊性常给正确应用带来的一定困难。提出从连续性方程本质出发确定压力修正方程边界条件的新方法。     

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

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

非结构网格在周期性充分发展流动数值模拟中的应用

为了解决周期性充分发展流动数值模拟的计算复杂性，该文提出了非结构网格上周期性边界进出口条件相互回代的方法。应用此方法对横掠叉排管束周期性充分发展层流流动进行了数值模拟，并将其结果与流动进口内8个管束层流流动的充分发展区域的数值模拟结果进行了比较，得到了较一致的结果。实验表明：单周期上进出口条件相互回代的方法可以替代多周期的数值计算，从而简化了计算过程，节省了计算时间。     

稀薄气体高超声流动的非结构DSMC的并行化计算

 采用非结构化三角网格为基本网格单元,在可变硬球（VHS）分子模型、Borgnakke-Larsen唯象模型、Bird的化学反应几率模型及壁面CLL反射模型的基础上,本文用Fortran语言编制了能够模拟内能松弛、热力学非平衡和化学非平衡的稀薄气体直接模拟Monte-Carlo（DSMC）源程序。针对多核计算机上进行并行计算实现技术,将并行OpenMP的模型应用于DSMC方法,编制了可在多核计算机中进行数值模拟的非结构DSMC并行程序。分别对不同稀薄领域的不同工况高超声速气体绕圆柱流动问题进行数值模拟,得到热非平衡态对飞行器流场的影响。通过数值结果的比较,验证了编制的DSMC并行程序的正确性和合理性,以双CPU、双核计算机为例,并行后的计算效率提高了近一倍。这些数值结果对飞行器流场特性分析和有效地完成飞行器热防护具有重要意义。     

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

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