后掠角对后掠机翼边界层稳定性及转捩的影响

后掠机翼边界层流动稳定性及转捩对翼型的设计及优化有着重要的参考价值,而机翼后掠角是引起后掠机翼边界层横流失稳的关键参数之一.以NACA0012翼型为研究对象,通过求解三维可压缩Navier-Stokes方程计算了展向无限长后掠机翼的基本流场;通过求解Orr-Sommerfeld方程得到了扰动Tollmien-Schishting波演化的中性曲线及幅值曲线,研究了后掠角对后掠机翼边界层流动稳定性的影响;最后采用eN方法进行了转捩预测.研究发现,随后掠角的增大,横流强度和扰动幅值放大指数n均先增加后减小,且后掠角在40°~50°之间横流强度达到最大值.当后掠角在50°左右时,用转捩预测eN方法计算的幅值增长指数N值最大,导致转捩发生所需的初始扰动幅值最小,转捩最易发生.      

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

基于流动控制的水下航行体流体动力特性分析

利用CFX对水下航行体流体动力特性进行了数值模拟研究.求解了混合介质的剪切压力输运湍流模型控制方程、雷诺-平均奈维尔-斯托克斯方程以及各相间的质量输运方程,采用三维数值模拟方法对比分析了不同流动控制方案的水下航行体空泡形态特征、表面压力分布和阻力系数变化特性,讨论了不同参数对减阻效果的影响.结果表明,对于有横流作用的水下垂直发射航行体,多相流动控制不仅可以降低通气空泡的不对称性和航行体阻力,同时可以均匀迎流面和背流面压力,从而实现航行体多相流空泡形态及表面压力特性等的控制.此外,通气口的位置对减阻效果具有显著影响.      

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

民用通风机气动设计中弯掠长叶片的数值模拟

采用数值模拟的方法,分别对弯掠叶片造型和直叶片造型两种情况下高性能民用通风机叶片内三维紊流流场进行了计算;计算程序中采用有限体积法,在任意非正交曲线坐标系下,用SIMPLE算法求解三维N-S方程和K-ε湍流方程.计算结果表明,弯掠叶片造型的民用通风机的流动参数和总性能参数更接近设计指标,而且采用弯掠叶片可以更有效地控制气流的流动,对减少二次流损失起了很大作用.      

飞行器表面三维流场与固壁温度场的耦合分析

分析了飞行器表面温度场形成的物理机制,建立了飞行器蒙皮温度场耦合计算的理论模型,较为完整的考虑了各种热源对飞行器蒙皮温度的作用.采用有限体积法和重整化群RNG(Renormalization Group)k-ε湍流模型对流场求解N-S方程,蒙皮固壁中采用热传导模型,采用离散坐标法DOM(Discrete Ordinates Method)计算蒙皮固体壁面对外界大气的辐射.通过数值计算研究了飞行高度、马赫数和内热源对蒙皮温度分布的影响,以及高度及马赫数对飞行器表面气动加热的作用.结果表明,该方法适合考虑飞行器外部流场、固壁等多种热源作用的温度场数值模拟.     

Rayleigh问题的Monte Carlo直接模拟

利用流场直接模拟的Ｍｏｎｔｅ Ｃａｒｌｏ方法数值模拟了Ｒａｙｌｅｉｇｈ问题的流场，获得了相应的结果，结果与解析解及经验解吻合较好，ＤＳＭＣ方法避免了求解复杂的控制方程，在流动比较复杂的情形下，有很大的优越性，是研究稀薄气体流动的一种好方法。     

脉冲射流强化喷流混合的数值模拟

采用高振幅低流量射流射入亚音剪切层以强化喷流混合,对利用脉冲射流强化喷流混合的流场进行了数值模拟.利用有限体积法和重整化群(RNG)k-ε湍流模型求解N-S方程.分别进行了低亚音马赫数及中等亚音马赫数喷流混合的数值模拟.数值模拟结果显示,无论是对于低亚音速热喷流还是中等亚音速热喷流,脉冲激励引起喷流在与激励源平行的平面迅速扩展,在与之垂直的平面收缩,当脉冲射流激励的斯德鲁哈尔数在0.2时,可以得到较好混合,而且少量的射流流量(占主流流量的1%)可以极大地增加速度及温度的衰减率,随着流量的增加(占主流流量的3%),喷流混合得到进一步的加强.     

浮空气囊定常及非定常气动特性计算

对于0~20 km范围内不同外形的浮空气囊,采用虚拟压缩方法求解不可压缩非定常N-S(Navier-Stokes)方程,数值模拟了5个高度下气囊的定常和非定常流场.定常状态时,随高度及气囊外形的变化,模拟了气囊气动力及背风区复杂流动分离.结果表明,当气囊厚度率逐渐增大时,其背风区存在复杂的空间流动分离.运用双时间步方法及动网格技术,成功实现了气囊平移或俯仰运动状态下非定常流场的数值模拟.考察了不同的平移方向及速度对气囊非定常气动特性的影响,得到了气囊平移过程中气动力及流场的变化情况.对于气囊的俯仰运动,获得了不同俯仰速度及不同最大俯仰角时气囊的气动力及流场.计算结果表明气囊的非定常运动对其流场和气动力影响显著.     

柔性杆子杆模型参数的遗传算法优化求解

研究了应用实数型的遗传算法解决柔性机器人子杆法建模中,子杆模型参数的优化求解问题.子杆模型参数形成高度非线性的方程组,它的求解是子杆法能否实现的关键.普通的牛顿迭代法及其变形,要求初始向量足够靠近解,一些所谓大范围收敛的算法又有方程组的Jacobian矩阵非奇异等要求,在此都难以应用.本文基于遗传算法的解法可求得模型参数的最优解,实例仿真显示了这一方法的可行性和高效率.      

卫星瞬态温度场的计算分析

本文讨论了地球同步通信卫星热设计中的几个典型瞬态溫度问题,采用Stiff稳定的自动积分法,对250阶的非线性微分方程组进行数值求解,在校正迭代计算中配合以稀疏矩阵的技术,获得了良好的效果。     

Newton-Kantorovich/Radau pseudospectral solution to rigid spacecraft time-optimal three-axis reorientation

The Newton-Kantorovich/Radau pseudospectral method is implemented in a homotopy scheme to solve the rigid spacecraft time-optimal three-axis reorientation problem. First, the reorientation differential equations associated with the current homotopic variable are linearized about the nominal (previous homotopic variable) solution. Then, the Radau pseudospectral method is used to solve the resulting linearized differential equations. By successively updating the nominal solution with the Newton-Kantorovich approach, the linearized solving procedure is iterated to achieve high-order precision. Because no numerical integration is needed and only linear algebraic equations are solved, the proposed method is computationally efficient. Finally, performance of the proposed method is illustrated through numerical examples.     

差分格式的优化组合

用偏微分方程组解实际的三维紊流流动问题时,如何选择各方程中对流项和扩散项的差分格式将对计算的稳定性、准确性与经济性影响很大.对一个实际的流动问题,计算机模拟验证了迎风差分、中心差分、混和格式、PLDS(Power Law Differencing Scheme)、QUICK(二次迎风插值)和OPTIMAL(Oscillation Preventing Technical algorithm based on the second-order Interpolation Method for Advection Linked convection)的各种组合对计算的稳定性、准确性与经济性的影响.对中心差分应用在雷诺时均方程中出现的数值解的振荡问题提出了在局部网格使用迎风差分的具体对策.对QUICK应用在κ,ε输运方程上出现的数值解的发散问题提出了QUICK和迎风差分格式交替使用的具体方法.      

矩阵带宽缩减技术在隐式间断有限元中的应用

为了数值求解二维Euler方程，以间断有限元方法作为空间离散、向后差分公式（BDF）作为时间离散。针对采用牛顿法求解源于隐式时间积分的非线性方程组，构造了相应的Jacobi矩阵，其具有阶数高、稀疏性强、数值非对称的特点。在每个时间步内，选择带预处理的广义极小残量（GMRES）方法求解线性方程组，预处理矩阵由不完全LU分解（ILU）方法构造。将矩阵带宽缩减技术应用于上述求解过程，无需额外的存储空间，就缩小了预处理矩阵与系数矩阵的差距，从而加快了GMRES方法的收敛、增大了可用的时间步长。通过求解典型的空气动力学问题，检验了该应用的有效性。      

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式,具有十分优异的分辨率特性。但在求解多维方程组时,HCW-R格式需要求解块状三对角方程组,因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P,构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式,其可沿迎风方向、由边界处向内推进求解,避免了处理三对角或块状三对角方程组,从而其计算代价与显式格式无异。虽然就分辨率而言,HCW-E稍逊于HCW-R,但前者的计算效率要显著高于后者。因此,当花费相同的计算代价,HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

带多个充液储箱航天器的耦合动力学建模方法

研究了充液航天器储箱内液体晃动与全星姿态运动的耦合动力学建模方法．分别基于单摆等效力学模型和质心面等效力学模型建立了全星姿态耦合动力学方程,并分别针对2个储箱串联布局、3个储箱并联布局的充液航天器在不同工况下的液体晃动问题进行了数值仿真研究;同时比较分析了两种储箱布局方案下液体晃动对航天器姿态运动的影响,为工程应用提供参考．     

对连续流式电泳简化数学模型的研究

用有限分析法对连续流式电泳的电流体力学过程进行了三维数值模拟。数值结果表明,Ｐｏｌｅｚｈａｅｖ（［１,２］）等人提出的电泳简化数学模型在处理重力条件下有回流问题时,所得结果与三维准确方程的计算结果严重不符。因此该模型不能用于重力有回流问题。为了克服上述缺陷,在解析分析的基础上,文章提出了一个新的电泳问题简化数学模型。数值计算显示,该模型所得结果与三维准确方程的计算结果十分吻合。     

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

叉流式热交换器动态特性的数值计算

建立了叉流式热交换器在热流体和冷流体都随时间变化的条件下的动态二维数学模型。使用一阶精度显式格式建立了热交换器的有限差分模型。该模型能够满足任意变化的温度及流量入口条件。经过实例计算表明,所得的理论热交换器出口平均温度变化曲线与相应的空中试飞测试曲线有较好的跟随性、吻合程度良好。