平板翼流体动力边界积分多重网格计算

讨论了平板翼流体动力边界积分法多重网格计算模型,给出了残差计算格式和函数修正值计算公式,推导出与步长无关的环量多重网格计算公式.计算结果表明,利用多重网格方法加快了收敛速度,其数值计算的流体动力与试验值吻合.     

面向复杂几何模型的多级并行四面体网格生成算法

高性能计算机的快速发展为数值模拟提供了必需的硬件环境,数值模拟领域对网格的需求已高达数十亿到数百亿量级,而网格生成作为数值模拟的重要一环,其发展则相对滞后,很难满足并行数值模拟求解器对大规模网格的需求。本文面向复杂几何模型提出一种多级并行四面体网格生成算法,该算法首先基于模型的几何特征建立网格的尺寸场,并基于尺寸场和几何实体间的邻接关系对几何实体进行分组,将分组后的几何实体分配到不同的计算节点,在计算节点间采用前沿推进法实现三角形面网格的并行生成,然后在计算节点内对三角形面网格进行二级区域分解,将分解后的子网格分配到各进程中,最后在进程内采用多线程并行方法实现四面体网格的并行生成。通过实际应用三峡大坝模型进行验证,该算法具有良好的并行效率和可扩展性,可以在数千处理器核上实现十亿规模高质量四面体网格的并行生成。     

坦克动力舱空气流动计算区域的网格生成

为了合理解决网格划分对坦克动力舱空气流动与传热三维数值计算效率和计算精度的影响,特提出动力舱空气流动计算区域的网格生成方法.对某型坦克动力舱的网格生成实例,显示了这种方法的实用性.     

一种解Euler方程的高效算法

本文介绍了目前解Euler（欧拉）方程的一般情况,认为采用多重网格对角隐格式,在块结构网格上解无粘可压流EUler方程是一种高效算法。引进多重网格循环的二种模式,可以考核收敛率、精度和效率;把这种算法运用在在并行机上,计算了绕翼型跨声速流场,结果所给出的压力分布和收敛历程验证了这种算法的精度、收敛率;本文还讨论了在并行运算中,用几个处理机所起的效果。     

滑移边界条件的收敛性分析及应用

用速度滑移与温度跳跃边界条件代替通常假定的无滑移边界条件,可有效地提高计算流体力学模型对高空滑移流区域流动的预测精度。应用Maxwell滑移边界条件时,通过直接计算速度梯度及温度梯度而得到速度滑移和温度跳跃量的处理方法在网格较密的时候会出现迭代计算发散的问题。理论分析表明,直接计算梯度的方法使边界条件的时间推进过程等价于雅克比迭代过程,因此必须满足相应的收敛性条件。为了消除收敛性条件的限制,给出了一种在任意网格密度下均收敛的边界条件处理方法并通过数值算例验证了该方法的正确性。针对高空高超声速流动,以空天飞机为例,对比了滑移/无滑移边界条件所得结果的差异,分析了滑移效应对飞行器气动特性及热环境的影响。     

复杂平面区域的三角网格生成算法

为生成复杂平面区域的有限元网格,提出了基于网格细化的三角网格生成算法.该算法首先采用耳尖移除法对区域边界做三角划分,得到粗略的初始网格.提出Delaunay优化平分方法,根据网格密度细化初始网格,该网格细化方法结合最长边平分技术与Delaunay边交换技术,可有效提高内点生成与单元细分的质量.实验表明,基于Delaunay 优化平分的三角网格生成算法可对任意平面域进行网格剖分,生成符合有限元计算要求的高质量三角网格.     

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

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

导管螺旋桨敞水性能的预报和比较

借助于计算流体力学软件,对导管螺旋桨的敞水性能进行了数值模拟,得到了在不同网格模型和湍流模型下导管螺旋桨的正车敞水性能曲线,通过与试验图谱的对比分析,发现Wilcoxk-ω模型更适合于导管螺旋桨的敞水性能计算,并且采用结构化网格和非结构化网格相结合的计算方法能够满足导管螺旋桨敞水性能预报的工程精度要求。同时,当网格中存在少量高度倾斜的网格单元时,在不影响计算收敛的情况下,仍能将计算误差控制在10%以内。此外,如果要获得更为精确的计算结果,应提高网格质量,尽量使用结构化网格,并将近壁面网格加密,合理控制壁面附近的Y+值。虽然网格数量的增加并不总是意味着计算误差的减少,但合理控制网格细密度能够获得更为可信的计算结果。     

跳跃-滑翔弹道扰动引力自适应网格快速赋值方法

为实现高超声速跳跃-滑翔弹道扰动引力的快速赋值,提出自适应网格赋值模型,并根据反距离加权理论,优化广义延拓逼近算法,对模型的逼近误差进行分析。该赋值模型的网格划分为两级,第一级网格根据标准弹道空域进行划分,第二级网格根据滑翔导弹实际弹道在线生成。根据一级网格节点数据,通过优化广义延拓逼近算法计算二级网格节点数据,最后根据二级单元内插计算实际弹道点的扰动引力值。仿真结果表明:在同等大小的网格划分下,优化广义延拓自适应网格模型的逼近精度高于一般赋值方法;在同等精度要求下,该赋值模型的最大单元格边长大于一般赋值方法,从而减少了单元格划分数量,进而降低弹上数据存储量;针对不同滑翔方向以及不同滑翔距离的跳跃-滑翔弹道,该模型逼近误差对应的落点偏差小于5 m,具有较好的适应性。该赋值模型在满足计算速度的前提下,提高了传统赋值方法的逼近精度,降低了弹上存储量,具有一定的工程应用价值。     

弹塑性流体计算中的一种实现方法

针对弹塑性流体计算问题，在ALE方法基础上提出了一种数值模拟的实现途径，主要包括含弹塑性的流体力学数学模型的建立，滑移面处理的方法，网格生成等，最后给出计算结果。     

影响柴油机喷雾数值模拟精度的若干因素分析

为探明影响柴油机喷雾数值模拟精度的若干因素,采用通用商业软件FLUENT对柴油机喷雾特性进行了CFD模拟计算,分别研究了网格尺寸、喷雾粒子数和最大时间步长对油束几何形状及油束贯穿距的影响.结果表明,网格划分的精度和喷雾粒子数量的多少对计算结果的正确性和油束的几何形状有直接的影响,网格尺寸过大或喷雾粒子数过少则喷雾几何形状失真.最大时间步长对数值模拟的精度在一定的范围内影响是比较小的.     

高超音速流场的隐式TVD格式并行计算

本文采用区域分割技术和拼接网格的并行策略，发展了一个适合于分布式存贮多机系统的ＴＶＤ隐式有限体积并行算法；并在ＰＶＭ并行环境下，对三维高超音速绕流流场实现了多机并行计算，通过负载平衡等方法得到了较高的加速比（在二处理机系统上加速比为１∶８４，在四处理机系统上为３∶４４）。     

水力旋流器湍流场数值模拟

应用CFD软件FLUENT对水力旋流器进行三维数值模拟，并对网格生成技术、湍流模型、离散方法及边界条件等问题进行了探讨。数值模拟所获得的轴截面上的压力场和速度分布规律与前人的研究结果基本吻合，验证了数值模拟的可行性。为进一步研究水力旋流器的结构优化和分离性能预测提供参考。     

Poisson方程生成数值网格方法中的源项计算法

网格生成是数值研究方法得以实现的基础。本文采用了Ｐｏｉｓｓｏｎ方程生成数值网格，并针对方程源项的选择问题，对两种形式的源项作了比较研究。     

迎风格式在接触间断的数值耗散及其诱导误差

在不同流场参数条件下对三种迎风格式在接触间断中的数值耗散问题进行了数值实验,并对数值耗散产生的机理进行了分析。数值计算结果和理论分析表明,矢通量分裂格式计算接触间断问题时,若流场静止或流场内存在亚声速区域,密度耗散的产生会诱导出以特征速度运动的数值扰动误差,该误差对数值耗散的大小无影响,但会影响流场的速度及压力参数分布,从而改变流场的结构。在二维问题中,诱导误差相互干扰会产生大量的复杂小尺度结构,给流场结构分析带来困难。同时,在密度参数线性分布的流场中,若空间离散格式重构的对象为对流通量,使用矢通量分裂格式计算流场会产生数值误差,使计算精度难以到达二阶。     

非结构网格上解二维Hamilton-Jacobi方程的自适应算法

提出了一种非结构网格上解Hamilton Jacobi方程的自适应局部加密方法。讨论了加密方法关键的一步:光滑度指标的确定,通过采取无振荡插值,得到了无振荡数值格式,典型算例表明,该方法只需增加少量的点和计算量却能够提高计算精度和对间断的分辨能力。     

全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢.原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩...     

自润滑弹带挤进过程中渗出油量的对比计算

从弹带的作用及特性出发, 通过数值计算和统计计算 2种方法计算了自润滑弹带挤进过程中的渗出油量,初步得出了渗出油量和行程之间的关系。这一问题的解决为自润滑弹带材料的选择提供了理论依据。通过对 2种算法的比较可知, 数值计算误差更小一些, 而统计计算对数学条件要求不太苛刻。