Abstract: In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function ψ (x,y)and the function φ (x,y)as independent variables where for irrotational flow φ (x,y)can be recognized as the velocity potential function, for rotational flow φ (x,y)ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on the finite difference scheme on a uniform mesh is employed. The technique described is capable of tackling the so-called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct geometries. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists.
Abstract: In this paper, a new approach is introduced to solve
Blasius equation using parameter identification of a nonlinear
function which is used as approximation function. Bees Algorithm
(BA) is applied in order to find the adjustable parameters of
approximation function regarding minimizing a fitness function
including these parameters (i.e. adjustable parameters). These
parameters are determined how the approximation function has to
satisfy the boundary conditions. In order to demonstrate the
presented method, the obtained results are compared with another
numerical method. Present method can be easily extended to solve a
wide range of problems.
Abstract: In the present work, study of the vibration of thin cylindrical shells made of a functionally gradient material (FGM) composed of stainless steel and nickel is presented. Material properties are graded in the thickness direction of the shell according to volume fraction power law distribution. The objective is to study the natural frequencies, the influence of constituent volume fractions and the effects of boundary conditions on the natural frequencies of the FG cylindrical shell. The study is carried out using third order shear deformation shell theory. The analysis is carried out using Hamilton's principle. The governing equations of motion of FG cylindrical shells are derived based on shear deformation theory. Results are presented on the frequency characteristics, influence of constituent volume fractions and the effects of free-free and clamped-clamped boundary conditions.
Abstract: In this paper, the dam-reservoir interaction is
analyzed using a finite element approach. The fluid is assumed to be
incompressible, irrotational and inviscid. The assumed boundary
conditions are that the interface of the dam and reservoir is vertical
and the bottom of reservoir is rigid and horizontal. The governing
equation for these boundary conditions is implemented in the
developed finite element code considering the horizontal and vertical
earthquake components. The weighted residual standard Galerkin
finite element technique with 8-node elements is used to discretize
the equation that produces a symmetric matrix equation for the damreservoir
system. A new boundary condition is proposed for
truncating surface of unbounded fluid domain to show the energy
dissipation in the reservoir, through radiation in the infinite upstream
direction. The Sommerfeld-s and perfect damping boundary
conditions are also implemented for a truncated boundary to compare
with the proposed far end boundary. The results are compared with
an analytical solution to demonstrate the accuracy of the proposed
formulation and other truncated boundary conditions in modeling the
hydrodynamic response of an infinite reservoir.
Abstract: Two-dimensional heat conduction within a composed solid material with a constant internal heat generation has been investigated numerically in a sector of the rotor a generator. The heat transfer between two adjacent materials is assumed to be purely conduction. Boundary conditions are assumed to be forced convection on the fluid side and adiabatic on symmetry lines. The control volume method is applied for the diffusion energy equation. Physical coordinates are transformed to the general curvilinear coordinates. Then by using a line-by-line method, the temperature distribution in a sector of the rotor has been determined. Finally, the results are normalized and the effect of cooling fluid on the maximum temperature of insulation is investigated.
Abstract: In this research, a 2-D computational analysis of
steady state free convection in a rectangular enclosure filled with an
electrically conducting fluid under Effect of Magnetic Field has been
performed. The governing equations (mass, momentum, and energy)
are formulated and solved by a finite volume method (FVM)
subjected to different boundary conditions. A parametric study has
been conducted to consider the influence of Grashof number (Gr),
Prantdl number (Pr) and the orientation of magnetic field on the flow
and heat transfer characteristics. It is observed that Nusselt number
(Nu) and heat flux will increase with increasing Grashof and Prandtl
numbers and decreasing the slope of the orientation of magnetic field.
Abstract: A numerical analysis used to simulate the effects of wavy surfaces and thermal radiation on natural convection heat transfer boundary layer flow over an inclined wavy plate has been investigated. A simple coordinate transformation is employed to transform the complex wavy surface into a flat plate. The boundary layer equations and the boundary conditions are discretized by the finite difference scheme and solved numerically using the Gauss-Seidel algorithm with relaxation coefficient. Effects of the wavy geometry, the inclination angle of the wavy plate and the thermal radiation on the velocity profiles, temperature profiles and the local Nusselt number are presented and discussed in detail.
Abstract: Natural convection heat transfer from a heated
horizontal semi-circular cylinder (flat surface upward) has been
investigated for the following ranges of conditions; Grashof number,
and Prandtl number. The governing partial differential equations
(continuity, Navier-Stokes and energy equations) have been solved
numerically using a finite volume formulation. In addition, the role of
the type of the thermal boundary condition imposed at cylinder
surface, namely, constant wall temperature (CWT) and constant heat
flux (CHF) are explored. Natural convection heat transfer from a
heated horizontal semi-circular cylinder (flat surface upward) has
been investigated for the following ranges of conditions; Grashof
number, and Prandtl number, . The governing partial differential
equations (continuity, Navier-Stokes and energy equations) have
been solved numerically using a finite volume formulation. In
addition, the role of the type of the thermal boundary condition
imposed at cylinder surface, namely, constant wall temperature
(CWT) and constant heat flux (CHF) are explored. The resulting flow
and temperature fields are visualized in terms of the streamline and
isotherm patterns in the proximity of the cylinder. The flow remains
attached to the cylinder surface over the range of conditions spanned
here except that for and ; at these conditions, a separated flow
region is observed when the condition of the constant wall
temperature is prescribed on the surface of the cylinder. The heat
transfer characteristics are analyzed in terms of the local and average
Nusselt numbers. The maximum value of the local Nusselt number
always occurs at the corner points whereas it is found to be minimum
at the rear stagnation point on the flat surface. Overall, the average
Nusselt number increases with Grashof number and/ or Prandtl
number in accordance with the scaling considerations. The numerical
results are used to develop simple correlations as functions of
Grashof and Prandtl number thereby enabling the interpolation of the
present numerical results for the intermediate values of the Prandtl or
Grashof numbers for both thermal boundary conditions.
Abstract: Dynamic shear test on simulated phantom can be used
to validate magnetic resonance elastography (MRE) measurements.
Phantom gel has been usually utilized for the cell culture of cartilage
and soft tissue and also been used for mechanical property
characterization using imaging systems. The viscoelastic property of
the phantom would be important for dynamic experiments and
analyses. In this study, An axisymmetric FE model is presented for
determining the dynamic shear behaviour of brain simulated phantom
using ABAQUS. The main objective of this study was to investigate
the effect of excitation frequencies and boundary conditions on shear
modulus and shear viscosity in viscoelastic media.
Abstract: In the present paper the results of a numerical study are presented, numerical models were developed to simulate the behaviour of vertical massive dikes. The proposed models were developed according to the geometry, boundary conditions, loading conditions and initial conditions of a physical model taken as reference. The results obtained were compared to the experimental data. As far as the overall behaviour, the displacements and the failure mechanisms of the dikes is concerned, the numerical results were in good agreement with the experimental results, which clearly indicates a good quality of numerical modelling. The validated numerical models were used in a parametric study were the displacements and failure mechanisms were fully investigated. Out of the results obtained, some conclusions and recommendations related to the design of massive dikes are proposed.
Abstract: A strip domain decomposition parallel algorithm for fast direct Poisson solver is presented on a 3D Cartesian staggered grid. The parallel algorithm follows the principles of sequential algorithm for fast direct Poisson solver. Both Dirichlet and Neumann boundary conditions are addressed. Several test cases are likewise addressed in order to shed light on accuracy and efficiency in the strip domain parallelization algorithm. Actually the current implementation shows a very high efficiency when dealing with a large grid mesh up to 3.6 * 109 under massive parallel approach, which explicitly demonstrates that the proposed algorithm is ready for massive parallel computing.
Abstract: In this paper, the melting of a semi-infinite body as a
result of a moving laser beam has been studied. Because the Fourier
heat transfer equation at short times and large dimensions does not
have sufficient accuracy; a non-Fourier form of heat transfer
equation has been used. Due to the fact that the beam is moving in x
direction, the temperature distribution and the melting pool shape are
not asymmetric. As a result, the problem is a transient threedimensional
problem. Therefore, thermophysical properties such as
heat conductivity coefficient, density and heat capacity are functions
of temperature and material states. The enthalpy technique, used for
the solution of phase change problems, has been used in an explicit
finite volume form for the hyperbolic heat transfer equation. This
technique has been used to calculate the transient temperature
distribution in the semi-infinite body and the growth rate of the melt
pool. In order to validate the numerical results, comparisons were
made with experimental data. Finally, the results of this paper were
compared with similar problem that has used the Fourier theory. The
comparison shows the influence of infinite speed of heat propagation
in Fourier theory on the temperature distribution and the melt pool
size.
Abstract: In this paper we present our results on the performance analysis of a multi-product manufacturing line. We study the influence of external perturbations, intermediate buffer content and the number of manufacturing stages on the production tracking error of each machine in the multi-product line operated under a surplusbased production control policy. Starting by the analysis of a single machine with multiple production stages (one for each product type), we provide bounds on the production error of each stage. Then, we extend our analysis to a line of multi-stage machines, where similarly, bounds on each production tracking error for each product type, as well as buffer content are obtained. Details on performance of the closed-loop flow line model are illustrated in numerical simulations.
Abstract: For many industrial applications plate heat
exchangers are demonstrating a large superiority over the
other types of heat exchangers. The efficiency of such a
device depends on numerous factors the effect of which needs
to be analysed and accurately evaluated.
In this paper we present a theoretical analysis of a cocurrent
plate heat exchanger and the results of its numerical
simulation.
Knowing the hot and the cold fluid streams inlet temperatures,
the respective heat capacities mCp
and the value of the
overall heat transfer coefficient, a 1-D mathematical model
based on the steady flow energy balance for a differential
length of the device is developed resulting in a set of N first
order differential equations with boundary conditions where N
is the number of channels.For specific heat exchanger
geometry and operational parameters, the problem is
numerically solved using the shooting method.
The simulation allows the prediction of the temperature
map in the heat exchanger and hence, the evaluation of its
performances. A parametric analysis is performed to evaluate
the influence of the R-parameter on the e-NTU values. For
practical purposes effectiveness-NTU graphs are elaborated
for specific heat exchanger geometry and different operating
conditions.
Abstract: To calculate the temperature distribution of the slab in
a hot rolled reheating furnace a mathematical model has been
developed by considering the thermal radiation in the furnace and
transient conduction in the slab. The furnace is modeled as radiating
medium with spatially varying temperature. Radiative heat flux within
the furnace including the effect of furnace walls, combustion gases,
skid beams and buttons is calculated using the FVM and is applied as
the boundary condition of the transient conduction equation of the
slab. After determining the slab emissivity by comparison between
simulation and experimental work, variation of heating characteristics
in the slab is investigated in the case of changing furnace temperature
with various time and the slab residence time is optimized with this
evaluation.
Abstract: The stability characteristics of water lubricated journal bearings having three axial grooves are obtained theoretically. In this lubricant (water) is fed under pressure from one end of the bearing, through the 3-axial grooves (groove angles may vary). These bearings can use the process fluid as the lubricant, as in the case of feed water pumps. The Reynolds equation is solved numerically by the finite difference method satisfying the boundary conditions. The stiffness and damping coefficient for various bearing number and eccentricity ratios, assuming linear pressure drop along the groove, shows that smaller groove angles better results.
Abstract: Calcium [Ca2+] dynamics is studied as a potential form
of neuron excitability that can control many irregular processes like
metabolism, secretion etc. Ca2+ ion enters presynaptic terminal and
increases the synaptic strength and thus triggers the neurotransmitter
release. The modeling and analysis of calcium dynamics in neuron
cell becomes necessary for deeper understanding of the processes
involved. A mathematical model has been developed for cylindrical
shaped neuron cell by incorporating physiological parameters like
buffer, diffusion coefficient, and association rate. Appropriate initial
and boundary conditions have been framed. The closed form solution
has been developed in terms of modified Bessel function. A computer
program has been developed in MATLAB 7.11 for the whole
approach.
Abstract: In this paper, an analytical approach is used to study the coupled lateral-torsional vibrations of laminated composite beam. It is known that in such structures due to the fibers orientation in various layers, any lateral displacement will produce a twisting moment. This phenomenon is modeled by the bending-twisting material coupling rigidity and its main feature is the coupling of lateral and torsional vibrations. In addition to the material coupling, the effects of shear deformation and rotary inertia are taken into account in the definition of the potential and kinetic energies. Then, the governing differential equations are derived using the Hamilton-s principle and the mathematical model matches the Timoshenko beam model when neglecting the effect of bending-twisting rigidity. The equations of motion which form a system of three coupled PDEs are solved analytically to study the free vibrations of the beam in lateral and rotational modes due to the bending, as well as the torsional mode caused by twisting. The analytic solution is carried out in three steps: 1) assuming synchronous motion for the kinematic variables which are the lateral, rotational and torsional displacements, 2) solving the ensuing eigenvalue problem which contains three coupled second order ODEs and 3) imposing different boundary conditions related to combinations of simply, clamped and free end conditions. The resulting natural frequencies and mode shapes are compared with similar results in the literature and good agreement is achieved.
Abstract: This paper is concerned with an investigation into the
localized non-stability of a thin elastic orthotropic semi-infinite plate.
In this study, a semi-infinite plate, simply supported on two edges
and different boundary conditions, clamped, hinged, sliding contact
and free on the other edge, are considered. The mathematical model
is used and a general solution is presented the conditions under which
localized solutions exist are investigated.
Abstract: An accurate procedure to determine free vibrations of
beams and plates is presented.
The natural frequencies are exact solutions of governing vibration
equations witch load to a nonlinear homogeny system.
The bilinear and linear structures considered simulate a bridge.
The dynamic behavior of this one is analyzed by using the theory of
the orthotropic plate simply supported on two sides and free on the
two others. The plate can be excited by a convoy of constant or
harmonic loads. The determination of the dynamic response of the
structures considered requires knowledge of the free frequencies and
the shape modes of vibrations. Our work is in this context. Indeed,
we are interested to develop a self-consistent calculation of the Eigen
frequencies.
The formulation is based on the determination of the solution of
the differential equations of vibrations. The boundary conditions
corresponding to the shape modes permit to lead to a homogeneous
system. Determination of the noncommonplace solutions of this
system led to a nonlinear problem in Eigen frequencies.
We thus, develop a computer code for the determination of the
eigenvalues. It is based on a method of bisection with interpolation
whose precision reaches 10 -12. Moreover, to determine the
corresponding modes, the calculation algorithm that we develop uses
the method of Gauss with a partial optimization of the "pivots"
combined with an inverse power procedure. The Eigen frequencies
of a plate simply supported along two opposite sides while
considering the two other free sides are thus analyzed. The results
could be generalized with the case of a beam by regarding it as a
plate with low width.
We give, in this paper, some examples of treated cases. The
comparison with results presented in the literature is completely
satisfactory.