• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.64-64
• /
• 2016

하천은 인간에게 용수의 이용 및 하천호안의 휴식처로써의 이용을 통해 직접적인 영향을 주고, 하천구조물의 심미적 영향, 랜드마크로써의 역할을 통해 간접적인 영향을 준다. 또한, 하천은 하천생태계에 서식하는 동 식물에게 영향을 준다. 그러나 하천유사로 인해 통수능이 감소하고, 하천구조물 주변에 침식을 야기할 뿐만 아니라, 댐과 저수지에 유사의 퇴적으로 저수용량의 감소시킨다. 그러므로 이를 예측하는 것은 경제적, 환경적으로 중요하다. 하상변동의 모의를 위해 기존의 2차원 모형은 만곡흐름에서 유동의 helical flow를 고려하지 않아 예측이 부정확하였다. 본 연구에서는 천수방정식을 이용한 하상변동 수치모의에 helical flow의 영향을 고려하였다. 하천과 같은 천수영역에서의 흐름 및 하상변동을 해석하기 위하여 수심평균 된 Navier-Stokes equations인 천수방정식을 이용하였다. 지배방정식은 곡선 좌표계에서 유한체적법으로 차분하였고, 비엇갈림격자를 사용하였다. 지배방정식의 닫힘 문제를 해결하기위해 0-방정식 난류모형을 사용하였고, "time marching" 기법의 적용을 위해 계산단계분할 방법을 이용하였다. 비엇갈림격자의 사용으로 인해 검사체적의 면에서의 유속이 필요하여 pressure-velocity coupling을 사용하여 유속의 진동을 줄였다. 또한, 만곡부의 helical flow를 모의하기위해 helical flow intensity model을 도입하였다. 앞에서 계산한 흐름을 바탕으로 유사량 산정공식과 Exner 방정식을 이용하여 하상변동을 모의하였다. 흐름의 검증, helical flow의 영향에 대한 확인, 하상변동의 적용을 위해 선행연구의 실험이 사용되었다.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.70-76
• /
• 1997

수치격자의 크기제한에 의한 자유표면 유동해석 문제를 효율적으로 다루기 위하여 자유표면의 모든 격자를 x 방향으로 4,8,12개로 등분할 하고, y 방향으로는 4개로 잘라서 계산하였다. 이중격자 또는 삼중격자로 Navier-Stokes 방정식의 각항에 크기가 다른 격자를 사용해 효율을 향상시키는 계산방법의 연장으로, 본 논문에서는 자유표면 방정식에 보다 세분화된 격자를 적용해, Marker Particle 이동 및 자유수면 형성에 효율향상을 줄수 있는 수치방법을 도입하였다. 계산결과에 의하면 초기사용 격자가 coarse한 경우가 본방법의 효과가 커짐을 알 수 있고 대상물로는 층류유동에서 Wigley모형과 낭류유동의 S103 모형이다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.813-820
• /
• 2014

This study presents a numerical modeling of mean flow and turbulence structures of partly-vegetated open-channel flows. For this, Reynolds-averaged Navier-Stokes equations with vegetation drag terms are solved numerically using the non-linear k-${\varepsilon}$ model. The numerical model is applied to laboratory experiments of Nezu and Onitsuka (2001), and simulated results are compared with data from measurement and computations by Kang and Choi's (2006) Reynolds stress model. The simulation results indicate that the proposed numerical model simulates the mean flow well. Twin vortices are found to be generated at the interface between vegetated and non-vegetated zones, where turbulence intensity and Reynolds stress show their maximums. The model simulates the pattern of the Reynolds stress well but under-predicts the intensity of Reynolds stress slightly.

• Journal of Korea Water Resources Association
• /
• v.56 no.10
• /
• pp.603-617
• /
• 2023

This study aims to provide a better understanding of the turbulent flow characteristics in swash zone. A double dam-break method is employed to generate the swash zone flow. Comparing with the conventional single dam-break method, a delay between two gate opening can be controlled to reproduce various interactions between uprush and backwash. For numerical simulations, overInterDyMFoam based on OpenFOAM is adopted. Using overInterDyMFoam, interface between two immiscible fluids having different densities (i.e., air and water phases) can be tracked in a moving mesh with multiple layers. Two-dimensional Reynolds-Averaged Navier-Stokes equations are solved with a standard 𝜅-𝜖 turbulence model for momentum and continuity. Numerical model results are validated with laboratory experiment data for the time series of water depth and streamwise velocity. Turbulent kinetic energy distribution is further investigated to identify the turbulence evolution for each flow regime (i.e., uprush, backwash, and swash-swash interaction).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.42-42
• /
• 2015

사행수로에서의 원심력은 비선형적인 압력분포를 야기하여 이차류, 편수위 등과 같은 불규칙하고 복잡한 흐름을 발생시킨다. 일반적으로 이들 흐름은 난류이고 매우 3차원적이며 자유수면과의 상호작용이 중요한 역할을 할 수도 있다. 환경, 유사이동, 지형 변화와 관련된 환경 수리학적 관점에서 사행수로에서의 흐름을 이해하고 설계하기 위해서는 이러한 복잡한 3차원 흐름을 정확하게 계산하는 것이 매우 중요하다. 이 연구에서는 유한차분법에 근거한 3차원 흐름해석 모형을 이용하여 사행수로에서의 난류 흐름을 모의하고자 한다. 지배방정식은 3차원 비정상 RANS(Reynolds averaged Navier-Stokes) 방정식이며, 난류 해석을 위해서 공학적으로 널리 이용되고 있는 난류 모형 중 k-omega 모형을 이용한다. 수치모형은 시간과 공간에 대해서 2차 정확도의 이산화 기법을 적용한다. 자유수면의 변동은 이상(two-phase) VOF (volume of fluid) 기법을 이용하여 계산한다. 수치모형의 적용 대상은 기존 문헌에서 제시되어 있는 키노시타 사인곡선을 이용하여 만든 폭 60cm의 사행수로에서 후르드수 0.23 그리고 레이놀즈수 41,700의 조건에서 발생시킨 난류 흐름이다. 적용한 난류모형들을 이용하여 해석한 결과들을 유속벡터분포와 수위의 항으로 비교분석하여 사행수로에서 발생되는 이차류와 편수위 변화 재현에 대한 수치모형의 적용성을 평가하고 각 난류모형들의 특성을 제시한다.

• 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 Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.3
• /
• pp.114-122
• /
• 2018

This paper introduces a non-hydrostatic wave model SWASH for simulating wave interactions with porous structures. This model calculates the flow in porous media based on volume-averaged Reynolds-averaged Navier-Stokes equations (VARANS) in ${\sigma}$-coordinate. The empirical coefficients of resistance used to account for the flow in a porous media often need to be measured or calibrated. In this study, the empirical resistance coefficients used in the model are calibrated and validated using laboratory experiments, involving dam-break flow through porous media, and solitary wave interactions with a porous structure. It is shown that the agreement between experimental and numerical results is generally satisfactory. It is also confirmed that non-hydrodynamic model, SWASH, is computationally much more efficient than the three-dimensional porous flow models based on VOF approach.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.

• Journal of the Korean Chemical Society
• /
• v.42 no.4
• /
• pp.398-410
• /
• 1998

For non-linear solution of the Boltzmann equation, the cumulant moment method has been studied. To apply the method to the normal shock wave problem, we restricted ourselves to the monatomic Maxwell molecular gases. The method is based on the iterative approach developed by Maxwell-Ikenberry-Truesdell (MIT). The original MIT approach employs the equilibrium distribution function for the initial values in beginning the iteration. In the present work, we use the Mott-Smith bimodal distribution function to calculate the initial values and follow the MIT iteration procedure. Calculations have been carried out up to the second iteration for the profiles of density, temperature, stress, heat flux, and shock thickness of strong shocks, including the weak shock thickness of Mach range less than 1.4. The first iteration gives a simple analytic expression for the shock profile, and the weak shock thickness limiting law which is in exact accord with the Navier-Stokes theory. The second iteration shows that the calculated strong shock profiles are consistent with the Monte Carlo values quantitatively.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6B
• /
• pp.379-385
• /
• 2012

Recent development of computing power and theoretical advances in computational fluid dynamics have made possible numerical simulations of water waves with full Navier-Stokes equations. In this study, an internal wave maker using the mass source function approach was combined with the level set finite element method for generation of waves. The model is first applied to the two-dimensional linear wave generation and propagation. Then, it is applied to the three-dimensional simulation of the same problem. To effectively utilize computational resources and enhance the speed of execution, parallel algorithms are developed and applied for the three-dimensional problem. The results of numerical simulations are compared with theoretical values and good agreements are observed.