Abstract: In this article, a simulation method called the Homotopy Perturbation Method (HPM) is employed in the steady flow of a Walter's B' fluid in a vertical channel with porous wall. We employed Homotopy Perturbation Method to derive solution of a nonlinear form of equation obtained from exerting similarity transforming to the ordinary differential equation gained from continuity and momentum equations of this kind of flow. The results obtained from the Homotopy Perturbation Method are then compared with those from the Runge–Kutta method in order to verify the accuracy of the proposed method. The results show that the Homotopy Perturbation Method can achieve good results in predicting the solution of such problems. Ultimately we use this solution to obtain the other terms of velocities and physical discussion about it.
Abstract: Both the minimum energy consumption and
smoothness, which is quantified as a function of jerk, are generally
needed in many dynamic systems such as the automobile and the
pick-and-place robot manipulator that handles fragile equipments.
Nevertheless, many researchers come up with either solely
concerning on the minimum energy consumption or minimum jerk
trajectory. This research paper proposes a simple yet very interesting
relationship between the minimum direct and indirect jerks
approaches in designing the time-dependent system yielding an
alternative optimal solution. Extremal solutions for the cost functions
of direct and indirect jerks are found using the dynamic optimization
methods together with the numerical approximation. This is to allow
us to simulate and compare visually and statistically the time history
of control inputs employed by minimum direct and indirect jerk
designs. By considering minimum indirect jerk problem, the
numerical solution becomes much easier and yields to the similar
results as minimum direct jerk problem.
Abstract: A Picard-Newton iteration method is studied to accelerate the numerical solution procedure of a class of two-dimensional nonlinear coupled parabolic-hyperbolic system. The Picard-Newton iteration is designed by adding higher-order terms of small quantity to an existing Picard iteration. The discrete functional analysis and inductive hypothesis reasoning techniques are used to overcome difficulties coming from nonlinearity and coupling, and theoretical analysis is made for the convergence and approximation properties of the iteration scheme. The Picard-Newton iteration has a quadratic convergent ratio, and its solution has second order spatial approximation and first order temporal approximation to the exact solution of the original problem. Numerical tests verify the results of the theoretical analysis, and show the Picard-Newton iteration is more efficient than the Picard iteration.
Abstract: The effect of a time delay on the transmission on
dengue fever is studied. The time delay is due to the presence of an
incubation period for the dengue virus to develop in the mosquito
before the mosquito becomes infectious. The conditions for the
existence of a Hopf bifurcation to limit cycle behavior are
established. The conditions are different from the usual one and they
are based on whether a particular third degree polynomial has
positive real roots. A theorem for determining whether for a given
set of parameter values, a critical delay time exist is given. It is
found that for a set of realistic values of the parameters in the model,
a Hopf bifurcation can not occur. For a set of unrealistic values of
some of the parameters, it is shown that a Hopf bifurcation can occur.
Numerical solutions using this last set show the trajectory of two of
the variables making a transition from a spiraling orbit to a limit
cycle orbit.
Abstract: Beginning from the creator of integro-differential
equations Volterra, many scientists have investigated these
equations. Classic method for solving integro-differential
equations is the quadratures method that is successfully applied up
today. Unlike these methods, Makroglou applied hybrid methods
that are modified and generalized in this paper and applied to the
numerical solution of Volterra integro-differential equations. The
way for defining the coefficients of the suggested method is also
given.
Abstract: Solution of some practical problems is reduced to the
solution of the integro-differential equations. But for the numerical
solution of such equations basically quadrature methods or its
combination with multistep or one-step methods are used. The
quadrature methods basically is applied to calculation of the integral
participating in right hand side of integro-differential equations. As
this integral is of Volterra type, it is obvious that at replacement with
its integrated sum the upper limit of the sum depends on a current
point in which values of the integral are defined. Thus we receive the
integrated sum with variable boundary, to work with is hardly.
Therefore multistep method with the constant coefficients, which is
free from noted lack and gives the way for finding it-s coefficients is
present.
Abstract: As we know, most differential equations concerning
physical phenomenon could not be solved by analytical method. Even if we use Series Method, some times we need an appropriate change of variable, and even when we can, their closed form solution may be
so complicated that using it to obtain an image or to examine the structure of the system is impossible. For example, if we consider Schrodinger equation, i.e.,
We come to a three-term recursion relations, which work with it takes, at least, a little bit time to get a series solution[6]. For this
reason we use a change of variable such as or when we consider the orbital angular momentum[1], it will be
necessary to solve. As we can observe, working with this equation is tedious. In this paper, after introducing Clenshaw method, which is a kind of Spectral method, we try to solve some of such equations.
Abstract: In the present analysis an unsteady laminar
forced convection water boundary layer flow is considered.
The fluid properties such as viscosity and Prandtl number are
taken as variables such that those are inversely proportional to
temperature. By using quasi-linearization technique the nonlinear
coupled partial differential equations are linearized and
the numerical solutions are obtained by using implicit finite
difference scheme with the appropriate selection of step sizes.
Non-similar solutions have been obtained from the starting
point of the stream-wise coordinate to the point where skin
friction value vanishes. The effect non-uniform mass transfer
along the surface of the cylinder through slot is studied on the
skin friction and heat transfer coefficients.
Abstract: A numerical solution of the initial boundary value
problem of the suspended string vibrating equation with the
particular nonlinear damping term based on the finite difference
scheme is presented in this paper. The investigation of how the
second and third power terms of the nonlinear term affect the
vibration characteristic. We compare the vibration amplitude as a
result of the third power nonlinear damping with the second power
obtained from previous report provided that the same initial shape
and initial velocities are assumed. The comparison results show that
the vibration amplitude is inversely proportional to the coefficient of
the damping term for the third power nonlinear damping case, while
the vibration amplitude is proportional to the coefficient of the
damping term in the second power nonlinear damping case.
Abstract: This manuscript presents a method for the numerical solution of the Cauchy type singular integral equations of the first kind, over a finite segment which is bounded at the end points of the finite segment. The Chebyshev polynomials of the second kind with the corresponding weight function have been used to approximate the density function. The force function is approximated by using the Chebyshev polynomials of the first kind. It is shown that the numerical solution of characteristic singular integral equation is identical with the exact solution, when the force function is a cubic function. Moreover, it also shown that this numerical method gives exact solution for other singular integral equations with degenerate kernels.
Abstract: A dual-reciprocity boundary element method is presented
for the numerical solution of a class of axisymmetric elastodynamic
problems. The domain integrals that arise in the integrodifferential
formulation are converted to line integrals by using the
dual-reciprocity method together suitably constructed interpolating
functions. The second order time derivatives of the displacement
in the governing partial differential equations are suppressed by
using Laplace transformation. In the Laplace transform domain, the
problem under consideration is eventually reduced to solving a system
of linear algebraic equations. Once the linear algebraic equations are
solved, the displacement and stress fields in the physical domain can
be recovered by using a numerical technique for inverting Laplace
transforms.
Abstract: This work examines thermal convection in two porous
layers. Flow in the upper layer is governed by Brinkman-s equations
model and in the lower layer is governed by Darcy-s model.
Legendre polynomials are used to obtain numerical solution when the
lower layer is heated from below.
Abstract: In this paper we investigate numerically positive solutions of the equation -Δu = λuq+up with Dirichlet boundary condition in a boundary domain ╬® for λ > 0 and 0 < q < 1 < p < 2*, we will compute and visualize the range of λ, this problem achieves a numerical solution.
Abstract: In this paper, a numerical solution based on sinc
functions is used for finding the solution of boundary value problems
which arise from the problems of calculus of variations. This
approximation reduce the problems to an explicit system of algebraic
equations. Some numerical examples are also given to illustrate the
accuracy and applicability of the presented method.
Abstract: This article is devoted to the numerical solution of
large-scale quadratic eigenvalue problems. Such problems arise in
a wide variety of applications, such as the dynamic analysis of
structural mechanical systems, acoustic systems, fluid mechanics,
and signal processing. We first introduce a generalized second-order
Krylov subspace based on a pair of square matrices and two initial
vectors and present a generalized second-order Arnoldi process for
constructing an orthonormal basis of the generalized second-order
Krylov subspace. Then, by using the projection technique and the
refined projection technique, we propose a restarted generalized
second-order Arnoldi method and a restarted refined generalized
second-order Arnoldi method for computing some eigenpairs of largescale
quadratic eigenvalue problems. Some theoretical results are also
presented. Some numerical examples are presented to illustrate the
effectiveness of the proposed methods.
Abstract: Semilinear elliptic equations are ubiquitous in natural sciences. They give rise to a variety of important phenomena in quantum mechanics, nonlinear optics, astrophysics, etc because they have rich multiple solutions. But the nontrivial solutions of semilinear equations are hard to be solved for the lack of stabilities, such as Lane-Emden equation, Henon equation and Chandrasekhar equation. In this paper, bifurcation method is applied to solving semilinear elliptic equations which are with homogeneous Dirichlet boundary conditions in 2D. Using this method, nontrivial numerical solutions will be computed and visualized in many different domains (such as square, disk, annulus, dumbbell, etc).
Abstract: This paper study the behavior of the solution at the crack edges for an elliptical crack with developing cusps, Ω in the plane elasticity subjected to shear loading. The problem of finding the resulting shear stress can be formulated as a hypersingular integral equation over Ω and it is then transformed into a similar equation over a circular region, D, using conformal mapping. An appropriate collocation points are chosen on the region D to reduce the hypersingular integral equation into a system of linear equations with (2N+1)(N+1) unknown coefficients, which will later be used in the determination of shear stress intensity factors and maximum shear stress intensity. Numerical solution for the considered problem are compared with the existing asymptotic solution, and displayed graphically. Our results give a very good agreement to the existing asymptotic solutions.
Abstract: Leptospirosis occurs worldwide (except the
poles of the earth), urban and rural areas, developed and
developing countries, especially in Thailand. It can be
transmitted to the human by rats through direct and indirect
ways. Human can be infected by either touching the infected rats
or contacting with water, soil containing urine from the infected
rats through skin, eyes and nose. The data of the people who
are infected with this disease indicates that most of the
patients are adults. The transmission of this disease is studied
through mathematical model. The population is separated into human
and rat. The human is divided into two classes, namely juvenile
and adult. The model equation is constructed for each class. The
standard dynamical modeling method is then used for
analyzing the behaviours of solutions. In addition, the
conditions of the parameters for the disease free and endemic
states are obtained. Numerical solutions are shown to support the
theoretical predictions. The results of this study guide the way to
decrease the disease outbreak.
Abstract: The reliable results of an insulated oval duct
considering heat radiation are obtained basing on accurate oval
perimeter obtained by integral method as well as one-dimensional
Plane Wedge Thermal Resistance (PWTR) model. This is an extension
study of former paper of insulated oval duct neglecting heat radiation.
It is found that in the practical situations with long-short-axes ratio a/b
4.5% while t/R2
Abstract: In this paper developed and realized absolutely new
algorithm for solving three-dimensional Poisson equation. This
equation used in research of turbulent mixing, computational fluid
dynamics, atmospheric front, and ocean flows and so on. Moreover in
the view of rising productivity of difficult calculation there was
applied the most up-to-date and the most effective parallel
programming technology - MPI in combination with OpenMP
direction, that allows to realize problems with very large data
content. Resulted products can be used in solving of important
applications and fundamental problems in mathematics and physics.