Faculty Publications
Permanent URI for this communityhttps://idr.nitk.ac.in/handle/123456789/18736
Publications by NITK Faculty
Browse
19 results
Search Results
Item Effect of CO2 based natural circulation loop for low temperature applications: CFD analysis(Toronto Metropolitan University, 2019) Wahidi, T.; Nagrani, P.P.; Yadav, A.K.Natural circulation loop (NCL) is a simple and economical heat transfer device in which flow occurs due to the buoyancy effect caused by thermally generated density gradient. In the present study, computational fluid dynamics (CFD) analyses are carried out to emphasize on the fluid ow and heat transfer characteristics of carbon dioxide (CO2 ) based NCL at low temperature (-38°C to 12°C). Studies are conducted in a three-dimensional (3-D) CFD model of NCL at different heat inputs i.e., 100W, 250W, 350W and 500W by keeping the loop fluid at pressure of 50 bar. Methanol is used as coolant in the heat exchanger at a fixed mass flow rate. Effect of loop operating pressure 50 bar on system performance is also investigated. Result are presented in the form of heat transfer rate, pressure drop, Reynolds number (Re) and temperature. Obtained results are validated with available correlations in the form of non-dimensional numbers, and found in good agreement. © 2019, Toronto Metropolitan University. All rights reserved.Item On the reliability of eddy viscosity based turbulence models in predicting turbulent flow past a circular cylinder using URANS approach(2012) Rajani, B.N.; Kandasamy, A.; Majumdar, S.Turbulent flow past circular cylinder at moderate to high Reynolds number has been analysed employing an secondorder time accurate pressure-based finite volume method solving two-dimensional Unsteady Reynolds Averaged Navier Stokes (URANS) equations for incompressible flow, coupled to eddy-viscosity based turbulence models. The major focus of the paper is to test the capabilities and limitations of the present turbulence model-based 2D URANS procedure to predict the phenomenon of Drag Crisis, usually manifested in reliable measurement data, as a sharp drop in the mean drag coefficient around a critical Reynolds number. The computation results are compared to corresponding measurement data for instantaneous aerodynamic coefficients and mean surface pressure and skin friction coefficients. Turbulence model-based URANS computations are in general found to be inadequate for correct prediction of the mean drag coefficients, the Strouhal number and also the coefficients of maximum fluctuating lift over the range of flow Reynolds number varying from 10 4 to 10 7.Item Transient analysis of subcritical/supercritical carbon dioxide based natural circulation loops with end heat exchangers: Numerical studies(Elsevier Ltd, 2014) Yadav, A.K.; Ram Gopal, M.; Bhattacharyya, S.Transient analysis of carbon dioxide based natural circulation loop (NCL) with end heat exchangers has been carried out. Subcritical and supercritical phases of CO2 are considered with operating pressures in the range of 50-100 bar for an operating temperature range of 323 K to 363 K. Studies are carried out for various loop tilt angles, different initial conditions, and different water mass flow rates. Results: are obtained for various inlet temperatures of water in the hot heat exchanger while keeping the inlet temperature of cooling water in the cold heat exchanger fixed. Effect of tilting the loop in XY and YZ planes on transient as well as steady state behaviour of loop are also studied. Validation of simulation results against experimental and numerical results reported in the literature in terms of modified Grashof number (Grm) and Reynolds number (Re) show good agreement. © 2014 Elsevier Ltd. All rights reserved.Item Effect of tilt angle on subcritical/supercritical carbon dioxide-based natural circulation loop with isothermal source and sink\(American Society of Mechanical Engineers (ASME) infocentral@asme.org, 2016) Yadav, A.K.; Ram Gopal, M.R.; Bhattacharyya, S.In recent years, a growing popularity of carbon dioxide (CO2) as a secondary fluid has been witnessed in both forced as well as in natural circulation loops (NCLs). This may be attributed to the favorable thermophysical properties of CO2 in addition to the environmental benignity of the fluid. However, an extensive literature review shows that studies on CO2-based NCLs are very limited. Also, most of the studies on NCLs do not consider the three-dimensional variation of the field variables. In the present work, threedimensional computational fluid dynamics (CFD) models of a NCL with isothermal source and sink have been developed to study the effect of tilt angle in different planes. Studies have been carried out employing subcritical (liquid and vapor) as well as supercritical phase of CO2 as loop fluid at different operating pressures and temperatures. Results are obtained for a range of tilt angles of the loop, and a significant effect is observed on heat transfer, mass flow rate, and stability of the loop. It was also found that changing the orientation of the loop could be an elegant and effective solution to the flow instability problem of NCLs.Item Effects of the Reynolds number on two-dimensional dielectrophoretic motions of a pair of particles under a uniform electric field(Korean Society of Mechanical Engineers, 2016) Kang, S.; Mannoor, M.; Maniyeri, R.This paper presents two-dimensional direct numerical simulations to explore the effect of the Reynolds number on the Dielectrophoretic (DEP) motion of a pair of freely suspended particles in an unbounded viscous fluid under an external uniform electric field. Accordingly, the electric potential is obtained by solving the Maxwell’s equation with a great sudden change in the electric conductivity at the particle-fluid interface and then the Maxwell stress tensor is integrated to determine the DEP force exerted on each particle. The fluid flow and particle movement, on the other hand, are predicted by solving the continuity and Navier-Stokes equations together with the kinetic equations. Numerical simulations are carried out using a finite volume approach, composed of a sharp interface method for the electric potential and a direct-forcing immersed-boundary method for the fluid flow. Through the simulations, it is found that both particles with the same sign of the conductivity revolve and eventually align themselves in a line with the electric field. With different signs, to the contrary, they revolve in the reverse way and eventually become lined up at a right angle with the electric field. The DEP motion also depends significantly on the Reynolds number defined based on the external electric field for all the combinations of the conductivity signs. When the Reynolds number is approximately below Recr ? 0.1, the DEP motion becomes independent of the Reynolds number and thus can be exactly predicted by the no-inertia solver that neglects all the inertial and convective effects. With increasing Reynolds number above the critical number, on the other hand, the particles trace larger trajectories and thus take longer time during their revolution to the eventual in-line alignment. © 2016, The Korean Society of Mechanical Engineers and Springer-Verlag Berlin Heidelberg.Item Numerical analysis of the buckling and recuperation dynamics of flexible filament using an immersed boundary framework(Elsevier B.V., 2019) Kanchan, M.; Maniyeri, R.The dynamics of flexible filaments in viscous shear flow is of interest to biologists and engineers in a wide variety of applications involving folding and unfolding sequence of long-chain biomolecules like DNA, non-motile sperm and microalgae. It is also helpful in understanding the deformation of natural and synthetic fibers which can be applied in areas such as biotechnology. In the present work, deformation and migration behavior of non-motile unicellular phytoplankton diatoms subjected to viscous shear flow are considered. These unicellular diatoms develop into colonies which are made up of linked chains. The complex fluid-structure interaction is solved by developing a two-dimensional numerical model with an immersed boundary framework. The simulation consists of suspending an elastic filament mimicking a diatom chain in a shear flow at low Reynolds number. The governing continuity and Navier–Stokes equations are solved on a Cartesian grid arranged in a staggered manner. A forcing term is added to the momentum equation that incorporates the presence of flexible filament in the fluid domain. The discretization of the governing equation is based on a finite volume method, and a SIMPLE algorithm is used to compute pressure and velocity. A computer code is developed to perform numerical simulations, and the model is first verified with the deformation study of a tethered flexible filament in uniform fluid flow. Next, the shape deformations for flexible filament placed freely in shear flow are compared with the studies of previous researchers. Further, the present results are validated with Jeffery's equation for particles immersed in shear flow along with classification plot for filament orbit regimes. All of these comparisons provide a reasonable validity for the developed model. The effect of bending rigidity and shear rate on the deformation and migration characteristics is ascertained with the help of parametric studies. A non-dimensional parameter called Viscous Flow Forcing value (VFF) is calculated to quantify the parametric results. An optimum Viscous Flow Forcing value is determined which indicates the transition of filaments exhibiting either a recuperative (regaining original shape past deformation) or non-recuperative (permanently deformed) behavior. The developed model is successful in capturing fluid motion, diatom buckling, shape recurrences and recuperation dynamics of diatom chains subjected to shear flow. Further, the developed computational model can successfully illustrate filament-fluid interaction for a wide variety of similar problems. © 2019 Elsevier Inc.Item Numerical simulation of flow in a wavy wall microchannel using immersed boundary method(Bentham Science Publishers, 2020) Kanchan, M.; Maniyeri, R.Background: Fluid flow in microchannels is restricted to low Reynolds number regimes and hence inducing chaotic mixing in such devices is a major challenge. Over the years, the Immersed Boundary Method (IBM) has proved its ability in handling complex fluid-structure interaction prob-lems. Objectives: Inspired by recent patents in microchannel mixing devices, we study passive mixing effects by performing two-dimensional numerical simulations of wavy wall in channel flow using IBM. Methods: The continuity and Navier-Stokes equations governing the flow are solved by fractional step based finite volume method on a staggered Cartesian grid system. Fluid variables are described by Eulerian coordinates and solid boundary by Lagrangian coordinates. A four-point Dirac delta function is used to couple both the coordinate variables. A momentum forcing term is added to the governing equation in order to impose the no-slip boundary condition between the wavy wall and fluid interface. Results: Parametric study is carried out to analyze the fluid flow characteristics by varying amplitude and wavelength of wavy wall configurations for different Reynolds number. Conclusion: Configurations of wavy wall microchannels having a higher amplitude and lower wavelengths show optimum results for mixing applications. © 2020 Bentham Science Publishers.Item Numerical simulation of buckling and asymmetric behavior of flexible filament using temporal second-order immersed boundary method(Emerald Publishing, 2020) Kanchan, M.; Maniyeri, R.Purpose: The purpose of this paper is to perform two-dimensional numerical simulation involving fluid-structure interaction of flexible filament. The filament is tethered to the bottom of a rectangular channel with oscillating fluid flow inlet conditions at low Reynolds number. The simulations are performed using a temporal second-order finite volume-based immersed boundary method (IBM). Further, to understand the relation between different aspect ratios i.e. ratio of filament length to channel height (Len/H) and fixed channel geometry ratio, i.e. ratio of channel height to channel length (H/Lc) on mixing and pumping capabilities. Design/methodology/approach: The discretization of governing continuity and Navier–Stokes equation is done by finite-volume method on a staggered Cartesian grid. SIMPLE algorithm is used to solve fluid velocity and pressure terms. Two cases of oscillatory flow conditions are used with the flexible filament tethered at the center of bottom channel wall. The first case is sinusoidal oscillatory flow with phase shift (SOFPS) and second case is sinusoidal oscillatory flow without phase shift (SOF). The simulation results are validated with filament dynamics studies of previous researchers. Further, parametric analysis is carried to study the effect of filament length (aspect ratio), filament bending rigidity and Reynolds number on the complex deformation and behavior of flexible filament interacting with nearby oscillating fluid motion. Findings: It is found that selection of right filament length and bending rigidity is crucial for fluid mixing scenarios. The phase shift in fluid motion is also found to critically effect filament displacement dynamics, especially for rigid filaments. Aspect ratio, suitable for mixing applications is dependent on channel geometry ratio. Symmetric deformation is observed for filaments subjected to SOFPS condition irrespective of bending rigidity, whereas medium and low rigidity filaments placed in SOF condition show severe asymmetric behavior. Two key findings of this study are: symmetric filament conformity without appreciable bending produces sweeping motion in fluid flow, which is highly suited for mixing application; and asymmetric behavior shown by the filament depicts antiplectic metachronism commonly found in beating cilia. As a result, it is possible to pin point the type of fluid motion governing fluid mixing and fluid pumping. The developed computational model can, thus, successfully demonstrate filament-fluid interaction for a wide variety of similar problems. Originality/value: The present study uses a temporal second-order finite volume-based IBM to examine flexible filament dynamics for various applications such as fluid mixing. Also, it highlights the relationship between channel geometry ratio and filament aspect ratio and its effect on filament sweep patterns. The study further reports the effect of filament displacement dynamics with or without phase shift for inlet oscillating fluid flow condition. © 2019, Emerald Publishing Limited.Item Experimental and numerical study of laminar separation bubble formation on low Reynolds number airfoil with leading-edge tubercles(Springer, 2020) Sreejith, B.K.; Sathyabhama, A.The present work reports the effect of leading-edge tubercles on aerodynamic performance and flow features of a cambered airfoil E216 at a Reynolds number of 100,000 and at various angles of attack in the pre-stall regime. Amplitude values of 2 mm, 4 mm and 8 mm and wavelength values of 15.5 mm, 31 mm and 62 mm are used for both experimental and simulation studies. The Transition-SST RANS model is used to simulate transition phenomenon (laminar separation bubble) and three-dimensional flow features over the airfoil. Wind tunnel experimental results are used for the performance analysis and the validation of the simulation methodology. The experimental values of Cl and Cd are 1.37 and 0.081, respectively, at a stall angle of 12 ? for the plain airfoil. The experimental results show that the lift generated by tubercled airfoils is higher than that produced by the plain airfoil in the pre-stall region but lower at the stall angle. A maximum benefit of 4.51% in Cl is obtained for the tubercled airfoil with the highest amplitude (8 mm) and wavelength (64 mm) at 6 ? angle of attack. A higher Cd is observed for all the tubercled airfoils than for the plain one. The simulation is mainly carried out to study the flow structure. Simulation results indicate the presence of laminar separation bubbles on the plain airfoil with a straight separation and reattachment line parallel to the trailing edge. The tubercles considerably altered the laminar separation bubble formation and the flow structure. A sinusoidal laminar separation bubble is formed on the tubercled airfoils with reduced bubble length. The laminar separation bubble along the trough is formed ahead of that at peak. Two pairs of counter-rotating vortices are formed on the airfoil surface along the trough at two different chord-wise locations which strongly alter the flow pattern over it. Prandtl’s secondary flow of the first kind is the key reason for the vortex formation. © 2020, The Brazilian Society of Mechanical Sciences and Engineering.