多重网格法求解离心压气机内部流场

本文提出了一种新的多重网格法,并将之与显式时间推进法、有限体积差分格式、粘性体积力法配合使用求解离心压气机内部流场。数值计算表明,作者多重网格法的计算结果与试验结果吻合比较好,同时该多重网格法的收敛速度比不用多重网格法加快了将近两倍。     

六阶精度的群速度直接控制紧致格式及其应用

差分数值解产生非物理振荡的直接原因是在于非均一的波群的群速度,因此本文利用直接群速度控制的方法重新构造具有六阶精度的紧致型差分格式,以达到改善激波数值解的目的。对所构造的格式的精度及数值解随控制参数变化的行为进行了分析。最后文中给出了一些典型算例,证明了文中所构造的格式具有方法简便,物理含义清楚,精度高,网格基架点小,和捕捉激波能力较强的优点。      

生成非结构网络的改进阵面推进法

通过对非结构化网格生成的阵面推进方法进行分析,对该算法进行了优化设计,采用笛卡尔直角坐标网格来代替背景非结构化网格,其上网格尺度参数通过控制体有限元法求解Possion方程的Dirchlet问题得到,提高了网格生成的效率和质量。同时利用Laplacian滤波算子和对角交换方法来改善三角形网格的质量,尤其是对角交换技术对消除畸形三角形单元效果明显。对固壁边界采用三次样条拟合后重新离散的方式形成初始阵面,改进了网格的贴体性,提高了网格的光滑性。在编制通用程序时利用链表结构和派生类型的数据变量,节省了计算内存,提高了网格生成效率。     

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

Comparison of autoignition characteristics of n-heptane, methylcyclohexane and toluene

The autoignition characteristics of three C7 hydrocarbon fuels,n-heptane,methylcyclohexane and toluene,were comparatively investigated.Ignitions were performed behind the reflected shock waves in a shock tube.The ignition delay times of these fuels were measured at the same igintion conditions with constant fuel mole fraction of 1.0%,equivalence ratio of 1.0,ignition pressure of 1.0×105 Pa(one more 2.0×105 Pa for n-heptane)and temperatures of 1 166-1 662 K.The correlation formula of ignition delay dependence of three fuels on igintion conditions was deduced separately.Results show that the ignition delay time of n-heptane is the shortest while that of toluene is the longest at the same ignition conditions.The ignition delay time of methylcyclohexane is most sensitive to the temperature while that of n-heptane is the least.The comparison of current ignition delay times with the predictions of available chemical kinetic reaction mechanisms has been presented to validate the reliability of mechanisms.The important chemical reactions during the ignition process have been obtained from the sensitivity analysis.     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

基于灵敏度的可靠性优化解耦方法

针对航空结构可靠性优化设计问题,提出了一种基于灵敏度的可靠性优化(RBO)解耦方法。首先将高效求解的可靠性灵敏度用于失效概率函数(FPF)的快速构建,其优点在于仅需要一次可靠性分析即可得到失效概率函数的局部近似,克服了常规求解方法中需要多次可靠性分析的缺点;然后将得出的FPF近似代入RBO模型中,将原问题解耦成确定性优化子问题,可用常规优化方法来求解。另外,采用了序列近似优化策略来保证可靠性优化解的正确性。文中给出了复合材料梁和机翼三盒段结构的优化求解算例来说明本文方法的可行性和正确性。     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

直齿锥齿轮低安装误差敏感性设计与实验验证

未修形直齿锥齿轮啮合时为线接触,为减小两齿面啮合对安装误差敏感性,对主动轮齿面采用鼓形修形,改变刨刀运动轨迹进行齿向抛物线修形、改变瞬时滚比进行齿廓修形,从而实现两齿面点接触啮合.为进一步降低安装误差的敏感性,以齿廓、齿向修形系数为优化变量,减小接触迹线上啮合点的差曲面高斯曲率波动,同时增大参考点的差曲面高斯曲率值,降低安装误差敏感性;为避免齿面接触应力过大,控制参考点的差曲面主曲率,保证瞬时接触椭圆的长度不小于齿宽的1/3.算例分析及加工和滚检实验显示:经过优化设计后的修形齿面安装误差敏感性较低,总轴向错位量和总轴向分离量分别达到法向模数的30%,齿面印痕仍具有较好的稳定性.