Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

Diffraction of elastic waves in the plane multiply-connected region and dynamic stress concentration

This paper deals with the problem of diffraction of elastic waves in the plane multiply-connected regions by the theory of complex functions. The complete function series which approach the solution of the problem and general expressions for boundary conditions are given.’ Then the problem is reduced to the solution to infinite series of algebraic equations and the solution can be directly obtained by using electronic computer. In particular, for the case of weak interaction, an asymptotic method is presented here, by which the problem ofp waves diffracted by a circular cavities is discussed in detail. Based on the solution of the diffracted wave field the general formulas for calculating dynamic stress concentration factor for a cavity of arbitrary shape in multiply-connected region are given.     

Weak shock waves and shear bands in compressible,inextensible thermoelastic solids

A study of weak shock waves propagating into a solid, which is compressible but temperature-dependent extensible in a specified direction is presented. The inextensible solid is also considered. The constitutive equations of constrained thermoelastic material are written as the summation of constrained and unconstrained counterparts of the relevant quantities. The equation of motion of weak shock waves, which is recovered by the theory of singular surfaces, reduces to an eigenvalue problem. The solution of this eigenvalue problem yields the speeds of propagation of weak shock waves. In the case of an undeformed solid, the speeds of these waves are explicitly expressed. Additionally, a discussion on the ductility limits of constrained thermoelastic material subjected to the uniaxial and biaxial extensions is presented.     

Regular reflection of a magnetohydrodynamic shock wave from a conducting wall

In two-dimensional supersonic gasdynamics, one of the classical steady-state problems, which include shock waves and other discontinuities, is the problem concerning the oblique reflection of a shock wave from a plane wall. It is well known [1–3] that two types of reflection are possible: regular and Mach. The problem concerning the regular reflection of a magnetohydrodynamic shock wave from an infinitely conducting plane wall is considered here within the scope of ideal magnetohydrodynamics [4]. It is supposed that the magnetic field, normal to the wall, is not equal to zero. The solution of the problem is constructed for incident waves of different types (fast and slow). It is found that, depending on the initial data, the solution can have a qualitatively different nature. In contrast from gasdynamics, the incident wave is reflected in the form of two waves, which can be centered rarefaction waves. A similar problem for the special case of the magnetic field parallel to the flow was considered earlier in [5, 6]. The normal component of the magnetic field at the wall was equated to zero, the solution was constructed only for the case of incidence of a fast shock wave, and the flow pattern is similar in form to that of gasdynamics. The solution of the problem concerning the reflection of a shock wave constructed in this paper is necessary for the interpretation of experiments in shock tubes [7–10].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 102–109, May–June, 1977.The author thanks A. A. Barmin, A. G. Kulikovskii, and G. A. Lyubimov for useful discussion of the results obtained.     

The structure and asymptotic behaviour of compression waves

In the study of weak solutions to nonlinear hyperbolic partial differential equations both rarefaction waves and compression waves arise. Although the behavior of rarefaction waves is known for all time, the characteristics that determine a compression wave intersect and hence the development of the wave is not easily determined. The purpose of this paper is to study compression waves. As a first step we consider the Cauchy problem for the nonlinear wave equation. We show that if the data outside some finite interval consist of constant states, then after finite time the solution involves the same states as does the solution to the Riemann problem determined by these constant states. This result is then applied to compression waves to obtain information on the shock that arises and on the steady-state solution. The region of interaction is also described. This information is obtained via a constructive procedure.     

Evolution of weak shocks in one dimensional planar and non-planar gasdynamic flows

Asymptotic decay laws for planar and non-planar shock waves and the first order associated discontinuities that catch up with the shock from behind are obtained using four different approximation methods. The singular surface theory is used to derive a pair of transport equations for the shock strength and the associated first order discontinuity, which represents the effect of precursor disturbances that overtake the shock from behind. The asymptotic behaviour of both the discontinuities is completely analysed. It is noticed that the decay of a first order discontinuity is much faster than the decay of the shock; indeed, if the amplitude of the accompanying discontinuity is small then the shock decays faster as compared to the case when the amplitude of the first order discontinuity is finite (not necessarily small). It is shown that for a weak shock, the precursor disturbance evolves like an acceleration wave at the leading order. We show that the asymptotic decay laws for weak shocks and the accompanying first order discontinuity are exactly the ones obtained by using the theory of non-linear geometrical optics, the theory of simple waves using Riemann invariants, and the theory of relatively undistorted waves. It follows that the relatively undistorted wave approximation is a consequence of the simple wave formalism using Riemann invariants.     

非线性岩土中球壳对冲击波的动力响应分析

采用小参数摄动法和Ｌａｐｌａｃｅ变换研究了非线性弹性岩土中球壳对冲击波的动力响应问题。系统的非线性方程由小参数摄动渐近展开后，利用Ｓｔｏｋｅｓ－Ｈｅｌｍｈｏｌｔｚ矢量分解定理把它们简化为一系列的线性波动方程，并由Ｌａｐｌａｃｅ变换和本征函数给出了各线性波动方程的求解。最后，对平面冲击波和球面冲击波给出了球壳动力响应的应力结果。     

Nonlinear dispersion of long internal waves

The propagation of long weakly nonlinear waves in an atmospheric waveguide is considered. A model system of Kadomtsev-Petviashvili equations [1], which describes the propagation of such waves, is derived. In the case of one excited wave mode the system of model equations goes over into the Kadomtsev-Petviashvili equation, in which, however, the variables x and t are interchanged. The reasons for this are clarified. In the two-dimensional case an approximate solution of the model equations is constructed, and steady nonlinear waves and their interaction in a collision are considered. The results of a numerical verification of the stability of the approximate steady solutions and of the solution to the problem of decay of the wave into quasisolitons are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 151–157, May–June, 1988.     

Three-dimensional propagation of the transmitted shock wave in a square cross-sectional chamber

An investigation into the three-dimensional propagation of the transmitted shock wave in a square cross-section chamber was described in this paper, and the work was carried out numerically by solving the Euler equations with a dispersion-controlled scheme. Computational images were constructed from the density distribution of the transmitted shock wave discharging from the open end of the square shock tube and compared directly with holographic interferograms available for CFD validation. Two cases of the transmitted shock wave propagating at different Mach numbers in the same geometry were simulated. A special shock reflection system near the corner of the square cross-section chamber was observed, consisting of four shock waves: the transmitted shock wave, two reflection shock waves and a Mach stem. A contact surface may appear in the four-shock system when the transmitted shock wave becomes stronger. Both the secondary shock wave and the primary vortex loop are three-dimensional in the present case due to the non-uniform flow expansion behind the transmitted shock.PACS: 43.40.Nm     

Solution of Riemann problem for dusty gas flow

A direct approach is used to solve the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady planar flow of an isentropic, inviscid compressible fluid in the presence of dust particles. The elementary wave solutions of the Riemann problem, that is, shock waves, rarefaction waves and contact discontinuities are derived and their properties are discussed for a dusty gas. The generalised Riemann invariants are used to find the solution between rarefaction wave and the contact discontinuity and also inside rarefaction fan. Unlike the ordinary gasdynamic case, the solution inside the rarefaction waves in dusty gas cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. Although the case of dusty gas is more complex than the ordinary gas dynamics case, all the parallel results for compressive waves remain identical. We also compare/contrast the nature of the solution in an ordinary gasdynamics and the dusty gas flow case.     

Theory of diffraction of weak shock waves

A numerical solution is considered to the universal nonlinear boundary-value diffraction problem which occurs in various problems of weak interaction [1, 2] in the asymptotic analysis of the flow in a region with large gradients of the parameters near the point of intersection of the incident, diffracted, and reflected waves. The analytical solutions to this type of problem usually approximately satisfy the conditions on the diffracted front, the position of which is not known beforehand, but is found along with the solution. In the present paper, the problem is solved by the numerical method of [3], which reduces the initial boundary-value problem for the system of short-wave equations with an unknown boundary to the solution of a series of boundary-value problems with a fixed boundary. The problem of the diffraction of a weak shock wave on a wedge with a finite apex angle is considered as an application of the solution. The data calculated by the asymptotic theory agree significantly better with the experimental data [5] than the theoretical data of [4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 176–178, November–December, 1984.     

Self-Switching of Waves in Materials

The problem is formulated on self-switching of acoustic waves in materials when the wavelength does not exceed considerably the representative dimension of the microstructure. Use is made of the microstructural theory of two-component mixtures, which is adapted well to study waves in composites. Truncated and evolutionary equations describing the interaction of two longitudinal plane waves — a strong pumping wave and a weak signal wave — are derived. A procedure of finding the exact solution to the evolutionary equations is described for the case where the wave numbers of the pumping and signal waves are equal. The solution is expressed in terms of the elliptic Jacobi function     

Growth and decay of weak shock waves in magnetogasdynamics

The purpose of the present study is to investigate the problem of the propagation of weak shock waves in an inviscid, electrically conducting fluid under the influence of a magnetic field. The analysis assumes the following two cases: (1) a planar flow with a uniform transverse magnetic field and (2) cylindrically symmetric flow with a uniform axial or varying azimuthal magnetic field. A system of two coupled nonlinear transport equations, governing the strength of a shock wave and the first-order discontinuity induced behind it, are derived that admit a solution that agrees with the classical decay laws for a weak shock. An analytic expression for the determination of the shock formation distance is obtained. How the magnetic field strength, whether axial or azimuthal, influences the shock formation is also assessed.     

Behaviour of the acceleration and shock waves in materials with internal state variables

The explicit expressions for the change in the amplitudes of one-dimensional acceleration and shock waves propagating through arbitrary homogeneous materials described by the strain and internal state variables/parameters/are derived. The existence of a critical amplitude β for the acceleration wave and a critical strain gradient λ for the shock wave is established. For an infinitesimal shock wave the general form of the solution of the governing differential equation is furnished. The differential equations for the amplitudes of these two kind of waves are applied to an elastic-viscoplastic material.     

General Relativistic Shock Waves that Extend the Oppenheimer-Snyder Model

In earlier work we constructed a class of spherically symmetric, fluid dynamical shock waves that satisfy the Einstein equations of general relativity. These shock waves extend the celebrated Oppenheimer-Snyder result to the case of non-zero pressure. Our shock waves are determined by a system of ordinary differential equations that describe the matching of a Friedmann-Robertson-Walker metric (a cosmological model for the expanding universe) to an Oppenheimer-Tolman metric (a model for the interior of a star) across a shock interface. In this paper we derive an alternate version of these ordinary differential equations, which are used to demonstrate that our theory generates a large class of physically meaningful (Lax-admissible) outgoing shock waves that model blast waves in a general relativistic setting. We also obtain formulas for the shock speed and other important quantities that evolve according to the equations. The resulting formulas are important for the numerical simulation of these solutions. (Accepted January 19, 1996)     

Solution to the Riemann problem for drift-flux model with modified Chaplygin two-phase flows

In this paper, we concern about the Riemann problem for compressible no-slip drift-flux model which represents a system of quasi-linear partial differential equations derived by averaging the mass and momentum conservation laws with modified Chaplygin two-phase flows. We obtain the exact solution of Riemann problem by elaborately analyzing characteristic fields and discuss the elementary waves namely, shock wave, rarefaction wave and contact discontinuity wave. By employing the equality of pressure and velocity across the middle characteristic field, two nonlinear algebraic equations with two unknowns as gas density ahead and behind the middle wave are formed. The Newton–Raphson method of two variables is applied to find the unknowns with a series of initial data from the literature. Finally, the exact solution for the physical quantities such as gas density, liquid density, velocity, and pressure are illustrated graphically.     

Simple stationary waves in an inclined magnetic field

One component of the solution to the problem of flow around a corner within the scope of magnetohydrodynamics, with the interception or stationary reflection of magnetohydrodynamic shock waves, and also steady-state problems comprising an ionizing shock wave, is the steady-state solution of the equations of magnetohydrodynamics, independent of length but depending on a combination of space variables, for example, on the angle. The flows described by these solutions are called stationary simple waves; they were considered for the first time in [1], where the behavior of the flow was investigated in stationary rotary simple waves, in which no change of density occurs. For a magnetic wave, of parallel velocity, the first integrals were found and the solution was reduced to a quadrature. The investigations and the applications of the solutions obtained for a qualitative construction of the problems of streamline flow were continued in [2–8]. In particular, problems were solved concerning flow around thin bodies of a conducting ideal gas. The general solution of the problem of streamline flow or the intersection of shock waves was not found because stationary simple waves with the magnetic field not parallel to the flow velocity were not investigated. The necessity for the calculation of such a flow may arise during the interpretation of the experimental results [9] in relation to the flow of an ionized gas. In the present paper, we consider stationary simple waves with the magnetic field not parallel to the flow velocity. A system of three nonlinear differential equations, describing fast and slow simple waves, is investigated qualitatively. On the basis of the pattern constructed of the behavior of the integral curves, the change of density, magnetic field, and velocity are found and a classification of the waves is undertaken, according to the nature of the change in their physical quantities. The relation between waves with outgoing and incoming characteristics is explained. A qualitative difference is discovered for the flow investigated from the flow in a magnetic field parallel to the flow velocity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 130–138, September–October, 1976.The author thanks A. A. Barmin and A. G. Kulikovskii for constant interest in the work and for valuable advice.     

Gas content factor in the interaction of shock waves of different intensity in a two-phase gas-liquid medium

The transition from regular to Mach interaction is investigated in connection with the interaction of two plane weak or moderate shock waves of different intensity in a two-phase gas-liquid medium over the entire range of gas contents. A nonmonotonic dependence of the transition limit and the flow parameters on the gas content is detected. The investigation extends the results of [1] corresponding to the reflection of a shock wave from a wall. At intermediate gas contents in the case of opposing shock waves, analogous to the normal reflection of a shock wave from a solid wall, the results are in agreement with [2]. In the case of weak shock waves non-linear asymptotic expansions [3] are employed. In the extreme cases of single-phase media the results coincide with the findings of [3, 4]. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 172–174, November–December, 1986.     

Some exact expressions for the temporal evolution of long Rossby waves on a beta-plane

Some exact expressions are derived to describe the temporal evolution of forced Rossby waves in a two-dimensional beta-plane configuration where the background flow has constant zonal-mean velocity. The meridional length scale of the problem is assumed to be small relative to the zonal length scale and so the long-wave limit of zero aspect ratio is taken. In the case where the background flow velocity is zero, an exact solution is obtained in terms of generalized hypergeometric functions. A late-time asymptotic approximation is obtained and it shows that the solution oscillates with time and its amplitude goes to zero in the limit of infinite time. In the case of a non-zero background flow velocity, the solution is evaluated using two different procedures which give two equivalent expressions in terms of different generalized hypergeometric functions. The late-time asymptotic behaviour is investigated and it is found that the solution approaches a steady state in the limit of infinite time.We also derive a solution in the form of an asymptotic series expansion for the more general situation where a Rossby wave packet is generated by a zonally-localized boundary condition comprising a continuous spectrum of wavenumbers or Fourier modes. The exact solutions found here can be used as leading-order solutions in weakly-nonlinear analyses and other studies involving more realistic configurations for time-dependent Rossby waves or wave packets.     

The flow field near the triple point in steady shock reflection

This paper investigates the flow field near three intersecting shock waves appearing in steady Mach reflection. Results of numerical computations reveal a “von Neumann Paradox”—like feature for weak shock waves, in which the flow field between the reflected and the Mach stem is smooth with no distinct slip flow region and changes rather smoothly. An analytical solution of the Navier–Stokes equations constructed using a polar–coordinate system gives a flow field with the same properties as the numerical simulation.