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: 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 potential flow model is used to study the unsteady
flow past two airfoils in configuration, each of which is suddenly set
into motion. The airfoil bound vortices are modeled using lumped
vortex elements and the wake behind the airfoil is modeled by discrete
vortices. This consists of solving a steady state flow problem at each
time-step where unsteadiness is incorporated through the “zero normal
flow on a solid surface" boundary condition at every time instant.
Additionally, along with the “zero normal flow on a solid surface"
boundary condition Kelvin-s condition is used to compute the strength
of the latest wake vortex shed from the trailing edge of the airfoil.
Location of the wake vortices is updated at each time-step to get the
wake shape at each time instant. Results are presented to show the
effect of airfoil-airfoil interaction and airfoil-wake interaction on the
aerodynamic characteristics of each airfoil.
Abstract: Unsteady natural convection and heat transfer in a square cavity partially filled with porous media using a thermal
non-equilibrium model is studied in this paper. The left vertical wall is
maintained at a constant hot temperature Th and the right vertical wall
is maintained at a constant cold temperature Tc, while the horizontal
walls are adiabatic. The governing equations are obtained by applying
the Darcy model and Boussinesq approximation. COMSOL’s finite
element method is used to solve the non-dimensional governing
equations together with specified boundary conditions. The governing
parameters of this study are the Rayleigh number (Ra = 10^5, and Ra = 10^6 ), Darcy namber (Da = 10^−2, and Da = 10^−3),
the modified thermal conductivity ratio (10^−1 ≤ γ ≤ 10^4), the inter-phase heat transfer coefficien (10^−1 ≤ H ≤ 10^3) and the
time dependent (0.001 ≤ τ ≤ 0.2). The results presented for
values of the governing parameters in terms of streamlines in both
fluid/porous-layer, isotherms of fluid in fluid/porous-layer, isotherms
of solid in porous layer, and average Nusselt number.
Abstract: The Navier–Stokes equations for unsteady, incompressible, viscous fluids in the axisymmetric coordinate system are solved using a control volume method. The volume-of-fluid (VOF) technique is used to track the free-surface of the liquid. Model predictions are in good agreement with experimental measurements. It is found that the dynamic processes after impact are sensitive to the initial droplet velocity and the liquid pool depth. The time evolution of the crown height and diameter are obtained by numerical simulation. The critical We number for splashing (Wecr) is studied for Oh (Ohnesorge) numbers in the range of 0.01~0.1; the results compares well with those of the experiments.
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: The paper reports on the results of experimental and
numerical study of nonstationary swirling flow in an isothermal
model of vortex burner. It has been identified that main source of the
instability is related to a precessing vortex core (PVC) phenomenon.
The PVC induced flow pulsation characteristics such as precession
frequency and its variation as a function of flowrate and swirl number
have been explored making use of acoustic probes. Additionally
pressure transducers were used to measure the pressure drops on the
working chamber and across the vortex flow. The experiments have
been included also the mean velocity measurements making use of a
laser-Doppler anemometry. The features of instantaneous flowfield
generated by the PVC were analyzed employing a commercial CFD
code (Star-CCM+) based on Detached Eddy Simulation (DES)
approach. Validity of the numerical code has been checked by
comparison calculated flowfield data with the obtained experimental
results. It has been confirmed particularly that the CFD code applied
correctly reproduces the flow features.
Abstract: Generally flow behavior in centrifugal fan is observed
to be in a state of instability with flow separation zones on suction
surface as well as near the front shroud. Overall performance of the
diffusion process in a centrifugal fan could be enhanced by
judiciously introducing the boundary layer suction slots. With easy
accessibility of CFD as an analytical tool, an extensive numerical
whole field analysis of the effect of boundary layer suction slots in
discrete regions of suspected separation points is possible. This paper
attempts to explore the effect of boundary layer suction slots
corresponding to various geometrical locations on the impeller with
converging configurations for the slots. The analysis shows that the
converging suction slots located on the impeller blade about 25%
from the trailing edge, significantly improves the static pressure
recovery across the fan. Also it is found that Slots provided at a
radial distance of about 12% from the leading and trailing edges
marginally improve the static pressure recovery across the fan.
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: 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 paper we propose, a Lagrangian method to solve unsteady gas equation which is a nonlinear ordinary differential equation on semi-infnite interval. This approach is based on Modified generalized Laguerre functions. This method reduces the solution of this problem to the solution of a system of algebraic equations. We also compare this work with some other numerical results. The findings show that the present solution is highly accurate.
Abstract: This paper analyses the unsteady, two-dimensional
stagnation point flow of an incompressible viscous fluid over a flat
sheet when the flow is started impulsively from rest and at the same
time, the sheet is suddenly stretched in its own plane with a velocity
proportional to the distance from the stagnation point. The partial
differential equations governing the laminar boundary layer forced
convection flow are non-dimensionalised using semi-similar
transformations and then solved numerically using an implicit finitedifference
scheme known as the Keller-box method. Results
pertaining to the flow and heat transfer characteristics are computed
for all dimensionless time, uniformly valid in the whole spatial region
without any numerical difficulties. Analytical solutions are also
obtained for both small and large times, respectively representing the
initial unsteady and final steady state flow and heat transfer.
Numerical results indicate that the velocity ratio parameter is found
to have a significant effect on skin friction and heat transfer rate at
the surface. Furthermore, it is exposed that there is a smooth
transition from the initial unsteady state flow (small time solution) to
the final steady state (large time solution).
Abstract: In the present study, a numerical analysis is carried
out to investigate unsteady MHD (magneto-hydrodynamic) flow and
heat transfer of a non-Newtonian second grade viscoelastic fluid
over an oscillatory stretching sheet. The flow is induced due to an
infinite elastic sheet which is stretched oscillatory (back and forth) in
its own plane. Effect of viscous dissipation and joule heating are
taken into account. The non-linear differential equations governing
the problem are transformed into system of non-dimensional
differential equations using similarity transformations. A newly
developed meshfree numerical technique Element free Galerkin
method (EFGM) is employed to solve the coupled non linear
differential equations. The results illustrating the effect of various
parameters like viscoelastic parameter, Hartman number, relative
frequency amplitude of the oscillatory sheet to the stretching rate and
Eckert number on velocity and temperature field are reported in
terms of graphs and tables. The present model finds its application in
polymer extrusion, drawing of plastic films and wires, glass, fiber
and paper production etc.
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: 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: 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: 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 considers the effect of heat generation
proportional l to (T - T∞ )p , where T is the local temperature and T∞
is the ambient temperature, in unsteady free convection flow near the
stagnation point region of a three-dimensional body. The fluid is
considered in an ambient fluid under the assumption of a step change
in the surface temperature of the body. The non-linear coupled partial
differential equations governing the free convection flow are solved
numerically using an implicit finite-difference method for different
values of the governing parameters entering these equations. The
results for the flow and heat characteristics when p ≤ 2 show that
the transition from the initial unsteady-state flow to the final steadystate
flow takes place smoothly. The behavior of the flow is seen
strongly depend on the exponent p.
Abstract: Flows over a harmonically oscillating NACA 0012
airfoil are simulated here using a two-dimensional, unsteady,
incompressibleNavier-Stokes solver.Both pure-plunging and
pitching-plunging combined oscillations are considered at a Reynolds
number of 5000. Special attention is paid to the vortex shedding and
interaction mechanism of the motions. For all the simulations
presented here, the reduced frequency (k) is fixed at a value of 2.5
and plunging amplitude (h) is selected to be in the range of 0.2-0.5.
The simulation results show that the interaction mechanism between
the leading and trailing edge vortices has a decisive effect on the
values of the resulting thrust and propulsive efficiency.
Abstract: The unsteady wake of an EPPLER 361 airfoil in
pitching motion has been investigated in a subsonic wind tunnel by
hot-wire anemometry. The airfoil was given the pitching motion
about the one-quarter chord axis at reduced frequency of 0182.
Streamwise mean velocity profiles (wake profiles) were investigated
at several vertically aligned points behind the airfoil at one-quarter
chord downstream distance from trailing edge. Oscillation amplitude
and mean angle of attack were varied to determine the effects on
wake profiles. When the maximum dynamic angle of attack was
below the static stall angle of attack, weak effects on wake were
found by increasing oscillation amplitude and mean angle of attack.
But, for higher angles of attack strong unsteady effects were
appeared on the wake.