
Two wellknown 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 flowsthey are {\em local} in the turbulent
flows, and {\em nonlocal} in laminar flows.

Discrete exterior calculus (DEC) is a structurepreserving 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 nonDelaunay 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 NavierStokes 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 structuredtriangular 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
shelltoshell energy transfers. These computations indicate that the magnetic
energy growth takes place mainly due to the energy transfers from largescale
velocity field to largescale magnetic field and that the magnetic energy flux
is forward. The steadystate 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 nonnormal leading to
the algebraic growth of the perturbations using nonmodal stability theory. The
typical time scales of the obtained instabilities are of the order of several
plasma periods. The firstorder 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
smallscale dynamo by studying energy transfers, mainly energy fluxes and
shelltoshell 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 largescale or forcingscale velocity field to smallscale
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 generalpurpose pseudospectral parallel code for simulating flows
involving fluids, magnetohydrodynamics, and RayleighB\'{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.