• 제목/요약/키워드: Third order singularly perturbed problems

검색결과 3건 처리시간 0.016초

AN OVERLAPPING SCHWARZ METHOD FOR SINGULARLY PERTURBED THIRD ORDER CONVECTION-DIFFUSION TYPE

  • ROJA, J. CHRISTY;TAMILSELVAN, A.
    • Journal of applied mathematics & informatics
    • /
    • 제36권1_2호
    • /
    • pp.135-154
    • /
    • 2018
  • In this paper, an almost second order overlapping Schwarz method for singularly perturbed third order convection-diffusion type problem is constructed. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the combination of classical finite difference scheme and central finite difference scheme on a uniform mesh while on the non-layer region we use the midpoint difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. We proved that, when appropriate subdomains are used, the method produces convergence of second order. Furthermore, it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme are it reduce iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

NUMERICAL METHOD FOR SINGULARLY PERTURBED THIRD ORDER ORDINARY DIFFERENTIAL EQUATIONS OF REACTION-DIFFUSION TYPE

  • ROJA, J. CHRISTY;TAMILSELVAN, A.
    • Journal of applied mathematics & informatics
    • /
    • 제35권3_4호
    • /
    • pp.277-302
    • /
    • 2017
  • In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

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
    • /
    • 제26권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