We present nsCouette, a highly scalable software tool to solve the Navier-Stokes equations for incompressible fluid flow between differentially heated and independently rotating, concentric cylinders. It is based on a pseudospectral spatial discretization and dynamic time-stepping. It is implemented in modern Fortran with a hybrid MPI-OpenMP parallelization scheme and thus designed to compute turbulent flows at high Reynolds and Rayleigh numbers. An additional GPU implementation (C-CUDA) for intermediate problem sizes and a basic version for turbulent pipe flow (nsPipe) are also provided.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Magnetorotational instability (MRI) is one of the fundamental processes in astrophysics, driving angular momentum transport and mass accretion in a wide variety of cosmic objects. Despite much theoretical/numerical and experimental efforts over the last decades, its saturation mechanism and amplitude, which sets the angular momentum transport rate, remains not well understood, especially in the limit of high resistivity, or small magnetic Prandtl numbers typical to interiors (dead zones) of protoplanetary disks, liquid cores of planets and liquid metals in laboratory. Using direct numerical simulations, in this paper we investigate the nonlinear development and saturation properties of the helical magnetorotational instability (HMRI) -- a relative of the standard MRI -- in a magnetized Taylor-Couette flow at very low magnetic Prandtl number (correspondingly at low magnetic Reynolds number) relevant to liquid metals. For simplicity, the ratio of azimuthal field to axial field is kept fixed. From the linear theory of HMRI, it is known that the Elsasser number, or interaction parameter determines its growth rate and plays a special role in the dynamics. We show that this parameter is also important in the nonlinear problem. By increasing its value, a sudden transition from weakly nonlinear, where the system is slightly above the linear stability threshold, to strongly nonlinear, or turbulent regime occurs. We calculate the azimuthal and axial energy spectra corresponding to these two regimes and show that they differ qualitatively. Remarkably, the nonlinear state remains in all cases nearly axisymmetric suggesting that HMRI turbulence is quasi two-dimensional in nature. Although the contribution of non-axisymmetric modes increases moderately with the Elsasser number, their total energy remains much smaller than that of the axisymmetric ones.

The magnetorotational instability (MRI) is thought to be a powerful source of turbulence in Keplerian accretion disks. Motivated by recent laboratory experiments, we study the MRI driven by an azimuthal magnetic field in an electrically conducting fluid sheared between two concentric rotating cylinders. By adjusting the rotation rates of the cylinders, we approximate angular velocity profiles $\omega \propto r^{q}$. We perform direct numerical simulations of a steep profile close to the Rayleigh line $q \gtrsim -2 $ and a quasi-Keplerian profile $q \approx -3/2$ and cover wide ranges of Reynolds ($Re\le 4\cdot10^4$) and magnetic Prandtl numbers ($0\le Pm \le 1$). In the quasi-Keplerian case, the onset of instability depends on the magnetic Reynolds number, with $Rm_c \approx 50$, and angular momentum transport scales as $\sqrt{Pm} Re^2$ in the turbulent regime. The ratio of Maxwell to Reynolds stresses is set by $Rm$. At the onset of instability both stresses have similar magnitude, whereas the Reynolds stress vanishes or becomes even negative as $Rm$ increases. For the profile close to the Rayleigh line, the instability shares these properties as long as $Pm\gtrsim0.1$, but exhibits a markedly different character if $Pm\rightarrow 0$, where the onset of instability is governed by the Reynolds number, with $Re_c \approx 1250$, transport is via Reynolds stresses and scales as $Re^2$. At intermediate $Pm=0.01$ we observe a continuous transition from one regime to the other, with a crossover at $Rm=\mathcal{O}(100)$. Our results give a comprehensive picture of angular momentum transport of the MRI with an imposed azimuthal field.

We report a direct-numerical-simulation study of Taylor-Couette flow in the quasi-Keplerian regime at shear Reynolds numbers up to $\mathcal{O}(10^5)$. Quasi-Keplerian rotating flow has been investigated for decades as a simplified model system to study the origin of turbulence in accretion disks that is not fully understood. The flow in this study is axially periodic and thus the experimental end-wall effects on the stability of the flow are avoided. Using optimal linear perturbations as initial conditions, our simulations find no sustained turbulence: the strong initial perturbations distort the velocity profile and trigger turbulence that eventually decays.

The protoMIRAX hard X-ray imaging telescope is a balloon-borne experiment developed as a pathfinder for the MIRAX satellite mission. The experiment consists essentially in a coded-aperture hard X-ray (30-200 keV) imager with a square array (13$\times$13) of 2mm-thick planar CZT detectors with a total area of 169 cm$^2$. The total, fully-coded field-of-view is $21^{\circ}\times 21^{\circ}$ and the angular resolution is 1$^{\circ}$43'. In this paper we describe the protoMIRAX instrument and all the subsystems of its balloon gondola, and we show simulated results of the instrument performance. The main objective of protoMIRAX is to carry out imaging spectroscopy of selected bright sources to demonstrate the performance of a prototype of the MIRAX hard X-ray imager. Detailed background and imaging simulations have been performed for protoMIRAX balloon flights. The 3$\sigma$ sensitivity for the 30-200 keV range is ~1.9 $\times$ 10$^{-5}$ photons cm$^{-2}$ s$^{-1}$ for an integration time of 8 hs at an atmospheric depth of 2.7 g cm$^{-2}$ and an average zenith angle of 30$^{\circ}$. We have developed an attitude control system for the balloon gondola and new data handling and ground systems that also include prototypes for the MIRAX satellite. We present the results of Monte Carlo simulations of the camera response at balloon altitudes, showing the expected background level and the detailed sensitivity of protoMIRAX. We also present the results of imaging simulations of the Crab region. The results show that protoMIRAX is capable of making spectral and imaging observations of bright hard X-ray source fields. Furthermore, the balloon observations will carry out very important tests and demonstrations of MIRAX hardware and software in a near space environment.

A hybrid-parallel direct-numerical-simulation method with application to turbulent Taylor-Couette flow is presented. The Navier-Stokes equations are discretized in cylindrical coordinates with the spectral Fourier-Galerkin method in the axial and azimuthal directions, and high-order finite differences in the radial direction. Time is advanced by a second-order, semi-implicit projection scheme, which requires the solution of five Helmholtz/Poisson equations, avoids staggered grids and renders very small slip velocities. Nonlinear terms are computed with the pseudospectral method. The code is parallelized using a hybrid MPI-OpenMP strategy, which is simpler to implement, reduces inter-node communications and is more efficient compared to a flat MPI parallelization. A strong scaling study shows that the hybrid code maintains very good scalability up to more than 20000 processor cores and thus allows to perform simulations at higher resolutions than previously feasible, and opens up the possibility to simulate turbulent Taylor-Couette flows at Reynolds numbers up to $\mathcal{O}(10^5)$. This enables to probe hydrodynamic turbulence in Keplerian flows in experimentally relevant regimes.