Stability of transonic shock fronts in two-dimensional Euler systems

We study the stability of stationary transonic shock fronts under two-dimensional perturbation in gas dynamics. The motion of the gas is described by the full Euler system. The system is hyperbolic ahead of the shock front, and is a hyperbolic-elliptic composed system behind the shock front. The stability of the shock front and the downstream flow under two-dimensional perturbation of the upstream flow can be reduced to a free boundary value problem of the hyperbolic-elliptic composed system. We develop a method to deal with boundary value problems for such systems. The crucial point is to decompose the system to a canonical form, in which the hyperbolic part and the elliptic part are only weakly coupled in their coefficients. By several sophisticated iterative processes we establish the existence and uniqueness of the solution to the described free boundary value problem. Our result indicates the stability of the transonic shock front and the flow field behind the shock.

     

三维弹性固体中冲击波传输方程的Lagrange描述

在Lagrange坐标中导出了三维非线性弹性固体中冲击波幅度在任意传播方向上的传输方程.导出的方程说明,冲击波的幅度在任意传播方向上随时间的变化率依赖于(i)冲击波阵面紧后方介质运动的一个未知量;(ii)冲击波阵面的两个主曲率;(iii)冲击波法向波速在阵面内的表面梯度;(iv)和冲击波前方介质运动有关的非齐次项,当前方介质处于均匀运动状态时此项为零.文中指出了适当选择传播矢量以简化传输方程的几种方法.我们还得到了一组与介质本构方程无关的、联系冲击波各跳跃量变化率的普适关系.     

Stability of transonic shocks in supersonic flow past a wedge

In this paper we study the stability of transonic shocks in steady supersonic flow past a wedge. We take the potential flow equation as the mathematical model to describe the compressible flow. It is known that in generic case such a problem admits two possible location of shock, connecting the flow ahead it and behind it. They can be distinguished as supersonic-supersonic shock and supersonic-subsonic shock (or transonic shock). Both these possible shocks satisfy the Rankine-Hugoniot conditions and entropy condition. In this paper we prove that the transonic shock is also stable under perturbation of the coming flow provided the pressure at infinity is well controlled.     

Stability of transonic shocks for the full Euler system in supersonic flow past a wedge

We study the stability of transonic shocks in steady supersonic flow past a wedge. It is known that in generic case such a problem admits two possible locations of the shock front, connecting the flow ahead of it and behind it. They can be distinguished as supersonic–supersonic shock and supersonic–subsonic shock (or transonic shock). Both these possible shocks satisfy the Rankine–Hugoniot conditions and the entropy condition. We prove that the transonic shock is conditionally stable under perturbation of the upstream flow or perturbation of wedge boundary. Copyright © 2005 John Wiley & Sons, Ltd.     

Ray expansions of solutions around fronts as a shock-fitting tool for shock loading simulation

In shock loading computations based on an implicit finite-difference scheme, the surfaces of velocity discontinuity and the discontinuity sizes are determined by computing an asymptotic (ray) expansion of the solution behind the front surfaces at every time step. The method for constructing ray expansions is based on a recurrence formulation of the geometric and kinematic consistency conditions for discontinuities of the derivatives of functions that are discontinuous on a moving surface. The algorithm is illustrated by computing a simple example of the out-of-plane motion of an incompressible elastic medium.     

Formulation of the problem of sonic boom by a maneuvering aerofoil as a one-parameter family of Cauchy problems

For the structure of a sonic boom produced by a simple aerofoil at a large distance from its source we take a physical model which consists of a leading shock (LS), a trailing shock (TS) and a one-parameter family of nonlinear wavefronts in between the two shocks. Then we develop a mathematical model and show that according to this model the LS is governed by a hyperbolic system of equations in conservation form and the system of equations governing the TS has a pair of complex eigenvalues. Similarly, we show that a nonlinear wavefront originating from a point on the front part of the aerofoil is governed by a hyperbolic system of conservation laws and that originating from a point on the rear part is governed by a system of conservation laws, which is elliptic. Consequently, we expect the geometry of the TS to be kink-free and topologically different from the geometry of the LS. In the last section we point out an evidence of kinks on the LS and kink-free TS from the numerical solution of the Euler’s equations by Inoue, Sakai and Nishida [5].     

On transonic shocks in a conic divergent nozzle with axi-symmetric exit pressures

In this paper, we establish the existence and stability of a 3-D transonic shock solution to the full steady compressible Euler system in a class of de Laval nozzles with a conic divergent part when a given variable axi-symmetric exit pressure lies in a suitable scope. Thus, for this class of nozzles, we have solved such a transonic shock problem in the axi-symmetric case described by Courant and Friedrichs (1948) in Section 147 of [8]: Given the appropriately large exit pressure pe(x), if the upstream flow is still supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle a shock front intervenes and the gas is compressed and slowed down to subsonic speed so that the position and the strength of the shock front are automatically adjusted such that the end pressure at the exit becomes pe(x).     

Propagation of an electromagnetic shock discontinuity in a non-linear isotropic material

The propagation of an electromagnetic shock discontinuity in a non-linear isotropic dielectric is studied. The velocity of propagation and the electromagnetic vectors behind the shock surface are determined in terms of the unit normal to the surface, the electric and magnetic induction fields ahead of the shock and their magnitudes behind the shock.

---

Zusammenfassung Es wird die Ausbreitung einer elektromagnetischen Stoßwelle in einem nichtlinear isotropen Dielektrikum untersucht. Die Fortpflanzungsgeschwindigkeit und die elektromagnetischen Vektoren hinter der Stoßfläche werden als Funktionen der Einheitsnormalen der Fläche, des elektrischen und magnetischen Induktionsfeldes vor dem Stoß und seiner Beträge hinter dem Stoß bestimmt.

     

Uniform asymptotic analysis for waves in an incompressible elastic rod II. Disturbances superimposed on a pre-stressed state

This paper studies the propagation of disturbances superimposedon a pre-stressed incompressible hyperelastic thin rod. Startingfrom the incremental equations given by Haughton and Ogden (1979,J. Mech. Phys. Solids 27, 179-212 and 489-512), we derive athree parameter-dependent one-dimensional rod equation as thegoverning equation. In particular, it is found that one parameterplays a crucial role. Depending on whether it is larger or smallerthan or equal to a critical value, the shear-wave velocity islarger or smaller than or equal to the bar-wave velocity. Inthe case that these two velocities are equal, there exist travelling-wavesolutions of arbitrary form. This implies that for this particularcase the initial disturbance would propagate along the rod withoutdistortion. To see the influence of the pre-stress in detail,we further consider an initial-value problem with an initialsingularity in the shear strain. The solutions are expressedin terms of integrals through the method of Fourier transform.We then conduct an asymptotic analysis for the solutions. Fora material point in a neighbourhood behind the shear-wave front,the phase function of these integrals has a stationary pointat infinity. Here, we use a technique of uniform asymtotic expansionto handle this case. An asymptotic expansion, correct up toorder O(t^-1), for the shear strain, which is uniformly validin a neighbourhood behind the shear-wave front, is derived.For material points in other spatial domains, the method ofstationary point is applicable, and asymptotic expansions (correctup to order O(t^-1)) are obtained. A novelty is that we are ableto deduce precise qualitative information about the waves inthe far field from our analytic results. Wave profiles for twoconcrete examples are also provided.     

Estimating the height of a tsunami wave propagating over a parabolic bottom in the ray approximation

In this paper, the kinematics of tsunami wave rays and wavefronts propagating over an uneven bottom is considered. Formulas to determine the wave height along a ray tube are obtained. An exact analytical solution for the trajectory of a wave ray over a parabolic bottom is derived. In the wave-ray approximation, this solution makes it possible to analytically determine the heights of tsunami waves over an area with a sloping bottom. The distribution of wave-height maxima over an area with a parabolic bottom is compared with that obtained by numerical computation with a shallow-water model.     

Bistable wavefronts in a diffusive and competitive Lotka-Volterra type system with nonlocal delays

A diffusive Lotka-Volterra type model with nonlocal delays for two competitive species is considered. The existence of a traveling wavefront analogous to a bistable wavefront for a single species is proved by transforming the system with nonlocal delays to a four-dimensional system without delay. Furthermore, in order to prove the asymptotic stability (up to translation) of bistable wavefronts of the system, the existence, regularity and comparison theorem of solutions of the corresponding Cauchy problem are first established for the systems on R by appealing to the theory of abstract functional differential equations. The asymptotic stability (up to translation) of bistable wavefronts are then proved by spectral methods. In particular, we also prove that the spreading speed is unique by upper and lower solutions technique. From the point of view of ecology, our results indicate that the nonlocal delays appeared in the interaction terms are not sensitive to the invasion of species of spatial isolation.     

Geometry and asymptotics of wavefronts

Methods of convex analysis and differential geometry are applied to the study of properties of nonconvex sets in the plane. Constructions of the theory of α-sets are used as a tool for investigation of problems of the control theory and the theory of differential games. The notions of the bisector and of a pseudovertex of a set introduced in the paper, which allow ones to study the geometry of sets and compute their measure of nonconvexity, are of independent interest. These notions are also useful in studies of evolution of sets of attainability of controllable systems and in constructing of wavefronts. In this paper, we develop a numerically-analytical approach to finding pseudovertices of a curve, computation of the measure of nonconvexity of a plane set, and constructing front sets on the basis these data.In the paper, we give the results of numerical constructing of bisectors and wavefronts for plane sets. We demonstrate the relation between nonsmoothness of wavefronts and singularity of the geometry of the initial set. We also single out a class of sets whose bisectors have an asymptote.     

积分中值ξ=ξ(x)的渐近性

对一类不满足g(n)≠0的函数g讨论了第一积分中值定理中ξ=ξ(x)在x→+∞时的渐近性质,并对第二积分中值定理的中值ξ=ξ(x)的渐近性进行了探讨,给出一些相关的结果.     

Study of the stability in the problem on flowing around a wedge. The case of strong wave

We consider the flow of an inviscid nonheatconducting gas in the thermodynamical equilibrium state around a plane infinite wedge and study the stationary solution to this problem, the so-called strong shock wave; the flow behind the shock front is subsonic.We find a solution to a mixed problem for a linear analog of the initial problem, prove that the solution trace on the shock wave is the superposition of direct and reflected waves, and, the main point, justify the Lyapunov asymptotical stability of the strong shock wave provided that the angle at the wedge vertex is small, the uniform Lopatinsky condition is fulfilled, the initial data have a compact support, and the solvability conditions take place if needed (their number depends on the class in which the generalized solution is found).     

A central limit theorem for random sums of random variables

Anscombe (1952) (also see Chung (1974)) has developed a central limit theoremof random sums of independent and identically distributed random variables. Applicability of this theorem in practice, however, is limited since the normalization requires random factors. In this paper we establish sufficient conditions under which the central limit theorem holds when such random factors are replaced by the underlying asymptotic mean and standard ddeviation. An application of this result in the context of shock models is also given.     

The flow of a supersonic ideal gas with “weak” and “strong” shocks over a wedge

The problem of the flow of a uniform supersonic ideal (inviscid and non-heat-conducting) gas over a wedge is considered. If the turning angle of the flow, which is equal to the angle of inclination of the generatrix of the wedge, is less than the maximum value, the problem has two solutions. In the solution with an oblique low-intensity (“weak”) shock, the uniform flow between the shock and the wedge is almost always supersonic. One exception is a small vicinity of the maximum turning angle. For an ideal gas this vicinity does not exceed a fraction of a degree at all Mach numbers. Behind a high-intensity (“strong”) shock, the flow of an ideal gas is always subsonic. “Weak” shocks are observed in all experiments with finite wedges. Some researchers attribute this preference to the “downstream” boundary conditions (“on the right at infinity” for a flow incident on the wedge from the left), and others attribute it to the instability (“Lyapunov” instability) of a flow with a strong shock when it flows over the wedge and to the stability of flow with a weak shock. The results presented below from calculations of the flows that occur for finite wedges within the two-dimensional unsteady Euler equations, when the parameters behind the strong shock are specified on the right-hand boundary, i.e., on the arc of a circle between the wedge and the shock, demonstrate the correctness of the conclusion of the first group of researchers and the incorrectness of the conclusion of the other group. In these calculations, after both small and fairly large perturbations, the flows investigated (which are, in fact, Lyapunov unstable!) return to the solution with a strong shock. In addition, the problem of steady flow over a wedge was regarded as the limit of the two-dimensional non-steady problems at infinite time. Simplification of one of them leads to the problem of the submerged over-expanded supersonic steady outflow. In the ideal gas model this problem is equivalent to flow over a wedge with both weak and strong shocks. All the solutions considered are stable.     

Travelling wavefronts in delayed lattice dynamical systems with global interaction

This paper is concerned with the travelling wavefronts of delayed lattice dynamical systems with global interaction. We establish the existence of the travelling wavefronts by upper–lower solutions technique and Schauder's fixed point theorem when the system satisfies the quasimonotone condition. The nonexistence of the travelling wavefronts of the system is considered by the comparison principle and the corresponding results of the scalar equation. Finally, we apply our main results to the Logistic model and Belousov–Zhabotinskii system on lattice. Our main finding here is that the global interaction can increase the minimal wave speed while the delay can decrease it.     

Exact solution of planar and nonplanar weak shock wave problem in gasdynamics

In the present paper, an analytical approach is used to determine a new exact solution of the problem of one dimensional unsteady adiabatic flow of planer and non-planer weak shock waves in an inviscid ideal fluid. Here it is assumed that the density ahead of the shock front varies according to the power law of the distance from the source of disturbance. The solution of the problem is presented in the form of a power in the distance and the time.     

Wavefronts of a Nonlocal State-dependent Time Delay Model

This paper is concerned with the travelling wavefronts of a nonlocal dispersal cooperation model with harvesting and state-dependent delay,which is assumed to be an increasing function of the population density with lower and upper bound.Especially,state-dependent delay is introduced into a nonlocal reaction-diffusion model.The conditions of Schauder's fixed point theorem are proved by constructing a reasonable set of functionsΓ(see Section 2)and a pair of upper-lower solutions,so the existence of traveling wavefronts is established.The present study is continuation of a previous work that highlights the Laplacian diffusion.     

On the genus of a maximal curve

The upper limit and the first gap in the spectrum of genera of -maximal curves are known, see [34], [16], [35]. In this paper we determine the second gap. Both the first and second gaps are approximately constant times , but this does not hold true for the third gap which is just 1 for while (at most) constant times q for This suggests that the problem of determining the third gap which is the object of current work on -maximal curves could be intricate. Here, we investigate a relevant related problem namely that of characterising those -maximal curves whose genus is equal to the third (or possible the forth) largest value in the spectrum. Our results also provide some new evidence on -maximal curves in connection with Castelnuovo's genus bound, Halphen's theorem, and extremal curves. Received: 1 January 2001 / Revised version: 30 July 2001 / Published online: 23 May 2002