飞行加速度对发动机内粒子轨迹影响研究

利用Euler-Lagrange模型方程，对发动机内的三维两相流动进行了数值模拟，分析了不同过载组合、不同粒子直径、不同壁面恢复系数对发动机燃烧室内粒子的聚集部位、浓度分布、冲刷速度、冲刷角度等的影响。计算结果表明：在纵向和横向过载相等的情形下，粒子的最大聚集密度随过载量的增加，呈先减小后增加的趋势，最大聚集点随过载量的增加向流动方向移动，过载量的大小对粒子的冲刷角度、冲刷速度影响较大；在相同工况下，随粒子直径的增加，最大聚集部位向流动方向移动，最大聚集密度增加；随粒子壁面恢复系数恢复得减小，最大聚集部位向发动机头部移动，最大聚集密度增加。     

局部扩张动脉管血液振荡流的切应力分布

本文求解局部缓慢扩张动脉管中血液振荡流的基本方程，得到血管内血液的流速与压力梯度的关系。通过导出压力梯度沿局部扩张管轴向的变化特性。建立利用扩张段上游血管均匀段中心流速波形确定局部扩张管中血液流的速度和切应力分布的方法，文章以人体颈动脉余弦扩张为例进行分析。详细讨论了局部扩张对血管壁切应力及其梯度分布的影响。数值结果表明，在与刚性均匀管中管壁切应力沿轴向保持不变不同，在局部扩张段，管壁切应力将随着血管半径的增大而减小，因而管壁切应力梯度一般不为零，甚至在某些位置达到相当大的数值。另外，随着血管扩张程度的增加，管壁切应力还将进一步减小，而且管壁切应力梯度也将进一步增大，血管扩张导致管壁切应力的这些变化将直接影响血管壁的结构和功能，使其产生适应性的变化。     

三爆轰管-脉冲爆轰发动机的近场噪声特性

为研究三爆轰管脉冲爆轰发动机噪声的形成机理和特性,设计了正三角形组合的三爆轰管PDE实验系统,分别在与中心轴线成0°、30°、60°及90°方向上,对8个不同距离处的声压进行了测量,结果表明:不同距离处噪声幅值最大值均在30°方向上;爆轰噪声的参考半径为3倍“名义管径”;爆轰噪声的A持续时间随着r(距离三爆轰管名义“中点”的直线距离)的增加近似线性减小,并随着角度的增大而减小。B持续时间与声压值的大小成反比关系,并且随着轴向距离的增加而增加。在不同的方位角上,峰值声压越大,B持续时间越小。     

基于流固耦合的燃气冲刷烧蚀内膛特性分析

火炮发射时，火药燃气与身管间发生剧烈的传热传质作用是导致身管烧蚀的重要因素。为了研究某155 mm火炮中高温高压高速的燃气流对身管的烧蚀特性，采用CFD流固耦合方法，建立了发射过程中的身管非稳态流动传热模型，并根据炮钢在不同温度下的烧蚀特点，将烧蚀过程分为热化学烧蚀和熔化烧蚀，建立了分段烧蚀模型。计算结果表明，身管内壁温度随时间的增加先迅速增大，随后逐渐降低。整体上，内壁温度随身管轴向距离的增大而逐渐降低。身管膛线起始区域的壁面温度最高，其烧蚀是熔化和热化学烧蚀共同导致的，而线膛部的大部分区域仅发生了热化学烧蚀。总烧蚀量随着身管轴向距离的增大而逐渐降低，膛线起始部的烧蚀最为严重，单发总烧蚀量（常温）为5.06μm。同时分析了不同工况对身管烧蚀特性的影响，发现最大烧蚀量与初始壁面温度呈现很强的正相关性，温度的升高会加剧身管的烧蚀。     

受限空间浮力羽流准稳态速度场的实验研究

利用作者自己发展的一种热球式低风速多点实时测量系统，测量了受限空间浮力羽流各横截面上竖直方向的速度分布。实验中改变火源相对于侧壁的距离以及火源的功率，发现火源距靠近的壁面越远，则浮力羽流横截面上的速度分布越趋于平缓，沿侧壁面Ｘ方向的速度分布是近似高斯分布，垂直壁面Ｙ方向的速度分布偏离高斯分布。羽流轴线向靠近的壁面偏斜，羽流轴线相对于火源中心线Ｚ轴之间的距离定义为轴线偏移距离δ。δ随火源与侧壁之间的距离Ｙ的增加先是增加然后下降；δ随火源功率的增加而增加。实验结果显示了受限空间浮力羽流场的基本特性。     

纳米粒子在方管中传输与沉降的研究

采用有限体积法离散并应用Simple方法对方截面弯曲管道内的纳米粒子传输和沉降进行了数值计算,结果表明Reynolds数和Schmidt数是影响纳米粒子传输和沉降的重要参数。粒子较小时,弯管中轴向速度较大的区域就是粒子的高浓度区域,沉降增强因子最大值出现在外弯侧的中心位置;粒子较大时,截面浓度的梯度值降低,沉降增强因子趋于平均,此时整个截面的粒子平均沉降。弯曲作用对于粒子较大且Dean数也较大时的影响更加明显。     

基于光滑粒子动力学的煤岩刨刀刨削过程数值模拟

为研究刨削参数对刨刀载荷特性的影响规律，对刨刀破煤过程进行仿真，为得到更可靠的仿真数据，采用光滑粒子动力学(SPH)与有限元进行耦合(FEM)对刨刀破煤过程进行模拟，并将仿真结果与传统有限元算法进行比较证明此算法的优越性。通过模拟不同刨削参数下的刨刀的破煤过程，得到刨刀载荷数据。研究结果表明：随着刨削速度的增加，刨刀在刨削过程中的整体受载变化不大，但刨刀所受载荷峰值和波动范围随着随刨削速度的增加而增大；随着刨削深度的增加，刨刀在刨削过程中的整体受载有着明显的变化，具体表现为：刨刀所受载荷峰值和均值随着随刨削深度的增加而增大，但刨削深度对刨刀载荷波动范围影响不明显。     

除尘管道内铝粉运移及沉积的数值研究

基于DPM(Discrete Phase Model)模型，研究了长直通风管道内粒径服从Rosin-Rammler 分布的铝粉的运移与沉积规律．基于颗粒与壁面的碰撞过程中的能量分析，建立了粉尘沉积-回弹模型，得出了粉尘沉积的判定准则及脱离壁面时的回弹速度．利用UDF 将沉积-回弹模型嵌入Fluent，完成了对管道内粉尘运移和沉积的数值模拟．粉尘沉积的数值结果与实验结果符合得较好，验证了所提模型的有效性．数值结果表明风速的增大使管道内粉尘浓度明显降低，管壁粉尘沉积率也降低；粒径的增大对粉尘浓度的大小影响不明显，主要影响粉尘浓度在管道内的分布情况，同时会增大粉尘在管壁的沉积率．     

水平流动边界层内气固相间作用的试验研究

应用三维粒子动态分析仪（threedimensionalparticledynamicsanalyzer）,测量了含有230μm颗粒的气固两相水平流的特性,特别是壁面边界层内的两相流动特性．结果表明颗粒载荷比（质量流率）对相间作用有较大影响,随颗粒流率的增加颗粒对气流平均速度和湍流的影响增大,颗粒使气流速度边界展变薄．颗粒和气流相互作用在不同方向上呈各向异性,颗粒对气流垂直方向的脉动影响较大．颗粒与湍流边界层气流的作用行为大致可以分成三个区域：贴壁区、中间区和外流区．     

壁面湍流标度律的PIV实验研究

利用高分辨率、高帧率PIV系统对湍流边界层内的自相似标度律在不同的壁面距离的具体表达形式进行了实验研究。在动量损失雷诺数（Reθ=2167）下测量平板湍流边界层内的二维瞬时速度场。应用小波分析和传统统计学方法,在垂直于平板的平面内,对不同法向位置的自相似的标度律的具体形式及其随尺度的变化情况进行考察,并与已知的She-leveque（简称SL）标度律进行比较分析。结果表明,越接近壁面,结构函数的概率密度曲线尾部越宽,表现出逐渐增强的间歇性;并且随着壁面距离的减小,ESS标度律曲线与SL标度律曲线之间的差别也就越大,而新形式标度律曲线则与SL标度律曲线的差别逐渐减小。在提取的结构尺度一样时,结构函数的标度指数随法向位置的不同变化不大,但是随着提取结构尺度的增加,两种形式的标度律适用的尺度范围也在逐渐增大。     

Numerical simulation of particle sedimentation in 3D rectangular channel

The 3D lattice Boltzmann method is used to simulate particle sedimentation in a rectangular channel. The results of single particle sedimentation indicate that the last position of the particle is along the center line of the channel regardless of the initial position, the particle diameter, and the particle Reynolds number. The wall effect on the terminal velocity is in good agreement with experimental results quantitatively. The drafting, kissing, and tumbling (DKT) process is reproduced and analyzed by simulating two-particle cluster sedimentation. The effects of the diameter ratio, initial position, and wall on the DKT process are investigated. When the two particles have equal diameter sediment in the rectangular channel, a periodical DKT process and the spiraling trajectory are found. The last equilibrium configuration is obtained from the simulation results. The interesting regular sedimentation phenomena are found when 49 particles fall down under gravity.     

Pipe flow of suspensions

This study shows that fully developed pipe flow of a particulate suspension is defined by four dimensionless parameters of particle-fluid interactions in addition to the Reynolds number. Effects accounted for include the Magnus effect due to fluid shear, electrostatic repulsion due to electric charges on the particles, and Brownian or turbulent diffusion. In the case of a laminar liquid-solid suspension, electrostatic effect is negligible but shear effect is prominent. Solution of the basic equations gives the density distribution of particles with a peak at the center (Einstein, Jeffery) or at other radii between the center and the pipe wall (Segré et al) depending on the magnitudes of the various flow parameters. In the case of a turbulent gas-solid suspension, the Magnus effect is significant only within the thickness of the laminar sublayer. However, charges induced on the particles by the impact of particles at the wall produce a higher density at the wall than at the center of the pipe. The velocity distribution of particles is characterized by a slip velocity at the wall and a lag in velocity in the core from the fluid phase. These results are verified by earlier measurements.     

Pipe flow of suspensions

This study shows that fully developed pipe flow of a particulate suspension is defined by four dimensionless parameters of particle-fluid interactions in addition to the Reynolds number. Effects accounted for include the Magnus effect due to fluid shear, electrostatic repulsion due to electric charges on the particles, and Brownian or turbulent diffusion. In the case of a laminar liquid-solid suspension, electrostatic effect is negligible but shear effect is prominent. Solution of the basic equations gives the density distribution of particles with a peak at the center (Einstein, Jeffery) or at other radii between the center and the pipe wall (Segré et al) depending on the magnitudes of the various flow parameters. In the case of a turbulent gas-solid suspension, the Magnus effect is significant only within the thickness of the laminar sublayer. However, charges induced on the particles by the impact of particles at the wall produce a higher density at the wall than at the center of the pipe. The velocity distribution of particles is characterized by a slip velocity at the wall and a lag in velocity in the core from the fluid phase. These results are verified by earlier measurements.     

An analytical approach for the closure equations of gas–solid flows with inter-particle collisions

This work examines the effect of inter-particle collisions on the motion of solid particles in two-phase turbulent pipe and channel flows. Two mechanisms for the particle–particle collisions are considered, with and without friction sliding. Based on these collision mechanisms, the correlations of the various velocity components of colliding particles are obtained analytically by using an averaging procedure. This takes into account three collision coordinates, two angles and the distance between the centers of colliding particles. The various stress tensor components are obtained and then introduced in the mass, linear momentum and angular momentum equations of the dispersed phase. The current approach applies to particle–particle collisions that result from both the average velocity difference and the turbulent velocity fluctuations. In order to close the governing equations of the dispersed phase, the pseudo-viscosity coefficients are defined and determined by the time of duration of the inter-particle collision process. The model is general enough to apply to both polydisperse and monodisperse particulate systems and has been validated by comparisons with experimental data.     

Particle behavior in a turbulent pipe flow with a flat bed

Particle behavior in a turbulent flow in a circular pipe with a bed height h = 0.5R is studied at Re_b = 40,000 and for two sizes of particles (5 μm and 50 μm) using large eddy simulation, one-way coupled with a Lagrangian particle tracking technique. Turbulent secondary flows are found within the pipe, with the curved upper wall affecting the secondary flow formation giving rise to a pair of large upper vortices above two smaller vortices close to the pipe floor. The behavior of the two sizes of particle is found to be quite different. The 50 μm particles deposit forming irregular elongated particle streaks close to the pipe floor, particularly at the center of the flow and the pipe corners due to the impact of the secondary flows. The deposition and resuspension rate of the 5 μm particles is high near the center of the floor and at the pipe corners, while values for the 50 μm particles are greatest near the corners. Near the curved upper wall of the pipe, the deposition rate of the 5 μm particles increases in moving from the wall center to the corners, and is greater than that for the larger particles due to the effects of the secondary flow. The maximum resuspension rate of the smaller particles occurs above the pipe corners, with the 50 μm particles showing their highest resuspension rate above and at the corners of the pipe.     

非球形固体颗粒对弯管内壁的磨损机制

为研究弯管内固体颗粒在液相夹带条件下的运动特性及颗粒对弯管内壁的磨损,采用计算流体动力学与离散元耦合的方法,建立数值模型,考虑固液两相之间的作用,对弯管内的固液两相流动进行数值模拟;通过软件的应用程序编程接口嵌入自编译磨损模型;借助试验结果,验证数值模型的有效性. 结果表明,所建立的数值模拟方案可以准确地模拟颗粒在管内的运动特征并能够预测弯管内壁的磨损位置以及磨损程度. 弯管内的二次流对颗粒运动有重要影响,弯管外侧壁面中心线附近的磨损较严重,磨损的形式以小角度划擦切削为主. 弯管磨损主要与颗粒对壁面的碰撞速度、碰撞角度及碰撞频率有关. 运动中的颗粒与壁面发生多次碰撞,碰撞角度逐渐减小. 随着颗粒球形度的增大,在相同碰撞条件下引起的磨损量变小,但是会降低颗粒的随流性. 颗粒形状影响颗粒在流场中的运动速度以及颗粒与壁面的碰撞. 随着颗粒球形度增大,严重磨损区域向弯管进口方向移动,壁面平均磨损量先减小后增大;当输送颗粒的球形度为0.91时,壁面磨损量最小.      

The effect of sphere–wall interactions on particle motion in a viscoelastic suspension of FENE dumbbells

A new method for simulating the motion of particles in viscoelastic Boger fluids is extended to problems with bounded geometries. Viscoelasticity is incorporated into the Stokesian dynamics method by modeling a viscoelastic fluid as a suspension of finite-extension nonlinear-elastic (FENE) dumbbells. Wall–particle and wall–bead interactions are included by using the image system method of Blake; particle–particle and particle–bead interactions are also modified by the presence of the wall. The method of incorporating sphere–wall interactions is verified by doing calculations for several problems involving particle–wall interactions in Newtonian fluids. The method is then used to study particle–wall interactions in viscoelastic dumbbell suspensions by examining several problems of interest: the sedimentation of a spherical particle near vertical and tilted walls; the sedimentation of a nonspherical particle between two flat plates; and the migration of a neutrally buoyant sphere in plane Poiseuille flow. We find that a single sphere falling near a wall moves toward the wall and exhibits anomalous rotation. When the wall is tilted by an amount less than a few degrees, the sphere still moves toward the wall, but tilting the wall greater than an angle of approximately 1.5° results in the sphere falling away from the wall. A nonspherical particle settling in a channel exhibits an oscillatory motion, but ultimately becomes centered in the channel with its long axis parallel to gravity. Finally, it is shown that a neutrally buoyant sphere in plane Poiseuille flow migrates to the channel center in wide channels, but migrates to the walls when the sphere is sufficiently large relative to the channel width.     

Modelling of gas–solid turbulent channel flow with non-spherical particles with large Stokes numbers

This paper describes a complete framework to predict the behaviour of interacting non-spherical particles with large Stokes numbers in a turbulent flow. A summary of the rigid body dynamics of particles and particle collisions is presented in the framework of Quaternions. A particle-rough wall interaction model to describe the collisions between non-spherical particles and a rough wall is put forward as well. The framework is coupled with a DNS-LES approach to simulate the behaviour of horizontal turbulent channel flow with 5 differently shaped particles: a sphere, two types of ellipsoids, a disc, and a fibre. The drag and lift forces and the torque on the particles are computed from correlations which are derived using true DNS.The simulation results show that non-spherical particles tend to locally maximise the drag force, by aligning their longest axis perpendicular to the local flow direction. This phenomenon is further explained by performing resolved direct numerical simulations of an ellipsoid in a flow. These simulations show that the high pressure region on the acute sides of a non-spherical particle result in a torque if an axis of the non-spherical particle is not aligned with the flow. This torque is only zero if the axis of the particle is perpendicular to the local direction of the flow. Moreover, the particle is most stable when the longest axis is aligned perpendicular to the flow.The alignment of the longest axis of a non-spherical particle perpendicular to the local flow leads to non-spherical particles having a larger average velocity compared to spherical particles with the same equivalent diameter. It is also shown that disc-shaped particles flow in a more steady trajectory compared to elongated particles, such as elongated ellipsoids and fibres. This is related to the magnitude of the pressure gradient on the acute side of the non-spherical particles. Finally, it is shown that the effect of wall roughness affects non-spherical particles differently than spherical particles. Particularly, a collision of a non-spherical particle with a rough wall induces a significant amount of rotational energy, whereas a corresponding collision with a spherical particle results in mostly a change in translational motion.     

Simulation of particles released near the wall in a turbulent boundary layer

A direct numerical simulation was used along with a Lagrangian particle tracking technique to study particle motion in a horizontal, spatially developing turbulent boundary layer along an upper-wall (with terminal velocity directed away from the wall). The objective of the research was to study particle diffusion, dispersion, reflection, and mean velocity in the context of two parametric studies: one investigated the effect of the drift parameter (the ratio of particle terminal velocity to fluid friction velocity) for a fixed and finite particle inertia, and the second varied the drift parameter and particle inertia by the same amount (i.e. for a constant Froude number). A range of drift parameters from 10⁻⁴ to 10⁰ were considered for both cases. The particles were injected into the simulation at a height of four wall units for several evenly distributed points across the span and a perfectly elastic wall collision was specified at one wall unit.Statistics collected along the particle trajectories demonstrated a transition in particle movement from one that is dominated by diffusion to one that is dominated by gravity. For small and intermediate sized particles (i.e. ones with outer Stokes numbers and drift parameters much less than unity) transverse diffusion away from the wall dominated particle motion. However, preferential concentration is seen near the wall for intermediate-sized particles due to inhomogeneous turbulence effects (turbophoresis), consistent with previous channel flow studies. Particle–wall collision statistics indicated that impact velocities tended to increase with increasing terminal velocity for small and moderate inertias, after which initial conditions become important. Finally, high relative velocity fluctuations (compared to terminal velocity) were found as particle inertia increased, and were well described with a quasi-one-dimensional fluctuation model.