Source author record

Ravi Samtaney

Ravi Samtaney appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2021arXiv

Linear Stability of an Impulsively Accelerated Density Interface in an Ideal Two-Fluid Plasma

We investigate the linear evolution of Richtmyer-Meshkov (RM) instability in the framework of an ideal two-fluid plasma model. The two-fluid plasma equations of motion are separated into a base state and a set of linearized equations governing the evolution of the perturbations. Different coupling regimes between the charged species are distinguished based on a non-dimensional Debye length parameter $d_{D,0}$. When $d_{D,0}$ is large, the coupling between ions and electrons is sufficiently small that the induced Lorentz force is very weak and the two species evolve as two separate fluids. When $d_{D,0}$ is small, the coupling is strong and the induced Lorentz force is strong enough that the difference between state of ions and electrons is rapidly decreased by the force. As a consequence, the ions and electrons are tightly coupled and evolve like one fluid. The temporal dynamics is divided into two phases: an early phase wherein electron precursor waves are prevalent, and a post ion shock-interface interaction phase during which the RM instability manifests itself. We also examine the effect of an initially applied magnetic field in the streamwise direction characterized by the non-dimensional parameter $β_0$. For a short duration after the ion shock-interface interaction, the growth rate is similar for different initial magnetic field strengths. As time progresses the suppression of the instability due to the magnetic field is observed. The growth rate shows oscillations with a frequency that is related to the ion or electron cyclotron frequency. The instability is suppressed due to the vorticity being transported away from the interface.

preprint2021arXiv

Magnetohydrodynamic shock refraction at an inclined density interface

Shock wave refraction at a sharp density interface is a classical problem in hydrodynamics. Presently, we investigate the strongly planar refraction of a magnetohydrodynamic (MHD) shock wave at an inclined density interface. A magnetic field is applied that is initially oriented either perpendicular and parallel to the motion of incident shock. We explore flow structure by varying the magnitude of the magnetic field governed by the non-dimensional parameter $β\in (0.5, 10^6)$ and the inclination angle of density interface $α\in (0.30, 1.52)$. The regular MHD shock refraction process results in a pair of outer fast shocks (reflected and transmitted) and a set of inner nonlinear magneto-sonic waves. By varying magnetic field (strength and direction) and inclination interface angle, the latter waves can be slow shocks, slow expansion fans, intermediate shocks or slow-mode compound waves. For a chosen incident shock strength, and density ratio, the MHD shock refraction transitions from regular (all nonlinear waves meeting at a single point) into irregular when the inclined density interface angle is less than a critical value. Since the MHD shock refraction is self-similar, we further explore by converting the initial value problem (IVP) into a boundary value problem (BVP) by a self-similar coordinate transformation. The self-similar solution to the BVP is numerically solved using an iterative method and implemented using the p4est adaptive mesh framework. The simulation shows that a Mach stem occurs in irregular MHD shock refraction, the flow structure can be an MHD equivalent to a single Mach reflexion irregular refraction $MRR$ and convex-forwards irregular refraction $CFR$ that occur in the hydrodynamic case.

preprint2020arXiv

A Primitive Variable Discrete Exterior Calculus Discretization of Incompressible Navier-Stokes Equations over Surface Simplicial Meshes

A conservative primitive variable discrete exterior calculus (DEC) discretization of the Navier-Stokes equations is performed. An existing DEC method (Mohamed, M. S., Hirani, A. N., Samtaney, R. (2016). Discrete exterior calculus discretization of incompressible Navier-Stokes equations over surface simplicial meshes. Journal of Computational Physics, 312, 175-191) is modified to this end, and is extended to include the energy-preserving time integration and the Coriolis force to enhance its applicability to investigate the late time behavior of flows on rotating surfaces, i.e., that of the planetary flows. The simulation experiments show second order accuracy of the scheme for the structured-triangular meshes, and first order accuracy for the otherwise unstructured meshes. The method exhibits second order kinetic energy relative error convergence rate with mesh size for inviscid flows. The test case of flow on a rotating sphere demonstrates that the method preserves the stationary state, and conserves the inviscid invariants over an extended period of time.

preprint2020arXiv

Enstrophy transfers in helical turbulence

In this paper we study the enstrophy transers in helical turbulence using direct numerical simulation. We observe that the helicity injection does not have significant effects on the inertial-range energy and helicity spectra ($\sim k^{-5/3}$) and fluxes (constants). We also calculate the separate contributions to enstrophy transfers via velocity to vorticity and vorticity to vorticity channels. There are four different enstrophy fluxes associated with the former channel or vorticity stretching, and one flux associated with the latter channel or vorticity advection. In the inertial range, the fluxes due to vorticity stretching are larger than that due to advection. These transfers too are insensitive to helicity injection.

preprint2020arXiv

Revisiting Reynolds and Nusselt numbers in turbulent thermal convection

In this paper, we extend Grossmann and Lohse's (GL) model [Phys. Rev. Lett. {\bf 86}, 3316 (2001)] for the predictions of Reynolds number (Re) and Nusselt number (Nu) in turbulent Rayleigh-Bénard convection (RBC). Towards this objective, we use functional forms for the prefactors of the dissipation rates in the bulk and the boundary layers. The functional forms arise due to inhibition of nonlinear interactions in the presence of walls and buoyancy compared to free turbulence, along with a deviation of viscous boundary layer profile from Prandtl-Blasius theory. We perform 60 numerical runs on a three-dimensional unit box for a range of Rayleigh numbers (Ra) and Prandtl numbers (Pr) and determine the aforementioned functional forms using machine learning. The revised predictions are in better agreement with the past numerical and experimental results than those of the GL model, especially for extreme Prandtl numbers.

preprint2020arXiv

The generation and sustenance of electric fields in sandstorms

Sandstorms are frequently accompanied by the generation of intense electric fields and lightning. In a very narrow region close to the ground level, sand particles undergo a charge exchange mechanism whereby larger (resp. smaller) sized sand grains become positively (resp. negatively) charged are then entrained by the turbulent fluid motion. Our central hypothesis is that differently sized sand particles get differentially transported by the turbulent flow resulting in a large-scale charge separation, and hence a large-scale electric field. We utilize our simulation framework, comprising of large-eddy simulation of the turbulent atmospheric boundary layer along with sand particle transport and an electrostatic Poisson solver, to investigate the physics of electric fields in sandstorms and thus, to confirm our hypothesis. We utilize the simulation framework to investigate electric fields in weak to strong sandstorms that are characterized by the number density of the sand particles. Our simulations reproduce observational measurements of both mean and RMS fluctuation values of the electric field. We propose a scaling law in which the electric field scales as the two-thirds power of the number density that holds for weak-to-medium sandstorms.

preprint2016arXiv

Discrete exterior calculus discretization of incompressible Navier-Stokes equations over surface simplicial meshes

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.

preprint2015arXiv

Energy transfers in dynamos with small magnetic Prandtl numbers

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.

preprint2015arXiv

Segmental Refinement: A Multigrid Technique for Data Locality

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.

preprint2014arXiv

Energy transfers and magnetic energy growth in small-scale dynamo

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.

preprint2014arXiv

Non-modal stability analysis and transient growth in a magnetized Vlasov plasma

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.

preprint2012arXiv

Benchmarking and scaling studies of a pseudospectral code Tarang for turbulence simulations

Tarang is a general-purpose pseudospectral parallel code for simulating flows involving fluids, magnetohydrodynamics, and Rayleigh-Bé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.