含尘气体中的非定常激波

对含尘气体中非定常激波进行了实验研究和数值模拟。实验在垂直两相激波管中实现。气体与含尘气体之间分界平面的形成以及分界面处粉尘浓度的瞬时测定为实验研究提供了有效的手段。在双流体模型基础上，应用时间分裂法，将ＴＶＤ格式和ＭａｃＣｏｒｍａｃｋ格式分别应用于气相和颗粒相方程，模拟了激波与含尘气体相互作用现象，讨论了含尘气体物理参数对激波传播行为的影响。计算结果与实验现象符合较好。     

气动阻尼对高层建筑横风向风振响应的影响

引入横风力谱模型，分析了某高宽比为6 的方形截面高层建筑在不同地貌和风速条件下的横风向风振响应. 同时考虑正气动阻尼的影 响，获得该建筑在不同自振基频下的风振响应及气动阻尼影响的规律. 分析结果表 明，对于处于低粗糙度地貌并受较高风速作用的低频建筑物，当计算其横风向风振响应时，应适当考虑 正气动阻尼的影响，使计算结果更具真实性.     

顶部重力加砂颗粒速度场及浓度场的实验研究

对砂尘环境试验设备顶部重力加砂气固两相流的颗粒速度场及浓度场进行了研究．利用中国科学院寒区旱区环境与工程研究所室内风沙环境风洞，采用激光数字式粒子成像测速技术(DPIV)，得到顺风向距加砂口不同水平距离截面位置颗粒相速度及空间分布规律．表明在顺风向距加砂口水平距离约3m左右处的位置，颗粒相速度才接近风洞风速，且粒子基本均匀扩散开，从而证明了国外砂尘实验设备的不合理性、为我国自行研制大型砂尘环境试验设备，确定加砂／尘方法及试验段参数等提供了依据．     

气固两相流场的湍流颗粒浓度理论模型

本文进行了气固两相流动颗粒湍流扩散现象的理论分析，提出了颗粒湍流扩散系数和气流弥散效应二个颗粒湍流模化新概念，在此基础上建立了气固两相流场湍流颗粒浓度模型。理论模型包括离心力和其它外加力场作用下颗粒运动和浓度分布的计算方法。运用湍流颗粒浓度模型，对直管气固两相流动、受限射流气固两相流动和９０°弯管气固两相流动等三种流动做了数值模拟，计算获得颗粒速度、颗粒浓度等主要流动参数。讨论了湍流颗粒浓度模型的适用性。     

结构横风向非平稳随机振动的矩分析

将双方向脉动风速引入到气动力中，从而得到结构横风向运动的非平稳随机振动方程．运用随机微分理论和矩方程理论，得到结构截面横风向运动时的一阶矩和二阶矩的微分方程．通过数值方法，得到结构截面横风向随机振动的时程曲线．通过研究不同地貌和不同方向脉动风速对结构响应的影响，得出脉动风速及其非线性耦合项对结构横向振动具有不可忽视的影响的结论。     

激波点燃粉尘的数值模拟研究

建立了用于模拟入射激波后可燃粉尘颗粒点火的一维非定常两相化学反应流模型，该模型考虑了气固两相间的相互作用、粉尘颗粒的加速、加热和化学反应。粉尘颗粒着火前的化学反应用发生在颗粒外表面和内孔表面的非均相反应描述，颗粒内部的温度变化用一含有化学反应源项的非稳态热传导方程来描述，以颗粒外表面温度的突跃上升作为可燃粉尘颗粒点燃的着火条件。我们用该模型和ＰＳＩＣ方法，对由中等强度激波从纯气相传入煤粉－氧气混合物而引起的非定常两相流动现象，包括气固两相间的相互作用、粉尘颗粒的加速、加热以及点火过程进行了数值研究，计算了对应于不同载荷比、马赫数为４～５的入射激波后煤尘颗粒的点火延迟时间，分析了由于可燃粉尘颗粒的存在，入射激波及波后气固两相流动参数的变化规律。数值计算结果与实验数据符合较好。文中建立的模型和所用的基于ＰＳＩＣ算法的数值方法，用最自然的方式描述气固两相流动，即用连续流模型（欧拉方程）描述输运相（气相）的流动，用轨道颗粒模型（拉格朗日方程）描述分散相（颗粒相）的运动。用这种方法模拟含尘介质中激波后颗粒的点火是很有效的，它可以清楚地确定哪一个颗粒群最先着火，它的初始位置以及在整个点火延迟时间内     

基于3自由度的新月形覆冰输电线舞动稳定性研究

针对覆冰输电线舞动问题提出了一种基于非对称空气动力系数矩阵的临界风速计算方法．基于拟静态理论得到覆冰输电线的气动载荷，该气动载荷考虑了横向运动以及扭转运动对相对风攻角的影响，最后建立等效的3自由度覆冰输电线舞动模型．在初始风攻角处对气动载荷进行泰勒展开，得到非对称的线性空气动力系数矩阵．结合3自由度振动方程以及非对称空气动力系数矩阵，采用Rourh-Hurwitz准则计算覆冰输电线舞动发生的临界风速．通过风洞实验测得新月型覆冰单导线的空气动力系数，根据本文提出的理论分析了竖向振动频率、面外振动频率以及扭转振动频率对临界风速的影响，最后与DenHartog理论得到的临界风速进行了对比．本文研究成果对于指导覆冰输电线路防舞设计具有理论意义．     

高超声速火星进入环境中颗粒运动特性研究

火星大气中会发生不同规模的沙尘暴,大气中蕴含的尘埃颗粒会对高速进入的火星探测器表面造成侵蚀并导致壁面热流增加,给探测器的热防护系统设计带来巨大挑战.文章针对高超声速火星进入环境两相流动问题,基于Euler-Lagrange框架建立非平衡流场与颗粒的单向耦合计算方法,采用模态半径为0.35μm的火星大气颗粒分布模型,研究不同尺寸颗粒在流场中的运动轨迹,获得高温相变模型对颗粒运动的影响以及不同粒径颗粒的撞击能量分布.结果表明,颗粒在高温流场中运动会吸热融化甚至蒸发,高温相变模型导致的颗粒直径减小对小尺寸颗粒运动轨迹有较大影响;当前计算状态下,直径3μm以上的颗粒具有较大的Stokes数且颗粒半径在运动过程中基本保持不变,其运动轨迹受流场影响较小,该尺寸颗粒的撞击分数均达95%以上,是造成壁面撞击的主要颗粒尺寸;撞击能量分数结果表明,直径3～10μm之间的颗粒是撞击能量的主要来源,约占总撞击能量的80%.     

水击驻波场中分散相颗粒的运动分析

根据水击在管段内形成的驻波场现象,分析了流体内分散相颗粒受到的驻波作用力;运用李雅普诺夫稳定判据研究了颗粒积聚与分离的机理;考虑到颗粒运动方程的严重刚性而很难进行数值求解,采用相空间和非对称分析方法获得了分散相颗粒的运动轨迹近似解,并进行了实验验证。结果表明:水击驻波场中分散相颗粒的受力方程中惯性项对颗粒初始运动速率的影响不可忽略;在水击驻波波节的±λ/4范围内,分散相颗粒经过一定的时间会发生积聚,其运动速度呈对称分布,最大速度出现在3λ/8位置处;随着分散相颗粒粒径和密度等物性参数以及水击驻波的频率和连续相初始速度的增大,颗粒达到平衡位置的时间呈减小趋势,且连续相的初始速度对颗粒到达波节时间的影响显著。     

低浓度固液两相流的颗粒相动理学模型

用广义Fokker-Planck扩散模型描述液相湍动对颗粒的挟带作用，用修正的BGK模型描述粒间碰撞效应，建立了封闭的颗粒相PDF输运方程．运用Chapman-Enskog迭代法求得方程的二阶近似解，获得颗粒相脉动速度二阶矩和三阶矩闭合关系．模型与颗粒流模型相容，与液相湍流闭合模型是否相容依赖于扩散模型的具体形式，并据此比较了不同的涡一颗粒作用模型．模型与二维明渠流轻质沙和天然沙试验资料符合很好．表明细小粒径颗粒能够充分跟随水流运动；大粒径颗粒的相间平均速度差和壁面滑移速度明显，近壁区内的颗粒沿流向和垂向脉动强度都可能大于水流，并存在一定程度的颗粒碰撞效应．     

激波诱导含灰气体边界层中的粒子分布

本文研究当激波沿着一个固体表面等速地穿越含灰气体运动时所诱导的层流边界层特性。考虑了作用在气体边界层中球形粒子的 Saffman 升力,建议了一种计算近壁区中弥散相密度剖面的方法,并给出了数值计算结果。本文结果表明:在激波后方存在着一个弯曲的薄层区域,其中的粒子密度可以比其波前原始值增加许多倍。这种粒子聚集效应对于工业中粉尘爆炸等实际问题具有重要意义。     

The distribution of particles in a shock-induced boundary layer of a dusty gas over a solid surface

The laminar boundary layer behind a constant-speed shock wave moving through a dusty gas along a solid surface is studied. The Saffman lift force acting on a spherical particle in a gas boundary layer is taken into account. A method for calculating the density profile of dispersed phase near the wall is proposed and some numerical results are given. It is shown that behind the shock wave, there exists a curved thin layer where the density of particles is many times higher than the original one. This dust collection effect may be of essential importance to the problem of dust explosion in industry.     

The flow structure of dilute gas-particle suspensions behind a shock wave moving along a flat surface

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented.The Saffman lift force acting on a particle in a shear flow istaken into account.It is shown that particle migration across the boundary layer leads tointersections of particle trajectories.The corresponding modification of dusty gas model isproposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at aconstant speed are obtained by using the method of matched asymptotic expansions.Themethod of the calculation of particle phase parameters in Lagrangian coordinates isdescribed in detail.Some numerical results for the case of small particle concentration aregiven.     

Boundary layer flow and heat transfer of a dusty fluid flow over a stretching sheet with non-uniform heat source/sink

The present paper deals with the analysis of boundary layer flow and heat transfer of a dusty fluid over a stretching sheet with the effect of non-uniform heat source/sink. Here we consider two types of heating processes namely (i) prescribed surface temperature and (ii) prescribed surface heat flux. The momentum and thermal boundary layer equations of motion are solved numerically using Runge Kutta Fehlberg fourth–fifth order method (RKF45 Method). The effects of fluid particle interaction parameter, Eckert number, Prandtl number, Number of dust particle and non-uniform heat generation/absorption parameter on temperature distribution are analyzed and also the effect of wall temperature gradient function and wall temperature function are tabulated and discussed.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Effect of viscous dissipation and heat source on flow and heat transfer of dusty fluid over unsteady stretching sheet

This paper investigates the problem of hydrodynamic boundary layer flow and heat transfer of a dusty fluid over an unsteady stretching surface.The study considers the effects of frictional heating(viscous dissipation) and internal heat generation or absorption.The basic equations governing the flow and heat transfer are reduced to a set of non-linear ordinary differential equations by applying suitable similarity transformations.The transformed equations are numerically solved by the Runge-Kutta-Fehlberg-45 order method.An analysis is carried out for two different cases of heating processes,namely,variable wall temperature(VWT) and variable heat flux(VHF).The effects of various physical parameters such as the magnetic parameter,the fluid-particle interaction parameter,the unsteady parameter,the Prandtl number,the Eckert number,the number density of dust particles,and the heat source/sink parameter on velocity and temperature profiles are shown in several plots.The effects of the wall temperature gradient function and the wall temperature function are tabulated and discussed.     

侧风风场特征对高速列车交会的影响研究

为了准确预测侧风环境下列车交会时的气动性能,考虑均匀风和大气底层速度边界型侧风,通过计算流体力学方法对比研究高速列车交会压力波特性和气动力特性. 结果表明,两种风场下交会压力波有所不同,但幅值差异不明显;采用均匀风场评估列车侧风环境下交会会高估列车所受气动力. 本文对高速列车行驶安全性评估和复杂运行场景的复杂流场的认识具有参考价值.     

Unsteady Boundary Layer on a Blunt Body in a Hypersonic Flow of Non-Uniform Dusty Gas

The unsteady heat transfer at the stagnation point on a blunt body traveling at hypersonic velocity through a layer of nonuniform dusty gas with low-inertia particles (not deposited on the body surface) is investigated. Using the matched asymptotic expansion method, the equations of the two-phase unsteady boundary layer near the symmetry axis of the body are derived with account for the polydispersity of the particles. The structure of the unsteady boundary layer and the variation of the friction and heat transfer coefficients at the stagnation point are studied numerically. Layered nonuniformities of the particle concentration and size are considered, the limits of variation of the thermal and mechanicals loads are found, and the effect of the dust polydispersity on the heat transfer is investigated.     

Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

A nonlinear model for aerodynamic configuration of wake behind horizontal-axis wind turbine

Determination of the aerodynamic configuration of wake is the key to analysis and evaluation of the rotor aerodynamic characteristics of a horizontal-axis wind turbine.According to the aerodynamic configuration, the real magnitude and direction of the onflow velocity at the rotor blade can be determined, and subsequently, the aerodynamic force on the rotor can be determined. The commonly employed wake aerodynamic models are of the cylindrical form instead of the actual expanding one. This is because the influence of the radial component of the induced velocity on the wake configuration is neglected. Therefore, this model should be called a "linear model". Using this model means that the induced velocities at the rotor blades and aerodynamic loads on them would be inexact. An approximately accurate approach is proposed in this paper to determine the so-called "nonlinear" wake aerodynamic configuration by means of the potential theory,where the influence of all three coordinate components of the induced velocity on wake aerodynamic configuration is taken into account to obtain a kind of expanding wake that approximately looks like an actual one. First, the rotor aerodynamic model composed of axial(central), bound, and trailing vortexes is established with the help of the finite aspect wing theory. Then, the Biot-Savart formula for the potential flow theory is used to derive a set of integral equations to evaluate the three components of the induced velocity at any point within the wake. The numerical solution to the integral equations is found,and the loci of all elementary trailing vortex filaments behind the rotor are determined thereafter. Finally, to formulate an actual wind turbine rotor, using the nonlinear wake model, the induced velocity everywhere in the wake, especially that at the rotor blade,is obtained in the case of various tip speed ratios and compared with the wake boundary in a neutral atmospheric boundary layer. Hereby, some useful and referential conclusions are offered for the aerodynamic computation and design of the rotor of the horizontal-axis wind turbine.