• Title/Summary/Keyword: Singularly perturbed problem

Search Result 42, Processing Time 0.044 seconds

Higher Order Uniformly Convergent Numerical Scheme for Singularly Perturbed Reaction-Diffusion Problems

  • Anilay, Worku Tilahun;Duressa, Gemechis File;Woldaregay, Mesfin Mekuria
    • Kyungpook Mathematical Journal
    • /
    • v.61 no.3
    • /
    • pp.591-612
    • /
    • 2021
  • In this paper, a uniformly convergent numerical scheme is designed for solving singularly perturbed reaction-diffusion problems. The problem is converted to an equivalent weak form and then a Galerkin finite element method is used on a piecewise uniform Shishkin mesh with linear basis functions. The convergence of the developed scheme is proved and it is shown to be almost fourth order uniformly convergent in the maximum norm. To exhibit the applicability of the scheme, model examples are considered and solved for different values of a singular perturbation parameter ε and mesh elements. The proposed scheme approximates the exact solution very well.

UNIFORMLY CONVERGENT NUMERICAL SCHEME FOR A SINGULARLY PERTURBED DIFFERENTIAL-DIFFERENCE EQUATIONS ARISING IN COMPUTATIONAL NEUROSCIENCE

  • DABA, IMIRU TAKELE;DURESSA, GEMECHIS FILE
    • Journal of applied mathematics & informatics
    • /
    • v.39 no.5_6
    • /
    • pp.655-676
    • /
    • 2021
  • A parameter uniform numerical scheme is proposed for solving singularly perturbed parabolic partial differential-difference convection-diffusion equations with a small delay and advance parameters in reaction terms and spatial variable. Taylor's series expansion is applied to approximate problems with the delay and advance terms. The resulting singularly perturbed parabolic convection-diffusion equation is solved by utilizing the implicit Euler method for the temporal discretization and finite difference method for the spatial discretization on a uniform mesh. The proposed numerical scheme is shown to be an ε-uniformly convergent accurate of the first order in time and second-order in space directions. The efficiency of the scheme is proved by some numerical experiments and by comparing the results with other results. It has been found that the proposed numerical scheme gives a more accurate approximate solution than some available numerical methods in the literature.

A NUMERICAL METHOD FOR SINGULARLY PERTURBED SYSTEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS OF CONVECTION DIFFUSION TYPE WITH A DISCONTINUOUS SOURCE TERM

  • Tamilselvan, A.;Ramanujam, N.
    • Journal of applied mathematics & informatics
    • /
    • v.27 no.5_6
    • /
    • pp.1279-1292
    • /
    • 2009
  • In this paper, a numerical method that uses standard finite difference scheme defined on Shishkin mesh for a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a discontinuous source term is presented. An error estimate is derived to show that the method is uniformly convergent with respect to the singular perturbation parameter. Numerical results are presented to illustrate the theoretical results.

  • PDF

AN ASYMPTOTIC FINITE ELEMENT METHOD FOR SINGULARLY PERTURBED HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS OF CONVECTION-DIFFUSION TYPE WITH DISCONTINUOUS SOURCE TERM

  • Babu, A. Ramesh;Ramanujam, N.
    • Journal of applied mathematics & informatics
    • /
    • v.26 no.5_6
    • /
    • pp.1057-1069
    • /
    • 2008
  • We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

  • PDF

AN SDFEM FOR A CONVECTION-DIFFUSION PROBLEM WITH NEUMANN BOUNDARY CONDITION AND DISCONTINUOUS SOURCE TERM

  • Babu, A. Ramesh;Ramanujam, N.
    • Journal of applied mathematics & informatics
    • /
    • v.28 no.1_2
    • /
    • pp.31-48
    • /
    • 2010
  • In this article, we consider singularly perturbed Boundary Value Problems(BVPs) for second order Ordinary Differential Equations (ODEs) with Neumann boundary condition and discontinuous source term. A parameter-uniform error bound for the solution is established using the Streamline-Diffusion Finite Element Method (SDFEM) on a piecewise uniform meshes. We prove that the method is almost second order of convergence in the maximum norm, independently of the perturbation parameter. Further we derive superconvergence results for scaled derivatives of solution of the same problem. Numerical results are provided to substantiate the theoretical results.

AN ASYMPTOTIC INITIAL VALUE METHOD FOR SECOND ORDER SINGULAR PERTURBATION PROBLEMS OF CONVECTION-DIFFUSION TYPE WITH A DISCONTINUOUS SOURCE TERM

  • Valanarasu, T.;Ramanujam, N.
    • Journal of applied mathematics & informatics
    • /
    • v.23 no.1_2
    • /
    • pp.141-152
    • /
    • 2007
  • In this paper a numerical method is presented to solve singularly perturbed two points boundary value problems for second order ordinary differential equations consisting a discontinuous source term. First, in this method, an asymptotic expansion approximation of the solution of the boundary value problem is constructed using the basic ideas of a well known perturbation method WKB. Then some initial value problems and terminal value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial value problems are happened to be singularly perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples provided to illustrate the method.

HYBRID DIFFERENCE SCHEMES FOR SINGULARLY PERTURBED PROBLEM OF MIXED TYPE WITH DISCONTINUOUS SOURCE TERM

  • Priyadharshini, R. Mythili;Ramanujam, N.;Valanarasu, T.
    • Journal of applied mathematics & informatics
    • /
    • v.28 no.5_6
    • /
    • pp.1035-1054
    • /
    • 2010
  • We consider a mixed type singularly perturbed one dimensional elliptic problem with discontinuous source term. The domain under consideration is partitioned into two subdomains. A convection-diffusion and a reaction-diffusion type equations are posed on the first and second subdomains respectively. Two hybrid difference schemes on Shishkin mesh are constructed and we prove that the schemes are almost second order convergence in the maximum norm independent of the diffusion parameter. Error bounds for the numerical solution and its numerical derivative are established. Numerical results are presented which support the theoretical results.

SPLINE DIFFERENCE SCHEME FOR TWO-PARAMETER SINGULARLY PERTURBED PARTIAL DIFFERENTIAL EQUATIONS

  • Zahra, W.K.;El-Azab, M.S.;Mhlawy, Ashraf M. El
    • Journal of applied mathematics & informatics
    • /
    • v.32 no.1_2
    • /
    • pp.185-201
    • /
    • 2014
  • In this paper, we construct a numerical method to solve singularly perturbed one-dimensional parabolic convection-diffusion problems. We use Euler method with uniform step size for temporal discretization and exponential-spline scheme on spatial uniform mesh of Shishkin type for full discretization. We show that the resulting method is uniformly convergent with respect to diffusion parameter. An extensive amount of analysis has been carried out to prove the uniform convergence with respect to the singular perturbation parameter. The obtained numerical results show that the method is efficient, stable and reliable for solving convection-diffusion problem accurately even involving diffusion parameter.

A FIFTH ORDER NUMERICAL METHOD FOR SINGULAR PERTURBATION PROBLEMS

  • Chakravarthy, P. Pramod;Phaneendra, K.;Reddy, Y.N.
    • Journal of applied mathematics & informatics
    • /
    • v.26 no.3_4
    • /
    • pp.689-706
    • /
    • 2008
  • In this paper, a fifth order numerical method is presented for solving singularly perturbed two point boundary value problems with a boundary layer at one end point. The two point boundary value problem is transformed into general first order ordinary differential equation system. A discrete approximation of a fifth order compact difference scheme is presented for the first order system. An asymptotically equivalent first order equation of the original singularly perturbed two point boundary value problem is obtained from the theory of singular perturbations. It is used in the fifth order compact difference scheme to get a two term recurrence relation and is solved. Several linear and non-linear singular perturbation problems have been solved and the numerical results are presented to support the theory. It is observed that the present method approximates the exact solution very well.

  • PDF

Robust Control for Singularly Perturbed Uncertain Systems with State Constraints

  • Lee, Sang-Yup;Kim, Eung-Ju;Kim, Beom-Soo;Lim, Myo-Taeg
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 2001.10a
    • /
    • pp.34.1-34
    • /
    • 2001
  • We deal with robust control problem for singularly perturbed linear systems with norm-bounded structured uncertainty under state constraints. We assume that the norm-bounded uncertainty is composed of repeated scalar-block and full-block forms. In the structured uncertainty, repeated scalar block forms account for uncertain physical parameter value and full-block forms may be some unknown nonlinear dynamics. In order deal with uncertainty and state constraints, we use LMI(Linear Matrix Inequality). The original problem is decomposed into two well behaved reduced order problems. Shinc two LMI problems are completely independent, each solution can be computed simultaneously and work in parallel.

  • PDF