Abstract: Heavy rainfall greatly affects the aerodynamic performance of the aircraft. There are many accidents of aircraft caused by aerodynamic efficiency degradation by heavy rain. In this Paper we have studied the heavy rain effects on the aerodynamic efficiency of NACA 64-210 & NACA 0012 airfoils. For our analysis, CFD method and preprocessing grid generator are used as our main analytical tools, and the simulation of rain is accomplished via two phase flow approach-s Discrete Phase Model (DPM). Raindrops are assumed to be non-interacting, non-deforming, non-evaporating and non-spinning spheres. Both airfoil sections exhibited significant reduction in lift and increase in drag for a given lift condition in simulated rain. The most significant difference between these two airfoils was the sensitivity of the NACA 64-210 to liquid water content (LWC), while NACA 0012 performance losses in the rain environment is not a function of LWC . It is expected that the quantitative information gained in this paper will be useful to the operational airline industry and greater effort such as small scale and full scale flight tests should put in this direction to further improve aviation safety.
Abstract: The electrokinetic flow resistance (electroviscous
effect) is predicted for steady state, pressure-driven liquid flow at
low Reynolds number in a microfluidic contraction of rectangular
cross-section. Calculations of the three dimensional flow are
performed in parallel using a finite volume numerical method. The
channel walls are assumed to carry a uniform charge density and the
liquid is taken to be a symmetric 1:1 electrolyte. Predictions are
presented for a single set of flow and electrokinetic parameters. It is
shown that the magnitude of the streaming potential gradient and the
charge density of counter-ions in the liquid is greater than that in
corresponding two-dimensional slit-like contraction geometry. The
apparent viscosity is found to be very close to the value for a
rectangular channel of uniform cross-section at the chosen Reynolds
number (Re = 0.1). It is speculated that the apparent viscosity for the
contraction geometry will increase as the Reynolds number is
reduced.
Abstract: The present microfluidic study is emphasizing the flow behavior within a Y shape micro-bifurcation in two similar flow configurations. We report here a numerical and experimental investigation on the velocity profiles evolution and secondary flows, manifested at different Reynolds numbers (Re) and for two different boundary conditions. The experiments are performed using special designed setup based on optical microscopic devices. With this setup, direct visualizations and quantitative measurements of the path-lines are obtained. A Micro-PIV measurement system is used to obtain velocity profiles distributions in a spatial evolution in the main flows domains. The experimental data is compared with numerical simulations performed with commercial computational code FLUENT in a 3D geometry with the same dimensions as the experimental one. The numerical flow patterns are found to be in good agreement with the experimental manifestations.
Abstract: Turbulence studies were made in the wake of a rotating
circular cylinder in a uniform free stream. The interest was to
examine the turbulence properties at the suppression of periodicity in
vortex formation process. An experimental study of the turbulent near
wake of a rotating circular cylinder was made at a Reynolds number
of 9000 for velocity ratios, λ between 0 and 2.7. Hot-wire
anemometry and particle image velocimetry results indicate that the
rotation of the cylinder causes significant changes in the vortical
activities. The turbulence quantities are getting smaller as λ increases
due to suppression of coherent vortex structures.
Abstract: As it is known, buoyancy and drag forces rule bubble's rise velocity in a liquid column. These forces are strongly dependent on fluid properties, gravity as well as equivalent's diameter. This study reports a set of bubble rising velocity experiments in a liquid column using water or glycerol. Several records of terminal velocity were obtained. The results show that bubble's rise terminal velocity is strongly dependent on dynamic viscosity effect. The data set allowed to have some terminal velocities data interval of 8.0 ? 32.9 cm/s with Reynolds number interval 1.3 -7490. The bubble's movement was recorded with a video camera. The main goal is to present an original set data and results that will be discussed based on two-phase flow's theory. It will also discussed, the prediction of terminal velocity of a single bubble in liquid, as well as the range of its applicability. In conclusion, this study presents general expressions for the determination of the terminal velocity of isolated gas bubbles of a Reynolds number range, when the fluid proprieties are known.
Abstract: This paper presents a numerical study on surface heat
transfer characteristics of laminar air flows in parallel-plate dimpled
channels. The two-dimensional numerical model is provided by
commercial code FLUENT and the results are obtained for channels
with symmetrically opposing hemi-cylindrical cavities onto both
walls for Reynolds number ranging from 1000 to 2500. The influence
of variations in relative depth of dimples (the ratio of cavity depth to
the cavity curvature diameter), the number of them and the thermophysical
properties of channel walls on heat transfer enhancement is
studied. The results are evident for existence of an optimum value for
the relative depth of dimples in which the largest wall heat flux and
average Nusselt number can be achieved. In addition, the results of
conjugation simulation indicate that the overall influence of the ratio
of wall thermal conductivity to the one of the fluid on heat transfer
rate is not much significant and can be ignored.
Abstract: A multiple-option analytical model for the evaluation of the energy performance and distribution of aerodynamic forces acting on a vertical-axis Darrieus wind turbine depending on both rotor architecture and operating conditions is presented. For this purpose, a numerical algorithm, capable of generating the desired rotor conformation depending on design geometric parameters, is coupled to a Single/Double-Disk Multiple-Streamtube Blade Element – Momentum code. Both single and double-disk configurations are analyzed and model predictions are compared to literature experimental data in order to test the capability of the code for predicting rotor performance. Effective airfoil characteristics based on local blade Reynolds number are obtained through interpolation of literature low-Reynolds airfoil databases. Some corrections are introduced inside the original model with the aim of simulating also the effects of blade dynamic stall, rotor streamtube expansion and blade finite aspect ratio, for which a new empirical relationship to better fit the experimental data is proposed. In order to predict also open field rotor operation, a freestream wind shear profile is implemented, reproducing the effect of atmospheric boundary layer.
Abstract: The dispersion of heavy particles line in an isotropic
and incompressible three-dimensional turbulent flow has been
studied using the Kinematic Simulation techniques to find out the
evolution of the line fractal dimension. In this study, the fractal
dimension of the line is found for different cases of heavy particles
inertia (different Stokes numbers) in the absence of the particle
gravity with a comparison with the fractal dimension obtained in the
diffusion case of material line at the same Reynolds number. It can
be concluded for the dispersion of heavy particles line in turbulent
flow that the particle inertia affect the fractal dimension of a line
released in a turbulent flow for Stokes numbers 0.02 < St < 2. At the
beginning for small times, most of the different cases are not affected
by the inertia until a certain time, the particle response time τa, with
larger time as the particles inertia increases, the fractal dimension of
the line increases owing to the particles becoming more sensitive to
the small scales which cause the change in the line shape during its
journey.
Abstract: The present study investigates numerically the
phenomenon of vortex-shedding and its suppression in twodimensional
mixed convective flow past a square cylinder under the
joint influence of buoyancy and free-stream orientation with respect
to gravity. The numerical experiments have been conducted at a
fixed Reynolds number (Re) of 100 and Prandtl number (Pr) of 0.71,
while Richardson number (Ri) is varied from 0 to 1.6 and freestream
orientation, α, is kept in the range 0o≤ α ≤ 90o, with 0o
corresponding to an upward flow and 90o representing a cross-flow
scenario, respectively. The continuity, momentum and energy
equations, subject to Boussinesq approximation, are discretized using
a finite difference method and are solved by a semi-explicit pressure
correction scheme. The critical Richardson number, leading to the
suppression of the vortex-shedding (Ric), is estimated by using
Stuart-Landau theory at various free-stream orientations and the
neutral curve is obtained in the Ri-α plane. The neutral curve
exhibits an interesting non-monotonic behavior with Ric first
increasing with increasing values of α upto 45o and then decreasing
till 70o. Beyond 70o, the neutral curve again exhibits a sharp
increasing asymptotic trend with Ric approaching very large values
as α approaches 90o. The suppression of vortex shedding is not
observed at α = 90o (cross-flow). In the unsteady flow regime, the
Strouhal number (St) increases with the increase in Richardson
number.
Abstract: The pressure drag from a cam shaped tube in cross flows have been investigated experimentally using pressure distribution measurement. The range of angle of attack and Reynolds number based on an equivalent circular tube are within 0≤α≤360° and 2×104< Reeq < 3.4 ×104, respectively. It is found that the pressure drag coefficient is at its highest at α=90° and 270° over the whole range of Reynolds number. Results show that the pressure drag coefficient of the cam shaped tube is lower than that of circular tube with the same surface area for more of the angles of attack. Furthermore, effects of the diameter ratio and finite length of the cam shaped tube upon the pressure drag coefficient are discussed.
Abstract: Addition of milli or micro sized particles to the heat
transfer fluid is one of the many techniques employed for improving
heat transfer rate. Though this looks simple, this method has
practical problems such as high pressure loss, clogging and erosion
of the material of construction. These problems can be overcome by
using nanofluids, which is a dispersion of nanosized particles in a
base fluid. Nanoparticles increase the thermal conductivity of the
base fluid manifold which in turn increases the heat transfer rate.
Nanoparticles also increase the viscosity of the basefluid resulting in
higher pressure drop for the nanofluid compared to the base fluid. So
it is imperative that the Reynolds number (Re) and the volume
fraction have to be optimum for better thermal hydraulic
effectiveness. In this work, the heat transfer enhancement using
aluminium oxide nanofluid using low and high volume fraction
nanofluids in turbulent pipe flow with constant wall temperature has
been studied by computational fluid dynamic modeling of the
nanofluid flow adopting the single phase approach. Nanofluid, up till
a volume fraction of 1% is found to be an effective heat transfer
enhancement technique. The Nusselt number (Nu) and friction factor
predictions for the low volume fractions (i.e. 0.02%, 0.1 and 0.5%)
agree very well with the experimental values of Sundar and Sharma
(2010). While, predictions for the high volume fraction nanofluids
(i.e. 1%, 4% and 6%) are found to have reasonable agreement with
both experimental and numerical results available in the literature.
So the computationally inexpensive single phase approach can be
used for heat transfer and pressure drop prediction of new nanofluids.
Abstract: Experiments have been carried out at sub-critical
Reynolds number to investigate free-to-roll motions induced by
forebody and/or wings complex flow on a 30° swept back nonslender
wings-slender body-model for static and dynamic (pitch-up)
cases. For the dynamic (pitch-up) case it has been observed that roll
amplitude decreases and lag increases with increase in pitching
speed. Decrease in roll amplitude with increase in pitch rate is
attributed to low disturbing rolling moment due to weaker interaction
between forebody and wing flow components. Asymmetric forebody
vortices dominate and control the roll motion of the model in
dynamic case when non-dimensional pitch rate ≥ 1x10-2.
Effectiveness of the active control scheme utilizing rotating nose with
artificial tip perturbation is observed to be low in the angle of attack
region where the complex flow over the wings has contributions from
both forebody and wings.
Abstract: A new design of a planar passive T-micromixer with fin-shaped baffles in the mixing channel is presented. The mixing efficiency and the level of pressure loss in the channel have been investigated by numerical simulations in the range of Reynolds number (Re) 1 to 50. A Mixing index (Mi) has been defined to quantify the mixing efficiency, which results over 85% at both ends of the Re range, what demonstrates the micromixer can enhance mixing using the mechanisms of diffusion (lower Re) and convection (higher Re). Three geometric dimensions: radius of baffle, baffles pitch and height of the channel define the design parameters, and the mixing index and pressure loss are the performance parameters used to optimize the micromixer geometry with a multi-criteria optimization method. The Pareto front of designs with the optimum trade-offs, maximum mixing index with minimum pressure loss, is obtained. Experiments for qualitative and quantitative validation have been implemented.
Abstract: The operation performance of a valveless micro-pump
is strongly dependent on the shape of connected nozzle/diffuser and
Reynolds number. The aims of present work are to compare the
performance curves of micropump with the original straight
nozzle/diffuser and contoured nozzle/diffuser under different back
pressure conditions. The tested valveless micropumps are assembled
of five pieces of patterned PMMA plates with hot-embracing
technique. The structures of central chamber, the inlet/outlet
reservoirs and the connected nozzle/diffuser are fabricated with laser
cutting machine. The micropump is actuated with circular-type PZT
film embraced on the bottom of central chamber. The deformation of
PZT membrane with various input voltages is measured with a
displacement laser probe. A simple testing facility is also constructed
to evaluate the performance curves for comparison.
In order to observe the evaluation of low Reynolds number
multiple vortex flow patterns within the micropump during suction
and pumping modes, the unsteady, incompressible laminar
three-dimensional Reynolds-averaged Navier-Stokes equations are
solved. The working fluid is DI water with constant thermo-physical
properties. The oscillating behavior of PZT film is modeled with the
moving boundary wall in way of UDF program. With the dynamic
mesh method, the instants pressure and velocity fields are obtained
and discussed.Results indicated that the volume flow rate is not
monotony increased with the oscillating frequency of PZT film,
regardless of the shapes of nozzle/diffuser. The present micropump
can generate the maximum volume flow rate of 13.53 ml/min when
the operation frequency is 64Hz and the input voltage is 140 volts.
The micropump with contoured nozzle/diffuser can provide 7ml/min
flow rate even when the back pressure is up to 400 mm-H2O. CFD
results revealed that the flow central chamber was occupied with
multiple pairs of counter-rotating vortices during suction and
pumping modes. The net volume flow rate over a complete
oscillating periodic of PZT
Abstract: The three-dimensional incompressible flow past a
rectangular open cavity is investigated, where the aspect ratio of the
cavity is considered as 4. The principle objective is to use large-eddy
simulation to resolve and control the large-scale structures, which are
largely responsible for flow oscillations in a cavity. The flow past an
open cavity is very common in aerospace applications and can be a
cause of acoustic source due to hydrodynamic instability of the shear
layer and its interactions with the downstream edge. The unsteady
Navier-stokes equations have been solved on a staggered mesh using
a symmetry-preserving central difference scheme. Synthetic jet has
been used as an active control to suppress the cavity oscillations in
wake mode for a Reynolds number of ReD = 3360. The effect of
synthetic jet has been studied by varying the jet amplitude and
frequency, which is placed at the upstream wall of the cavity. The
study indicates that there exits a frequency band, which is larger than
a critical value, is effective in attenuating cavity oscillations when
blowing ratio is more than 1.0.
Abstract: An original Direct Numerical Simulation (DNS) method to tackle the problem of particulate flows at moderate to high concentration and finite Reynolds number is presented. Our method is built on the framework established by Glowinski and his coworkers [1] in the sense that we use their Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) formulation and their operator-splitting idea but differs in the treatment of particle collisions. The novelty of our contribution relies on replacing the simple artificial repulsive force based collision model usually employed in the literature by an efficient Discrete Element Method (DEM) granular solver. The use of our DEM solver enables us to consider particles of arbitrary shape (at least convex) and to account for actual contacts, in the sense that particles actually touch each other, in contrast with the simple repulsive force based collision model. We recently upgraded our serial code, GRIFF 1 [2], to full MPI capabilities. Our new code, PeliGRIFF 2, is developed under the framework of the full MPI open source platform PELICANS [3]. The new MPI capabilities of PeliGRIFF open new perspectives in the study of particulate flows and significantly increase the number of particles that can be considered in a full DNS approach: O(100000) in 2D and O(10000) in 3D. Results on the 2D/3D sedimentation/fluidization of isometric polygonal/polyedral particles with collisions are presented.
Abstract: The measurement of aerodynamic forces and moments
acting on an aircraft model is important for the development of wind
tunnel measurement technology to predict the performance of the full
scale vehicle. The potentials of an aircraft model with and without
winglet and aerodynamic characteristics with NACA wing No. 65-3-
218 have been studied using subsonic wind tunnel of 1 m × 1 m
rectangular test section and 2.5 m long of Aerodynamics Laboratory
Faculty of Engineering (University Putra Malaysia). Focusing on
analyzing the aerodynamic characteristics of the aircraft model, two
main issues are studied in this paper. First, a six component wind
tunnel external balance is used for measuring lift, drag and pitching
moment. Secondly, Tests are conducted on the aircraft model with
and without winglet of two configurations at Reynolds numbers
1.7×105, 2.1×105, and 2.5×105 for different angle of attacks. Fuzzy
logic approach is found as efficient for the representation,
manipulation and utilization of aerodynamic characteristics.
Therefore, the primary purpose of this work was to investigate the
relationship between lift and drag coefficients, with free-stream
velocities and angle of attacks, and to illustrate how fuzzy logic
might play an important role in study of lift aerodynamic
characteristics of an aircraft model with the addition of certain
winglet configurations. Results of the developed fuzzy logic were
compared with the experimental results. For lift coefficient analysis,
the mean of actual and predicted values were 0.62 and 0.60
respectively. The coreelation between actual and predicted values
(from FLS model) of lift coefficient in different angle of attack was
found as 0.99. The mean relative error of actual and predicted valus
was found as 5.18% for the velocity of 26.36 m/s which was found to
be less than the acceptable limits (10%). The goodness of fit of
prediction value was 0.95 which was close to 1.0.
Abstract: In this work, we perform numerical simulation of fluid
mixing in a floor-grooved micro-channel with wavy sidewalls which
may impose perturbation on the helical flow induced by the slanted
grooves on the channel floor. The perturbation is caused by separation
vortices in the recesses of the wavy-walled channel as the Reynolds
number is large enough. The results show that the effects of the wavy
sidewalls of the present micromixer on the enhancement of fluid
mixing increase with the increase of Reynolds number. The degree of
mixing increases with the increase of the corrugation angle, until the
angle is greater than 45 degrees. Besides, the pumping pressure of the
micromixer increases with the increase of the corrugation angle
monotonically. Therefore, we would suggest setting the corrugation
angle of the wavy sidewalls to be 45 degrees.
Abstract: The effects of dynamic subgrid scale (SGS) models are
investigated in variational multiscale (VMS) LES simulations of bluff
body flows. The spatial discretization is based on a mixed finite
element/finite volume formulation on unstructured grids. In the VMS
approach used in this work, the separation between the largest and the
smallest resolved scales is obtained through a variational projection
operator and a finite volume cell agglomeration. The dynamic version
of Smagorinsky and WALE SGS models are used to account for
the effects of the unresolved scales. In the VMS approach, these
effects are only modeled in the smallest resolved scales. The dynamic
VMS-LES approach is applied to the simulation of the flow around a
circular cylinder at Reynolds numbers 3900 and 20000 and to the flow
around a square cylinder at Reynolds numbers 22000 and 175000. It
is observed as in previous studies that the dynamic SGS procedure
has a smaller impact on the results within the VMS approach than in
LES. But improvements are demonstrated for important feature like
recirculating part of the flow. The global prediction is improved for
a small computational extra cost.
Abstract: Experiments have been performed to investigate the radiation effects on mixed convection heat transfer for thermally developing airflow in vertical ducts with two differentially heated isothermal walls and two adiabatic walls. The investigation covers the Reynolds number Re = 800 to Re = 2900, heat flux varied from 256 W/m2 to 863 W/m2, hot wall temperature ranges from 27°C to 100 °C, aspect ratios 1 & 0.5 and the emissivity of internal walls are 0.05 and 0.85. In the present study, combined flow visualization was conducted to observe the flow patterns. The effect of surface temperature along the walls was studied to investigate the local Nusselt number variation within the duct. The result shows that flow condition and radiation significantly affect the total Nusselt number and tends to reduce the buoyancy condition.