Abstract: A new mechanism responsible for structural life
consumption due to resonant fatigue in turbine blades, or vanes, is
presented and explained. A rotating blade or vane in a gas turbine can
change its contour due to erosion and/or material build up, in any of
these instances, the surface pressure distribution occurring on the
suction and pressure sides of blades-vanes can suffer substantial
modification of their pressure and temperatures envelopes and flow
characteristics. Meanwhile, the relative rotation between the blade
and duct vane while the pressurized gas flows and the consequent
wake crossings, will induce a fluctuating thrust force or lift that will
excite the blade.
An actual totally used up set of vane-blade components in a HP
turbine power stage in a gas turbine is analyzed. The blade suffered
some material erosion mostly at the trailing edge provoking a
peculiar surface pressure envelope which evolved as the relative
position between the vane and the blade passed in front of each other.
Interestingly preliminary modal analysis for this eroded blade
indicates several natural frequencies within the aeromechanic power
spectrum, moreover, the highest frequency component is 94% of one
natural frequency indicating near resonant condition.
Independently of other simultaneously occurring fatigue cycles
(such as thermal, centrifugal stresses).
Abstract: The problem of incompressible steady flow simulation around an airfoil is discussed. For some simplest airfoils (circular, elliptical, Zhukovsky airfoils) the exact solution is known from complex analysis. It allows to compute the intensity of vortex layer which simulates the airfoil. Some modifications of the vortex element method are proposed and test computations are carried out. It-s shown that the these approaches are much more effective in comparison with the classical numerical scheme.
Abstract: In the present analysis an unsteady laminar
forced convection water boundary layer flow is considered.
The fluid properties such as viscosity and Prandtl number are
taken as variables such that those are inversely proportional to
temperature. By using quasi-linearization technique the nonlinear
coupled partial differential equations are linearized and
the numerical solutions are obtained by using implicit finite
difference scheme with the appropriate selection of step sizes.
Non-similar solutions have been obtained from the starting
point of the stream-wise coordinate to the point where skin
friction value vanishes. The effect non-uniform mass transfer
along the surface of the cylinder through slot is studied on the
skin friction and heat transfer coefficients.
Abstract: This is a study on numerical simulation of the convection-diffusion transport of a chemical species in steady flow through a small-diameter tube, which is lined with a very thin layer made up of retentive and absorptive materials. The species may be subject to a first-order kinetic reversible phase exchange with the wall material and irreversible absorption into the tube wall. Owing to the velocity shear across the tube section, the chemical species may spread out axially along the tube at a rate much larger than that given by the molecular diffusion; this process is known as dispersion. While the long-time dispersion behavior, well described by the Taylor model, has been extensively studied in the literature, the early development of the dispersion process is by contrast much less investigated. By early development, that means a span of time, after the release of the chemical into the flow, that is shorter than or comparable to the diffusion time scale across the tube section. To understand the early development of the dispersion, the governing equations along with the reactive boundary conditions are solved numerically using the Flux Corrected Transport Algorithm (FCTA). The computation has enabled us to investigate the combined effects on the early development of the dispersion coefficient due to the reversible and irreversible wall reactions. One of the results is shown that the dispersion coefficient may approach its steady-state limit in a short time under the following conditions: (i) a high value of Damkohler number (say Da ≥ 10); (ii) a small but non-zero value of absorption rate (say Γ* ≤ 0.5).
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: In this study, aeroelastic response and performance
analyses have been conducted for a 5MW-Class composite wind
turbine blade model. Advanced coupled numerical method based on
computational fluid dynamics (CFD) and computational flexible
multi-body dynamics (CFMBD) has been developed in order to
investigate aeroelastic responses and performance characteristics of
the rotating composite blade. Reynolds-Averaged Navier-Stokes
(RANS) equations with k-ω SST turbulence model were solved for
unsteady flow problems on the rotating turbine blade model. Also,
structural analyses considering rotating effect have been conducted
using the general nonlinear finite element method. A fully implicit
time marching scheme based on the Newmark direct integration
method is applied to solve the coupled aeroelastic governing equations
of the 3D turbine blade for fluid-structure interaction (FSI) problems.
Detailed dynamic responses and instantaneous velocity contour on the
blade surfaces which considering flow-separation effects were
presented to show the multi-physical phenomenon of the huge rotating
wind- turbine blade model.
Abstract: In this paper we study the transformation of Euler equations 1 , u u u Pf t (ρ ∂) + ⋅∇ = − ∇ + ∂ G G G G ∇⋅ = u 0, G where (ux, t) G G is the velocity of a fluid, P(x, t) G is the pressure of a fluid andρ (x, t) G is density. First of all, we rewrite the Euler equations in terms of new unknown functions. Then, we introduce new independent variables and transform it to a new curvilinear coordinate system. We obtain the Euler equations in the new dependent and independent variables. The governing equations into two subsystems, one is hyperbolic and another is elliptic.
Abstract: The present study investigates numerically the
phenomenon of vortex-shedding and its suppression in twodimensional
mixed convective flow past a square cylinder under the
joint influence of buoyancy and free-stream orientation with respect
to gravity. The numerical experiments have been conducted at a
fixed Reynolds number (Re) of 100 and Prandtl number (Pr) of 0.71,
while Richardson number (Ri) is varied from 0 to 1.6 and freestream
orientation, α, is kept in the range 0o≤ α ≤ 90o, with 0o
corresponding to an upward flow and 90o representing a cross-flow
scenario, respectively. The continuity, momentum and energy
equations, subject to Boussinesq approximation, are discretized using
a finite difference method and are solved by a semi-explicit pressure
correction scheme. The critical Richardson number, leading to the
suppression of the vortex-shedding (Ric), is estimated by using
Stuart-Landau theory at various free-stream orientations and the
neutral curve is obtained in the Ri-α plane. The neutral curve
exhibits an interesting non-monotonic behavior with Ric first
increasing with increasing values of α upto 45o and then decreasing
till 70o. Beyond 70o, the neutral curve again exhibits a sharp
increasing asymptotic trend with Ric approaching very large values
as α approaches 90o. The suppression of vortex shedding is not
observed at α = 90o (cross-flow). In the unsteady flow regime, the
Strouhal number (St) increases with the increase in Richardson
number.
Abstract: An experimental investigation was performed on pulp
liquid flow in straight ducts with a square cross section. Fully
developed steady flow was visualized and the fiber concentration was
obtained using a light-section method developed by the author et al.
The obtained results reveal quantitatively, in a definite form, the
distribution of the fiber concentration. From the results and
measurements of pressure loss, it is found that the flow characteristics
of pulp liquid in ducts can be classified into five patterns. The
relationships among the distributions of mean and fluctuation of fiber
concentration, the pressure loss and the flow velocity are discussed,
and then the features for each pattern are extracted. The degree of
nonuniformity of the fiber concentration, which is indicated by the
standard deviation of its distribution, is decreased from 0.3 to 0.05
with an increase in the velocity of the tested pulp liquid from 0.4 to
0.8%.
Abstract: Cardiovascular diseases, principally atherosclerosis, are responsible for 30% of world deaths. Atherosclerosis is due to the formation of plaque. The fatty plaque may be at risk of rupture, leading typically to stroke and heart attack. The plaque is usually associated with a high degree of lumen reduction, called a stenosis.It is increasingly recognized that the initiation and progression of disease and the occurrence of clinical events is a complex interplay between the local biomechanical environment and the local vascular biology. The aim of this study is to investigate the flow behavior through a stenosed artery. A physical experiment was performed using an artery model and blood analogue fluid. An axisymmetric model constructed consists of contraction and expansion region that follow a mathematical form of cosine function. A 30% diameter reduction was used in this study. The flow field was measured using particle image velocimetry (PIV). Spherical particles with 20μm diameter were seeded in a water-glycerol-NaCl mixture. Steady flow Reynolds numbers are 250. The area of interest is the region after the stenosis where the flow separation occurs. The velocity field was measured and the velocity gradient was investigated. There was high particle concentration in the recirculation zone. High velocity gradient formed immediately after the stenosis throat created a lift force that enhanced particle migration to the flow separation area.
Abstract: This paper presents a numerical investigation of the
unsteady flow around an American 19th century vertical-axis
windmill: the Stevens & Jolly rotor, patented on April 16, 1895. The
computational approach used is based on solving the complete
transient Reynolds-Averaged Navier-Stokes (t-RANS) equations: a
full campaign of numerical simulation has been performed using the
k-ω SST turbulence model. Flow field characteristics have been
investigated for several values of tip speed ratio and for a constant
unperturbed free-stream wind velocity of 6 m/s, enabling the study of
some unsteady flow phenomena in the rotor wake. Finally, the global
power generated from the windmill has been determined for each
simulated angular velocity, allowing the calculation of the rotor
power-curve.
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: 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.
Abstract: The rheological properties of light crude oil and its mixture with water were investigated experimentally. These rheological properties include steady flow behavior, yield stress, transient flow behavior, and viscoelastic behavior. A RheoStress RS600 rheometer was employed in all of the rheological examination tests. The light crude oil exhibits a Newtonian and for emulsion exhibits a non-Newtonian shear thinning behavior over the examined shear rate range of 0.1–120 s-1. In first time, a series of samples of crude oil from the Algerian Sahara has been tested and the results expressed in terms of τ=f(γ) have demonstrated their Newtonian character for the temperature included in [20°C, 70°C]. In second time and at T=20°C, the oil-water emulsions (30%, 50% and 70%) by volume of water), thermodynamically stable, have demonstrated a non-Newtonian rheological behavior that is to say, Herschel-Bulkley and Bingham types. For each type of crude oil-water emulsion, the rheological parameters are calculated by numerical treatment of results.