Abstract: In order to study the aerodynamic performance of a
semi-flexible membrane wing, Fluid-Structure Interaction simulations
have been performed. The fluid problem has been modeled using
two different approaches which are the vortex panel method and the
numerical solution of the Navier-Stokes equations. Nonlinear analysis
of the structural problem is performed using the Finite Element
Method. Comparison between the two fluid solvers has been made.
Aerodynamic performance of the wing is discussed regarding its
lift and drag coefficients and they are compared with those of the
equivalent rigid wing.
Abstract: The computational fluid dynamics (CFD) study of
stirred tank with the air-water interface are carried out in the presence
of different types of the impeller and with or without baffles. A
multiple reference frame (MRF) approach with the volume of fluid
(VOF) method is used to capture the air-water interface. The RANS
(Reynolds Averaged Navier-Stokes) equations with k-ε turbulence
model are solved to predict the flow behavior of water and air phase
which are treated as a different phases. The predicted results have
shown that the VOF method is able to capture the interface in the
unbaffled tank. While, the VOF method is showing an unfeasible
results in the baffled tank with high rotational impeller speed. For
continuous stirred tank, the air-water interface is disturbed by the
inflow and the level of water is also increased with time.
Abstract: This study presents the numerical simulation of three-dimensional incompressible steady and laminar fluid flow and conjugate heat transfer of a trapezoidal microchannel heat sink using water as a cooling fluid in a silicon substrate. Navier-Stokes equations with conjugate energy equation are discretized by finite-volume method. We perform numerical computations for a range of 50 ≦ Re ≦ 600, 0.05W ≦ P ≦ 0.8W, 20W/cm2 ≦q"≦ 40W/cm2. The present study demonstrates the numerical optimization of a trapezoidal microchannel heat sink design using the response surface methodology (RSM) and the genetic algorithm method (GA). The results show that the average Nusselt number increases with an increase in the Reynolds number or pumping power, and the thermal resistance decreases as the pumping power increases. The thermal resistance of a trapezoidal microchannel is minimized for a constant heat flux and constant pumping power.
Abstract: The aim of this work is to analyze a viscous flow
around the axisymmetric blunt body taken into account the mesh size
both in the free stream and into the boundary layer. The resolution of
the Navier-Stokes equations is realized by using the finite volume
method to determine the flow parameters and detached shock
position. The numerical technique uses the Flux Vector Splitting
method of Van Leer. Here, adequate time stepping parameter, CFL
coefficient and mesh size level are selected to ensure numerical
convergence. The effect of the mesh size is significant on the shear
stress and velocity profile. The best solution is obtained with using a
very fine grid. This study enabled us to confirm that the
determination of boundary layer thickness can be obtained only if the
size of the mesh is lower than a certain value limits given by our
calculations.
Abstract: The typical insects employ a flapping-wing mode of flight. The numerical simulations on free flight of a model fruit fly (Re=143) including hovering and are presented in this paper. Unsteady aerodynamics around a flapping insect is studied by solving the three-dimensional Newtonian dynamics of the flyer coupled with Navier-Stokes equations. A hybrid-grid scheme (Generalized Finite Difference Method) that combines great geometry flexibility and accuracy of moving boundary definition is employed for obtaining flow dynamics. The results show good points of agreement and consistency with the outcomes and analyses of other researchers, which validate the computational model and demonstrate the feasibility of this computational approach on analyzing fluid phenomena in insect flight. The present modeling approach also offers a promising route of investigation that could complement as well as overcome some of the limitations of physical experiments in the study of free flight aerodynamics of insects. The results are potentially useful for the design of biomimetic flapping-wing flyers.
Abstract: The external flow simulation of a flying car at take off phase is a daunting task owing to the fact that the prediction of the transient unsteady flow features during its deployment phase is very complex. In this paper 3D numerical simulations of external flow of Ferrari F430 proposed flying car with different NACA 9618 rectangular wings have been carried. Additionally, the aerodynamics characteristics have been generated for optimizing its geometry for achieving the minimum take off velocity with better overall performance in both road and air. The three-dimensional standard k-omega turbulence model has been used for capturing the intrinsic flow physics during the take off phase. In the numerical study, a fully implicit finite volume scheme of the compressible, Reynolds-Averaged, Navier-Stokes equations is employed. Through the detailed parametric analytical studies we have conjectured that Ferrari F430 flying car facilitated with high wings having three different deployment histories during the take off phase is the best choice for accomplishing its better performance for the commercial applications.
Abstract: The practical application of the Computational Fluid Dynamics (CFD), for predicting the flow pattern around Multipurpose Amphibious Vehicle (MAV) hull has made much progress over the last decade. Today, several of the CFD tools play an important role in the land and water going vehicle hull form design. CFD has been used for analysis of MAV hull resistance, sea-keeping, maneuvering and investigating its variation when changing the hull form due to varying its parameters, which represents a very important task in the principal and final design stages. Resistance analysis based on CFD (Computational Fluid Dynamics) simulation has become a decisive factor in the development of new, economically efficient and environmentally friendly hull forms. Three-dimensional finite volume method (FVM) based on Reynolds Averaged Navier-Stokes equations (RANS) has been used to simulate incompressible flow around three types of MAV hull bow models in steady-state condition. Finally, the flow structure and streamlines, friction and pressure resistance and velocity contours of each type of hull bow will be compared and discussed.
Abstract: The effect of an axial electric field on the capillary instability of a cylindrical interface in the presence of heat and mass transfer 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, heat transfer capillary number, conductivity ratio, permittivity ratio on the stability of the system. It has been observed that the axial electric field and heat and mass transfer both have stabilizing effect on the stability of the system.
Abstract: To solve saddle point systems efficiently, several preconditioners have been published. There are many methods for constructing preconditioners for linear systems from saddle point problems, for instance, the relaxed dimensional factorization (RDF) preconditioner and the augmented Lagrangian (AL) preconditioner are used for both steady and unsteady Navier-Stokes equations. In this paper we compare the RDF preconditioner with the modified AL (MAL) preconditioner to show which is more effective to solve Navier-Stokes equations. Numerical experiments indicate that the MAL preconditioner is more efficient and robust, especially, for moderate viscosities and stretched grids in steady problems. For unsteady cases, the convergence rate of the RDF preconditioner is slightly faster than the MAL perconditioner in some circumstances, but the parameter of the RDF preconditioner is more sensitive than the MAL preconditioner. Moreover the convergence rate of the MAL preconditioner is still quite acceptable. Therefore we conclude that the MAL preconditioner is more competitive than the RDF preconditioner. These experiments are implemented with IFISS package.
Abstract: In this work we study the thermodynamic behavior of some ventilated facades under summer operating conditions in Southern Spain. Under these climatic conditions, indoor comfort implies a high energetic demand due to high temperatures that usually are reached in this season in the considered geographical area.
The aim of this work is to determine if during summer operating conditions in Southern Spain, ventilated façades provide some energy saving compared to the non-ventilated façades and to deduce their behavior patterns in terms of energy efficiency.
The modelization of the air flow in the channel has been performed by using Navier-Stokes equations for thermodynamic flows. Numerical simulations have been carried out with a 2D Finite Element approach.
This way, we analyze the behavior of ventilated façades under different weather conditions as variable wind, variable temperature and different levels of solar irradiation.
CFD computations show the combined effect of the shading of the external wall and the ventilation by the natural convection into the air gap achieve a reduction of the heat load during the summer period. This reduction has been evaluated by comparing the thermodynamic performances of two ventilated and two unventilated façades with the same geometry and thermophysical characteristics.
Abstract: This work presents a numerical simulation of the interaction of an incident shock wave propagates from the left to the right with a cone placed in a tube at shock. The Mathematical model is based on a non stationary, viscous and axisymmetric flow. The Discretization of the Navier-stokes equations is carried out by the finite volume method in the integral form along with the Flux Vector Splitting method of Van Leer. Here, adequate combination of time stepping parameter, CFL coefficient and mesh size level is selected to ensure numerical convergence. The numerical simulation considers a shock tube filled with air. The incident shock wave propagates to the right with a determined Mach number and crosses the cone by leaving behind it a stationary detached shock wave in front of the nose cone. This type of interaction is observed according to the time of flow.
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.
Abstract: A three-dimensional and pulsatile blood flow in the left ventricle of heart model has been studied numerically. The geometry was derived from a simple approximation of the left ventricle model and the numerical simulations were obtained using a formulation of the Navier-Stokes equations. In this study, simulation was used to investigate the pattern of flow velocity in 3D model of heart with consider the left ventricle based on critical parameter of blood under steady condition. Our results demonstrate that flow velocity focused from mitral valve channel and continuous linearly to left ventricle wall but this skewness progresses into outside wall in atrium through aortic valve with random distribution that is irregular due to force subtract from ventricle wall during cardiac cycle. The findings are the prediction of the behavior of the blood flow velocity pattern in steady flow condition which can assist the medical practitioners in their decision on the patients- treatments.
Abstract: Acoustic function plays an important role in
aerodynamic mechanical engineering. It can classify the kind of
air-vehicle such as subsonic or supersonic. Acoustic velocity
relates with velocity and Mach number. Mach number relates
again acoustic stability or instability condition. Mach number
plays an important role in growth or decay in energy system.
Acoustic is a function of temperature and temperature is directly
proportional to pressure. If we control the pressure, we can control
acoustic function. To get pressure stability condition, we apply
Navier-Stokes equations.
Abstract: The motion of a sphere moving along the axis of a
rotating viscous fluid is studied at high Reynolds numbers and
moderate values of Taylor number. The Higher Order Compact
Scheme is used to solve the governing Navier-Stokes equations. The
equations are written in the form of Stream function, Vorticity
function and angular velocity which are highly non-linear, coupled
and elliptic partial differential equations. The flow is governed by
two parameters Reynolds number (Re) and Taylor number (T). For
very low values of Re and T, the results agree with the available
experimental and theoretical results in the literature. The results are
obtained at higher values of Re and moderate values of T and
compared with the experimental results. The results are fourth order
accurate.
Abstract: The simulation of extrusion process is studied widely
in order to both increase products and improve quality, with broad
application in wire coating. The annular tube-tooling extrusion was
set up by a model that is termed as Navier-Stokes equation in
addition to a rheological model of differential form based on singlemode
exponential Phan-Thien/Tanner constitutive equation in a twodimensional
cylindrical coordinate system for predicting the
contraction point of the polymer melt beyond the die. Numerical
solutions are sought through semi-implicit Taylor-Galerkin pressurecorrection
finite element scheme. The investigation was focused on
incompressible creeping flow with long relaxation time in terms of
Weissenberg numbers up to 200. The isothermal case was considered
with surface tension effect on free surface in extrudate flow and no
slip at die wall. The Stream Line Upwind Petrov-Galerkin has been
proposed to stabilize solution. The structure of mesh after die exit
was adjusted following prediction of both top and bottom free
surfaces so as to keep the location of contraction point around one
unit length which is close to experimental results. The simulation of
extrusion process is studied widely in order to both increase products
and improve quality, with broad application in wire coating. The
annular tube-tooling extrusion was set up by a model that is termed
as Navier-Stokes equation in addition to a rheological model of
differential form based on single-mode exponential Phan-
Thien/Tanner constitutive equation in a two-dimensional cylindrical
coordinate system for predicting the contraction point of the polymer
melt beyond the die. Numerical solutions are sought through semiimplicit
Taylor-Galerkin pressure-correction finite element scheme.
The investigation was focused on incompressible creeping flow with
long relaxation time in terms of Weissenberg numbers up to 200. The
isothermal case was considered with surface tension effect on free
surface in extrudate flow and no slip at die wall. The Stream Line
Upwind Petrov-Galerkin has been proposed to stabilize solution. The
structure of mesh after die exit was adjusted following prediction of
both top and bottom free surfaces so as to keep the location of
contraction point around one unit length which is close to
experimental results.
Abstract: Computational simulation of steam flow and heat transfer in power plant condensers on the basis of the threedimensional mathematical model for the flow through porous media is presented. In order to solve the mathematical model of steam flow and heat transfer in power plant condensers, the Streamline Upwind Petrov-Galerkin finite element method is applied. By comparison of the results of simulation with experimental results about an experimental condenser, it is confirmed that SUPG finite element method can be successfully applied for solving the three-dimensional mathematical model of steam flow and heat transfer in power plant condensers.
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 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.