Abstract: The temperature distribution and the heat transfer
rates through a multi-layer door of a furnace were investigated. The
inside of the door was in contact with hot air and the other side of the
door was in contact with room air. Radiation heat transfer from the
walls of the furnace to the door and the door to the surrounding area
was included in the problem. This work is a two dimensional steady
state problem. The Churchill and Chu correlation was used to find
local convection heat transfer coefficients at the surfaces of the
furnace door. The thermophysical properties of air were the functions
of the temperatures. Polynomial curve fitting for the fluid properties
were carried out. Finite difference method was used to discretize for
conduction heat transfer within the furnace door. The Gauss-Seidel
Iteration was employed to compute the temperature distribution in
the door.
The temperature distribution in the horizontal mid plane of the
furnace door in a two dimensional problem agrees with the one
dimensional problem. The local convection heat transfer coefficients
at the inside and outside surfaces of the furnace door are exhibited.
Abstract: Saccharomyces cerevisiae (baker-s yeast) can exhibit
sustained oscillations during the operation in a continuous bioreactor
that adversely affects its stability and productivity. Because of
heterogeneous nature of cell populations, the cell population balance
models can be used to capture the dynamic behavior of such cultures.
In this paper an unstructured, segregated model is used which is
based on population balance equation(PBE) and then in order to
simulation, the 4th order Rung-Kutta is used for time dimension and
three methods, finite difference, orthogonal collocation on finite
elements and Galerkin finite element are used for discretization of the
cell mass domain. The results indicate that the orthogonal collocation
on finite element not only is able to predict the oscillating behavior of
the cell culture but also needs much little time for calculations.
Therefore this method is preferred in comparison with other methods.
In the next step two controllers, a globally linearizing control (GLC)
and a conventional proportional-integral (PI) controller are designed
for controlling the total cell mass per unit volume, and performances
of these controllers are compared through simulation. The results
show that although the PI controller has simpler structure, the GLC
has better performance.
Abstract: In this paper we use quintic non-polynomial
spline functions to develop numerical methods for approximation
to the solution of a system of fourth-order boundaryvalue
problems associated with obstacle, unilateral and contact
problems. The convergence analysis of the methods has been
discussed and shown that the given approximations are better
than collocation and finite difference methods. Numerical
examples are presented to illustrate the applications of these
methods, and to compare the computed results with other
known methods.
Abstract: An analysis is performed to study the influence of nonuniform double slot suction on a steady laminar boundary layer flow over a rotating sphere when fluid properties such as viscosity and Prandtl number are inverse linear functions of temperature. Nonsimilar solutions have been obtained from the starting point of the streamwise co-ordinate to the exact point of separation. The difficulties arising at the starting point of the streamwise co-ordinate, at the edges of the slot and at the point of separation have been overcome by applying an implicit finite difference scheme in combination with the quasi-linearization technique and an appropriate selection of the finer step sizes along the stream-wise direction. The present investigation shows that the point of ordinary separation can be delayed by nonuniform double slot suction if the mass transfer rate is increased and also if the slots are positioned further downstream. In addition, the investigation reveals that double slot suction is found to be more effective compared to a single slot suction in delaying ordinary separation. As rotation parameter increase the point of separation moves upstream direction.
Abstract: Based on Rayleigh beam theory, the sub-impacts of a
free-free beam struck horizontally by a round-nosed rigid mass is
simulated by the finite difference method and the impact-separation
conditions. In order to obtain the sub-impact force, a uniaxial
compression elastic-plastic contact model is employed to analyze the
local deformation field on contact zone. It is found that the horizontal
impact is a complicated process including the elastic plastic
sub-impacts in sequence. There are two sub-zones of sub-impact. In
addition, it found that the elastic energy of the free-free beam is more
suitable for the Poisson collision hypothesis to explain compression
and recovery processes.
Abstract: Many computational techniques were applied to
solution of heat conduction problem. Those techniques were the
finite difference (FD), finite element (FE) and recently meshless
methods. FE is commonly used in solution of equation of heat
conduction problem based on the summation of stiffness matrix of
elements and the solution of the final system of equations. Because
of summation process of finite element, convergence rate was
decreased. Hence in the present paper Cellular Automata (CA)
approach is presented for the solution of heat conduction problem.
Each cell considered as a fixed point in a regular grid lead to the
solution of a system of equations is substituted by discrete systems of
equations with small dimensions. Results show that CA can be used
for solution of heat conduction problem.
Abstract: This paper deals with a numerical analysis of the
transient response of composite beams with strain rate dependent
mechanical properties by use of a finite difference method. The
equations of motion based on Timoshenko beam theory are derived.
The geometric nonlinearity effects are taken into account with von
Kármán large deflection theory. The finite difference method in
conjunction with Newmark average acceleration method is applied to
solve the differential equations. A modified progressive damage
model which accounts for strain rate effects is developed based on
the material property degradation rules and modified Hashin-type
failure criteria and added to the finite difference model. The
components of the model are implemented into a computer code in
Mathematica 6. Glass/epoxy laminated composite beams with
constant and strain rate dependent mechanical properties under
dynamic load are analyzed. Effects of strain rate on dynamic
response of the beam for various stacking sequences, load and
boundary conditions are investigated.
Abstract: This paper proposes the numerical simulation of the
investment casting of gold jewelry. It aims to study the behavior of
fluid flow during mould filling and solidification and to optimize the
process parameters, which lead to predict and control casting defects
such as gas porosity and shrinkage porosity. A finite difference
method, computer simulation software FLOW-3D was used to
simulate the jewelry casting process. The simplified model was
designed for both numerical simulation and real casting production.
A set of sensor acquisitions were allocated on the different positions
of the wax tree of the model to detect filling times, while a set of
thermocouples were allocated to detect the temperature during
casting and cooling. Those detected data were applied to validate the
results of the numerical simulation to the results of the real casting.
The resulting comparisons signify that the numerical simulation can
be used as an effective tool in investment-casting-process
optimization and casting-defect prediction.
Abstract: Dynamic characteristics of a four-lobe journal bearing
of micropolar fluids are presented. Lubricating oil containing
additives and contaminants is modelled as micropolar fluid. The
modified Reynolds equation is obtained using the micropolar
lubrication theory and solving it by using finite difference technique.
The dynamic characteristics in terms of stiffness, damping
coefficients, the critical mass and whirl ratio are determined for
various values of size of material characteristic length and the
coupling number. The results show compared with Newtonian fluids,
that micropolar fluid exhibits better stability.