动力刚化问题的实验研究

应用频散可控耗散格式对环形激波在圆柱形激波管内绕射、反射和聚焦 的问题进行了数值模拟研究. 研究结果表明环形激波形成强烈聚焦的关键因素是环形激波在 圆柱形管道中向对称轴运动时，绕射激波就不断加速而不作通常情况下的衰减；不同马赫数 的环形激波绕射也产生不同马赫数及形状的准柱形激波，导致聚焦效果和位置的差异；另外， 环形激波聚焦于一个点而圆柱形激波聚焦于一条线,两者有本质不同.     

激波与堆积粉尘相互作用的数值模拟

基于双流体模型和测定的堆积粉尘的本构方程 ,利用AUSM+ 格式 ,对激波与堆积粉尘的相互作用进行了数值模拟。计算所反映的流场结构与实验图像一致。此外还对激波强度 ,颗粒材料密度等对流场的影响进行了讨论。     

环形激波绕射, 反射和聚焦的数值模拟研究

应用频散可控耗散格式对环形激波在圆柱形激波管内绕射、反射和聚焦的问题进行了数值模拟研究，研究结果表明环形激波形成强烈聚焦的关键因素是环形激波在圆柱形管道中向对称轴运动时，绕射激波就不断加速而不作通常情况下的衰减；不同马赫数的环形激波绕射也产生不同马赫数及形状的准柱形激波，导致聚焦效果和位置的差异；另外，环形激波聚焦于一个点而圆柱形激波聚焦于一条线，两者有本质不同。     

激波和湍流相互作用的数值模拟

本文综述了关于激波和湍流相互作用数值模拟的近期研究进展, 主要包括激波和均匀各向同性湍流、激波和湍流边界层、激波和射流以及激波和尾迹的相互作用. 激波和湍流相互作用特性受到诸多因素的影响,如激波的强度、位置、形状和流动边界以及来流的湍流状态和可压缩性等. 激波和湍流的相互作用会引起流场结构、激波特性和湍流统计特性的显著变化. 最后简要讨论了激波和湍流相互作用数值研究需要关注的一些问题.      

可压缩流向涡与激波轴对称干扰的数值模拟1)

用NS方程数值模拟了可压缩流向涡和激波轴对称相互作用现象.数值模拟包括定常和非定常两种情况,计算结果分别与相应的实验进行了比较.结果表明数值模拟成功地捕捉到了激波和旋涡相互作用过程中发生的激波波面变形,激波振荡,涡核变大以及激波波后出现驻点、回流区等流场特征.提出了判断流向涡与运动激波相互作用中旋涡破碎的准则.     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

一种改进的中心-WENO混合格式

基于中心差分与WENO格式混合可以改善WENO格式耗散特性的思想,在理论推导的基础上,给出了一种用于激波捕捉计算的守恒型中心-WENO混合格式,该混合格式可视为三阶WENO格式和二阶中心差分格式的加权平均。在数值研究现有加权函数的基础上,给出了适用于该混合格式的加权函数,使其能够自适应地调整数值耗散以捕捉激波间断。数值结果表明:与三阶WENO格式相比,混合格式HY3_4能够降低数值耗散,更陡峭地捕捉间断,对复杂流场结构具有较高的分辨率;混合格式HY3_5对于包含高压比激波间断流场结构,能给出无振荡、低耗散的结果。     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全.针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见.开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据.采用直接数值模拟方法对来流马赫数6.0,34?压缩拐角内激波与湍流边界层的干扰问题进行了研究.基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性.通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制.采用动态模态分解方法,探讨了干扰流场的非定常运动历程.研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性.干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力-膨胀项,而对膨胀-耗散项影响较小.高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.     

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

反射激波冲击重气柱的RM不稳定性数值研究

数值研究了二维气柱在入射激波以及反射激波作用下的Richtmyer-Meshkov(RM)不稳定性发展规律, 采用有限体积法结合网格自 适应技术的VAS2D程序, 精确刻画激波和界面的演化. 入射平面激波的马赫数为1.2, 气柱界面内气体为六氟化硫(SF₆), 环境气体为空气, 激波管的尾端为固壁. 通过改变气柱与尾端之间的距离调节反射激波再次作用已经变形的气柱的时间, 获得不同时刻下已经变形的气柱形态、界面尺寸以及环量演化受到反射激波的影响. 结果表明, 反射激波再次作用气柱时, 气柱所处发展阶段不同, 界面演化规律以及环量随时间的变化也不相同, 反射激波与气柱相互作用过程中的涡量产生和分布与无反射情况差异较大, 揭示了不同情况下界面演化的物理机理.     

Numerical simulation of fast granular flow facing obstacles on steep terrains

The interaction between fast shallow granular flow and obstacles on steep terrain is an important aspect of granular mechanics and defending against geological hazards. In this study, we used a depth-averaged model for granular flow facing obstacles on steep terrains in a bed-fitted coordinate system where the obstacle system is treated as a local bed deviation term. A second-order Riemann-free scheme is extended to compute the depth-averaged model with a wetting–drying technique, which is verified by several granular flow cases, such as aluminum bar collapse and granular flow runout on a steep slope. Numerical simulations were performed for the case of granular flow facing a (i) single hemispherical obstacle and (ii) system of three hemispherical obstacles to produce a dynamical process and deposit profile, and show good agreement with experimental results. Granular flow facing a single obstacle on a concave plane produces a detached shock wave that moves upstream and a tailing rapid transition zone that moves down, which will merge to form a new shock for deposition. Granular flows facing a three-hemisphere obstacle system produce a tailing rapid transition zone that moves downstream and a downstream wavy shock that results from the interaction of three bow shocks in front of each obstacle. The downstream wavy shock moves upstream and merges with the upstream transition zone to form a new curved shock, which later relaxes to a deposit owing to bed friction. These findings provide some supplemental understandings of flow structures of fast granular flow facing obstacles.     

Variable density bore interaction with block obstacles

A horizontally variable density flow model is used to simulate hydraulic bore interactions with idealised urban obstacles. The 2D non-linear shallow water equations are solved using a second-order Monotonic Upstream-centered Schemes for Conservation Laws-Hancock Godunov-type HLLC approximate Riemann scheme. Validation test results are reported for wave propagation over a hump, a constant-density circular dam break and two 1D dam breaks involving different spatial distributions of solute concentration. Detailed parameter studies are then considered for hydraulic bore interactions with single and multiple-square obstacles under subcritical, critical and supercritical flow conditions. In all cases, reflected and diffracted wave patterns are generated immediately after the bore impacts the obstacle(s). Later, the incident bore reconstitutes itself downstream of the obstacle(s). Variable density flows are also considered, with the upstream volumetric concentrations set to values corresponding to water–sediment mixture densities of 1165 and 1495 kg/m³. It is found that the upstream Froude number, gap spacing between obstacles and upstream to downstream density difference influence the strength of the bore–structure interaction, run-up at the front face of the obstacle(s), and subsequent wave–wave interactions.     

Impingement of a supersonic underexpanded jet onto obstacles with different permeabilities

Results of an experimental study of the flow structure and characteristics of flow oscillations induced by impinging of a supersonic underexpanded jet onto obstacles made of permeable and impermeable materials are presented. Spectral characteristics of the acoustic field in the regime of unsteady interaction between the jet and the obstacle are determined. The influence of permeability of porous obstacles on these characteristics is identified.     

Active bogies and chassis levelling for a vehicle operating in rough terrain

Four axle vehicles with bogies can adapt the position of the wheels to follow irregularities in the terrain, having an obstacle surpassing ability far greater than conventional 2-axle vehicles. Still, the ability to overcome discrete obstacles on a steep slope is very different depending on the wheel that is facing the obstacle. A possible solution to diminish this variation can be found if the vehicle is able to actively redistribute the load on each wheel. One strategy is to design the suspension mechanism so it can regulate its height, being able to level the chassis. Also, an active torque on the pin join between the bogie and the chassis can be applied with the same goal, adopting a system of active bogies. Both solutions have been parametrically studied in a bi-dimensional multibody model of a 4-axle vehicle with double bogies. The results show an improvement independent of obstacle position and terrain angle when using active bogies. With height regulation, this improvement is limited to the rear bogie wheels, but the obstacle surmounting capacity of the vehicle as a whole can be considerably increased if the optimal regulation point is found. Possible applications for such enhanced vehicles with bogies are performing different tasks in forest areas with obstacles on steep slopes or unstructured terrain exploration.     

Free surface Stokes flows obstructed by multiple obstacles

Gravity‐driven Stokes flow down an inclined plane over and around multiple obstacles is considered. The flow problem is formulated in terms of a boundary integral equation and solved using the boundary element method. A Hermitian radial basis function (RBF) is used for the interpolation of the free surface, generation of the unit normal and curvature, and to prescribe the far‐field conditions. For flow over an obstacle, hemispheres are taken. For flow around an obstacle, circular cylinders are modelled and the contact angle condition on the obstacle/free surface intersection specified using the RBF formulation. Explicit profiles are produced for flow over and around two obstacles placed in various locations relative to one another. Interaction due to two obstacles is given by comparisons made with the profiles for flow over and around individual obstacles. In general, when the obstacles are separated by a sufficiently large distance the flow profiles are identical to a single obstacle analysis. For flow over and around two obstacles in‐line with the incident flow, effects of the governing parameters are examined, with variations in plane inclination angle, Bond number, obstacle size, and in the case of obstacles intersecting the free surface, static contact angle is considered. Finally flows over and around three obstacles are modelled. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of snow cover on obstacle performance of vehicles

Trafficability of terrain is a function of soft soil, hard or rough ground, geometric obstacles, vegetation, and the riverine environment. All of these terrain aspects are altered by cold temperatures and snow cover. This paper examines the effect of snow cover on obstacle crossing performance of vehicles. The mathematical expressions describing step negotiation, trench crossing, and slope climbing on snow-covered obstacles are given in terms of tracked vehicle, obstacle, and snow parameters.Tests of two tracked vehicles on snow-covered slopes, stream crossings, steps, and trenches were conducted, and some of the results were compared with computed values. Differences between computed and experimental values are attributed to neglecting slip-sinkage and track deflection in the computations.     

NUMERICAL STUDY ON PROPERTIES OF WEAK INTERLAYER OF LAMINATED COMPOSITE

<正>With assumption of material inhomogeneity on meso level,a three-point bending beam model has been adopted to simulate crack propagation in two dimensional laminated composite and to verify the toughening mechanism of energy dissipation and crack deflection along the weak interlayer.Moreover,the effect of strength,elastic modulus and thickness of the weak interlayer on both strength and toughness of the laminated composite have also been investigated in this paper.     

Dynamic and mass transfer characteristics of the flow issued from a bent chimney around buildings

The present work consists in an experimental investigation of the flow issuing from a bent chimney over a downstream obstacle. Our purpose is to explore the resulting flow field and its different characterizing features. These features were captured by means of the Particle Image Velocimetry technique. A numerical simulation of the problem has also been carried out and validated after comparison of the corresponding results to the experimental data. A good level of agreement was achieved between the experiments and the calculations. Then, we tried to upgrade our model by adopting large (real) scale dimensions. Our purpose consisted mainly in the observation and evaluation of the behavior of the incoming flow in presence of a double tandem obstacle. In a second step, we proposed to increase the number of the placed obstacles to four. The results given by the three-dimensional model are likely to highlight the dynamic features of the established field as well as the resulting mass transfer. Finally, we tried to evaluate the effect of further parameters on the characterizing features of the resulting flow filed such as the velocity ratio, the obstacles’ gap, the arrangement of the obstacles and the obstacles’ geometry.     

Towards the numerical simulation of the horizontal slug front

Steps towards the numerical simulation of the flow behind the slug front in horizontal slug flow performed with a streamfunction-vorticity representation of the mean flow and an energy dissipation model for the turbulence are discussed. The flow field consists of two vortices, one saddle point and four stagnation regions. Attention is focused on the following boundary conditions: moving wall jet, moving wall, free jet velocity discontinuity and vertical liquid-gas open surface. A dissipation flux boundary condition is suggested to simulate the interaction of the turbulent eddies with the open surface. A method to assess the necessity to use a transport model equation for the dissipation rather than a geometric specification of a length is suggested. Three different ways to characterize the mixing zone length are proposed.     

颗粒介质固-流态转变的理论分析及实验研究

颗粒介质由大量离散的颗粒聚集而成,因而与传统固体和流体不同,运动过程中的颗粒介质中可能同时存在多种流态及其相互间复杂的转换过程. 颗粒介质弹性失稳机理、不可恢复应变量化是研究颗粒介质固态和流态及固-流态转变的关键. 在前期建立的双颗粒温度热力学(two-granular-temperature, TGT) 理论基础上,确定了颗粒介质的弹性稳定性条件,建立了不可恢复应变流动法则,搭建了描述颗粒固态-液态及其相互转化的简单模型. 颗粒堆积体坍塌过程是典型的颗粒介质固态和流态及其转变过程,因此本文首先开展了25 167 个陶颗粒堆积体坍塌过程的实验研究,并使用基于TGT 理论的物质点方法和离散元方法对物理实验进行了模拟. 结果表明,模型数值结果与物理实验在颗粒堆坍塌过程中的形态、速度分布等细节上吻合很好,同时也发现了现阶段所使用的物质点方法和TGT 理论的不足. 初步说明TGT 理论可以实现颗粒介质固态和流态,以及状态转变的描述.