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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.1 no.3
• /
• pp.141-145
• /
• 1977

In this paper, approximate solutions of the von Karman equations for the free flexural vibration of a transversely isotropic thin rectangular plate with two simply supported edges and two clamped edges are obtained. Applying one term Ritz-Galerkin procedure, the spatial dependent part of the equation is separated and time dependent function is found to be the Duffing's equation. Then the relation between nonlinear period and amplitude of the vibration is obtained by using averaging method which is a method of the perturbation procedure. It can be seen that averaging method is easy and agrees well with prior results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2167-2173
• /
• 1996

An approximate weight function technique using the indirect boundary integral equation has been presented for the analysis of stress intensity foactors(SIFs) of a thick pipe. One-term boundary integral was introduced to represent the crack surface displacement field for the displacement based weight function technique. An explicit closed-form SIF solution applicable to symmetric cracked pipes without any modification of the solution including both circumferential and radial cracks has been derived. The necessary information in the analysis is two or three reference SIFs. In most cases the SIF solution were in good agreement with those available in the literature.

• International Journal of Safety
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2003

The present study deals with an approximate integral equation approach to finite deflection of elastic plates with arbitrary plane form. An integral formulation leads to a system of boundary integral equations involving values of deflection, slope, bending moment and transverse shear force along the edge. The basic principles of the development of boundary element technique are reviewed. A computer program for solving for stresses and deflections in a isotropic, homogeneous, linear and elastic bending plate is developed. The fundamental solution of deflection and moment is employed in this program. The deflections and moments are assumed constant within the quadrilateral element. Numerical solutions for sample problems, obtained by the direct boundary element method, are presented and results are compared with known solutions.

• Journal of information and communication convergence engineering
• /
• v.12 no.1
• /
• pp.39-45
• /
• 2014

A numerical method for solving Hamiltonian equations is said to be symplectic if it preserves the symplectic structure associated with the equations. Various symplectic methods are widely used in many fields of science and technology. A symplectic method preserves an approximate Hamiltonian perturbed from the original Hamiltonian. It theoretically supports the effectiveness of symplectic methods for long-term integration. Although it is also related to long-term integration, numerical stability of symplectic methods have received little attention. In this paper, we consider explicit symplectic methods defined for Hamiltonian equations with Hamiltonians of the special form, and study their numerical stability using the harmonic oscillator as a test equation. We propose a new stability criterion and clarify the stability of some existing methods that are visually based on the criterion. We also derive a new method that is better than the existing methods with respect to a Courant-Friedrichs-Lewy condition for hyperbolic equations; this new method is tested through a numerical experiment with a nonlinear wave equation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.252-259
• /
• 1999

One of the most important factors for a proper design of a slender compression member may be the exact determination of the elastic critical load of that member. In the cases of non-prismatic compression member, however, there are times when the exact critical load becomes impossible to determinate if one relies on the neutral equilibrium method or energy principle. Here in this paper, the approximate critical loads of symmetrically or non-symmetrically tapered members are computed by finite element method. The two parameters considered in this numerical analysis are the taper parameter, $\alpha$ and the sectional property parameters, m. The computed results for each sectional property parameter, m are presented in an algebraic equation which agrees with those by F.E.M The algebraic equation can be easily used by structural engineers, who are engaged in structural analysis and design of non-prismatic compression member.

• 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.24 no.5 s.176
• /
• pp.1166-1174
• /
• 2000

In this paper, we introduce a positioner based on magnetic levitation to eliminate the friction which is the most severe effect to limit high resolution on the micro level. Differently from existing electromagnetic device, the proposed positioner consists of air core solenoid and permanent magnet. Although the combination produces small magnetic force, it is suitable for realizing micro motion repeatedly without the accumulation of error because there is no hysteresis caused by ferromagnetic materials, no eddy current loss, no flux saturation. First, the approximate modeling of stiffness and damping effects between the magnetic elements is made and verified experimentally. Then, we have formulated the dynamic equation of one d.o.f magnetic levitation positioner using linear perturbation method and discussed the necessity of optimization for the chief design parameters to maximize the stability performance.

• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.227-233
• /
• 2004

Simplified model of slide-fed curtain coating flow has been developed and tested in this study. It rests on the sheet profile equations for curtain thickness in curtain flow and its trajectory derived by the integral momentum balance approach of Higgins and Scriven (1979) and Kistler (1983). It also draws on the film profile equation of film thickness variation in flow down a slide. The equations have been solved in finite difference approximation by Newton iteration with continuation. The results show that how inertia (Rey­nolds number), surface tension (capillary number), inclination angle of the slide, and air pressure difference across the curtain affect sheet trajectory and thickness profile. It has been revealed that approximate models can be useful to easily analyze coating flow dynamics without complex computations, giving qualitative agreement with full theory and with experiment.

• Journal of the Korean Mathematical Society
• /
• v.51 no.2
• /
• pp.251-266
• /
• 2014

In this paper, we use Newton's method to approximate a locally unique solution of an equation in Banach spaces and introduce recurrent functions to provide a weaker semilocal convergence analysis for Newton's method than before [1]-[13], in some interesting cases, provided that the Fr$\acute{e}$chet-derivative of the operator involved is p-H$\ddot{o}$lder continuous (p${\in}$(0, 1]). Numerical examples involving two boundary value problems are also provided.