Abstract: The significant effects of the interactions between the
system boundaries and the near wall molecules in miniaturized
gaseous devices lead to the formation of the Knudsen layer in which
the Navier-Stokes-Fourier (NSF) equations fail to predict the correct
associated phenomena. In this paper, the well-known lattice
Boltzmann method (LBM) is employed to simulate the fluid flow and
heat transfer processes in rarefied gaseous micro media. Persuaded
by the problematic deficiency of the LBM in capturing the Knudsen
layer phenomena, present study tends to concentrate on the effective
molecular mean free path concept the main essence of which is to
compensate the incapability of this mesoscopic method in dealing
with the momentum and energy transport within the above mentioned
kinetic boundary layer. The results show qualitative and quantitative
accuracy comparable to the solutions of the linearized Boltzmann
equation or the DSMC data for the Knudsen numbers of O (1) .
Abstract: Natural convection heat transfer from a heated
horizontal semi-circular cylinder (flat surface upward) has been
investigated for the following ranges of conditions; Grashof number,
and Prandtl number. The governing partial differential equations
(continuity, Navier-Stokes and energy equations) have been solved
numerically using a finite volume formulation. In addition, the role of
the type of the thermal boundary condition imposed at cylinder
surface, namely, constant wall temperature (CWT) and constant heat
flux (CHF) are explored. Natural convection heat transfer from a
heated horizontal semi-circular cylinder (flat surface upward) has
been investigated for the following ranges of conditions; Grashof
number, and Prandtl number, . The governing partial differential
equations (continuity, Navier-Stokes and energy equations) have
been solved numerically using a finite volume formulation. In
addition, the role of the type of the thermal boundary condition
imposed at cylinder surface, namely, constant wall temperature
(CWT) and constant heat flux (CHF) are explored. The resulting flow
and temperature fields are visualized in terms of the streamline and
isotherm patterns in the proximity of the cylinder. The flow remains
attached to the cylinder surface over the range of conditions spanned
here except that for and ; at these conditions, a separated flow
region is observed when the condition of the constant wall
temperature is prescribed on the surface of the cylinder. The heat
transfer characteristics are analyzed in terms of the local and average
Nusselt numbers. The maximum value of the local Nusselt number
always occurs at the corner points whereas it is found to be minimum
at the rear stagnation point on the flat surface. Overall, the average
Nusselt number increases with Grashof number and/ or Prandtl
number in accordance with the scaling considerations. The numerical
results are used to develop simple correlations as functions of
Grashof and Prandtl number thereby enabling the interpolation of the
present numerical results for the intermediate values of the Prandtl or
Grashof numbers for both thermal boundary conditions.
Abstract: We have studied the migration of a charged permeable aggregate in electrolyte under the influence of an axial electric field and pressure gradient. The migration of the positively charged aggregate leads to a deformation of the anionic cloud around it. The hydrodynamics of the aggregate is governed by the interaction of electroosmotic flow in and around the particle, hydrodynamic friction and electric force experienced by the aggregate. We have computed the non-linear Nernest-Planck equations coupled with the Dracy- Brinkman extended Navier-Stokes equations and Poisson equation for electric field through a finite volume method. The permeability of the aggregate enable the counterion penetration. The penetration of counterions depends on the volume charge density of the aggregate and ionic concentration of electrolytes at a fixed field strength. The retardation effect due to the double layer polarization increases the drag force compared to an uncharged aggregate. Increase in migration sped from the electrophretic velocity of the aggregate produces further asymmetry in charge cloud and reduces the electric body force exerted on the particle. The permeability of the particle have relatively little influence on the electric body force when Double layer is relatively thin. The impact of the key parameters of electrokinetics on the hydrodynamics of the aggregate is analyzed.
Abstract: The performance of Advection Upstream Splitting
Method AUSM schemes are evaluated against experimental flow
fields at different Mach numbers and results are compared with
experimental data of subsonic, supersonic and hypersonic flow fields.
The turbulent model used here is SST model by Menter. The
numerical predictions include lift coefficient, drag coefficient and
pitching moment coefficient at different mach numbers and angle of
attacks. This work describes a computational study undertaken to
compute the Aerodynamic characteristics of different air vehicles
configurations using a structured Navier-Stokes computational
technique. The CFD code bases on the idea of upwind scheme for the
convective (convective-moving) fluxes. CFD results for GLC305
airfoil and cone cylinder tail fined missile calculated on above
mentioned turbulence model are compared with the available data.
Wide ranges of Mach number from subsonic to hypersonic speeds are
simulated and results are compared. When the computation is done
by using viscous turbulence model the above mentioned coefficients
have a very good agreement with the experimental values. AUSM
scheme is very efficient in the regions of very high pressure gradients
like shock waves and discontinuities. The AUSM versions simulate
the all types of flows from lower subsonic to hypersonic flow without
oscillations.
Abstract: A numerical study is made of laminar, unsteady flow
behind a rotationally oscillating circular cylinder using a recently
developed higher order compact (HOC) scheme. The stream function
vorticity formulation of Navier-Stokes (N-S) equations in cylindrical
polar coordinates are considered as the governing equations. The
temporal behaviour of vortex formation and relevant streamline
patterns of the flow are scrutinized over broad ranges of two
externally specified parameters namely dimensionless forced
oscillating frequency Sf and dimensionless peak rotation rate αm for
the Reynolds-s number Re = 200. Excellent agreements are found
both qualitatively and quantitatively with the existing experimental
and standard numerical results.
Abstract: The development of biomimetic micro-aerial-vehicles
(MAVs) with flapping wings is the future trend in military/domestic
field. The successful flight of MAVs is strongly related to the
understanding of unsteady aerodynamic performance of low Reynolds
number airfoils under dynamic flapping motion. This study explored
the effects of flapping frequency, stroke amplitude, and the inclined
angle of stroke plane on lift force and thrust force of a bio-inspiration
corrugated airfoil with 33 full factorial design of experiment and
ANOVA analysis. Unsteady vorticity flows over a corrugated thin
airfoil executing flapping motion are computed with time-dependent
two-dimensional laminar incompressible Reynolds-averaged
Navier-Stokes equations with the conformal hybrid mesh. The tested
freestream Reynolds number based on the chord length of airfoil as
characteristic length is fixed of 103. The dynamic mesh technique is
applied to model the flapping motion of a corrugated airfoil. Instant
vorticity contours over a complete flapping cycle clearly reveals the
flow mechanisms for lift force generation are dynamic stall, rotational
circulation, and wake capture. The thrust force is produced as the
leading edge vortex shedding from the trailing edge of airfoil to form a
reverse von Karman vortex. Results also indicated that the inclined
angle is the most significant factor on both the lift force and thrust
force. There are strong interactions between tested factors which mean
an optimization study on parameters should be conducted in further
runs.
Abstract: Channel junctions can be analyzed in two ways of
division (lateral intake) and combined flows (confluence). The
present paper investigates 3D flow pattern at lateral intake using
Navier-Stokes equation and κ -ε (RNG) turbulent model. The
equations are solved by Finite-Volume Method (FVM) and results
are compared with the experimental data of (Barkdoll, B.D., 1997)
to test the validity of the findings. Comparison of the results with
the experimental data indicated a close proximity between the two
sets of data which suggest a very close simulation. Results further
indicated an inverse relation between the effects of discharge ratio
( r Q ) on the length and width of the separation zone. In other words,
as the discharge ration increases, the length and width of separation
zone decreases.
Abstract: The major focus of this work was to characterize hydrodynamics in a packed-bed with and without static mixer by using Computational Fluid Dynamic (CFD). The commercial software: COMSOL MULTIPHYSICSTM Version 3.3 was used to simulate flow fields of mixed-gas reactants i.e. CO and H2. The packed-bed was a single tube with the inside diameter of 0.8 cm and the length of 1.2 cm. The static mixer was inserted inside the tube. The number of twisting elements was 1 with 0.8 cm in diameter and 1.2 cm in length. The packed-bed with and without static mixer were both packed with approximately 700 spherical structures representing catalyst pellets. Incompressible Navier-Stokes equations were used to model the gas flow inside the beds at steady state condition, in which the inlet Reynolds Number (Re) was 2.31. The results revealed that, with the insertion of static mixer, the gas was forced to flow radially inward and outward between the central portion of the tube and the tube wall. This could help improving the overall performance of the packed-bed, which could be utilized for heterogeneous catalytic reaction such as reforming and Fischer- Tropsch reactions.
Abstract: Solution to unsteady Navier-Stokes equation by Splitting method in physical orthogonal algebraic curvilinear coordinate system, also termed 'Non-linear grid system' is presented. The linear terms in Navier-Stokes equation are solved by Crank- Nicholson method while the non-linear term is solved by the second order Adams-Bashforth method. This work is meant to bring together the advantage of Splitting method as pressure-velocity solver of higher efficiency with the advantage of consuming Non-linear grid system which produce more accurate results in relatively equal number of grid points as compared to Cartesian grid. The validation of Splitting method as a solution of Navier-Stokes equation in Nonlinear grid system is done by comparison with the benchmark results for lid driven cavity flow by Ghia and some case studies including Backward Facing Step Flow Problem.
Abstract: The RANS method with Saffman-s turbulence model
was employed to solve the time-dependent turbulent Navier-Stokes
and energy equations for oscillating pipe flows. The method of
partial sums of the Fourier series is used to analyze the harmonic
velocity and temperature results. The complete structures of the
oscillating pipe flows and the averaged Nusselt numbers on the tube
wall are provided by numerical simulation over wide ranges of ReA
and ReR. Present numerical code is validated by comparing the
laminar flow results to analytic solutions and turbulence flow results
to published experimental data at lower and higher Reynolds
numbers respectively. The effects of ReA and ReR on the velocity,
temperature and Nusselt number distributions have been di scussed.
The enhancement of the heat transfer due to oscillating flows has
also been presented. By the way of analyzing the overall Nusselt
number over wide ranges of the Reynolds number Re and Keulegan-
Carpenter number KC, the optimal ratio of the tube diameter over
the oscillation amplitude is obtained based on the existence of a
nearly constant optimal KC number. The potential application of the
present results in sea water cooling has also been discussed.
Abstract: Shape optimization of the airfoil with high aspect ratio
of long endurance unmanned aerial vehicle (UAV) is performed by the
multi-objective optimization technology coupled with computational
fluid dynamics (CFD). For predicting the aerodynamic characteristics
around the airfoil the high-fidelity Navier-Stokes solver is employed
and SMOGA (Simple Multi-Objective Genetic Algorithm), which is
developed by authors, is used for solving the multi-objective
optimization problem. To obtain the optimal solutions of the design
variable (i.e., sectional airfoil profile, wing taper ratio and sweep) for
high performance of UAVs, both the lift and lift-to-drag ratio are
maximized whereas the pitching moment should be minimized,
simultaneously. It is found that the lift force and lift-to-drag ratio are
linearly dependent and a unique and dominant solution are existed.
However, a trade-off phenomenon is observed between the lift-to-drag
ratio and pitching moment. As the result of optimization, sixty-five
(65) non-dominated Pareto individuals at the cutting edge of design
spaces that is decided by airfoil shapes can be obtained.
Abstract: Direct numerical simulation (DNS) is used to study the evolution of a boundary layer that was laminar initially followed by separation and then reattachment owing to generation of turbulence. This creates a closed region of recirculation, known as the laminar-separation bubble. The present simulation emulates the flow environment encountered in a modern LP turbine blade, where a laminar separation bubble may occur on the suction surface. The unsteady, incompressible three-dimensional (3-D) Navier-Stokes (NS) equations have been solved over a flat plate in the Cartesian coordinates. The adverse pressure gradient, which causes the flow to separate, is created by a boundary condition. The separated shear layer undergoes transition through appearance of ╬ø vortices, stretching of these create longitudinal streaks. Breakdown of the streaks into small and irregular structures makes the flow turbulent downstream.
Abstract: Axial Flux Permanent Magnet (AFPM) Machines require effective cooling due to their high power density. The detrimental effects of overheating such as degradation of the insulation materials, magnets demagnetization, and increase of Joule losses are well known. This paper describes the CFD simulations performed on a test rig model of an air cooled Axial Flux Permanent Magnet (AFPM) generator built at Durham University to identify the temperatures and heat transfer coefficient on the stator. The Reynolds Averaged Navier-Stokes and the Energy equations are solved and the flow pattern and heat transfer developing inside the machine are described. The Nusselt number on the stator surfaces has been found. The dependency of the heat transfer on the flow field is described temperature field obtained. Tests on an experimental are undergoing in order to validate the CFD results.
Abstract: A numerical simulation of vortex-induced vibration of
a 2-dimensional elastic circular cylinder with two degree of freedom
under the uniform flow is calculated when Reynolds is 200.
2-dimensional incompressible Navier-Stokes equations are solved
with the space-time finite element method, the equation of the cylinder
motion is solved with the new explicit integral method and the mesh
renew is achieved by the spring moving mesh technology. Considering
vortex-induced vibration with the low reduced damping parameter, the
variety trends of the lift coefficient, the drag coefficient, the
displacement of cylinder are analyzed under different oscillating
frequencies of cylinder. The phenomena of locked-in, beat and
phases-witch were captured successfully. The evolution of vortex
shedding from the cylinder with time is discussed. There are very
similar trends in characteristics between the results of the one degree
of freedom cylinder model and that of the two degree of freedom
cylinder model. The streamwise vibrations have a certain effect on the
lateral vibrations and their characteristics.
Abstract: Numerical integration of initial boundary problem for advection equation in 3 ℜ is considered. The method used is
conditionally stable semi-Lagrangian advection scheme with high order interpolation on unstructured mesh. In order to increase time step integration the BFECC method with limiter TVD correction is used. The method is adopted on parallel graphic processor unit environment using NVIDIA CUDA and applied in Navier-Stokes solver. It is shown that the calculation on NVIDIA GeForce 8800
GPU is 184 times faster than on one processor AMDX2 4800+ CPU. The method is extended to the incompressible fluid dynamics solver. Flow over a Cylinder for 3D case is compared to the experimental data.
Abstract: In this paper, numerical simulations are performed to investigate the effect of disturbance block on flow field of the classical square lid-driven cavity. Attentions are focused on vortex formation and studying the effect of block position on its structure. Corner vortices are different upon block position and new vortices are produced because of the block. Finite volume method is used to solve Navier-Stokes equations and PISO algorithm is employed for the linkage of velocity and pressure. Verification and grid independency of results are reported. Stream lines are sketched to visualize vortex structure in different block positions.
Abstract: We study bifurcation structure of the zonal jet flow the
streamfunction of which is expressed by a single spherical harmonics
on a rotating sphere. In the non-rotating case, we find that a steady
traveling wave solution arises from the zonal jet flow through Hopf
bifurcation. As the Reynolds number increases, several traveling
solutions arise only through the pitchfork bifurcations and at high
Reynolds number the bifurcating solutions become Hopf unstable. In
the rotating case, on the other hand, under the stabilizing effect of
rotation, as the absolute value of rotation rate increases, the number
of the bifurcating solutions arising from the zonal jet flow decreases
monotonically. We also carry out time integration to study unsteady
solutions at high Reynolds number and find that in the non-rotating
case the unsteady solutions are chaotic, while not in the rotating cases
calculated. This result reflects the general tendency that the rotation
stabilizes nonlinear solutions of Navier-Stokes equations.
Abstract: A two dimensional numerical simulation has been
performed for incompressible and compressible fluid flow through
microchannels in slip flow regime. The Navier-Stokes equations have
been solved in conjunction with Maxwell slip conditions for
modeling flow field associated with slip flow regime. The wall
roughness is simulated with triangular microelements distributed on
wall surfaces to study the effects of roughness on fluid flow. Various
Mach and Knudsen numbers are used to investigate the effects of
rarefaction as well as compressibility. It is found that rarefaction has
more significant effect on flow field in microchannels with higher
relative roughness. It is also found that compressibility has more
significant effects on Poiseuille number when relative roughness
increases. In addition, similar to incompressible models the increase
in average fRe is more significant at low Knudsen number flows but
the increase of Poiseuille number duo to relative roughness is sharper
for compressible models. The numerical results have also validated
with some available theoretical and experimental relations and good
agreements have been seen.
Abstract: Due to the three- dimensional flow pattern interacting with bed material, the process of local scour around bridge piers is complex. Modeling 3D flow field and scour hole evolution around a bridge pier is more feasible nowadays because the computational cost and computational time have significantly decreased. In order to evaluate local flow and scouring around a bridge pier, a completely three-dimensional numerical model, SSIIM program, was used. The model solves 3-D Navier-Stokes equations and a bed load conservation equation. The model was applied to simulate local flow and scouring around a bridge pier in a large natural river with four piers. Computation for 1 day of flood condition was carried out to predict the maximum local scour depth. The results show that the SSIIM program can be used efficiently for simulating the scouring in natural rivers. The results also showed that among the various turbulence models, the k-ω model gives more reasonable results.
Abstract: The effect of porous medium on the capillary instability of a cylindrical interface in the presence of axial electric field has been investigated using viscous potential flow theory. In viscous potential flow, the viscous term in Navier-Stokes equation vanishes as
vorticity is zero but viscosity is not zero. Viscosity enters through normal stress balance in the viscous potential flow theory and tangential stresses are not considered. A dispersion relation that accounts for the growth of axisymmetric waves is derived and stability is discussed theoretically as well as numerically. Stability criterion is given by critical value of applied electric field as well as critical wave number. Various graphs have been drawn to show the effect of various physical parameters such as electric field, viscosity ratio, permittivity ratio on the stability of the system. It has been observed that the axial electric field and porous medium both have stabilizing effect on the stability of the system.