BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ashoke De

Ashoke De contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
46works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

46 published item(s)

preprint2023arXiv

Understanding the liquid jet break-up in various regimes at elevated pressure using a compressible VOF-LPT coupled framework

The present work develops a compressible VOF-LPT coupled solver in OpenFOAM and utilizes it to simulate a LJICF numerically. This methodology helps accurately predict a complex primary breakup in the Eulerian framework and the secondary atomization of spherical droplets using a computationally efficient LPT method. The coupled solver with AMR is rigorously validated for a liquid jet in crossflow at varying operating conditions. We have further carried out a thorough investigation to study the effect of momentum flux ratio and weber number on the various flow features and liquid jet break-up phenomenon in a crossflow while identifying the stream-wise location of the liquid jet breakup region. At low momentum flux ratios in the bag breakup regime, the predictions reveal that the liquid jet breakup occurs due to the growth of similar instability as usually observed in the high-speed liquid sheet atomization. The short wavelength assumption of the inviscid dispersion relation resembles the Kelvin-Helmholtz type instability observed in this case, as opposed to Rayleigh-Taylor instability at high momentum flux ratio in the surface breakup regime. It is also proposed that the shear breakup along the transverse edges of the liquid column occurs due to the shear layer instability of the air passing around the liquid column. The simulation wavelength closely matches the Williamson correlation for shear layer instability around cylinders: a shape similar to the cross-section of the bottom of the liquid column. The results show a distinct streamer or bifurcation phenomenon at low momentum flux ratios and moderate weber numbers. Further investigation suggests that the internal liquid boundary layer and the three-dimensional flow field behind the liquid jet are responsible for streamer formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamics of Vortex-Induced-Vibrations of a Slit-Offset Circular Cylinder for Energy Harvesting at Low Reynolds Number

Vortex-Induced Vibrations (VIV) offer a safe, renewable, and environmentally friendly energy source for energy harvesting. To enhance the energy harvesting capability of the circular cylinder-based devices, the authors explore the placement of the normal slit by determining the most effective slit offset location from the cylinder's center. Using the open-source Computational Fluid Dynamics (CFD) program OpenFOAM, a series of numerical simulations are conducted to determine the utility of the slit placement on the circular cylinder and how it influences (positively/negatively) the harvester's performance for a 1-degree-of-freedom (1-DOF) VIV system. The study shows that the slit-cylinder displays the three-branch VIV response (i.e. initial branch (IB), upper branch (UB), and lower branch (LB)). At a Reynolds number of 150, the cylinder with no slit exhibits the two-branch VIV response (i.e. IB and LB) but lacks the upper branch. This is the first time the Upper Branch has been acquired for such low Re flows. The results indicate that adding a normal slit to the middle of the cylinder improves the alternating suction and blowing phenomena. Offsetting it from the center towards the back stagnation point suppresses the VIV and makes it unsuitable for energy harvesting applications. While positioning the slit toward the front stagnation point improves aerodynamic lift, and the cylinder sheds vortex closer to each other. It enhances the amplitude of transverse oscillation and, consequently, the power extraction. In addition, the peak energy transfer ratio for these scenarios is comparable to that of the no-slit case but with a larger range of peak energy transfer ratio values. It makes it suitable for energy harvesting applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flow Characteristics of Elastically Mounted Slit Cylinder at Sub-critical Reynolds Number

The present work numerically investigates vortex-induced vibrations (VIV) of a two-dimensional circular cylinder with an axisymmetric slit at Reynolds number 500. The study examines the effects of slit shape (i.e., converging, diverging, and parallel slits), the effect of slit-area ratios, slit-width, and its angle with the freestream velocity on aerodynamic forces, vibration response, and associated flow characteristics. The results demonstrate that the addition of the slit assists the VIV suppression by adding an extra amount of flow to the main flow. It results in the stabilized wake with the pressure recovery downstream of the cylinder and causes a reduction in the lift force over the cylinder. Also, there exist different shedding patterns associated with different slit shapes. Proper Orthogonal Decomposition (POD) of the flow field suggests that among all the three slits, the parallel slit heavily modifies the flow behind the cylinder by distributing the amount of energy to a large number of modes as compared to other slits (converging slit and diverging slit), where the most of the energy is contained in a couple of modes. Further, increasing the slit angle (for parallel slit) with respect to the freestream increases the slit effectiveness up to a specific value of slit-angle and beyond that starts to affect the VIV suppression adversely. Observations for the effect of slit-width are also reported from the perspective of suppression of VIV.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical analysis of dilute methanol spray flames in vitiated coflow using extended Flamelet Generated Manifold model

The present work focuses on the large eddy simulation (LES) and the study of turbulent dilute methanol spray flames in vitiated coflow using the secondary oxidizer Flamelet Generated Model (FGM). The modified FGM model uses an additional secondary oxidizer parameter in addition to the three other parameters previously used for spray flames progress variable, mixture fraction, and enthalpy. The results for gas phase and droplet properties are validated against the dilute methanol spray flame database for varying fuel injection amounts. The droplet statistics and the lift-off flame heights are accurately captured for all the cases. A proper orthogonal decomposition (POD) of the scalar fluctuating hydroxyl radical (OH) field and velocity-temperature field capture the flame structures in the downstream region of ignition kernels. The detailed POD analysis reveals that the base frequency of the dominant OH field equals that of the dominant vortical structure of 67.3 Hz. The flame propagation happens around these dominant vortical structures because of the less-strained fluid mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal Cooling Enhancements in a Heated Channel Using Flow-Induced Motion

The paper presents the comparative study of the vortex-induced cooling of a heated channel for the four different cross-sections of the rigid cylinder, i.e., circular, square, semi-circular, and triangular, with or without the rigid/flexible splitter plate at the Reynolds number (based on the hydraulic diameter) of 200. The study presents a comprehensive analysis of the flow and thermal performance for all the cases. For flexible plate cases, a partitioned approach is invoked to solve the coupled fluid-structure-convection problem. The simulations show the reduction in the thermal boundary layer thickness at the locations of the vortices resulting in the improved Nusselt number. Further, the thin plate's flow-induced motion significantly increases the vorticity field inside the channel, resulting in improved mixing and cooling. It is observed that the plate-motion amplitude is maximum when the plate is attached to the cylinder with the triangular cross-section. The power requirement analysis shows that the flexible plate reduces the power required to pump the channel's cold fluid. Thus, based on the observations of the present study, the authors recommend using the flexible plate attached to the cylinder for improved convective cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Novel Sharp Interface Immersed Boundary Framework for Viscous Flow Simulations at Arbitrary Mach Number Involving Complex and Moving Boundaries

This work presents a robust and efficient sharp interface immersed boundary (IBM) framework, which is applicable for all-speed flow regimes and is capable of handling arbitrarily complex bodies (stationary or moving). The work deploys an in-house, parallel, multi-block structured finite volume flow solver, which employs a 3D unsteady Favre averaged Navier Stokes equations in a generalized curvilinear coordinate system; while we employ a combination of HCIB (Hybrid Cartesian Immersed boundary) method and GC(Ghost-cell) for solution reconstruction near immersed boundary interface. A significant difficulty for these sharp interface approaches is of handling sharp features/edges of complex geometries. In this study, we observe that apart from the need for robust node classification strategy and higher order boundary formulations, the direction in which the reconstruction procedures are performed plays an important role in handling sharp edges. Taking this into account we present a versatile interface tracking procedure based on ray tracing algorithm and a novel three step solution reconstruction procedure that computes pseudo-normals in the regions where the normal is not well-defined and reconstructs the flow field along those directions. We demonstrate that this procedure enables solver to efficiently handle and accurately represent sharp-edged regions. A fifth-order weighted essentially non-oscillatory (WENO) scheme is used for capturing shock-induced discontinuities and complex fluid-solid interactions with high resolution. The developed IBM framework is applied to a wide range of flow phenomena encompassing all-speed regimes (M=0.001 to M = 2.0). A total of seven benchmark cases (three stationary and four moving bodies) are presented involving various geometries (cylinder, airfoil, wedge) and the predictions are found to be in excellent agreement with the published results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Experimental and Computational Study of a Swirl-Stabilized Premixed Flame

An unconfined strongly swirled flow is investigated for different Reynolds numbers using particle image velocimetry (PIV) and Large Eddy Simulation (LES) with a Thickened Flame (TF) model. Both reacting and non-reacting flow results are presented. In the LES-TF approach, the flame front is resolved on the computational grid through artificial thickening and the individual species transport equations are directly solved with the reaction rates specified using Arrhenius chemistry. Good agreement is found when comparing predictions with the experimental data. Also the predicted RMS fluctuations exhibit a double-peak profile with one peak in the burnt and the other in the un-burnt region. The measured and predicted heat release distributions are in qualitative agreement with each other and exhibit the highest values along the inner edge of the shear layer. The precessing vortex core (PVC) is clearly observed in both the non-reacting and reacting cases. However, it appears more axially-elongated for the reacting cases and the oscillations in the PVC are damped with reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assessment of Soot Formation Models in Lifted Ethylene/Air Turbulent Diffusion Flame

In the present study, soot formation in the turbulent lifted diffusion flame, consisting of ethylene-air is numerically investigated using three different soot modeling approaches and is comprehensively reported. For turbulence-chemistry interaction, Flamelet generated manifold (FGM) model is used. A detailed kinetics is used which is represented through POLIMI mechanism (Ranzi et al. 2012). Soot formation is modeled using two different approaches, semi-empirical two-equation approach and Quadrature methods of moments approach, where both the approaches consider various sub-processes such as nucleation, coagulation, surface growth and oxidation. The radiation heat transfer is taken into account considering four fictitious gasses in conjunction with the weighted-sum-of-gray gas (WSSGM) approach for modeling absorption coefficient. The experimental and earlier published numerical data from Köhler et al. (2012) and Blacha et al. (2011) are used for assessment of different soot modeling approaches. The discrepancies between numerical and experimental data are observed due to under-prediction of OH radicals concentration and poor fuel-air mixing ratios in the vicinity of the fuel jet region leading to early soot formation and the trends are unaffected after invoking radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assessment of turbulence-chemistry interaction models in MILD combustion regime

The present paper reports on the assessment of different turbulence-chemistry interaction closures for modeling turbulent combustion in the Moderate and Intense Low oxygen Dilution combustion regime. 2D RANS simulations have been carried out to model flames issuing from two burners DJHC burner and Adelaide JHC burner which imitate the MILD combustion. In order to model these flames, two different approaches of turbulence-chemistry interaction models are considered; while in the PDF based modeling, two different variants are invoked to understand the applicability of the PDF based models in the MILD regime: one is based on presumed shape PDF approach and the other one is transported PDF approach. A comprehensive study has been carried out by comparing the results obtained from these different models. For the DJHC burner, the computations are carried out for a jet speed corresponding to Reynolds numbers of Re=4100, whereas the Adelaide JHC burner computations are performed for a jet speed corresponding to Reynolds number of Re=10000. The effects of molecular diffusion on the flame characteristics are also studied by using different micro-mixing models. In the case of DJHC burner, it has been observed that the mean axial velocity and the turbulent kinetic energy profiles are in good agreement with the measurements. However, the temperature profiles are over-predicted in the downstream region by both EDC and the PDF based models. In the context of Adelaide JHC burner, the profiles of the temperature and the mass fraction of major species (CH4, H2, O2, H2O, CO2) obtained using LPDF approach are in better agreement with the measurements compared to those obtained using EDC model; although, both the solution approaches fail to capture CO and OH radical profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

CFD-CAA simulation of flow acoustic coupling in a half-dump combustor

This paper reports on the investigation of combustion instabilities in a Methane air nonpremixed half dump combustor for different flow Reynolds number using Computational aero acoustic simulation. In order to simulate the flow physics under turbulent flow conditions, Detached Eddy Simulation in conjunction with generalized eddy dissipation model is adopted, while a CAA formulation based on Lighthills acoustic analogy is used for computing the acoustic field. It is observed that the unsteady pressure signals either predominantly arise from the natural acoustic modes of the duct or the local flow fluctuations in the vortex shedding process downstream of the dump plane, giving rise to different dominant frequencies at different spatial locations at lower Re and a single dominant lockedon frequency at higher Re. The nondimensional numbers (Helmholtz and Strouhal numbers) are used to characterize the duct acoustic modes from the vortex shedding modes. Under reacting flow conditions, unsteady heat release and pressure oscillations are monitored to compute Rayleigh index in order to verify if the instability is driving or damping at different frequency levels. Moreover, the predicted Helmholtz and Strouhal numbers are found to be in excellent agreement with the experimental data available in open literature for wide range of Re.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Characterization of turbulent supersonic flow over a backward facing step through POD

The present work reports on the flow physics of turbulent supersonic flow over backward facing step (BFS) at Mach 2 using LES methodology where the dynamic Smagorinsky model is used for SGS modeling, while POD is invoked to identify the coherent structures present in the flow. The mean data obtained through the computations is in good agreement with the experimental measurements, while the iso-surfaces of Q-criterion at different time instants show the complex flow structures. The presence of counter rotating vortex pair in the shear layer along with the complex shock wave/boundary layer interaction leading to the separation of boundary layer is also evident from the contours of both Q and the modulus of vorticity. Further, the POD analysis reveals the presence of coherent structures, where the first and second modes confirm the vortical structures near the step as well as along the shear layer in the downstream region; while the second, third and fourth modes confirm the presence of vortices along the shear layer due to Kelvin-Helmholtz (K-H) instability. Moreover, POD as well as frequency analysis is extended at different planes to extract the detailed flow features.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamics of upstream flame propagation in a hydrogen-enriched premixed flame

An unconfined strongly swirled flow is investigated to study the effect of hydrogen addition on upstream flame propagation in a methane-air premixed flame using Large Eddy Simulation (LES) with a Thickened Flame (TF) model. A laboratory-scale swirled premixed combustor operated under atmospheric conditions for which experimental data for validation is available has been chosen for the numerical study. In the LES-TF approach, the flame front is resolved on the computational grid through artificial thickening and the individual species transport equations are directly solved with the reaction rates specified using Arrhenius chemistry. Good agreement is found when comparing predictions with the published experimental data including the predicted RMS fluctuations. Also, the results show that the initiation of upstream flame propagation is associated with balanced maintained between hydrodynamics and reaction. This process is associated with the upstream propagation of the center recirculation bubble, which pushes the flame front in the upstream mixing tube. Once the upstream movement of the flame front is initiated, the hydrogen-enriched mixture exhibits more unstable behavior; while in contrast, the CH4 flame shows stable behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effect of grid sensitivity on the performance of wall adapting SGS models for LES of swirling and separating reattaching flows

The present study assesses the performance of the Wall Adapting SGS models along with the Dynamic Smagorinsky model for flows involving separation, reattachment and swirl. Due to the simple geometry and wide application in a variety of engineering systems, the Backward-Facing Step (BFS) geometry and Confined Swirling Flow (CSF) geometry are invoked in the present case. The calculation of the SGS stresses employs three models, namely, the Dynamic Smagorinsky model, the Wall Adapting Local Eddy viscosity (WALE) model and the Dynamic WALE model. For studying the effect of the grid sensitivity, the simulations are performed over two sets of grids with different resolutions based on the non-dimensional wall distance parameter ( y+ ). Grids corresponding to y+ = 70 and y+ = 20 are employed for the subsonic flow over the BFS while grids corresponding to y+ = 40 and y+ = 20 are employed for supersonic flow over the BFS and for confined swirling flow geometry. The validation against the experimental results includes the mean flow fields and the turbulent stresses obtained for each case. The results reveal that for the fine grid (y+ = 20), the near wall eddy viscosity profile for the WALE model is better than both the Dynamic WALE and the Dynamic Smagorinsky model. The difference between the predictions of the coarse and fine grids for Dynamic Smagorinsky and the WALE model is high whereas, the Dynamic WALE model is almost insensitive to the grid resolutions considered for the present case. The mean velocity and pressure values as well as the turbulent quantities predicted by the Dynamic WALE model are closest to the experimental values for all the cases.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effect of Leading Edge Tubercles on Compressor Cascade Performance

Tubercles are modifications to the leading edge of an airfoil in the form of blunt wave-like serrations. Several studies on the effect of tubercles on isolated airfoils have shown a beneficial effect in the post-stall regime, as reduced drag and increased lift, leading to a delay of stall. The prospect of delaying stall is particularly attractive to designers of axial compressors in gas turbines, as this leads to designs with higher loading and therefore higher pressure rise with fewer number of stages. In the present study, experiments were performed on a cascade of airfoils with NACA 65209 profile with different tubercle geometries. The measurements were made over an exit plane using a five-hole probe to compare the cascade performance parameters. Additionally, hot-wire measurements were taken near the blade surface to understand the nature of the flow in the region close to the tubercles. Oil-flow visualization on the cascade end wall reveal the flow through the passage of blades with and without tubercles. For the cascade considered, the estimated stall angle for the best performing set of blades is found to increase up to 8.6° from that of the unmodified blade of 6.0°. Application of such structures in axial compressor blades may well lead to suppression of stall in axial compressors and extend the operating range.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effect of Precursors and Radiation on Soot Formation in Turbulent Diffusion Flame

Soot formation in Delft flame III, a pilot stabilized turbulent diffusion flame burning natural gas/air, is investigated using ANSYS FLUENT by considering two different approaches for soot inception. In the first approach soot inception is based on the formation rate of acetylene, while the second approach considers the formation rate of two and three-ringed aromatics to describe the soot inception [1]. Transport equations are solved for soot mass fraction and radical nuclei concentration to describe inception, coagulation, surface growth, and oxidation processes. The turbulent-chemistry interactions and soot precursors are described by the steady laminar flamelet model (SLFM). Two chemical mechanisms GRI 3.0 [2] and POLIMI [3] are used to represent the effect of species concentration on soot formation. The radiative properties of the medium are included based on the non-gray modeling approach by considering four factious gases; the weighted sum of gray gas (WSGGM) approach is used to model the absorption coefficient. The effect of soot on radiative transfer is modeled in terms of effective absorption coefficient of the medium. A beta probability density function (\b{eta}-PDF) in terms of normalized temperature is used to describe the effect of turbulence on soot formation. The results clearly elucidate the strong effect of radiation and species concentration on soot volume fraction predictions. Due to increase in radiative heat loss with soot, flame temperature decreases slightly. The inclusion of ethylene has less synergic effect than that of both benzene and ethylene. Both cases have less impact on the nucleation of soot. The increase in soot volume fraction with soot-turbulence interaction is in consistence with the DNS predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Flow Past Stationary and Oscillating Airfoil at Low Reynolds number Using Sharp Interface Immersed Boundary Approach

The present study reports on flow past airfoils (stationary and moving) using sharp interface immersed-boundary approach. Non-boundary conforming approach like immersed-boundary method offers a viable alternative over traditional boundary conforming approach by allowing us to model flow past arbitrarily complex shapes, by eliminating the need to re-grid the flow domain as the body exhibits motion. We present flow past a NACA 0012 airfoil at stationary conditions as well as exhibiting pitching motion. Evolution of vortex dynamics and wake structures are presented to show that the developed sharp interface immersed-boundary approach captures the flow physics of dynamic stall accurately. Moving body problems involving immersed-boundary approach usually encounter the issues of spurious oscillations and mass conservation. This is handled through a field extension strategy based on ghost cell approach, which allows for extrapolating the flow field value onto the ghost nodes, ensuring smooth temporal transition as the immersed surface moves through time. The results presented here show excellent agreement with the experimental results found in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of asymmetrically pitching airfoil at high reduced frequency

The expanding application in Micro-Air Vehicles has encouraged many researchers to understand the unsteady flow around a flapping foil at a low Reynolds number. We numerically investigate an incompressible unsteady flow around a two-dimensional pitching airfoil (SD7003) at high reduced frequency in the laminar regime. This study interrogates the effect of different unsteady parameters, namely amplitude (A), reduced frequency (k), Reynolds number (Re), and asymmetry parameter (S) for pitching motion on the force coefficients. The inviscid theoretical model is utilized to calculate the lift coefficient for sinusoidal motion in the viscous regime, and a comparison is made with the numerical results. The theoretical analysis identifies the influence of the non-circulatory lift over circulatory lift at a high reduced frequency. Further, the results indicate that the reduced frequency (k) and asymmetry parameter (S) have a significant impact on the instantaneous and time-averaged force coefficients as well as on the vortex structure in the wake. Finally, the Fast Fourier Transformation (FFT) analysis is carried out over a simulated case with fixed amplitude and Reynolds number for distinct k and S values. The findings confirm that the dominant frequency in the flow (k*) has a direct correlation to the airfoil pitching frequency (k).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of Flow Characteristics inside a Dual Bell Nozzle with and without Film Cooling

In this study, we perform a two-dimensional axisymmetric simulation to assess the flow characteristics and understand the film cooling process in a dual bell nozzle. The secondary stream with low temperature is injected at three different axial locations on the nozzle wall, and the simulations are carried out to emphasize the impact of injection location (secondary flow) on film cooling of the dual bell nozzle. The cooling effect is demonstrated through the temperature and pressure distributions on the nozzle wall or, in-turn, the separation point movement. Downstream of the injection point, the Mach number and temperature profiles document the mixing of the main flow and secondary flow. The inflection region is observed to be the most promising location for the injection of the secondary flow. We have further investigated the effect of Mach number of the secondary stream. The current study demonstrates that one can control the separation point in a dual bell nozzle with the help of secondary injection (Mach number) so that an optimum amount of thrust can be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of flow field of clap & fling motion using immersed boundary coupled lattice Boltzmann method

This paper deals with the investigation of flow field due to clap and fling mechanism using immersed boundary coupled with lattice Boltzmann method. The lattice Boltzmann method (LBM), an alternative to Navier Stokes solver, is used because of its simplicity and computational efficiency in solving complex moving boundary problems. Benchmark problems are simulated to validate the code, which is then used for simulating flow over two elliptic wing of aspect ratio 5 performing clap and fling flapping motion for different flow parameters such as Reynolds number (Re = 75,100,150), advance ratio (J = 10E-3,10E-2,2) and frequency (f = 0.05Hz,0.25Hz). Numerical simulation is able to capture typical low Reynolds number unsteady phenomena such as, &#39;wake vortex wing interaction&#39;, &#39;Kramer effect&#39; and &#39;delayed stall&#39;. The results are both qualitatively and quantitatively consistent with experimental observation. The parametric study involving different combinations of Re, f and J depict distinctly different aerodynamic performances providing physical insights into the flow physics. It is observed that a combination of low f, low J and high Re flow results in better aerodynamic performance. Pronounced lift enhancement via leading edge vortices are obtained in unsteady regime (J<1) compared to quasi-steady regime (J>1). The role of leading edge vortices in enhancing lift are investigated by studying the size and strength of these vortices for different flow conditions. For a given Re, the magnitude of maximum lift coefficient decreases with increasing f irrespective of the value of J; while the same is enhanced with the increasing Re.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of flow structures in a turbulent separating flow using hybrid RANS-LES model

Hybrid Reynolds-averaged Navier Stokes large eddy simulation (RANS LES) methods have become popular for simulation of massively separated flows at high Reynolds numbers due to their reduced computational cost and good accuracy. The current study aims to examine the performance of LES and hybrid RANS-LES model for a given grid resolution. For better assessment and contrast of model performance, both mean and instantaneous flow fields have been investigated. For studying instantaneous flow, proper orthogonal decomposition has been used. Current analysis shows that hybrid RANS-LES is capable of achieving similar accuracy in prediction of both mean and instantaneous flow fields at a very coarse grid as compared to LES. Focusing mostly on the practical applications of computation, most of the attention has been given to the prediction of one-point flow statistics and little consideration has been put to two-point statistics. Here, two-point statistics has been considered using POD to investigate unsteady turbulent flow.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of Mixing Characteristics in Strut Injectors using Modal Decomposition

Effect of large-scale vortical structure on mixing and spreading of shear layer is numerically investigated. Two strut configurations namely Straight & Tapered strut at two convective Mach numbers (Mc = 1.4 & 0.37) for two jet heights (0.6 & 1mm) are investigated. Hydrogen jet is injected through a two-dimensional slot in oncoming coflow at Mach 2. Excellent agreement between simulated and experimental data is witnessed, whereas the instantaneous data reveal the presence of various large-scale structures in the flow field. From the instantaneous field, it becomes apparent that both the geometries have different vortical breakdown locations. It is also noticed that an early onset of vortex breakdown manifests itself into the mixing layer thickness enhancement, effect of which is reflected in overall mixing characteristics. It becomes evident that the shear strength plays an important role in the near field mixing. The higher shear strength promotes the generation of large vortices. The analysis shows that SS-0.6 case offers highest mixing efficiency being dominated by relatively large-scale structures. Eigenmodes obtained through Proper Orthogonal Decomposition (POD) confirm the presence of dominating structures and shed light into the series of events involved in vortex pairing/merging and breakdown. Dynamic modal decomposition (DMD) also strengthens the observation made through the POD.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of NOx in piloted stabilized methane-air diffusion flames using Finite-rate and Infinitely-fast chemistry based combustion models

The present work reports on the numerical investigation of NOx in three turbulent piloted diffusion flames of different levels of extinction. The study involves two-dimensional axisymmetric modeling of combustion in these flames with fairly detailed chemistry, i.e. GRI 3.0 mechanism. The main focus of the study is to analyze the effects of the two different combustion model approaches, such as infinitely fast chemistry based unsteady flamelet and finite rate chemistry based EDC, in predicting the NOx formation in three piloted methane jet flames (Sandia D, E, and F). The EDC approach is able to predict the passive scalar quantities but shows over-prediction in the reactive scalar quantities and NO prediction, while the unsteady flamelet modeling is found to be essential in predicting the accurate formation of slow kinetic species like NOx. The inability of flamelet and EDC approach in capturing localized flame extinction is observed, which lead to an over-prediction of NOx at larger downstream locations. Further, the dominance of NOx formation pathways is investigated in all three flames.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of passive flow control over an airfoil using leading edge tubercles

The present paper numerically investigates flow control over NACA0021 (National Advisory Committee for Aeronautics) airfoil by applying a leading-edge tubercle as a passive flow control device at transitional Reynolds number (Re) regime. The study includes simulations over the unmodified airfoil and modified airfoil operating at Re=120000 using kT-kL-ω based Delayed Detached Eddy Simulation (DDES) model. The calculations include both the pre-stall and post-stall regions, while the performance of the mean flow aerodynamic properties like pressure, lift, skin-friction, and drag coefficients are also investigated in the presence of tubercles. The spanwise flow over the modified airfoil interacts with the primary flow. The turbulent separation gets delayed compared to the unmodified airfoil, thereby significantly improving the stall behavior and flow behavior in the post-stall region. The modified airfoil provides the gradual stall behavior against the steep stall behavior for the unmodified airfoil. To optimize tubercle configuration, the study considers two modified airfoils with different combinations of amplitude and wavelength. The unsteady analysis, including the reduced-order modeling, i.e. 3D Proper orthogonal decomposition (POD), characterizes the dominant vortical structure of the flow to evaluate the improvement in the aerodynamic performance of the modified airfoil. The model decomposition for flow over airfoil helps to provide a deeper understanding of the flow control phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of Shock Wave Interactions involving Stationary and Moving Wedges

The present study investigates the shock wave interactions involving stationary and moving wedges using a sharp interface immersed boundary method combined with a fifth order weighted essentially non oscillatory (WENO) scheme. Inspired by Schardins problem, which involves moving shock interaction with a finite triangular wedge, we study influences of incident shock Mach number and corner angle on the resulting flow physics in both stationary and moving conditions. The present study involves three incident shock Mach numbers (1.3, 1.9, 2.5) and three corner angles (60deg, 90deg, 120deg), while its impact on the vorticity production is investigated using vorticity transport equation, circulation, and rate of circulation production. Further, the results yield that the generation of the vorticity due to the viscous effects are quite dominant compared to baroclinic or compressibility effects. The moving cases presented involve shock driven wedge problem. The fluid and wedge structure dynamics are coupled using the Newtonian equation. These shock driven wedge cases show that wedge acceleration due to the shock results in a change in reflected wave configuration from Single Mach Reflection (SMR) to Double Mach Reflection (DMR). The intermediary state between them, the Transition Mach Reflection (TMR), is also observed in the process. The effect of shock Mach number and corner angle on Triple Point (TP) trajectory, as well as on the drag coefficient, is analyzed in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigation of strut-ramp injector in a Scramjet combustor: effect of strut geometry, fuel and jet diameter on mixing characteristics

The strut based injector has been found to be one of the most promising injector designs for supersonic combustor, offering en-hanced mixing of fuel and air. The mixing and flow field characteristics of the straight (SS) & tapered strut (TS), with fixed ramp an-gle and height at freestream Mach number 2 in conjunction with fuel injection at Mach 2.3 have been investigated numerically and reported. In the present investigation, hydrogen (H2) and ethylene (C2H4) are injected in oncoming supersonic flow from the back of the strut, where jet to freestream momentum ratio is maintained at 0.79 and 0.69 for H2 & C2H4, respectively. The predicted wall static pressure and species mole fractions at various downstream locations are compared with the experimental data for TS case with 0.6 mm jet diameter and found to be in good agreement. Further, the effect of jet diameter and strut geometry on the near field mixing in strut ramp configuration is discussed for both the fuel. The numerical results are assessed based on various parameters for the performance evaluation of different strut ramp configurations. The SS configuration for both the injectant is found to be an optimum candidate, also it is observed that for higher jet diameter larger combustor length is required to achieve satisfactory near field mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Eddy Simulation of a Premixed Bunsen flame using a modified Thickened-Flame model at two Reynolds number

A modified Thickened Flame (TF) model based on Large Eddy Simulation methodology is used to investigate premixed combustion and the model predictions are evaluated by comparing with the piloted premixed stoichiometric methane-air flame data (Chen et al., 1996) for Reynolds numbers Re = 24,000 (flame F3) and Re=52,000 (flame F1). The basic idea of Thickened-Flame approach is that the flame front is artificially thickened to resolve on the computational LES grid while keeping the laminar flame speed ( ) constant. The artificially thickening of the flame front is obtained by enhancing the molecular diffusion and decreasing the pre-exponential factor of the Arrhenius law. Since the flame front is artificially thickened, the response of the thickened flame to turbulence is affected and taken care of by incorporating an efficiency function (E) in the governing equations. The efficiency function (E) in the modified TF model is proposed based on the direct numerical simulations (DNS) data set of flame-vortex interactions (Colin et al., 2000). The predicted simulation results are compared with the experimental data and with computations reported using RANS based probability distribution function (PDF) modeling approach (Lindstedt, R. P. and Vaos, E. M., 2006, Transported PDF modeling of high-Reynolds-number premixed turbulent flames, Combustion and Flame, 145, 495) and RANS based G-equation approach (Herrmann, M., 2006). It is shown that the results with the modified TF model are generally in good agreement with the data, with the TF predictions consistently comparable to the PDF model predictions and superior to the results with the G-equation approach.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Eddy Simulation of MILD combustion using PDF based turbulence-chemistry interaction models

The present work reports on LES of turbulent reacting flow of the Delft-Jet-in-Hot-Coflow (DJHC) burner, emulating MILD combustion, with transported PDF based combustion models using ANSYS FLUENT 13.0. Two different eddy viscosity models for LES (Dynamic Smagorinsky and Kinetic Energy Transport) along with two solution approaches for PDF transport equation, i.e. Eulerian and Lagrangian, have been used in the present study. Moreover, the effects of chemical kinetics and the micro-mixing models have also been investigated for two different fuel jet Reynolds number (Re = 4,100 and Re = 8,800). The mean velocity and turbulent kinetic energy predicted by the different models are in good agreement with experimental data. Both the composition PDF models predict an early ignition resulting in higher radial mean temperature predictions at burner exit. The models, however, correctly predict the formation mechanism of ignition kernels and the decreasing trend of the lift-off height with increasing jet Reynolds number, as observed experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large Eddy Simulation of Premixed Combustion with a Thickened Flame Approach

A Thickened Flame (TF) modeling approach is combined with a Large Eddy Simulation (LES) methodology to model premixed combustion and the accuracy of these model predictions is evaluated by comparing with the piloted premixed stoichiometric methane-air flame data of Chen et al. [Combust. Flame 107 (1996) 233-244] at a Reynolds number Re = 24,000. In the TF model, the flame front is artificially thickened to resolve it on the computational LES grid and the reaction rates are specified using reduced chemistry. The response of the thickened flame to turbulence is taken care of by incorporating an efficiency function in the governing equations. The efficiency function depends on the characteristics of the local turbulence and on the characteristics of the premixed flame such as laminar flame speed and thickness. Three variants of the TF model are examined: the original Thickened Flame model, the Power-law flame wrinkling model, and the dynamically modified TF model. Reasonable agreement is found when comparing predictions with the experimental data and with computations reported using a probability distribution function (PDF) modeling approach. The results of the TF model are in better agreement with data when compared with the predictions of the G-equation approach

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Leading Edge Vortex Dynamics of Airfoils, Pitching Continuously at High Amplitudes

Airfoils pitching in the stalled regime have been of keen interest in recent years due to their desirable aerodynamic force characteristics. In this study, we are numerically investigating the unsteady flow past a NACA 0012 airfoil under sinusoidally pitching motion, using a finite volume based sharp interface immersed boundary solver. The flow is investigated in the low Reynolds number regime (Re=3000) for reduced frequencies of 0.1 and 0.5, at three different pivot locations (c/3, c/2 and 2c/3 from the leading edge). The airfoil is subjected to sinusoidal oscillations with its incidence angle varying from 15o to 45o. Leading edge vortices (LEVs) that are formed during the pitching motion dictate the transient aerodynamic characteristics. The flow field data is used to identify individual LEVs in the flow field. The spatio-temporal evolution of LEVs in terms of their strengths are traced throughout the pitching cycle to obtain a quantitative estimate of its evolution. The effect of pivot location and pitching frequency on the vortex evolution and aerodynamic forces is investigated. Change in pitching frequency is found to have a drastic effect on the vortex dynamics. LEV formation is delayed to higher fractional times and LEVs are found to interact and merge into each other at higher frequencies. The maximum lift and drag are also found to increase with pitching frequency. Pivot location is found to have a higher influence on the magnitude and phase of aerodynamic coefficients at higher frequencies. However, at both frequencies, moving the pivot axis aftward is found to delay LEV evolution and cause a phase lag on the aerodynamic forces.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Modal Decomposition of Turbulent Supersonic Cavity

Self-sustained oscillation in Mach 3 supersonic cavity with a length-to-depth ratio of 3 is investigated using wall modeled Large Eddy Simulation (LES) methodology for ReD = 3.39*10^5. The unsteady data obtained through computation is utilized to investigate the spatial and temporal evolution of the flow field, especially second invariant of velocity tensor; while the phase averaged data is analyzed over a feedback cycle to study the spatial structures. This analysis is accompanied by the Proper Orthogonal Decomposition (POD) data, which reveals the presence of discrete vortices along the shear layer. The POD analysis is performed in both the spanwise and streamwise planes to extract the coherence in flow structures. Finally, Dynamic Mode Decomposition (DMD) is performed on the data sequence to get the dynamic information and deeper insight into the self-sustained mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Investigation of Coaxial GCH4/LOx Combustion at Supercritical Pressures

This article aims to numerically investigate the combustion phenomenon of coaxial gaseous CH4 LOx at supercritical pressures. The choice of turbulence model, real gas model, and chemical kinetics model are the critical parameters in numerical simulations of cryogenic combustion at high pressure. At this supercritical operating pressure, the ideal gas law does not remain valid for such cases. Therefore, we have systematically carried out a comparative study to analyze the importance of real gas models, turbulence parameters, and chemical kinetics at such conditions. The comparison of real gas models with the NIST database reveals better conformity of SRK (Soave Redlich Kwong Equation of State (EoS)) model predictions with the database. Further, the computed results indicate that the Standard k-e turbulence model with modified constant captures the better flame shape and temperature peak position compared to other RANS based turbulence models while invoking the non-premixed steady b-PDF flamelet model for simulating the combustion process. Furthermore, a comparative study comparing two different chemical kinetics models indicates that the reduced Jones Lindstedt mechanism can accurately predict the flame characteristics with the least computational cost. Finally, we have studied the effect of chamber pressure and LOx inlet temperature on the flame characteristics. The flame characteristics exhibit a strong sensitivity towards the chamber pressure due to the weakening of the pseudo-boiling effect with an increase in pressure. As a consequence of lower turbulent rates of energy and mass transfer through the transcritical mixing layer, the flame spreading becomes narrower at elevated pressure and temperature, thereby yielding an increased flame length at transcritical conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Investigation of Flow Structures Around a Cylindrical Afterbody Under Supersonic Condition

Large Eddy Simulation (LES) with dynamic Smagorinsky model has been applied to numerically investigate the complicated flow structures that evolve in the near wake of a cylindrical after body aligned with a uniform Mach 2.46 flow. Mean flow field properties obtained from numerical simulations, such as axial velocity, pressure on base surface, have been compared with the experimental measurements as well as with other published results. It has been found that standard k-epsilon model fails to predict the flow properties in the recirculation region where better agreement has been observed between the data obtained from LES and experimental measurements. Flow Statistics like turbulent kinetic energy and primary Reynolds stress have also been calculated and compared with the results obtained from experiments in order to quantitatively assess the ability of LES technique to predict the turbulence properties of flow field in the highly compressible shear layer region. The data obtained from LES has been further analyzed to understand the evolution of coherent structures in the flow field. Proper Orthogonal Decomposition (POD) of the data obtained from central plane in the wake region has been performed in order to reveal the most energetic structures present in the flow field.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical investigation of MILD combustion using Multi-Environment Eulerian Probability Density Function modeling

In the present paper, the flames imitating Moderate and Intense Low Oxygen Dilution (MILD) combustion are studied using the Probability Density Function (PDF) modeling approach. Two burners which imitate MILD combustion are considered for the current study: one is Adelaide Jet-in-Hot-Coflow (JHC) burner and the other one is Delft-Jet-In-Hot-Coflow (DJHC) burner. 2D RANS simulations have been carried out using Multi-environment Eulerian Probability Density Function (MEPDF) approach along with the Interaction-by-Exchange-with-Mean (IEM) micro-mixing model. A quantitative comparison is made to assess the accuracy and predictive capability of the MEPDF model in the MILD combustion regime. The computations are performed for two different jet speeds corresponding to Reynolds numbers of Re=4100 and Re=8800 for DJHC burner, while Re=10000 is considered for the Adelaide burner. In the case of DJHC burner, for Re=4100, it has been observed that the mean axial velocity profiles and the turbulent kinetic energy profiles are in good agreement with the experimental database while the temperature profiles are slightly over-predicted in the downstream region. For the higher Reynolds number case (Re=8800), the accuracy of the predictions is found to be reduced. Whereas in the case of Adelaide burner, the computed profiles of temperature and the mass fraction of major species (CH4, H2, N2, O2) are found to be in excellent agreement with the measurements while the discrepancies are observed in the mass fraction profiles of CO2 and H2O. In addition, the effects of differential diffusion are observed due to the presence of H2 in the fuel mixture.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Investigation of Soot Formation in Turbulent Diffusion Flame with Strong Turbulence Chemistry Interaction

The present work is aimed at examining the ability of different models in predicting soot formation in Delft flame III, which is a non-premixed pilot stabilized natural gas flame. The turbulence-chemistry interactions are modeled using a laminar steady flamelet model (SLFM). One step and two step models are used to describe the formation, growth, and oxidation of soot particles. One Step is an empirical model which solves the soot mass fraction equation. The two-step models are semi empirical models, where the soot formation is modelled by solving the governing transport equations for the soot mass fraction and normalized radical nuclei concentration. The effect of radiative heat transfer due to gas and soot particulates is included using P1 approximation. The absorption coefficient of the mixture is modeled using the weighted sum of gray gases model (WSGGM). The turbulence-chemistry interaction effects on soot formation are studied using a single variable probability density function (PDF) in terms of a normalized temperature or mixture fraction. The results shown in this work clearly elucidate the effect of radiation and turbulence-chemistry interaction on soot formation. The soot volume fraction decreases with the introduction of radiation interactions, which is consistence with the theoretical predictions. It has also been observed in the current work that the soot volume fraction is sensitive to the variable used in the PDF to incorporate the turbulence interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical simulation of Delft-jet-in-hot-coflow (DJHC) flames using the Eddy Dissipation Concept model for turbulence-chemistry interaction

In this paper we report results of a numerical investigation of turbulent natural gas combustion for jet in a coflow of lean combustion products in the Delft-Jet-in-Hot-Coflow (DJHC) burner which emulates MILD (Moderate and Intense Low Oxygen Dilution) combustion behavior. The focus is on assessing the performance of the Eddy Dissipation Concept (EDC) model in combination with two-equation turbulence models and chemical kinetic schemes for about 20 species (Correa mechanism and DRM19 mechanism) by comparing predictions with experimental measurements. We study two different flame conditions corresponding to two different oxygen levels (7.6% and 10.9% by mass) in the hot coflow, and for two jet Reynolds number (Re=4100 and Re=8800). The mean velocity and turbulent kinetic energy predicted by different turbulence models are in good agreement with data without exhibiting large differences among the model predictions. The realizable k-e model exhibits better performance in the prediction of entrainment. The EDC combustion model predicts too early ignition leading to a peak in the radial mean temperature profile at too low axial distance. However the model correctly predicts the experimentally observed decreasing trend of lift-off height with jet Reynolds number. A detailed analysis of the mean reaction rate of the EDC model is made and as possible cause for the deviations between model predictions and experiments a low turbulent Reynolds number effect is identified. Using modified EDC model constants prediction of too early ignition can be avoided. The results are weakly sensitive to the sub-model for laminar viscosity and laminar diffusion fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Study of Bifurcating Flow through Sudden Expansions: Effect of divergence and geometric asymmetry

A numerical study of laminar flow through symmetric and slightly asymmetric sudden expansion, of expansion ratio 1:3, in channels with increasing cross section, is carried out using two different approaches - Conventional CFD and Lattice Boltzmann Method. The effect of divergence of walls of the channel, after a sudden expansion, on the symmetry of flow and recirculation is studied for various Reynolds numbers. It is seen that the angles of the walls play an active role in disrupting the symmetry of flow. For non-parallel walls, the symmetry breaking bifurcation phenomenon no longer exists and the loss of symmetry is a gradual process. The effect of asymmetry of geometry on flow is also studied by considering two types of asymmetry - First type is at the plane of expansion where the steps on either side are of unequal heights, while the second one deals the walls of the channel are at a different inclination with the direction of inflow. The present study reveals that small asymmetry impedes the sharp transition and thus smoothens the pitchfork bifurcation. Increased asymmetry finally leads to the unique solution of the governing Equations, indicating complete loss of bifurcation pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical study of flow physics in supersonic base-flow with mass bleed

Large Eddy Simulation (LES) with dynamic sub-grid scale eddy viscosity model has been applied to numerically investigate the evolution of complicated flow structures in supersonic base flow with mass bleed. Mean flow properties obtained from numerical simulations, such as axial velocity, pressure on the base surface, have been compared with the experimental measurements to show that LES is able to predict the mean flow properties with acceptable accuracy. The data obtained from LES has been further analyzed to understand the evolution of coherent structures in the flow field. Periodical shedding of vortical structures from the outer shear layer has been observed and it has also been found that this vortex shedding is associated with the flapping of the outer shear layer. The frequency of flapping of the outer shear layer has been found out and the phase-averaged streamlines have been analyzed to further study the evolution of vortical structures associated with this flapping. The phase-averaged streamline plots clearly elucidate the evolution of vortical structures along the outer shear layer. Further, the study of these structures is investigated by performing Proper Orthogonal Decomposition (POD) analysis of the data, obtained along the central plane in the wake region. The POD results also seem to agree well with the observations made in the phase averaged streamline plots, as the concentrated energy and enstropy are observed in the outer shear layer with fewer POD modes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical Study of Trailing and Leading Vortex Dynamics in a Forced Jet with Coflow

The interaction of trailing and leading vortex structures in low Reynolds number forced (varicose) jet with coflow is studied using Direct Numerical Simulation (DNS), where we report the dynamics of these vortices for a range of varying parameters. In a circular forced jet with Strouhal number 0.2, continuous formation of these vortices is observed which undergo pairing, tearing and disintegration. The forced jet is analyzed for varying frequency of forced perturbations, coflow temperature, turbulence intensity, momentum thickness, coflow intensity and Reynolds number at a constant forcing amplitude. The observation reveals that any reduction in momentum thickness increases the circulation as well as energy capacity of the leading and the trailing vortices, thereby enhancing the vortex pairing time. However, the hot coflow increases the strength of the leading vortex to such an extent that the generated trailing vortex is weak and dissipates even before it undergoes vortex pairing, while an incomplete disintegration of the leading vortex occurs in this case. Moreover, the strength of leading and trailing vortices is found to decrease with an increase in coflow power, thereby leading to a weaker vortex pairing. The decrease in Reynolds number also reduces the intensity of trailing vortices; whereas, at a lower Reynolds Number (Re) the trailing vortices dissipate soon after they form. The findings are achieved based on the analysis done by Proper Orthogonal Decomposition (POD), Dynamic Mode Decomposition (DMD), fast Fourier transform (FFT), Q-criterion and vortex tracking method to understand the dynamic interaction and behavior of the vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Parametric study of upstream flame propagation in hydrogen-enriched premixed combustion: effects of swirl, geometry and premixedness

The effect of swirl, premixedness and geometry has been investigated for hydrogen enriched premixed flame using Large Eddy Simulation (LES) with a Thickened Flame (TF) model. Swirl strength has been varied to study the effects of swirl on flame behavior in a laboratory-scale premixed combustor operated under atmospheric conditions. In addition, the levels of premixedness and geometry have also been changed to study the role of these quantities on flame behavior. The turbulent flow field and the chemistry are coupled through TF model. In the LES-TF approach, the flame front is resolved on the computational grid through artificial thickening and the individual species transport equations are directly solved with the reaction rates specified using Arrhenius chemistry. Good agreement is found when comparing predictions with the published experimental data including the predicted RMS fluctuations. Also, the results show that higher swirl strength and increase in level of premixedness make the system more susceptible to upstream flame movement due to higher combustibility of hydrogen, which increases the reaction along the flame front, thereby raises temperature in the reaction zone and leads to combustion induced vortex breakdown (CIVB). Moreover, upstream flame movement is always observed at higher swirl strength irrespective of level of premixedness and burner geometry, whereas the premixed systems exhibit stable behavior while operating at low swirl.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Passive control of the onset of vortex shedding in flow past a circular cylinder using Slit

The present study investigates the passive flow control phenomena over a two-dimensional circular cylinder using numerical simulations in the laminar regime. The aim is to explore one of the passive control techniques, which involves the introduction of a slit to the geometry of the cylinder. The two parameters, slit width ratio S/D (slit width/diameter) and slit angle (measured with respect to the incoming flow direction), play an essential role in determining the trend of critical Reynolds number (Rec). Most of the analysis invokes flow visualization and saturation amplitude methods to obtain the critical Reynolds number (indicative of the onset of vortex shedding) for different cases. Further, Hopf bifurcation analysis using Stuart-Landau equation and global stability analysis confirm the accuracy and consistency in the predicted solutions. The additional amount of flow through slit increases the pressure downstream of the cylinder, which consequently leads to an increase in Rec of the modified cylinder. The critical Reynolds number increases with S/D of the modified cylinder at 0 deg slit angle (as an additional amount of flow grows with S/D). The critical Reynolds number shows an increasing trend with the slit angle in the range of S/D (0.05-0.15) as the fluctuation intensity reduces with slit angle in this range. For S/D= 0.15-0.25, the extra amount of flow through slit induces the instability in wake which causes a decrease in Rec with slit angle. A correlation is obtained that estimates the critical Reynolds number for a given S/D and slit angle.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pore-scale Simulation of Shear-Thinning Fluid Flow using Lattice Boltzmann Method

Present work attempts to identify the roles of flow- and geometric-variables on the scaling factor which is a necessary parameter for modeling the apparent viscosity of non-Newtonian fluid in porous media. While idealizing the porous media microstructure as arrays of circular and square cylinders, present study uses multi-relaxation time lattice Boltzmann method to conduct pore-scale simulation of shear thinning non-Newtonian fluid flow. Variation in the size and inclusion ratio of the solid cylinders generates wide range of porous media with varying porosity and permeability. Present study also used stochastic reconstruction technique to generate realistic, random porous microstructures. For each case, pore-scale fluid flow simulation enables the calculation of equivalent viscosity based on the computed shear rate within the pores. It is observed that the scaling factor has strong dependence on porosity, permeability, tortuosity and the percolation threshold, while approaching the maximum value at the percolation threshold porosity. Present investigation quantifies and proposes meaningful correlations between the scaling factor and the macroscopic properties of the porous media.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Prediction of Separation Induced Transition on Thick Airfoil Using Nonlinear URANS Based Turbulence Model

Most of the turbulence models in practice are based on the assumption of a linear relation between Reynolds stresses and mean flow strain rates which generally provides a good approximation in case of attached and fully turbulent flows. A two dimensional numerical study has been carried out over NACA 0021 with k-ωSST model with non-linear correction at Re = 120,000 for various angles of attack which experiences the formation of a laminar separation bubble (LSB). A correct prediction of LSB requires an accurate resolution of anisotropy in Reynolds stresses. For comparison with other linear models, the simulations are also performed with k-kl-ω, k-ωSST and Spalart Allmaras. The performance of these models is assessed through aerodynamic lift, drag, pressure and friction coefficients. It is found that the non-linear k-ωSST and k-kl-ωtransition model provide comparable quality of prediction in lift and drag coefficients (in spite of the fact that non-linear k-ωSST involves solving less number of transport equation than the transition model) as observed in the experiments whereas k-ωSST and SA models under predict the drag coefficient value at low angle of attack due to inability to capture the separation induced transition. It is also observed that the location of laminar separation bubble is captured accurately when nonlinear or transition model is used as opposed to the SA or linear SST models, which lack in the ability to predict the same.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Role of Jet Spacing and Strut Geometry on the formation of Large Scale Structures & Mixing Characteristics

The present study primarily focuses on the effect of jet spacing and strut geometry on the evolution and structure of the large-scale vortices which play a key role in mixing characteristics in turbulent supersonic flows. Numerically simulated results corresponding to varying parameters such as strut geometry and jet spacing (Xn=nDj such that n=2, 3 & 5) for square jet of height Dj=0.6 mm is presented in the current study, while the work also investigates the presence of the local quasi-two-dimensionality for the X2(2Dj) jet spacing; however, the same is not true for higher jet spacing. Further, the tapered strut (TS) section is modified to the straight strut (SS) for investigation, where the remarkable difference in flow physics is unfolded between the two configurations for similar jet spacing (X2: 2Dj). The instantaneous density and vorticity contours reveal the structures of varying scales undergoing different evolution for the different configurations. The effect of local spanwise rollers is clearly manifested in the mixing efficiency and the jet spreading rate. SS configuration exhibits excellent near field mixing behavior amongst all the arrangements. However, in case of TS cases, only X2(2Dj) configuration performs better due to the presence of local spanwise rollers. The qualitative and quantitative analysis reveals that near-field mixing is strongly affected by the two-dimensional rollers, while the early onset of wake mode is another crucial parameter to have improved mixing. Modal decomposition performed for SS arrangement sheds light into the spatial and temporal coherence of the structures, where the most dominant structures are found to be the Von-Karman street vortices in the wake region.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Simulation of Coupled Heat and Mass Transport With Reaction in PEM Fuel Cell Cathode using Lattice Boltzmann Method

Fluid, heat and species transport and oxygen reduction in the cathode of a PEM fuel cell are simulated using multi relaxation lattice Boltzmann method. Heat generation due to oxygen reduction and its effects on transport and reaction are considered. Simulations for various cell operating voltages, temperatures and flow rates with various values of porous media properties, namely, permeability, porosity, and effective porous media diffusion coefficient, are performed to study transport and operating characteristics in the electrode. It is seen that maximum output power density achievable is limited by the mass transport rate. A small increase in current density is obtained by increasing the operating temperature. However, this results in an increase in the rate of heat generation. Permeability and porosity of the gas diffusion layer do not show a significant impact on the performance in the range of values presently simulated. Higher permeability, in turn, resulted in enhancement of thermal gradients in the porous layer. A significant increase in the maximum current density obtainable is observed with increase in flow rates. The higher convection associated with high flow rate facilitates better transport of species and heat at the catalyst layer resulting in larger current density with a lesser chance of hotspot formation. Increased species diffusion coefficient also resulted in increasing the power output substantially. In addition, the fuel utilization is also improved at high diffusion rates in the porous media. The study analyses and shows the impact of various operating and material parameters affecting the performance of a PEM fuel cell with special attention on enhancing the maximum power density attainable.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stabilization of lifted hydrogen jet diffusion flame in a vitiated co-flow: effects of jet and coflow velocities, coflow temperature and mixing

The present paper reports on the numerical investigation of lifted turbulent jet flames with H2/N2 fuel issuing into a vitiated coflow of lean combustion products of H2/air using conditional moment closure method (CMC). A 2D axisymmetric formulation has been used for the predictions of fluid flow, while CMC equations are solved with detailed chemistry to represent the turbulence-chemistry interaction. Simulations are carried out for different coflow temperatures, jet and coflow velocities in order to investigate the impact on the flame lift-off height as well as on the flame stabilization. Furthermore, the role of conditional velocity models on the flame has also been investigated. In addition, the effect of mixing is investigated over a range of coflow temperatures and the stabilization mechanism is determined from the analysis of the transport budgets. It is found that the lift-off height is highly sensitive to the coflow temperature, while the predicted lift-off height using the mixing model constant, i.e., CΦ=4, is found to be the closest to the experimental results. For all the coflow temperatures, the balance is found between the chemical, axial convection and molecular diffusion terms while the contribution from axial and radial diffusion is negligible, thus indicating auto-ignition as the flame stabilization mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Vortex-induced Vibrations of a Confined Circular Cylinder for Efficient Flow Power Extraction

A simple method to increase the flow power extraction efficiency of a circular cylinder, undergoing vortex-induced vibration (VIV), by confining it between two parallel plates is proposed. A two-dimensional numerical study was performed on VIV of a circular cylinder inside a parallel plate channel of height H at Reynolds number 150 to quantify the improvement. The cylinder is elastically mounted with a spring such that it is only free to vibrate in the direction transverse to the channel flow and has a fixed mass ratio (m*) of 10. The energy extraction process is modelled as a damper, with spatially constant damping ration ((), attached to the cylinder. The simulations are performed by varying the reduced velocity for a set of fixed mass-damping (α = m*() values ranging between 0 to 1. The blockage ratio (b = D/H) is varied from 0.25 to 0.5 by changing the channel height. The quasi-periodic initial branch found for the unconfined cylinder shrinks with the increasing blockage. The extracted power is found to increase rapidly with the blockage. For maximum blockage (b = 0.2), the maximum flow power extracted by the cylinder is an order of magnitude larger as compared to what it would extract in an open domain with free stream velocity equal to the channel mean velocity. The optimal mass-damping (αc ) for extracting maximum power is found to lie between 0.2 to 0.3. An expression is derived to predict the maximum extracted power from the undamped response of a confined/unconfined cylinder. With the assumption αc = 0.25, the derived expression can predict the maximum power extraction within +-20% of the actual values obtained from present and previous numerical and experimental studies.

0 citations0 reviewsNo code linkedOpen access