Abstract: Exploratory missions to Mars and Titan have increased recently with various endeavors to find an alternate home to humankind. The use of surface rovers has its limitations due to rugged and uneven surfaces of these planetary bodies. The use of aerial robots requires the complete aerodynamic characterization of these vehicles in the atmospheric conditions of these planetary bodies. The dynamic stall phenomenon is extremely important for rotary wings performance under low Reynolds number that can be encountered in Martian and Titan’s atmosphere. The current research focuses on the aerodynamic characterization and exploration of the dynamic stall phenomenon of two different airfoils viz. E387 and Selig-Donovan7003 in Martian and Titan’s atmosphere at low Reynolds numbers of 10000 and 50000. The two-dimensional numerical simulations are conducted using commercially available finite volume solver with multi-species non-reacting mixture of gases as the working fluid. The k-epsilon (k-ε) turbulence model is used to capture the unsteady flow separation and the effect of turbulence. The dynamic characteristics are studied at a fixed different constant rotational extreme of angles of attack. This study of airfoils at different low Reynolds number and atmospheric conditions on Mars and Titan will be resulting in defining the aerodynamic characteristics of these airfoils for unmanned aerial missions for outer space exploration.
Abstract: Biopolymers are popular in many areas, like petrochemicals, food industry and agriculture due to their favorable properties like environment-friendly, availability, and cost. In this study, a biopolymer gum Arabic was used to find its effect on the pressure drop at various concentrations (100 ppm – 300 ppm) with various Reynolds numbers (10000 – 45000). A rheological study was also done by using the same concentrations to find the effect of the shear rate on the shear viscosity. Experiments were performed to find the effect of injection of gum Arabic directly near the boundary layer and to investigate its effect on the maximum possible drag reduction. Experiments were performed on a test section having i.d of 19.50 mm and length of 3045 mm. The polymer solution was injected from the top of the test section by using a peristaltic pump. The concentration of the polymer solution and the Reynolds number were used as parameters to get maximum possible drag reduction. Water was circulated through a centrifugal pump having a maximum 3000 rpm and the flow rate was measured by using rotameter. Results were validated by using Virk's maximum drag reduction asymptote. A maximum drag reduction of 62.15% was observed with the maximum concentration of gum Arabic, 300 ppm. The solution was circulated in the closed loop to find the effect of degradation of polymers with a number of cycles on the drag reduction percentage. It was observed that the injection of the polymer solution in the boundary layer was showing better results than premixed solutions.
Abstract: An essential component of a finite volume method (FVM) is the advection scheme that estimates values on the cell faces based on the calculated values on the nodes or cell centers. The most widely used advection schemes are upwind schemes. These schemes have been developed in FVM on different kinds of structured and unstructured grids. In this research, the physical influence scheme (PIS) is developed for a cell-centered FVM that uses an implicit coupled solver. Results are compared with the exponential differencing scheme (EDS) and the skew upwind differencing scheme (SUDS). Accuracy of these schemes is evaluated for a lid-driven cavity flow at Re = 1000, 3200, and 5000 and a backward-facing step flow at Re = 800. Simulations show considerable differences between the results of EDS scheme with benchmarks, especially for the lid-driven cavity flow at high Reynolds numbers. These differences occur due to false diffusion. Comparing SUDS and PIS schemes shows relatively close results for the backward-facing step flow and different results in lid-driven cavity flow. The poor results of SUDS in the lid-driven cavity flow can be related to its lack of sensitivity to the pressure difference between cell face and upwind points, which is critical for the prediction of such vortex dominant flows.
Abstract: Sudden change of bed condition is frequent in open channel flow. Change of bed condition affects the turbulence characteristics in both streamwise and wall-normal direction. Understanding the turbulence intensity in open channel flow is of vital importance to the modeling of sediment transport and resuspension, bed formation, entrainment, and the exchange of energy and momentum. A comprehensive study was carried out to understand the extent of the effect of Reynolds number and bed roughness on different turbulence characteristics in an open channel flow. Four different bed conditions (impervious smooth bed, impervious continuous rough bed, pervious rough sand bed, and impervious distributed roughness) and two different Reynolds numbers were adopted for this cause. The effect of bed roughness on different turbulence characteristics is seen to be prevalent for most of the flow depth. Effect of Reynolds number on different turbulence characteristics is also evident for flow over different bed, but the extent varies on bed condition. Although the same sand grain is used to create the different rough bed conditions, the difference in turbulence characteristics is an indication that specific geometry of the roughness has an influence on turbulence characteristics. Roughness increases the contribution of the extreme turbulent events which produces very large instantaneous Reynolds shear stress and can potentially influence the sediment transport, resuspension of pollutant from bed and alter the nutrient composition, which eventually affect the sustainability of benthic organisms.
Abstract: Since graphene nanoplatelet (GNP) is a promising material due to desirable thermal properties, this paper is related to the thermophysical and heat transfer performance of covalently functionalized GNP-based water/ethylene glycol nanofluid through an annular channel. After experimentally measuring thermophysical properties of prepared samples, a computational fluid dynamics study has been carried out to examine the heat transfer and pressure drop of well-dispersed and stabilized nanofluids. The effect of concentration of GNP and Reynolds number at constant wall temperature boundary condition under turbulent flow regime on convective heat transfer coefficient has been investigated. Based on the results, for different Reynolds numbers, the convective heat transfer coefficient of the prepared nanofluid is higher than that of the base fluid. Also, the enhancement of convective heat transfer coefficient and thermal conductivity increase with the increase of GNP concentration in base-fluid. Based on the results of this investigation, there is a significant enhancement on the heat transfer rate associated with loading well-dispersed GNP in base-fluid.
Abstract: A two-dimensional numerical study for flow past a square cylinder in presence of flat plate both at upstream and downstream position is carried out using the single-relaxation-time lattice Boltzmann method for gap spacing 0.5 and 1. We select Reynolds numbers from 80 to 200. The wake structure mechanism within gap spacing and near wake region, vortex structures around and behind the main square cylinder in presence of flat plate are studied and compared with flow pattern around a single square cylinder. The results are obtained in form of vorticity contour, streamlines, power spectra analysis, time trace analysis of drag and lift coefficients. Four different types of flow patterns were observed in both configurations, named as (i) Quasi steady flow (QSF), (ii) steady flow (SF), (iii) shear layer reattachment (SLR), (iv) single bluff body (SBB). It is observed that upstream flat plate plays a vital role in significant drag reduction. On the other hand, rate of suppression of vortex shedding is high for downstream flat plate case at low Reynolds numbers. The reduction in mean drag force and root mean square value of drag force for upstream flat plate case are89.1% and 86.3% at (Re, g) = (80, 0.5d) and (120, 1d) and reduction for downstream flat plate case for mean drag force and root mean square value of drag force are 11.10% and 97.6% obtained at (180, 1d) and (180, 0.5d).
Abstract: Optimizations of Plate Heat Exchangers (PHS) have received great attention in the past decade. In this study, heat transfer and pressure drop coefficients are compared for rectangular and circular PHS employing numerical simulations. Plates are designed to have equivalent areas. Simulations were implemented to investigate the efficiency of PHSs considering heat transfer, friction factor and pressure drop. Amount of heat transfer and pressure drop was obtained for different range of Reynolds numbers. These two parameters were compared with aim of F "weighting factor correlation". In this comparison, the minimum amount of F indicates higher efficiency. Results reveal that the F value for rectangular shape is less than circular plate, and hence using rectangular shape of PHS is more efficient than circular one. It was observed that, the amount of friction factor is correlated to the Reynolds numbers, such that friction factor decreased in both rectangular and circular plates with an increase in Reynolds number. Furthermore, such simulations revealed that the amount of heat transfer in rectangular plate is more than circular plate for different range of Reynolds numbers. The difference is more distinct for higher Reynolds number. However, amount of pressure drop in circular plate is less than rectangular plate for the same range of Reynolds numbers which is considered as a negative point for rectangular plate efficiency. It can be concluded that, while rectangular PHSs occupy more space than circular plate, the efficiency of rectangular plate is higher.
Abstract: This paper investigates the dynamic behavior of laminar water flows inside superhydrophobic micro-pipes patterned with square micro-posts features under different operating conditions. It also investigates the effects of air fraction and Reynolds number on the frictional performance of these pipes. Rather than modeling the air-water interfaces of superhydrophobic as a flat inflexible surface, a transient, incompressible, three-dimensional, volume-of-fluid (VOF) methodology has been employed to continuously track the air–water interface shape inside micro-pipes. Also, the entrance effects on the flow field have been taken into consideration. The results revealed the strong dependency of the frictional performance on the air fractions and Reynolds number. The frictional resistance reduction becomes increasingly more significant at large air fractions and low Reynolds numbers. Increasing Reynolds number has an adverse effect on the frictional resistance reduction.
Abstract: The effectiveness of microchannels in enhancing heat
transfer has been demonstrated in the semiconductor industry. In
order to tap the microscale heat transfer effects into macro
geometries, overcoming the cost and technological constraints,
microscale passages were created in macro geometries machined
using conventional fabrication methods. A cylindrical insert was
placed within a pipe, and geometrical profiles were created on the
outer surface of the insert to enhance heat transfer under steady-state
single-phase liquid flow conditions. However, while heat transfer
coefficient values of above 10 kW/m2·K were achieved, the heat
transfer enhancement was accompanied by undesirable pressure drop
increment. Therefore, this study aims to address the high pressure
drop issue using Constructal theory, a universal design law for both
animate and inanimate systems. Two designs based on Constructal theory were developed to study
the effectiveness of Constructal features in reducing the pressure drop
increment as compared to parallel channels, which are commonly
found in microchannel fabrication. The hydrodynamic and heat
transfer performance for the Tree insert and Constructal fin (Cfin)
insert were studied using experimental methods, and the underlying
mechanisms were substantiated by numerical results. In technical
terms, the objective is to achieve at least comparable increment in
both heat transfer coefficient and pressure drop, if not higher
increment in the former parameter. Results show that the Tree insert improved the heat transfer
performance by more than 16 percent at low flow rates, as compared
to the Tree-parallel insert. However, the heat transfer enhancement
reduced to less than 5 percent at high Reynolds numbers. On the
other hand, the pressure drop increment stayed almost constant at 20
percent. This suggests that the Tree insert has better heat transfer
performance in the low Reynolds number region. More importantly,
the Cfin insert displayed improved heat transfer performance along
with favourable hydrodynamic performance, as compared to Cfinparallel
insert, at all flow rates in this study. At 2 L/min, the
enhancement of heat transfer was more than 30 percent, with 20
percent pressure drop increment, as compared to Cfin-parallel insert.
Furthermore, comparable increment in both heat transfer coefficient
and pressure drop was observed at 8 L/min. In other words, the Cfin
insert successfully achieved the objective of this study. Analysis of the results suggests that bifurcation of flows is
effective in reducing the increment in pressure drop relative to heat
transfer enhancement. Optimising the geometries of the Constructal
fins is therefore the potential future study in achieving a bigger stride
in energy efficiency at much lower costs.
Abstract: In this paper, the 2-D unsteady viscous flow around
two cam shaped cylinders in tandem arrangement is numerically
simulated in order to study the characteristics of the flow in turbulent
regimes. The investigation covers the effects of high subcritical and
supercritical Reynolds numbers and L/D ratio on total drag
coefficient. The equivalent diameter of cylinders is 27.6 mm The
space between center to center of two cam shaped cylinders is define
as longitudinal pitch ratio and it varies in range of 1.5< L/D
Abstract: An experimental study with four different types of bed
conditions was carried out to understand the effect of roughness in
open channel flow at two different Reynolds numbers. The bed
conditions include a smooth surface and three different roughness
conditions, which were generated using sand grains with a median
diameter of 2.46 mm. The three rough conditions include a surface
with distributed roughness, a surface with continuously distributed
roughness and a sand bed with a permeable interface. A commercial
two-component fibre-optic LDA system was used to conduct the
velocity measurements. The variables of interest include the mean
velocity, turbulence intensity, correlation between the streamwise and
the wall normal turbulence, Reynolds shear stress and velocity triple
products. Quadrant decomposition was used to extract the magnitude
of the Reynolds shear stress of the turbulent bursting events. The
effect of roughness was evident throughout the flow depth. The
results show that distributed roughness has the greatest roughness
effect followed by the sand bed and the continuous roughness.
Compared to the smooth bed, the streamwise turbulence intensity
reduces but the vertical turbulence intensity increases at a location
very close to the bed due to the introduction of roughness. Although
the same sand grain is used to create the three different rough bed
conditions, the difference in the turbulence intensity is an indication
that the specific geometry of the roughness has an influence on
turbulence structure.
Abstract: The hydrodynamic and thermal behaviors of fluid
flow in wavy microchannel are investigated numerically. Effects of
Reynolds number on the hydrodynamics and thermal behaviors are
investigated. Three cases of Reynolds number (580, 1244, and 1910)
are adopted in this study. It is found that the separation zone begin
appears when Reynolds number is greater than 1910 at the endsection
of the wave. Also it is found that dimensionless maximum
velocity at the mid-section of the wave decreases and becomes as a
turbulent behavior as Reynolds numbers increases. The maximum
temperature at the center line at the mid-section of the wave increases
as Reynolds number increases until it reaches the turbulent behavior
when Reynolds number is equal or greater than 1244, while this
behavior will be achieved at very high velocities at the end section of
the wave.
Abstract: In this numerical work, mixed convection and entropy
generation of Cu–water nanofluid in a lid-driven square cavity have
been investigated numerically using the Lattice Boltzmann Method.
Horizontal walls of the cavity are adiabatic and vertical walls have
constant temperature but different values. The top wall has been
considered as moving from left to right at a constant speed, U0. The
effects of different parameters such as nanoparticle volume
concentration (0–0.05), Rayleigh number (104–106) and Reynolds
numbers (1, 10 and 100) on the entropy generation, flow and
temperature fields are studied. The results have shown that addition
of nanoparticles to the base fluid affects the entropy generation, flow
pattern and thermal behavior especially at higher Rayleigh and low
Reynolds numbers. For pure fluid as well as nanofluid, the increase
of Reynolds number increases the average Nusselt number and the
total entropy generation, linearly. The maximum entropy generation
occurs in nanofluid at low Rayleigh number and at high Reynolds
number. The minimum entropy generation occurs in pure fluid at low
Rayleigh and Reynolds numbers. Also at higher Reynolds number,
the effect of Cu nanoparticles on enhancement of heat transfer was
decreased because the effect of lid-driven cavity was increased. The
present results are validated by favorable comparisons with
previously published results. The results of the problem are presented
in graphical and tabular forms and discussed.
Abstract: The reliability of the filtered HVBK model is now
investigated via some large eddy simulations (LES) of freely
decaying isotropic superfluid turbulence. For homogeneous
turbulence at very high Reynolds numbers, comparison of the terms
in the spectral kinetic energy budget equation indicates, in the
energy-containing range, that the production and energy transfer
effects become significant except for dissipation. In the inertial range,
where the two fluids are perfectly locked, the mutual friction maybe
neglected with respect to other terms. Also, the LES results for the
other terms of the energy balance are presented.
Abstract: A finite difference/front tracking method is used to
study the motion of three-dimensional deformable drops suspended in
plane Poiseuille flow at non-zero Reynolds numbers. A parallel
version of the code was used to study the behavior of suspension on a
reasonable grid resolution (grids). The viscosity and density of drops
are assumed to be equal to that of the suspending medium. The effect
of the Reynolds number is studied in detail. It is found that drops
with small deformation behave like rigid particles and migrate to an
equilibrium position about half way between the wall and the
centerline (the Segre-Silberberg effect). However, for highly
deformable drops there is a tendency for drops to migrate to the
middle of the channel, and the maximum concentration occurs at the
centerline. The effective viscosity of suspension and the fluctuation
energy of the flow across the channel increases with the Reynolds
number of the flow.
Abstract: In this study the augmentation of heat transfer in a rectangular channel with triangular vortex generators is evaluated. The span wise averaged Nusselt number, mean temperature and total heat flux are compared with and without vortex generators in the channel at a blade angle of 30° for Reynolds numbers 800, 1200, 1600, and 2000. The use of vortex generators increases the span wise averaged Nusselt number compared to the case without vortex generators considerably. At a particular blade angle, increasing the Reynolds number results in an enhancement in the overall performance and span wise averaged Nusselt number was found to be greater at particular location for larger Reynolds number. The total heat flux from the bottom wall with vortex generators was found to be greater than that without vortex generators and the difference increases with increase in Reynolds number.
Abstract: This paper reports an experimental investigation of the energy spectrum of turbulent velocity fields at low Reynolds numbers in grid turbulence. Hot wire measurements are carried out in grid turbulence with subjected to a 1.36:1 contraction of the wind tunnel. Three different grids are used: (i) large square perforated grid (mesh size 43.75mm), (ii) small square perforated grid (mesh size 14. and (iii) woven mesh grid (mesh size 5mm). The results indicate that the energy spectrum at small Reynolds numbers does not follow Kolmogorov’s universal scaling. It is further found that the critical Reynolds number, below which the scaling breaks down, is around 25.
Abstract: A method for simulating flow around the solid bodies has been presented using hybrid meshfree and mesh-based schemes. The presented scheme optimizes the computational efficiency by combining the advantages of both meshfree and mesh-based methods. In this approach, a cloud of meshfree nodes has been used in the domain around the solid body. These meshfree nodes have the ability to efficiently adapt to complex geometrical shapes. In the rest of the domain, conventional Cartesian grid has been used beyond the meshfree cloud. Complex geometrical shapes can therefore be dealt efficiently by using meshfree nodal cloud and computational efficiency is maintained through the use of conventional mesh-based scheme on Cartesian grid in the larger part of the domain. Spatial discretization of meshfree nodes has been achieved through local radial basis functions in finite difference mode (RBF-FD). Conventional finite difference scheme has been used in the Cartesian ‘meshed’ domain. Accuracy tests of the hybrid scheme have been conducted to establish the order of accuracy. Numerical tests have been performed by simulating two dimensional steady and unsteady incompressible flows around cylindrical object. Steady flow cases have been run at Reynolds numbers of 10, 20 and 40 and unsteady flow problems have been studied at Reynolds numbers of 100 and 200. Flow Parameters including lift, drag, vortex shedding, and vorticity contours are calculated. Numerical results have been found to be in good agreement with computational and experimental results available in the literature.
Abstract: An experimental investigation was conducted to study the effect of surface roughness on friction factor and heat transfer characteristics in single-phase fluid flow in a stainless steel micro-tube having diameter of 0.85 mm and average internal surface roughness of 1.7 μm with relative surface roughness of 0.002. Distilled water and R134a liquids were used as the working fluids and testing was conducted with Reynolds numbers ranging from 100 to 10,000 covering laminar, transition and turbulent flow conditions. The experiments were conducted with the micro-tube oriented horizontally with uniform heat fluxes applied at the test section. The results indicated that the friction factor of both water and R134a can be predicted by the Hagen-Poiseuille equation for laminar flow and the modified Miller correlation for turbulent flow and early transition from laminar to turbulent flows. The heat transfer results of water and R134a were in good agreement with the conventional theory in the laminar flow region and lower than the Adam’s correlation for turbulent flow region which deviates from conventional theory.
Abstract: The migration of a deformable drop in simple shear
flow at finite Reynolds numbers is investigated numerically by
solving the full Navier-Stokes equations using a finite
difference/front tracking method. The objectives of this study are to
examine the effectiveness of the present approach to predict the
migration of a drop in a shear flow and to investigate the behavior of
the drop migration with different drop sizes and non-unity viscosity
ratios. It is shown that the drop deformation depends strongly on the
capillary number, so that; the proper non-dimensional number for the
interfacial tension is the capillary number. The rate of migration
increased with increasing the drop radius. In other words, the
required time for drop migration to the centreline decreases. As the
viscosity ratio increases, the drop rotates more slowly and the
lubrication force becomes stronger. The increased lubrication force
makes it easier for the drop to migrate to the centre of the channel.
The migration velocity of the drop vanishes as the drop reaches the
centreline under viscosity ratio of one and non-unity viscosity ratios.
To validate the present calculations, some typical results are
compared with available experimental and theoretical data.