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.     

Efficient finite-difference multi-scheme approach to the simulation of seismic waves in anisotropic media

This paper presents an original multi-scheme approach to the numerical simulation of seismic wave propagation in models with anisotropic formations. To simulate wave propagation in the anisotropic parts of the model, the Lebedev scheme is used. This scheme is rather universal, but highly expensive in terms of computational efficiency. In the main part of the model, a highly efficient standard staggered grid scheme is proposed. The two schemes are coupled to ensure convergence of the reflection/propagation coefficients with a prescribed order. The algorithm combines the universality of the Lebedev scheme and the efficiency of the standard staggered grid scheme.     

On the application of the nonstationary form of the Tappert equation as an artificial boundary condition

In this paper, a nonstationary analog of the range refraction parabolic equation is derived. A new approach to the derivation of Tappert’s operator asymptotic formula with the use of noncommutative analysis is presented. The obtained nonstationary equation is proposed as an artificial boundary condition for the wave equation in underwater acoustics. This form of artificial boundary condition has low computational cost and systematically takes into account variations of sound speed. This is confirmed by various numerical experiments, including propagation of normal modes and wave fields produced by point source.     

A parallel adaptive grid algorithm for computational shock hydrodynamics

Adaptive Mesh Refinement (AMR) algorithms that dynamically match the local resolution of the computational grid to the numerical solution being sought have emerged as powerful tools for solving problems that contain disparate physical scales. In particular several workers have demonstrated the effectiveness of employing an adaptive, hierarchical block-structured grid system for time-accurate simulations of complex shock wave phenomena. Here we present an overview of one such block-structured AMR algorithm. Our formulation has progressed far beyond the development stage to become a reliable numerical tool for performing detailed investigations of complex flows. While our refinement machinery is not tied to a specific application and is intended for general use, in this paper we adopt detonation phenomena as a theme so as to provide a sense of purpose.     

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.     

流体饱和多孔隙介质波动方程小波有限差分法

研究流体饱和多孔隙介质中波动方程的数值模拟．针对求解二维弹性波方程问题,提出小波有限差分法．该方法综合了小波多分辨分析计算灵活、计算效率高特性和有限差分易于实现的优点．数值模拟的结果显示,此方法对于求解流体饱和多孔隙介质方程的数值模拟是有效稳定的．     

Modulational instability in a full‐dispersion shallow water model

We propose a shallow water model that combines the dispersion relation of water waves and Boussinesq equations, and that extends the Whitham equation to permit bidirectional propagation. We show that its sufficiently small and periodic traveling wave is spectrally unstable to long wavelength perturbations if the wave number is greater than a critical value, like the Benjamin‐Feir instability of a Stokes wave. We verify that the associated linear operator possesses infinitely many collisions of purely imaginary eigenvalues, but they do not contribute to instability to the leading order in the amplitude parameter. We discuss the effects of surface tension. The results agree with those from a formal asymptotic expansion and a numerical computation for the physical problem.     

Numerical simulation of spinning detonation in circular section channels

Numerical simulation of three-dimensional structures of gas detonation in circular section channels that emerge due to the instability when the one-dimensional flow is initiated by energy supply at the closed end of the channel is performed. It is found that in channels with a large diameter, an irregular three-dimensional cellular detonation structure is formed. Furthermore, it is found that in channels with a small diameter circular section, the initially plane detonation wave is spontaneously transformed into a spinning detonation wave, while passing through four phases. A critical value of the channel diameter that divides the regimes with the three-dimensional cellular detonation and spinning detonation is determined. The stability of the spinning detonation wave under perturbations occurring when the wave passes into a channel with a greater (a smaller) diameter is investigated. It is found that the spin is preserved if the diameter of the next channel (into which the wave passes) is smaller (respectively, greater) than a certain critical value. The computations were performed on the Lomonosov supercomputer using from 0.1 to 10 billions of computational cells. All the computations of the cellular and spinning detonation were performed in the whole long three-dimensional channel (up to 1 m long) rather than only in its part containing the detonation wave; this made it possible to adequately simulate and investigate the features of the transformation of the detonation structure in the process of its propagation.     

Numerical modeling of seismic and acoustic-gravity waves propagation in an “Earth-Atmosphere” model in the presence of wind in the air

In this paper, an effective numerical algorithm for 2.5D seismic and acoustic-gravitational wave propagation is applied to a combined “Earth-Atmosphere” model in the presence of wind in the air. Seismic wave propagation in an elastic half-space is described by a system of first-order dynamic equations of elasticity theory. The propagation of acoustic-gravitational waves in the atmosphere in the presence of wind is described by the linearized Navier-Stokes equations. The algorithm is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform with respect to a spatial coordinate combined with a finite difference method for the reduced problem.     

On the local fractional wave equation in fractal strings

The key aim of the present study is to attain nondifferentiable solutions of extended wave equation by making use of a local fractional derivative describing fractal strings by applying local fractional homotopy perturbation Laplace transform scheme. The convergence and uniqueness of the obtained solution by using suggested scheme is also examined. To determine the computational efficiency of offered scheme, some numerical examples are discussed. The results extracted with the aid of this technique verify that the suggested algorithm is suitable to execute, and numerical computational work is very interesting.     

Properties of Detonation Waves

The phenomenon of unphysical wave propagation speeds sometimes occurs in numerical computations of detonation waves on coarse grids. The strong detonation wave splits into two parts, a weak detonation which travels with the speed of one cell per time step and an ordinary shock wave. We analyse a simplified set of equations and look for travelling wave solutions. It is shown that the solution depends on the dimensionless number Kr = μK/Qρ₁. Here μ is the viscosity, K is the rate of reaction, Q is the heat release available in the process and ρ₁ is the density at the unburnt state. It is shown that the density peak of the travelling wave depends on Kr and also, that if Kr is sufficiently large there is no travelling wave solution. The erroneous behaviour above is explained as an effect of the artificial viscosity necessarily inherent in the numerical methods when coarse grids are used. To prevent this unphysical behaviour we suggest the use of an ‘artificial rate of reaction’ such that the actual value of Kr used in the numerical method retains its correct physical value.     

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

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

A quadratic b-spline based isogeometric analysis of transient wave propagation problems with implicit time integration method

A uniform quadratic b-spline isogeometric element is exclusively considered for wave propagation problem with the use of desirable implicit time integration scheme. A generalized numerical algorithm is proposed for dispersion analysis of one-dimensional (1-D) and two-dimensional (2-D) wave propagation problems where the quantified influence of the defined CFL number on wave velocity error is analyzed and obtained. Meanwhile, the optimal CFL (Courant–Friedrichs–Lewy) number for the proposed 1-D and 2-D problems is suggested. Four representative numerical simulations confirm the effectiveness of the proposed method and the correctness of dispersion analysis when appropriate spatial element size and time increment are adopted. The desirable computation efficiency of the proposed isogeometric method was confirmed by conducting time cost and calculation accuracy analysis of a 2-D numerical example where the referred FEM was also tested for comparison.     

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

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

Control of space-time chaos in a system of equations of the FitzHugh–Nagumo type

We perform an analytic and numerical study of a system of partial differential equations that describes the propagation of nerve impulses in the heart muscle. We show that, for fixed parameter values, the system has infinitely many distinct stable wave solutions running along the spatial axis at arbitrary velocities and infinitely many distinct modes of space-time chaos, where the bifurcation parameter is the velocity of running wave propagation along the spatial axis, which does not explicitly occur in the original system of equations. We suggest an algorithm for controlling the space-time chaos in the system, which permits one to stabilize any of its unstable periodic running waves.     

A mathematical description of the acoustic coupling of the mass/spring model

This paper describes hybrid mathematical model which couples the mechanics of the mass/spring model to the acoustic wave propagation model for use in generating the acoustic signal emitted by complex structures of paper fibres under strain. A discussion of the coupling method is presented including remarks on the errors encountered intrinsic to the discretisation scheme. The numerical results of a vibrating rubber band and a vibrating paper fibre are compared to their experimental counterparts. The fundamental frequencies of the acoustic signals are compared showing a close agreement between the experimental and numerical results.     

Numerical identification of the leading coefficient of a parabolic equation

For a multidimensional parabolic equation, we study the problem of finding the leading coefficient, which is assumed to depend only on time, on the basis of additional information about the solution at an interior point of the computational domain. For the approximate solution of the nonlinear inverse problem, we construct linearized approximations in time with the use of ordinary finite-element approximations with respect to space. The numerical algorithm is based on a special decomposition of the approximate solution for which the transition to the next time level is carried out by solving two standard elliptic problems. The capabilities of the suggested numerical algorithm are illustrated by the results of numerical solution of a model inverse two-dimensional problem.     

Front propagation for discrete periodic monostable equations

This paper deals with front propagation for discrete periodic monostable equations. We show that there is a minimal wave speed such that a pulsating traveling front solution exists if and only if the wave speed is above this minimal speed. Moreover, in comparing with the continuous case, we prove the convergence of discretized minimal wave speeds to the continuous minimal wave speed.     

Dynamics of coal dust flame acceleration: A feedback control model for unsteady flow

Several modeling concepts borrowed from control theory are employed to develop an algebraic and ordinary differential equations model for the dynamics of unsteady coal dust flame acceleration in a constant area duct closed at one end, e.g., in a coal mine tunnel. We are particularly concerned with modeling the feedback mechanisms which cause a coal dust flame to accelerate, leading to detonation. Previous experimental studies have been conducted on both coal dust flame propagation and on individual coal particle combustion. Based on the results, a physical model is proposed in which coal dust flame acceleration is entirely controlled, in a feedback fashion, by volatiles emission and their reaction. A control system model is developed that employs five well-stirred reactor subsystems with three feedback interaction mechanisms. The model consists of a leading shock wave, followed by a variable length volatiles emission region ahead of the flame, a fixed length burning region immediately behind the flame front, and a variable length exhaust region extending back to the closed end of the duct. The feedback mechanisms incorporated into the model include heat transfer and pressurization from the burning region to the volatiles emission region, and pressurization from the volatiles emission region to the turbulent mixing region behind the shock wave. Each well-stirred reactor is described by a system of algebraic and ordinary differential equations for the rate of change of conditions inside the reactor. Numerical simulation results reveal that, despite far-reaching simplifications (ordinary instead of partial differential equations, ideal gases insteady of two-phase flow, separation of volatiles emission and combustion, neglection of char burning), the model exhibits the fundamental dynamic properties of the flame propagation process. The model agrees with qualitative photographic experimental results and is applicable to both the case where the flame accelerates to detonation and to the case where the combustion process dies out.     

直线爆炸所形成的柱面波的近似分析解

由直线爆炸所形成的柱面波,文献[1]得到了数值解。本文对于绝热指数γ=C_p/C_v较大的情形,取ε=1/γ2为小参数,用伸缩坐标法得到了问题的一级近似分析解。算例(γ=3)表明,近似解与准确解(数值解)十分符合。