梯形翼风洞试验模型数值模拟技术

基于雷诺平均Navier-Stokes (RANS)方程和结构网格技术,采用二阶空间离散精度的MUSCL格式,并结合k-ω剪切应力输运(SST)两方程湍流模型和γ-Re_θ转捩模型,研究了梯形翼风洞试验模型中前缘缝翼、后缘襟翼连接装置对气动特性的影响。简要介绍了本文采用的计算方法;介绍了梯形翼的风洞试验模型及风洞试验结果;在网格收敛性研究的基础上,采用"全湍流"方式和转捩模型研究了梯形翼试验模型连接装置对气动特性的影响。通过与不带连接装置的计算结果的对比,采用"全湍流"模拟方式,计算模型中考虑试验模型的连接装置引起升力系数下降、阻力系数下降、低头力矩减小以及失速迎角提前;通过与试验数据的对比,进一步考虑转捩影响可以提高梯形翼风洞试验模型气动特性的计算结果与试验结果的吻合程度,梯形翼风洞试验模型失速迎角附近的气动特性数值模拟技术还需要进一步的研究。     

第一届航空CFD可信度研讨会总结

总结了第一届航空可信度研讨会(AeCW-1)的数值模拟结果。AeCW-1会议由中国空气动力学会计算空气动力学专委会主办,并在第一届空气动力学大会上做了专题研讨。本次研讨会吸引了14家单位参与,共提交了26组计算结果。活动的主要目的是针对典型运输机巡航构型,评估国内CFD当前技术状态,探索CFD进一步发展的研究方向,促进CFD验证和确认工作的开展,为大型客机的研制提供技术支撑。AeCW-1选择的研究模型为中国空气动力研究与发展中心(CARDC)设计的单通道运输机模型(CHN-T1),CHN-T1包括机身/机翼/平尾/立尾等部件,分别在CRADC的2.4m×2.4m风洞(FL-26)、欧洲的2.4m×2.0m风洞中(European Transonic Wind Tunnel,ETW)开展了风洞试验。AeCW-1组委会提供了系列结构网格和非结构网格,研究工况主要包括网格收敛性研究、抖振特性研究、雷诺数影响研究,风洞模型支撑装置和静气动弹性变形对数值模拟结果的影响是重要的研究内容。本文对提交组委会的数值模拟结果进行统计分析,并与FL-26及部分ETW的风洞试验结果进行了对比,对未来的研究工作提出了意见和建议。     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

离心泵双密封环平衡系统泄漏特性预测方法研究

为了准确预测离心泵双密封环平衡系统的泄漏特性,通过分析该平衡系统的几何特征和液体流动特点,对其主要过流部件进行简化,将泵腔液体流动简化为由圆周剪切流和径向压差流叠加而成的轴对称二维粘性层流运动,同时将密封环间隙液体流动简化为两平行平板间的定常不可压缩粘性流体层流运动.然后,利用N-S方程分别建立了泵腔和密封环间隙的液体...     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

柔性机器人动力学建模的一种方法

提出一种建立具有柔性关节的多柔杆机器人动力学方程的方法.首先给出了柔性关节及柔性臂杆的简化模型,然后利用Kane方法和假设模态法推导出完整的动力学方程并给出了递推公式.在此基础上利用一个算例研究了臂杆柔性与关节柔性对机器人动力学响应的影响.结果表明:所建立的动力学模型更接近于实际的柔性机器人系统;臂杆柔性与关节柔性都起着非常重要的作用,因此在柔性机器人的动力学建模与控制中应充分考虑它们的影响.      

基于N-S方程的尾迹面法翼型气动阻力计算

研究了计算翼型阻力改进的方法——尾迹面法及尾迹面积分的位置和相应的积分技术.在亚跨声速时分别采用了表面积分和尾迹积分求得给定翼型RAE2822的阻力.结果显示,2种方法具有相同的结果,表明尾迹面积分方法是有效的,与实验值吻合较好.尾迹面法中,尾迹面位置应处于离后缘相对弦长距离0.6~1.0之间.尾迹面积分方法中积分结果不依赖于物体的详细几何外形,可以预计对曲率变化大的三维复杂外形,该方法有更大的优势.      

Modelling and observing Jovian electron propagation times in the inner heliosphere

The propagation of Jovian electrons in interplanetary space was modelled by solving the relevant transport equation numerically through the use of stochastic differential equations. This approach allows us to calculate, for the first time, the propagation time of Jovian electrons from the Jovian magnetosphere to Earth. Using observed quiet-time increases of electron intensities at Earth, we also derive values for this quantity. Comparing the modelled and observed propagation times we can gauge the magnitude of the transport parameters sufficiently to place a limit on the 6 MeV Jovian electron flux reaching Earth. We also investigate how the modelled propagation time, and corresponding Jovian electron flux, varies with the well-known ∼13 month periodicity in the magnetic connectivity of Earth and Jupiter. The results show that the Jovian electron intensity varies by a factor of ∼10 during this cycle of magnetic connectivity.     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。