• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.57-62
• /
• 2008

Finite difference method using not general SSG(standard staggered grid) but RSG(rotated staggered grid) was applied to simulation of elastic wave propagation. Special free surface boundary condition such as imaging method is needed in finite difference method using SSG in elastic wave propagation but free surface boundary condition in finite difference method using RSG is easily solved with adding air layer. Recently PML(Perfectly Matched layer) is widely used to eliminate artificial reflection waves from finite boundary because of its' greate efficiency. Absorbing ability of CPML(convolutional Perfectly Matched Layer) that is more efficient than that of PML was applied to FDM using RSG in this study. The results of CPML eliminated artificial boundary waves very effectively in FDM using RSG in being compared with that of Cerjan's absorbing method.

• Geophysics and Geophysical Exploration
• /
• v.12 no.2
• /
• pp.183-191
• /
• 2009

Finite difference method using not general SSG (standard staggered grid) but RSG (rotated staggered grid) was applied to simulation of elastic wave propagation. Special free surface boundary condition such as imaging method is needed in finite difference method using SSG in elastic wave propagation. But free surface boundary condition in finite difference method using RSG is easily solved with adding air layer or vacuum layer. Recently PML (Perfectly Matched layer) is widely used to eliminate artificial reflection waves from finite boundary because of its' greate efficiency. Absorbing ability of CPML (convolutional Perfectly Matched Layer) that is more efficient than that of PML and CPML that don't use splitting of wave equation that should be adapted to PML was applied to FDM using RSG in this study. Frequency absorbing characteristic and energy absorbing ability in CPML layer were investigated and CPML eliminated artificial boundary waves very effectively in FDM using RSG in being compared with that of Cerjan's absorbing method. CPML method also diminished amplitude of waves in boundary layer of solid-liquid model very well.

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.1
• /
• pp.212-218
• /
• 1995

Turbulent flow and heat transfer characteristics around staggered tube bundles were studied using a non-orthogonal boundary fitted coordinate system and the low Reynolds .kappa. - .epsilon. turbulence model suggested by Lam and Bremhorst. The predicted flow characteristics for two tube pitches and tube arrangement showed good agreement with the experimental data except the strongly curved region. The predicted Nusselt number was compared with measurements obtained in the staggered rough bundles and it revealed the similar trend to measurements, but the location of the maximum and minimum heat transfer differed somewhat from the measurements.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.20-28
• /
• 2016

Recently many sinkholes have appeared in urban areas of Korea, threatening public safety. To predict the occurrence of sinkholes, it is necessary to investigate the existence of cavity under urban roads. Ground-penetrating radar (GPR) has been recognized as an effective means for detecting underground cavity in urban areas. In order to improve the understanding of the governing physical processes associated with GPR wave propagation, and interpret underground cavity effectively, a theoretical approach using numerical modeling is required. We have developed an algorithm employing a three-dimensional (3D) staggered-grid finite-difference time-domain (FDTD) method. This approach allows us to model the full electromagnetic wavefield associated with GPR surveys. We examined the GPR response for a simple cavity model, and the modeling results showed that our 3D FDTD modeling algorithm is useful to assess the underground cavity under urban roads.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.8
• /
• pp.2611-2625
• /
• 1996

A numerical method on multi-block technique by Bi-CGSTAB(Bi-Conjugate Gradient STABilized) solver has been proposed. The present multi-block technique can reduce the numerical manipulation greatly because the common regions at the interface of each block are not necessary. In order to test the computational performance of present multi-block technique, the flow characteristics in a T type duct system and a N type duct system have been investigated by three kinds of methods such as the single-block method, the previous multi-block technique and the multi-block technique with Bi-CGSTAB solver. The results indicated that the required CPU time by present multi block technique was shorter than that of other two numerical methods and the convergency history was shown very stable at the present multi-block technique.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2001

Wave-induced currents drive nearshore transport processes, and hence an accurate understanding of wave-current interaction is required for proper management of coastal zone. This paper presents details of an adaptive quadtree grid based numerical model of the coupled wave climate and depth-averaged current field. The model accounts for wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, mixing processes, bottom friction effects, and movement of land-water interface at the shoreline. The wave period- and depth-averaged governing equations arc discrctized explicitly by means of an Adarns¬Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Results from the numerical model are in reasonable agreement with the laboratory data of longshore current generated by oblique waves on a plane beach (Visser 1980, 1991).

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.127-132
• /
• 1998

The present study is aimed to investigate flow characteristics of Two dimensional backward-facing step by numerical approach. A convection conservative difference scheme based upon SOLA algorithm is used for the solution of the two-dimensional incompressible Navier-Stokes equations to simulate the laminar, transitional and turbulent flow conditions at which the experimental data can be available for the backward-facing step. The twenty kinds of Reynolds number are used for the calculations. In an effort to demonstrate that the reported solutions are dependent on the mesh refinement, computations are performed on seven different meshes of uniformly increasing refinement. Also to investigate the result of inflow dependence, two kinds of the inflow profile are chosen for the laminar flow. As criterion of benchmarking the result of numerical simulation, reattachment length is used for the selected Reynolds numbers.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2006.03a
• /
• pp.25-32
• /
• 2006

We present a imaging method of seismic sources by time reversal propagation of seismic waves. Time-reversal wave propagation is actively used in medical imaging, non destructive testing and waveform tomography. Time-reversal wave propagation is based on the time-reversal invariance and the spatial reciprocity of the wave equation. A signal is recorded by an array of receivers, time-reversed and then back-propagated into the medium. The time-reversed signal propagates back into the same medium and the energy refocuses back at the source location. The increasing power of computers and numerical methods makes it possible to simulate more accurately the propagation of seismic waves in heterogenous media. In this work, a staggered-grid finite-difference solution of the elastic wave equation is employed for the wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.2
• /
• pp.70-76
• /
• 2000

A numerical study on the viscous flow around a 2-dimensional circulation control foil is carried out for application on the field of naval architecture and ocean engineering. The governing equations are the RANS and the continuity equations. The equations are discretized by finite difference method and MAC method and the pressure poisson equation is calculate by a SOR method and an O-type non-staggered boundary fitted coordinate system which is overlapped near the slot is used to improve the numerical accuracy. Turbulence is approximated by a modified Baldwin-Lomax turbulence model. In the present paper, the Coanda effect on a 2-dimensional foil of a 20% thickness ellipse with modified rounded trailing edge has been numerically studied. The change in drag and lift of the foil with various jet momentums are calculated and compared to the experimental results to show good agreements.