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: 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.
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.
Abstract: The operation performance of a valveless micro-pump
is strongly dependent on the shape of connected nozzle/diffuser and
Reynolds number. The aims of present work are to compare the
performance curves of micropump with the original straight
nozzle/diffuser and contoured nozzle/diffuser under different back
pressure conditions. The tested valveless micropumps are assembled
of five pieces of patterned PMMA plates with hot-embracing
technique. The structures of central chamber, the inlet/outlet
reservoirs and the connected nozzle/diffuser are fabricated with laser
cutting machine. The micropump is actuated with circular-type PZT
film embraced on the bottom of central chamber. The deformation of
PZT membrane with various input voltages is measured with a
displacement laser probe. A simple testing facility is also constructed
to evaluate the performance curves for comparison.
In order to observe the evaluation of low Reynolds number
multiple vortex flow patterns within the micropump during suction
and pumping modes, the unsteady, incompressible laminar
three-dimensional Reynolds-averaged Navier-Stokes equations are
solved. The working fluid is DI water with constant thermo-physical
properties. The oscillating behavior of PZT film is modeled with the
moving boundary wall in way of UDF program. With the dynamic
mesh method, the instants pressure and velocity fields are obtained
and discussed.Results indicated that the volume flow rate is not
monotony increased with the oscillating frequency of PZT film,
regardless of the shapes of nozzle/diffuser. The present micropump
can generate the maximum volume flow rate of 13.53 ml/min when
the operation frequency is 64Hz and the input voltage is 140 volts.
The micropump with contoured nozzle/diffuser can provide 7ml/min
flow rate even when the back pressure is up to 400 mm-H2O. CFD
results revealed that the flow central chamber was occupied with
multiple pairs of counter-rotating vortices during suction and
pumping modes. The net volume flow rate over a complete
oscillating periodic of PZT
Abstract: This paper analytically investigates the 3D flow
pattern at the confluences of two rectangular channels having 900
angles using Navier-Stokes equations based on Reynolds Stress
Turbulence Model (RSM). The equations are solved by the Finite-
Volume Method (FVM) and the flow is analyzed in terms of steadystate
(single-phased) conditions. The Shumate experimental findings
were used to test the validity of data. Comparison of the simulation
model with the experimental ones indicated a close proximity
between the flow patterns of the two sets. Effects of the discharge
ratio on separation zone dimensions created in the main-channel
downstream of the confluence indicated an inverse relation, where a
decrease in discharge ratio, will entail an increase in the length and
width of the separation zone. The study also found the model as a
powerful analytical tool in the feasibility study of hydraulic
engineering projects.
Abstract: The three-dimensional incompressible flow past a
rectangular open cavity is investigated, where the aspect ratio of the
cavity is considered as 4. The principle objective is to use large-eddy
simulation to resolve and control the large-scale structures, which are
largely responsible for flow oscillations in a cavity. The flow past an
open cavity is very common in aerospace applications and can be a
cause of acoustic source due to hydrodynamic instability of the shear
layer and its interactions with the downstream edge. The unsteady
Navier-stokes equations have been solved on a staggered mesh using
a symmetry-preserving central difference scheme. Synthetic jet has
been used as an active control to suppress the cavity oscillations in
wake mode for a Reynolds number of ReD = 3360. The effect of
synthetic jet has been studied by varying the jet amplitude and
frequency, which is placed at the upstream wall of the cavity. The
study indicates that there exits a frequency band, which is larger than
a critical value, is effective in attenuating cavity oscillations when
blowing ratio is more than 1.0.
Abstract: Incompressible Navier-Stokes equations are reviewed
in this work. Three-dimensional Navier-Stokes equations are solved
analytically. The Mathematical derivation shows that the solutions
for the zero and constant pressure gradients are similar. Descriptions
of the proposed formulation and validation against two laminar
experiments and three different turbulent flow cases are reported in
this paper. Even though, the analytical solution is derived for nonreacting
flows, it could reproduce trends for cases including
combustion.
Abstract: This research proposes an algorithm for the simulation
of time-periodic unsteady problems via the solution unsteady Euler
and Navier-Stokes equations. This algorithm which is called Time
Spectral method uses a Fourier representation in time and hence
solve for the periodic state directly without resolving transients
(which consume most of the resources in a time-accurate scheme).
Mathematical tools used here are discrete Fourier transformations. It
has shown tremendous potential for reducing the computational cost
compared to conventional time-accurate methods, by enforcing
periodicity and using Fourier representation in time, leading to
spectral accuracy. The accuracy and efficiency of this technique is
verified by Euler and Navier-Stokes calculations for pitching airfoils.
Because of flow turbulence nature, Baldwin-Lomax turbulence
model has been used at viscous flow analysis. The results presented
by the Time Spectral method are compared with experimental data. It
has shown tremendous potential for reducing the computational cost
compared to the conventional time-accurate methods, by enforcing
periodicity and using Fourier representation in time, leading to
spectral accuracy, because results verify the small number of time
intervals per pitching cycle required to capture the flow physics.
Abstract: A new numerical method for solving the twodimensional,
steady, incompressible, viscous flow equations on a
Curvilinear staggered grid is presented in this paper. The proposed
methodology is finite difference based, but essentially takes
advantage of the best features of two well-established numerical
formulations, the finite difference and finite volume methods. Some
weaknesses of the finite difference approach are removed by
exploiting the strengths of the finite volume method. In particular,
the issue of velocity-pressure coupling is dealt with in the proposed
finite difference formulation by developing a pressure correction
equation in a manner similar to the SIMPLE approach commonly
used in finite volume formulations. However, since this is purely a
finite difference formulation, numerical approximation of fluxes is
not required. Results obtained from the present method are based on
the first-order upwind scheme for the convective terms, but the
methodology can easily be modified to accommodate higher order
differencing schemes.
Abstract: A 3D simulation study for an incompressible
slip flow around a spherical aerosol particle was performed.
The full Navier-Stokes equations were solved and the velocity
jump at the gas-particle interface was treated numerically by
imposition of the slip boundary condition. Analytical solution
to the Stokesian slip flow past a spherical particle was used as
a benchmark for code verification, and excellent agreement
was achieved. The Simulation results showed that in addition
to the Knudsen number, the Reynolds number affects the slip
correction factor. Thus, the Cunningham-based slip corrections
must be augmented by the inclusion of the effect of
Reynolds number for application to Lagrangian tracking of
fine particles. A new expression for the slip correction factor
as a function of both Knudsen number and Reynolds number
was developed.
Abstract: This is the second part of the paper. It, aside from the
core subroutine test reported previously, focuses on the simulation of
turbulence governed by the full STF Navier-Stokes equations on a
large scale. Law of the wall is found plausible in this study as a model
of the boundary layer dynamics. Model validations proceed to
include velocity profiles of a stationary turbulent Couette flow, pure
sloshing flow simulations, and the identification of water-surface
inclination due to fluid accelerations. Errors resulting from the
irrotational and hydrostatic assumptions are explored when studying
a wind-driven water circulation with no shakings. Illustrative
examples show that this numerical strategy works for the simulation
of sloshing-shear mixed flow in a 3-D rigid rectangular base tank.
Abstract: In this paper, the differential quadrature method is applied to simulate natural convection in an inclined cubic cavity using velocity-vorticity formulation. The numerical capability of the present algorithm is demonstrated by application to natural convection in an inclined cubic cavity. The velocity Poisson equations, the vorticity transport equations and the energy equation are all solved as a coupled system of equations for the seven field variables consisting of three velocities, three vorticities and temperature. The coupled equations are simultaneously solved by imposing the vorticity definition at boundary without requiring the explicit specification of the vorticity boundary conditions. Test results obtained for an inclined cubic cavity with different angle of inclinations for Rayleigh number equal to 103, 104, 105 and 106 indicate that the present coupled solution algorithm could predict the benchmark results for temperature and flow fields. Thus, it is convinced that the present formulation is capable of solving coupled Navier-Stokes equations effectively and accurately.
Abstract: A parallel computational fluid dynamics code has been
developed for the study of aerodynamic heating problem in hypersonic
flows. The code employs the 3D Navier-Stokes equations as the basic
governing equations to simulate the laminar hypersonic flow. The cell
centered finite volume method based on structured grid is applied for
spatial discretization. The AUSMPW+ scheme is used for the inviscid
fluxes, and the MUSCL approach is used for higher order spatial
accuracy. The implicit LU-SGS scheme is applied for time integration
to accelerate the convergence of computations in steady flows. A
parallel programming method based on MPI is employed to shorten
the computing time. The validity of the code is demonstrated by
comparing the numerical calculation result with the experimental data
of a hypersonic flow field around a blunt body.
Abstract: A general purpose viscous flow solver Ansys CFX
was used to solve the unsteady three-dimensional (3D) Reynolds
Averaged Navier-Stokes Equation (RANSE) for simulating a 3D
numerical viscous wave tank. A flap-type wave generator was
incorporated in the computational domain to generate the desired
incident waves. Authors have made effort to study the physical
behaviors of Flap type wave maker with governing parameters.
Dependency of the water fill depth, Time period of oscillations and
amplitude of oscillations of flap were studied. Effort has been made
to establish relations between parameters. A validation study was
also carried out against CFD methodology with wave maker theory.
It has been observed that CFD results are in good agreement with
theoretical results. Beaches of different slopes were introduced to
damp the wave, so that it should not cause any reflection from
boundary. As a conclusion this methodology can simulate the
experimental wave-maker for regular wave generation for different
wave length and amplitudes.
Abstract: In this study, the numerical solution of unsteady flow
between two concentric rotating spheres with suction and blowing at
their boundaries is presented. The spheres are rotating about a
common axis of rotation while their angular velocities are constant.
The Navier-Stokes equations are solved by employing the finite
difference method and implicit scheme. The resulting flow patterns
are presented for various values of the flow parameters including
rotational Reynolds number Re , and a blowing/suction Reynolds
number Rew . Viscous torques at the inner and the outer spheres are
calculated, too. It is seen that increasing the amount of suction and
blowing decrease the size of eddies generated in the annulus.