• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제5권2호
• /
• pp.25-37
• /
• 2001

We consider the finite element method applied to nonlinear Sobolev equation with smooth data and demonstrate for arbitrary order ($k{\geq}2$) finite element spaces the optimal rate of convergence in $L_{\infty}\;W^{1,{\infty}}({\Omega})$ and $L_{\infty}(L_{\infty}({\Omega}))$ (quasi-optimal for k = 1). In other words, the nonlinear Sobolev equation can be approximated equally well as its linear counterpart. Furthermore, we also obtain superconvergence results in $L_{\infty}(W^{1,{\infty}}({\Omega}))$ for the difference between the approximate solution and the generalized elliptic projection of the exact solution.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2000년도 추계 학술발표회 논문집 Proceedings of EESK Conference-Fall 2000
• /
• pp.81-88
• /
• 2000

The elastic wave equation is solved using the finite-difference method in 3D space to simulate the seismic wave propagation. It is based on the velocity-stress formulation of the equation of motion on a staggered grid. The nonreflecting boundary conditions are used to attenuate the wave field close to the numerical boundary. To satisfy the stress-free conditions at the free-surface boundary, a new formulation combining the zero-stress formalism with the vacuum one is applied. The effective media parameters are employed to satisfy the traction continuity condition across the media interface. With use of the moment-tensor components, the wide range of source mechanism parameters can be specified. The numerical experiments are carried out in order to test the applicability and accuracy of this scheme and to understand the fundamental features of the wave propagation under the generalized elastic media structure. Computational results show that the scheme is sufficiently accurate for modeling wave propagation in 3D elastic media and generates all the possible phases appropriately in under the given heterogeneous velocity structure. Also the characteristics of the ground motion in an sedimentary basin such as the amplification, trapping, and focusing of the elastic wave energy are well represented. These results demonstrate the use of this simulation method will be helpful for modeling the ground motion of seismological and engineering purpose like earthquake hazard assessment, seismic design, city planning, and etc..

• Journal of applied mathematics & informatics
• /
• 제25권1_2호
• /
• pp.435-443
• /
• 2007

Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

• Tribology and Lubricants
• /
• 제24권1호
• /
• pp.14-20
• /
• 2008

In this paper, the characteristics of trapezoidal air screws is studied. These screws can e applied to the elevator having column driving shaft. The generalized coordinate transformation technique is used to solve incompressible Reynolds' equation because the air lubricated plane is twisted. The transformed equation is discretized by the base of Finite difference method. Using Visual C++ language, a GUI program which can calculate he load carrying capacity for this kind of air screw is developed, then the design variables for these air screws is studied.

• 대한기계학회논문집B
• /
• 제22권12호
• /
• pp.1796-1804
• /
• 1998

A numerical investigation has been performed to discuss the radiation-affected steady-laminar natural convection induced by a hot inner cylinder under a large temperature difference in the annuli filled with a gray gas. To examine the effects of thermal radiation on thermo-fluid dynamic behaviors in the eccentric geometry, the generalized body-fitted coordinate system is introduced while the finite volume method (FVM) is used for solving the radiative transport equation. After validating the numerical results for the case without radiation, the detailed radiation effect has been discussed. Based on the results of this study, when there exists a large temperature difference between two cylinders, the existence of radiatively participating medium is found to incur a distinct difference in fluid dynamic as well as thermal behavior.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2009년도 춘계 학술발표회
• /
• pp.434-443
• /
• 2009

Slope stability analysis is an essential part of rock slope design. For highly fractured rock, the limit equilibrium method (LEM) based slope stability analysis with a circular failure surface is often carried out assuming the rock mass behaves more or less as a continuum. This paper examines first, the applicability of the finite-element method (FEM) based shear strength reduction (SSR) technique for highly fractured rock slope, and second the use of Mohr-Coulomb (MC) failure criterion in conjunction with generalized Hoek-Brown (HB) failure criterion. The numerical results on a number of cases are compared in terms of the factor of safety (FS). The results indicated that the FEM-based SSR technique yields almost the same FSs from LEM, and that the MC and HB failure criteria yield almost identical FSs when the strength parameters for MC failure criterion are obtained based on the modified HB failure criterion if and only if value of the Hoek-Brown constant $m_i$ is smaller than 10 and slope angle is smaller than 1:1, otherwise MC failure criteria over-estimate the factor of safety.

• 한국전산유체공학회지
• /
• 제10권1호
• /
• pp.31-38
• /
• 2005

Three-dimensional flow over the sand dunes have been studied numerically by using Large-Eddy Simulation (LES) method. In the direction of initial flow and span direction cyclic boundary conditions are imposed for velocity and pressure. The movement of the sand dune which is formed by converging wind direction has been investigated. The numerical method employed in this study can be divided into three parts: (i) calculation of the air flow over the sand dune using standard MAC method with a generalized coordinate system; (ii) estimation of the sand transfer caused by the flow through the friction; (iii) determination of the shape of the sand surface. Since the computational area has been changed due to step (iii), (i)-(iii) are repeated. The simulated dune, which has initially elliptic cross section, extends at the converging direction, which is known as linear dunes.

• 대한의용생체공학회:학술대회논문집
• /
• 대한의용생체공학회 1997년도 추계학술대회
• /
• pp.272-275
• /
• 1997

Noninvasive multifrequency microwave radiometry using coaxial waveguide antenna has been investigated for a homogeneous and our layer human body model. We derived finite-difference time-domain(FDTD) algorithm and equation of MUR and generalized perfectly matched layer(GPML) absorbing boundary conditions(ABCs) in cylindrical coordination. The coupling between coaxial waveguide antenna and a biological object was analyzed by use of the FDTD method using MUR and GPML ABCs to obtain the absorbed power patterns in the media. The specific absorption rates(SAR) distribution which was corresponding to the temperature distribution was calculated in each region by use of the steady-state response in FDTD method. The SAR patterns of FDTD method using MUR ABCs was compared with those of FDTD method using GPML ABCs.

• Korean Journal of Hydrosciences
• /
• 제5권
• /
• pp.85-97
• /
• 1994

A hybrid finite difference method for the longitudinal dispersion equation, which is based on combining the Holly-Preissmann scheme with fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme, is described and comparatively evaluated with other characteristics-based numerical methods. Longitudinal dispersion of an instantaneously-loaded pollutant source is simulated, and computational results are compared with the exact solution. The present method is free from wiggles regardless of the Courant number, and exactly reproduces the location of the peak concentration. Overall accuracy of the computation increases for smaller value of the weighting factor, $\theta$of the model. Larger values of $\theta$ overestimates the peak concentration. Smaller Courant number yields better accuracy, in general, but the sensitivity is very low, especially when the value of $\theta$ is small. From comparisons with the hybrid method using cubic interpolating polynomial and with splitoperator methods, the present method shows the best performance in reproducing the exact solution as the advection becomes more dominant.

• Geomechanics and Engineering
• /
• 제16권1호
• /
• pp.59-69
• /
• 2018

This study focused on the mechanical and hydraulic characteristics of underwater tunnels based on Mohr-Coulomb (M-C), Hoek-Brown (H-B) and generalized H-B failure criteria. An improved approach for calculating stress, displacement and plastic radius of the circular tunnel considering hydraulic-mechanical coupling was developed. The innovation of this study was that the radius-incremental-approach was reconstructed (i.e., the whole plastic zone is divided into a finite number of concentric annuli by radius), stress and displacement of each annulus were determined in terms of numerical method and Terzaghi's effective stress principle. The validation of the proposed approach was conducted by comparing with the results in Brown and Bray (1982) and Park and Kim (2006). In addition, the Rp-pin curve (plastic radius-internal supporting pressure curve) was obtained using the numerical iterative method, and the plastic radius of the deep-buried tunnel could be obtained by interpolation method in terms of the known value of internal supporting pressure pin. Combining with the theories in Carranza and Fairhurst (2000), the improved technique for assessing the reliability of the tunnel support was proposed.