• Journal of applied mathematics & informatics
• /
• v.27 no.1_2
• /
• pp.441-452
• /
• 2009

In this paper, a fifth order numerical method is presented for solving singularly perturbed differential-difference equations with negative shift. In recent papers the term negative shift has been using for delay. Similar boundary value problems are associated with expected first exit time problem of the membrane, potential in models for neuron and in variational problems in control theory. In the numerical treatment for such type of boundary value problems, first we use Taylor approximation to tackle terms containing small shifts which converts it to a boundary value problem for singularly perturbed differential equation. 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 and is solved using the boundary conditions. Several numerical examples are solved and compared with exact solution. It is observed that present method approximates the exact solution very well.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.331-337
• /
• 2018

This paper presents dynamic algorithm of the MLS(moving least squares) difference method using first order differential Approximation. The governing equations are only discretized by the first order MLS derivative approximation. The system equation consists of an assembly of the approximate function, so the shape of system equation is similar to FEM(finite element method). The CDM(central difference method) is used for time integration of dynamic equilibrium equation. The natural frequency analyses of the MLS difference method and FEM are performed, and two analysis results are compared. Also, the accuracy of the proposed numerical method is verified by displaying the dynamic analysis results together with the results by the existing second order differential approximation. In the process of assembling the first order MLS derivative approximation, the oscillation error was suppressed and the stress distribution was interpreted as relatively uniform.

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

• Kyungpook Mathematical Journal
• /
• v.56 no.1
• /
• pp.221-233
• /
• 2016

The present paper proposes a new monotone iteration method for existence as well as approximation of the coupled solutions for a coupled periodic boundary value problem of first order ordinary nonlinear differential equations. A new hybrid coupled fixed point theorem involving the Dhage iteration principle is proved in a partially ordered normed linear space and applied to the coupled periodic boundary value problems for proving the main existence and approximation results of this paper. An algorithm for the coupled solutions is developed and it is shown that the sequences of successive approximations defined in a certain way converge monotonically to the coupled solutions of the related differential equations under some suitable mixed hybrid conditions. A numerical example is also indicated to illustrate the abstract theory developed in the paper.

• Honam Mathematical Journal
• /
• v.37 no.3
• /
• pp.299-315
• /
• 2015

In this paper, we study the first-order system least-squares (FOSLS) spectral method for parabolic partial differential equations. There were lots of least-squares approaches to solve elliptic partial differential equations using finite element approximation. Also, some approaches using spectral methods have been studied in recent. In order to solve the parabolic partial differential equations in parallel, we consider a parallel numerical method based on a hybrid method of the frequency-domain method and first-order system least-squares method. First, we transform the parabolic problem in the space-time domain to the elliptic problems in the space-frequency domain. Second, we solve each elliptic problem in parallel for some frequencies using the first-order system least-squares method. And then we take the discrete inverse Fourier transforms in order to obtain the approximate solution in the space-time domain. We will introduce such a hybrid method and then present a numerical experiment.

• 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.26 no.1_2
• /
• pp.135-149
• /
• 2008

Many deterministic systems are described by Ordinary differential equations and can often be improved by including stochastic effects, but numerical methods for solving stochastic differential equations(SDEs) are required, and work in this area is far less advanced than for deterministic differential equations. In this paper,first we follow [7] to describe Runge-Kutta methods with order 2 from Taylor approximations in the weak sense and present two well known Runge-Kutta methods, RK2-TO and RK2-PL. Then we obtain a new 3-stage explicit Runge-Kutta with order 2 in weak sense and compare the numerical results among these three methods.

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

• 전기의세계
• /
• v.22 no.1
• /
• pp.28-34
• /
• 1973

This paper treats with the sub-optimal control problem of noninear systems by approximation method. This method involves the approximation by linearization which provides the sub-optimal solution of non-linear control problems. The result of this work shows that, in the problem in which the controlled plant is characterized by an ordinary differential equation of first order, the solution obtained by this method coincides with the exact solution of problem. In of case of the second or higher order systems, it is proved analytically that this method of linearization produces the sub-optimal solution of the given problem. It is also shown that the sub-optimality of solution by the method can be evaluated by introducing the upper and lower bounded performance indices. Discussion is made on the procedure with some illustrative examples whose performance indices are given in the quadratic forms.

• Structural Engineering and Mechanics
• /
• v.22 no.3
• /
• pp.331-349
• /
• 2006

The stresses and deflections in a laminated rectangular plate under thermal vibration are determined by using the moving least squares differential quadrature (MLSDQ) method based on the first order shear deformation theory. The weighting coefficients used in MLSDQ approximation are obtained through a fast computation of the MLS shape functions and their partial derivatives. By using this method, the governing differential equations are transformed into sets of linear homogeneous algebraic equations in terms of the displacement components at each discrete point. Boundary conditions are implemented through discrete grid points by constraining displacements, bending moments and rotations of the plate. Solving this set of algebraic equations yields the displacement components. Then substituting these displacements into the constitutive equation, we obtain the stresses. The approximate solutions for stress and deflection of laminated plate with cross layer under thermal load are obtained. Numerical results show that the MLSDQ method provides rapidly convergent and accurate solutions for calculating the stresses and deflections in a multi-layered plate of cross ply laminate subjected to thermal vibration of sinusoidal temperature including shear deformation with a few grid points.