• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.822-826
• /
• 2009

홍수터에 식재된 수목에 의한 운동량 교환은 흐름의 추가적인 흐름저항으로 작용한다. 따라서 수목으로 인한 통수능의 감소로 부분적인 유속 상승이 발생하기도 하지만, 전반적으로 평균유속이 감소하며 수위는 상승하게 된다. 실무적 사용목적으로 홍수터 수목의 저항은 통합 매개변수인 합성조도계수가 사용되고 있으나, 정량적인 조도계수의 추정은 매우 어려운 문제이다. 본 연구에서는 경험적으로 추정되는 비정량적 합성조도 계수의 한계를 극복하고, 실무적으로 사용이 용이하도록 수목에 의한 흐름저항 항을 분리한 1차원 상류이송 음해기법이 개발되었다. 개발된 기법의 평가를 위하여 균일한 단면에서 수행된 수리실험에 적용하였으며, 다양한 유량규모 및 식재조건에서 제안된 기법의 적용성 및 정확성을 확인하였다. 불규칙한 단면형상을 가진 실제 하천의 홍수터에 식재된 수목군의 경우를 가정하여 1/25 정상축척의 수리실험을 실시하였고, 100년 빈도 홍수량에서 수목식재 전후의 실측결과를 개발된 수치기법으로 계산하여 비교하였다. 개발된 기법은 모든 경우에서 사용자의 주관적인 판단 없이 정확하면서 안정적인 해를 제공하였으므로, 흐름의 2차원 효과가 상쇄될 수 있는 광폭의 자연하천에서 수목의 식재기준 분석을 위한 대안으로 사용될 수 있을 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.604-604
• /
• 2012

국지성 집중호우에 따른 도심지 내수 침수 피해의 주원인으로 하수관거의 설계기준을 초과하는 강우가 침수피해의 주요 원인이며, 도심화로 인해 불투수 면적이 증가함에 따라 유출되는 시간이 짧아 저지대의 피해는 불가피하다. 2010년과 2011년에 100년 이상의 강우사상이 서울시에 연이어 나타나면서 집중호우로 인한 피해지역이 유사하게 나타났으며, 광화문 거리의 연이은 침수는 현재 서울시의 하수관거의 용량과 빗물펌프장 및 저류조 시설로 구성된 기존 수방대책의 한계점을 보이고 있다. 이에 서울시는 광화문 일대의 배수능력을 향상시키기 위하여 효자배수분구 빗물배수터널을 계획하고 있다. 일본, 미국 및 유럽 등지에서는 대심도 지하수로 시설에 대한 수리실험 및 수치 연구를 바탕으로 다양한 지하방수로가 건설되어 국지성 집중 강우에 대해 적절히 대응하고 있으나, 국내의 경우에는 대심도 지하방수로 시설에 대한 연구가 미비하여 지하방수로 설계 지침 및 기술적 자료가 부족한 실정이다. 그러므로 대심도 빗물배수터널 시설에서의 흐름특성 분석에 관한 수리실험 및 수치해석 등의 구체적인 연구가 필요하다고 판단된다. 본 연구에서는 수리모형 실험의 물질적 및 시간적 한계를 극복하기 위하여 일반적으로 3차원 유체거동의 특성분석에 많이 사용되는 Fluent 6.3 모형을 이용하여 대심도 빗물배수터널 시설의 접선식 유입구에 대한 흐름특성을 수치모의 하였다. 접선식 유입구 및 수직갱(drop shaft)에 대한 기하 모형의 격자망은 수치해석의 안정성 확보를 위하여 그림 1과 같이 6면체 격자로 구성하였다. 맨홀 내의 다상유동을 고려하기 위하여 VOF(Volume of Fluid) Scheme을 적용하였으며, 수치해석 방법으로는 비정상류, 1st order implicit method를 사용하였다. Fluent에서의 난류 흐름을 계산하는 방법에는 난류 운동에너지와 난류 에너지 소산율 $\epsilon$의 전달 방정식을 도입한 k-$\epsilon$ 난류 모형을 채택하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.103-112
• /
• 1991

In this study, a numerical model was established and applied to simulate the steady-state groundwater and heat flow in an isotropic, heterogeneous, three dimensional aquifer system with uniform thermal properties and no change of state. This model was developed as an aid in screening large groundwater-flow systems as prospects for underground waste storage. Driving forces on the system are external hydrologic conditions of recharge from precipitation and fixed hydraulic head boundaries. Heat flux includes geothermal heat-flow, conduction to the land surface, advection from recharge, and advection to or from fixed-head boundaries. The model uses an iterative procedure that alternately solves the groundwater-flow and heat-flow equations, updating advective flux after solution of the groundwater-flow equation, and updating hydraulic conductivity after solution of the heat-flow equation. Dierect solution is used for each equation. Travel time is determined by particle tracking through the modeled space. Velocities within blocks are linear interpolations of velocities at block faces. Applying this model to the groundwater-flow system located in Jigyung-ri. Songla-myun, Youngil-gun. Kyungsangbuk-do, the groundwater-flow system including distribution of head, temperature and travel time and flow line, is analyzed.

• Journal of the Korean Geotechnical Society
• /
• v.26 no.4
• /
• pp.47-58
• /
• 2010

Seepage rate through the core zone of rockfill dam, estimated from graphical technique and the equation by Sakamoto (1998), is different from the real condition because of neglecting unsaturated flow. With existing method to estimate total seepage rate, it is difficult to understand the tendency of total seepage rate changes by reservoir water level change. Steady state seepage rate and the factors affecting the time needed to attain to changes of reservoir water level and saturated hydraulic conductivity and unsaturated hydraulic properties of core material are analysed thorough the 2-D steady and unsteady state seepage analyses of Soyanggang dam. Numerical results revealed that the seepage rate can be expressed by the linear equation form and the value of unsaturated soil parameter n is the most important factor affecting the seepage rate and the time needed to attain steady state. The estimation method presented in this study can be used by the designer and the personnel of dam safety for convenient estimation of seepage rate and quantitative analysis of measured seepage rate without 2-D and 3-D numerical analyses.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.37-45
• /
• 1997

벽면 근처에 가까이 위치하는 회전와류쌍을 음원으로 갖는 비정상 유동장에서 벽면이 음장에 미치는 효과를 알아보기 위해 이차원 음장 수치해석을 시도하였다. 비압축성 유동장에 대한 비정상 수력정보를 기반으로 오일러식에서 교란 압축성 소음항을 도출하였다. 원거리 자유 경계면은 비반사 경계조건을 이용하였으며, 벽면에서는 벽면 효과를 음향장에 고려하였다. 자유흐름장에 놓인 와류쌍이 대칭인 나선팔을 갖는 반면에, 벽면이 있는 경우엔 음파가 전달되는 경로를 따라 방향성이 존재함을 알 수 있었다. 본 연구를 통하여 벽면이 존재하는 경우에 비정상 수력정보를 이용하여 근거리와 원거리 음장을 동시에 수행할 수 있음을 알아내었다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.3B
• /
• pp.293-303
• /
• 2011

An efficient diagonalized approximate factorization algorithm (DAF) is developed for the solution of three-dimensional incompressible viscous flows. The pressure-based, artificial compressibility (AC) method is used for calculating steady incompressible Navier-Stokes equations. The AC form of the governing equations is discretized in space using a second-order-accurate finite volume method. The present DAF method is applied to derive a second-order accurate splitting of the discrete system of equations. The primary objective of this study is to investigate the computational efficiency of the present DAF method. The solutions of the DAF method are evaluated relative to those of well-known four-stage Runge-Kutta (RK4) method for fully developed and developing laminar flows in curved square ducts and a laminar flow in a cavity. While converged solutions obtained by DAF and RK4 methods on the same computational meshes are essentially identical because of employing the same discrete schemes in space, both algorithms shows significant discrepancy in the computing efficiency. The results reveal that the DAF method requires substantially at least two times less computational time than RK4 to solve all applied flow fields. The increase in computational efficiency of the DAF methods is achieved with no increase in computational resources and coding complexity.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.3
• /
• pp.146-159
• /
• 2016

This study used Navier-Stokes Solver(LES-WASS-2D) for analyzing hydrodynamic characteristics with high order in order to analyze self-burial mechanism of pipeline with spoiler under steady flow. For the validity and effectiveness of numerical model used, it was compared and analyzed with the experiment to show flow characteristics around the pipeline with and without the spoiler. And the hydraulic(flow, vortex, and pressure) and force characteristics were numerically analyzed around the pipeline according to the incident velocity, and shape and arrangement of spoiler. Primarily, if the spoiler is attached to the pipeline, the projected area is increased resulting in higher flow velocity toward the back and strong vortex caused by wake stream in the back. Secondly, the spoiler causes vertically asymmetric flow and vorticity fields and thus asymmetric pressure field. It increases the asymmetry of force on the pipe and thus develops large downward fluid force. Both of them are the causes of selfburying of the pipeline with spoiler.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.11
• /
• pp.73-82
• /
• 2017

Recently, offshore bucket cut-off walls were developed to solve several problems in conventional offshore cut-off walls. In this study, a numerical analysis was carried out to investigate the seepage stability of offshore bucket cut-off walls. The ground was assumed as uniform homogeneous sand and steady state flow conditions were applied. The flow condition was compared among 2-dimensional flow (2-D), 2-dimensional concentrated flow (2-DC), and axisymmetric flow. The analysis results showed that the seepage velocities in axisymmetric flow were about 1.5 and 2 times larger than those of 2-DC and 2-D flow conditions, respectively. Thereafter, the axisymmetric flow condition was applied because the seepage flow was concentrated toward the center of the circular-shaped wall. A parametric study was performed varying bucket radius, penetration depth, total head difference between in and outside of the wall. The exit gradient, which used for the calculation of piping stability, decreased with increase of the penetration depth and bucket radius. Design charts were proposed to estimate the factor of safety and the exit gradient at various analysis conditions. Finally, the design equation was proposed to calculate the exit gradient for the preliminary design of the bucket cut-off wall.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.12
• /
• pp.1097-1106
• /
• 2011

Unsteady numerical analysis was performed to simulate air-launched missiles from a complete helicopter in hover by using an unstructured overset mesh flow solver coupled with a module of six degree-of-freedom motion of equations. The unsteady computations have been performed to obtain flow fields around the complete helicopter including main rotor, tail rotor, and fuselage equipped with multiple missiles, and six-DOF simulation has been performed to predict the behavior of the air-launched missile. The effects of the launching position and the missile thrust on the trajectory of the missile were investigated as well as the aerodynamic interference of the air-launched missile under the unsteady downwash produced by main rotor.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.249-262
• /
• 2011

In this study, in order to reduce the numerical oscillation due to the unbalance between source and flux terms as the HLLC scheme is applied to the flow analysis on the irregular bed topography, a unstructured finite volume model based on the well-balanced HLLC scheme and the shallow water equations is developed and applied to problems of dam-break waves. The well-balanced HLLC scheme considers directly the gradient of bed topography as the flux terms is calculated. This scheme provides the good numerical balance between the source and flux terms in the case of the application to the steady-state transcritical flow. To verify the numerical model developed in this study, it is applied to three cases of hydraulic model experiments and a field case study of Mapasset dam failure (France). As a result of the verification, the predicted numerical results agree relatively well with available laboratory and field measurements. The model provides slightly more accurate results compared with the existing models.