• 한국지구물리탐사학회:학술대회논문집
• /
• 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.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.218-226
• /
• 2008

A component-scale two-phase analysis code has been developed for a realistic simulation of two-phase flow transients in a light water nuclear reactor component. In the code, a two-fluid three-field model is adopted and the governing equations are solved on an unstructured mesh. For the numerical solution scheme, the semi-implicit method used in the RELAP5 code was selected, which has been proved to be very stable and accurate for most of practical applications. However, some modifications were needed for its application to an unstructured non-staggered grid. This paper presents the modified semi-implicit numerical method for unstructured grid and the preliminary results of the calculations.

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

• The Journal of the Acoustical Society of Korea
• /
• v.23 no.8
• /
• pp.576-582
• /
• 2004

A computer code that is based on the Pseudo-spectral finite difference algorithm using staggered grid is developed for the wave propagation modeling in the time domain. The advantage of a finite difference approximation is that any geometrically complicated media can be modeled. Staggered grids are advantageous as it provides much more accuracy than using a regular grid. Pseudo-spectral methods are those that evaluate spatial derivatives by multiplying a wavenumber by the Fourier transform of a pressure wave-field and performing the inverse Fourier transform. This method is very stable and reduces memory and the number of computations. The synthetic results by this algorithm agree with the analytic solution in the infinite and half space. The time domain modeling was implemented in various models. such as half-space. Pekeris waveguide, and range dependent environment. The snapshots showing the total wave-field reveals the Propagation characteristic or the acoustic waves through the complex ocean environment.

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.1-8
• /
• 2007

A code is developed using the hybrid Cartesian/immersed boundary method and it is applied to simulate flows around a three-dimensional deforming body. A new criterion is suggested to distribute the immersed boundary nodes based on edges crossing a body boundary. Velocities are reconstructed at the immersed boundary nodes using the interpolation along a local normal line to the boundary. Reconstruction of the pressure at the immersed boundary node is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other experimental and numerical results for the velocity profiles around a circular cylinder under the forced in-line oscillation and the pressure coefficient distribution on a sphere. The code is applied to simulate the flow fields around a plate whose tail is periodically flapping under a translation. The effects of the velocity and acceleration due to the deformation on the periodic shedding of pairs of tip vortices are investigated.

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

Various discretiation methods of Laplacian operator in the Pressure-Poisson equation are investigated for the solution of incompressible Navier-Stokes equations using the non-staggered grid. Laplacian operators previously proposed by other researchers are applied to a Driven-Cavity problem. The computational results are compared with those of Ghia. The results show the characteristics of the discrete Laplacian operators.

• Journal of computational fluids engineering
• /
• v.7 no.1
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.640-653
• /
• 2011

We simulate the propagation of earthquake waves in the continental margin of Antarctica using the elastic wave modeling algorithm, which is modified to be suitable for acoustic-elastic coupled media and earthquake source. To simulate the various types of earthquake source, the staggered-grid finite-difference method, which is composed of velocity-stress formulae, can be more appropriate to use than the conventional, displacement-based, finite-difference method. We simulate the elastic wave propagation generated by earthquakes combining 3D staggered-grid finite-difference algorithm composed of displacement-velocity-stress formulae with double couple mechanisms for earthquake source. Through numerical tests for left-lateral strike-slip fault, normal fault and reverse fault, we could confirm that the first arrival of P waves at the surface is in a good agreement with the theoretically-predicted results based on the focal mechanism of an earthquake. Numerical results for a model made after the subduction zone in the continental margin of Antarctica showed that earthquake waves, generated by the reverse fault and propagating through the continental crust, the oceanic crust and the ocean, are accurately described.