散心柱面胞格爆轰演化数值研究

采用有限体积方法,在自适应非结构网格上求解二维含化学反应Euler方程,数值研究了柱面胞格爆轰波演化现象．化学反应计算采用单步可逆总包反应模型．数值结果演示了散心柱面胞格爆轰波演化过程中胞格结构的分裂现象,获得了与实验结果定性一致的结果．胞格结构的分裂演化在点火区近场和远场显示了不同的特点,其中爆轰波传播过程中波阵面当地曲率的变化是控制胞格分裂演化行为的关键因素．数值结果也显示胞格结构的分裂现象来自于爆轰波前锋结构中横波的自组织行为,即沿爆轰波波面传播的小扰动发展成为横波的过程,这种现象与胞格爆轰波的不稳定性密切相关．     

炸药爆轰物质点法三维数值模拟

炸药爆轰以及多点起爆所产生的爆轰波汇聚问题很难应用有限元法进行模拟分析,尤其当网格发生畸变时,导致有限元法计算效率和数值精度严重下降,甚至无法得到正确结果.为此,该文应用显式积分算法的物质点法对炸药两点起爆和按时间序列的多点起爆的爆轰过程进行数值模拟,与炸药爆轰的理论计算结果相吻合.物质点法不但可以有效地避免网格畸变问题,而且为炸药爆轰的数值模拟提供了新的思路.     

膛口流场动力学过程数值研究

采用基于ALE方法的动网格及嵌入网格技术,运用有限体积方法,结合二阶精度Roe格式,对弹丸由高压气体驱动从静止状态加速至超音速,射出膛口到完全飞离初始流场的整个过程进行了数值模拟．根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程,揭示了在这一变化过程中激波与激波、激波与漩涡、激波与弹丸等的相互作用以及激波衍射、聚焦等对弹丸加速的影响．     

基于外场实验的风力机叶片三维效应研究

针对一台33 kW水平轴风电机组开展了外场实验,得到其叶片7个断面翼型的压力分布曲线;基于求解时均N-S方程对风轮进行三维数值模拟,以及将叶片各断面作为二维翼型进行数值计算,分别得到各断面翼型的压力分布曲线及升阻力系数．通过将外场实验、三维和二维数值计算所得压力分布曲线及升阻力系数进行对比分析,研究了三维效应对风力机气动性能的影响．研究表明,从叶尖到叶根各断面翼型的压差先增大后逐渐减小,叶片表面压力分布曲线比较明显地反映了从叶尖到叶根流动分离的变化;叶片表面压力分布的三维数值计算结果较二维计算结果更加接近于外场实验值;风力机叶片表面的三维流动对叶片的气动性能影响较大,在叶尖和叶根部分尤为突出.     

层流管流入口压力突然升高引起的水力瞬变的理论分析

研究了层流状态下管道入口压力突然升高引起的水力瞬变过程,建立了瞬态压力分布的偏微分方程和初边值条件,用分离变量法求得了压力的理论解．根据压力和流量间的约束关系,得到了关于流量的偏微分方程和初边值条件,用分离变量求得了瞬变过程流量分布理论解．最后,用特征线法(MOC)对该问题进行了数值求解,理论解和数值解吻合很好．     

锥形血管入口区域内管壁应力分析

本文对锥形血管入口区域的流动进行了探讨，导出了压力分布、轴向和径向的速度分布以及流场的切应力分布、管壁应力分布等公式，进行了相应的数值算例的研究和分析，还着重讨论了血管锥度角对管壁应力、压力分布等的影响。     

基于映射函数的三阶WENO改进格式及其应用

低耗散、高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义.在传统三阶WENO格式（WENO-JS3）和三阶WENO Z格式（WENO-Z3）基础上,基于映射函数,给出WENO-M3、WENO MZ3格式.选用Sod激波管、激波与熵波相互作用、双爆轰波碰撞及双Mach(马赫)反射等经典算例,考察上述格式的计算性能.数值结果表明,WENO-MZ3格式相较其他格式具有耗散低、对流场结构分辨率高的特性.为了进一步扩展WENO-MZ3格式的应用范围,采用该格式数值研究封闭方形舱室内柱形高压、高密度气体爆炸波传播过程,波系演化规律以及壁面典型测点压力载荷.数值计算结果表明WENO-MZ3格式能够较好地模拟包含高压比、高密度比的爆炸波且给出数值耗散较小的壁面压力载荷.     

隐式TVD格式在气液两相爆轰数值模拟中的应用

为研究FAE(燃料-空气炸药)爆炸参数规律,运用二维轴对称气液两相方程组模型,针对其中的气相方程组采用高分辨率的隐式TVD格式,液相方程组采用MacCormack格式,较好地对FAE气液两相爆轰产生的爆轰波的发展和传播过程进行了数值模拟,计算结果与国内外的研究结果符合良好.     

多翼裂缝压裂偏心井半解析模型及其渗流特征

考虑压裂多翼裂缝偏心井的实际情况,建立了多翼裂缝偏心井的数学模型.采用Laplace变换和压降叠加原理得到Laplace空间多翼裂缝压裂偏心井井底压力的半解析解.采用非均匀流量法,对井底压力的半解析解进行离散.结合Stehfest数值反演获得实空间井底压力的数值解和产量分布.借助SAPHIR试井分析软件建立了储层的数值试井模型并进行了数值离散计算.将计算结果与该文的半解析模型计算结果进行了对比,验证了该文模型的正确性.结果表明,多翼裂缝压裂偏心井井底压力变化可划分为8个主要流动阶段.最后讨论了裂缝的无因次导流能力、裂缝的不对称因子和井的偏心距对井底压力变化和产量分布特征的影响.     

轧板过程中单位压力的分布规律

轧板时单位压力的分布规律不仅用于计算总轧制力的数值,而且提供计算宽展、设计合理辊型的依据.然而迄今文献上所进行的分析,绝大部分只限于一维表达,不能反映单位压力沿变形区宽度上的变化.本文从变分法的原理出发,得出单位压力分布规律的二维表达式.     

High resolution multi-moment finite volume method for supersonic combustion on unstructured grids

In this study, we present a novel numerical model for simulating detonation waves on unstructured grids. In contrast to the conventional finite volume method (FVM), two types of moment comprising the volume-integrated average (VIA) and the point value (PV) at the cell vertex are treated as the evolution variables for the reacting Euler equations. The VIA is computed based on a finite volume formulation of the flux form where the conventional Riemann problem is solved by the HLLC Riemann solver. The PV is updated in a point-wise manner by using the differential formulation where the Roe solver is used to compute the differential Riemann problems. In order to increase the accuracy around discontinuities, numerical oscillations and dissipations are reduced using the boundary variation diminishing algorithm. Convergence tests demonstrated that the proposed model could achieve third-order accuracy with unstructured grids for reacting Euler equations. The high resolution property of the proposed method was verified based on simulations of several detonation wave propagation problems in two and three dimensions. In particular, the current model could resolve the cellular structures with fewer degrees of freedom for the unstable oblique detonation wave problem. These fine structures may be smoothed out by the conventional FVM due to the excessive amount of numerical dissipation errors. Importantly, a simulation of stiff detonation waves showed that the proposed method could capture the correct position of the reaction front whereas the conventional FVMs produced spurious phenomena. Thus, the proposed model can obtain highly accurate solutions for detonation problems on unstructured grids, which is highly advantageous for real applications involving complex geometrical configurations.     

陶瓷/金属复合靶板受可变形弹体撞击的理论分析模型

针对弹体撞击陶瓷/金属复合靶板的问题,将弹体的变形、陶瓷面板的碎裂和金属背板的变形结合起来,建立了新的可变形弹体垂直撞击陶瓷/金属靶板的理论分析模型．模型中计入了弹体刚性区长度和运动速度、塑性变形区长度、横截面积和运动速度的变化以及弹体对靶板的侵入速度和深度;对陶瓷面板考虑了陶瓷锥体积和抗压强度的变化;对金属背板的变形,根据其塑性变形功、外力功及其动能守恒原理,得到金属背板的运动方程．最后对具体算例进行了分析,得到了各物理量随时间的变化,给出了一些有价值的规律．结果表明,模型能较好地描述撞击过程中的有关规律;与实验结果和数值模拟结果进行对比,吻合较好,说明了模型的有效性．     

装药驱动飞片引爆炸药性能影响参数分析

扩爆装药结构对爆轰波传播、飞片驱动过程以及对主炸药引爆性能有直接的影响。为分析装药结构对飞片威力参数的影响,针对装药直径、飞片厚度、飞片拱起高度等主要结构参数,利用正交实验原理设计了数值实验方案,并采用动力学有限差分程序建立了相应的数值模拟模型。通过对数值实验结果的对比和统计分析,获得了影响飞片速度、动量、比动能等引爆炸药威力指标的主要装药结构参数及其影响规律。其结果可为相关扩爆装药设计提供理论依据。     

Numerical Simulation of Detonation Excitation at Shock Wave Interaction with Dust Cloud in Channel

The study of detonation ability of reactive particle gas mixtures is necessary to prevent industrial explosions in industries where dispersed powders are used. The present paper focuses on numerical simulation of the shock wave interaction with semiinfinite aluminum dust cloud, which is situated inside a plane channel. The cloud fills entirely or partly the channel cross‐section and has initially a rectangular shape. The scenarios of detonation initiation in the cloud are determined depending on the incident shock wave amplitude values. The processes of transformation and spreading of finite width clouds under weak incident shock wave action (when the particles do not ignite) are investigated. The types of an oblique shock wave reflection from the plane of symmetry in the cloud are analyzed. The processes of particle ignition and detonation structure formation at strong incident shock wave action are investigated. Nonstationary periodic fuctuations take place in the detonation flow due to transversal wave effect. Nevertheless the detonation structure established propagates in quasistationary regime. If the incident shock wave is attenuated with a rarefaction wave then the detonation formation fails at clouds of insufficient width.     

爆炸冲击载荷下航空铝合金平板动态响应数值分析方法

为得到适用于爆炸冲击载荷下航空铝合金平板动态响应的数值分析方法,采用LS-DYNA显式动力学分析软件对爆炸冲击载荷下的铝合金平板进行数值仿真计算.主要研究了不同的任意Lagrange-Euler(拉格朗日-欧拉)网格（ALE）输运步算法、流固耦合方式、流固耦合点数量、网格尺寸、有限元单元类型对计算结果的影响.通过计算结果与实验结果的分析对比,表明采用van Leer+HIS输运步算法、罚函数耦合方式、在流体网格与结构网格之间采用3个耦合点、结构网格尺寸与空气域网格尺寸比例设为2∶1、结构单元采用163号壳单元时可以较为准确地计算航空铝合金平板在爆炸冲击载荷下的动态响应,并且能提高计算效率,节约计算时间.     

An analytical solution for conduction-dominated unidirectional solidification of binary mixtures

In this paper, a fully analytical solution technique is established for the solution of unidirectional, conduction-dominated, alloy solidification problems. By devising appropriate averaging techniques for temperature and phase-fraction gradients, governing equations inside the mushy region are made inherently homogeneous. The above formulation enables one to obtain complete analytical solutions for solid, liquid and mushy regions, without resorting to any numerical iterative procedure. Due considerations are given to account for variable properties and different microscopic models of alloy solidification (namely, equilibrium and non-equilibrium models) in the two-phase domain. The results are tested for the problem of solidification of a NH₄Cl–H₂O solution, and compared with those from existing analytical models as well as with the corresponding results from a fully numerical simulation. The effects of different microscopic models on solidification behaviour are illustrated, and transients in temperature and heat flux distribution are also analysed. A good agreement between the present solutions and results from computational simulation is observed.     

Numerical and empirical approach in predicting the penetration of a concrete target by an ogive-nosed projectile

This paper demonstrates the application of both numerical simulation and empirical equation in predicting the penetration of a concrete target by an ogive-nosed projectile. The results from the experiment performed by Gran and Frew [In-target radial stress measurements from penetration experiments into concrete by ogive-nose steel projectiles, Int. J. Impact Eng. 19 (8) (1997) 715–726] are used as a benchmark for comparison. In the numerical simulations a 3.0-caliber radius-head steel ogival-nose projectile with a mass of 2.3 kg is fired against cylindrical concrete target with a striking velocity of 315 m/s. The simulation, performed using AUTODYN 2-D, assesses three numerical schemes, namely Langrange, Euler–Lagrange coupling and smooth particles hydrodynamics SPH–Lagrange coupling, in predicting the maximum depth of penetration and the radial stress–time response of the concrete target. When assessing the three solution techniques we hypothesize that the effect of strain rate on strength for the concrete target does not adversely affect the prediction on the maximum depth of penetration and the radial stress–time response of the concrete target. In the empirical approach the penetration equation developed by Forrestal et al. [An empirical equation for penetration depth of ogive-nose projectiles into concrete targets, Int. J. Impact Eng. 15 (4) (1994) 395–405] is used to determine the maximum depth of penetration and the deceleration–time response. The deceleration–time response for the projectile using the empirical approach is compared with those obtained from the numerical simulations. Results from both the numerical and empirical approaches are consistent. The calculated depth of penetration from both approaches yield relatively good agreement with that obtained from the experiment. The numerical simulations using each of the three numerical schemes are also able to reproduce the profiles from the radial stress measurements. Simulations using the SPH numerical scheme give the best overall agreement. The good overall agreement with the experimental radial stress measurements and consistent results between both empirical and numerical approach, enhanced the confidence in engineers and ballisticians when using these two approaches in complementing full-scale testing.     

DEPENDENCE OF QUALITATIVE BEHAVIOR OF THE NUMERICAL SOLUTIONS ON THE IGNITION TEMPERATURE FOR A COMBUSTION MODEL

We study the dependence of qualitative behavior of the numerical solutions （obtained by a projective and upwind finite difference scheme） on the ignition temperature for a combustion model problem with general initial condition. Convergence to weak solution is proved under the Courant-Friedrichs-Lewy condition. Some condition on the ignition temperature is given to guarantee the solution containing a strong detonation wave or a weak detonation wave. Finally, we give some numerical examples which show that a strong detonation wave can be transformed to a weak detonation wave under some well-chosen ignition temperature.