两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

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

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

区域分裂法在湍流流场数值解中的研究

应用区域分裂算法(简称DDM)对湍流流场进行数值研究。对湍流流动控制方程建立一种通用的DDM-1计算格式S_{I,ω,给出重叠区域内变量处理方法。DDM-1的SI,ω计算格式的适用性由突扩圆管分离流流场数值计算所验证,数值计算结果与实验数据比较表明:采用SI,ω计算格式所预测的再附点位置及平均轴向速度、湍动能强度比在整个计算区域内采用单一模型的方法能够更好地与实验值吻合。}      

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

半模试验翼身组合体垫块的数值模拟分析

采用数值模拟方法研究存在风洞侧壁边界层影响时垫块厚度和修型对翼身组合体气动特性的影响.本文采用Jameson中心有限体积、多步Runge-Kutta时间推进格式求解三维可压、雷诺平均Navier-Stokes方程,引入Baldwin-Lomax湍流模型,分析了垫块厚度和修型对DLR-F4翼身组合体半模型气动力和局部流场的影响.通过本文计算得出的结论对于跨声速半模型风洞试验垫块优化设计具有重要的指导意义,能更好地提高风洞试验数据的精准度和试验技术水平.     

绕翼型低雷诺数流动的数值分析研究

分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非定常流动,并对两方程SSTk-ω湍流模型,代数B-L模型和层流的计算结果作了比较。N-S方程和两方程湍流模型的控制方程非耦合求解,时间推进均采用近似因子方法(AF);空间离散无粘项采用ROE格式,粘性项用中心差分方法计算;计算非定常流场时,采用伪时间子迭代(τt-s)方法保证二阶时间精度,湍流模型计算时都是在固定点转捩。最后分析了转捩对低雷诺数流动的影响、分离泡的存在引起的流动不稳定性和周期性的脱出涡,数值计算给出的时间平均升力系数、阻力系数和压力分布与实验结果比较吻合。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

基于非结构/混合网格的脱体涡模拟算法

为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与脱体涡模拟(DES)方法相匹配的二阶混合耗散自适应格式,根据流场信息自动调节格式耗散,并分别基于Spalart-Allmaras一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,发展了基于非结构/混合网格的DES方法。采用该方法计算了雷诺数为3 900的圆柱绕流和NACA 0021翼型60°大迎角分离流两个典型算例,通过与试验数据以及其他数值结果的对比验证了该方法的可行性。同时开展了不同数值格式、湍流模型的对比分析,研究结果表明:采用混合格式的DES算法能够解析更小尺度的湍流涡结构、计算数据更接近试验值;本文的DES类算法受其基准湍流模型影响较小。     

基于模型的飞行仿真系统设计

采用基于模型的仿真软件设计思想,研究了飞行仿真系统的设计方法,针对具体的任务和功能 要求,建立了仿真系统的结构体系。研究了仿真软件框架的设计方法和结构层次的关系,并对各个层次功能 要求进行了介绍。对飞行仿真软件系统的模块设计进行了研究,确立了模型的描述形式和接口设计规范。     

风洞试验模型技术新发展

风洞模型试验是航空航天飞行器研制的重要环节之一。试验模型的设计制造关系到风洞试验的数据质量、效率、周期和成本。本文归纳了近年来国外风洞模型技术的最新发展,分析了快速成型技术在风洞试验模型制造中的发展和应用;阐述了欧、美遥控风洞模型技术的发展理念、关键技术和应用研究;概述了风洞试验模型采用的新材料、抑振和变形测量技术。     

基于混沌理论的滚动轴承振动信号融合模型预测

提出融合算法模型,在混沌理论的基础上对滚动轴承振动信号进行预测。基于相图法、最大Lyapunov指数法和关联维数法对滚动轴承振动信号进行混沌判别,证明其混沌性。以预测值和真值间差值范数最小为目标导向优化出Kriging模型、最小二乘支持向量机(LSSVM)模型和极端学习机(ELM)模型的权重,加权法构建融合算法模型。相空间重构法构建滚动轴承振动信号预测的训练样本,并对融合模型、Kriging模型、LSSVM模型和ELM模型进行训练,训练好的模型用于振动轴承振动信号混沌预测。以案例1和案例2共两个实验的滚动轴承振动信号为对象进行验证,两案例的最大Lyapunov指数大于0,从而判断这两个案例的轴承振动信号呈现混沌特性。另外,从方均误差、方均根误差和平均绝对误差指标来评价,融合算法模型的指标值均小于单一模型算法,融合算法模型的预测精度优于单一模型算法。     

集群多机器人系统建模研究的发展与展望

凭借高效、鲁棒、应用广泛的特性,集群多机器人系统已经成为当今研究的热点课题之一,具有重要的实用价值。首先简述了自顶而下和自底而上两种多机器人系统研究思路的当前研究概况。然后从拟生物集群系统模型引入,进而引出一致性系统模型,针对低阶、高阶、异质、时延等一致性模型进行分析总结,单独阐述了多智能体强化学习系统模型的研究情况,并分别讨论了三种集群多机器人系统自组织建模方法的研究现状与各自存在的问题,总结与分析了集群多机器人系统运动的发生机理。最后,分析了现有集群多机器人系统模型尚待解决的关键问题和面临的挑战,并对其未来发展进行了展望。     

基于滞回碰摩力模型的转子系统碰摩响应研究

转/静子碰摩是影响旋转机械的稳定性与安全性的重要因素之一,而碰摩力的表征则是预测转子系统动力学响应的关键问题。讨论了已有转/静子碰摩力模型的适用范围,重点阐述滞回碰摩力模型的物理意义,并与几种经典的碰摩力模型进行比较;以Jeffcott转子作为典型例子,基于滞回碰摩力模型给出系统的碰摩响应,并与采用线性碰摩力模型所得到的碰摩响应进行比较;利用数值仿真的结果分析静子刚度、转子偏心距等因素对采用滞回碰摩力模型的转子碰摩响应的影响。数值仿真的结果表明:采用滞回碰摩力模型得出的失稳速度(1226 rad/s)低于采用线性碰摩力模型得出的失稳速度(1476 rad/s),且系统的动力学特性更为复杂,因此,可以更真实地反应碰摩的物理实质。     

航空发动机参数不确定模型辨识方法

针对航空发动机线性模型在工况点附近使用范围较小的问题,提出一种航空发动机不确定性模型辨识方法。该方法使用非线性规划处理航空发动机模型辨识问题,求解过程考虑线性模型中矩阵参数不确定性的影响,以期获得一个具有适用范围大、形式简单的航空发动机模型。使用该方法对DGEN380发动机在某一稳态点进行辨识,定义实际工况偏差参数和模型最大偏差参数分析DGEN380发动机参数不确定模型的误差范围。使用实验数据与参数不确定模型仿真结果进行比较,结果表明油门杆角度变化小于22%时,参数不确定模型与实际发动机状态的偏差量较小,能够在1%误差范围内模拟实际发动机状态。     

飞机几何外形三维数模构建的一种新算法

提出了一种飞机几何外形三维数字模型构建的新方法。该方法从飞机外形的工业设计效果图和飞机外形图片入手，利用图像处理技术，进行图像三维外形恢复，利用CAD三维平台进行数字化主模构建，同时对构建过程中的关键技术进行了分析研究，并提出了可行性的解决方案。     

飞行控制模型管理系统

模型管理是决策支持系统的核心，将模型管理应用于飞行控制系统是现代试飞管理的一个重要手段，也是一个新的探索领域。本文从模型及模型库管理系统的概念出发，分析传统的及近年来进入市场的主流数据库技术，研究对飞控系统进行模型管理的可行性。首先提出对飞控模型利用对象关系数据库进行管理的思想并将模型对象化；其次是模型管理概念的提出及模型库的数据模型，模型库组织和存储；最后是飞控系统模型管理的实现。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。