• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.10
• /
• pp.10-16
• /
• 2005

In this paper, a three-dimensional Euler aeroelastic analysis program is developed with a second-order staggered algorithm to reduce the lagging errors between the fluid and structural solvers. In the unsteady aerodynamic analysis, a dual-time stepping method based on the diagonalized-ADI algorithm is adopted to improve the time accuracy and a parallelized multi-grid method is used to save the computing time. The aeroelastic analyses of AGARD 445.6 wing model have been performed to verify the Euler aeroelastic analysis code. The analysis results are compared with the experimental data and other computational results. The results show comparatively good correlation when they are compared with other references.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2005.06a
• /
• pp.231-238
• /
• 2005

On the basis of flow resistance theory the conceptual model and related mathematical descriptions is proposed for resistance modeling of groundwater flow in CPM(continuum Porous medium), DFN(discrete fracture network) and fractured-porous medium. The proposed model is developed on the basis of finite volume method assuming steady-state, constant density groundwater flow. The basic approach of the method is to evaluate inter-block flow resistance values for a staggered grid arrangement, i.e. fluxes are stored at cell walls and scalars at cell centers. The balance of forces, i.e. the Darcy law, is utilized for each control volume centered around the point where the velocity component is stored. The transmissivity (or permeability) at the interface is assumed to be the harmonic average of neighboring blocks. Flow resistance theory was utilized to relate the fluxes between the grid blocks with residual pressures. The flow within porous medium is described by three dimensional equations and that within an individual fracture is described by a two dimensional equivalent of the flow equations for a porous medium. Newly proposed models would contribute to develop flow simulation techniques with various matrix characteristics.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.99-104
• /
• 2009

We propose a waveform inversion method for SH-wave data obtained in a shallow seismic refraction survey, to determine a 2D inhomogeneous S-wave profile of shallow soils. In this method, a 2.5D equation is used to simulate SH-wave propagation in 2D media. The equation is solved with the staggered grid finite-difference approximation to the 4th-order in space and 2nd-order in time, to compute a synthetic wave. The misfit, defined using differences between calculated and observed waveforms, is minimised with a hybrid heuristic search method. We parameterise a 2D subsurface structural model with blocks with different depth boundaries, and S-wave velocities in each block. Numerical experiments were conducted using synthetic SH-wave data with white noise for a model having a blind layer and irregular interfaces. We could reconstruct a structure including a blind layer with reasonable computation time from surface seismic refraction data.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.2
• /
• pp.159-164
• /
• 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 second 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 in Yangsan. The numerical result for urban district was visualization for 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.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.5 s.149
• /
• pp.538-549
• /
• 2006

A code is developed to simulate a viscous flow field around a deformable body using the hybrid Cartesian/immersed boundary method. In this method, the immersed boundary(IB) nodes are defined near the body boundary then velocities at the IB nodes are reconstructed based on the interpolation along the normal direction to the body surface. A new method is suggested to define the IB nodes so that a closed fluid domain is guaranteed by a set of IB nodes and the method is applicable to a zero-thickness body such as a sail. To validate the developed code, the vorticity fields are compared with other recent calculations where a cylinder orbits and moves into its own wake. It is shown the code can handle a sharp trailing edge at Reynolds number of $10^5$ under moderate requirements on girds. Finally the developed code is applied to simulate the vortex shedding behind a deforming foil with flapping tail like a fish. It is shown that the acceleration of fluids near the flapping tail contributes to the generation of the thrust for propulsion.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.181-186
• /
• 2003

An experimental study has been conducted to measure the temperature fields for two and three staggered rows of the rectangular shaped-holes with high turbulence intensity. 10 % turbulence intensity is obtained by installation of two kinds of grids which have different shapes. One grid which is installed at 30d upstream from center of 1st row of holes is composed of vertical cylinders of which diameter is 10 mm and center to center distance is 18 mm. The other installed 15d apart to upstream from center of 1st row of holes which has square pattern is constructed of 3 rum square bars and bar spacing is 25 mm. Temperature fields are measured by using a thermocouple rake which is attached on three-axis traversing system. The results show that the overall values are decreased and the thicker film of coolant is fanned downstream of rows of holes for high mainstream turbulence intensity.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.379-388
• /
• 2008

A hybrid Cartesian/immersed boundary code is expanded to simulate flow field around a three-dimensional body which undergoes large dynamic deformation. Immersed boundary nodes are automatically distributed based on the edges crossing triangles on body boundary. Velocity vectors are reconstructed at those immersed boundary nodes along local normal lines to the boundary. The reconstruction of pressure is avoided using the hybrid staggered/non-staggered grid method. The developed code is validated through comparisons with other results on laminar flow over a sphere. The code is applied to simulate flow around a foil which is attached to a body of revolution and rotates under periodic deformation. The periodic variation of the tip vortex is observed and the effects of the deformation on hydrodynamic force acting on the body are investigated.

• Journal of computational fluids engineering
• /
• v.7 no.2
• /
• pp.8-16
• /
• 2002

A numerical study on the use of the momentum interpolation method for flows with a large body force is presented. The inherent problems of the momentum interpolation method are discussed first. The origins of problems of the momentum interpolation methods are the validity of linear assumptions employed for the evaluation of the cell-face velocities, the enforcement of mass conservation for the cell-centered velocities and the specification of pressure and pressure correction at the boundary. Numerical experiments are performed for a typical flow involving a large body force. The numerical results are compared with those by the staggered grid method. The fact that the momentum interpolation method may result in physically unrealistic solutions is demonstrated. Numerical experiments changing the numerical grid have shown that a simple way of removing the physically unrealistic solution is a proper grid refinement where there is a large pressure gradient. An effective way of specifying the pressure and pressure correction at the boundary by a local mass conservation near the boundary is proposed, and it is shown that this method can effectively remove the inherent problem of the specification of pressure and pressure correction at the boundary when one uses the momentum interpolation method.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.395-406
• /
• 2011

In this research a numerical simulation method is developed for moving body in free surface flows using fixed staggered rectangular grid system. The non-linear free surface near the body is defined by marker-density method. The body boundary is defined by line segment connecting the points where the body surface and grid line meet. Continuity equation and Navier-Stokes equations are used as governing equations and the equations are coupled with two-step projection method. The velocities and pressures of body boundary and free surface cells are calculated with simultaneous iterative method. To treat a body movement in a fixed grid system, the volume displaced by moving body is added to the divergence of the body boundary cell. For the verification of the present numerical method. vortex shedding period of advancing cylinder is calculated and the period is compared with existing experiment results. Moreover, added mass and damping coefficients of a vertically excited box are calculated and the computed results are compared with published experiment results. Impulsive pressure and water level variation due to sloshing phenomenon are simulated and the results are compared with published experiment results. Varying the plunger shape, the waves generated by plunging type wave maker are compared with the 2nd order Stokes wave theory The plunger shape generating the wave that shows the best agreement with the theory is represented.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.