High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability: Comparison to experimental data and to amplitude growth model predictions

The reshocked single-mode Richtmyer-Meshkov instability is simulated in two spatial dimensions using the fifth- and ninth-order weighted essentially nonoscillatory shock-capturing method with uniform spatial resolution of 256 points per initial perturbation wavelength. The initial conditions and computational domain are modeled after the single-mode, Mach 1.21 air(acetone)/SF6 shock tube experiment of Collins and Jacobs [J. Fluid Mech. 464, 113 (2002)]. The simulation densities are shown to be in very good agreement with the corrected experimental planar laser-induced fluorescence images at selected times before reshock of the evolving interface. Analytical, semianalytical, and phenomenological linear and nonlinear, impulsive, perturbation, and potential flow models for single-mode Richtmyer-Meshkov unstable perturbation growth are summarized. The simulation amplitudes are shown to be in very good agreement with the experimental data and with the predictions of linear amplitude growth models for small times, and with those of nonlinear amplitude growth models at later times up to the time at which the driver-based expansion in the experiment (but not present in the simulations or models) expands the layer before reshock. The qualitative and quantitative differences between the fifth- and ninth-order simulation results are discussed. Using a local and global quantitative metric, the prediction of the Zhang and Sohn [Phys. Fluids 9, 1106 (1997)] nonlinear Padé model is shown to be in best overall agreement with the simulation amplitudes before reshock. The sensitivity of the amplitude growth model predictions to the initial growth rate from linear instability theory, the post-shock Atwood number and amplitude, and the velocity jump due to the passage of the shock through the interface is also investigated numerically

Stokes' first problem for some non-Newtonian fluids: Results and mistakes

The well-known problem of unidirectional plane flow of a fluid in a half-space due to the impulsive motion of the plate it rests upon is discussed in the context of the second-grade and the Oldroyd-B non-Newtonian fluids. The governing equations are derived from the conservation laws of mass and momentum and three correct known representations of their exact solutions given. Common mistakes made in the literature are identified. Simple numerical schemes that corroborate the analytical solutions are constructed.Comment: 10 pages, 2 figures; accepted for publication in Mechanics Research Communications; v2 corrects a few typo

A survey on stationary problems, Green's functions and spectrum of Sturm–Liouville problem with nonlocal boundary conditions

In this paper, we present a survey of recent results on the Green's functions and on spectrum for stationary problems with nonlocal boundary conditions. Results of Lithuanian mathematicians in the field of differential and numerical problems with nonlocal boundary conditions are described. *The research was partially supported by the Research Council of Lithuania (grant No. MIP-047/2014)

Fractional Calculus - Theory and Applications

In recent years, fractional calculus has led to tremendous progress in various areas of science and mathematics. New definitions of fractional derivatives and integrals have been uncovered, extending their classical definitions in various ways. Moreover, rigorous analysis of the functional properties of these new definitions has been an active area of research in mathematical analysis. Systems considering differential equations with fractional-order operators have been investigated thoroughly from analytical and numerical points of view, and potential applications have been proposed for use in sciences and in technology. The purpose of this Special Issue is to serve as a specialized forum for the dissemination of recent progress in the theory of fractional calculus and its potential applications

Global three-dimensional flow of a neutron superfluid in a spherical shell in a neutron star

We integrate for the first time the hydrodynamic Hall-Vinen-Bekarevich-Khalatnikov equations of motion of a $^{1}S_{0}$-paired neutron superfluid in a rotating spherical shell, using a pseudospectral collocation algorithm coupled with a time-split fractional scheme. Numerical instabilities are smoothed by spectral filtering. Three numerical experiments are conducted, with the following results. (i) When the inner and outer spheres are put into steady differential rotation, the viscous torque exerted on the spheres oscillates quasiperiodically and persistently (after an initial transient). The fractional oscillation amplitude ($\sim 10^{-2}$) increases with the angular shear and decreases with the gap width. (ii) When the outer sphere is accelerated impulsively after an interval of steady differential rotation, the torque increases suddenly, relaxes exponentially, then oscillates persistently as in (i). The relaxation time-scale is determined principally by the angular velocity jump, whereas the oscillation amplitude is determined principally by the gap width. (iii) When the mutual friction force changes suddenly from Hall-Vinen to Gorter-Mellink form, as happens when a rectilinear array of quantized Feynman-Onsager vortices is destabilized by a counterflow to form a reconnecting vortex tangle, the relaxation time-scale is reduced by a factor of $\sim 3$ compared to (ii), and the system reaches a stationary state where the torque oscillates with fractional amplitude $\sim 10^{-3}$ about a constant mean value. Preliminary scalings are computed for observable quantities like angular velocity and acceleration as functions of Reynolds number, angular shear, and gap width. The results are applied to the timing irregularities (e.g., glitches and timing noise) observed in radio pulsars.Comment: 6 figures, 23 pages. Accepted for publication in Astrophysical Journa