Abstract: Cavitation inside diesel injector nozzle is investigated using Reynolds-Stress-Navier stokes equations. Schnerr-Sauer cavitation model is used for modeling cavitation inside diesel injector nozzle. The carrying fluid utilized in the current study is diesel fuel. The flow is verified at the beginning by comparing with the previous experimental data and it was found that K-Epsilon turbulent model could lead to a better accuracy comparing to K-Omega turbulent model. Moreover, mass flow rate obtained numerically is compared with the experimental value and discrepancy was found to be less than 5% - which shows the accuracy of the current results. Finally, a real-size four-hole nozzle is investigated and the flow inside it is visualized based on velocity profile, discharge coefficient and cavitation number. It was found that the mesh density could be reduced significantly by utilizing periodic boundary condition. Velocity contour at the mid nozzle showed that maximum value of velocity occurs at the end of the needle before entering the orifice area. Last but not least, at the same boundary conditions, when different needle heights were utilized, it was found that as needle height increases with an increase in cavitation number, discharge coefficient increases, while the mentioned increases is more tangible at smaller values of needle heights.
Abstract: In this work, small-scale experiments of fire whirl were conducted to study the spinning fire phenomenon and to gain comprehensive understandings of fire tornadoes and the factors that affect their behavior. High speed imaging was used to track the flames at both temporal and spatial scales. This allowed us to better understand the role of the near-ground height in creating a boundary layer flow profile that, in turn contributes to formation of vortices around the fire, and consequent fire whirls. Based on the results obtained from these observations, we were able to spot the differences in the fuel burning rate of the fire itself as a function of a newly defined specific non-dimensional near-ground height. Based on our observations, there is a cutoff non-dimensional height, beyond which a normal fire can be turned into a fire whirl. Additionally, the results showed that the fire burning rate decreases by moving the fire to a height higher than the ground level. These effects were justified by the interactions between vortices formed by, the back pressure and the boundary layer velocity profile, and the vortices generated by the fire itself.
Abstract: In this paper, an enhancement of the heat transfer using non-Newtonian nanofluids by magnetohydrodynamic (MHD) mixed convection along stretching sheets embedded in an isotropic porous medium is investigated. Case of the Maxwell nanofluids is studied using the two phase mathematical model of nanofluids and the Darcy model is applied for the porous medium. Important effects are taken into account, namely, non-linear thermal radiation, convective boundary conditions, electromagnetic force and presence of the heat source/sink. Suitable similarity transformations are used to convert the governing equations to a system of ordinary differential equations then it is solved numerically using a fourth order Runge-Kutta method with shooting technique. The main results of the study revealed that the velocity profiles are decreasing functions of the Darcy number, the Deborah number and the magnetic field parameter. Also, the increase in the non-linear radiation parameters causes an enhancement in the local Nusselt number.
Abstract: The Oscillatory electroosmotic flow (OEOF) in power
law fluids through a microchannel is studied numerically. A
time-dependent external electric field (AC) is suddenly imposed
at the ends of the microchannel which induces the fluid motion.
The continuity and momentum equations in the x and y direction
for the flow field were simplified in the limit of the lubrication
approximation theory (LAT), and then solved using a numerical
scheme. The solution of the electric potential is based on the
Debye-H¨uckel approximation which suggest that the surface potential
is small,say, smaller than 0.025V and for a symmetric (z : z)
electrolyte. Our results suggest that the velocity profiles across
the channel-width are controlled by the following dimensionless
parameters: the angular Reynolds number, Reω, the electrokinetic
parameter, ¯κ, defined as the ratio of the characteristic length scale
to the Debye length, the parameter λ which represents the ratio
of the Helmholtz-Smoluchowski velocity to the characteristic length
scale and the flow behavior index, n. Also, the results reveal that
the velocity profiles become more and more non-uniform across the
channel-width as the Reω and ¯κ are increased, so oscillatory OEOF
can be really useful in micro-fluidic devices such as micro-mixers.
Abstract: The earliest theories of sloshing waves and solitary waves based on potential theory idealisations and irrotational flow have been extended to be applicable to more realistic domains. To this end, the computational fluid dynamics (CFD) methods are widely used. Three-dimensional CFD methods such as Navier-Stokes solvers with volume of fluid treatment of the free surface and Navier-Stokes solvers with mappings of the free surface inherently impose high computational expense; therefore, considerable effort has gone into developing depth-averaged approaches. Examples of such approaches include Green–Naghdi (GN) equations. In Cartesian system, GN velocity profile depends on horizontal directions, x-direction and y-direction. The effect of vertical direction (z-direction) is also taken into consideration by applying weighting function in approximation. GN theory considers the effect of vertical acceleration and the consequent non-hydrostatic pressure. Moreover, in GN theory, the flow is rotational. The present study illustrates the application of GN equations to propagation of sloshing waves and solitary waves. For this purpose, GN equations solver is verified for the benchmark tests of Gaussian hump sloshing and solitary wave propagation in shallow basins. Analysis of the free surface sloshing of even harmonic components of an initial Gaussian hump demonstrates that the GN model gives predictions in satisfactory agreement with the linear analytical solutions. Discrepancies between the GN predictions and the linear analytical solutions arise from the effect of wave nonlinearities arising from the wave amplitude itself and wave-wave interactions. Numerically predicted solitary wave propagation indicates that the GN model produces simulations in good agreement with the analytical solution of the linearised wave theory. Comparison between the GN model numerical prediction and the result from perturbation analysis confirms that nonlinear interaction between solitary wave and a solid wall is satisfactorilly modelled. Moreover, solitary wave propagation at an angle to the x-axis and the interaction of solitary waves with each other are conducted to validate the developed model.
Abstract: Turbulence modelling is still evolving, and efforts are on to improve and develop numerical methods to simulate the real turbulence structures by using the empirical and experimental information. The monotonically integrated large eddy simulation (MILES) is an attractive approach for modelling turbulence in high Re flows, which is based on the solving of the unfiltered flow equations with no explicit sub-grid scale (SGS) model. In the current work, this approach has been used, and the action of the SGS model has been included implicitly by intrinsic nonlinear high-frequency filters built into the convection discretization schemes. The MILES solver is developed using the opensource CFD OpenFOAM libraries. The role of flux limiters schemes namely, Gamma, superBee, van-Albada and van-Leer, is studied in predicting turbulent statistical quantities for a fully developed channel flow with a friction Reynolds number, ReT = 180, and compared the numerical predictions with the well-established Direct Numerical Simulation (DNS) results for studying the wall generated turbulence. It is inferred from the numerical predictions that Gamma, van-Leer and van-Albada limiters produced more diffusion and overpredicted the velocity profiles, while superBee scheme reproduced velocity profiles and turbulence statistical quantities in good agreement with the reference DNS data in the streamwise direction although it deviated slightly in the spanwise and normal to the wall directions. The simulation results are further discussed in terms of the turbulence intensities and Reynolds stresses averaged in time and space to draw conclusion on the flux limiter schemes performance in OpenFOAM context.
Abstract: The aim of this paper is to report the different experimental studies, conducted in the laboratory, dealing with the flow in the presence of an obstacle lying in a rectangular hydraulic channel. Both subcritical and supercritical regimes are considered. Generally, when considering the theoretical problem of the free-surface flow, in a fluid domain of finite depth, due to the presence of an obstacle, we suppose that the water is an inviscid fluid, which means that there is no sheared velocity profile, but constant upstream. In a hydraulic channel, it is impossible to satisfy this condition. Indeed, water is a viscous fluid and its velocity is null at the bottom. The two configurations are presented, i.e. a flow over an obstacle and a towed obstacle in a resting fluid.
Abstract: The paper aims to extend the knowledge about jet behavior and jet interaction between five plane unventilated jets with large aspect ratio (AR). The distance between the single plane jets is two times the channel height. The experimental investigation applies 2D Particle Image Velocimetry (PIV) and static pressure measurements. Our study focuses on the influence of two different outlet nozzle geometries (triangular shape with 2 x 7.5° and blunt geometry) with respect to variation of Reynolds number from 5500 - 12000. It is shown that the outlet geometry has a major influence on the jet formation in terms of uniformity of velocity profiles downstream of the sudden expansion. Furthermore, we describe characteristic regions like converging region, merging region and combined region. The triangular outlet geometry generates most uniform velocity distributions in comparison to a blunt outlet nozzle geometry. The blunt outlet geometry shows an unstable behavior where the jets tend to attach to one side of the walls (ceiling) generating a large recirculation region on the opposite side. Static pressure measurements confirm the observation and indicate that the recirculation region is connected to larger pressure drop.
Abstract: Transportation of unrefined crude oil from the production unit to a refinery or large storage area by a pipeline is difficult due to the different properties of crude in various areas. Thus, the design of a crude oil pipeline is a very complex and time consuming process, when considering all the various parameters. There were three very important parameters that play a significant role in the transportation and processing pipeline design; these are: viscosity profile, temperature profile and the velocity profile of waxy crude oil through the crude oil pipeline. Knowledge of the Rheological computational technique is required for better understanding the flow behavior and predicting the flow profile in a crude oil pipeline. From these profile parameters, the material and the emulsion that is best suited for crude oil transportation can be predicted. Rheological computational fluid dynamic technique is a fast method used for designing flow profile in a crude oil pipeline with the help of computational fluid dynamics and rheological modeling. With this technique, the effect of fluid properties including shear rate range with temperature variation, degree of viscosity, elastic modulus and viscous modulus was evaluated under different conditions in a transport pipeline. In this paper, two crude oil samples was used, as well as a prepared emulsion with natural and synthetic additives, at different concentrations ranging from 1,000 ppm to 3,000 ppm. The rheological properties was then evaluated at a temperature range of 25 to 60 °C and which additive was best suited for transportation of crude oil is determined. Commercial computational fluid dynamics (CFD) has been used to generate the flow, velocity and viscosity profile of the emulsions for flow behavior analysis in crude oil transportation pipeline. This rheological CFD design can be further applied in developing designs of pipeline in the future.
Abstract: The heat and mass transfer characteristics of
contaminants in groundwater subjected to a biodegradation reaction
is analyzed by taking into account the thermal diffusion (Soret)
effects. This phenomenon is modulated mathematically by a
system of partial differential equations which govern the motion
of fluid (groundwater) and solid (contaminants) particles. The
numerical results are presented graphically for different values of
the parameters entering into the problem on the velocity profiles of
fluid, contaminants, temperature and concentration profile.
Abstract: In recent years, increasing the usage of railway transportations especially in developing countries caused more attention to control systems railway vehicles. Consequently, designing and implementing the modern control systems to improve the operating performance of trains and locomotives become one of the main concerns of researches. Dynamic braking systems is an important safety system which controls the amount of braking torque generated by traction motors, to keep the adhesion coefficient between the wheel-sets and rail road in optimum bound. Adhesion force has an important role to control the braking distance and prevent the wheels from slipping during the braking process. Cascade control structure is one of the best control methods for the wide range of industrial plants in the presence of disturbances and errors. This paper presents cascade control structure based on two forward simple controllers with two feedback loops to control the slip ratio and braking torque. In this structure, the inner loop controls the angular velocity and the outer loop control the longitudinal velocity of the locomotive that its dynamic is slower than the dynamic of angular velocity. This control structure by controlling the torque of DC traction motors, tries to track the desired velocity profile to access the predefined braking distance and to control the slip ratio. Simulation results are employed to show the effectiveness of the introduced methodology in dynamic braking system.
Abstract: In this study, laboratory experiments in open channel flows over a sand bed were conducted. A porous bed (sand bed) with porosity of ε=0.70 and porous thickness of s΄=3 cm was tested. Vertical distributions of velocity were evaluated by using a two-dimensional (2D) Particle Image Velocimetry (PIV). Velocity profiles are measured above the impermeable bed and above the sand bed for the same different total water heights (h= 6, 8, 10 and 12 cm) and for the same slope S=1.5. Measurements of mean velocity indicate the effects of the bed material used (sand bed) on the flow characteristics (Velocity distribution and Reynolds number) in comparison with those above the impermeable bed.
Abstract: The exploitation of flow pulsation in micro- and
mini-channels is a potentially useful technique for enhancing cooling
of high-end photonics and electronics systems. It is thought that
pulsation alters the thickness of the hydrodynamic and thermal
boundary layers, and hence affects the overall thermal resistance
of the heat sink. Although the fluid mechanics and heat transfer
are inextricably linked, it can be useful to decouple the parameters
to better understand the mechanisms underlying any heat transfer
enhancement. Using two-dimensional, two-component particle image
velocimetry, the current work intends to characterize the heat transfer
mechanisms in pulsating flow with a mean Reynolds number of
48 by experimentally quantifying the hydrodynamics of a generic
liquid-cooled channel geometry. Flows circulated through the test
section by a gear pump are modulated using a controller to achieve
sinusoidal flow pulsations with Womersley numbers of 7.45 and
2.36 and an amplitude ratio of 0.75. It is found that the transient
characteristics of the measured velocity profiles are dependent on the
speed of oscillation, in accordance with the analytical solution for
flow in a rectangular channel. A large velocity overshoot is observed
close to the wall at high frequencies, resulting from the interaction
of near-wall viscous stresses and inertial effects of the main fluid
body. The steep velocity gradients at the wall are indicative of
augmented heat transfer, although the local flow reversal may reduce
the upstream temperature difference in heat transfer applications.
While unsteady effects remain evident at the lower frequency, the
annular effect subsides and retreats from the wall. The shear rate at
the wall is increased during the accelerating half-cycle and decreased
during deceleration compared to steady flow, suggesting that the flow
may experience both enhanced and diminished heat transfer during
a single period. Hence, the thickness of the hydrodynamic boundary
layer is reduced for positively moving flow during one half of the
pulsation cycle at the investigated frequencies. It is expected that the
size of the thermal boundary layer is similarly reduced during the
cycle, leading to intervals of heat transfer enhancement.
Abstract: The construction of most coastal infrastructure developments around the world are usually made considering wave height, current velocities and river discharges; however, little effort has been paid to surveying sediment transport during dredging or the modification to currents outside the ports or marinas during and after the construction. This study shows a complete survey during the construction of one of the largest ports of the Gulf of Mexico. An anchored Acoustic Doppler Current Velocity profiler (ADCP), a towed ADCP and a combination of model outputs were used at the Veracruz port construction in order to describe the hourly sediment transport and current modifications in and out of the new port. Owing to the stability of the system the new port was construction inside Vergara Bay, a low wave energy system with a tidal range of up to 0.40 m. The results show a two-current system pattern within the bay. The north side of the bay has an anticyclonic gyre, while the southern part of the bay shows a cyclonic gyre. Sediment transport trajectories were made every hour using the anchored ADCP, a numerical model and the weekly data obtained from the towed ADCP within the entire bay. The sediment transport trajectories were carefully tracked since the bay is surrounded by coral reef structures which are sensitive to sedimentation rate and water turbidity. The survey shows that during dredging and rock input used to build the wave breaker sediments were locally added (< 2500 m2) and local currents disperse it in less than 4 h. While the river input located in the middle of the bay and the sewer system plant may add more than 10 times this amount during a rainy day or during the tourist season. Finally, the coastal line obtained seasonally with a drone suggests that the southern part of the bay has not been modified by the construction of the new port located in the northern part of the bay, owing to the two subsystem division of the bay.
Abstract: The aim of this study is to examine, through
experimentation in the laboratory, the supercritical flow in the
presence of an obstacle in a rectangular channel. The supercritical
regime in the whole hydraulic channel is achieved by adding a
convergent. We will observe the influence of the obstacle shape and
dimension on the characteristics of the supercritical flow, mainly the
free-surface elevation and the velocity profile. The velocity
measurements have been conducted with the one dimension laser
anemometry technique.
Abstract: Geological and tectonic framework indicates that
Bangladesh is one of the most seismically active regions in the world.
The Bengal Basin is at the junction of three major interacting plates:
the Indian, Eurasian, and Burma Plates. Besides there are many
active faults within the region, e.g. the large Dauki fault in the north.
The country has experienced a number of destructive earthquakes due
to the movement of these active faults. Current seismic provisions of
Bangladesh are mostly based on earthquake data prior to the 1990.
Given the record of earthquakes post 1990, there is a need to revisit
the design provisions of the code. This paper compares the base shear
demand of three major cities in Bangladesh: Dhaka (the capital city),
Sylhet, and Chittagong for earthquake scenarios of magnitudes
7.0MW, 7.5MW, 8.0MW, and 8.5MW using a stochastic model. In
particular, the stochastic model allows the flexibility to input region
specific parameters such as shear wave velocity profile (that were
developed from Global Crustal Model CRUST2.0) and include the
effects of attenuation as individual components. Effects of soil
amplification were analysed using the Extended Component
Attenuation Model (ECAM). Results show that the estimated base
shear demand is higher in comparison with code provisions leading to
the suggestion of additional seismic design consideration in the study
regions.
Abstract: The understanding of geotechnical characteristics of
near-surface material and the effects of the groundwater is very
important problem in such as site studies. For showing the relations
between seismic data and groundwater, we selected about 25 km2 as
the study area. It has been presented which is a detailed work of
seismic data and groundwater depths of Gokpinar Damp area.
Seismic waves velocity (Vp and Vs) are very important parameters
showing the soil properties. The seismic records were used the
method of the multichannel analysis of surface waves near area of
Gokpinar Damp area. Sixty sites in this area have been investigated
with survey lines about 60 m in length. MASW (Multichannel
analysis of surface wave) method has been used to generate onedimensional
shear wave velocity profile at locations. These shear
wave velocities are used to estimate equivalent shear wave velocity in
the study area at every 2 and 5 m intervals up to a depth of 45 m.
Levels of equivalent shear wave velocity of soil are used the
classified of the study area. After the results of the study, it must be
considered as components of urban planning and building design of
Gokpinar Damp area, Denizli and the application and use of these
results should be required and enforced by municipal authorities.
Abstract: The steady flow of a second order fluid through
constricted tube with slip velocity at wall is modeled and analyzed
theoretically. The governing equations are simplified by implying no
slip in radial direction. Based on Karman Pohlhausen procedure
polynomial solution for axial velocity profile is presented.
Expressions for pressure gradient, shear stress, separation and
reattachment points, and radial velocity are also calculated. The
effect of slip and no slip velocity on magnitude velocity, shear stress,
and pressure gradient are discussed and depicted graphically. It is
noted that when Reynolds number increases magnitude velocity of
the fluid decreases in both slip and no slip conditions. It is also found
that the wall shear stress, separation, and reattachment points are
strongly affected by Reynolds number.
Abstract: Due to the interference effects, the intrinsic
aerodynamic parameters obtained from the individual component
testing are always fundamentally different than those obtained for
complete model testing. Consideration and limitation for such testing
need to be taken into account in any design work related to the
component buildup method. In this paper, the scaled model of a
straight rectangular canard of a hybrid buoyant aircraft is tested at 50
m/s in IIUM-LSWT (Low Speed Wind Tunnel). Model and its
attachment with the balance are kept rigid to have results free from
the aeroelastic distortion. Based on the velocity profile of the test
section’s floor; the height of the model is kept equal to the
corresponding boundary layer displacement. Balance measurements
provide valuable but limited information of overall aerodynamic
behavior of the model. Zero lift coefficient is obtained at -2.2o and
the corresponding drag coefficient was found to be less than that at
zero angle of attack. As a part of the validation of low fidelity tool,
plot of lift coefficient plot was verified by the experimental data and
except the value of zero lift coefficients, the overall trend has under
predicted the lift coefficient. Based on this comparative study, a
correction factor of 1.36 is proposed for lift curve slope obtained
from the panel method.
Abstract: MHD chemically reacting viscous fluid flow towards
a vertical surface with slip and convective boundary conditions has
been conducted. The temperature and the chemical species
concentration of the surface and the velocity of the external flow are
assumed to vary linearly with the distance from the vertical surface.
The governing differential equations are modeled and transformed
into systems of ordinary differential equations, which are then solved
numerically by a shooting method. The effects of various parameters
on the heat and mass transfer characteristics are discussed. Graphical
results are presented for the velocity, temperature, and concentration
profiles whilst the skin-friction coefficient and the rate of heat and
mass transfers near the surface are presented in tables and discussed.
The results revealed that increasing the strength of the magnetic field
increases the skin-friction coefficient and the rate of heat and mass
transfers toward the surface. The velocity profiles are increased
towards the surface due to the presence of the Lorenz force, which
attracts the fluid particles near the surface. The rate of chemical
reaction is seen to decrease the concentration boundary layer near the
surface due to the destructive chemical reaction occurring near the
surface.