Abstract: The radiative exchange method is introduced as a
numerical method for the simulation of radiative heat transfer in an
absorbing, emitting and isotropically scattering media. In this
method, the integro-differential radiative balance equation is solved
by using a new introduced concept for the exchange factor. Even
though the radiative source term is calculated in a mesh structure that
is coarser than the structure used in computational fluid dynamics,
calculating the exchange factor between different coarse elements by
using differential integration elements makes the result of the method
close to that of integro-differential radiative equation. A set of
equations for calculating exchange factors in two and threedimensional
Cartesian coordinate system is presented, and the
method is used in the simulation of radiative heat transfer in twodimensional
rectangular case and a three-dimensional simple cube.
The result of using this method in simulating different cases is
verified by comparing them with those of using other numerical
radiative models.
Abstract: Catalytic converters are used for minimizing the release of pollutants to the atmosphere. It is during the warm-up period that hydrocarbons are seen to be released in appreciable quantities from these converters. In this paper the conversion of a fast oxidizing hydrocarbon propylene is analysed using two numerical methods. The quasi steady state method assumes the accumulation terms to be negligible in the gas phase mass and energy balance equations, however this term is present in the solid phase energy balance. The unsteady state model accounts for the accumulation term to be present in the gas phase mass and energy balance and in the solid phase energy balance. The results derived from the two models for gas concentration, gas temperature and solid temperature are compared.
Abstract: In the present paper some recommendations for the
use of software package “Mathematica" in a basic numerical analysis
course are presented. The methods which are covered in the course
include solution of systems of linear equations, nonlinear equations
and systems of nonlinear equations, numerical integration,
interpolation and solution of ordinary differential equations. A set of
individual assignments developed for the course covering all the
topics is discussed in detail.
Abstract: The Eulerian numerical method is proposed to analyze
the explosion in tunnel. Based on this method, an original software
M-MMIC2D is developed by Cµ program language. With this
software, the explosion problem in the tunnel with three
expansion-chambers is numerically simulated, and the results are
found to be in full agreement with the observed experimental data.
Abstract: This study presents a hybrid neural network and Gravitational Search Algorithm (HNGSA) method to solve well known Wessinger's equation. To aim this purpose, gravitational search algorithm (GSA) technique is applied to train a multi-layer perceptron neural network, which is used as approximation solution of the Wessinger's equation. A trial solution of the differential equation is written as sum of two parts. The first part satisfies the initial/ boundary conditions and does not contain any adjustable parameters and the second part which is constructed so as not to affect the initial/boundary conditions. The second part involves adjustable parameters (the weights and biases) for a multi-layer perceptron neural network. In order to demonstrate the presented method, the obtained results of the proposed method are compared with some known numerical methods. The given results show that presented method can introduce a closer form to the analytic solution than other numerical methods. Present method can be easily extended to solve a wide range of problems.
Abstract: In contrast to existing methods which do not take into account multiconnectivity in a broad sense of this term, we develop mathematical models and highly effective combination (BIEM and FDM) numerical methods of calculation of stationary and quasistationary temperature field of a profile part of a blade with convective cooling (from the point of view of realization on PC). The theoretical substantiation of these methods is proved by appropriate theorems. For it, converging quadrature processes have been developed and the estimations of errors in the terms of A.Ziqmound continuity modules have been received. For visualization of profiles are used: the method of the least squares with automatic conjecture, device spline, smooth replenishment and neural nets. Boundary conditions of heat exchange are determined from the solution of the corresponding integral equations and empirical relationships. The reliability of designed methods is proved by calculation and experimental investigations heat and hydraulic characteristics of the gas turbine first stage nozzle blade.
Abstract: As is known, one of the priority directions of research
works of natural sciences is introduction of applied section of
contemporary mathematics as approximate and numerical methods to
solving integral equation into practice. We fare with the solving of
integral equation while studying many phenomena of nature to whose
numerically solving by the methods of quadrature are mainly applied.
Taking into account some deficiency of methods of quadrature for
finding the solution of integral equation some sciences suggested of
the multistep methods with constant coefficients. Unlike these papers,
here we consider application of hybrid methods to the numerical
solution of Volterra integral equation. The efficiency of the suggested
method is proved and a concrete method with accuracy order p = 4
is constructed. This method in more precise than the corresponding
known methods.
Abstract: In contrast to existing methods which do not take into account multiconnectivity in a broad sense of this term, we develop mathematical models and highly effective combination (BIEM and FDM) numerical methods of calculation of stationary and quasi-stationary temperature field of a profile part of a blade with convective cooling (from the point of view of realization on PC). The theoretical substantiation of these methods is proved by appropriate theorems. For it, converging quadrature processes have been developed and the estimations of errors in the terms of A.Ziqmound continuity modules have been received. For visualization of profiles are used: the method of the least squares with automatic conjecture, device spline, smooth replenishment and neural nets. Boundary conditions of heat exchange are determined from the solution of the corresponding integral equations and empirical relationships. The reliability of designed methods is proved by calculation and experimental investigations heat and hydraulic characteristics of the gas turbine first stage nozzle blade.
Abstract: Impinging jets are used in various industrial areas as a cooling and drying technique. The current research is concerned with the means of improving the heat transfer for configurations with a minimum distance of the nozzle to the impingement surface. The impingement heat transfer is described using numerical methods over a wide range of parameters for an array of planar jets. These parameters include varying jet flow speed, width of nozzle, distance of nozzle, angle of the jet flow, velocity and geometry of the impingement surface. Normal pressure and shear stress are computed as additional parameters. Using dimensionless characteristic numbers the parameters and the results are correlated to gain generalized equations. The results demonstrate the effect of the investigated parameters on the flow.
Abstract: In this work the numerical simulation of transient heat
transfer in a cylindrical probe is done. An experiment was conducted
introducing a steel cylinder in a heating chamber and registering its
surface temperature along the time during one hour. In parallel, a
mathematical model was solved for one dimension transient heat
transfer in cylindrical coordinates, considering the boundary
conditions of the test. The model was solved using finite difference
method, because the thermal conductivity in the cylindrical steel bar
and the convection heat transfer coefficient used in the model are
considered temperature dependant functions, and both conditions
prevent the use of the analytical solution. The comparison between
theoretical and experimental results showed the average deviation is
below 2%. It was concluded that numerical methods are useful in
order to solve engineering complex problems. For constant k and h,
the experimental methodology used here can be used as a tool for
teaching heat transfer in mechanical engineering, using mathematical
simplified models with analytical solutions.
Abstract: A numerical method for solving the time-independent Schrödinger equation of a particle moving freely in a three-dimensional
axisymmetric region is developed. The boundary of the region
is defined by an arbitrary analytic function. The method uses a
coordinate transformation and an expansion in eigenfunctions. The
effectiveness is checked and confirmed by applying the method to a
particular example, which is a prolate spheroid.
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: To compute dynamic characteristics of nonlinear viscoelastic springs with elastic structures having huge degree-of-freedom, Yamaguchi proposed a new fast numerical method using finite element method [1]-[2]. In this method, restoring forces of the springs are expressed using power series of their elongation. In the expression, nonlinear hysteresis damping is introduced. In this expression, nonlinear complex spring constants are introduced. Finite element for the nonlinear spring having complex coefficients is expressed and is connected to the elastic structures modeled by linear solid finite element. Further, to save computational time, the discrete equations in physical coordinate are transformed into the nonlinear ordinary coupled equations using normal coordinate corresponding to linear natural modes. In this report, the proposed method is applied to simulation for impact responses of a viscoelastic shock absorber with an elastic structure (an S-shaped structure) by colliding with a concentrated mass. The concentrated mass has initial velocities and collides with the shock absorber. Accelerations of the elastic structure and the concentrated mass are measured using Levitation Mass Method proposed by Fujii [3]. The calculated accelerations from the proposed FEM, corresponds to the experimental ones. Moreover, using this method, we also investigate dynamic errors of the S-shaped force transducer due to elastic mode in the S-shaped structure.
Abstract: In recent years, a new numerical method has been
developed, the extended finite element method (X-FEM). The
objective of this work is to exploit the (X-FEM) for the treatment of
the fracture mechanics problems on 3D geometries, where we
showed the ability of this method to simulate the fatigue crack
growth into two cases: edge and central crack. In the results we
compared the six first natural frequencies of mode shapes uncracking
with the cracking initiation in the structure, and showed the stress
intensity factor (SIF) evolution function as crack size propagation
into structure, the analytical validation of (SIF) is presented. For to
evidence the aspects of this method, all result is compared between
FEA and X-FEM.
Abstract: A prototype model of an emulsion separator was
designed and manufactured. Generally, it is a cylinder filled with
different fractal modules. The emulsion was fed into the reactor by a
peristaltic pump through an inlet placed at the boundary between the
two phases. For hydrodynamic design and sizing of the reactor the
assumptions of the theory of filtration were used and methods to
describe the separation process were developed. Based on this
methodology and using numerical methods and software of Autodesk
the process is simulated in different operating modes. The basic
hydrodynamic characteristics - speed and performance for different
types of fractal systems and decisions to optimize the design of the
reactor were also defined.
Abstract: Seismic design may require non-conventional
concept, due to the fact that the stiffness and layout of the structure
have a great effect on the overall structural behaviour, on the seismic
load intensity as well as on the internal force distribution. To find an
economical and optimal structural configuration the key issue is the
optimal design of the lateral load resisting system. This paper focuses
on the optimal design of regular, concentric braced frame (CBF)
multi-storey steel building structures. The optimal configurations are
determined by a numerical method using genetic algorithm approach,
developed by the authors. Aim is to find structural configurations
with minimum structural cost. The design constraints of objective
function are assigned in accordance with Eurocode 3 and Eurocode 8
guidelines. In this paper the results are presented for various building
geometries, different seismic intensities, and levels of energy
dissipation.
Abstract: The ultimate goal of this article is to develop a robust and accurate numerical method for solving hyperbolic conservation laws in one and two dimensions. A hybrid numerical method, coupling a cheap fourth order total variation diminishing (TVD) scheme [1] for smooth region and a Robust seventh-order weighted non-oscillatory (WENO) scheme [2] near discontinuities, is considered. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme, while the smooth regions are computed with fourth order total variation diminishing (TVD). For time integration, we use the third order TVD Runge-Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.
Abstract: Intravitreal injection (IVI) is the most common treatment for eye posterior segment diseases such as endopthalmitis, retinitis, age-related macular degeneration, diabetic retinopathy, uveitis, and retinal detachment. Most of the drugs used to treat vitreoretinal diseases, have a narrow concentration range in which they are effective, and may be toxic at higher concentrations. Therefore, it is critical to know the drug distribution within the eye following intravitreal injection. Having knowledge of drug distribution, ophthalmologists can decide on drug injection frequency while minimizing damage to tissues. The goal of this study was to develop a computer model to predict intraocular concentrations and pharmacokinetics of intravitreally injected drugs. A finite volume model was created to predict distribution of two drugs with different physiochemical properties in the rabbit eye. The model parameters were obtained from literature review. To validate this numeric model, the in vivo data of spatial concentration profile from the lens to the retina were compared with the numeric data. The difference was less than 5% between the numerical and experimental data. This validation provides strong support for the numerical methodology and associated assumptions of the current study.