Abstract: Flow through micro and mini channels requires relatively
high driving pressure due to the large fluid pressure drop
through these channels. Consequently the forces acting on the walls of
the channel due to the fluid pressure are also large. Due to these forces
there are displacement fields set up in the solid substrate containing
the channels. If the movement of the substrate is constrained at some
points, then stress fields are established in the substrate. On the other
hand, if the deformation of the channel shape is sufficiently large
then its effect on the fluid flow is important to be calculated. Such
coupled fluid-solid systems form a class of problems known as fluidstructure
interactions. In the present work a co-located finite volume
discretization procedure on unstructured meshes is described for
solving fluid-structure interaction type of problems. A linear elastic
solid is assumed for which the effect of the channel deformation
on the flow is neglected. Thus the governing equations for the
fluid and the solid are decoupled and are solved separately. The
procedure is validated by solving two benchmark problems, one from
fluid mechanics and another from solid mechanics. A fluid-structure
interaction problem of flow through a U-shaped channel embedded
in a plate is solved.
Abstract: Calcium [Ca2+] is an important second messenger
which plays an important role in signal transduction. There are
several parameters that affect its concentration profile like buffer
source etc. The effect of stationary immobile buffer on Ca2+
concentration has been incorporated which is a very important
parameter needed to be taken into account in order to make the
model more realistic. Interdependence of all the important parameters
like diffusion coefficient and influx over [Ca2+] profile has been
studied. Model is developed in the form of advection diffusion
equation together with buffer concentration. A program has been
developed using finite volume method for the entire problem and
simulated on an AMD-Turion 32-bit machine to compute the
numerical results.
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: In this paper, the melting of a semi-infinite body as a
result of a moving laser beam has been studied. Because the Fourier
heat transfer equation at short times and large dimensions does not
have sufficient accuracy; a non-Fourier form of heat transfer
equation has been used. Due to the fact that the beam is moving in x
direction, the temperature distribution and the melting pool shape are
not asymmetric. As a result, the problem is a transient threedimensional
problem. Therefore, thermophysical properties such as
heat conductivity coefficient, density and heat capacity are functions
of temperature and material states. The enthalpy technique, used for
the solution of phase change problems, has been used in an explicit
finite volume form for the hyperbolic heat transfer equation. This
technique has been used to calculate the transient temperature
distribution in the semi-infinite body and the growth rate of the melt
pool. In order to validate the numerical results, comparisons were
made with experimental data. Finally, the results of this paper were
compared with similar problem that has used the Fourier theory. The
comparison shows the influence of infinite speed of heat propagation
in Fourier theory on the temperature distribution and the melt pool
size.
Abstract: Fluid flow in cylinders of elliptic cross-section was investigated. Fluid used is Liquefied petroleum gas (LPG). LPG found in Nigeria contains majorly butane with percentages of propane. Commercial available code FLUENT which uses finite volume method was used to solve fluid flow governing equations. There has been little attention paid to fluid flow in cylindrical elliptic pipes. The present work aims to predict the LPG gas flow in cylindrical pipes of elliptic cross-section. Results of flow parameters of velocity and pressure distributions are presented. Results show that the pressure drop in elliptic pipes is higher than circular pipe of the same cross-sectional area. This is an important result as the pressure drop is related to the pump power needed to drive the flow. Results show that the velocity increases towards centre of the pipe as the flow moves downstream, and also increases towards the outlet of the pipe.
Abstract: In this paper, a two-dimensional (2D) numerical
model for the tidal currents simulation in Persian Gulf is presented.
The model is based on the depth averaged equations of shallow water
which consider hydrostatic pressure distribution. The continuity
equation and two momentum equations including the effects of bed
friction, the Coriolis effects and wind stress have been solved. To
integrate the 2D equations, the Alternative Direction Implicit (ADI)
technique has been used. The base of equations discritization was
finite volume method applied on rectangular mesh. To evaluate the
model validation, a dam break case study including analytical
solution is selected and the comparison is done. After that, the
capability of the model in simulation of tidal current in a real field is
represented by modeling the current behavior in Persian Gulf. The
tidal fluctuations in Hormuz Strait have caused the tidal currents in
the area of study. Therefore, the water surface oscillations data at
Hengam Island on Hormoz Strait are used as the model input data.
The check point of the model is measured water surface elevations at
Assaluye port. The comparison between the results and the
acceptable agreement of them showed the model ability for modeling
marine hydrodynamic.
Abstract: In this work, axisymetric CFD simulation of fixed bed
GTL reactor has been conducted, using computational fluid dynamics
(CFD). In fixed bed CFD modeling, when N (tube-to-particle
diameter ratio) has a large value, it is common to consider the packed
bed as a porous media. Synthesis gas (a mixture of predominantly
carbon monoxide and hydrogen) was fed to the reactor. The reactor
length was 20 cm, divided to three sections. The porous zone was in
the middle section of the reactor. The model equations were solved
employing finite volume method. The effects of particle diameter,
bed voidage, fluid velocity and bed length on pressure drop have
been investigated. Simulation results showed these parameters could
have remarkable impacts on the reactor pressure drop.
Abstract: Mathematical and computational modeling of calcium
signalling in nerve cells has produced considerable insights into how
the cells contracts with other cells under the variation of biophysical
and physiological parameters. The modeling of calcium signaling in
astrocytes has become more sophisticated. The modeling effort has
provided insight to understand the cell contraction. Main objective
of this work is to study the effect of voltage gated (Operated)
calcium channel (VOC) on calcium profile in the form of advection
diffusion equation. A mathematical model is developed in the form
of advection diffusion equation for the calcium profile. The model
incorporates the important physiological parameter like diffusion
coefficient etc. Appropriate boundary conditions have been framed.
Finite volume method is employed to solve the problem. A program
has been developed using in MATLAB 7.5 for the entire problem
and simulated on an AMD-Turion 32-bite machine to compute the
numerical results.
Abstract: The Spalart and Allmaras turbulence model has been
implemented in a numerical code to study the compressible turbulent
flows, which the system of governing equations is solved with a
finite volume approach using a structured grid. The AUSM+ scheme
is used to calculate the inviscid fluxes. Different benchmark
problems have been computed to validate the implementation and
numerical results are shown. A special Attention is paid to wall jet
applications. In this study, the jet is submitted to various wall
boundary conditions (adiabatic or uniform heat flux) in forced
convection regime and both two-dimensional and axisymmetric wall
jets are considered. The comparison between the numerical results
and experimental data has given the validity of this turbulence model
to study the turbulent wall jets especially in engineering applications.
Abstract: Turbulent forced convection flow in a 2-dimensional channel over periodic grooves is numerically investigated. Finite volume method is used to study the effect of turbulence model. The range of Reynolds number varied from 10000 to 30000 for the ribheight to channel-height ratio (B/H) of 2. The downstream wall is heated by a uniform heat flux while the upstream wall is insulated. The investigation is analyzed with different types of nanoparticles such as SiO2, Al2O3, and ZnO, with water as a base fluid are used. The volume fraction is varied from 1% to 4% and the nanoparticle diameter is utilized between 20nm to 50nm. The results revealed 114% heat transfer enhancement compared to the water in a grooved channel by using SiO2 nanoparticle with volume fraction and nanoparticle diameter of 4% and 20nm respectively.
Abstract: The aim of this work is to analyze a viscous flow in
the axisymmetric nozzle 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 supersonic flow parameters at the exit of convergingdiverging
nozzle. The numerical technique uses the Flux Vector
Splitting method of Van Leer. Here, adequate time stepping
parameter, along with CFL coefficient and mesh size level is selected
to ensure numerical convergence. The effect of the boundary layer
thickness is significant at the exit of the nozzle. The best solution is
obtained with using a very fine grid, especially near the wall, where
we have a strong variation of velocity, temperature and shear stress.
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: Flow around a flat tube is studied numerically. Reynolds number is defined base on equivalent circular tube and it is varied in range of 100 to 300. Equations are solved by using finite volume method and results are presented in form of drag and lift coefficient. Results show that drag coefficient of flat tube is up to 66% lower than circular tube with equivalent diameter. In addition, by increasing l/D from 1 to 2, the drag coefficient of flat tube is decreased about 14-27%.
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: In mechanical and environmental engineering, mixed
convection is a frequently encountered thermal fluid phenomenon
which exists in atmospheric environment, urban canopy flows, ocean
currents, gas turbines, heat exchangers, and computer chip cooling
systems etc... . This paper deals with a numerical investigation of
mixed convection in a vertical heated channel. This flow results from
the mixing of the up-going fluid along walls of the channel with the
one issued from a flat nozzle located in its entry section. The fluiddynamic
and heat-transfer characteristics of vented vertical channels
are investigated for constant heat-flux boundary conditions, a
Rayleigh number equal to 2.57 1010, for two jet Reynolds number
Re=3 103 and 2104 and the aspect ratio in the 8-20 range. The system
of governing equations is solved with a finite volumes method and an
implicit scheme. The obtained results show that the turbulence and
the jet-wall interaction activate the heat transfer, as does the drive of
ambient air by the jet. For low Reynolds number Re=3 103, the
increase of the aspect Ratio enhances the heat transfer of about 3%,
however; for Re=2 104, the heat transfer enhancement is of about
12%. The numerical velocity, pressure and temperature fields are
post-processed to compute the quantities of engineering interest such
as the induced mass flow rate, and average Nusselt number, in terms
of Rayleigh, Reynolds numbers and dimensionless geometric
parameters are presented.
Abstract: Hypersonic flows around spatial vehicles during their
reentry phase in planetary atmospheres are characterized by intense
aerothermal phenomena. The aim of this work is to analyze high
temperature flows around an axisymmetric blunt body taking into
account chemical and vibrational non-equilibrium for air mixture
species. For this purpose, a finite volume methodology is employed
to determine the supersonic flow parameters around the axisymmetric
blunt body, especially at the stagnation point and along the wall of
spacecraft for several altitudes. This allows the capture shock wave
before a blunt body placed in supersonic free stream. The numerical
technique uses the Flux Vector Splitting method of Van Leer. Here,
adequate time stepping parameter, along with CFL coefficient and
mesh size level are selected to ensure numerical convergence, sought
with an order of 10-8
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: Predictions of flow and heat transfer characteristics and shape optimization in internally finned circular tubes have been performed on three-dimensional periodically fully developed turbulent flow and thermal fields. For a trapezoidal fin profile, the effects of fin height h, upper fin widths d1, lower fin widths d2, and helix angle of fin ? on transport phenomena are investigated for the condition of fin number of N = 30. The CFD and mathematical optimization technique are coupled in order to optimize the shape of internally finned tube. The optimal solutions of the design variables (i.e., upper and lower fin widths, fin height and helix angle) are numerically obtained by minimizing the pressure loss and maximizing the heat transfer rate, simultaneously, for the limiting conditions of d1 = 0.5~1.5 mm, d2 = 0.5~1.5 mm, h= 0.5~1.5mm, ? = 10~30 degrees. The fully developed flow and thermal fields are predicted using the finite volume method and the optimization is carried out by means of the multi-objective genetic algorithm that is widely used in the constrained nonlinear optimization problem.
Abstract: We numerically study the three-dimensional
magnetohydrodynamics (MHD) stability of oscillatory natural
convection flow in a rectangular cavity, with free top surface, filled
with a liquid metal, having an aspect ratio equal to A=L/H=5, and
subjected to a transversal temperature gradient and a uniform
magnetic field oriented in x and z directions. The finite volume
method was used in order to solve the equations of continuity,
momentum, energy, and potential. The stability diagram obtained in
this study highlights the dependence of the critical value of the
Grashof number Grcrit , with the increase of the Hartmann number
Ha for two orientations of the magnetic field. This study confirms
the possibility of stabilization of a liquid metal flow in natural
convection by application of a magnetic field and shows that the
flow stability is more important when the direction of magnetic field
is longitudinal than when the direction is transversal.
Abstract: It is essential to have a uniform and calm flow field
for a settling tank to have high performance. In general, the
recirculation zones always occurred in sedimentation tanks. The
presence of these regions may have different effects. The nonuniformity
of the velocity field, the short-circuiting at the surface and
the motion of the jet at the bed of the tank that occurs because of the
recirculation in the sedimentation layer, are affected by the geometry
of the tank. There are some ways to decrease the size of these dead
zones, which would increase the performance. One of the ways is to
use a suitable baffle configuration. In this study, the presence of
baffle with different position has been investigated by a finite volume
method, with VOF (Volume of Fluid) model. Besides, the k-ε
turbulence model is used in the numerical calculations. The results
indicate that the best position of the baffle is obtained when the
volume of the recirculation region is minimized or is divided to
smaller part and the flow field trend to be uniform in the settling
zone.
Abstract: Solid fuel transient burning behavior under oxidizer
gas flow is numerically investigated. It is done using analysis of the
regression rate responses to the imposed sudden and oscillatory
variation at inflow properties. The conjugate problem is considered
by simultaneous solution of flow and solid phase governing
equations to compute the fuel regression rate. The advection
upstream splitting method is used as flow computational scheme in
finite volume method. The ignition phase is completely simulated to
obtain the exact initial condition for response analysis. The results
show that the transient burning effects which lead to the combustion
instabilities and intermittent extinctions could be observed in solid
fuels as the solid propellants.