The parallelization of numerical codes on unstructured grids still is a demanding programming task, which is often performed in an application-speciic manner. However, many of the algorithms used in this eld exhibit strong structural similarities. A concept exploiting the latter is developed in this paper. Furthermore, it will be shown how the abstract concepts can be implemented in a reusable ...

Ordinary differential equations are ubiquitous in science and engineering: in geometry and mechanics from the first examples onwards (Newton, Leibniz, Euler, Lagrange), in chemical reaction kinetics, molecular dynamics, electronic circuits, population dynamics, and many more application areas. They also arise, after semidiscretization in space, in the numerical treatment of time-dependent parti...

A new algorithm is proposed for solving the inverse Sturm–Liouville problem of reconstructing a symmetric potential from eigenvalues. It uses Numerov’s method instead of the second order method of the related algorithm of Fabiano, Knobel and Lowe. An extension by Andrew and Paine of the asymptotic correction technique of Paine, de Hoog and Anderssen is the key to the success of the new algorith...

This work presents two high-order, semi-discrete, central-upwind schemes for computing approximate solutions of 1D systems of conservation laws. We propose a central weighted essentially non-oscillatory (CWENO) reconstruction, also we apply a fourth-order reconstruction proposed by Peer et al., and afterwards, we combine these reconstructions with a semi-discrete central-upwind numerical flux ...

Many time-varying phenomena of various fields in science and engineering can be modeled as a stochastic differential equations, so investigation of conditions for existence of solution and obtain the analytical and numerical solutions of them are important. In this paper, the Adomian decomposition method for solution of the stochastic differential equations are improved.  Uniqueness and converg...

in the present article, a numerical method is proposed for the numerical solution of thekdv equation by using a new approach by combining cubic b-spline functions. in this paper we convert the kdv equation to system of two equations. the method is shown to be unconditionally stable using von-neumann technique. to test accuracy the error norms2l, ∞l are computed. three invariants of motion are p...

the aim of this paper is to present a numerical method for singularly perturbed convection-diffusion problems with a delay. the method is a combination of the asymptotic expansion technique and the reproducing kernel method (rkm). first an asymptotic expansion for the solution of the given singularly perturbed delayed boundary value problem is constructed. then the reduced regular delayed diffe...

in this paper , the quintic b-spline collocation scheme is employed to approximate numerical solution of the kdv-like rosenau equation . this scheme is based on the crank-nicolson formulation for time integration and quintic b-spline functions for space integration . the unconditional stability of the present method is proved using von- neumann approach . since we do not know the exact solution...

substrates' concentration profile was studied in a porous matrix containing immobilized amyloglucosidase for glucose production. this analysis was performed by using an analytical method called least square method, and the results were compared with numerical solution. effects of effective diffusivity, michael's constant, maximum reaction rate and initial substrate concentration were studied on...