Abstract: In this work, we present an efficient approach for
solving variable-order time-fractional partial differential equations,
which are based on Legendre and Laguerre polynomials. First, we
introduced the pseudo-operational matrices of integer and variable
fractional order of integration by use of some properties of
Riemann-Liouville fractional integral. Then, applied together with
collocation method and Legendre-Laguerre functions for solving
variable-order time-fractional partial differential equations. Also, an
estimation of the error is presented. At last, we investigate numerical
examples which arise in physics to demonstrate the accuracy of the
present method. In comparison results obtained by the present method
with the exact solution and the other methods reveals that the method
is very effective.
Abstract: The magnetohydrodynamic (MHD) Falkner-Skan
equations appear in study of laminar boundary layers flow over
a wedge in presence of a transverse magnetic field. The partial
differential equations of boundary layer problems in presence of
a transverse magnetic field are reduced to MHD Falkner-Skan
equation by similarity solution methods. This is a nonlinear ordinary
differential equation. In this paper, we solve this equation via
spectral collocation method based on Bessel functions of the first
kind. In this approach, we reduce the solution of the nonlinear
MHD Falkner-Skan equation to a solution of a nonlinear algebraic
equations system. Then, the resulting system is solved by Newton
method. We discuss obtained solution by studying the behavior
of boundary layer flow in terms of skin friction, velocity, various
amounts of magnetic field and angle of wedge. Finally, the results
are compared with other methods mentioned in literature. We can
conclude that the presented method has better accuracy than others.
Abstract: In this paper, we introduce a generalized Chebyshev
collocation method (GCCM) based on the generalized Chebyshev
polynomials for solving stiff systems. For employing a technique
of the embedded Runge-Kutta method used in explicit schemes, the
property of the generalized Chebyshev polynomials is used, in which
the nodes for the higher degree polynomial are overlapped with those
for the lower degree polynomial. The constructed algorithm controls
both the error and the time step size simultaneously and further
the errors at each integration step are embedded in the algorithm
itself, which provides the efficiency of the computational cost. For
the assessment of the effectiveness, numerical results obtained by the
proposed method and the Radau IIA are presented and compared.
Abstract: In this paper, numerical solution of system of
Fredholm and Volterra integral equations by means of the Spline
collocation method is considered. This approximation reduces the
system of integral equations to an explicit system of algebraic
equations. The solution is collocated by cubic B-spline and the
integrand is approximated by the Newton-Cotes formula. The error
analysis of proposed numerical method is studied theoretically. The
results are compared with the results obtained by other methods to
illustrate the accuracy and the implementation of our method.
Abstract: The generalized wave equation models various
problems in sciences and engineering. In this paper, a new three-time
level implicit approach based on cubic trigonometric B-spline for the
approximate solution of wave equation is developed. The usual finite
difference approach is used to discretize the time derivative while
cubic trigonometric B-spline is applied as an interpolating function in
the space dimension. Von Neumann stability analysis is used to
analyze the proposed method. Two problems are discussed to exhibit
the feasibility and capability of the method. The absolute errors and
maximum error are computed to assess the performance of the
proposed method. The results were found to be in good agreement
with known solutions and with existing schemes in literature.
Abstract: In this paper, the numerical solution of optimal control problem (OCP) for systems governed by Volterra integro-differential
(VID) equation is considered. The method is developed by means
of the Legendre wavelet approximation and collocation method. The
properties of Legendre wavelet together with Gaussian integration
method are utilized to reduce the problem to the solution of nonlinear
programming one. Some numerical examples are given to confirm the
accuracy and ease of implementation of the method.
Abstract: In this paper, we propose Hermite collocation method for solving Thomas-Fermi equation that is nonlinear ordinary differential equation on semi-infinite interval. This method reduces the solution of this problem to the solution of a system of algebraic equations. We also present the comparison of this work with solution of other methods that shows the present solution is more accurate and faster convergence in this problem.
Abstract: In this paper the Kuramoto-Sivashinsky equation is solved numerically by collocation method. The solution is approximated as a linear combination of septic B-spline functions. Applying the Von-Neumann stability analysis technique, we show that the method is unconditionally stable. The method is applied on some test examples, and the numerical results have been compared with the exact solutions. The global relative error and L∞ in the solutions show the efficiency of the method computationally.
Abstract: In this paper, numerical solutions of the nonlinear Benjamin-Bona-Mahony-Burgers (BBMB) equation are obtained by a method based on collocation of cubic B-splines. Applying the Von-Neumann stability analysis, the proposed method is shown to be unconditionally stable. The method is applied on some test examples, and the numerical results have been compared with the exact solutions. The L∞ and L2 in the solutions show the efficiency of the method computationally.
Abstract: A generalized Dirichlet to Neumann map is
one of the main aspects characterizing a recently introduced
method for analyzing linear elliptic PDEs, through which it
became possible to couple known and unknown components
of the solution on the boundary of the domain without
solving on its interior. For its numerical solution, a well conditioned
quadratically convergent sine-Collocation method
was developed, which yielded a linear system of equations
with the diagonal blocks of its associated coefficient matrix
being point diagonal. This structural property, among others,
initiated interest for the employment of iterative methods for
its solution. In this work we present a conclusive numerical
study for the behavior of classical (Jacobi and Gauss-Seidel)
and Krylov subspace (GMRES and Bi-CGSTAB) iterative
methods when they are applied for the solution of the Dirichlet
to Neumann map associated with the Laplace-s equation
on regular polygons with the same boundary conditions on
all edges.
Abstract: Sinc-collocation scheme is one of the new techniques
used in solving numerical problems involving integral equations. This
method has been shown to be a powerful numerical tool for finding
fast and accurate solutions. So, in this paper, some properties of the
Sinc-collocation method required for our subsequent development
are given and are utilized to reduce integral equation of the first
kind to some algebraic equations. Then convergence with exponential
rate is proved by a theorem to guarantee applicability of numerical
technique. Finally, numerical examples are included to demonstrate
the validity and applicability of the technique.
Abstract: In this paper we study some numerical methods to solve a model one-dimensional convection–diffusion equation. The semi-discretisation of the space variable results into a system of ordinary differential equations and the solution of the latter involves the evaluation of a matrix exponent. Since the calculation of this term is computationally expensive, we study some methods based on Krylov subspace and on Restrictive Taylor series approximation respectively. We also consider the Chebyshev Pseudospectral collocation method to do the spatial discretisation and we present the numerical solution obtained by these methods.
Abstract: This article presents a numerical study of the doublediffusive
mixed convection in a vertical channel filled with porous
medium by using non-equilibrium model. The flow is assumed
fully developed, uni-directional and steady state. The controlling
parameters are thermal Rayleigh number (RaT ), Darcy number (Da),
Forchheimer number (F), buoyancy ratio (N), inter phase heat transfer
coefficient (H), and porosity scaled thermal conductivity ratio
(γ). The Brinkman-extended non-Darcy model is considered. The
governing equations are solved by spectral collocation method. The
main emphasize is given on flow profiles as well as heat and solute
transfer rates, when two diffusive components in terms of buoyancy
ratio are in favor (against) of each other and solid matrix and fluid
are thermally non-equilibrium. The results show that, for aiding flow
(RaT = 1000), the heat transfer rate of fluid (Nuf ) increases upto a
certain value of H, beyond that decreases smoothly and converges
to a constant, whereas in case of opposing flow (RaT = -1000),
the result is same for N = 0 and 1. The variation of Nuf in (N,
Nuf )-plane shows sinusoidal pattern for RaT = -1000. For both cases
(aiding and opposing) the flow destabilize on increasing N by inviting
point of inflection or flow separation on the velocity profile. Overall,
the buoyancy force have significant impact on the non-Darcy mixed
convection under LTNE conditions.
Abstract: The problem addressed herein is the efficient management of the Grid/Cluster intense computation involved, when the preconditioned Bi-CGSTAB Krylov method is employed for the iterative solution of the large and sparse linear system arising from the discretization of the Modified Helmholtz-Dirichlet problem by the Hermite Collocation method. Taking advantage of the Collocation ma-trix's red-black ordered structure we organize efficiently the whole computation and map it on a pipeline architecture with master-slave communication. Implementation, through MPI programming tools, is realized on a SUN V240 cluster, inter-connected through a 100Mbps and 1Gbps ethernet network,and its performance is presented by speedup measurements included.
Abstract: Recently, it is found that telegraph equation is more suitable than ordinary diffusion equation in modelling reaction diffusion for such branches of sciences. In this paper, a numerical solution for the one-dimensional hyperbolic telegraph equation by using the collocation method using the septic splines is proposed. The scheme works in a similar fashion as finite difference methods. Test problems are used to validate our scheme by calculate L2-norm and L∞-norm. The accuracy of the presented method is demonstrated by two test problems. The numerical results are found to be in good agreement with the exact solutions.
Abstract: In this paper, the telegraph equation is solved numerically by cubic B-spline quasi-interpolation .We obtain the numerical scheme, by using the derivative of the quasi-interpolation to approximate the spatial derivative of the dependent variable and a low order forward difference to approximate the temporal derivative of the dependent variable. The advantage of the resulting scheme is that the algorithm is very simple so it is very easy to implement. The results of numerical experiments are presented, and are compared with analytical solutions by calculating errors L2 and L∞ norms to confirm the good accuracy of the presented scheme.
Abstract: In this paper, we explore the applicability of the Sinc-
Collocation method to a three-dimensional (3D) oceanography model.
The model describes a wind-driven current with depth-dependent
eddy viscosity in the complex-velocity system. In general, the
Sinc-based methods excel over other traditional numerical methods
due to their exponentially decaying errors, rapid convergence and
handling problems in the presence of singularities in end-points.
Together with these advantages, the Sinc-Collocation approach that
we utilize exploits first derivative interpolation, whose integration
is much less sensitive to numerical errors. We bring up several
model problems to prove the accuracy, stability, and computational
efficiency of the method. The approximate solutions determined by
the Sinc-Collocation technique are compared to exact solutions and
those obtained by the Sinc-Galerkin approach in earlier studies. Our
findings indicate that the Sinc-Collocation method outperforms other
Sinc-based methods in past studies.