TY - JOUR
T1 - Spontaneous singularity formation in converging cylindrical shock waves
AU - Mostert, W.
AU - Pullin, D. I.
AU - Samtaney, R.
AU - Wheatley, V.
N1 - Publisher Copyright:
© 2018 American Physical Society.
PY - 2018/7
Y1 - 2018/7
N2 - We develop a nonlinear, Fourier-based analysis of the evolution of a perturbed, converging cylindrical strong shock using the approximate method of geometrical shock dynamics (GSD). This predicts that a singularity in the shock-shape geometry, corresponding to a change in Fourier-coefficient decay from exponential to algebraic, is guaranteed to form prior to the time of shock impact at the origin, for arbitrarily small, finite initial perturbation amplitude. Specifically for an azimuthally periodic Mach-number perturbation on an initially circular shock with integer mode number q and amplitude proportional to ϵ1, a singularity in the shock geometry forms at a mean shock radius Ru,c∼(q2ϵ)-1/b1, where b1(γ)<0 is a derived constant and γ the ratio of specific heats. This requires q2ϵ1, q≫1. The constant of proportionality is obtained as a function of γ and is independent of the initial shock Mach number M0. Singularity formation corresponds to the transition from a smooth perturbation to a faceted polygonal form. Results are qualitatively verified by a numerical GSD comparison.
AB - We develop a nonlinear, Fourier-based analysis of the evolution of a perturbed, converging cylindrical strong shock using the approximate method of geometrical shock dynamics (GSD). This predicts that a singularity in the shock-shape geometry, corresponding to a change in Fourier-coefficient decay from exponential to algebraic, is guaranteed to form prior to the time of shock impact at the origin, for arbitrarily small, finite initial perturbation amplitude. Specifically for an azimuthally periodic Mach-number perturbation on an initially circular shock with integer mode number q and amplitude proportional to ϵ1, a singularity in the shock geometry forms at a mean shock radius Ru,c∼(q2ϵ)-1/b1, where b1(γ)<0 is a derived constant and γ the ratio of specific heats. This requires q2ϵ1, q≫1. The constant of proportionality is obtained as a function of γ and is independent of the initial shock Mach number M0. Singularity formation corresponds to the transition from a smooth perturbation to a faceted polygonal form. Results are qualitatively verified by a numerical GSD comparison.
UR - http://www.scopus.com/inward/record.url?scp=85051138830&partnerID=8YFLogxK
U2 - 10.1103/PhysRevFluids.3.071401
DO - 10.1103/PhysRevFluids.3.071401
M3 - Article
AN - SCOPUS:85051138830
SN - 2469-990X
VL - 7
JO - Physical Review Fluids
JF - Physical Review Fluids
IS - 3
M1 - 071401
ER -