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: 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: Fossil fuels are the major source to meet the world
energy requirements but its rapidly diminishing rate and adverse
effects on our ecological system are of major concern. Renewable
energy utilization is the need of time to meet the future challenges.
Ocean energy is the one of these promising energy resources. Threefourths
of the earth-s surface is covered by the oceans. This enormous
energy resource is contained in the oceans- waters, the air above the
oceans, and the land beneath them. The renewable energy source of
ocean mainly is contained in waves, ocean current and offshore solar
energy. Very fewer efforts have been made to harness this reliable
and predictable resource. Harnessing of ocean energy needs detail
knowledge of underlying mathematical governing equation and their
analysis. With the advent of extra ordinary computational resources
it is now possible to predict the wave climatology in lab simulation.
Several techniques have been developed mostly stem from numerical
analysis of Navier Stokes equations. This paper presents a brief over
view of such mathematical model and tools to understand and
analyze the wave climatology. Models of 1st, 2nd and 3rd generations
have been developed to estimate the wave characteristics to assess the
power potential. A brief overview of available wave energy
technologies is also given. A novel concept of on-shore wave energy
extraction method is also presented at the end. The concept is based
upon total energy conservation, where energy of wave is transferred
to the flexible converter to increase its kinetic energy. Squeezing
action by the external pressure on the converter body results in
increase velocities at discharge section. High velocity head then can
be used for energy storage or for direct utility of power generation.
This converter utilizes the both potential and kinetic energy of the
waves and designed for on-shore or near-shore application. Increased
wave height at the shore due to shoaling effects increases the
potential energy of the waves which is converted to renewable
energy. This approach will result in economic wave energy
converter due to near shore installation and more dense waves due to
shoaling. Method will be more efficient because of tapping both
potential and kinetic energy of the waves.
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: Numerical studies have been carried out using a two
dimensional code to examine the influence of pressure / thrust
transient of solid propellant rockets at liftoff. This code solves
unsteady Reynolds-averaged thin-layer Navier–Stokes equations by
an implicit LU-factorization time-integration method. The results
from the parametric study indicate that when the port is narrow there
is a possibility of increase in pressure / thrust-rise rate due to
relatively high flame spread rate. Parametric studies further reveal
that flame spread rate can be altered by altering the propellant
properties, igniter jet characteristics and nozzle closure burst pressure
without altering the grain configuration and/or the mission
demanding thrust transient. We observed that when the igniter
turbulent intensity is relatively low the vehicle could liftoff early due
to the early flow choking of the rocket nozzle. We concluded that the
high pressurization-rate has structural implications at liftoff in
addition to transient burning effect. Therefore prudent selection of the
port geometry and the igniter, for meeting the mission requirements,
within the given envelop are meaningful objectives for any designer
for the smooth liftoff of solid propellant rockets.
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: 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: The Navier Stokes Equations (NSE) for an incompressible fluid of variable viscosity in the presence of an unknown external force in Von-Mises system x,\ are transformed, and some new exact solutions for a class of flows characterized by equation y f x a\b for an arbitrary state equation are determined, where f x is a function, \ the stream function, a z 0 and b are the arbitrary constants. In three, out of four cases, the function f x is arbitrary, and the solutions are the solutions of the flow equations for all the flows characterized by the equationy f x a\b. Streamline patterns for some forms of f x in unbounded and bounded regions are given.
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: 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.