We develop a minimal model for \textit{pulsar glitches} by introducing a solid-crust potential in the three-dimensional (3D) Gross-Pitaevskii-Poisson equation (GPPE), which we have used earlier to study gravitationally bound Bose-Einstein Condensates (BECs), i.e., bosonic stars. In the absence of the crust potential, we show that, if we rotate such a bosonic star, it is threaded by vortices. We then show, via extensive direct numerical simulations (DNSs), that the interaction of these vortices with the crust potential yields (a) stick-slip dynamics and (b) dynamical glitches. We demonstrate that, if enough momentum is transferred to the crust from the bosonic star, then the vortices are expelled from the star and the crust's angular momentum $J_c$ exhibits features that can be interpreted naturally as glitches. From the time series of $J_c$, we compute the cumulative probability distribution functions (CPDFs) of event sizes, event durations, and waiting times. We show that these CPDFs have signatures of self-organized criticality (SOC), which have been seen in observations on pulsar glitches.

We investigate the clustering properties of radio sources in the Alternative Data Release 1 of the TIFR GMRT Sky Survey (TGSS), focusing on large angular scales, where previous analyses have detected a large clustering signal. After appropriate data selection, the TGSS sample we use contains ~110,000 sources selected at 150 MHz over ~70% of the sky. The survey footprint is largely superimposed on that of the NRAO VLA Sky Survey (NVSS) with the majority of TGSS sources having a counterpart in the NVSS sample. These characteristics make TGSS suitable for large-scale clustering analyses and facilitate the comparison with the results of previous studies. In this analysis we focus on the angular power spectrum, although the angular correlation function is also computed to quantify the contribution of multiple-component radio sources. We find that on large angular scales, corresponding to multipoles $2 \leq \ell \leq 30$, the amplitude of the TGSS angular power spectrum is significantly larger than that of the NVSS. We do not identify any observational systematic effects that may explain this mismatch. We have produced a number of physically motivated models for the TGSS angular power spectrum and found that all of them fail to match observations, even when taking into account observational and theoretical uncertainties. The same models provide a good fit to the angular spectrum of the NVSS sources. These results confirm the anomalous nature of the TGSS large-scale power, which has no obvious physical origin and seems to indicate that unknown systematic errors are present in the TGSS dataset.

We show that a shell-model version of the three-dimensional Hall-magnetohydrodynamic (3D Hall-MHD) equations provides a natural theoretical model for investigating the multiscaling behaviors of velocity and magnetic structure functions. We carry out extensive numerical studies of this shell model, obtain the scaling exponents for its structure functions, in both the low-$k$ and high-$k$ power-law ranges of 3D Hall-MHD, and find that the extended-self-similarity (ESS) procedure is helpful in extracting the multiscaling nature of structure functions in the high-$k$ regime, which otherwise appears to display simple scaling. Our results shed light on intriguing solar-wind measurements.

We present a natural framework for studying the persistence problem in two-dimensional fluid turbulence by using the Okubo-Weiss parameter $\Lambda$ to distinguish between vortical and extensional regions. We then use a direct numerical simulation (DNS) of the two-dimensional, incompressible Navier--Stokes equation with Ekman friction to study probability distribution functions (PDFs) of the persistence times of vortical and extensional regions by employing both Eulerian and Lagrangian measurements. We find that, in the Eulerian case, the persistence-time PDFs have exponential tails; by contrast, this PDF for Lagrangian particles, in vortical regions, has a power-law tail with an exponent $\theta=2.9\pm0.2$.

We carry out systematic and high-resolution studies of dynamo action in a shell model for magnetohydrodynamic (MHD) turbulence over wide ranges of the magnetic Prandtl number $Pr_{\rm M}$ and the magnetic Reynolds number $Re_{\rm M}$. Our study suggests that it is natural to think of dynamo onset as a nonequilibrium, first-order phase transition between two different turbulent, but statistically steady, states. The ratio of the magnetic and kinetic energies is a convenient order parameter for this transition. By using this order parameter, we obtain the stability diagram (or nonequilibrium phase diagram) for dynamo formation in our MHD shell model in the $(Pr^{-1}_{\rm M}, Re_{\rm M})$ plane. The dynamo boundary, which separates dynamo and no-dynamo regions, appears to have a fractal character. We obtain hysteretic behavior of the order parameter across this boundary and suggestions of nucleation-type phenomena.

We derive relations for the decay of the kinetic and magnetic energies and the growth of the Taylor and integral scales in unforced, incompressible, homogeneous and isotropic three-dimensional magnetohydrodynamic (3DMHD) turbulence with power-law initial energy spectra. We also derive bounds for the decay of the cross- and magnetic helicities. We then present results from systematic numerical studies of such decay both within the context of an MHD shell model and direct numerical simulations (DNS) of 3DMHD. We show explicitly that our results about the power-law decay of the energies hold for times $t<t_*$, where $t_*$ is the time at which the integral scales become comparable to the system size. For $t<t_*$, our numerical results are consistent with those predicted by the principle of `permanence of large eddies'.