Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
The von Neumann Triple Point Paradox Richard Sanders 1∗ and Allen M. Tesdall 2† 1 Department of Mathematics, University of Houston, Houston, TX 77204, USA sanders@math.uh.edu 2 Fields Institute, Toronto, ON M5T 3J1, Canada and Department of Mathematics, University of Houston, Houston, TX 77204, USA atesdall@fields.utoronto.ca Summary. We describe the problem of weak shock reflection off a wedge and discuss the triple point paradox that arises. When the shock is sufficiently weak and the wedge is thin, Mach reflection appears to be observed but is impossible according to what von Neumann originally showed in 1943. We summarize some recent numerical results for weak shock reflection problems for the unsteady transonic small disturbance equations, the nonlinear wave system, and the Euler equations. Rather than finding a standard but mathematically inadmissible Mach reflection with a shock triple point, the solutions contain a complex structure: there is a sequence of triple points and supersonic patches in a tiny region behind the leading triple point, with an expansion fan originating at each triple point. The sequence of patches may be infinite, and we refer to this structure as Guderley Mach reflection. The presence of the expansion fans at the triple points resolves the paradox. We describe some recent experimental evidence which is consistent with these numerical findings. Key words: self-similar solutions, two-dimensional Riemann problems, triple point paradox 1 Introduction Consider a planar normal shock in an inviscid compressible and calorically perfect gas which impinges on a fixed wedge with apex half angle θ w ,see Figure 1. Given an upstream state with density ρ = ρ r , velocity u = v =0 and pressure p = p r , one calculates that downstream of a fast (i.e., u + c) shock ∗ Research supported by the National Science Foundation, Grant DMS 03-06307. † Research supported by the National Science Foundation, Grant DMS 03-06307, NSERC grant 312587-05, and the Fields Institute.
114 R. Sanders and A.M. Tesdall I R θ w R I S M (a) (b) Fig. 1. A planar shock moving from left to right impinges on a wedge. After contact, I indicates the incident shock and R indicates the reflected shock. On the right, the dotted line S indicates a slip line and M is the Mach stem. Regular reflection is depicted on the left. Irregular reflection is depicted on the right. U l R I R U l I ? U r U U r S ? M (a) (b) Fig. 2. A blow-up of the incident and reflected shock intersection. Regular reflection is on the left and irregular on the right. The constant states upstream and downstream of the incident shock are denoted by U r and U l . Whether or not constant states indicated by the question marks exist depends on the strength of I. p l = 2γ p r γ +1 M 2 − γ − 1 γ +1 , ρ l ρ r = u l = 2 ( M − 1 ) , c r γ +1 M (1) 2 (γ +1)M 2+(γ − 1) M 2 , where γ denotes the ratio of specific heats, and M > 1 denotes the shock Mach number defined as the Rankine–Hugoniot shock speed divided by the upstream speed of sound c r = √ γp r /ρ r . Following interaction, a number of self-similar (with respect to the wedge apex) reflection patterns are possible, depending on the values of M and θ w . This wedge reflection problem has a rich history, experimentally, analytically, and numerically. Probably the earliest and most significant analytical result was found by von Neumann [Neu43]. In this work were first formulated the equations which describe two and three planar shocks meeting at a point separated by constant states, see Figure 2. The two shock theory leads to what is known as regular reflection. The three shock theory leads to Mach reflection. For supersonic regular reflection, state U immediately behind the reflected shock R is supersonic and becomes subsonic across a sonic line downstream (toward the wedge’s apex). When the incident shock angle is increased
- Page 69 and 70: 58 E.J. Dean and R. Glowinski 7 Num
- Page 71 and 72: 60 E.J. Dean and R. Glowinski Fig.
- Page 73 and 74: 62 E.J. Dean and R. Glowinski Assum
- Page 75 and 76: Higher Order Time Stepping for Seco
- Page 77 and 78: u n+1 h Optimal Higher Order Time D
- Page 79 and 80: Optimal Higher Order Time Discretiz
- Page 81 and 82: Optimal Higher Order Time Discretiz
- Page 83 and 84: Optimal Higher Order Time Discretiz
- Page 85 and 86: Optimal Higher Order Time Discretiz
- Page 87 and 88: Optimal Higher Order Time Discretiz
- Page 89 and 90: Optimal Higher Order Time Discretiz
- Page 91 and 92: Optimal Higher Order Time Discretiz
- Page 93 and 94: Optimal Higher Order Time Discretiz
- Page 95 and 96: Optimal Higher Order Time Discretiz
- Page 97 and 98: Optimal Higher Order Time Discretiz
- Page 99 and 100: Optimal Higher Order Time Discretiz
- Page 101 and 102: Optimal Higher Order Time Discretiz
- Page 103 and 104: 96 I. Sazonov et al. To provide a p
- Page 105 and 106: 98 I. Sazonov et al. In the first s
- Page 107 and 108: 100 I. Sazonov et al. Fig. 1. An ex
- Page 109 and 110: 102 I. Sazonov et al. H z 1 exact F
- Page 111 and 112: 104 I. Sazonov et al. (a) (b) Fig.
- Page 113 and 114: 106 I. Sazonov et al. Scattering Wi
- Page 115 and 116: 108 I. Sazonov et al. 6.4 Scatterin
- Page 117 and 118: 110 I. Sazonov et al. (a) (b) Fig.
- Page 119: 112 I. Sazonov et al. [MHP96] K. Mo
- Page 123 and 124: 116 R. Sanders and A.M. Tesdall imp
- Page 125 and 126: 118 R. Sanders and A.M. Tesdall alo
- Page 127 and 128: 120 R. Sanders and A.M. Tesdall (a)
- Page 129 and 130: 122 R. Sanders and A.M. Tesdall D C
- Page 131 and 132: 124 R. Sanders and A.M. Tesdall 8.6
- Page 133 and 134: 126 R. Sanders and A.M. Tesdall 0.3
- Page 135 and 136: 128 R. Sanders and A.M. Tesdall [TR
- Page 137 and 138: 132 S. Lapin et al. Ω R γ Ω 2 Γ
- Page 139 and 140: 134 S. Lapin et al. ∫ ∂ 2 ∫
- Page 141 and 142: 136 S. Lapin et al. Ω R γ Fig. 3.
- Page 143 and 144: 138 S. Lapin et al. 4 Energy Inequa
- Page 145 and 146: 140 S. Lapin et al. 5 Numerical Exp
- Page 147 and 148: 142 S. Lapin et al. Fig. 6. Contour
- Page 149 and 150: 144 S. Lapin et al. Fig. 9. Obstacl
- Page 151 and 152: Domain Decomposition and Electronic
- Page 153 and 154: Domain Decomposition Approach for C
- Page 155 and 156: Domain Decomposition Approach for C
- Page 157 and 158: Domain Decomposition Approach for C
- Page 159 and 160: Domain Decomposition Approach for C
- Page 161 and 162: Domain Decomposition Approach for C
- Page 163 and 164: Domain Decomposition Approach for C
- Page 165 and 166: Domain Decomposition Approach for C
- Page 167 and 168: Domain Decomposition Approach for C
- Page 169 and 170: Numerical Analysis of a Finite Elem
The von Neumann Triple Point Paradox<br />
Richard S<strong>and</strong>ers 1∗ <strong>and</strong> Allen M. Tesdall 2†<br />
1 Department of Mathematics, University of Houston, Houston, TX 77204, USA<br />
s<strong>and</strong>ers@math.uh.edu<br />
2 Fields Institute, Toronto, ON M5T 3J1, Canada <strong>and</strong> Department of<br />
Mathematics, University of Houston, Houston, TX 77204, USA<br />
atesdall@fields.utoronto.ca<br />
Summary. We describe the problem of weak shock reflection off a wedge <strong>and</strong> discuss<br />
the triple point paradox that arises. When the shock is sufficiently weak <strong>and</strong> the<br />
wedge is thin, Mach reflection appears to be observed but is impossible according to<br />
what von Neumann originally showed in 1943. We summarize some recent numerical<br />
results for weak shock reflection problems for the unsteady transonic small disturbance<br />
equations, the nonlinear wave system, <strong>and</strong> the Euler equations. Rather than<br />
finding a st<strong>and</strong>ard but mathematically inadmissible Mach reflection with a shock<br />
triple point, the solutions contain a complex structure: there is a sequence of triple<br />
points <strong>and</strong> supersonic patches in a tiny region behind the leading triple point, with<br />
an expansion fan originating at each triple point. The sequence of patches may be<br />
infinite, <strong>and</strong> we refer to this structure as Guderley Mach reflection. The presence<br />
of the expansion fans at the triple points resolves the paradox. We describe some<br />
recent experimental evidence which is consistent with these numerical findings.<br />
Key words: self-similar solutions, two-dimensional Riemann problems, triple<br />
point paradox<br />
1 Introduction<br />
Consider a planar normal shock in an inviscid compressible <strong>and</strong> calorically<br />
perfect gas which impinges on a fixed wedge with apex half angle θ w ,see<br />
Figure 1. Given an upstream state with density ρ = ρ r , velocity u = v =0<br />
<strong>and</strong> pressure p = p r , one calculates that downstream of a fast (i.e., u + c)<br />
shock<br />
∗ Research supported by the National Science Foundation, Grant DMS 03-06307.<br />
† Research supported by the National Science Foundation, Grant DMS 03-06307,<br />
NSERC grant 312587-05, <strong>and</strong> the Fields Institute.