• Two well-known turbulence models that describe the energy spectrum in the inertial and dissipative ranges simultaneously are by Pao~(1965) and Pope~(2000). In this paper, we compute the energy spectrum $E(k)$ and energy flux $\Pi(k)$ using direct numerical simulations on grids up to $4096^3$, and show consistency between the numerical results and the predictions by the aforementioned models. We also construct a model for laminar flows that predicts $E(k)\sim k^{-1} \exp(-k)$ and $\Pi(k)\sim k \exp(-k)$. Our model predictions match with the numerical results. We emphasize differences on the energy transfers in the two flows---they are {\em local} in the turbulent flows, and {\em nonlocal} in laminar flows.
  • Discrete exterior calculus (DEC) is a structure-preserving numerical framework for partial differential equations solution, particularly suitable for simplicial meshes. A longstanding and widespread assumption has been that DEC requires special (Delaunay) triangulations, which complicated the mesh generation process especially on curved surfaces. This paper presents numerical evidences demonstrating that this restriction is unnecessary. Convergence experiments are carried out for various physical problems using both Delaunay and non-Delaunay triangulations. Signed diagonal definition for the key DEC operator (Hodge star) is adopted. The errors converge as expected for all considered meshes and experiments. This relieves the DEC paradigm from unnecessary triangulation limitation.
  • A conservative discretization of incompressible Navier-Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.
  • We investigate a domain decomposed multigrid technique, segmental refinement, for solving general nonlinear elliptic boundary value problems. Brandt and Diskin first proposed this method in 1994; we continue this work by analytically and experimentally investigating its complexity. We confirm that communication of traditional parallel multigrid can be eliminated on fine grids with modest amounts of extra work and storage while maintaining the asymptotic exactness of full multigrid, although we observe a dependence on an additional parameter not considered in the original analysis. We present a communication complexity analysis that quantifies the communication costs ameliorated by segmental refinement and report performance results with up to 64K cores of a Cray XC30.
  • We perform numerical simulation of dynamo with magnetic Prandtl number $\mathrm{Pm} =0.2$ on $1024^3$ grid, and compute the energy fluxes and the shell-to-shell energy transfers. These computations indicate that the magnetic energy growth takes place mainly due to the energy transfers from large-scale velocity field to large-scale magnetic field and that the magnetic energy flux is forward. The steady-state magnetic energy is much smaller than the kinetic energy, rather than equipartition; this is because the magnetic Reynolds number is near the dynamo transition regime. We also contrast our results with those for dynamo with $\mathrm{Pm} =20$ and decaying dynamo.
  • Collisionless plasmas, such as those encountered in tokamaks, exhibit a rich variety of instabilities. The physical origin, triggering mechanisms and fundamental understanding of many plasma instabilities, however, are still open problems. We investigate the stability properties of a collisionless Vlasov plasma in a stationary homogeneous magnetic field. We narrow the scope of our investigation to the case of Maxwellian plasma. For the first time using a fully kinetic approach we show the emergence of the local instability, a transient growth, followed by classical Landau damping in a stable magnetized plasma. We show that the linearized Vlasov operator is non-normal leading to the algebraic growth of the perturbations using non-modal stability theory. The typical time scales of the obtained instabilities are of the order of several plasma periods. The first-order distribution function and the corresponding electric field are calculated and the dependence on the magnetic field and perturbation parameters is studied. Our results offer a new scenario of the emergence and development of plasma instabilities on the kinetic scale.
  • In this letter we investigate the dynamics of magnetic energy growth in small-scale dynamo by studying energy transfers, mainly energy fluxes and shell-to-shell energy transfers. We perform dynamo simulations for magnetic Prandtl number $\mathrm{Pm}=20$ on $1024^3$ grid using pseudospectral method. We demonstrate that the magnetic energy growth is caused by nonlocal energy transfers from the large-scale or forcing-scale velocity field to small-scale magnetic field. The peak of these energy transfers move towards lower wavenumbers as dynamo evolves, which is the reason why the integral scale of the magnetic field increases with time. The energy transfers $U2U$ (velocity to velocity) and $B2B$ (magnetic to magnetic) are forward and local.
  • Tarang is a general-purpose pseudospectral parallel code for simulating flows involving fluids, magnetohydrodynamics, and Rayleigh-B\'{e}nard convection in turbulence and instability regimes. In this paper we present code validation and benchmarking results of Tarang. We performed our simulations on $1024^3$, $2048^3$, and $4096^3$ grids using the {\em HPC system} of IIT Kanpur and {\em Shaheen} of KAUST. We observe good "weak" and "strong" scaling for Tarang on these systems.