激波诱导两相流中影响阻力系数的特性参数研究

基于双流体模型 ,利用Euler Lagrange组合方法 ,对激波诱导的气固两相流场进行了数值计算 ,系统研究了影响颗粒群阻力系数的几个重要特性参数。结果表明 :目前采用激波管技术研究非定常条件下颗粒群阻力系数时界定这些因素的影响程度是必要的。     

Motion of tracer particles in supersonic flows

 Factors that may act on particle motion in high-speed flow are investigated. The classical expressions of drag coefficient C _D for a sphere are reviewed. Then, a drag expression is proposed, extending Cunningham’s method to higher velocities and Knudsen numbers. This law, valid from continuum to free molecule conditions, for Re≲200 and M≲1 (where Re and M are, respectively, the Reynolds and Mach numbers based on relative velocity), is used to compare calculated and experimental values of the drag coefficient, as well as the particle velocities across an oblique shock wave. Calculated results are found to be in agreement with experiments. Received: 3 June 1997/Accepted: 16 August 1998     

气体动理学统一算法中相容性条件不满足引起的数值误差及其影响研究

求解玻尔兹曼(Boltzmann) 模型方程的气体动理学统一算法(unified gas kinetic scheme,UGKS) 是为模拟存在显著稀薄气体效应流动而建立的. 在该方法中,如果速度空间离散采用传统的离散速度坐标法(discreteordinate method,DOM),将会导致相容性条件得不到严格满足,从而引入数值误差. 本文从理论分析及数值试验两方面说明了该数值误差,正比于来流马赫数,反比于来流努森数. 引入了守恒型的离散速度坐标法(conservativediscrete ordinate method,CDOM),在离散层面上确保了相容性条件得到严格满足. 圆柱绕流计算结果表明,来流马赫数较高、努森数较小时,相容性条件满足与否对计算结果影响较大,采用CDOM 可以在较稀的速度空间网格上得到网格无关解,缩减计算量最大可达2/3.      

New formula for drag coefficient of cylindrical particles

The drag force on a cylindrical particle is calculated using lattice Boltzmann method.The results show that the drag coefficient of a particle with different orientation angles decreases with increasing Reynolds number.When the principal axis of the particle is parallel to flow,the drag coefficient is much larger than that of others and decreases fastest with increasing Reynolds number,which becomes more obvious with increasing particle aspect ratio.When the principal axis of the particle is inclined to flo...     

Aerodynamic coefficients of a spinning sphere in a rarefied-gas flow

A three-dimensional rarefied-gas flow past a spinning sphere in the transitional and near-continuum flow regimes is studied numerically. The rarefaction and compressibility effects on the lateral (Magnus) force and the aerodynamic torque exerted on the sphere are investigated for the first time. The coefficients of the drag force, the Magnus force, and the aerodynamic torque are found for Mach numbers ranging from 0.1 to 2 and Knudsen numbers ranging from 0.05 to 20. In the transitional regime, at a certain Knudsen number depending on the Mach number the Magnus force direction changes. This change is attributable to the increase in the role of normal stresses and the decrease in the contribution of the shear stresses to the Magnus force with decrease in the Knudsen number. A semi-empirical formula for the calculation of the Magnus force coefficient in the transitional flow regime is proposed.     

Shock tunnel study of spiked aerodynamic bodies flying at hypersonic Mach numbers

A three-component accelerometer balance system is used to study the drag reduction effect of an aerodisc on large angle blunt cones flying at hypersonic Mach numbers. Measurements in a hypersonic shock tunnel at a freestream Mach number of 5.75 indicate more than 50% reduction in the drag coefficient for a 120° apex angle blunt cone with a forward facing aerospike having a flat faced aerodisc at moderate angles of attack. Enhancement of drag has been observed for higher angles of attack due to the impingement of the flow separation shock on the windward side of the cone. The flowfields around the large angle blunt cone with aerospike assembly flying at hypersonic Mach numbers are also simulated numerically using a commercial CFD code. The pressure and density levels on the model surface, which is under the aerodynamic shadow of the flat disc tipped spike, are found very low and a drag reduction of 64.34% has been deduced numerically.     

Numerical simulation of sedimentation of rectangular particle in Newtonian fluid

The sedimentation of a rectangular particle falling in a two-dimensional channel filled with Newtonian fluid was simulated with finite element arbitrary Lagrangian–Eulerian domain method. The numerical procedure was validated by comparison of the simulation results with existing numerical work. Moreover, good agreement was obtained between the simulation results and experimental measurements performed in the current study. The equilibrium position, stable orientation and drag coefficient of a rectangular particle for different particle Reynolds numbers (Re_p) were studied. The results show that there is a critical particle Reynolds number for the preferred orientation of a rectangular particle falling in a Newtonian fluid. When Re_p is smaller than the critical value, the particle falls with its long side parallel to gravity; otherwise the particle falls with its long side perpendicular to gravity. The critical particle Reynolds number is a decreasing function of the blockage ratio and aspect ratio. The distributions of pressure and shear stress on rectangular particle surface were analyzed. Moreover, the drag coefficient of the rectangular particle decreases as Re_p or the blockage ratio increases; however, it appears to be independent of aspect ratio.     

乘波体前体/进气道一体化设计研究

为了研究乘波体几何外形参数和飞行参数对前体/进气道一体化设计的影响,采用理论分析和数值模拟相结合的方法,以马赫数Ma=6和攻角α=0为设计状态、进气道总压恢复系数和前体阻力系数为目标函数,对乘波体前体/进气道进行了优化设计,并在此基础上研究了攻角、马赫数、前缘半径、前体宽度对气动参数的影响。结果表明:该乘波体前体/进气道构型具有良好的攻角特性,总压恢复系数比基准构型提高17.79%,阻力系数比基准构型降低78.5%,符合高超声速飞行器高升力、低阻力的要求,且非常适合小攻角高超声速巡航飞行;为了得到较高升阻比的前体,在前缘半径R≤2mm的范围内进行流场反设计时,可以将设计马赫数的取值比预期低一些。     

Shock tube spherical particle accelerating study for drag coefficient determination

An original particle accelerating technique has been developed for a shock tube. The trajectories of calibrated spherical particles and in diameter have been measured by the multiple exposure shadowgraph technique coupled with a high speed drum camera. Both particle velocity and acceleration, deduced from the experimental trajectories, allow the determination of the drag coefficients for different, subsonic and supersonic, flow regimes for the particle Reynolds numbers from to and the particle Mach numbers from 0.6 to 1.2. The drag coefficient values have been compared with different correlations found in the literature. Received 8 April 2002/ Accepted 17 June 2002 Published online 19 December 2002 Correspondence to: L. Houas (e-mail: Lazhar.Houas@polytech.univ-mrs.fr)     

Drag forces of interacting spheres in power-law fluids

Drag forces of interacting particles suspended in power-law fluid flows were investigated in this study. The drag forces of interacting spheres were directly measured by using a micro-force measuring system. The tested particles include a pair of interacting spheres in tandem and individual spheres in a cubic matrix of multi-sphere in flows with the particle Reynolds number from 0.7 to 23. Aqueous carboxymethycellulose (CMC) solutions and glycerin solutions were used as the fluid media in which the interacting spheres were suspended. The range of power-law index varied from 0.6 to 1.0. In conjunction to the drag force measurements, the flow patterns and velocity fields of power-law flows over a pair of interacting spheres were also obtained from the laser assisted flow visualization and numerical simulation.

Both experimental and computational results suggest that, while the drag force of an isolated sphere depends on the power-index, the drag coefficient ratio of an interacting sphere is independent from the power-law index but strongly depends on the separation distance and the particle Reynolds number. Our study also shows that the drag force of a particle in an assemblage is strongly positions dependent, with a maximum difference up to 38%.     

Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

Numerical simulation of the accumulation of heavy particles in a circular bounded vortex flow

In present work, an Eulerian–Lagrangian CFD model based on the discrete element method (DEM) and immersed boundary method (IBM) has been developed, validated and used to investigate the accumulation of heavy particles in a circular bounded viscous vortex flow. The inter-particle and particle-wall collisions are resolved by a hard-sphere model. Effects of one-way and two-way coupling, Reynolds number, and particle diameter are systematically explored. Results show that, in case of one-way coupling, the majority of particles will spiral into an accumulation point located near the stagnation point of the flow field. The accumulation point represents a stable equilibrium point as the drag created by the flow field balances the destabilizing centrifugal force on the particle. However, in case of two-way coupling, there does not exist a stable accumulation point due to the strong interaction between the particles and fluid dynamics. Instead most particles are expelled from the circular domain and accumulate on the confining wall. The percentage of accumulated particles on the wall increases with increasing Reynolds number and particle diameter. Moreover, influence of three well-known drag models is also studied and they give consistent results on the particle accumulation behavior, although small quantitative differences can still be discerned.     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.     

基于重叠网格与结构网格的圆柱绕流数值模拟

为研究重叠网格与结构网格在圆柱绕流数值模拟中的区别,以二维圆柱为例,利用有限元分析软件ANSYS 19.2中的DM与Mesh建立模型并划分重叠网格,利用ANSYS 19.2中的ICEM建立模型并划分结构网格。采用FLUENT 19.2中laminar模型模拟分析系统中的平均升力系数、平均阻力系数、斯特劳哈尔数St等流体动力特性。通过改变流体流速得到两种不同网格下各6组雷诺数Re,这6组雷诺数在60~160之间。结果表明:结构网格与重叠网格的St都随着Re的增加而增加,但相同雷诺数下重叠网格对应的St数值更大,St的增长速度更快;重叠网格与结构网格的平均升力系数与平均阻力系数随着Re的增加趋于稳定的速度都加快,但结构网格的平均升力系数与平均阻力系数趋于稳定的速度更快,且两种网格的平均升力系数与平均阻力系数趋于稳定速度的差距逐渐缩小,当Re=160时,两种网格的平均升力系数与平均阻力系数趋于稳定的速度几乎相同;当雷诺数在60~160之间时,采用重叠网格计算出来的斯特劳哈尔数比结构网格更加接近理论值;从升力功率谱密度分布曲线中可以看出,随着雷诺数的增加,两种网格下的频率逐渐变大,并且相同雷诺数下重叠网格的频率比结构网格大。     

真实气体效应对返回舱气动力特性的影响

采用7 组元6 反应化学动力学模型,通过数值方法研究了真实气体效应对阿波罗(Apollo) 返回舱流场及气动力特性的影响. 并利用典型弹道点的飞行和实验数据验证了化学非平衡流计算程序的可靠性. 计算结果表明:真实气体效应主要发生在物面附近很薄的激波层内,真实气体效应使得激波脱体距离减小;真实气体效应使阻力系数和升力系数增加,且在小攻角时增加幅度最大;真实气体效应产生附加的低头力矩,使压心位置后移. 真实气体效应的影响随着马赫数的增加变化不明显.     

Hydrodynamic and thermal fields analysis in gas-solid two-phase flow

The present work aims to investigate numerically the flowfield and heat transfer process in gas-solid suspension in a vertical pneumatic conveying pipe. The Eulerian-Lagrangian model is used to simulate the flow of the two-phases. The gas phase is simulated based on Reynolds Average Navier-Stokes equations (RANS) with low Reynolds number k-ε model, while particle tracking procedure is used for the solid phase. An anisotropic model is used to calculate the Reynolds stresses and the turbulent Prandtl number is calculated as a function of the turbulent viscosity. The model takes into account the lift and drag forces and the effect of particle rotation as well as the particles dispersion by turbulence effect. The effects of inter-particles collisions and turbulence modulation by the solid particles, i.e. four-way coupling, are also included in the model. Comparisons between different models for turbulence modulation with experimental data are carried out to select the best model. The model is validated against published experimental data for velocities of the two phases, turbulence intensity, solids concentration, pressure drop, heat transfer rates and Nusselt number distribution. The comparisons indicate that the present model is able to predict the complex interaction between the two phases in non-isothermal gas-solid flow in the tested range. The results indicate that the particle-particle collision, turbulence dispersion and lift force play a key role in the concentration distribution. In addition, the heat transfer rate increases as the mass loading ratio increases and Nusselt number increases as the pipe diameter increases.     

微尺度气体流动

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和Millikan的小球阻力系数的测量,实验结果揭示了稀薄气体效应即尺度效应对气体运动的重要影响.由于流动特征长度很小,微尺度气流经常处于滑流区甚至过渡领域,流动的相似参数为Knudsen数和Mach数.因此可以考虑利用相似准则,通过增大几何尺寸、减小压力的途径,解决微机电系统实验观测遇到的困难.为解决直接模拟Monte Carlo方法分析微机电系统中低速稀薄气流遇到的统计涨落困难,我们提出了信息保存法(IP),该方法能够有效克服统计散布,并已成功用于多种微尺度气流.      

Computation of wave interference and relaxation of particles after passing of a shock wave

Interaction of a shock wave with a system of motionless or relaxing particles is numerically simulated. Regimes of the gas flow around these particles are described, and the influence of the initial parameters of the examined phenomenon on the flow pattern is analyzed. The drag coefficient of particles is calculated as a function of the Mach number behind the shock wave at a fixed Reynolds number. The dynamics of heat exchange for particles of different sizes (10 μm–1 mm) is determined, and the laws of thermal relaxation after passing of a shock wave over the system of particles are found. The times of thermal and velocity relaxation of particles are estimated as functions of the Reynolds number, and the predicted relaxation time is compared with the corresponding empirical dependences.     

Hyperbolic Approximation of the Navier-Stokes Equations for Viscous Mixed Flows

Simplified two-dimensional Navier-Stokes equations of the hyperbolic type are derived for viscous mixed (with transition through the sonic velocity) internal and external flows as a result of a special splitting of the pressure gradient in the predominant flow direction into hyperbolic and elliptic components. The application of these equations is illustrated with reference to the calculation of Laval nozzle flows and the problem of supersonic flow past blunt bodies. The hyperbolic approximation obtained adequately describes the interaction between the stream and surfaces for internal and external flows and can be used over a wide Mach number range at moderate and high Reynolds numbers. Examples of the calculation of viscous mixed flows in a Laval nozzle with large longitudinal throat curvature and in a shock layer in the neighborhood of a sphere and a large-aspect-ratio hemisphere-cylinder are given. The problem of determining the drag coefficient of cold and hot spheres is solved in a new formulation for supersonic air flow over a wide range of Reynolds numbers. In the case of low and moderate Reynolds numbers a drag reduction effect is detected when the surface of the sphere is cooled.     

Regimes of sedimentation of a consolidated system of solid spherical particles

The results of the experimental investigation of gravitational sedimentation of a consolidated system of solid monodisperse spherical particles in a viscous liquid are represented over wide ranges of the particle number density and the Reynolds and Stokes numbers. Empirical dependences of the velocity of sedimentation of a particle aggregate and the drag coefficient of a system of particles as functions of the initial volume number density are obtained. The boundary values of the particle number density separating the sedimentation regimes are determined.