基于动态嵌套网格方法的摆动翼型粘性绕流数值模拟

该文根据摆动翼型的流场特性运用动态嵌套网格技术将其划分为不同的子区,在不同的子区中使用不同的流动控制方程,以提高计算效率.对翼型生成一个随其摆动的相对较小较密的贴体网格,在它上面采用Navier Stokes方程;对整个流场生成了一个大的相对较稀的静止背景网格,在它上面采用Euler方程.两个网格间的流场信息交换通过动态嵌套网格技术来实现.通过与单块翼型网格的计算过程和计算结果的对比分析,表明该方法在不损失计算精度的前提下计算非定常问题时具有相对较高的计算效率.     

基于预处理技术的扰动风对多段翼型非定常绕流作用的数值模拟

运用双时间推进法和动态嵌套网格技术对扰动风作用下的多段翼型非定常绕流进行了数值模拟. 首先, 结合预处理技术对麻省理工学院(Massachusetts Institude of Technology, MIT) 的翼型模型在扰动风作用下的流动进行数值模拟, 并运用其实验结果验证程序的正确性; 然后, 对多段翼型在周期性的扰动风作用下的非定常粘性绕流进行数值模拟, 得到非定常情况下气动力随扰动风变化的迟滞规律. 计算范围内的结果表明: 扰动风的频率越低, 幅度越大, 多段翼型整体上的气动力越大; 扰动风的频率越高, 幅值越大, 对多段翼型的气动力影响越大.     

预处理算法在非定常流动计算中的应用

介绍了一种可应用于全速度范围的非定常N-S方程数值求解方法,并全面考核该方法的适用性。该方法通过引进伪时间导数项,对伪时间导数进行预处理以解决低速时方程系数矩阵刚性过大的问题;同时用双时间步推进算法保证时间精度,将该方法推广应用到非定常问题的求解。文中给出了若干算例,覆盖了从极低马赫数到超音速、无粘/层流/湍流、二维/三维、定常/非定常情况。算例结果表明,该方法在较宽的使用范围内均能得到理想的计算结果。     

基于CFD技术的螺旋桨非定常流场数值模拟

针对单独螺旋桨模型,基于动态结构网格搭接技术,采用CFD方法对其非定常流场进行了数值模拟。首先针对全部外形生成多块点搭接和面搭接两套计算网格,通过分析两套网格的计算结果,验证动态网格搭接技术的可行性和正确性,得出螺旋桨性能参数随转速的变化规律;在此基础上分析了雷诺平均NS(RANS)方程和Euler方程的数值模拟结果;最后研究了攻角对单独螺旋桨性能参数的影响。研究结果表明：动态网格搭接技术能够很好地模拟非定常螺旋桨的性能参数;螺旋桨的拉力系数和功率系数都随着转速的增加而增大;对螺旋桨非定常流场进行数值模拟时,采用Euler方程就能够满足工程精度的需求;攻角0度时单片桨叶的拉力系数和转矩系数都是恒值,而在攻角10度时却都呈现出周期性的正弦曲线。     

粘性层流流速度分布的数值计算方法

介绍了如何运用数值计算方法来求解粘性层流流体的速度分布,通过坐标系转化理论、连续性方程、Navier-Stokes方程导出有关速度的非线性微分方程,运用数值计算方法分别对对称性速度流体和非对称性速度流体进行数值求解,其结果对深入了解粘性流体运动的基本规律有着重要的现实意义.     

二阶迎风有限体积法方腔流数值模拟

介绍了一种基于非结构网格的二阶迎风有限体积离散格式，在对流项的离散过程中，为了达到二阶精度，在界面上对物理量咖作Taylor展开，在处理展开项中的梯度时经入NND格式的优点，克服了中心差分格式不稳定的缺点，并对方腔流动进行了系统的数值模拟．计算网格采用三角形网格，节点数为12960，单元数为25600．压强修正基于SIMPLEC方法．给出了雷诺数Re达到5000时的定常流动结果，与以往方腔流计算的标准解非常吻合，但所用的网格数要少．计算结果表明该二阶迎风有限体积算法有很好的收敛性和稳定性，且松弛因子对计算结果影响很小．     

多段翼型无网格算法及其布点技术

用显式无网格算法实现了多段翼型的数值模拟。与传统的网格算法不同,求解区域用“点云”离散代替通常的网格划分。基于当地点云离散结构,用二次极小曲面逼近计算空间导数。在研究该算法的基础上,给出了Euler方程无网格离散形式,运用Runge-Kutta显式时间推进格式推进求解。此外,还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了复杂多段翼型的绕流。     

基于动态嵌套网格技术的飞行器导弹发射的数值模拟

该文采用基于动态嵌套网格技术上的双时间法,模拟了飞行器导弹发射时的粘性流场和导弹气动性能的变化,湍流模型为K-ω两方程模型.该文研究的是一个简单的二维模型,但所得的结论有助于理解和掌握激波对飞行器发射导弹的飞行速度、飞行稳定性和飞行轨迹的影响,并具有一定的参考价值.     

低雷诺数下翼型粘性绕流主动流动控制的数值模拟

通过数值求解Navier-Stokes方程模拟零质量射流对低雷诺数翼型绕流的控制作用.空间离散采用中心有限体积格式,时间推进为双时间推进方法.计算结果表明,采用零质量射流进行主动流动控制能有效抑制大攻角下大尺度的流动分离,改善翼型的气动特性 ,起到显著的增升减阻作用.     

Benard对流问题五模类Lorenz方程组混沌吸引子的数值模拟

考虑底部加热的两平行板间流层的Benard对流问题,该问题可以用Navier-Stokes方程与热传导方程来描述,对其速度场和温度场等相关变量进行二维傅立叶展开,经过复杂运算得到新五模类Lorenz方程组,对此方程组的动力学行为进行了数值模拟     

基于动态嵌套网格技术的微型共轴式双旋翼流场的数值模拟

运用动态嵌套网格技术和双时间推进算法,对微型共轴式双旋翼的非定常粘性绕流进行数值模拟和研究.针对上、下两层旋翼间距离较小和流动粘性影响较大的特点,在动态嵌套网格技术的基础上引入滑移网格技术,保证了流动信息交换的准确性.首先,模拟低雷诺数下单旋翼的流动,结果与实验数据吻合较好,验证了算法的有效性;然后,利用该方法实现对微型共轴式双旋翼流场的数值模拟.结果表明,微型共轴式双旋翼流场中,桨尖涡起主要作用,其以螺旋方式向下运动,同时与周期性运动的旋翼相互作用,使得作用在旋翼上的总体拉力呈现出周期性的变化规律.     

一种求解不可压N-S方程的非结构化网格方法

将Date压力修正算法推广到非结构化网格中来处理压力与速度的耦合,用具有二阶精度的迎风最小二乘格式离散对流项,用梯度投影法离散扩散项,得到了一种在非结构化网格上求解N-S方程的同位网格法.它充分利用了单元梯度,避免了利用顶点信息计算交叉扩散通量的麻烦,所有界面速度都采用线性插值,克服了Rhie＆Chow动量插值法的缺点.算例表明,该算法精度较高且健壮稳定,收敛特性也比较好.     

风对建筑物体绕流的数值模拟

利用有限分析法数值模拟了风对建筑物体的绕流，对高雷诺数问题引进了一个经验公式较好地克服了实际中的困难，对方形、圆形及两头尖等各种形状物体绕流的计算结果表明数值方法在准确模拟各种形状建筑物体的绕流流场方面是成功的。     

三维流体诱发振动的理论模型与数值模拟

建立了三维流体诱发振动的理论及数值模型，采用任意拉格朗日-欧拉框架下的三维非定常Navier-Stokes方程组及波动方程组，分别描述了流体的流动及结构振动．流动及振动控制方程组分别采用基于非结构化运动网格上的有限体积法和有限差分法求解．对刚性、受限振动及自由振动的三维圆柱体进行了数值计算．结果表明：圆柱的自由振动基本上为Strouhal型涡诱发振动，即每个截面按“8”字型轨迹运动，且横向振幅远大于流向振幅，计算结果与现有的实验及数值解吻合较好．首次采用流体诱发振动模型对极限环振动模式进行了定量分析，发现在圆柱中间截面附近高频气动力对振动起阻尼作用，而在端面附近气动力与位移间的稳定相位差对振动起激励作用．     

冲刷坑冲刷过程的数值模拟

为了探索用数值模拟的方法来模拟冲刷坑形成过程及预测冲刷深度的可行性,在贴体非正交曲线坐标系中,采用控制体积法对二维散粒体河床的冲刷进行了数值模拟,根据冲刷过程中坑底动水压强和冲刷深度之间的平衡关系,来控制坑底边界网格移动进而模拟冲刷坑的冲刷过程及流场变化过程,实现了对冲坑底边界动态发展过程的自动模拟.计算得到的控制参数变化规律,表明用此方法模拟冲坑发展过程的可行性;计算结果与试验吻合较好,说明本方法可以有效地用来模拟冲坑发展过程,并能正确的反映挑流冲刷中流场结构和速度、压力等参数的变化.     

液滴溅落的数值模拟

将二阶投影方法和水平集技巧相结合求解二维Navier-Stokes方程,以模拟水溅落的过程。与其他方法不同,本方法不需额外的技巧,就能自然地处理自由界面出现的尖点,角点等。数值计算过程中无重新初始化步骤。