• Journal of applied mathematics & informatics
• /
• v.35 no.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.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.6
• /
• pp.1542-1550
• /
• 1994

For a class of singularly perturbed uncertain system, an output feedback control law is designed. The controller structure is designed based on the uncertain reduced-order system, and the controller parameters are determined by information on the reduced-order and full-order systems. It has been shown that the reduces-order system with the designed controller possesses a stability property(specifically, a global uniform attractivity). Furthermore, the stability property of this control scheme is robust with respect to singular perturbation ; i.e., the full-order system, subject to the same controller, possesses the global uniform attractivity, provided the singular perturbation parameter $\mu<\mu^{*}$, where a threshold value $\mu^{*}$ can be computed from the information available on the full-order system.

• Journal of applied mathematics & informatics
• /
• v.33 no.5_6
• /
• pp.635-648
• /
• 2015

We consider a weakly coupled system of singularly perturbed convection-diffusion equations with discontinuous source term. The diffusion term of each equation is associated with a small positive parameter of different magnitude. Presence of discontinuity and different parameters creates boundary and weak interior layers that overlap and interact. A numerical method is constructed for this problem which involves an appropriate piecewise uniform Shishkin mesh. The numerical approximations are proved to converge to the continuous solutions uniformly with respect to the singular perturbation parameters. Numerical results are presented which illustrates the theoretical results.

• 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.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.109-130
• /
• 2010

In this paper, a numerical method for a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with the mixed type boundary conditions is presented. Parameter-uniform error bounds for the numerical solution and also to numerical derivative are established. Numerical results are provided to illustrate the theoretical results.

• Proceedings of the KIEE Conference
• /
• 1987.07b
• /
• pp.1225-1228
• /
• 1987

It is shown that the singularly perturbed continuous-time system is led to two different discrete versions according to slow or fast sampling rates. The design of stabilizing feedback control of singularly perturbed discrete-time stochastic system is decomposed into the design of slow and fast controllers, which is combined to form the composite control. Composite control law is derived for the case of both single rate measurement and multirate measurement.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1001-1015
• /
• 2009

In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.49-65
• /
• 2006

A robust numerical method for a singularly perturbed second-order ordinary differential equation having two parameters with a discontinuous source term is presented in this article. Theoretical bounds are derived for the derivatives of the solution and its smooth and singular components. An appropriate piecewise uniform mesh is constructed, and classical upwind finite difference schemes are used on this mesh to obtain the discrete system of equations. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are provided to illustrate the convergence of the numerical approximations.

• Nuclear Engineering and Technology
• /
• v.21 no.2
• /
• pp.99-110
• /
• 1989

The optimal control of xenon concentration in a nuclear reactor is posed as a linear quadratic regulator problem with state feedback control. Since it is not possible to measure the state variables such as xenon and iodine concentrations directly, implementation of the optimal state feedback control law requires estimation of the unmeasurable state variables. The estimation method used is based on the Luenberger observer. The set of the reactor kinetics equations is a stiff system. This singularly perturbed system arises from the interaction of slow dynamic modes (iodine and xenon concentrations) and fast dynamic modes (neutron flux, fuel and coolant temperatures). The singular perturbation technique is used to overcome this stiffness problem. The observer-based controller of the original system is effected by separate design of the observer and controller of the reduced subsystem and the fast subsystem. In particular, since in the reactor kinetics control problem analyzed in the study the fast mode dies out quickly, we need only design the observer for the reduced slow subsystem. The results of the test problems demonstrated that the state feedback control of the xenon oscillation can be accomplished efficiently and without sacrificing accuracy by using the observer combined with the singular perturbation method.