• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.5
• /
• pp.461-471
• /
• 2008

An analytical solution of one dimensional mild slope equation is derived by use of the WKB method, which has a form similar to Porter's solution(2003). The present solution is so general in the sense of application that it is comparable to the corresponding numerical solutions. In the derivation we also presented the solution of refraction equation in terms of surface displacement. Some numerical results of the present solution by use of Bremmer's method are presented which agree with existing numerical solutions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• Tunnel and Underground Space
• /
• v.21 no.3
• /
• pp.216-224
• /
• 2011

In this study, Class II behavior of rock failure process under uniaxial and biaxial compression has been numerically simulated using bonded particle model. Class II behavior of rock was simulated by radial strain controlled uniaxial and biaxial compression tests using a suggested method of ISRM. Micro-parameters used in the simulation were determined based on the laboratory uniaxial compression tests carried out at ${\"{A}}sp{\"{o}}$ Hard Rock Laboratory, Sweden. Class II behavior of ${\"{A}}sp{\"{o}}$ rock was effectively simulated using newly proposed numerical technique in this study, and the results of numerical simulations show good similarity with the complete stress-strain curves for Class II behavior obtained from the laboratory tests.

• Journal of the Korean Geotechnical Society
• /
• v.17 no.6
• /
• pp.245-255
• /
• 2001

Existing cut-and-cover tunnels are designed regardless of cut-slope under the assumption that the overburden weight of backfill soil acts on tunnel arch and the earth pressure at rest acts on tunnel walls. However, actual earth pressures acting on the tunnel lining depend on open-cut size composed of cut-slope and cut-width, and thus the tunnel lining shows a different structural behavior. This study investigated the effect of cut-slope on structural behavior of the cut-and-cover tunnel lining as follows; Firstly, a comprehensive numerical analysis method using FLAC2D code was used and verified by field measurements of tunnel profile. Secondly, based on the verified numerical analysis technique, earth pressure acting on the lining, and displacement and sectional force developed on the lining were estimated with various shapes of cut-slopes$30^{\circ}\;, 456{\circ},\; 60^{\circ},\; and\;75^{\circ}%). Numerical analysis results indicate that the steeper cut-slope shows the more displacement and moment of the tunnel lining.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.1
• /
• pp.133-138
• /
• 2008

To analyze incident waves traveling from the deep ocean is very important in that it is based on resolving problems occurred in coastal areas. In general, numerical models and analytical solutions are used to analyze wave transformation. Although a numerical model can be applied to various bottoms and wave conditions, it may have some cumbersome numerical errors. On the other hand, an analytical solution has an advantage of obtaining the solution quickly and accurately without numerical errors. The analytical solution can, however, be utilized only for specific conditions. In this study, the analytical solution of the mild-slope equation has been developed. It can be applied to various conditions combing a numerical technique and an analytical approach while minimizing the numerical errors. As a result of comparing the obtained solutions in this study with those of the previously developed numerical model, A good agreement was observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.4
• /
• pp.359-370
• /
• 2015

The present study aims to numerically investigate the behavior of liquid droplet driven by capillarity force imbalance on horizontal surfaces ranging from hydrophilic to hydrophobic, under various conditions. The droplet behavior has been simulated using an in-house solution code(PowerCFD), which employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The detailed droplet behavior was obtained under various conditions for droplets with different initial shapes, contact angles and surface tension forces(or Bond number). The mechanism of droplet transport was examined using the numerical results on the droplet shapes.

• The Journal of the Acoustical Society of Korea
• /
• v.25 no.5
• /
• pp.229-238
• /
• 2006

This paper shows the analytic error caused by the inconsistency of the approximation order between the non local boundary condition (NLBC) and the parabolic governing equation. To obtain the analytic error, we first transform the NLBC to the half space domain using plane wave analysis. Then, the analytic error is derived on the boundary between the true numerical domain and the half space domain equivalent to the NLBC. The derived analytic error is physically expressed as the artificial reflection. We examine the characteristic of the analytic error for the grazing angle, the approximation order of the PE or the NLBC. Our main contribution is to present the analytic method of error estimation and the application limit for the high order parabolic equation and the NLBC.

• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.751-760
• /
• 2003

Error analysis of finite difference equation on the Muskingum-Cunge flood routing method with free time and space weighting factor was carried out. The error analysis shows that the numerical solution of the Muskingum-Cunge method becomes diverged with time when the sum of time weighting factor and space weighting factor is greater than 1.0. Numerical diffusion increases when the sum of time weighting factor and space weighting factor decreases. Numerical diffusion and numerical oscillation increase when the grid resolution is coarse. Numerical experiments and field applications show that the Muskingum-Cunge method with free space weighting factor is more effective for simulating the flood routing with great peak diminution than conventional Muskingum-Cunge method with fixed space weighting factor, 0.5.

• Water for future
• /
• v.21 no.3
• /
• pp.291-291
• /
• 1988

The two dimensional time dependent flow past a circular cylinder is analyzed numerically. In the analysis, equations of conservation of mass and momentum are transformed to equations of stream function-vorticity and vorticity transport, and nondimensionalized by nondimensional parameters representing flow characteristics, The resulting stream function-vorticity equstion and vorticity transport equation are solved by successive over relaxation scheme and alternating direct implicit scheme. Numerical experments are performed for the flow in the range of Reynolds number 125 to 275. The time dependent streamlines, vorticities, pressure on cylinder surface, separation angle, and drag and lift coefficients are calculated, and the method for estimation of pressure on cylinder surface and the outer boundary limit are developed.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.