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

In this paper, a numerical method for singular perturbation problems arising in chemical reactor theory for general singularly perturbed two point boundary value problems with boundary layer at one end(left or right) of the underlying interval is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1273-1287
• /
• 2008

In this paper, the numerical integration method for general singularly perturbed two point boundary value problems with mixed boundary conditions of both left and right end boundary layer is presented. The original second order differential equation is replaced by an approximate first order differential equation with a small deviating argument. By using the trapezoidal formula we obtain a three term recurrence relation, which is solved using Thomas Algorithm. To demonstrate the applicability of the method, we have solved four linear (two left and two right end boundary layer) and one nonlinear problems. From the results, it is observed that the present method approximates the exact or the asymptotic expansion solution very well.

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