Source author record

Sanjiva K. Lele

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

ResearcherUnclaimed source record

physics.flu-dyn Distributed, Parallel, and Cluster Computing physics.comp-ph

Catalog footprint

What is connected

7works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Aeroacoustics of Twin Rectangular Jets Including Screech: Large-eddy Simulations with Experimental Validation

High-fidelity large-eddy simulations (LES) are performed to investigate aeroacoustic characteristics of jets issuing from twin rectangular nozzles with an aspect ratio of 2:1 at two over-expanded conditions and the design condition. For all three jet conditions simulated, LES predicts qualitatively similar near-field flow statistics to those measured at the University of Cincinnati. Using the Ffowcs Williams-Hawkings method, LES captures the fundamental screech tone and its harmonics fairly well at multiple observer locations in the far-field. Intense jet flapping motions in the near-field along the minor axis, which are influenced by jet-to-jet interactions, are found to correspond to those frequencies. Moreover, the predicted overall sound pressure levels are within 1-2 dB of the experimental measurements. However, the screech tones appear to be intermittent, as the twin-jet interaction pattern varies irregularly. To extract dominant flow structures at the screech frequencies and identify the twin-jet coupling modes, spectral proper orthogonal decomposition (SPOD) analysis is used. SPOD analysis recovers energetic peaks at the screech frequencies, and the corresponding leading modes indicate strong upstream radiation originating from the fifth/sixth shock-cells. For the two over-expanded conditions, the leading modes show anti-symmetric coupling in the minor axis at the fundamental screech frequencies. In contrast, the two jets behave symmetrically with respect to each other in the major axis, in line with the absence of jet flapping in this direction. Furthermore, the leading SPOD eigenvalues turn out to be, at least, two orders of magnitude larger than higher-order eigenvalues, suggesting potential of reduced-order models for the twin-jet screech.

preprint2022arXiv

Analysis of second-moments and their budgets for Richtmyer-Meshkov instability and variable-density turbulence induced by re-shock

Nonlinear Richtmyer--Meshkov instability and mixing transition induced by a Mach 1.45 shock and subsequent re-shock at an interface between two ideal gases (sulfur hexafluoride and air) with high Atwood number are studied with second-moment analysis using data from high-resolution compressible Navier--Stokes simulations. The analysis first addresses the importance of two second-order moments: turbulent mass flux and density-specific-volume covariance, together with their transport equations. These quantities play an essential role in the development of Favre-averaged Reynolds stress and turbulent kinetic energy in this variable-density flow. Then, grid sensitivities and the time evolution of the turbulent quantities which include the second-moments are investigated, followed by a detailed study of the transport equations for the second-moments including the Reynolds stress and the turbulent kinetic energy with well-resolved data before re-shock. After re-shock, budgets of the same but large-scale turbulent quantities are studied with the effects of the subfilter-scale stress taken into account. The budgets of these large-scale quantities are shown to have an insignificant influence from the numerical regularization. Finally, the effects of the subfilter-scale stress on the budgets of the large-scale turbulent quantities with different degrees of filtering are also examined.

preprint2021arXiv

Scalable Parallel Linear Solver for Compact Banded Systems on Heterogeneous Architectures

A scalable algorithm for solving compact banded linear systems on distributed memory architectures is presented. The proposed method factorizes the original system into two levels of memory hierarchies, and solves it using parallel cyclic reduction on both distributed and shared memory. This method has a lower communication footprint across distributed memory partitions compared to conventional algorithms involving data transpose or re-partitioning. The algorithm developed in this work is generalized to cyclic compact banded systems with flexible data decompositions. For cyclic compact banded systems, the method is a direct solver with a deterministic operation and communication counts depending on the matrix size, its bandwidth, and the partition strategy. The implementation and runtime configuration details are discussed for performance optimization. Scalability is demonstrated on the linear solver as well as on a representative fluid mechanics application problem, in which the dominant computational cost is solving the cyclic tridiagonal linear systems of compact numerical schemes on a 3D periodic domain. The algorithm is particularly useful for solving the linear systems arising from the application of compact finite difference operators to a wide range of partial differential equation problems, such as but not limited to the numerical simulations of compressible turbulent flows, aeroacoustics, elastic-plastic wave propagation, and electromagnetics. It alleviates obstacles to their use on modern high performance computing hardware, where memory and computational power are distributed across nodes with multi-threaded processing units.

preprint2020arXiv

Drag of a Heated Sphere at Low Reynolds Numbers in the Presence of Buoyancy

Fully resolved simulations are used to quantify the effects of heat transfer in the presence of buoyancy on the drag of a spatially fixed heated spherical particle at low Reynolds numbers ($Re$) in the range $10^{-3} \le Re \le 10$ in a variable property fluid. The amount of heat addition from the sphere encompasses both, the heating regime where the Boussinesq approximation holds and the regime where it breaks down. The particle is assumed to have a low Biot number which means that the particle is uniformly at the same temperature and has no internal temperature gradients. Scaling buoyancy with inertial and viscous forces yields two related non-dimensional quantities, called Buoyancy Induced Viscous Reynolds Number ($Re_{BV}$) and Buoyancy Induced Inertial Reynolds Number ($Re_{BI}$). For ideal gases, $Re_{BV}$ is analogous to the Grashof number ($Gr$). No assumptions are made on the magnitude of $Re_{BI}$ (or equivalently $Re_{BV}$). The effects of the orientation of gravity relative to the free-stream velocity are examined. Large deviations in the value of the drag coefficient are observed when the Froude number ($Fr$) decreases and/or the temperature of the sphere increases. Under appropriate constraints on $Re_{BI}$ and $Re$, the total drag on a heated sphere in a low $Re$ flow in the presence of buoyancy (mixed convection) is shown to be, within 10% error, the linear superposition of the drag computed in two canonical setups: one being the drag on a steadily moving heated sphere in the absence of buoyancy (forced convection) and the other being natural convection. However, the effect of temperature variation on the drag of a sphere in both, forced and natural convection, is significant.

preprint2014arXiv

Linear and Nonlinear Modeling of a Traveling-Wave Thermoacoustic Heat Engine

We have carried out three-dimensional Navier-Stokes simulations, from quiescent conditions to the limit cycle, of a traveling-wave thermoacoustic heat engine (TAE) composed of a long variable-area resonator shrouding a smaller annular tube, which encloses the hot (HHX) and ambient (AHX) heat-exchangers, and the regenerator (REG). Simulations are wall-resolved, with no-slip and adiabatic conditions enforced at all boundaries, while the heat transfer and drag due to the REG and HXs are modeled. HHX temperatures have been investigated in the range 440K - 500K with AHX temperature fixed at 300K. The initial exponential growth of acoustic energy is due to a network of traveling waves amplified by looping around the REG/HX unit in the direction of the imposed temperature gradient. A simple analytical model demonstrates that such thermoacoustic instability is a Lagrangian thermodynamic process resembling a Stirling cycle. A system-wide linear stability model based on Rott's theory is able to accurately predict the frequency and the growth rate, exhibiting properties consistent with a supercritical Hopf bifurcation. The limit cycle is governed by acoustic streaming -- a rectified flow resulting from high-amplitude acoustics. Its key features are explained with an axially symmetric incompressible model driven by the wave-induced stresses extracted from the compressible calculations. These features include Gedeon streaming and mean recirculations due to flow separation. The first drives the mean advection of hot fluid from the HHX to a secondary heat exchanger (AHX2), in the thermal buffer tube (TBT), necessary to achieve saturation of the acoustic energy growth. The direct evaluation of the energy fluxes reveals that the efficiency of the device deteriorates with the drive ratio and that the acoustic power in the TBT is balanced primarily by the mean advection and thermoacoustic heat transport.

preprint2013arXiv

The Shock and the Turbulence: the story of an interaction

Video highlighting multiple features of the interaction between isotropic turbulence and a nominally planar shock wave, normal to the mean flow direction. The visualizations showcased in this video were obtained from a Direct Numerical Simulation.

preprint2011arXiv

Large-Eddy Simulation of an over-expanded planar nozzle

This fluid dynamics video shows visualizations of a Large-Eddy Simulation (LES) of an over-expanded planar nozzle. This configuration represents the experimental setup of Papamoschou et. al. which found the position of the internal shock to be unstable. Our LES calculations confirm this instability and offer a vibrant and dynamic view of the underlying flow physics. The interaction between shock and turbulent boundary layer is shown as is the subsequent separation region downstream. Numerical Schlieren provide a glimpse of the low frequency shock motion and suggest potential mechanisms for the instability. Key features include the asymmetry of the shock structure (with large and small lambda shocks), compression and expansion waves downstream of the shock and large scale flow reversal. Full details of the experiment and the calculation can be found in the references.