• Journal of applied mathematics & informatics
• /
• 제27권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.

• 한국전산구조공학회논문집
• /
• 제31권6호
• /
• pp.331-337
• /
• 2018

본 논문은 MLS(moving least squares) 차분법의 1차 미분 근사함수를 바탕으로 시간에 따른 수치해석이 가능한 해석기법을 제시한다. 오직 1차 미분 근사함수로만 지배방정식을 이산화했으며, 근사함수를 조립하는 형태로 전체 시스템 방정식을 구성하여 차분법으로 이산화된 운동방정식이 유한요소법(finite element method)과 유사한 모습을 갖게 되었다. 운동방정식을 시간적분하기 위해서 중앙차분법(central difference method)을 사용하였다. 유한요소 알고리즘을 통해서 MLS 차분법과 유한요소법의 고유진동 해석을 수행하였으며, 두 해석결과를 비교하였다. 또한, 동적해석 결과를 기존의 2차 미분 근사함수를 활용한 해석결과와 함께 도시함으로써 제안된 수치기법의 정확성을 검증하였다. 1차 미분 근사함수를 조립하는 과정에서 해석결과의 떨림현상이 억제되었으며 상대적으로 균일한 응력분포를 구할 수 있었다.

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

• Kyungpook Mathematical Journal
• /
• 제56권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.

• 호남수학학술지
• /
• 제37권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
• /
• 제26권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
• /
• 제26권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
• /
• 제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.

• 전기의세계
• /
• 제22권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
• /
• 제22권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.