The Richtmyer–Meshkov instability in magnetohydrodynamics

In ideal magnetohydrodynamics (MHD), the Richtmyer–Meshkov instability can be suppressed by the presence of a magnetic field. The interface still undergoes some growth, but this is bounded for a finite magnetic field. A model for this flow has been developed by considering the stability of an impulsively accelerated, sinusoidally perturbed density interface in the presence of a magnetic field that is parallel to the acceleration. This was accomplished by analytically solving the linearized initial value problem in the framework of ideal incompressible MHD. To assess the performance of the model, its predictions are compared to results obtained from numerical simulation of impulse driven linearized, shock driven linearized, and nonlinear compressible MHD for a variety of cases. It is shown that the analytical linear model collapses the data from the simulations well. The predicted interface behavior well approximates that seen in compressible linearized simulations when the shock strength, magnetic field strength, and perturbation amplitude are small. For such cases, the agreement with interface behavior that occurs in nonlinear simulations is also reasonable. The effects of increasing shock strength, magnetic field strength, and perturbation amplitude on both the flow and the performance of the model are investigated. This results in a detailed exposition of the features and behavior of the MHD Richtmyer–Meshkov flow. For strong shocks, large initial perturbation amplitudes, and strong magnetic fields, the linear model may give a rough estimate of the interface behavior, but it is not quantitatively accurate. In all cases examined the accuracy of the model is quantified and the flow physics underlying any discrepancies is examine

Power-law versus log-law in wall-bounded turbulence: A large-eddy simulation perspective

The debate whether the mean streamwise velocity in wall-bounded turbulent flows obeys a log-law or a power-law scaling originated over two decades ago, and continues to ferment in recent years. As experiments and direct numerical simulation can not provide sufficient clues, in this study we present an insight into this debate from a large-eddy simulation (LES) viewpoint. The LES organically combines state-of-the-art models (the stretched-vortex model and inflow rescaling method) with a virtual-wall model derived under different scaling law assumptions (the log-law or the power-law by George and Castillo [“Zero-pressure-gradient turbulent boundary layer,” Appl. Mech. Rev.50, 689 (1997)]). Comparison of LES results for Re θ ranging from 10^5 to 10^(11) for zero-pressure-gradient turbulent boundary layer flows are carried out for the mean streamwise velocity, its gradient and its scaled gradient. Our results provide strong evidence that for both sets of modeling assumption (log law or power law), the turbulence gravitates naturally towards the log-law scaling at extremely large Reynolds numbers

Magnetic reconnection and stochastic plasmoid chains in high-Lundquist-number plasmas

A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is \normEeff\sim 0.02 independently of $S$ for $S\gg10^4$. The plasmoid flux ($\Psi$) and half-width ($w_x$) distribution functions scale as $f(\Psi)\sim \Psi^{-2}$ and $f(w_x)\sim w_x^{-2}$. The joint distribution of $\Psi$ and $w_x$ shows that plasmoids populate a triangular region $w_x\gtrsim\Psi/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic "monster" plasmoids with $w_x\sim 10$% of the system size are shown to emerge in just a few Alfv\'en times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.Comment: 5 pages, 6 figures, submitted for publicatio

Large-eddy simulation of separation and reattachment of a flat plate turbulent boundary layer

We present large-eddy simulations (LES) of separation and reattachment of a flat-plate turbulent boundary-layer flow. Instead of resolving the near wall region, we develop a two-dimensional virtual wall model which can calculate the time- and space-dependent skin-friction vector field at the wall, at the resolved scale. By combining the virtual-wall model with the stretched-vortex subgrid-scale (SGS) model, we construct a self-consistent framework for the LES of separating and reattaching turbulent wall-bounded flows at large Reynolds numbers. The present LES methodology is applied to two different experimental flows designed to produce separation/reattachment of a flat-plate turbulent boundary layer at medium Reynolds number Re_θ based on the momentum boundary-layer thickness θ. Comparison with data from the first case at Re_θ=2000 demonstrates the present capability for accurate calculation of the variation, with the streamwise co-ordinate up to separation, of the skin friction coefficient, Re_θ, the boundary-layer shape factor and a non-dimensional pressure-gradient parameter. Additionally the main large-scale features of the separation bubble, including the mean streamwise velocity profiles, show good agreement with experiment. At the larger Re_θ=11000 of the second case, the LES provides good postdiction of the measured skin-friction variation along the whole streamwise extent of the experiment, consisting of a very strong adverse pressure gradient leading to separation within the separation bubble itself, and in the recovering or reattachment region of strongly-favourable pressure gradient. Overall, the present two-dimensional wall model used in LES appears to be capable of capturing the quantitative features of a separation-reattachment turbulent boundary-layer flow at low to moderately large Reynolds numbers

Magnetohydrodynamic implosion symmetry and suppression of Richtmyer-Meshkov instability in an octahedrally symmetric field

We present numerical simulations of ideal magnetohydrodynamics showing suppression of the Richtmyer-Meshkov instability in spherical implosions in the presence of an octahedrally symmetric magnetic field. This field configuration is of interest owing to its high degree of spherical symmetry in comparison with previously considered dihedrally symmetric fields. The simulations indicate that the octahedral field suppresses the instability comparably to the other previously considered candidate fields for light-heavy interface accelerations while retaining a highly symmetric underlying flow even at high field strengths. With this field, there is a reduction in the root-mean-square perturbation amplitude of up to approximately 50% at representative time under the strongest field tested while maintaining a homogeneous suppression pattern compared to the other candidate fields

Effects of seed magnetic fields on magnetohydrodynamic implosion structure and dynamics

The effects of various seed magnetic fields on the dynamics of cylindrical and spherical implosions in ideal magnetohydrodynamics are investigated. Here, we present a fundamental investigation of this problem utilizing cylindrical and spherical Riemann problems under three seed field configurations to initialize the implosions. The resulting flows are simulated numerically, revealing rich flow structures, including multiple families of magnetohydrodynamic shocks and rarefactions that interact non-linearly. We fully characterize these flow structures, examine their axi- and spherisymmetry-breaking behaviour, and provide data on asymmetry evolution for different field strengths and driving pressures for each seed field configuration. We find that out of the configurations investigated, a seed field for which the implosion centre is a saddle point in at least one plane exhibits the least degree of asymmetry during implosion

Converging cylindrical shocks in ideal magnetohydrodynamics

We consider a cylindrically symmetrical shock converging onto an axis within the framework of ideal, compressible-gas non-dissipative magnetohydrodynamics (MHD). In cylindrical polar co-ordinates we restrict attention to either constant axial magnetic field or to the azimuthal but singular magnetic field produced by a line current on the axis. Under the constraint of zero normal magnetic field and zero tangential fluid speed at the shock, a set of restricted shock-jump conditions are obtained as functions of the shock Mach number, defined as the ratio of the local shock speed to the unique magnetohydrodynamic wave speed ahead of the shock, and also of a parameter measuring the local strength of the magnetic field. For the line current case, two approaches are explored and the results compared in detail. The first is geometrical shock-dynamics where the restricted shock-jump conditions are applied directly to the equation on the characteristic entering the shock from behind. This gives an ordinary-differential equation for the shock Mach number as a function of radius which is integrated numerically to provide profiles of the shock implosion. Also, analytic, asymptotic results are obtained for the shock trajectory at small radius. The second approach is direct numerical solution of the radially symmetric MHD equations using a shock-capturing method. For the axial magnetic field case the shock implosion is of the Guderley power-law type with exponent that is not affected by the presence of a finite magnetic field. For the axial current case, however, the presence of a tangential magnetic field ahead of the shock with strength inversely proportional to radius introduces a length scale R = √μ_0/p_0I/(2π) where I is the current, μ_0 is the permeability, and p_0 is the pressure ahead of the shock. For shocks initiated at r ≫ R, shock convergence is first accompanied by shock strengthening as for the strictly gas-dynamic implosion. The diverging magnetic field then slows the shock Mach number growth producing a maximum followed by monotonic reduction towards magnetosonic conditions, even as the shock accelerates toward the axis. A parameter space of initial shock Mach number at a given radius is explored and the implications of the present results for inertial confinement fusion are discussed

Linear simulations of the cylindrical Richtmyer-Meshkov instability in magnetohydrodynamics

Numerical simulations and analysis indicate that the Richtmyer-Meshkov instability (RMI) is suppressed in ideal magnetohydrodynamics (MHD) in Cartesian slab geometry. Motivated by the presence of hydrodynamic instabilities in inertial confinement fusion and suppression by means of a magnetic field, we investigate the RMI via linear MHD simulations in cylindrical geometry. The physical setup is that of a Chisnell-type converging shock interacting with a density interface with either axial or azimuthal (2D) perturbations. The linear stability is examined in the context of an initial value problem (with a time-varying base state) wherein the linearized ideal MHD equations are solved with an upwind numerical method. Linear simulations in the absence of a magnetic field indicate that RMI growth rate during the early time period is similar to that observed in Cartesian geometry. However, this RMI phase is short-lived and followed by a Rayleigh-Taylor instability phase with an accompanied exponential increase in the perturbation amplitude. We examine several strengths of the magnetic field (characterized by β = 2p/B ) and observe a significant suppression of the instability for β ≤ 4. The suppression of the instability is attributed to the transport of vorticity away from the interface by Alfvén fronts

Converging cylindrical magnetohydrodynamic shock collapse onto a power-law-varying line current

We investigate the convergence behaviour of a cylindrical, fast magnetohydrodynamic (MHD) shock wave in a neutrally ionized gas collapsing onto an axial line current that generates a power law in time, azimuthal magnetic field. The analysis is done within the framework of a modified version of ideal MHD for an inviscid, non-dissipative, neutrally ionized compressible gas. The time variation of the magnetic field is tuned such that it approaches zero at the instant that the shock reaches the axis. This configuration is motivated by the desire to produce a finite magnetic field at finite shock radius but a singular gas pressure and temperature at the instant of shock impact. Our main focus is on the variation with shock radius r, as r→0, of the shock Mach number M(r) and pressure behind the shock p(r) as a function of the magnetic field power-law exponent μ ⩾ 0, where μ = 0 gives a constant-in-time line current. The flow problem is first formulated using an extension of geometrical shock dynamics (GSD) into the time domain to take account of the time-varying conditions ahead of the converging shock, coupled with appropriate shock-jump conditions for a fast, symmetric MHD shock. This provides a pair of ordinary differential equations describing both M(r) and the time evolution on the shock, as a function of r, constrained by a collapse condition required to achieve tuned shock convergence. Asymptotic, analytical results for M(r) and p(r) are obtained over a range of μ for general γ, and for both small and large r. In addition, numerical solutions of the GSD equations are performed over a large range of r, for selected parameters using γ=5/3. The accuracy of the GSD model is verified for some cases using direct numerical solution of the full, radially symmetric MHD equations using a shock-capturing method. For the GSD solutions, it is found that the physical character of the shock convergence to the axis is a strong function of μ. For 0 ⩽ μ < 4/13, M and p both approach unity at shock impact r=0 owing to the dominance of the strong magnetic field over the amplifying effects of geometrical convergence. When μ⩾0.816 (for γ=5/3), geometrical convergence is dominant and the shock behaves similarly to a converging gas dynamic shock with singular M(r) and p(r), r→0. For 4/13 < μ ⩽ 0.816 three distinct regions of M(r) variation are identified. For each of these p(r) is singular at the axis