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: For many chemical and biological processes, the understanding of the mixing phenomenon and flow behavior in a stirred tank is of major importance. A three-dimensional numerical study was performed using the software Fluent, to study the flow field in a stirred tank with a Rushton turbine. In this work, we first studied the flow generated in the tank with a Rushton turbine. Then, we studied the effect of the variation of turbine’s submergence on the thermodynamic quantities defining the flow field. For that, four submergences were considered, while maintaining the same rotational speed (N =250rpm). This work intends to optimize the aeration performances of a Rushton turbine in a stirred tank.
Abstract: The incorporation of computational fluid dynamics in the design of modern hydraulic turbines appears to be necessary in order to improve their efficiency and cost-effectiveness beyond the traditional design practices. A numerical optimization methodology is developed and applied in the present work to a Turgo water turbine. The fluid is simulated by a Lagrangian mesh-free approach that can provide detailed information on the energy transfer and enhance the understanding of the complex, unsteady flow field, at very small computing cost. The runner blades are initially shaped according to hydrodynamics theory, and parameterized using Bezier polynomials and interpolation techniques. The use of a limited number of free design variables allows for various modifications of the standard blade shape, while stochastic optimization using evolutionary algorithms is implemented to find the best blade that maximizes the attainable hydraulic efficiency of the runner. The obtained optimal runner design achieves considerably higher efficiency than the standard one, and its numerically predicted performance is comparable to a real Turgo turbine, verifying the reliability and the prospects of the new methodology.
Abstract: Technology of thin film deposition is of interest in
many engineering fields, from electronic manufacturing to corrosion
protective coating. A typical deposition process, like that developed
at the University of Eindhoven, considers the deposition of a thin,
amorphous film of C:H or of Si:H on the substrate, using the
Expanding Thermal arc Plasma technique. In this paper a computing
procedure is proposed to simulate the flow field in a deposition
chamber similar to that at the University of Eindhoven and a
sensitivity analysis is carried out in terms of: precursor mass flow
rate, electrical power, supplied to the torch and fluid-dynamic
characteristics of the plasma jet, using different nozzles. To this
purpose a deposition chamber similar in shape, dimensions and
operating parameters to the above mentioned chamber is considered.
Furthermore, a method is proposed for a very preliminary evaluation
of the film thickness distribution on the substrate. The computing
procedure relies on two codes working in tandem; the output from
the first code is the input to the second one. The first code simulates
the flow field in the torch, where Argon is ionized according to the
Saha-s equation, and in the nozzle. The second code simulates the
flow field in the chamber. Due to high rarefaction level, this is a
(commercial) Direct Simulation Monte Carlo code. Gas is a mixture
of 21 chemical species and 24 chemical reactions from Argon plasma
and Acetylene are implemented in both codes. The effects of the
above mentioned operating parameters are evaluated and discussed
by 2-D maps and profiles of some important thermo-fluid-dynamic
parameters, as per Mach number, velocity and temperature. Intensity,
position and extension of the shock wave are evaluated and the
influence of the above mentioned test conditions on the film
thickness and uniformity of distribution are also evaluated.
Abstract: The flow field around a flat plate of infinite span has
been investigated for several values of the angle of attack. Numerical
predictions have been compared to experimental measurements, in
order to examine the effect of turbulence model and grid resolution
on the resultant aerodynamic forces acting on the plate. Also the
influence of the free-stream turbulence intensity, at the entrance of
the computational domain, has been investigated. A full campaign of
simulations has been conducted for three inclination angles (9°, 15°
and 30°), in order to obtain some practical guidelines to be used for
the simulation of the flow field around inclined plates and discs.
Abstract: Polymer melt compressibility and mold surface roughness, which are generally ignored during the filling stage of the conventional injection molding, may become increasingly significant in micro injection molding where the parts become smaller. By employing the 2.5D generalized Hele-Shaw model, we presented here the effects of polymer compressibility and mold surface roughness on mold-filling in a micro-thickness cavity. To elucidate the effects of surface roughness, numerical investigations were conducted using a cavity flat plate which has two halves with different surface roughness. This allows the comparison of flow field on two different halves under identical processing conditions but with different roughness. Results show that polymer compressibility and mold surface roughness have effects on mold filling in micro injection molding. There is in shrinkage reduction as the density is increased due to polymer melt compressibility during the filling stage.
Abstract: Non-uniform current distribution in polymer
electrolyte membrane fuel cells results in local over-heating,
accelerated ageing, and lower power output than expected. This
issue is very critical when fuel cell experiences water flooding. In
this work, the performance of a PEM fuel cell is investigated under
cathode flooding conditions. Two-dimensional partially flooded
GDL models based on the conservation laws and electrochemical
relations are proposed to study local current density distributions
along flow fields over a wide range of cell operating conditions.
The model results show a direct association between cathode inlet
humidity increases and that of average current density but the
system becomes more sensitive to flooding. The anode inlet
relative humidity shows a similar effect. Operating the cell at
higher temperatures would lead to higher average current densities
and the chance of system being flooded is reduced. In addition,
higher cathode stoichiometries prevent system flooding but the
average current density remains almost constant. The higher anode
stoichiometry leads to higher average current density and higher
sensitivity to cathode flooding.
Abstract: The two-phase flow field and the motion of the free
surface in an oscillating channel are simulated numerically to assess
the methodology for simulating nuclear reacotr thermal hydraulics
under seismic conditions. Two numerical methods are compared: one
is to model the oscillating channel directly using the moving grid of
the Arbitrary Lagrangian-Eulerian method, and the other is to simulate
the effect of channel motion using the oscillating acceleration acting
on the fluid in the stationary channel. The two-phase flow field in the
oscillating channel is simulated using the level set method in both
cases. The calculated results using the oscillating acceleration are
found to coinside with those using the moving grid, and the theoretical
back ground and the limitation of oscillating acceleration are discussed.
It is shown that the change in the interfacial area between liquid and
gas phases under seismic conditions is important for nuclear reactor
thermal hydraulics.
Abstract: The flow field over a three dimensional pole barn
characterized by a cylindrical roof has been numerically investigated.
Wind pressure and viscous loads acting on the agricultural building
have been analyzed for several incoming wind directions, so as to
evaluate the most critical load condition on the structure. A constant
wind velocity profile, based on the maximum reference wind speed in
the building site (peak gust speed worked out for 50 years return
period) and on the local roughness coefficient, has been simulated.
In order to contemplate also the hazard due to potential air
wedging between the stored hay and the lower part of the ceiling, the
effect of a partial filling of the barn has been investigated.
The distribution of wind-induced loads on the structure have been
determined, allowing a numerical quantification of the effect of wind
direction on the induced stresses acting on a hemicylindrical roof.
Abstract: The aerodynamic noise radiation from a side view mirror (SVM) in the high-speed airflow is calculated by the combination of unsteady incompressible fluid flow analysis and acoustic analysis. The transient flow past the generic SVM is simulated with variable turbulence model, namely DES Detached Eddy Simulation and LES (Large Eddy Simulation). Detailed velocity vectors and contour plots of the time-varying velocity and pressure fields are presented along cut planes in the flow-field. Mean and transient pressure are also monitored at several points in the flow field and compared to corresponding experimentally data published in literature. The acoustic predictions made using the Ffowcs-Williams-Hawkins acoustic analogy (FW-H) and the boundary element (BEM).
Abstract: Future space vehicles will require the use of non-toxic, cryogenic propellants, because of the performance advantages over the toxic hypergolic propellants and also because of the environmental and handling concerns. A prototypical capillary flow liquid acquisition device (LAD) for cryogenic propellants was fabricated with a mesh screen, covering a rectangular flow channel with a cylindrical outlet tube, and was tested with liquid oxygen (LOX). In order to better understand the performance in various gravity environments and orientations with different submersion depths of the LAD, a series of computational fluid dynamics (CFD) simulations of LOX flow through the LAD screen channel, including horizontally and vertically submersions of the LAD channel assembly at normal gravity environment was conducted. Gravity effects on the flow field in LAD channel are inspected and analyzed through comparing the simulations.
Abstract: The flow field and the motion of the free surface in an
oscillating container are simulated numerically to assess the numerical
approach for studying two-phase flows under oscillating conditions.
Two numerical methods are compared: one is to model the oscillating
container directly using the moving grid of the ALE method, and the
other is to simulate the effect of container motion using the oscillating
body force acting on the fluid in the stationary container. The
two-phase flow field in the container is simulated using the level set
method in both cases. It is found that the calculated results by the body
force method coinsides with those by the moving grid method and the
sloshing behavior is predicted well by both the methods. Theoretical
back ground and limitation of the body force method are discussed,
and the effects of oscillation amplitude and frequency are shown.
Abstract: This paper presents a numerical analysis of the
performance of a three-bladed Darrieus vertical-axis wind turbine
based on the DU91-W2-250 airfoil. A complete campaign of 2-D
simulations, performed for several values of tip speed ratio and based
on RANS unsteady calculations, has been performed to obtain the
rotor torque and power curves. Rotor performances have been
compared with the results of a previous work based on the use of the
NACA 0021 airfoil. Both the power coefficient and the torque
coefficient have been determined as a function of the tip speed ratio.
The flow field around rotor blades has also been analyzed. As a final
result, the performance of the DU airfoil based rotor appears to be
lower than the one based on the NACA 0021 blade section. This
behavior could be due to the higher stall characteristics of the NACA
profile, being the separation zone at the trailing edge more extended
for the DU airfoil.
Abstract: In this research a comparison between k-epsilon and
LES model for a deoiling hydrocyclone is conducted. Flow field of
hydrocyclone is obtained by three-dimensional simulations with
OpenFOAM code. Potential of prediction for both methods of this
complex swirl flow is discussed. Large eddy simulation method
results have more similarity to experiment and its results are
presented in figures from different hydrocyclone cross sections.
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: This paper investigates the nature of the development
of two-dimensional laminar flow of an incompressible fluid at the
reversed stagnation-point. ". In this study, we revisit the problem
of reversed stagnation-point flow over a flat plate. Proudman and
Johnson (1962) first studied the flow and obtained an asymptotic
solution by neglecting the viscous terms. This is no true in neglecting
the viscous terms within the total flow field. In particular it is pointed
out that for a plate impulsively accelerated from rest to a constant
velocity V0 that a similarity solution to the self-similar ODE is
obtained which is noteworthy completely analytical.
Abstract: A new conserving approach in the context of Immersed Boundary Method (IBM) is presented to simulate one dimensional, incompressible flow in a moving boundary problem. The method employs control volume scheme to simulate the flow field. The concept of ghost node is used at the boundaries to conserve the mass and momentum equations. The Present method implements the conservation laws in all cells including boundary control volumes. Application of the method is studied in a test case with moving boundary. Comparison between the results of this new method and a sharp interface (Image Point Method) IBM algorithm shows a well distinguished improvement in both pressure and velocity fields of the present method. Fluctuations in pressure field are fully resolved in this proposed method. This approach expands the IBM capability to simulate flow field for variety of problems by implementing conservation laws in a fully Cartesian grid compared to other conserving methods.
Abstract: In the present work steady inviscid hypersonic flows
are calculated by approximate Method. Maslens' inverse method is the chosen approximate method. For the inverse problem, parabolic
shock shape is chosen for the two-dimensional flow, and the body shape and flow field are calculated using Maslen's method. For the axisymmetric inverse problem paraboloidal shock is chosen and the
surface distribution of pressure is obtained.
Abstract: This paper looks into detailed investigation of
thermal-hydraulic characteristics of the flow field in a fuel rod
model, especially near the spacer. The area investigate represents a
source of information on the velocity flow field, vortex, and on the
amount of heat transfer into the coolant all of which are critical for
the design and improvement of the fuel rod in nuclear power plants.
The flow field investigation uses three-dimensional Computational
Fluid Dynamics (CFD) with the Reynolds stresses turbulence model
(RSM). The fuel rod model incorporates a vertical annular channel
where three different shapes of spacers are used; each spacer shape is
addressed individually. These spacers are mutually compared in
consideration of heat transfer capabilities between the coolant and
the fuel rod model. The results are complemented with the calculated
heat transfer coefficient in the location of the spacer and along the
stainless-steel pipe.
Abstract: We present a Large-Eddy simulation of a vortex cell
with circular shaped. The results show that the flow field can be sub
divided into four important zones, the shear layer above the cavity,
the stagnation zone, the vortex core in the cavity and the boundary
layer along the wall of the cavity. It is shown that the vortex core
consits of solid body rotation without much turbulence activity. The
vortex is mainly driven by high energy packets that are driven into the
cavity from the stagnation point region and by entrainment of fluid
from the cavity into the shear layer. The physics in the boundary
layer along the cavity-s wall seems to be far from that of a canonical
boundary layer which might be a crucial point for modelling this
flow.