Dynamics of the rear part of a long-rod penetrator’s core during firing

A problem of the dynamics of long-rod penetrators loaded with the gunpowder combustion products has been solved in closed form by using the wave method (d’Alambert). The analytic recurrence formulae have been derived to define the fields of longitudinal stress and particle velocity in the rod penetrator. The problem of estimating the axial stress field and particle velocity caused by the gunpowder combustion products for a subcalibre projectile penetrator has been solved by using the substitute dynamic theoretical model. The reason for using the substitute dynamic theoretical model is that the pressure on the rear part of the projectile is controlled by gas-dynamic processes and it is not possible to readily estimate the loading. The method of characteristics net gets very complicated.      

PREDICTION OF THE VISCOELASTIC PROPERTIES OF THE EQUIVALENT PARTICLE FOR THE INTERCALATED MULTI-LAYER STACK OF NANOPLASTICS

The aim of this paper is to apply the asymptotic homogenization method to deter- mining analytically and numerically the transversely isotropic viscoelastic relaxation moduli of the equivalent particle for the intercalated multi-layer stack of intercalated type nanoplastics.A two-phase multilayered material containing n layers is considered.The matrix is assumed to be an isotropic viscoelastic standard linear body and the reinforcement is assumed to be an isotropic elastic body.Final explicit analytical formulae for the effective elastic moduli of the multilay- ered material are derived first;and then the correspondence principle is employed to obtain the homogenized relaxation moduli of the equivalent intercalated particle.A numerical example is given.Final explicit analytical formulae in the time domain derived here make it convenient to estimate the influence of all the particle parameters of micro-structural details on the effective properties of the equivalent intercalated particle.The results of this paper can also be applied to multi-layer composites.     

Kinetic modeling of particles in stratified flow – Evaluation of dispersion tensors in inhomogeneous turbulence

The continuum equations for a dilute particle distribution in inhomogeneous turbulence are tested against results from a Langevin particle tracking simulation. Reeks’ version of the kinetic theory is used to generate the mass, momentum and kinetic stress equations for the particle distribution. The particle tracking data are used to directly evaluate the dispersion tensors λ and μ which serve as closure relations for the continuum equations. These exact forms are compared to approximate, local forms. Even for low Stokes numbers (corresponding to low particle inertia defined by τ/τ_p ? 1), the tensor λ is strongly affected by the inhomogeneity and depends on turbulence parameters in the volume corresponding to the particle path dispersion over the particle Lagrangian integral timescale τ. In contrast, the locally homogeneous form of the velocity dispersion tensor μ is a sufficient approximation, since it depends on the dispersion volume over the much smaller particle relaxation time τ_p. It is demonstrated that the body force due to the dispersion vector γ cannot be neglected. In the limit of passive tracers (zero stopping distance), γ is equal to the gradient of λ, if the physical setting is such that we can invoke constant tracer density in this limit.     

D'alembert principle in the velocity space

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｒｅｓｅａｒｃｈｅｒｓｈａｖｅｂｅｅｎｍａｄｅｏｎｔｈｅｖａｒｉｏｕｓｆｏｒｍｕｌａｅｏｆＤ’Ａｌｅｍｂｅｒｔｐｒｉｎｃｉｐｌｅ,ａｔｔｅｎｔｉｏｎｈａｓｂｅｅｎｇｉｖｅｎｔｏｍｅｃｈａｎｉｃａｌｐｒｉｎｃｉｐｌｅｓｉｎｔｈｅｖｅｌｏｃｉｔｙｓｐａｃｅａｎｄｔｈｅｖａｒｉａｔｉｏｎａｌｐｒｉｎｃｉｐｌｅｉｎｔｈｅｖｅｌｏｃｉｔｙｓｐａｃｅ,ａｎｄｔｈｅｖｅｌｏｃｉｔｙｓｐａｃｅｏｆｒｅｌａｔｉｖｉｓｔｉｃｍｅｃｈａｎｉｃｓｈａｓｂｅｅｎｍａｄｅ…     

基于离散单元法和人工神经网络的近壁颗粒动力学特征研究

颗粒与壁面的相互作用往往对颗粒流动具有显著影响. 为研究颗粒与壁面作用机理, 对滚筒内颗粒流动过程进行离散单元法(DEM)数值模拟. 基于模拟结果统计分析靠近壁面处颗粒的运动特征, 结果表明, 小摩擦系数时颗粒平动和旋转速度均近似满足正态分布, 但由于壁面影响, 摩擦系数增大时颗粒沿滚筒轴向的旋转速度偏离正态分布, 颗粒动力学理论推导壁面边界条件时应考虑速度正态分布的修正及速度脉动的各向异性. 采用人工神经网络(ANN)构建了颗粒无因次旋转温度、滑移速度和平动温度之间的函数模型, 进而可以在常规双流模型壁面边界条件中考虑颗粒旋转的影响. 基于DEM模拟及结果分析可以为壁面边界条件的理论构造和半经验修正提供基础数据和封闭模型.      

Semi-theoretical analyses of the concrete plate perforated by a rigid projectile

Based on the three-stage perforation model, a semi-theoretical analysis is conducted for the ballistic performances of a rigid kinetic projectile impacting on concrete plates. By introducing the projectile resistance coefficients, dimensionless formulae are proposed for depth of penetration (DOP), perforation limit thickness, ballistic limit velocity, residual velocity and perforation ratio, with the projectile nosed geometries and projectile-target interfacial friction taken into account. Based on the proposed formula for DOP and lots of penetration tests data of normal and high strength concrete targets, a new expression is obtained for target strength parameter. By comparisons between the results of the proposed formulae and existing empirical formulae and large amount of projectile penetration or perforation tests data for monolithic and segmented concrete targets, the validations of the proposed formulae are verified. It is found that the projectile-target interfacial friction can be neglected in the predictions of characteristic ballistic parameters. The dimensionless DOP for low-to-mid speed impacts of non-flat nosed projectiles increases almost linearly with the impact factor by a coefficient of 2/(πS). The anti-perforation ability of the multilayered concrete plates is dependent on both the target plate thickness and the projectile impact velocity. The variation range of the perforation ratio is 1–3.5 for concrete targets.     

柱形长杆弹侵彻的界面击溃分析

在Alekseevski-Tate模型的基础上,分析了柱形长杆弹的界面击溃过程,给出了弹体速度下降及质 量侵蚀的计算公式;讨论了弹体速度下降及质量侵蚀对动能损失的影响;特别针对柱形长杆弹在界面击溃过 程中弹体速度准定常小量变化的特点,近似给出了弹体速度、弹体质量随时间变化的简化解析表达式,为工程 应用提供便利。     

A numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 John Wiley & Sons, Ltd.     

一种基于湍动能方程的转捩判定方法简

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.     

高压驱动研究方法的改进及其对质点运动的研究

从分析二维运动的测试原理入手 ,利用高速狭缝扫描摄影技术建立了测试系统 ,该系统能获得二维运动的轨迹和速度 ,从而把高速狭缝扫描摄影的研究领域由一维拓展到了二维 ,并对不同种测试系统进行了比较。利用新的测试系统获得了物体运动的照片 ,研究了质点高压下的运动 ,同时 ,获得了变形速度和变形角。     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

Lagrangian Covariance Analysis of β-Plane Turbulence

The effects of Rossby wave–turbulence interactions on particle dispersion are investigated in a Lagrangian analysis of the potential vorticity equation. The analysis produces several exact statistical results for fluid particle dispersion in barotropic turbulence on a β-plane. In the inviscid problem the first integral time scale of the meridional velocity is found to be zero, as might occur in pure wave processes, and the meridional particle dispersion is bounded. The second integral time scale, which determines the magnitude of the bound, is shown to depend explicitly on β, the enstrophy and the energy of the meridional velocity. Expressions relating the autocorrelation of the vorticity to the autocorrelation of the meridional velocity are derived and the Lagrangian integral time scale of the relative vorticity is diagnostically related to the meridional velocity correlation. The applicability of these predictions is verified in a series of numerical simulations. For a range of β values, the meridional extent of quasisteady alternating zonally averaged jets occurring in the numerical solutions scales with a length scale given by the the standard deviation of the meridional particle dispersion. Received 8 March 1999 and accepted 12 December 1999     

Relaxation of a fluidized bed to the state of a continuous medium

The kinetic equation proposed in [1,2] for describing the behavior of a system of particles in a gas flow differs from the usual Boltzmann equation with respect to the additional terms that take into account random variations of the particle velocity under the influence of the flow. As shown in [2], the collision operator and the Brownian-type operator in the starting kinetic equation describe essentially different simultaneous physical processes of change of state of the particle system: equalization of the mean kinetic energy of the particles and change of energy due to the action of the viscous forces associated with the suspending flow. Therefore the method of solving the kinetic equation used in [2], a direct generalization of the Chapman-Enskog method of solving the kinetic equation it is necessary to investigate method of solving the kinetic equation it is necessar y to investigate the relaxation processes in the system. Moreover, the relaxation of systems of the fluidized-bed type to the continuum state is also of independent interest in connection with the analysis of fast processes in the system, i.e., processes with a characteristic duration of the order of the mean free time.     

Measurement of instantaneous Eulerian acceleration fields by particle image accelerometry: method and accuracy

Acceleration is a fundamental quantity in fluid mechanics because it reflects the sum of all forces (pressure and viscous) present within the flow. However, measurements of acceleration have been difficult to achieve relative to the ease with which fluid velocity can be measured. A particle image accelerometer (PIA) has been developed to measure Eulerian acceleration fields by time-differencing successive measurements of the Eulerian velocity field as measured by particle image velocimetry (PIV). The measurements can also be made in uniformly translating frames. With current video camera technology, it is often not possible to measure the two velocity fields with a time separation sufficiently small enough to permit accurate finite difference approximation of the time derivative. A two-CCD-camera system has been developed to alleviate this limitation. Polarization filtering is utilized to separate the particle images viewed by each camera. The polarization filtering is achieved using cross-polarized light-sheets and a polarization filter just upstream of the imaging optics of the cameras. In this manner, PIV measurements can be achieved easily at time delays several orders of magnitude smaller than the shutter-time of the CCD cameras. The accuracy of the acceleration measurements is determined by numerical finite differencing errors and random noise and bias errors associated with the measurement of velocity. These errors, and methods of compensating for them, are studied.     

Experiments and simulations of gas-solid flow in an airlift loop reactor

The hydrodynamics in a gas-solid airlift loop reactor was investigated systematically using experimental measurements and CFD simulation.In the experiments,the time averaged parameters,such as solid fraction and particle velocity,were measured by optical fiber probe.In the simulation,the modified Gidaspow drag model accounting for the interparticles clustering was incorporated into the Eulerian-Eulerian CFD model with particulate-phase kinetic theory.Predicted values of solid fraction and particle velocity ...     

Theory on transverse migration of particles in a turbulent two-phase suspension flow due to turbulent diffusion—I

A new theoretical model has been developed to explain the behavior of transverse particle transport in turbulent flow of a dilute two-phase suspension due to turbulent diffusion. This model is based on the ability of a particle to respond to surrounding fluid motion and depends on particle size and density relative to the carrier fluid, the fractional variation in particle concentration in the transverse direction as well as the existing turbulence structure of the surrounding fluid. The model developed in this investigation has been formulated by dividing the transverse fluid velocity, as seen by a particular particle, into two superimposed components representing, respectively, the transverse turbulent fluid fluctuations and an apparent transverse local fluid drifting velocity due to the effect on the transverse oscillatory component of fluid motion by the transverse concentration distribution of particles. A subsequent paper will show that the theory (together with other new results on the concentration effects on particle drag and lift and fluid turbulence properties) can help to explain the phenomena measured previously.     

Statistical models for predicting the effect of bidisperse particle collisions on particle velocities and stresses in homogeneous anisotropic turbulent flows

The purpose of this paper is to present and compare two statistical models for predicting the effect of collisions on particle velocities and stresses in bidisperse turbulent flows. These models start from a kinetic equation for the probability density function (PDF) of the particle velocity distribution in a homogeneous anisotropic turbulent flow. The kinetic equation describes simultaneously particle–turbulence and particle–particle interactions. The paper is focused on deriving the collision terms in the governing equations of the PDF moments. One of the collision models is based on a Grad-like expansion for the PDF of the velocity distributions of two particles. The other model stems from a Grad-like expansion for the joint fluid–particle PDF. The validity of these models is explored by comparing with Lagrangian simulations of particle tracking in uniformly sheared and isotropic turbulent flows generated by LES. Notwithstanding the fact that the fluid turbulence may be isotropic, the particle velocity fluctuations are anisotropic due to the impact of gravitational settling. Comparisons of the model predictions and the numerical simulations show encouraging agreement.     

Head-on Collision of a Detonation with a Planar Shock Wave

The phenomenon that occurs when a Chapman–Jouguet (CJ) detonation collides with a shock wave is discussed. Assuming a one-dimensional steady wave configuration analogous to a planar shock–shock frontal interaction, analytical solutions of the Rankine–Hugoniot relationships for the transmitted detonation and the transmitted shock are obtained by matching the pressure and particle velocity at the contact surface. The analytical results indicate that there exist three possible regions of solutions, i.e. the transmitted detonation can have either strong, weak or CJ solution, depending on the incident detonation and shock strengths. On the other hand, if we impose the transmitted detonation to have a CJ solution followed by a rarefaction fan, the boundary conditions are also satisfied at the contact surface. The existence of these multiple solutions is verified by an experimental investigation. It is found that the experimental results agree well with those predicted by the second wave interaction model and that the transmitted detonation is a CJ detonation. Unsteady numerical simulations of the reactive Euler equations with both simple one-step Arrhenius kinetic and chain-branching kinetic models are also carried out to look at the transient phenomena and at the influence of a finite reaction thickness of a detonation wave on the problem of head-on collision with a shock. From all the computational results, a relaxation process consisting of a quasi-steady period and an overshoot for the transmitted detonation subsequent to the head-on collisions can be observed, followed by the asymptotic decay to a CJ detonation as predicted theoretically. For unstable pulsating detonations, it is found that, due to the increase in the thermodynamic state of the reactive mixture caused by the shock, the transmitted pulsating detonation can become more stable with smaller amplitude and period oscillation. These observations are in good agreement with experimental evidence obtained from smoked foils where there is a significant decrease in the detonation cell size after a region of relaxation when the detonation collides head-on with a shock wave.     

一种基于湍动能方程的转捩判定方法

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.