• 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 the Korean Mathematical Society
• /
• v.51 no.4
• /
• pp.679-702
• /
• 2014

Investigation of the numerical solution of singularly perturbed turning point problems dates back to late 1970s. However, due to the presence of layers, not many high order schemes could be developed to solve such problems. On the other hand, one could think of applying the convergence acceleration technique to improve the performance of existing numerical methods. However, that itself posed some challenges. To this end, we design and analyze a novel fitted operator finite difference method (FOFDM) to solve this type of problems. Then we develop a fitted mesh finite difference method (FMFDM). Our detailed convergence analysis shows that this FMFDM is robust with respect to the singular perturbation parameter. Then we investigate the effect of Richardson extrapolation on both of these methods. We observe that, the accuracy is improved in both cases whereas the rate of convergence depends on the particular scheme being used.

• Kyungpook Mathematical Journal
• /
• v.60 no.3
• /
• pp.629-645
• /
• 2020

A uniformly convergent numerical method is developed for solving singularly perturbed 1-D parabolic convection-diffusion problems. The developed method applies a non-standard finite difference method for the spatial derivative discretization and uses the implicit Runge-Kutta method for the semi-discrete scheme. The convergence of the method is analyzed, and it is shown to be first order convergent. To validate the applicability of the proposed method two model examples are considered and solved for different perturbation parameters and mesh sizes. The numerical and experimental results agree well with the theoretical findings.

• International Journal of Control, Automation, and Systems
• /
• v.6 no.1
• /
• pp.76-85
• /
• 2008

Motivated by the fact that in many industrial robots the joint velocity is estimated from position measurements, the trajectory tracking of robot manipulators with output feedback is addressed in this paper. The fact that robot actuators have limited power is also taken into account. Let us notice that few solutions for the torque-bounded output feedback tracking control problem have been proposed. In this paper we contribute to this subject by presenting a theoretical reexamination of a known controller, by using the theory of singularly perturbed systems. Motivated by this analysis, a redesign of that controller is introduced. As another contribution, we present an experimental evaluation in a two degrees-of-freedom revolute-joint direct-drive robot, confirming the practical feasibility of the proposed approach.

• Communications of the Korean Mathematical Society
• /
• v.37 no.3
• /
• pp.939-955
• /
• 2022

In this paper, we study a first-order non-linear singularly perturbed Volterra integro-differential equation (SPVIDE). We discretize the problem by a uniform difference scheme on a Bakhvalov-Shishkin mesh. The scheme is constructed by the method of integral identities with exponential basis functions and integral terms are handled with interpolating quadrature rules with remainder terms. An effective quasi-linearization technique is employed for the algorithm. We establish the error estimates and demonstrate that the scheme on Bakhvalov-Shishkin mesh is O(N^-1) uniformly convergent, where N is the mesh parameter. The numerical results on a couple of examples are also provided to confirm the theoretical analysis.

• Journal of applied mathematics & informatics
• /
• v.40 no.3_4
• /
• pp.491-505
• /
• 2022

This paper is concerned with singularly perturbed convection-diffusion parabolic partial differential equations which have time-delayed. We used the Crank-Nicolson(CN) scheme to build a fitted operator to solve the problem. The underling method's stability is investigated, and it is found to be unconditionally stable. We have shown graphically the unstableness of CN-scheme without fitting factor. The order of convergence of the present method is shown to be second order both in space and time in relation to the perturbation parameter. The efficiency of the scheme is demonstrated using model examples and the proposed technique is more accurate than the standard CN-method and some methods available in the literature, according to the findings.

• Interaction and multiscale mechanics
• /
• v.4 no.4
• /
• pp.291-311
• /
• 2011

To simulate the self excited torsional vibrations of rotating drill strings (DSs) in vertical bore-holes, the nonlinear wave models of homogeneous and sectional torsional pendulums are formulated. The stated problem is shown to be of singularly perturbed type because the coefficient appearing before the second derivative of the constitutive nonlinear differential equation is small. The diapasons ${\omega}_b\leq{\omega}\leq{\omega}_l$ of angular velocity ${\omega}$ of the DS rotation are found, where the torsional auto-oscillations (of limit cycles) of the DS bit are generated. The variation of the limit cycle states, i.e. birth (${\omega}={\omega}_b$), evolution (${\omega}_b<{\omega}<{\omega}_l$) and loss (${\omega}={\omega}_l$), with the increase in angular velocity ${\omega}$ is analyzed. It is observed that firstly, at birth state of bifurcation of the limit cycle, the auto-oscillation generated proceeds in the regime of fast and slow motions (multiscale motion) with very small amplitude and it has a relaxation mode with nearly discontinuous angular velocities of elastic twisting. The vibration amplitude increases as ${\omega}$ increases, and then it decreases as ${\omega}$ approaches ${\omega}_l$. Sectional drill strings are also considered, and the conditions of the solution at the point of the upper and lower section joints are deduced. Besides, the peculiarities of the auto-oscillations of the sectional DSs are discussed.

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

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.937-948
• /
• 2011

The method of Asymptotic Inner Boundary Condition for Singularly Perturbed Two-Point Boundary value Problems is presented. By using a terminal point, the original second order problem is divided in to two problems namely inner region and outer region problems. The original problem is replaced by an asymptotically equivalent first order problem and using the stretching transformation, the asymptotic inner condition in implicit form at the terminal point is determined from the reduced equation of the original second order problem. The modified inner region problem, using the transformation with implicit boundary conditions is solved and produces a condition for the outer region problem. We used Chawla's fourth order method to solve both the inner and outer region problems. The proposed method is iterative on the terminal point. Some numerical examples are solved to demonstrate the applicability of the method.

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

In this paper, we apply a new decomposition method for solving initial and boundary value problems, which is due to Noor and Noor [18]. The analytical results are calculated in terms of convergent series with easily computable components. The diagonal Pade approximants are applied to make the work more concise and for the better understanding of the solution behavior. The proposed technique is tested on boundary layer problem; Thomas-Fermi, Blasius and sixth-order singularly perturbed Boussinesq equations. Numerical results reveal the complete reliability of the suggested scheme. This new decomposition method can be viewed as an alternative of Adomian decomposition method and homotopy perturbation methods.