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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1436-1444
• /
• 1996

The problem ofinstability of laminated circular cylindrical shell under the action of torsio or lateral pressure is investigated. The analysis is based on the Sander's theory for finite deformations of thin shell. The buckling is elastic for thin compoisite shell nad the geometry is assumed to be free of initial imperfections. The equilibrium equations are obrained by usitn the p[erturbation technique. Solution procedure is based on the Galerkin mehtod. The computer program for numerical results is made for several stacking sequence, length-to-radius ratio, and radius-to-thickness ratio. The numerical results of buckling load are present.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2001.04a
• /
• pp.17-20
• /
• 2001

Combustion characteristics of KSR-III liquid rocket engine are investigated numerically in the standpoints of engine performance and acoustic instability. In the present calculation, engine performance for design and off-design conditions is estimated effectively with reasonable error. Numerical results of acoustic instability show that engine operation for the design condition has sufficient stability margin, but for a certain off-design condition, acoustic instability can be triggered by artificial pressure perturbation. The present results are in a good agreement with the available experimental results and can be adopted for the prediction of engine performance and stability, depending on the specific operating condition.

• Journal of Institute of Control, Robotics and Systems
• /
• v.14 no.12
• /
• pp.1270-1277
• /
• 2008

A Haar wavelet based numerical method for solving singularly perturbed linear time invariant system is presented in this paper. The reduced pure slow and pure fast subsystems are obtained by decoupling the singularly perturbed system and differential matrix equations are converted into algebraic Sylvester matrix equations via Haar wavelet technique. The operational matrix of integration and its inverse matrix are utilized to reduce the computational time to the solution of algebraic matrix equations. Finally a numerical example is given to demonstrate the validity and applicability of the proposed method.

• Journal of KSNVE
• /
• v.8 no.3
• /
• pp.408-419
• /
• 1998

Dynamic characteristics are investigated when a nonlinear system showing periodic and chaotic responses under harmonic excitation is exposed to random perturbation. Approach for both qulitative and quantitative analysis of the noise effect in a nonlinear system under harmonic excitation is presented. For the qualitative analysis, Lyapunov exponents are calculated and Poincar map is illustrated. For the quatitative analysis. Fokker-Planck equatin is solved numerical by means of a Path-integral solution procedure. Eigenvalue problem obtained from the numerical caculation is solved and the relation of eigenvalue, eigenvector and chaotic motion is investigated.

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

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.418-423
• /
• 2008

Among various sensitivity evaluation techniques, semi-analytical method is quite popular since this method is more advantageous than analytical method and global finite difference method. However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified for individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, the adjoint variable method combined with complex variable is proposed to obtain the shape and size sensitivity for structural optimization. The complex variable can present accurate results regardless of the perturbation size as well as easy to be implemented. Through a few numerical examples of the static problem for the structural sensitivity, the efficiency and reliability of the adjoint variable method combined with complex variable is demonstrated.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.585-592
• /
• 2003

Among various sensitivity evaluation techniques, semi-analytical method(SAM) is quite popular since this method is more advantageous than analytical method(AM) and global finite difference method(FDM). However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified fur individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, an iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes and the error of SAM caused by numerical difference scheme is alleviated by using a Von Neumann series approximation considering the higher order terms for the sensitivity derivatives.

• Nuclear Engineering and Technology
• /
• v.35 no.1
• /
• pp.45-54
• /
• 2003

A hyperbolic two-phase, two-fluid equation system developed in the previous work has been implemented in an existing nuclear safety analysis code, MARS. Although the implicit treatment of interfacial pressure force term introduced in momentum equation of the hyperbolic equation system is required to enhance the numerical stability, it is very difficult to implement in the code because it is not possible to maintain the existing numerical solution structure. As an alternative, two-step approach with stabilizer momentum equations has been selected. The results of a linear stability analysis by Von-Neumann method show the equivalent stability improvement with fully-implicit solution method. To illustrate the applicability, the new solution scheme has been implemented into the best-estimate thermal-hydraulic analysis code, MARS. This paper also includes the comparisons of the simulation results for the perturbation propagation and water faucet problems using both two-step method and the original solution scheme.

• Bulletin of the Korean Chemical Society
• /
• v.15 no.6
• /
• pp.488-496
• /
• 1994

An approximate nonlinear theory of the Oregonator model is obtained with the aid of an ordinary perturbation method when the system is perturbed by some kinds of external input. The effects of internal and external parameters on the oscillations are discussed in detail by taking specific values of the parameters. A simple approximate solution for the Oregonator model under the influence of a constant input is obtained and the result is compared with the numerical result. For other types of external inputs the approximate solutions up to the fourth order expansion are compared with the numerical results. For a periodic input, we found that the entrainment depends crucially on the difference between the internal and external frequencies near the bifurcation point.