钝感炸药圆筒试验与爆轰产物JWL状态方程研究

 对以TATB为主的钝感炸药JB-9014进行了25 mm和50 mm两种装药直径的圆筒试验，测试了TU1圆筒在爆轰产物驱动下的膨胀过程R(t)关系。对圆筒试验进行了综合分析和二维数值模拟计算，综合评估了JB-9014的作功能力和圆筒试验的相似性。通过二维流体动力学数值模拟，确定得到了JB-9014炸药爆轰产物JWL状态方程参数。经对JB-9014二维平面滑移爆轰驱动试验的数值模拟检验，证明确定的JWL状态方程参数是可靠的，具有较高的精度和普适性。     

ON THE STRUCTURE OF STEADY FLOW THROUGH DUAL-INTAKE ENGINE PORTS

The influence of intake port design on the flow field in a dual-intake valve engine was investigated using computational fluid dynamics, in order to study the effect of inlet port design on the in-cylinder flow. A detailed 3D computational grid incorporating all the features of the Ford Zetec production engine inlet ports, valves and cylinder head was initially created and the flow structure modelled at 5 and 10 mm valve lifts under steady flow conditions. Comparisons of computational results with experimental data obtained by laser Doppler anemometry indicate that the flow characteristics have been predicted well in most regions. Flow generated by different intake port designs was also simulated by introducing air into the cylinder at different directions to the inlet valve axes and the effects of port deactivation, throttling and exhaust gas recirculation were examined. The implications of the results for intake port design are discussed.     

Effects of Thermal Radiation and Slip Mechanism on Mixed Convection Flow of Williamson Nanofluid Over an Inclined Stretching Cylinder

The current investigation highlights the mixed convection slip flow and radiative heat transport of uniformly electrically conducting Williamson nanofluid yield by an inclined circular cylinder in the presence of Brownian motion and thermophoresis parameter.A Lorentzian magnetic body force model is employed and magnetic induction effects are neglected.The governing equations are reduced to a system of nonlinear ordinary differential equations with associated boundary conditions by applying scaling group transformations.The reduced nonlinear ordinary differential equations are then solved numerically by Runge-Kutta-Fehlberg fifth-order method with shooting technique.The effects of magnetic field,Prandtl number,mixed convection parameter,buoyancy ratio parameter,Brownian motion parameter,thermophoresis parameter,heat generation/absorption parameter,mass transfer parameter,radiation parameter and Schmidt number on the skin friction coefficient and local Nusselt are analyzed and discussed.It is found that the velocity of the fluid decreases with decrease in curvature parameter,whereas it increases with mixed convection parameter.Further,the local Nusselt number decreases with an increase in the radiation parameter.The numerical comparison is also presented with the existing published results and found that the present results are in excellent agreement which also confirms the validity of the present methodology.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.     

内部爆炸作用下多层钢筒的动态响应

为评估内部爆炸作用下多层钢筒结构的防护效果,考察多层钢筒结构动态响应和变形吸能特征,采用两端开口、总厚度为50 mm的4层圆柱形Q345钢筒,在8.90~18.18 kg TNT药量下进行爆炸实验,并在容器外壁进行应变电测。实验后钢筒结构爆心局部发生塑性变形,内层钢筒变形最大,但未发生破坏。根据研究得到初步认识:采用爆心单位环面变形吸能的设计方法,可以较好地预估给定药量下所需钢筒的厚度;不同药量下,轴向距离超过多层钢筒结构的1/4内径后,其外壁环向变形峰值约减小为爆心截面环向应变峰值的1/2。     

LATTICE BOLTZMANN SIMULATION OF POWER-LAW FLOW PAST AN ELLIPTIC CYLINDER

基于格子波尔兹曼方法(lattice Boltzmann method) 和幂律流体本构方程, 建立二维流动模型, 将充分发展的速度分布与理论解进行对比, 吻合良好. 对幂律流体的圆柱绕流进行模拟, 采用了反弹格式的无滑移圆柱边界, 并使用应力积分法计算阻力系数, 分析了稳态圆柱绕流时, 阻力系数随幂律指数n 以及雷诺数Re 的变化规律. 分析了椭圆横轴/纵轴长度比和幂律指数n, 对压力系数C_p 和黏度系数C_v 的影响. 得到的变化规律与有限元方法规律一致, 验证了格子波尔兹曼模型的可行性.     

Large deflections of a beam subject to three-point bending

In this paper, a solution for the equilibrium configuration of an elastic beam subject to three-point bending is given in terms of Jacobi elliptical functions. General equations are derived, and the domain of the solution is established. Several examples that illustrate a use of the solution are discussed. The obtained numerical results are compared with the results of other authors. An approximation formula by which the beam load is given as a polynomial function of beam deflection is also derived. The range of applicability of the approximation is illustrated by numerical examples.     

Applications of scale-adaptive simulation technique based on one-equation turbulence model

A modified scale-adaptive simulation (SAS) technique based on the Spalart-Allmaras (SA) model is proposed. To clarify its capability in prediction of the complex turbulent flow, two typical cases are carried out, i.e., the subcritical flow past a circular cylinder and the transonic flow over a hemisphere cylinder. For comparison, the same cases are calculated by the detached-eddy simulation (DES), the delayed-detached eddy simulation (DDES), and the XY-SAS approaches. Some typical results including the mean pressure coefficient, velocity, and Reynolds stress profiles are obtained and compared with the experiments. Extensive calculations show that the proposed SAS technique can give better prediction of the massively separated flow and shock/turbulent-boundary-layer interaction than the DES and DDES methods. Furthermore, by the comparison of the XY-SAS model with the present SAS model, some improvements can be obtained.     

Influence of aspect ratio on the flow above the free end of a surface-mounted finite cylinder

The flow above the free end of a surface-mounted finite-height cylinder was studied in a low-speed wind tunnel using particle image velocimetry (PIV). Velocity measurements were made in vertical and horizontal measurement planes above the free end of finite cylinders of aspect ratios AR = 9, 7, 5 and 3, at a Reynolds number of Re = 4.2 × 10⁴. The relative thickness of the boundary layer on the ground plane was δ/D = 1.7. Flow separating from the leading edge formed a prominent recirculation zone on the free-end surface. The legs of the mean arch vortex contained within the recirculation zone terminate on the free-end surface on either side of the centreline. Separated flow from the leading edge attaches onto the upper surface of the cylinder along a prominent attachment line. Local separation downstream of the leading edge is also induced by the reverse flow and arch vortex circulation within the recirculation zone. As the cylinder aspect ratio is lowered from AR = 9 to AR = 3, the thickness of the recirculation zone increases, the arch vortex centre moves downstream and higher above the free-end surface, the attachment position moves downstream, and the termination points of the arch vortex move upstream. A lowering of the aspect ratio therefore results in accentuated curvature of the arch vortex line. Changes in aspect ratio also influence the vorticity generation in the near-wake region and the shape of the attachment line.