The study of radially symmetric motion is important for the theory of explosion waves. We construct rigorously self-similar entropy solutions to Riemann initial-boundary value problems for the radially symmetric relativistic Euler equations. We use the assumption of self-similarity to reduce the relativistic Euler equations to a system of nonlinear ordinary differential equations, from which we obtain detailed structures of solutions besides their existence. For the ultra-relativistic Euler equations, we also obtain the uniqueness of the self-similar entropy solution to the Riemann initial-boundary value problems.

The equation of state (EOS) plays a crucial role in hyperbolic conservation laws for the compressible fluid. Whereas, the solid constitutive model with elastic-plastic phase transition makes the analysis of the solid Riemann problem more difficult. In this paper, one-dimensional elastic-perfectly plastic solid Riemann problem is investigated and its exact Riemann solver is proposed. Different from previous works treating the elastic and plastic phases integrally, we resolve the elastic wave and plastic wave separately to understand the complicate nonlinear waves within the solid and then assemble them together to construct the exact Riemann solver for the elastic-perfectly plastic solid. After that, the exact solid Riemann solver is associated with the fluid Riemann solver to decouple the fluid-solid multi-material interaction. Numerical tests, including gas-solid, water-solid high-speed impact problems are simulated by utilizing the modified ghost fluid method (MGFM).

In this paper, we extend using the Runge-Kutta discontinuous Galerkin method together with the front tracking method to simulate the compressible two-medium flow on unstructured meshes. A Riemann problem is constructed in the normal direction in the material interfacial region, with the goal of obtaining a compact, robust and efficient procedure to track the explicit sharp interface precisely. Extensive numerical tests including the gas-gas and gas-liquid flows are provided to show the proposed methodologies possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow and the interfacial vicinities of the two-medium flow in many occasions.

We develop a numerical method to simulate a two-phase compressible flow with sharp phase interface on Eulerian grids. The scheme makes use of a levelset to depict the phase interface numerically. The overall scheme is basically a finite volume scheme. By approximately solving a two-phase Riemann problem on the phase interface, the normal phase interface velocity and the pressure are obtained, which is used to update the phase interface and calculate the numerical flux between the flows of two different phases. We adopt an aggregation algorithm to build cell patches around the phase interface to remove the numerical instability due to the breakdown of the CFL constraint by the cell fragments given by the phase interface depicted using the levelset function. The proposed scheme can handle problems with tangential sliping on the phase interface, topological change of the phase interface and extreme contrast in material parameters in a natural way. Though the perfect conservation of the mass, momentum and energy in global is not achieved, it can be quantitatively identified in what extent the global conservation is spoiled. Some numerical examples are presented to validate the numerical method developed.

A front tracking method combined with the real ghost fluid method (RGFM) is proposed for simulations of fluid interfaces in two-dimensional compressible flows. In this paper the Riemann problem is constructed along the normal direction of interface and the corresponding Riemann solutions are used to track fluid interfaces. The interface boundary conditions are defined by the RGFM, and the fluid interfaces are explicitly tracked by several connected marker points. The Riemann solutions are also used directly to update the flow states on both sides of the interface in the RGFM. In order to validate the accuracy and capacity of the new method, extensive numerical tests including the bubble advection, the Sod tube, the shock-bubble interaction, the Richtmyer-Meshkov instability and the gas-water interface, are simulated by using the Euler equations. The computational results are also compared with earlier computational studies and it shows good agreements including the compressible gas-water system with large density differences.

We construct an approximate Riemann solver for the isentropic Baer−Nunziato two-phaseflow model, that is able to cope with arbitrarily small values of the statistical phasefractions. The solver relies on a relaxation approximation of the model for which theRiemann problem is exactly solved for subsonic relative speeds. In an original manner, theRiemann solutions to the linearly degenerate relaxation system are allowed to dissipatethe total energy in the vanishing phase regimes, thereby enforcing the robustness andstability of the method in the limits of small phase fractions. The scheme is proved tosatisfy a discrete entropy inequality and to preserve positive values of the statisticalfractions and densities. The numerical simulations show a much higher precision and a morereduced computational cost (for comparable accuracy) than standard numerical schemes usedin the nuclear industry. Finally, two test-cases assess the good behavior of the schemewhen approximating vanishing phase solutions.

We formulate a one-dimensional time-dependent non-linear mathematical model for some types of physiological fluid flow in collapsible tubes with discontinuous material properties. The resulting 6 × 6 hyperbolic system is analysed and the associated Riemann problem is solved exactly. Although the solution algorithm deals with idealised cases, it is nonetheless uniquely well-suited for assessing the performance of numerical methods intended for simulating more general situations. Moreover, our model may be a useful starting point for numerical calculations of realistic flows involving rapid and discontinuous material property variations. One important example in mind is the simulation of blood flow in medium-to-large veins in humans. Finally, we also discuss some peculiarities of the model regarding the loss of strict hyperbolicity and uniqueness. In particular we show an example in which the solution of the Riemann problem is non unique.

Existence and admissibility of δ-shock solutions is discussed for the non-convex strictly hyperbolic system of equations

The system is fully nonlinear, i.e. it is nonlinear with respect to both unknowns, and it does not admit the classical Lax-admissible solution for certain Riemann problems. By introducing complex-valued corrections in the framework of the weak asymptotic method, we show that a compressive δ-shock solution resolves such Riemann problems. By letting the approximation parameter tend to zero, the corrections become real valued, and the solutions can be seen to fit into the framework of weak singular solutions defined by Danilov and Shelkovich. Indeed, in this context, we can show that every 2 × 2 system of conservation laws admits δ-shock solutions.

We study the theoretical and numericalcoupling of two hyperbolic systems of conservation laws at a fixed interface. As already proven in the scalar case, the couplingpreserves in a weak sense the continuity of the solution at the interfacewithout imposing the overall conservativity of the coupled model. We develop a detailed analysis of the coupling inthe linear case. In the nonlinear case, we either use a linearized approach or a coupling method based on the solution of a Riemann problem. We discuss both approaches in the case of the coupling of two fluid models at a material contact discontinuity, the models being the usual gas dynamics equations with different equations ofstate. We also study the coupling of two-temperature plasma fluid models and illustrate the approach by numericalsimulations.

A definition of a shock layer of thickness $\delta$ is proposed when a parabolic perturbation is applied to a scalar conservation law. The asymptotic equality $\delta\asymp\sqrt{\varepsilon}$ is established, where $\varepsilon$ denotes viscosity. This equality is proved to be optimal. Nevertheless, the equality $\delta\asymp\varepsilon$ is also proved to be valid for a class of shocks in accordance with the Mises conjecture.

AMS 2000 Mathematics subject classification: Primary 35L67. Secondary 35K65