Computational modelling of multi-material energetic materials and systems

ABSTRACT

In this paper, we present a systematic roadmap for developing a robust and parallel multi-material reactive hydrodynamic solver that integrates historically stable algorithms with new and current modern methods to solve explosive system design problems. The Ghost Fluid Method and Riemann solvers were used to enforce appropriate interface boundary conditions. Improved performance in terms of computational work and convergence properties was achieved by modifying a local node sorting strategy that decouples ghost nodes, allowing us to set material boundary conditions via an explicit procedure, removing the need to solve a coupled system of equations numerically. The locality and explicit nature of the node sorting concept allows for greater levels of parallelism and lower computational cost when populating ghost nodes. Non-linear numerical issues endemic to the use of real Equations of State in hydro-codes were resolved by using more thermodynamically consistent forms allowing us to accurately resolve large density gradients associated with high energy detonation problems at material interfaces. Pre-computed volume tables were implemented adding to the robustness of the solver base.     

生产率差异与企业对外直接投资进入模式选择

基于企业异质性视角,借鉴演化博弈理论探究企业通过OFDI“走出去”时的跨国并购方和标的方在全要素生产率上的差异对企业OFDI进入策略选择的影响机制。利用Hotelling模型构造资本密集型企业与劳动密集型企业OFDI进入策略选择博弈模型,讨论两类企业OFDI进入方式决策在并购双方生产率差异变化之下的变迁机理。研究发现:当并购方和被并购方生产率差异较小时,市场达到两类企业都倾向于跨国并购策略的演化稳定状态;当并购双方生产率差异较大时,企业的生产要素密集度及其结构对其OFDI进入模式决策具有重要的作用,知识或技术密集度较高的资本密集型企业与综合实力强势的部分劳动密集型企业采取跨国并购战略,而另一部分劳动密集型企业选择绿地新建方式进入国际市场。即企业OFDI进入策略选择因其所属的行业生产要素密集度及其结构不同而具有差异性。研究结论在一定程度上弥补了现有研究的不足,为以后的研究工作提供一定的理论参考。     

CH4/空气简化动力学机理数值验证

选择绕圆柱预混燃烧算例，验证CH₄/空气三种简化动力学机理（16s41r、15s19r和53s325r）.考虑均匀来流，忽略湍流和湍流与燃烧相互作用以及燃料扩散效应，假设层流有限反应速率，采用保自由流5阶WENO格式求解多组分Euler方程组，得到CH₄/空气预混燃烧流场温度等值线、沿驻点线压力和温度及其CH₄、CO和CO₂质量百分数分布.结果表明：三种简化动力学机理给出的流场均出现弓形激波和火焰面，弓形激波和火焰驻点距离及其形状、诱导区宽度和简化动力学机理相关.当圆柱直径增大，弓形激波和火焰向圆柱上游移动，对应的驻点距离均增大，诱导区宽度变短，点火延时变小，但火焰和弓形激波位置次序未变化.53s325r模型要比16s41r模型和15s19r模型精度要高，点火延时覆盖的压力和温度范围也较宽，所有简化机理均未完全反应，在较大圆柱直径下游达到化学平衡.     

An improved generalized predictive control algorithm based on the difference equation CARIMA model for the SISO system with known strong interference

ABSTRACT

The single input single output (SISO) system with known strong interference is widely used in various occasions. Due to its strong interference, the control accuracy is hard to guarantee. To solve this problem, an improved generalized predictive control (IGPC) algorithm is developed. The IGPC firstly builds the difference equation CARIMA (Controlled Auto-Regressive Integrated Moving-Average) model of the SISO system and then treats the system as a two input single output (TISO) system and calculates its predictive vector, then transforms it into a SISO system and uses the TISO system predictive vector to calculate the SISO system control increment. A new parameter called phase coefficient is added to inhibit the control lag. Simulations are performed to make the comparison among the traditional GPC, PID control, velocity synchronization control (VSC), fuzzy adaptive PID control (FAPID), model-based robust PID control (BPID) and the IGPC. Results show that IGPC has best performance compared to the others. Finally, experiments are developed which proved that the IGPC algorithm has a higher accuracy in the SISO system with known strong interference than that of VSC.     

A three‐dimensional Cartesian cut cell method for incompressible viscous flow with irregular domains

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L₁ and L _∞ norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

Implementation of ADER scheme for a bore on an unsaturated permeable slope

The paper details the implementation of the Godunov‐type finite volume Arbitrary high order schemes using Derivatives (ADER) scheme for the case of a large source term in the continuity equation of the nonlinear shallow water equations. The particular application is the movement of a bore on a highly permeable slope. The large source term is caused by the infiltration into the initially unsaturated slope material. Infiltration is modelled as vertical downwards piston‐like flow with Forchheimer quadratic parameterisation of the resistance law. The corresponding ODE is solved using the fourth‐order Runge–Kutta method. The surface and subsurface flow models have been tested by comparison with analytical solutions. Example predictions of surface bore propagation and wetting front propagation are presented for a range of slope permeabilities. The effects of permeability on bore run‐up, water depths and velocities are illustrated. The ADER scheme is capable of handling the source term, including the extreme case when this term dominates the volume balance. Copyright © 2011 John Wiley & Sons, Ltd.     

Unsteady hydromagnetic flow of a reactive variable viscosity third‐grade fluid in a channel with convective cooling

This paper examines the combined effects of a transverse magnetic field and variable viscosity on unsteady flow of a reactive third‐grade electrically conducting fluid and heat transfer in a channel with convective cooling at the surface. It is assumed that the fluid has small electrical conductivity and the electromagnetic force produced is very small. The coupled nonlinear partial differential equations governing the problem are derived and solved numerically using a semi‐implicit finite‐difference scheme. Both numerical and graphical results are presented and physical aspects of the problem are discussed with respect to various parameters embedded in the system. It is in general noted that those parameters that increase/decrase one flow quantity (velocity or temperature) also lead to the increase/decrease respectively of the other quantity. Copyright © 2011 John Wiley & Sons, Ltd.     

Combining GPS carrier phase and Doppler observations for precise velocity determination

The truncation error and propagation error are analyzed for velocity determination through differential GPS carrier phase observations,and an approach for the choice of the best number of points for the central difference method is developed.In order to overcome the disadvantages of existing GPS velocity determination methods,a new velocity determination algorithm is presented,based on combining carrier phase and Doppler observations.The basic idea is that two types of observation are combined by adding the...     

Approximant crystals of icosahedral Mg–(Ga,Al)–Zn alloys

The zero field cooled (ZFC) and field cooled (FC) low-field magnetic moment m of a dense frozen ferrofluid containing Fe₅₅Co₄₅ particles of size 4.6nm in hexane exhibits irreversibility at temperatures T?T _b≈ 30?K. FC in μ ₀ H ≤ 1?T gives rise to shifted minor hysteresis loops below T _b. At T _c≈ 10?K, sharp peaks of m ^ZFC and of the ac susceptibility χ ′, a kink of the thermoremanent magnetic moment m ^TRM, a sizeable reduction of the coercive field H _c, and the appearance of a spontaneous moment m ^SFM indicate a phase transition with near mean-field critical behaviour of both m ^SFM and χ ′ . These features are explained within a core-shell model of nanoparticles, whose strongly disordered shells gradually become blocked below T _b, while their soft ferromagnetic cores couple dipolarly and become superferromagnetic (SFM) below T _c.