内外环对高速球轴承中球体拖曳阻力系数的影响

在液体润滑的高速球轴承中,球体公转运动会受到润滑流体的拖曳阻力,对球轴承的动力及产热等性能造成显著影响。针对球轴承中球形滚动体公转的流体拖曳阻力问题,建立平板夹层单球体绕流模型,通过数值模拟和实验研究的方法,研究内外环对球体绕流的流动特征及压力分布的影响,分析单球体绕流的拖曳阻力系数变化规律,并与开放空间球绕流模型进行对比。研究表明,建立的数值计算方法和实验结果吻合良好,误差不超过10%。在相同雷诺数下,由于内外壁面的限流作用,球体与壁面接触区域压差阻力明显增大,导致夹层空间球绕流中球体拖曳阻力系数明显高于开放空间球绕流,在103Re105内约为开放空间球体阻力系数的两倍。     

N-S方程的非定常绕流计算——有迎角Joukowsky翼型的粘性绕流

本文利用N-S方程计算了绕Joukowsky翼型的粘性流动,迎角为0°和42°,对称翼型的厚度为25.8％,无量纲时间间隔为t=0到t=10。计算结果包括:流函数分布,翼面压力分布和升力系数以及阻力系数。     

固壁夹道中多体Stokes绕流的数值分析

本文利用格林函数方法,数值求解了固壁夹道中多体的Stokes绕流问题,计算了两平板间的二圆柱体的绕流,三圆柱体的绕流,三正方体的绕流以及五个不规则形状物体的绕流,对其流动特性进行了分析。表明格林函数方法是目前对这类流动问题进行数值分析的最有效方法。     

多孔介质空心圆柱绕流的摄动解

本文主要讨论了绕多孔介质空心圆柱的有势流动问题。运用摄动方法,给出了速度势的四阶近似,然后得到了圆柱表面的无量纲的压力系数和阻力系数。结果发现在势流的情况下多孔介质圆柱受到附加阻力作用,并且在圆柱的喷射或吸收流体的特征速度与来流特征速度的比值是一个小量的时候,阻力系数与特征速度比成简单的线性关系。     

不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

斜、直圆柱绕流的CFD模拟

采用CFD软件CFX5.5对直圆柱和斜圆柱绕流进行数值模拟,计算了直、斜置圆柱的阻力系数、升力系数、周向平均压力系数及其脉动压力系数的均方根值,并与试验结果进行了比较.计算表明,斜圆柱绕流与直圆柱绕流有很大的差异:斜圆柱绕流的平均阻力系数减小,平均升力系数不等于零,升力系数的频域除包含有传统的卡门涡脱频率外,还包含有许多低频成份.     

尖前缘平板的高超声速粘性干扰效应引起的压力分布规律研究

采用CFD方法计算了多种来流条件下高超声速平板绕流的粘性干扰现象,计算与实验结果吻合良好。针对理论分析在平板前缘失1效的情况,结合理论分析、CFD结果和试验结果对原有沿流向壁面压力分布关系式进行了修正,使得修正后关系式在整个平板上与相关数据吻合,显著增加了压力分布关系式的应用范围。同时根据CFD计算结果,分析了平板附面层内压力等流动参数的分布规律,发现附面层内法向压力分布存在峰值,壁面压力分布大于来流静压,并给出了峰值压力及其位置的拟合式。     

谱单元法及其在多圆柱绕流分析中的应用

系统阐明谱单元方法,基于谱单元方法对低雷诺数Re=200时不同间距下的顺排两圆柱和Re=150正方形排列的四圆柱绕流及其阻力系数、升力系数等进行数值模拟。研究比较不同间距比L/D(两圆柱圆心距离与圆柱直径之比)对两圆柱和四圆柱绕流的影响,计算分析了涡量图分布、平均阻力系数和斯托罗哈数随间距比的变化。研究表明,间距比对顺排两圆柱和正方形四圆柱绕流影响显著;顺排两圆柱绕流存在临界间距比,在Re=200时临界间距比约为3.6。正方形排列四圆柱存在三种流态。当流场从一种流动形态变成另外一种流动形态时,力学参数发生显著变化,在某些间距比区间内出现骤升或骤降现象。     

均匀磁场中高超声速弱电离气体流动数值模拟

采用7组元化学模型并应用组分公式计算电导率,通过求解黏性MHD(magneto-hydro-dy-namics)方程组,研究了不同强度均匀磁场对三维钝头体高超声速绕流化学非平衡流动的影响.结果表明,随着磁场的增强,激波脱体距离逐渐增加;总阻力系数和壁面温度逐渐减小.在B_y=0.03 T磁场作用下,与无磁场的结果相比,化学非平衡流中的激波脱体距离增加约7%,总阻力系数减小约5%,局部壁面温度最大降低74%;而冻结流中的激波脱体距离增加约43%,总阻力系数减小约6.9%,局部壁面温度最大降低18%.      

机器学习辅助湍流建模在分离流预测中的应用

数据驱动湍流建模是近年来发展的提高雷诺平均N-S方程预测精度的有效手段,通过机器学习算法能够从高置信度数据中自动提取特征,建立准确的从平均流动特征到雷诺应力的预测模型。针对高雷诺数下积冰翼型绕流这一类典型的复杂流动分离问题,基于此前研究者提出的机器学习预测框架,从输入输出特征选择和翼型绕流中数据分布特性两个方面出发,对机器学习预测结果的光滑性和准确性进行改善。提出了基于雷诺应力张量分析和流动特征辨识的输入特征选择准则;提出了局部区域建模方法以及基准模型和机器学习预测模型混合的代入计算方法。将改进方法应用于积冰翼型绕流问题之中,结果表明改进的方法能够准确给出训练集和预测集上的雷诺应力结果,并且代入平均流计算可以得到和真实分布更加接近的流动和机翼表面压力分布。     

基于2D-LDV技术的船舶螺旋桨三维流场结构分析(英文)

利用2D-LDV设备,通过优化布置测试点,获得了均匀来流下螺旋桨尾流场三维流动数据。展示了宏观的流场动态以及微观的流动信息,如尾涡和梢涡在流场中的分布情况、尾涡部流动的剪切效应、梢涡部的绕流及流动分离状态。计算了水动力螺距角以及螺旋桨附着涡的总环量。试验显示了LDV技术在流动细节测试方面的优势,获得的定量信息可以为相关研究人员提供技术支持。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

稀薄流到连续流的气体运动论统一算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数的各流域均适用的气体分子速度分布函数简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而在速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。发展可用于速度空间宏观取矩的离散速度数值积分方法，获取物理空间各点的流动参数，由此发展一套能有效模拟各流域三维绕流问题的气体运动论统一算法。研究气体运动论数值算法所适合的并行方案，基于统一算法的HPF并行实现，建立一套能有效模拟不同流域复杂外形体绕流的HPF并行算法软件。通过对不同Knudsen数的一维、二维、三维气体绕流问题进行数值计算表明，计算结果与有关实验数据及其它途径得到的研究结果吻合较好，证实了本文发展的统一算法在求解稀薄流到连续流不同流域复杂绕流问题方面的可行性。     

空腔流动特性数值模拟研究

本文采用数值方法模拟研究了空腔的长深比（L／D）对流动特性的影响。控制方程采用Navier-Stokes方程，计算方法采用有限体积法，湍流模型采用Baldwin-Lomax模型。通过分析三种典型空腔的流线图谱、压力云图、底部压力分布曲线，得到空腔流动特性与其几何参数和来流参数间的相应关系。从数值模拟的结果来看开式空腔流动的空腔内压力分布均匀，流场结构相对稳定，有利于内埋式武器安全平稳地分离。     

用二维N—S方程计算汽车流场的地面效应

采用有限体积方法求解了Navier-Stokes方程,计算了绕二维类汽车物体的绕流并研究了地面效应。在非正交网格系统上以强守恒形式构造出有限体积方程,采用压力校正算法来求解所得到的方程。对于湍流情况,采用了双方程模型的κ_ε模型。 为了验证采用本程序计算移动壁面效应的正确性,本文给出了Couette流动的计算结果。 本文也给出了汽车绕流的计算结果,它包括层流情况,从汽车头部开始由层流转为湍流的混合流动情况及湍流情况。对于每一种情况都比较了移动壁面的计算结果和固定壁面的计算结果。     

有限体积格子Boltzmann方法用于近空间连续流区绕流模拟

针对航天器再入解体形成残骸碎片的近空间绕流计算问题,拓展了格子Boltzmann方法在可压缩流动模拟的能力。引入有限体积隐式格式求解耦合双分布函数格子Boltzmann模型方程及采用圆函数为基础构造的D2Q13离散速度模型;引入IMEX-RK格式进行了时间项离散解决模型方程的源项刚性问题;对Riemann问题、平板双马赫反射问题、RAE2822翼型跨声速绕流等近空间连续流区、可压缩典型案例进行了数值模拟。通过比较分析,初步验证了耦合双分布函数有限体积格子Boltzmann方法对连续流区、可压缩流动的模拟能力;进一步开展了方柱形解体残骸的超声速绕流模拟,得到与N-S方程计算结果吻合一致的绕流结果,证实经改进的耦合双分布函数有限体积格子Boltzmann方法对解体残骸碎片绕流问题具有较好的模拟能力。     

过渡区球锥气动力特性的桥函数估算与DSMC仿真研究

采用DSMC(DirectSimulationMonte-Carlo)方法模拟球锥体过渡区流动。论文给出了不同Knudsen数球锥绕流流场的温度与压力分布图,以及在不同Knudsen数、不同马赫数条件下球锥阻力系数和底部阻力系数的分布,并将阻力系数模拟结果与桥函数一体化估算结果进行了比较分析。     

穹顶结构表面平均压力分布及绕流结构数值模拟

通过商用软件Fluent求解Realizable形式的时均k-ε湍流模型,得到了穹顶结构煤场的绕流结构及表面平均压力分布,并结合流动分离、旋涡结构的特点对其平均风压分布规律进行了解释。最后通过参数变化,研究了来流边界条件、煤场钢混结构高度以及煤场侧窗开启状态对煤场顶部穹顶结构表面平均风压分布的影响。     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.