Shock-resolved Navier–Stokes simulation of the Richtmyer–Meshkov instability start-up at a light–heavy interface

The single-mode Richtmyer–Meshkov instability is investigated using a first-order perturbation of the two-dimensional Navier–Stokes equations about a one-dimensional unsteady shock-resolved base flow. A feature-tracking local refinement scheme is used to fully resolve the viscous internal structure of the shock. This method captures perturbations on the shocks and their influence on the interface growth throughout the simulation, to accurately examine the start-up and early linear growth phases of the instability. Results are compared to analytic models of the instability, showing some agreement with predicted asymptotic growth rates towards the inviscid limit, but significant discrepancies are noted in the transient growth phase. Viscous effects are found to be inadequately predicted by existing models

A study of planar Richtmyer-Meshkov instability in fluids with Mie-Grüneisen equations of state

We present a numerical comparison study of planar Richtmyer-Meshkov instability with the intention of exposing the role of the equation of state. Results for Richtmyer-Meshkov instability in fluids with Mie-Grüneisen equations of state derived from a linear shock-particle speed Hugoniot relationship (Jeanloz, J. Geophys. Res. 94, 5873, 1989; McQueen et al., High Velocity Impact Phenomena (1970), pp. 294–417; Menikoff and Plohr, Rev. Mod. Phys. 61(1), 75 1989) are compared to those from perfect gases under nondimensionally matched initial conditions at room temperature and pressure. The study was performed using Caltech’s Adaptive Mesh Refinement, Object-oriented C++ (AMROC) (Deiterding, Adaptive Mesh Refinement: Theory and Applications (2005), Vol. 41, pp. 361–372; Deiterding, “Parallel adaptive simulation of multi-dimensional detonation structures,” Ph.D. thesis (Brandenburgische Technische Universität Cottbus, September 2003)) framework with a low-dissipation, hybrid, center-difference, limiter patch solver (Ward and Pullin, J. Comput. Phys. 229, 2999 (2010)). Results for single and triple mode planar Richtmyer-Meshkov instability when a reflected shock wave occurs are first examined for mid-ocean ridge basalt (MORB) and molybdenum modeled by Mie-Grüneisen equations of state. The single mode case is examined for incident shock Mach numbers of 1.5 and 2.5. The planar triple mode case is studied using a single incident Mach number of 2.5 with initial corrugation wavenumbers related by k_1 = k_2+k_3. Comparison is then drawn to Richtmyer-Meshkov instability in perfect gases with matched nondimensional pressure jump across the incident shock, post-shock Atwood ratio, post-shock amplitude-to-wavelength ratio, and time nondimensionalized by Richtmyer’s linear growth time constant prediction. Differences in start-up time and growth rate oscillations are observed across equations of state. Growth rate oscillation frequency is seen to correlate directly to the oscillation frequency for the transmitted and reflected shocks. For the single mode cases, further comparison is given for vorticity distribution and corrugation centerline shortly after shock interaction. Additionally, we examine single mode Richtmyer-Meshkov instability when a reflected expansion wave is present for incident Mach numbers of 1.5 and 2.5. Comparison to perfect gas solutions in such cases yields a higher degree of similarity in start-up time and growth rate oscillations. The formation of incipient weak waves in the heavy fluid driven by waves emanating from the perturbed transmitted shock is observed when an expansion wave is reflected

Interaction Between Convection and Pulsation

This article reviews our current understanding of modelling convection dynamics in stars. Several semi-analytical time-dependent convection models have been proposed for pulsating one-dimensional stellar structures with different formulations for how the convective turbulent velocity field couples with the global stellar oscillations. In this review we put emphasis on two, widely used, time-dependent convection formulations for estimating pulsation properties in one-dimensional stellar models. Applications to pulsating stars are presented with results for oscillation properties, such as the effects of convection dynamics on the oscillation frequencies, or the stability of pulsation modes, in classical pulsators and in stars supporting solar-type oscillations.Comment: Invited review article for Living Reviews in Solar Physics. 88 pages, 14 figure

Chaotic and pseudochaotic attractors of perturbed fractional oscillator

We consider a nonlinear oscillator with fractional derivative of the order alpha. Perturbed by a periodic force, the system exhibits chaotic motion called fractional chaotic attractor (FCA). The FCA is compared to the ``regular'' chaotic attractor. The properties of the FCA are discussed and the ``pseudochaotic'' case is demonstrated.Comment: 20 pages, 7 figure

Stochastic Perturbations of Periodic Orbits with Sliding

Vector fields that are discontinuous on codimension-one surfaces are known as Filippov systems and can have attracting periodic orbits involving segments that are contained on a discontinuity surface of the vector field. In this paper we consider the addition of small noise to a general Filippov system and study the resulting stochastic dynamics near such a periodic orbit. Since a straight-forward asymptotic expansion in terms of the noise amplitude is not possible due to the presence of discontinuity surfaces, in order to quantitatively determine the basic statistical properties of the dynamics, we treat different parts of the periodic orbit separately. Dynamics distant from discontinuity surfaces is analyzed by the use of a series expansion of the transitional probability density function. Stochastically perturbed sliding motion is analyzed through stochastic averaging methods. The influence of noise on points at which the periodic orbit escapes a discontinuity surface is determined by zooming into the transition point. We combine the results to quantitatively determine the effect of noise on the oscillation time for a three-dimensional canonical model of relay control. For some parameter values of this model, small noise induces a significantly large reduction in the average oscillation time. By interpreting our results geometrically, we are able to identify four features of the relay control system that contribute to this phenomenon.Comment: 44 pages, 9 figures, submitted to: J Nonlin. Sc

Dark solitons in cigar-shaped Bose-Einstein condensates in double-well potentials

We study the statics and dynamics of dark solitons in a cigar-shaped Bose-Einstein condensate confined in a double-well potential. Using a mean-field model with a non-cubic nonlinearity, appropriate to describe the dimensionality crossover regime from one to three dimensional, we obtain branches of solutions in the form of single- and multiple-dark soliton states, and study their bifurcations and stability. It is demonstrated that there exist dark soliton states which do not have a linear counterpart and we highlight the role of anomalous modes in the excitation spectra. Particularly, we show that anomalous mode eigenfrequencies are closely connected to the characteristic soliton frequencies as found from the solitons' equations of motion, and how anomalous modes are related to the emergence of instabilities. We also analyze in detail the role of the height of the barrier in the double well setting, which may lead to instabilities or decouple multiple dark soliton states.Comment: 35 pages, 12 figure