不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A)扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.     

一种结合壁面模型模拟湍流流动的间断界面浸入边界法

提出了一种具有鲁棒性的间断界面浸入边界法，用来模拟不可压湍流的物体绕流问题。开发了一种新的OpenFOAM求解器，该求解器结合了压力隐式分裂算法和k-ωSST湍流模型求解Navier-Stokes方程。为了减少对近壁面网格的要求，在实施边界条件时，采用壁面模型计算得到的壁面剪切力将速度和湍流变量修正联系起来。壁面切应力采用远离壁面的镜像点求得，镜像点上流场变量通过反距离线性插值由周围流场信息得到。使用该求解器模拟了平板绕流、NACA0012表面压强和切应力分布、Buice扩散管、圆柱绕流、方柱绕流5个湍流流动问题，本文数值模拟结果与文献结果有较好的一致性，验证了该方法模拟湍流流动的有效性和准确性。     

用三维N—S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松驰迭代来实现方程的求解，湍流的计算采用了标准的双方程模型，本文给出了钝体和类小汽车外形的计算结果，钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

可压缩粘性流动笛卡尔网格虚拟单元方法研究

以二维高雷诺数可压缩粘性流动问题为背景,提出了一种全新的笛卡尔网格虚拟单元方法。基于壁面函数基本假设,构造了壁面函数-虚拟单元方法(WF-GCM),用于定义湍流壁面边界条件。引入参考点的概念计算虚拟单元上的基本变量与湍流变量值,定义了"非贴体"笛卡尔网格下的湍流壁面边界条件,并通过壁面函数模型修正近壁面单元与界面单元。基于自适应笛卡尔网格体系,采用发展的具有二阶精度的格心格式有限体积求解器,数值模拟了跨音速RAE2822翼型绕流问题与超音速圆柱绕流问题,计算结果与实验值吻合良好,显示了WF-GCM对高雷诺数可压缩粘性问题是有效的。     

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

湍流附面层方程的有限差分解

本文用三点隐式有限差分法数值求解了轴对称流动的层流、转捩和湍流可压缩附面层方程,湍流附面层用两层概念处理,而每一层用适当的涡粘性模型代替Reynolds应力项,用常值湍流Prandtl数将湍流热流项与Reynolds应力关联起来。支配方程用迭代法求解。计算结果表明,采用本文改进的涡粘性模型和涡粘性分布的平滑办法,对文中所计算的几种情况,都给出了满意的结果。     

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

扰流柱形状对层板流阻特性的影响研究

应用FLUENT软件,采用隐式有限体积法求解雷诺时均N-S方程,对几种不同扰流柱形状的层板冷却结构的内部流动进行了数值模拟.湍流模型采用Realizable κ-ε双方程模型,近壁面采用壁面函数法处理,利用SIMPLE算法求解速度与压力的耦合.计算结果表明,层板内部流场结构十分复杂.同时,还对几种不同扰流柱形状的层板进行了流阻试验,采用沿程阻力关系式及流量关系式得出了其流阻特性.试验结果与计算得出的流阻特性符合较好.     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

μ理论在电液负载模拟器中的应用

 研究μ理论在电液负载模拟器中的应用，提出电液负载模拟器的鲁棒控制策略。电液负载模拟器是一个复杂的机-电-液复合系统，且是一个强耦合、时变受控对象。综合考虑参数变化、模型变动和干扰等不确定性的影响，利用μ综合控制理论设计电液负载模拟器的鲁棒力控制器，并通过μ分析使用鲁棒力控制器时系统的鲁棒稳定性和鲁棒性能。给出使用鲁棒力控制器和经典力控制器时的实验结果，实验结果表明所设计鲁棒力控制器的有效性和优越性。     

基于圆型限制性三体问题模型的行星卫星逃逸能量研究

 基于由飞行器、行星及其卫星组成圆型限制性三体问题模型，通过庞加莱映射的方法，研究了飞行器从行星卫星附近逃逸的问题。在Jacobi常数确定的前提下，通过逆向积分，飞行器从L1或L2点附近返回近月点，得到近月点速度出发速度。研究结果表明绕飞L1点和L2点逃逸行星卫星需要的最低能量是不同的，从月球表面逃逸所需的速度脉冲分别比开普勒算法节省46.5m/s和42.3m/s，且均小于Villac等人在Hill模型下得到38.9m/s，从而改进了Villac等人的相关工作，同时也给出了从太阳系主要行星卫星表面逃逸所需的最小能量。     

Model Development and Adaptive Imbalance Vibration Control of Magnetic Suspended System

对星载磁悬浮单框架控制力矩陀螺,建立其转子的平动与转动的动力学方程。所建模型反映了转子的动、静不平衡特性,以及转子运动与框架、星体运动的耦合关系。模型是主动磁悬浮控制与动框架效应分析的基础。对动、静不平衡量未知的转子,设计了自适应对中控制器。在分散PD、比例交叉反馈控制器与高低通滤波器的基础上,通过动、静不平衡量的辨识与补偿控制,使转子绕其惯量主轴旋转,衰减角动量交换执行机构高速转子因不平衡而对基座产生的干扰力,消除星上主要颤振源,达到降低卫星姿态抖动的目的。     

Love Waves in Layered Graded Composite Structures with Imperfectly Bonded Interface

Theoretical analysis and numerical calculations of Love wave propagation in layered graded composites with imperfectly bonded interface are described in this paper. On the basis of WKB method, the approximate analytic solutions for Love waves are obtained. By the interface shear spring model, the dispersion relations for Love waves in layered graded composite structures with rigid, slip, and imperfectly bonded interfaces are given, and the effects of the interface conditions on the phase velocities of Love waves in SiC/Al lay- ered graded composites are discussed. Numerical analysis shows that the phase velocity decreases when the defined flexibility parameter is greater. For the general imperfectly bonded interface, the phase velocity changes in the range of the velocities for the rigid and slip interface conditions.     

Experimental Investigation on Performance of Pulse Detonation Rocket Engine Model

The PDRE test model used in these experiments utilized kerosene as the fuel, oxygen as oxidizer, and nitrogen as purge gas. The solenoid valves were employed to control intermittent supplies of kerosene, oxygen and purge gas. PDRE test model was 50 mm in inner diameter by 1.2 m long. The DDT （deflagration to detonation transition） enhancement device Shchelkin spiral was used in the test model. The effects of detonation frequency on its time-averaged thrust and specific impulse were experimentally investigated. The obtained results showes that the time-averaged thrust of PDRE test model was approximately proportional to the detonation frequency. For the detonation frequency 20 Hz, the time-averaged thrust was around 107 N, and the specific impulse was around 125 s. The nozzle experiments were conducted using PDRE test model with three traditional nozzles. The experimental results obtained demonstrated that all of those nozzles could augment the thrust and specific impulse. Among those three nozzles, the convergent nozzle had the largest increased augmentation, which was approximately 18%, under the specific condition of the experiment.     

Compression Creep Behavior of High Volume Fraction of SiC Particles Reinforced Al Composite Fabricated by Pressureless Infiltration

The compression creep deformation of the high volume fraction of SiC particles reinforced Al-Mg-Si composite fabricated by pressure-less infiltration was investigated. The experimental results show that the creep stress exponents are very high at temperatures of 673 K, 723 K and 773 K, and if taking the threshold stress into account, the true stress exponent of minimum creep strain rate is still approximately 5, although the volume fraction of reinforcements is very high. The creep strain rate in the high volume fraction rein- forced aluminum alloy matrix composites is controlled by matrix lattice diffusion. It is found that the creep-strengthening effect of high volume fraction of silicon carbide particles is significant, although the particles do not form effective obstacles to dislocation motion.     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

Adaptive Control of Helicopter Ground Resonance with Magnetorheological Damper

The methodology for adaptive control of helicopter ground resonance with magnetorheological （MR） damper is presented. The adaptive inverse control method is used to control the output damping force of MR damper and the range of the damping force is given. Through the adaptive inverse control, the damping force of MR damper is fit to a desired damping force. With the background of applying MR damper to control of helicopter ground resonance, a model of loss force and an adaptive arithmetic for stabilization of the coupled rotor/fuselage system are presented. The simulation shows that the controller presented in this paper can stabilize the rotor/fuselage coupling system quickly and control the helicopter ground resonance effectively.     

The Full Flowpath Analysis of a Hypersonic Vehicle

对一种类X43-A吸气式高超飞行器全流道开展了M7一级的三维数值仿真研究,深入探索了进气道处于起动状态和不起动状态时全流道的冷流流场结构和和气动力特性,且部分结果与实验数据进行了对比。研究结果表明:(1)前体横截面上存在显著的展向压强梯度,使得经过预压缩的气流偏离了进气道进口,但同时也减少了进入内通道的边界层气流,提高了进口流场的品质;(2)后体喷流股的膨胀过程受到了周围外流的显著干扰,因而沿流动方向其截面形状不断发生变化,如在喷口附近为近似矩形,而在后体末端附近则演化为近似三角形;(3)当进气道处于不起动状态时,其全流道流动结构发生了显著变化,进气道的外部压缩波系往复振荡,尾喷管出口的喷流股也在不停的膨胀和收缩,具有极强的非定常特征;(4)当进气道处于不起动状态时,全机的气动力特性呈周期性变化,升阻比的变化幅度较大,在最大、最小值分别可达起动状态的2倍和1/4倍。另外,升力、阻力系数的变化曲线之间存在一定的相位差;(5)与实验结果的对照表明,所采用的数值仿真方法具有较高的精度。