Journal Articles
Permanent URI for this collectionhttps://idr.nitk.ac.in/handle/123456789/19884
Browse
16 results
Search Results
Item Numerical simulation of laminar flow past a circular cylinder(2009) Rajani, B.N.; Kandasamy, A.; Majumdar, S.The present paper focuses on the analysis of two- and three-dimensional flow past a circular cylinder in different laminar flow regimes. In this simulation, an implicit pressure-based finite volume method is used for time-accurate computation of incompressible flow using second order accurate convective flux discretisation schemes. The computation results are validated against measurement data for mean surface pressure, skin friction coefficients, the size and strength of the recirculating wake for the steady flow regime and also for the Strouhal frequency of vortex shedding and the mean and RMS amplitude of the fluctuating aerodynamic coefficients for the unsteady periodic flow regime. The complex three dimensional flow structure of the cylinder wake is also reasonably captured by the present prediction procedure. © 2008 Elsevier Inc. All rights reserved.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 An investigation on the effect of pitchwise endwall design in a turbine cascade at different incidence angles(Elsevier Masson SAS 62 rue Camille Desmoulins Issy les Moulineaux Cedex 92442, 2017) Kiran, K.N.; Anish, S.This paper describes the effects of non-axisymmetric endwall profiling on the aerodynamic performance of a linear turbine cascade at different incidence angles. The sinusoidal profiling is carried out with constant profile curvature along the mean streamline path. Three different profiles, with varying hump to dip height, are analyzed numerically and the performances are compared with the planar profile. Reynolds Averaged Navier Stokes (RANS) equations are solved in their conservative form using Finite Volume Method with SST turbulence model. The calculated results indicate that the profiled endwall minimizes the lateral movement of weaker boundary layer fluid from the hub-pressure side corner. In comparison with planar case, the flow deviations are largely contained with endwall profiling but closer to the endwall it enhances the overturning and secondary flow kinetic energy. The reduction in loss coefficient is estimated to be 1.3%, 8.7% and 38% for incidence angles of ?10°, nominal and +15° respectively. The sinusoidal profiling has brought down the pitch averaged flow deviation and secondary flow kinetic energy at nominal and positive incidence angles but the impact is insignificant at negative incidence. Profiling minimizes the rolling up of the passage vortex and makes the passage vortex to migrate closer to the endwall. This flow modification brings down the losses in the core flow but enhances the losses near the endwall. © 2017 Elsevier Masson SASItem Inverse approach for estimating boundary properties in a transient fin problem(Springer, 2018) Gnanasekaran, N.; Balaji, S.A solution methodology is proposed for an inverse estimation of boundary conditions from the knowledge of transient temperature data. A forward model based on prevalent time-dependent heat conduction fin equation is solved using a fully implicit finite volume method. First, the inverse model is formulated and accomplished for time-invariant heat flux at the fin base, and later extended to transient heat flux, base temperature and average heat transfer coefficient. Secondly, the Nusselt number is then replaced with Rayleigh number in the forward model to realistically estimate the base temperature, which varies with respect to time, based on in-house transient fin heat transfer experiments. This scenario further corroborates the validation of the proposed inverse approach. The experimental set-up consists of a mild steel 250×150×6mm3 fin mounted centrally on an aluminium base 250×150×8mm3 plate. The base is attached to an electrical heater and insulated with glass-wool to prevent heat loss to surroundings. Five calibrated K-type thermocouples are used to measure temperature along the fin. The functional form of the unknown parameters is not known beforehand; sensitivity studies are performed to determine suitability of the estimation and location of sensors for the inverse approach. Conjugate gradient method with adjoint equation is chosen as the inverse technique and the study is performed as a numerical optimization; subsequently, the estimates show satisfactory results. © 2018, Indian Academy of Sciences.Item Numerical simulation of oscillating lid driven square cavity(Elsevier B.V., 2018) Indukuri, J.V.; Maniyeri, R.This paper aim to develop a two-dimensional computational model to study the fluid dynamic behaviour in a square cavity driven by an oscillating lid using staggered grid based finite volume method. Firstly the developed computational model is validated with that of other researcher's results for the case of finite wall motion. Later the numerical simulations are performed for the case of top wall oscillations for various combinations of Reynolds number and frequencies. From these simulations an optimum frequency is chosen and then with the optimum frequency the simulations are carried out to explore the vortex behaviour for the cases of parallel wall oscillations (both top and bottom walls moving in the same direction) and anti-parallel wall oscillations (both top and bottom walls moving in the opposite direction). From these simulations it may be concluded that Re = 1000 is medium range of operation for better mixing inside the cavity for the cases of parallel and anti-parallel wall oscillations. © 2017 Faculty of Engineering, Alexandria UniversityItem Numerical study of forced convection heat transfer in an oscillating lid driven cavity with heated top wall(International Information and Engineering Technology Association info@iieta.org, 2018) Indukuri, J.V.; Maniyeri, R.The present work is aimed to study the fluid flow and heat transfer behaviour in an oscillating lid-driven cavity using finite volume method by developing a two-dimensional computational model. Firstly, the developed computational model is validated by comparing our numerical results with that of the other researcher's results for the case of wall moving with finite motion. Next, the simulations are conducted for oscillating cavity problem with top wall oscillation for Reynolds number (Re =5 00) and frequency (?=2?/6). Later, the simulations are carried out for cases of oscillating parallel wall (upper and lower walls oscillating with sync) and oscillating anti-parallel wall (upper and lower walls oscillating with reverse sync) with the same optimum frequency and fixed Reynolds number (Re = 500). Secondly, the same optimum frequency is used to study the heat transfer characteristics in an oscillating lid-driven square cavity with heated top wall and lower cold wall for various Reynolds numbers (Re = 100-1000) and Prandtl numbers (Pr = 0.2 to 1.0). From this study, it is found that for high Prandtl number case (Pr = 1.0) the flow of high temperature isotherms inside the cavity is more when compared with low Prandtl number cases due to increase in molecular diffusion of momentum. © 2018 International Information and Engineering Technology Association.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 Numerical simulation and prediction model development of multiple flexible filaments in viscous shear flow using immersed boundary method and artificial neural network techniques(IOP Publishing Ltd custserv@iop.org, 2020) Kanchan, M.; Maniyeri, R.Many chemical and biological systems have applications involving fluid-structure interaction (FSI) of flexible filaments in viscous fluid. The dynamics of single- and multiple-filament interaction are of interest to engineers and biologists working in the area of DNA fragmentation, protein synthesis, polymer segmentation, folding-unfolding analysis of natural and synthetic fibers, etc. To perform numerical simulation of the above-mentioned FSI applications is challenging. In this direction, methods like the immersed boundary method (IBM) have been quite successful. We simulate the dynamics of multiple flexible filaments subjected to planar shear flow at low Reynolds number using the finite volume method-based IBM. The governing continuity and Navier-Stokes equations are solved by the SIMPLE algorithm on a staggered Cartesian grid system. The validation of the developed model is done using previous works. The length of the filament, its bending rigidity and fluid shear rate are taken as parametric variables and numerical simulations are carried out. Viscous flow forcing and fractional contraction terms are incorporated so as to effectively categorize filament motion into various deformation regimes. The effects of tumbling motion on the filament migration and recuperative aspects are studied. The mutual interaction of two filaments placed side by side is thus observed. Finally, an artificial neural network model is developed from the IBM simulation results to predict tumbling counts for different filament parameters. © 2020 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.