Abstract: The method of moments combined with the method
of symmetrical components is used for the analysis of interstitial
hyperthermia applicators. The basis and testing functions are both
piecewise sinusoids, qualifying our technique as a Galerkin one. The
dielectric coatings are modeled by equivalent volume polarization
currents, which are simply related to the conduction current
distribution, avoiding in that way the introduction of additional
unknowns or numerical integrations. The results of our method
for a four dipole circular array, are in agreement with those
already published in literature for a same hyperthermia configuration.
Apart from being accurate, our approach is more general, more
computationally efficient and takes into account the coupling between
the antennas.
Abstract: In this paper, a new trend for improvement in semianalytical
method based on scale boundaries in order to solve the 2D
elastodynamic problems is provided. In this regard, only the
boundaries of the problem domain discretization are by specific subparametric
elements. Mapping functions are uses as a class of higherorder
Lagrange polynomials, special shape functions, Gauss-Lobatto-
Legendre numerical integration, and the integral form of the weighted
residual method, the matrix is diagonal coefficients in the equations
of elastodynamic issues. Differences between study conducted and
prior research in this paper is in geometry production procedure of
the interpolation function and integration of the different is selected.
Validity and accuracy of the present method are fully demonstrated
through two benchmark problems which are successfully modeled
using a few numbers of DOFs. The numerical results agree very well
with the analytical solutions and the results from other numerical
methods.
Abstract: A new basis function neural network algorithm is proposed for numerical integration. The main idea is to construct neural network model based on spline basis functions, which is used to approximate the integrand by training neural network weights. The convergence theorem of the neural network algorithm, the theorem for numerical integration and one corollary are presented and proved. The numerical examples, compared with other methods, show that the algorithm is effective and has the characteristics such as high precision and the integrand not required known. Thus, the algorithm presented in this paper can be widely applied in many engineering fields.
Abstract: Consider the Gregory integration (G) formula
with end corrections where h Δ is the forward difference operator with step size h. In this study we prove that can be optimized by
minimizing some of the coefficient k a in the remainder term by
particle swarm optimization. Experimental tests prove that can be rendered a powerful formula for library use.
Abstract: In general dynamic analyses, lower mode response is
of interest, however the higher modes of spatially discretized
equations generally do not represent the real behavior and not affects
to global response much. Some implicit algorithms, therefore, are
introduced to filter out the high-frequency modes using intended
numerical error. The objective of this study is to introduce the
P-method and PC α-method to compare that with dissipation method
and Newmark method through the stability analysis and numerical
example. PC α-method gives more accuracy than other methods
because it based on the α-method inherits the superior properties of the
implicit α-method. In finite element analysis, the PC α-method is more
useful than other methods because it is the explicit scheme and it
achieves the second order accuracy and numerical damping
simultaneously.
Abstract: Nonlinear response behaviour of a cracked RC beam under harmonic excitation is analysed to investigate various instability phenomena like, bifurcation, jump phenomena etc. The nonlinearity of the system arises due to opening and closing of the cracks in the RC beam and is modelled as a cubic polynomial. In order to trace different branches at the bifurcation point on the response curve (amplitude versus frequency of excitation plot), an arc length continuation technique along with the incremental harmonic balance (IHBC) method is employed. The stability of the solution is investigated by the Floquet theory using Hsu-s scheme. The periodic solutions obtained by the IHBC method are compared with these obtained by the numerical integration of the equation of motion. Characteristics of solutions fold bifurcation, jump phenomena and from stable to unstable zones are identified.
Abstract: Carriers scattering in the inversion channel of n-
MOSFET dominates the drain current. This paper presents an effective
electron mobility model for the pocket implanted nano scale
n-MOSFET. The model is developed by using two linear pocket
profiles at the source and drain edges. The channel is divided into
three regions at source, drain and central part of the channel region.
The total number of inversion layer charges is found for these three
regions by numerical integration from source to drain ends and the
number of depletion layer charges is found by using the effective
doping concentration including pocket doping effects. These two
charges are then used to find the effective normal electric field,
which is used to find the effective mobility model incorporating the
three scattering mechanisms, such as, Coulomb, phonon and surface
roughness scatterings as well as the ballistic phenomena for the
pocket implanted nano-scale n-MOSFET. The simulation results show
that the derived mobility model produces the same results as found
in the literatures.
Abstract: Numerical integration of initial boundary problem for advection equation in 3 ℜ is considered. The method used is
conditionally stable semi-Lagrangian advection scheme with high order interpolation on unstructured mesh. In order to increase time step integration the BFECC method with limiter TVD correction is used. The method is adopted on parallel graphic processor unit environment using NVIDIA CUDA and applied in Navier-Stokes solver. It is shown that the calculation on NVIDIA GeForce 8800
GPU is 184 times faster than on one processor AMDX2 4800+ CPU. The method is extended to the incompressible fluid dynamics solver. Flow over a Cylinder for 3D case is compared to the experimental data.
Abstract: The Goursat partial differential equation arises in
linear and non linear partial differential equations with mixed
derivatives. This equation is a second order hyperbolic partial
differential equation which occurs in various fields of study such as
in engineering, physics, and applied mathematics. There are many
approaches that have been suggested to approximate the solution of
the Goursat partial differential equation. However, all of the
suggested methods traditionally focused on numerical differentiation
approaches including forward and central differences in deriving the
scheme. An innovation has been done in deriving the Goursat partial
differential equation scheme which involves numerical integration
techniques. In this paper we have developed a new scheme to solve
the Goursat partial differential equation based on the Adomian
decomposition (ADM) and associated with Boole-s integration rule to
approximate the integration terms. The new scheme can easily be
applied to many linear and non linear Goursat partial differential
equations and is capable to reduce the size of computational work.
The accuracy of the results reveals the advantage of this new scheme
over existing numerical method.
Abstract: The present study concentrates on solving the along wind oscillation problem of a tall square building from first principles and across wind oscillation problem of the same from empirical relations obtained by experiments. The criterion for human comfort at the worst condition at the top floor of the building is being considered and a limiting value of height of a building for a given cross section is predicted. Numerical integrations are carried out as and when required. The results show severeness of across wind oscillations in comparison to along wind oscillation. The comfort criterion is combined with across wind oscillation results to determine the maximum allowable height of a building for a given square cross-section.
Abstract: Gas hydrates can agglomerate and block multiphase oil and gas pipelines when water is present at hydrate forming conditions. Using "Cold Flow Technology", the aim is to condition gas hydrates so that they can be transported as a slurry mixture without a risk of agglomeration. During the pipeline shut down however, hydrate particles may settle in bends and build hydrate plugs. An experimental setup has been designed and constructed to study the flow of such plugs at start up operations. Experiments have been performed using model fluid and model hydrate particles. The propagations of initial plugs in a bend were recorded with impedance probes along the pipe. The experimental results show a dispersion of the plug front. A peak in pressure drop was also recorded when the plugs were passing the bend. The evolutions of the plugs have been simulated by numerical integration of the incompressible mass balance equations, with an imposed mixture velocity. The slip between particles and carrier fluid has been calculated using a drag relation together with a particle-fluid force balance.
Abstract: Stick models are widely used in studying the
behaviour of straight as well as skew bridges and viaducts subjected
to earthquakes while carrying out preliminary studies. The
application of such models to highly curved bridges continues to
pose challenging problems. A viaduct proposed in the foothills of the
Himalayas in Northern India is chosen for the study. It is having 8
simply supported spans @ 30 m c/c. It is doubly curved in horizontal
plane with 20 m radius. It is inclined in vertical plane as well. The
superstructure consists of a box section. Three models have been
used: a conventional stick model, an improved stick model and a 3D
finite element model. The improved stick model is employed by
making use of body constraints in order to study its capabilities. The
first 8 frequencies are about 9.71% away in the latter two models.
Later the difference increases to 80% in 50th mode. The viaduct was
subjected to all three components of the El Centro earthquake of May
1940. The numerical integration was carried out using the Hilber-
Hughes-Taylor method as implemented in SAP2000. Axial forces
and moments in the bridge piers as well as lateral displacements at
the bearing levels are compared for the three models. The maximum
difference in the axial forces and bending moments and
displacements vary by 25% between the improved and finite element
model. Whereas, the maximum difference in the axial forces,
moments, and displacements in various sections vary by 35%
between the improved stick model and equivalent straight stick
model. The difference for torsional moment was as high as 75%. It is
concluded that the stick model with body constraints to model the
bearings and expansion joints is not desirable in very sharp S curved
viaducts even for preliminary analysis. This model can be used only
to determine first 10 frequency and mode shapes but not for member
forces. A 3D finite element analysis must be carried out for
meaningful results.
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 flow and heat transfer characteristics for natural
convection along an inclined plate in a saturated porous medium with
an applied magnetic field have been studied. The fluid viscosity has
been assumed to be an inverse function of temperature. Assuming
temperature vary as a power function of distance. The transformed
ordinary differential equations have solved by numerical integration
using Runge-Kutta method. The velocity and temperature profile
components on the plate are computed and discussed in detail for
various values of the variable viscosity parameter, inclination angle,
magnetic field parameter, and real constant (λ). The results have also
been interpreted with the aid of tables and graphs. The numerical
values of Nusselt number have been calculated for the mentioned
parameters.