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: In this paper, based on flume experimental data, the velocity distribution in open channel flows is re-investigated. From the analysis, it is proposed that the wake layer in outer region may be divided into two regions, the relatively weak outer region and the relatively strong outer region. Combining the log law for inner region and the parabolic law for relatively strong outer region, an explicit equation for mean velocity distribution of steady and uniform turbulent flow through straight open channels is proposed and verified with the experimental data. It is found that the sediment concentration has significant effect on velocity distribution in the relatively weak outer region.
Abstract: This paper presents a procedure of forming the
mathematical model of radial electric power systems for simulation
of both transient and steady-state conditions. The research idea has
been based on nodal voltages technique and on differentiation of
Kirchhoff's current law (KCL) applied to each non-reference node of
the radial system, the result of which the nodal voltages has been
calculated by solving a system of algebraic equations. Currents of the
electric power system components have been determined by solving
their respective differential equations. Transforming the three-phase
coordinate system into Cartesian coordinate system in the model
decreased the overall number of equations by one third. The use of
Cartesian coordinate system does not ignore the DC component
during transient conditions, but restricts the model's implementation
for symmetrical modes of operation only. An example of the input
data for a four-bus radial electric power system has been calculated.
Abstract: The steady-state temperature for one-dimensional transpiration cooling system has been conducted experimentally and numerically to investigate the heat transfer characteristics of combined convection and radiation. The Nickel –Chrome (Ni-Cr) open-cellular porous material having porosity of 0.93 and pores per inch (PPI) of 21.5 was examined. The upper surface of porous plate was heated by the heat flux of incoming radiation varying from 7.7 - 16.6 kW/m2 whereas air injection velocity fed into the lower surface was varied from 0.36 - 1.27 m/s, and was then rearranged as Reynolds number (Re). For the report of the results in the present study, two efficiencies including of temperature and conversion efficiency were presented. Temperature efficiency indicating how close the mean temperature of a porous heat plate to that of inlet air, and increased rapidly with the air injection velocity (Re). It was then saturated and had a constant value at Re higher than 10. The conversion efficiency, which was regarded as the ability of porous material in transferring energy by convection after absorbed from heat radiation, decreased with increasing of the heat flux and air injection velocity. In addition, it was then asymptotic to a constant value at the Re higher than 10. The numerical predictions also agreed with experimental data very well.
Abstract: Multiport diffusers are the effective engineering
devices installed at the modern marine outfalls for the steady
discharge of effluent streams from the coastal plants, such as
municipal sewage treatment, thermal power generation and seawater
desalination. A mathematical model using a two-dimensional
advection-diffusion equation based on a flat seabed and incorporating
the effect of a coastal tidal current is developed to calculate the
compounded concentration following discharges of desalination
brine from a sea outfall with multiport diffusers. The analytical
solutions are computed graphically to illustrate the merging of
multiple brine plumes in shallow coastal waters, and further
approximation will be made to the maximum shoreline's
concentration to formulate dilution of a multiport diffuser discharge.
Abstract: This paper proposes an effective adaptation learning
algorithm based on artificial neural networks for speed control of an
induction motor assumed to operate in a high-performance drives
environment. The structure scheme consists of a neural network
controller and an algorithm for changing the NN weights in order that
the motor speed can accurately track of the reference command. This
paper also makes uses a very realistic and practical scheme to
estimate and adaptively learn the noise content in the speed load
torque characteristic of the motor. The availability of the proposed
controller is verified by through a laboratory implementation and
under computation simulations with Matlab-software. The process is
also tested for the tracking property using different types of reference
signals. The performance and robustness of the proposed control
scheme have evaluated under a variety of operating conditions of the
induction motor drives. The obtained results demonstrate the
effectiveness of the proposed control scheme system performances,
both in steady state error in speed and dynamic conditions, was found
to be excellent and those is not overshoot.
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: Dynamic models of power converters are normally
time-varying because of their switching actions. Several approaches
are applied to analyze the power converters to achieve the timeinvariant
models suitable for system analysis and design via the
classical control theory. The paper presents how to derive dynamic
models of the power system consisting of a three-phase controlled
rectifier feeding an uncontrolled buck converter by using the
combination between the well known techniques called the DQ and
the generalized state-space averaging methods. The intensive timedomain
simulations of the exact topology model are used to support
the accuracies of the reported model. The results show that the
proposed model can provide good accuracies in both transient and
steady-state responses.
Abstract: A biophysically based multilayer continuum model of the facial soft tissue composite has been developed for simulating wrinkle formation. The deformed state of the soft tissue block was determined by solving large deformation mechanics equations using the Galerkin finite element method. The proposed soft tissue model is composed of four layers with distinct mechanical properties. These include stratum corneum, epidermal-dermal layer (living epidermis and dermis), subcutaneous tissue and the underlying muscle. All the layers were treated as non-linear, isotropic Mooney Rivlin materials. Contraction of muscle fibres was approximated using a steady-state relationship between the fibre extension ratio, intracellular calcium concentration and active stress in the fibre direction. Several variations of the model parameters (stiffness and thickness of epidermal-dermal layer, thickness of subcutaneous tissue layer) have been considered.
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: The authors present an optimization algorithm for order reduction and its application for the determination of the relative mapping errors of linear time invariant dynamic systems by the simplified models. These relative mapping errors are expressed by means of the relative integral square error criterion, which are determined for both unit step and impulse inputs. The reduction algorithm is based on minimization of the integral square error by particle swarm optimization technique pertaining to a unit step input. The algorithm is simple and computer oriented. It is shown that the algorithm has several advantages, e.g. the reduced order models retain the steady-state value and stability of the original system. Two numerical examples are solved to illustrate the superiority of the algorithm over some existing methods.
Abstract: Steady state experiments have been conducted for
natural and mixed convection heat transfer, from five different sized
protruding discrete heat sources, placed at the bottom position on a
PCB and mounted on a vertical channel. The characteristic length (
Lh ) of heat sources vary from 0.005 to 0.011 m. The study has been
done for different range of Reynolds number and modified Grashof
number. From the experiment, the surface temperature distribution
and the Nusselt number of discrete heat sources have been obtained
and the effects of Reynold number and Richardson number on them
have been discussed. The objective is to find the rate of heat
dissipation from heat sources, by placing them at the bottom position
on a PCB and to compare both modes of cooling of heat sources.
Abstract: Wind catchers are traditional natural ventilation
systems attached to buildings in order to ventilate the indoor air. The
most common type of wind catcher is four sided one which is
capable to catch wind in all directions. CFD simulation is the perfect
way to evaluate the wind catcher performance. The accuracy of CFD
results is the issue of concern, so sensitivity analyses is crucial to
find out the effect of different settings of CFD on results. This paper
presents a series of 3D steady RANS simulations for a generic
isolated four-sided wind catcher attached to a room subjected to wind
direction ranging from 0º to 180º with an interval of 45º. The CFD
simulations are validated with detailed wind tunnel experiments. The
influence of an extensive range of computational parameters is
explored in this paper, including the resolution of the computational
grid, the size of the computational domain and the turbulence model.
This study found that CFD simulation is a reliable method for wind
catcher study, but it is less accurate in prediction of models with non
perpendicular wind directions.
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: In order to obtain an accurate result of the heat transfer
of the rib in the internal cooling Rectangular channel, using separation
of variables, analytical solutions of three dimensional steady-state heat
conduction in rectangular ribs are given by solving three dimensional
steady-state function of the rectangular ribs. Therefore, we can get
solution of three dimensional temperature field in the rib. Based on the
solution, we can get how the Bi number affected on heat transfer.
Furthermore, comparisons of the analytical and numerical results
indicate agreement on temperature field in the rib.
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: Saturated two-phase fluid flows are often subject to
pressure induced oscillations. Due to compressibility the vapor
bubbles act as a spring with an asymmetric non-linear characteristic.
The volume of the vapor bubbles increases or decreases differently if
the pressure fluctuations are compressing or expanding;
consequently, compressing pressure fluctuations in a two-phase pipe
flow cause less displacement in the direction of the pipe flow than
expanding pressure fluctuations. The displacement depends on the
ratio of liquid to vapor, the ratio of pressure fluctuations over average
pressure and on the exciting frequency of the pressure fluctuations.
In addition, pressure fluctuations in saturated vapor bubbles cause
condensation and evaporation within the bubbles and change
periodically the ratio between liquid to vapor, and influence the
dynamical parameters for the oscillation. The oscillations are
conforming to an isenthalpic process at constant enthalpy with no
heat transfer and no exchange of work.
The paper describes the governing non-linear equation for twophase
fluid oscillations with condensation and evaporation, and
presents steady state approximate solutions for free and for pressure
induced oscillations. Resonance criteria and stability are discussed.
Abstract: In this paper, three dimensional flow characteristic was
presented by a revision of an impeller of an axial turbo fan for
improving the airflow rate and the static pressure. TO consider an
incompressible steady three-dimensional flow, the RANS equations
are used as the governing equations, and the standard k-ε turbulence
model is chosen. The pitch angles of 44°, 54°, 59°, and 64° are
implemented for the numerical model. The numerical results show that
airflow rates of each pitch angle are 1,175 CMH, 1,270 CMH, 1,340
CMH, and 800 CMH, respectively. The difference of the static
pressure at impeller inlet and outlet are 120 Pa, 214 Pa, 242 Pa, and 60
Pa according to respective pitch angles. It means that the 59° of the
impeller pitch angle is optimal to improve the airflow rate and the
static pressure.
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.