• 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.9 no.2
• /
• pp.139-147
• /
• 2006

We had conducted a numerical modeling to investigate seismic properties of gas hydrate with field parameters acquired over the East sea in 1998. We used a 2-D staggered grid finite difference method to generate synthetic elastic seismograms for multi-channel seismic survey, OBC (Ocean Bottom Cable) survey and VCS (Vertical Cable Seismic) survey. The results of this study showed that the method using staggered grid yielded stable results and could be used to seismic imaging. We could find out the high amplitude anomaly and the phase reversal phenomenon of reflection wave at interface between the gas hydrate layer and free gas layer such a BSR (Bottom Simulating Reflector) which is the evidence for existence of gas hydrate in seismic reflection data. And we computed the reflection coefficients at the incident angles corresponding to offset distance with the synthetic seismograms. The reflection coefficients acquired from the numerical modeling were nearly consistent with the reflection coefficient computed by Shuey's equation.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.3
• /
• pp.66-81
• /
• 2021

In this paper, we review the lowest order staggered discontinuous Galerkin methods on polygonal meshes in 2D. The proposed method offers many desirable features including easy implementation, geometrical flexibility, robustness with respect to mesh distortion and low degrees of freedom. Discrete function spaces for locally H¹ and H(div) spaces are considered. We introduce special properties of a sub-mesh from a given star-shaped polygonal mesh which can be utilized in the construction of discrete spaces and implementation of the staggered discontinuous Galerkin method. For demonstration purposes, we consider the lowest case for the Poisson equation. We emphasize its efficient computational implementation using only geometrical properties of the underlying mesh.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.179-182
• /
• 2008

This study simulated the flood inundations of the Nakdong River catchment running through Yangsan, a small city located in the south eastern area of Korea by using the depth averaged two-dimensional hydrodynamic numerical model. The numerical model employs the staggered grid system including moving boundary and a finite different method to solve the Saint-Venant equations. A second order upwind scheme is used to discretize the nonlinear convection terms of the momentum equations, whereas linear terms are discretized by a first order leap-frog scheme(Cho and Yoon, 1998). The numerical model was applied to a real topography to simulate the flood inundation of the Yangsan basin. The numerical results for urban district are visualized in three dimension. These results can be essentially utilized to construct the three dimensional inundation map after building the GIS-based database in local public organizations in order to protect the life and property safely.

• 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.29 no.3
• /
• pp.59-70
• /
• 1992

A Finite volume method for the computation of the two-dimensional, incompressible, steady, laminar Navier-Stokes equation is developed using a non-staggered grid system in a general curvilinear coordinate. The numerical pressure fluctuations, usually encountered when the non-staggered grid system is used, is suppressed by the momentum interpolation method. Flows around a NACA0012 foil section have been computed by the present method and the results show good agreements with other experimental and numerical ones.

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

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

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