Numerical solutions of transonic two-dimensional flows at a ninety-degree wedge

We present numerical results on self-similar two-dimensional Riemann problems governed by the compressible Euler system and the nonlinear wave system, which give rise to a transonic shock. We consider a configuration for a vertical incident shock moving to the right above a rectangular object. The incident shock then interacts with a sonic circle soon after it moves beyond the object, and creates a transonic region. We implement Lax–Liu positive schemes and Strang splitting, and obtain linear correlations of the incident shock strength and the shock strength at the vertical wall. We further implement Roe average methods and finite volume methods on quadrilateral grids to capture a contact discontinuity of the Euler system near the corner of the object. The contact discontinuity creates a new supersonic state and a transonic shock inside the transonic region.     

Transonic shock and rarefaction wave interactions of two-dimensional Riemann problems for the self-similar nonlinear wave system

This paper addresses the self-similar transonic irrotational flow in gas dynamics in two space dimensions.We consider a configuration that the incident shock becomes a transonic shock as it enters the sonic circle, interacts with the rarefaction wave downstream, and then becomes sonic. The rarefaction wave further downstream becomes sonic (degenerate) creating an unknown boundary for the governing system. We present the Riemann data for this configuration, provide the characteristic decomposition of the system, and formulate the boundary value problem for this configuration. The numerical results are presented, and a method to establish the existence result is briefly discussed.     

High‐resolution relaxation scheme for the two‐dimensional Riemann problems in gas dynamics

We present a new relaxation method for the numerical approximation of the two‐dimensional Riemann problems in gas dynamics. The novel feature of the technique proposed here is that it does not require either a Riemann solver or a characteristics decomposition. The high resolution of the method is achieved by using a third‐order reconstruction for the space discretization and a third‐order TVD Runge‐Kutta scheme for the time integration. Numerical experiments, using several configurations of Riemann problems in gas dynamics, are included to confirm the high resolution of the new relaxation scheme. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Transonic shock reflection problems for the self-similar two-dimensional nonlinear wave system

We discuss the regular transonic shock reflections for a model problem of multidimensional conservation laws, the nonlinear wave system. We consider two incident shocks that create reflected shocks, and the state behind the reflected shocks becomes subsonic. We present the existence of the global solution to this configuration, and provide an analysis to handle a degeneracy occurred in the problem and Lipschitz estimates near the sonic boundary. We further implement Lax–Liu positive schemes and Strang splitting, and obtain linear correlations of the incident shock strength and the reflected shock strength. The result obtained in this paper develops a mathematical theory of transonic shock reflection problems. Furthermore the numerical result provides better understanding of the solution structure. This paper provides an application of an important physical problem.     

A global subsonic solution to an interacting transonic shock for the self-similar nonlinear wave equation

We present the existence of the subsonic solution to a two-dimensional Riemann problem governed by a self-similar nonlinear wave equation where the boundary of the subsonic region consists of a transonic shock and the sonic circle. Thus the governing equation becomes a free boundary problem on the transonic shock and degenerates on the sonic circle. By utilizing the barrier methods and iterative methods, we show the well-posedness of the transonic shock in the entire subsonic region and thus establish the global solution. This result does not rely on any smallness of Riemann data.     

The ultra‐relativistic Euler equations

We study the ultra‐relativistic Euler equations for an ideal gas, which is a system of nonlinear hyperbolic conservation laws. We first analyze the single shocks and rarefaction waves and solve the Riemann problem in a constructive way. Especially, we develop an own parametrization for single shocks, which will be used to derive a new explicit shock interaction formula. This shock interaction formula plays an important role in the study of the ultra‐relativistic Euler equations. One application will be presented in this paper, namely, the construction of explicit solutions including shock fronts, which gives an interesting example for the non‐backward uniqueness of our hyperbolic system. Copyright © 2014 John Wiley & Sons, Ltd.     

An Interaction of a Rarefaction Wave and a Transonic Shock for the Self-Similar Two-Dimensional Nonlinear Wave System

We study a two dimensional Riemann problem for the self-similar nonlinear wave system which gives rise to an interaction of a transonic shock and a rarefaction wave. The interesting feature of this problem is that the governing equation changes its type from supersonic in the far field to subsonic near the origin. The subsonic region is then bounded above by the sonic line (degenerate) and below by the transonic shock (free boundary). Furthermore due to the rarefaction wave in the downstream, which interacts with the transonic shock, the problem becomes inhomogeneous and degenerate. We establish the existence result of the global solution to this configuration, and present analysis to understand the solution structure of this problem.     

A MUSCL smoothed particles hydrodynamics for compressible multi‐material flows

By incorporating the Monotone Upwind Scheme of Conservation Law (MUSCL) scheme into the smoothed particles hydrodynamics (SPH) method and making use of an interparticle contact algorithm, we present a MUSCL–SPH scheme of second order for multifluid computations, which extends the Riemann‐solved‐based SPH method. The numerical tests demonstrate high accuracy and resolution of the scheme for both shocks, contact discontinuities, and rarefaction waves in the one‐dimensional shock tube problem. For the two‐dimensional cylindrical Noh and shock‐bubble interaction problems, the MUSCL–SPH scheme can resolve shocks well. Copyright © 2015 John Wiley & Sons, Ltd.     

Bilinear Form and N‐Shock‐Wave Solutions for a (2+1)‐Dimensional Breaking Soliton Equation in Certain Fluids with the Bell Polynomials and Auxiliary Function

In this paper, we will investigate a (2+1)‐dimensional breaking soliton (BS) equation for the (2+1)‐dimensional collision of a Riemann wave with a long wave in certain fluids. Using the Bell polynomials and an auxiliary function, we derive a new bilinear form for the (2+1)‐dimensional BS equation, which is different from those in the previous literatures. One‐, two‐ and N‐shock‐wave solutions are obtained with the Hirota method and symbolic computation. One shock wave is found to be able to stably propagate. Two shock waves are observed to have the parallel collision, oblique collision, and stable propagation of the V‐type structure. In addition, we present the collision between one shock wave and V‐type structure, and the collision between two V‐type structures.     

Undercompressive shock waves in Hall‐magnetohydrodynamics

We study the propagation of nonlinear waves in a Hall‐magnetohydrodynamic model. An asymptotic method is used to derive the Gardner‐Burgers equation for fast magnetosonic waves; here, the flux function is nonconvex with both quadratic and cubic nonlinearities, and the evolution equation involves both second‐ and third‐order derivatives representing diffusion and dispersion terms, respectively. Effects of Hall parameter are discussed on the evolution of waves and their interaction by solving a pair of Riemann problems both analytically and numerically. It is shown that the Hall parameter is responsible for shock splitting—a phenomenon that is completely absent in ideal magnetohydrodynamic; indeed, the Hall parameter plays a significant role in deciding about the structure of the solution that involves undercompressive shocks and their interaction with refracted waves and the Lax shocks. It is found that increasing Hall parameter means increasing dispersion that triggers the physical mechanism causing speed and strength of an undercompressive shock to increase and the wave‐fan width to decrease; numerical solutions substantiate these features predicted by the analytical solution.     

Zu einigen mathematischen problemen in der theorie der transonischen strömung

The work deals with a definition of a weak solution of steady plane transonic flows past a thin profile, with the properties of the solution across a shock wave, and with a derivation of a conservative difference scheme suitable for numerical solution of the above mentioned problem by a finite difference method. The work presents several examples of numerical solution of transonic flows past a profile, through a plane cascade and some three-dimensional results. The numerical results presented are compared with experimental results or with numerical results by other authors.     

Nonlinear Steepest Descent and Numerical Solution of Riemann‐Hilbert Problems

The effective and efficient numerical solution of Riemann‐Hilbert problems has been demonstrated in recent work. With the aid of ideas from the method of nonlinear steepest descent for Riemann‐Hilbert problems, the resulting numerical methods have been shown, in practice, to retain accuracy as values of certain parameters become arbitrarily large. Remarkably, this numerical approach does not require knowledge of local parametrices; rather, the deformed contour is scaled near stationary points at a specific rate. The primary aim of this paper is to prove that this observed asymptotic accuracy is indeed achieved. To do so, we first construct a general theoretical framework for the numerical solution of Riemann‐Hilbert problems. Second, we demonstrate the precise link between nonlinear steepest descent and the success of numerics in asymptotic regimes. In particular, we prove sufficient conditions for numerical methods to retain accuracy. Finally, we compute solutions to the homogeneous Painlevé II equation and the modified Korteweg–de Vries equation to explicitly demonstrate the practical validity of the theory. © 2014 Wiley Periodicals, Inc.     

An Efficient Method for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations

We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography,and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results in a well-balanced scheme that can capture the steady state solution with a remarkable accuracy. This method approximates the source terms via the numerical fluxes produced by an (approximate) Riemann solver for the homogeneous hyperbolic systems with slight additional computation complexity using Newton‘s iterations and numerical integrations. This method solves well the subor super-critical flows, and with a transonic fix, also handles well the transonic flows over the concentration. Numerical examples provide strong evidence on the effectiveness of this new method for both unsteady and steady state calculations.     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

Shatalov‐Sternin's construction of complex WKB solutions and the choice of integration paths

We re‐examine Shatalov‐Sternin's proof of existence of resurgent solutions of a linear ODE. In particular, we take a closer look at the “Riemann surface” (actually, a two‐dimensional complex manifold) whose existence, endless continuability and other properties are claimed by those authors. We present a detailed argument for a part of the “Riemann surface” relevant for the exact WKB method.     

Riemann problem and elementary wave interactions in isentropic magnetogasdynamics

In this paper, we consider the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady simple wave flow of an isentropic, inviscid and perfectly conducting compressible fluid, subjected to a transverse magnetic field. This class of equations includes, as a special case, the equations of isentropic gasdynamics. We study the shock and rarefaction waves and their properties, and discuss the geometry of shock curves using the Riemann invariant coordinates. Under certain conditions, we show the existence and uniqueness of the solution to the Riemann problem for arbitrary initial data, and then discuss the vacuum state in isentropic magnetogasdynamics. Finally, we discuss numerical results for different initial data, and discuss all possible interactions of elementary waves. It is noticed that although the magnetogasdynamic system is more complex than the corresponding gasdynamic system, all the parallel results remain identical. However, unlike the ordinary gasdynamic case, the solution inside rarefaction waves in magnetogasdynamics cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. It is also observed that the presence of a magnetic field makes both the shock and rarefaction stronger compared to what they would have been in the absence of a magnetic field.     

A note on Riemann problems for single‐periodic polyanalytic functions

In this article, a new canonical function has been established to deal with Riemann boundary‐value problem of periodic analytic functions discussed in 16 . In comparison with the corresponding result in 16 , the expression of solution obtained here is much simpler. Then, we demonstrate the equivalence of solutions for the homogeneous Riemann problem. What's more, we obtain the precise rank of matrix of coefficients for the system of linear algebraic equations (4.35) in 16 . Those results can simplify the discussion of Riemann problem of single‐periodic polyanalytic functions in 16 .     

Delta-shocks and vacuums in zero-pressure gas dynamics by the flux approximation

In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.     

Two‐grid methods for time‐harmonic Maxwell equations

In this paper, we develop several two‐grid methods for the Nédélec edge finite element approximation of the time‐harmonic Maxwell equations. We first present a two‐grid method that uses a coarse space to solve the original problem and then use a fine space to solve a corresponding symmetric positive definite problem. Then, we present two types of iterative two‐grid methods, one is to add the kernel of the curl ‐operator in the fine space to a coarse mesh space to solve the original problem and the other is to use an inner iterative method for dealing with the kernel of the curl ‐operator in the fine space and the coarse space, separately. We provide the error estimates for the first two methods and present numerical experiments to show the efficiency of our methods.Copyright © 2012 John Wiley & Sons, Ltd.     

Compact finite difference scheme for the solution of a time fractional partial integro‐differential equation with a weakly singular kernel

In this paper, we develop a high‐order finite difference scheme for the solution of a time fractional partial integro‐differential equation with a weakly singular kernel. The fractional derivative is used in the Riemann‐Liouville sense. We prove the unconditional stability and convergence of scheme using energy method and show that the convergence order is . We provide some numerical experiments to confirm the efficiency of suggested scheme. The results of numerical experiments are compared with analytical solutions to show the efficiency of proposed scheme. It is illustrated that the numerical results are in good agreement with theoretical ones.