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

• 유체기계공업학회:학술대회논문집
• /
• 2002.12a
• /
• pp.453-458
• /
• 2002

The feasibility study on supporting a turbocharger rotor on air foil bearing is investigated. Based on finite difference method and Newton-Raphson method, the static equilibrium position of a turbocharger rotor is predicted. And using finite difference method and perturbation method, dynamic characteristics of air foil bearings are calculated. Rotordynamic analysis is performed by finite element method, with collaboration of calculated stiffness and damping of foil bearing. The effect of compliance and clearance of bump foil bearing on the oil-free turbocharger is investigated in terms of rotordynamics. And the critical speeds, eccentricity ratio, vibration amplitude, and stability are considered. It is demonstrated that foil bearings offer a rlausible replacement for oil-lubricated bearings in turbocharger.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.3
• /
• pp.175-187
• /
• 2010

This paper presents the numerical valuation of the two-asset step-down equitylinked securities (ELS) option by using the operator-splitting method (OSM). The ELS is one of the most popular financial options. The value of ELS option can be modeled by a modified Black-Scholes partial differential equation. However, regardless of whether there is a closedform solution, it is difficult and not efficient to evaluate the solution because such a solution would be represented by multiple integrations. Thus, a fast and accurate numerical algorithm is needed to value the price of the ELS option. This paper uses a finite difference method to discretize the governing equation and applies the OSM to solve the resulting discrete equations. The OSM is very robust and accurate in evaluating finite difference discretizations. We provide a detailed numerical algorithm and computational results showing the performance of the method for two underlying asset option pricing problems such as cash-or-nothing and stepdown ELS. Final option value of two-asset step-down ELS is obtained by a weighted average value using probability which is estimated by performing a MC simulation.

• Tunnel and Underground Space
• /
• v.12 no.2
• /
• pp.83-91
• /
• 2002

Stability analysis of rocky embankment slopes is done by both the limit equilibrium method and the finite difference method. The height or the rocky embankment is approximately 40 m and the side slope is 1 vertical to 1.5 horizontal. The cohesion and internal friction angle of rock debris are assumed zero and 43$^{\circ}$, respectively. For finite difference analysis, strength reduction method is used to calculate the saft factor of the slope. As a result, the safety factor of the slope is discovered to be 1.4 by using either methods. Considering that the design criteria of the safety factor is 1.3, it can be judged that the rock fragments embankment slope is in a stable state.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.195-199
• /
• 2005

Wave breaking phenomenon near the fore body of a ship is numerically simulated. The ship advance with uniform velocity in calm water. For the simulation, incompressible Navier-Stokes equations and continuity equation are adopted as governing equations. The simulation is carried out in staggered variable mesh system with finite difference method. Marker and Cell(MAC) method and Marker-Density method are employed to track the free surface. Body boundary conditions are satisfied with the adoption of porosity method and no-slip condition on the hull surface. The ship model has a wedge type fore-body, and the computational domain is an appropriate region around the fore-body. The computation results are compared with some experimental results. Also the difference of the free surface tracking methods are discussed.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.4
• /
• pp.5-10
• /
• 2006

A special finite difference method for nonlinear dynamic response analysis of semi-infinite foundation soil using mapping which transforms semi-infinite domain into finite domain is presented here. For the region of engineering interest, mapping is isometric, and fur far field, shrink mapping which transforms infinite interval into finite interval is adopted. At first, the responses of semi-infinite foundation soil with linear constituting model are computed, and compared with theoretical results and those of existing method. Good agreements are obtained among the results of the proposed method, Lamb's theory and FEM with extensive mesh model. Then the responses of infinite foundation soil are computed by the present method, using small and large mesh model. The results of small and large mesh models agree well with each other, demonstrating the effectiveness of the proposed method.

• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.9
• /
• pp.110-121
• /
• 2001

During hot pipe bending using induction heating, the wall of bending outside is thinned by tensile stress. In design requirement, the reduction of wall thickness is not allowed to exceed 12.5%. So in this study, two methods of bending, one is loading of reverse moment and the other is loading of temperature gradient, have been investigated to design pipe bending process that satisfy design requirements. For this purpose, finite element analysis with a bending radius 2Do(outer diameter of pipe) has been performed to calculate proper reverse moment and temperature gradient to be applied. Induction heating process has been analyzed to estimate influence of heating process parameters on heating characteristic by finite difference method. Then pipe bending experiments have been performed for verification of finite element and finite difference analysis results. Experimental results are in good agreement with the results of simulations.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.121-128
• /
• 2017

Interpretation of time-domain electromagnetic (TEM) data has been made mostly based on one-dimensional (1-D) inversion scheme in Korea. A proper interpretation of TEM data should employ 3-D TEM forward and inverse modeling algorithms. This study developed a 3-D TEM modeling algorithm using a finite difference time-domain (FDTD) method with staggered grid. In numerically solving Maxwell equations, fictitious displacement current is included based on an explicit FDTD method using a central difference approximation scheme. The developed modeling algorithm simulated a small-coil source configuration to be verified against analytic solutions for homogeneous half-space models. Further, TEM responses for a 3-D anomaly are modeled and analyzed. We expect that it will contribute greatly to the precise interpretation of TEM data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5A
• /
• pp.457-465
• /
• 2009

This study presents a moving least squares (MLS) finite difference method for solving elastic crack problems with stress singularity at the crack tip. Near-tip functions are intrinsically employed in the MLS approximation to model near-tip field inducing singularity in stress field. employment of the functions does not lose the merit of the MLS Taylor polynomial approximation which approximates the derivatives of a function without actual differentiating process. In the formulation of crack problem, computational efficiency is considerably improved by taking the strong formulation instead of weak formulation involving time consuming numerical quadrature Difference equations are constructed on the nodes distributed in computational domain. Numerical experiments for crack problems show that the intrinsically enriched MLS finite difference method can sharply capture the singular behavior of near-tip stress and accurately evaluate stress intensity factors.