原始变量法计算封闭腔内自然对流

在文献[1，2]的基础上，利用原始变量法求解封闭腔内自然对流换热问题．对求解域采用非等距交错网格剖分，利用泰勒级数于网格点展开取二阶精度进行(u,v,θ)方程各项的离散，采用SIMPLE方法对压力及压力修正进行求解．不同瑞利数(Ra)条件下的数值试验显示，该数值方法物理概念清晰，计算稳定,可通过调节网格的疏密等方法，达到改善求解的收敛及提高求解的精度等目的．     

湖泊整治工程二维水沙数值模拟

应用破开算子法原理,将平面二维不恒定流及泥沙扩散基本方程按不同的物理意义破成对流、扩散和传播三个子方程．基于三角形网格划分,利用有限体积守恒格式求解连续方程,用特征线方法求解对流方程并采用有限差分法求解扩散和传播方程,经数值实验及实际工程验证,该算法计算简便、边界符合良好,数值模拟精度也较高．     

同位网格下平面二维水流数学模型的有限分析解法

为避免在求解平面二维水流数学模型方程时,流速、水位分别采用有限分析法和有限差分法求解而导致运动方程和连续性方程的计算精度不匹配,进而降低模型整体的计算精度和收敛速度.参照SIMPLE算法,推导出了能够在同位网格下采用有限分析法求解的一般曲线坐标系中的水位控制方程,实现在同位网格下采用有限分析方法求解全部控制方程的同精度算法,避免了交错网格,提高了计算精度,加快了收敛速度.运用于实际工程,计算结果与实测值基本吻合.     

非结构网格特定条件下热传导的数值计算

用基元有限容积方法,在非结构化网格中离散和求解热传导方程,二阶迎风格式和全程隐式迭代求解稳态导热问题.在验证求解方法的算例中,结构化和非结构化网格同时采用,并且提供了3种非结构化网格(三角形网格、四边形网格和混合网格)的计算结果,同时用精确解与得到的数值解相比较.虽然两种网格形式都能得到满意的精确度结果,但对于非正交的非结构化网格,二次扩散项对提高精度是十分重要的.     

高阶间断有限元方法求解Euler方程时的数值积分问题研究

文章对非结构网格上高阶间断有限元方法求解Euler方程时的数值积分方式进行了研究。首先根据间断有限元方法的基本原理推导出了有限元离散控制方程中线积分和体积分项的数值积分精度要求。然后给出了采用Gauss-Legendre和Gauss-Lobatto积分公式处理线积分项,以及采用Guass积分公式和重构积分方法处理体积分项的情况下,为满足积分精度要求所需使用的最少积分节点数目。最后,通过具体算例对上述积分精度要求进行了验证,并考察了不同数值积分方法对于求解效率和精度的影响。     

一维对流扩散方程的迎风变换及其有限元解

提出了两种迎风变换函数．利用这种变换函数将一维对流扩散方程进行变换，得到等价的扩散方程．然后使用有限元（ＦＥＭ）求解等价的扩散方程，得到对流扩散方程的数值解．计算成果表明，用这种方法计算．可突破网格Ｐｅ≤２的界限，可求得较大和更大的Ｐｅ数情况下的数值解．解的精度和收敛性是令人满意的．     

一种隐式时间算法在非结构混合网格粘性流动计算中的应用

进行了一种隐式时间算法在非结构混合网格求解粘性流场中的应用研究，采用计算量较小的中心有限体积格式对N-S方程进行空间离散，运用通量线性化假设和最大特征值方法进行雅可比矩阵分裂来实现隐式高斯-赛德尔时间迭代，并利用作者提出的自适应当地时间步长和隐式残值光顺等措施加速收敛。以M6机翼和DLR-F6带发动机短舱的翼身组合体2种外形的跨音速流场为数值算例，计算结果表明：改进的非结构混合网格粘性流场隐式求解算法具有精度好、效率高的特点，可用于激波——附面层干扰、飞行器部件干扰等多种复杂流动的计算。     

湖泊整治工程二维水沙数值模拟

应用破开算子法原理,交坪面二维不恒定流及尼潲扩散基本方程按不同的物理意义破成对流、扩散和传播三个子方程,基于三角形网格划分,利用有限体积守恒格式求解连续方程,用特征线方法求解对流方程并采用有限差分法求解扩散和传播方程,经数值实验及实际工程验证,该算法计算简便、边界符合良好、数值模拟精度也较高。     

非结构网格迎风格式中的一种重构方法研究

研究了一种基于网格节点重构的二阶迎风格式在非结构网格无粘流场计算中的应用，提出了一种新的网格节点流场参数的重构方法。并通过ONERA M6机翼、带副翼的直机翼、翼身组合体以及翼身组合体带增升装置无粘绕流流动计算对该重构方法进行了验证。算例表明新方法可以在确保格式精度和不显著增加流场计算时间的前提下，明显改善迎风格式的稳定性、收敛性和鲁棒性，减小了网格质量对流场求解精度的影响，因此十分适合非结构网格计算复杂外形流场时使用。     

新型低阻导线气动特性数值模拟

采用流体力学计算软件FLUENT6．3,对新型低阻导线LP810导线进行雷诺数76000时气动力特性以及流场的数值模拟。计算采用有限体积法和SIMPLEC算法求解均匀黏性不可压缩流体的NaV1er—St0keS方程（N—S方程）,采用三阶精度二次迎风插值（OUICK）离散N—S方程,以降低求解过程的数值振荡。计算区域网格划分采用具有良好形状适应性的非结构化三角形网格,贴壁网格尺寸取为约0．0005倍的导线直径,网格尺寸在导线周围直径4倍于导线直径的圆形加密区内沿半径方向按照1．009的比例增大。计算结果表明,计算方法具有良好的计算精度,与试验结果吻合较好。新型低阻导线LP810导线的特殊截面设计,使得导线在高雷诺数时的分离点显著推后,降低了在高雷诺数时的阻力系数。     

基于直角叉树切割网格的复杂外形N-S方程数值模拟

对于粘性绕流的数值模拟，在直角叉树切割网格的基础上，结合三角形非结构网格和矩形结构化网格，利用其各自的优势和特点，提出了一种生成混合染交网格的思路和方法，在物面附近生成了适合粘性流计算的大长宽比矩形网格，地远场分布直角叉树切割网格，快速离散计算空间。对于复杂的多体问题，采用三角形网格来连接 各体网格，并运用网格合并的方法，保证了各网格之间的光滑过渡过连接，提高了网格的质量，针对复杂多段翼型绕流问题，在上述网格基础上，采用B－L代数湍流模型和中心有限体积法，完成了Navier-Stokes方程（N－S方程）的数值模拟，计算结果表明，网格生成和流场计算都是正确的。     

贴体网格的生成方法及其在天然河道数值模拟中的应用

由于生成网格质量的优劣直接影响流场计算的精度和效率,在对目前常用的几种网格生成方法进行比较的基础上,提出了代数法和偏微分方程相结合的一类新的混合方法。该方法能有效地控制所生成网格的疏密,且在边界处有良好的正交性。通过具体算例,说明该生成方法具有良好的可靠性和较强的适应性等特点,可以用于解决天然河道、水库等不规则边界网格生成问题。     

复杂边界正交曲线网格生成技术

在采用Poisson方程进行曲线网格生成时,如何确定合适的调节因子P、Q函数是一个重要的研究内容.分析了Thompson的曲线网格生成方法中控制网格分布的调节函数的缺陷,经过严格推导得出了一组新的调节函数P、Q的表达式,并采用新调节函数生成曲线网格.网格实例表明,该调节函数能够对复杂边界的单连通域或多连通域生成理想的正交曲线网格,网格的尺寸与形状连续光滑地变化,网格线的光滑性和间距得到了有效的控制,网格分布能够适应物理量场变化的情形;正交曲线网格能够在实际水利工程流场数值模拟中准确地使用边界条件,从而提高其求解的精确度.     

不同网格加密方法在河流回补地下水模拟中的对比

河流回补地下水是线状回补过程,河道附近地下水位高精度模拟是回补效果准确评估的关键。基于非结构化网格加密方法,设定4种不同的网格剖分方案(两种结构化网格和两种局部加密的非结构化网格),在不同边界条件下对河流回补地下水过程进行模拟。通过对末流场、距离河道不同范围回补前后水位变差、地下水均衡和模型运行效率的对比分析,阐明线状回补模拟中不同网格加密方法的优劣。结果表明:与粗结构网格模型相比,两种局部加密的非结构化网格模型以及细结构网格模型的模拟精度高,能更准确刻画因河道回补导致的地下水位变化; 两种局部加密的非结构化网格模型比细结构网格模型运行时间短,在没有大幅降低模拟精度的条件下,运行效率更高; 四叉树网格模型相对于嵌套网格模型具有长期模拟的优势,在河流回补期末,模拟精度要稍高于嵌套网格模型; 同时,四叉树网格模型具有多级加密的优势,加密级次越多,剖分的网格数量相对越少,运行时间相对更短。综上所述,四叉树网格加密方法既保证了模拟精度,又一定程度上提高了模型运行效率,是值得推广的一种局部网格加密方法。     

渗流计算中的三维非结构化网格下的自由面适应网格方法

在2维正交网格系统下自由面适应网格方法的基础上,提出了3维非结构化网格下的自由面适应网格方法.该方法在迭代求解渗流自由面时,无需改动网格,而是用映射的方法去让自由面适应既定网格系统,也不需要对出逸边界进行特殊处理,其方法简单,易于程序实现;通过将数值计算结果与甘油试验及室内砂槽模拟结果相对比,表明了该方法在3维非结构化网格系统中捕捉浸润面的精度较高,拓宽了自由面适应网格方法在3维带自由面渗流问题中的应用范围.     

基于非均匀网格的高效高精度洪涝过程模拟

随着计算机技术的发展,网格划分技术日趋成熟.针对变化环境下洪涝灾害频发以及实际地形范围广泛,且通常存在狭长沟道与广阔泛洪区等问题,本文建立了一套具有层次拓扑关系的结构化非均匀网格划分模型,并与基于GPU加速的高分辨率数值模型相结合模拟地表水流动过程.高质量网格影响着模型计算精度与计算效率,非均匀网格划分技术以地形高程的...     

复杂几何模型数值计算网格划分及独立性考核

以开孔折流板矩形通道为例,对流体流经复杂几何通道进行数值计算时网格的划分方式进行了探讨,研究了对于同一计算模型,在不同工况下采用相同的网格划分时,网格独立性的考核标准。通过设置有两种类型开孔折流板通道的数值计算表明:对于同一计算模型,在不同工况下采用相同网格进行数值计算时,网格独立性的考核应该在高Re数工况下进行.     

阶背式MIRA模型气动阻力数值模拟网格无关性研究

对汽车气动阻力数值模拟结果影响较大的网格划分形式和方法进行了对比研究.利用现有充分的计算资源,建立阶背式MIRA模型,通过划分常见的四面体网格和混合网格,分别对两种网格形式在整车周围和车身尾部对网格进行加密,并保证足够多的网格数量,得到了数值稳定的网格无关性的阻力系数结果,总结出了阻力系数随网格加密的变化趋势,并与现有风洞试验值进行对比分析.结果表明：网格一旦加密后阻力系数即有显著下降,随着网格数的不断增加,阻力系数均呈现出明显的下降趋势,且下降幅度逐渐减小,最后都能在3000万网格左右得到网格无关性的结果.最终四种方案的阻力系数值均小于风洞试验值,且与之存在一定的误差.     

非结构网格生成技术中一种直接提供背景信息的方法

本文利用点源或线源所影响的势函数分布规律,结合计算流体力学中外部绕流的特点,发展了一种直接提供背景信息的方法,用来控制非结构网格生成过程中的网格分布;并且通过对部分点源或线源的影响势函数施加方向控制开关,使所生成的网格按期望的规律分布。使用证明,本文方法是一种有效、快速的处理方法,与传统用背景网格提供背景信息的方法相比,在计算量、程序量、光滑性和灵活性方面都能表现出良好的品质。     

翼身组合体的网格生成及EULER方程数值模拟

采用代数方法中的广义无限插值法生成翼身组合体三维Ｏ－Ｃ型贴体网格。该方法使用双边界法和引入了物面外法向导数以保证物面附近网格正交的处理方法，提高了代数方法生成翼身组合体网格的质量。并用中心格式有限体积法求解跨声速ＥＵＬＥＲ方程，以某翼身组合体模型为例，计算结果与实验符合良好。