scholarly journals An exact Riemann-solver-based solution for regular shock refraction

2009 ◽  
Vol 627 ◽  
pp. 33-53 ◽  
Author(s):  
P. DELMONT ◽  
R. KEPPENS ◽  
B. VAN DER HOLST

We study the classical problem of planar shock refraction at an oblique density discontinuity, separating two gases at rest. When the shock impinges on the density discontinuity, it refracts, and in the hydrodynamical case three signals arise. Regular refraction means that these signals meet at a single point, called the triple point. After reflection from the top wall, the contact discontinuity becomes unstable due to local Kelvin–Helmholtz instability, causing the contact surface to roll up and develop the Richtmyer–Meshkov instability (RMI). We present an exact Riemann-solver-based solution strategy to describe the initial self-similar refraction phase, by which we can quantify the vorticity deposited on the contact interface. We investigate the effect of a perpendicular magnetic field and quantify how its addition increases the deposition of vorticity on the contact interface slightly under constant Atwood number. We predict wave-pattern transitions, in agreement with experiments, von Neumann shock refraction theory and numerical simulations performed with the grid-adaptive code AMRVAC. These simulations also describe the later phase of the RMI.

2001 ◽  
Vol 449 ◽  
pp. 395-411 ◽  
Author(s):  
LUCIANO REZZOLLA ◽  
OLINDO ZANOTTI

A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general, when solving the Riemann problem numerically, the determination of the specific wave pattern produced is obtained through some initial guess which can be successively discarded or improved. We here discuss a new procedure, suitable for implementation in an exact Riemann solver in one dimension, which removes the initial ambiguity in the wave pattern. In particular we focus our attention on the relativistic velocity jump between the two initial states and use this to determine, through some analytic conditions, the wave pattern produced by the decay of the initial discontinuity. The exact Riemann problem is then solved by means of calculating the root of a nonlinear equation. Interestingly, in the case of two rarefaction waves, this root can even be found analytically. Our procedure is straightforward to implement numerically and improves the efficiency of numerical codes based on exact Riemann solvers.


2015 ◽  
Vol 773 ◽  
pp. 366-394 ◽  
Author(s):  
Xisheng Luo ◽  
Minghu Wang ◽  
Ting Si ◽  
Zhigang Zhai

The interaction of a planar shock wave ($M\approx 1.2$) with an $\text{SF}_{6}$ polygonal inhomogeneity surrounded by air is experimentally investigated. Six polygons including a square, two types of rectangle, two types of triangle, and a diamond are generated by the soap film technique developed in our previous work, in which thin pins are used as angular vertexes to avoid the pressure singularities caused by the surface tension. The evolutions of the shock-accelerated $\text{SF}_{6}$ polygons are captured by a high-speed schlieren system from which wave systems and the interface characteristics can be clearly identified. Both regular and irregular refraction phenomena are observed outside the volume, and more complex wave patterns, including transmitted shock, refracted shock, Mach stem and the interactions between them, are found inside the volume. Two typical irregular refraction phenomena (free precursor refraction, FPR, and free precursor von Neumann refraction, FNR) are observed and analysed, and the transition from FPR to FNR is found, providing the experimental evidence for the transition between different wave patterns numerically found in the literature. Combined with our previous work (Zhai et al., J. Fluid Mech., vol. 757, 2014, pp. 800–816), the reciprocal transitions between FPR and FNR are experimentally confirmed. The velocities and trajectories of the triple points are further measured and it is found that the motions of the triple points are self-similar or pseudo-stationary. Besides the shock dynamics phenomena, the evolutions of these shocked heavy polygonal volumes, which are quite different from the light ones, are captured and found to be closely related to their initial shapes. Specifically, for square and rectangular geometries, the different width–height ratios result in different behaviours of shock–shock interaction inside the volume, and subsequently different features for the outward jet and the interface. Quantitatively, the time-variations of the interface scales, such as the width and the normalized displacements of the edges, are obtained and compared with those from previous work. The comparison illustrates the superiority of the interface formation method and the significant effect of the initial interface shape on the interface features. Furthermore, the characteristics of the vortex core, including the velocity and vortex spacing, are experimentally measured, and the vortex velocity is compared with those from some circulation models to check the validity of the models. The results in the present work enrich understanding of the shock refraction phenomenon and the database of research into Richtmyer–Meshkov instability (RMI).


1994 ◽  
Vol 98 (979) ◽  
pp. 325-339 ◽  
Author(s):  
E. F. Toro ◽  
A. Chakraborty

Abstract An improved version (HLLC) of the Harten, Lax, van Leer Riemann solver (HLL) for the steady supersonic Euler equations is presented. Unlike the HLL, the HLLC version admits the presence of the slip line in the structure of the solution. This leads to enhanced resolution of computed slip lines by Godunov type methods. We assess the HLLC solver in the context of the first order Godunov method and the second order weighted average flux method (WAF). It is shown that the improvement embodied in the HLLC solver over the HLL solver is virtually equivalent to incorporating the exact Riemann solver.


2001 ◽  
Vol 437 ◽  
pp. 255-282 ◽  
Author(s):  
O. IGRA ◽  
X. WU ◽  
J. FALCOVITZ ◽  
T. MEGURO ◽  
K. TAKAYAMA ◽  
...  

The complex flow and wave pattern following an initially planar shock wave transmitted through a double-bend duct is studied experimentally and theoretically/numerically. Several different double-bend duct geometries are investigated in order to assess their effects on the accompanying flow and shock wave attenuation while passing through these ducts. The effect of the duct wall roughness on the shock wave attenuation is also studied. The main flow diagnostic used in the experimental part is either an interferometric study or alternating shadow–schlieren diagnostics. The photos obtained provide a detailed description of the flow evolution inside the ducts investigated. Pressure measurements were also taken in some of the experiments. In the theoretical/numerical part the conservation equations for an inviscid, perfect gas were solved numerically. It is shown that the proposed physical model (Euler equations), which is solved by using the second-order-accurate, high-resolution GRP (generalized Riemann problem) scheme, can simulate such a complex, time-dependent process very accurately. Specifically, all wave patterns are numerically simulated throughout the entire interaction process. Excellent agreement is found between the numerical simulation and the experimental results. The efficiency of a double-bend duct in providing a shock wave attenuation is clearly demonstrated.


2010 ◽  
Vol 6 (S274) ◽  
pp. 441-444
Author(s):  
Rony Keppens ◽  
Peter Delmont ◽  
Zakaria Meliani

AbstractThe interaction of (strong) shock waves with localized density changes is of particular relevance to laboratory as well as astrophysical research. Shock tubes have been intensively studied in the lab for decades and much has been learned about shocks impinging on sudden density contrasts. In astrophysics, modern observations vividly demonstrate how (even relativistic) winds or jets show complex refraction patterns as they encounter denser interstellar material.In this contribution, we highlight recent insights into shock refraction patterns, starting from classical up to relativistic hydro and extended to magnetohydrodynamic scenarios. Combining analytical predictions for shock refraction patterns exploiting Riemann solver methodologies, we confront numerical, analytical and (historic) laboratory insights. Using parallel, grid-adaptive simulations, we demonstrate the fate of Richtmyer-Meshkov instabilities when going from gaseous to magnetized plasma scenarios. The simulations invoke idealized configurations closely resembling lab analogues, while extending them to relativistic flow regimes.


2008 ◽  
Vol 103 (1) ◽  
pp. 77
Author(s):  
Junsheng Fang ◽  
Don Hadwin ◽  
Xiujuan Ma

We compute spectra and Brown measures of some non self-adjoint operators in $(M_2(\mathsf {C}), \frac{1}{2}\mathrm{Tr})*(M_2(\mathsf{C}), \frac{1}{2}\mathrm{Tr})$, the reduced free product von Neumann algebra of $M_2(\mathsf {C})$ with $M_2(\mathsf {C})$. Examples include $AB$ and $A+B$, where $A$ and $B$ are matrices in $(M_2(\mathsf {C}), \frac{1}{2}\mathrm{Tr})*1$ and $1*(M_2(\mathsf {C}), \frac{1}{2}\mathrm{Tr})$, respectively. We prove that $AB$ is an R-diagonal operator (in the sense of Nica and Speicher [12]) if and only if $\mathrm{Tr}(A)=\mathrm{Tr}(B)=0$. We show that if $X=AB$ or $X=A+B$ and $A,B$ are not scalar matrices, then the Brown measure of $X$ is not concentrated on a single point. By a theorem of Haagerup and Schultz [9], we obtain that if $X=AB$ or $X=A+B$ and $X\neq \lambda 1$, then $X$ has a nontrivial hyperinvariant subspace affiliated with $(M_2(\mathsf{C}), \frac{1}{2}\mathrm{Tr})*(M_2(\mathsf{C}), \frac{1}{2}\mathrm{Tr})$.


2009 ◽  
Vol 410-411 ◽  
pp. 373-379 ◽  
Author(s):  
Camille Robert ◽  
Lanouar Ben Ayed ◽  
Arnaud Delamézière ◽  
Phillippe dal Santo ◽  
J.L. Batoz

The numerical simulation of the Single Point Incremental Forming process (SPIF) is time consuming due to the necessity to take into account various non-linearity such as the material behaviour, large strain deformation and the evolution of the tool-flange contact. Classical contact algorithms give good agreement with experimental results, but are time consuming. In this paper, we investigate the development of a procedure to simplify the management of the contact interface between the tool and the sheet. Nodes with imposed displacements are determined by a geometrical approximation of the deformed sheet. In order to have a better approximation of the local stresses in the flange, a pressure is applied on the tool side of the elements in the contact zone. The pressure value is obtained by an analytical model. A classical contact algorithm and the present simplified approach are compared in terms of an incremental forming benchmark. It has been shown that, for the benchmark problem studied here, a CPU time reduction of approximately 65% can be achieved while at the same time have good simulation results.


2000 ◽  
Vol 14 (2) ◽  
pp. 117-131 ◽  
Author(s):  
EMMANUELLE DECLERCQ ◽  
ALAIN FORESTIER ◽  
JEAN-MARC HÉRARD ◽  
XAVIER LOUIS ◽  
GÉRARD POISSANT

Sign in / Sign up

Export Citation Format

Share Document