• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.313-325
• /
• 2011

Finite element grid system using h-refinement on the Yellow Sea was constructed based on previous study (Suh, 1999b) from 14 K to 210 K and special attention was concentrated on refining the coastal zone. In grid generation, depth change between adjacent points and non-dimensional tidal wave length ratio were considered. As a result approximately a quarter of the total nodes are located nearby 5 m of shallow area. Accurate bathymetry data using 30's and ETOPO1 with open boundary conditions of 8 major tidal constituents extracted automatically from FES2004 have been applied. In tidal simulation a 3-dimensional nonlinear harmonic model was setup and tidal amplification due to changes in vertical turbulent and bottom friction were simulated. In this study not only 8 major tidal constituents but also nonlinear shallow tides $M_4,$, $MS_4$ and long period $M_f,$, $M_{sf}$ were reproduced. It is found that implication of spatial variation of friction coefficient plays a very important role in reproduction of astronomical and shallow tides which are computed by iterative computation of nonlinear terms. Also it should be considered differently with respect to tidal periods. To understand the distribution of tidal asymmetry, amplitude ratio of $M_4/M_2$ and phase differences $2g(M_2)-g(M_4)$ were calculated. Tidal distortion ratio marks up to 0.2 on the west coast showing shallow coastal characteristics and somewhat wide range of ebb-dominances in front of Mokpo area are reproduced.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.1
• /
• pp.39-45
• /
• 2005

In the situation which Russia's ratification of the Kyoto protocol at February,2004, ANNEX I nations must reduce GHG(Green House Gas) discharge rate from 2008 by 2012 to the reduction level at 1990. We introduce the CO₂ ocean sequestration that is one of promising method for getting the stable CO₂ concentration in the atmosphere. There are four categories : ocean transportation technique, ocean initial dissolution technique, ocean deep current evaluation technique, and ocean biological evaluation technique. In this paper, we carried out the fundamental numerical study on the ocean initial dissolution technique, when the Liquidized CO₂ is emitted at the deep ocean, It is very important to the dissolution rate of movable CO₂ interface because it Is directly impact to the ocean organism. In order to investigate the relation of the interface movement and rate of the dissolution, we develope CR(Computational Fluid Dynamics) code that was constructed by the finite volume method based on the unstructured mesh, and a droplet's boundary surface can move and one direction dissolution from disperse phase into continuous phase adopted as its physics be. This study clarifies hydrodynamic relation between solubility and movement of the droplet through the verification of the Cm code.

• Journal of Korea Water Resources Association
• /
• v.56 no.spc1
• /
• pp.1059-1069
• /
• 2023

Dam-break flow occurs when an elevated dam suddenly collapses, resulting in the catastrophic release of rapid and uncontrolled impounded water. This study compares laminar and turbulent closure models for simulating three-dimensional dam-break flows using OpenFOAM. The Reynolds-Averaged Navier-Stokes (RANS) model, specifically the k-ε model, is employed to capture turbulent dissipation. Two scenarios are evaluated based on a laboratory experiment and a modified multi-layered block obstacle scenario. Both models effectively represent dam-break flows, with the turbulent closure model reducing oscillations. However, excessive dissipation in turbulent models can underestimate water surface profiles. Improving numerical schemes and grid resolution enhances flow recreation, particularly near structures and during turbulence. Model stability is more significantly influenced by numerical schemes and grid refinement than the use of turbulence closure. The k-ε model's reliance on time-averaging processes poses challenges in representing dam-break profiles with pronounced discontinuities and unsteadiness. While simulating turbulence models requires extensive computational efforts, the performance improvement compared to laminar models is marginal. To achieve better representation, more advanced turbulence models like Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) are recommended, necessitating small spatial and time scales. This research provides insights into the applicability of different modeling approaches for simulating dam-break flows, emphasizing the importance of accurate representation near structures and during turbulence.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.60-67
• /
• 2008

In this paper we consider an application of the method of electromagnetic (EM) migration to the interpretation of a typical marine controlled-source (MCSEM) survey consisting of a set of sea-bottom receivers and a moving electrical bipole transmitter. Three-dimensional interpretation of MCSEM data is a very challenging problem because of the enormous number of computations required in the case of the multi-transmitter and multi-receiver data acquisition systems used in these surveys. At the same time, we demonstrate that the MCSEM surveys with their dense system of transmitters and receivers are extremely well suited for application of the migration method. In order to speed up the computation of the migration field, we apply a fast form of integral equation (IE) solution based on the multigrid quasi-linear (MGQL) approximation which we have developed. The principles of migration imaging formulated in this paper are tested on a typical model of a sea-bottom petroleum reservoir.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.04a
• /
• pp.147-152
• /
• 1995

실험에 의한 모우드 해석 방법들은 1980년대부터 활발히 연구되어 많은 새로운 방법들이 개발되어 발표되었다. 그러나 개발된 대부분의 방법들은 측정된 데이타를 일괄처리하는 밸치(또는 off-line) 방법들이다. 최근에는 시간에 따라서 변하는 구조물의 동특성을 규명하는 분야에 모우드 해석 방법이 응용되어 사용되고 있다. 이러한 응용분야에서는 모우드 변수들의 변화되는 값을 새로운 데이타가 샘플링 될 때마다 그 값들을 수정하면서 추정할 수 있는 회귀적인(recursive 또는 on-line) 방법을 사용하여야 한다. Davies와 Hammond[1]는 회귀적 선형 자승법(Recursive Least Squares : RLS)을 이용하여 모우드 변수를 구하고 이를 벧치방법인 Instrumental Variable 방법과 Fourier 방법의 결과와 비교하였다. 그러나, 그 결과에서 보여준것처럼 RLS 방법은 잡음 대 시호비가 낮을 때에만 모우드 변수 값들을 정확하게 추정할 수 있었다. Sundararajan과 Montgomrey[2]는 회귀적 선형 최소자승 격자필터(lattice filter)를 이용하여 구조물의 차수(order)와 고유진동형, 그리고 진폭을 결정한 후 이를 토대로 회귀적 gradient형태의 방정식 오차 규명 방법(equation-error identification algorithm)에 의하여 모우드 변수들을 추정하였다. 이 방법은 2차원 격자구조물의 모우드 변수 추정에 사용되었으며, 또한 적응모우드제어에도 성공적으로 이용되었다. 그러나, 이 방법도 잡음 대 신호비가 낮은 환경에서만 사용할 수 있다는 단점이 있다. 위에서 언급한 방법들은 모두 RLS 방법을 기초로 하여 개발되었으나, RLS 방법은 전형적인 결정적(deterministic)방법으로서 잡음이 섞인 데이타를 처리하기에는 부적절한 방법임이 널리 알려진 사실이다[3]. 최근에 Ben Mrad와 Fassois[4]는 신호에 잡음이 존재하여도 이를 잘 처리할 수 있는 확률적(stochastic) 방법을 개발하여 기존의 결정적 방법들과 그 결과를 비교하였다. 그러나, 개발된 방법은 응답 신호에 백색잡음(white noise)이 섞이는 특수한 경우에만 사용할 수 있게 만들어져서 이 방법의 실질적인 적용에는 어려움이 있다. 본 연구에서는 기존의 방법들의 단점을 극복할 수 있는 새로운 회귀적 모우드 변수 규명 방법을 개발하였다. 이는 Fassois와 Lee가 ARMAX모델의 계수를 효율적으로 추정하기 위하여 개발한 뱉치방법인 Suboptimum Maximum Likelihood 방법[5]를 기초로 하여 개발하였다. 개발된 방법의 장점은 응답 신호에 유색잡음이 존재하여도 모우드 변수들을 항상 정확하게 구할 수 있으며, 또한 알고리즘의 안정성이 보장된 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.165-165
• /
• 2011

현재 계획 또는 설계 단계에서 수행되고 있는 관거 시설의 수리계산에는 연결관 내에서의 마찰손실만을 감안하여 수행하고 있으며, 맨홀에서의 에너지 손실은 고려되지 않는 실정이다. 그러나 연결관 내부와 맨홀의 내부는 여러 가지 수리학적 조건이 다르므로 에너지 손실이 발생하게 된다(최원석과 송호면, 2002). 더욱이 직선으로 연결된 중간맨홀보다 유입관과 유출관이 $90^{\circ}$의 각도로 접합된 합류맨홀은 연결 구조상 유수교란에 의한 에너지 손실이 커질 것으로 예상됨에도 불구하고 현재 실무에서 우수 배수시설의 설계 시 직선 연결맨홀과 $90^{\circ}$ 접합맨홀의 손실을 구별하지 않고 사용하고 있는 실정이다. 그러므로 $90^{\circ}$ 접합맨홀에서 우수관거 시스템의 우수 배제 능력을 증가시켜 도심지의 침수를 방지하기 위한 관거시설의 적정 설계 기준이 필요하며, 합리적인 설계 기준을 제시하기 위하여 $90^{\circ}$ 접합맨홀 내에서의 수두 손실을 분석할 필요가 있다. 본 연구에서는 수리모형 실험의 물질적 및 시간적 한계를 극복하기 위하여 일반적으로 3차원 유체거동의 특성분석에 많이 사용되는 Fluent 6.3 모형을 이용하여 과부하 $90^{\circ}$ 접합맨홀에서의 흐름특성을 수치모의 하였으며, 맨홀 내 손실수두의 변화를 계산하여 손실계수를 산정하였다. 맨홀 및 접합 관거의 기하 모형의 격자망은 수치해석의 안정성 확보를 위하여 그림 1과 같이 6면체 격자로 구성하였다. 또한 $90^{\circ}$ 접합맨홀에서 급격한 와류에 의해 발생하는 에너지 손실을 저감하기 위하여 $90^{\circ}$ 접합맨홀의 내부 형상 및 접합 조건을 변화시켜 손실계수를 산정하였다. 수치모형의 적용 결과 맨홀 내에서의 유속변화, 수심변화 및 압력변화에 대해서는 수리모형 실험 결과와 유사한 경향을 나타내고 있으며, 수치모형에 의하여 산정된 $90^{\circ}$ 접합맨홀에서 에서의 손실계수 값과 수리모형에 의하여 산정된 손실계수 값이 거의 유사하게 나타났다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.7
• /
• pp.37-42
• /
• 2004

In order to monitor the impact locations in smart structures, multipoint ultrasonic sensors are to be employed. In this study, a multiplexing demodulator with wide dynamic range was proposed to detect the impact locations using FBG sensors, and a stabilization controlling system was also developed for the maintenance of maximum sensitivity of sensors. Two FBG sensors were attached on the bottom side of the aluminum beam specimen and low velocity impact tests were performed to detect the one-dimensional impact locations. As a result, multiplexed in-line FBG sensors could detect the moment of impact precisely, and found the impact locations with the average location error below 0.58mm.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.