• Journal of Korea Water Resources Association
• /
• v.52 no.5
• /
• pp.361-372
• /
• 2019

The purpose of this study is to develop a two-dimensional land surface flood analysis model based on uniform square grid using the governing equations except for the convective acceleration term in the momentum equation. Finite volume method and implicit method were applied to spatial and temporal discretization. In order to reduce the execution time of the model, parallel computation techniques using CPU were applied. To verify the developed model, the model was compared with the analytical solution and the behavior of the model was evaluated through numerical experiments in the virtual domain. In addition, inundation analyzes were performed at different spatial resolutions for the domestic Janghowon area and the Sebou river area in Morocco, and the results were compared with the analysis results using the CAESER-LISFLOOD (CLF) model. In model verification, simulation results were well matched with the analytical solution, and the flow analyses in the virtual domain were also evaluated to be reasonable. The results of inundation simulations in the Janghowon and the Sebou river area by this study and CLF model were similar with each other and for Janghowon area, the simulation result was also similar to the flooding area of flood hazard map. The different parts in the simulation results of this study and the CLF model were compared and evaluated for each case. The results of this study suggest that the model proposed in this study can simulate the flooding well in the floodplain. However, in case of flood analysis using the model presented in this study, the characteristics and limitations of the model by domain composition method, governing equation and numerical method should be fully considered.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.4
• /
• pp.3-11
• /
• 2004

A finite volume model is developed to simulate the surface irrigation at a paddy field. The model's capabilities are validated through comparison with the simulafed results and the observed data obtained by various experimental tests, and the simulated results are in good agreement with the observed pending depth. The result of surface irrigation simulation shows that the longer the paddy field's the length of long-sided becomes, the longer the advance and storage time is taken. To analyze surface irrigation performance with variable inflow rate, three patterns of flow variation-constant rate, initially high then low, and initially low then high-were studied. The results show that at the pattern with initially high followed by low during the latter half of the irrigation the advance time is shortest, but the pending depth of irrigation completion and irrigation effiency are the little difference between irrigation patterns.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.761-772
• /
• 2008

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.3
• /
• pp.69-78
• /
• 2007

In this study, a time-dependent aspect of an embankment failure is considered to simulate a flood inundation map and calculate overtopping discharge induced by an embankment failure. A numerical model has been developed by solving the two dimensional nonlinear shallow water equations with a finite volume method on unstructured grids. To analyze a Riemann problem, the HLLC approximate Riemann solver and the Weighted Averaged Flux method are employed by using a TVD limiter and the source term treatment is also employed by using the operator splitting method. Firstly, the numerical model is applied to a dam break problem and a sloping seawall. Obtained numerical results show good agreements with experimental data. Secondly, the model is applied to a flow induced by an embankment failure by assuming that the width and elevation of embankment are varied with time-dependent functions. As a result of the comparison with each numerical overtopping discharge, established flood inundation discharges in the previous studies are overestimated than the result of the present numerical model.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.2
• /
• pp.39-47
• /
• 1997

In the present paper, viscous flow fields around a wigley hull with appendages are analysed to study interactions between the hull and appendages. Navier-Stokes and continuity equations are solved by a finite volume method in a body-fitted coordinate system which conforms three dimensional ship geometries with appendages. A Sub-Grid Scale(SGS) turbulent model is used for a calculation of high Reynolds number flow. Numerical computations has been done for a Wigley hull form at $Rn=1.0{\times}10^6$. The results show that the present approach can predict, at least in qualitative sense, the influence of the appendages upon the flow field around a ship.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.55-63
• /
• 2016

Nowadays, large container ships are continually developed and that's why the bow and stern structural stability problems by slamming become a significant more and more. However, due to the complexity of slamming, it is difficult to consider those problems at the design stage. For this reason, we attempt numerical analysis through SNUFOAM by generating the bow and stern two-dimensional cross-sectional grid in WILS JIP experiment at KRISO. Unlike the conventional method for the computation time saving, by setting the inlet flow conditions referred to the model test, we analyzed the slamming without applying the grid deformation method. As a result, when the stern model, as in the previous studies, it was possible to obtain quantitatively the fluid impulse is close to the experimental results. When the bow model, we can found the change by the position of force sensors which are derived for the bulbous bow and obtained fluid impulse and flow shape at slamming similar to the model test.

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

식생된 개수로에서 식생의 영향을 파악하기 위해 k-$\in$ 난류 모형을 이용하여 수치모의를 하였다. 식생의 영향을 고려하기 위해 항력항을 추가한 지배방정식을 구성하였으며, 지배방정식을 해석하기 위하여 유한체적법을 사용하였다. 수치모의에서 구한 식생된 개수로의 흐름구조를 기존의 수리실험 결과와 비교하여 비교적 잘 일치함을 확인할 수 있다. 난류의 생성과 소멸을 수치모의한 결과, 부분구간 식생된 경우 식생높이 보다 낮은 구간에서는 후류에 의한 난류 생성이 지배적이며, 식생높이보다 높은 구간에서는 주로 마찰에 의한 난류 생성이 지배적임을 보였다. 기존의 연구들은 식생의 영향을 고려하여 개수로의 흐름을 연구한 예는 드물며, 현재까지 진행되어진 국내의 연구는 난류모형을 이용하여 식생된 개수로에서의 흐름 구조를 모의하였다. 따라서 난류흐름을 모의하는데 가장 보편적인 k-$\in$ 난류모형을 이용하여 식생된 개수로에서 수직방향으로의 흐름구조와 식생의 영향을 해석하는 것은 그 자체로도 의미 있는 연구이며, 앞으로의 환경수리 문제를 해결하기 위해 선행되어야 하는 연구이다. 식생된 개수로에서의 난류구조와 부유사 이동에 대한 식생의 영향을 비정상 1차원 수직모형으로 해석하였으며, 폐합문제를 위해 2-방정식인 k-$\in$ 난류모형을 사용하였다. k-$\in$ 난류모형에 식생에 의한 항력항을 더하여 지배방정식을 구성하였다. 수직방향에 대해 흐름방향 유속 u, 난류에너지 k, 그리고 난류에너지 소산율 $\in$의 분포를 구하고, 부유사에 대한 수송방정식을 풀었다. 식생된 개수로와 식생되지 않은 개수로에서의 유속분포, 난류강도, 레이놀즈 응력 분포와 난류의 생성과 소멸을 구하여 식생이 난류흐름에 미치는 영향을 분석하였다.

• Journal of Korea Water Resources Association
• /
• v.42 no.12
• /
• pp.1079-1089
• /
• 2009

The objective of this study is to develop the accurate, robust and high resolution two-dimensional numerical model that solves the computationally difficult hydraulic problems, including the wave front propagation over dry bed and abrupt change in bathymetry. The developed model in this study solves the conservative form of the two-dimensional shallow water equations using an unsplit finite volume scheme and HLLC approximate Riemann solvers to compute the interface fluxes. Bed-slope term is discretized by the divergence theorem in the framework of FVM for application of unsplit scheme. Accurate and stable SGM, in conjunction with the MUSCL which is second-order-accurate both in space and time, is adopted to balance with fluxes and source terms. The exact C-property is shown to be satisfied for balancing the fluxes and the source terms. Since the spurious oscillations in second-order schemes are inherent, an efficient slope limiting technique is used to supply TVD property. The accuracy, conservation property and application of developed model are verified by comparing numerical solution with analytical solution and experimental data through the simulations of one-dimensional dam break flow without bed slope, steady transcritical flow over a hump and two-dimensional dam break flow with a constriction.

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

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.3
• /
• pp.3-14
• /
• 2004

A two-dimensional numerical model based on a finite volume method was formulated to solve the shallow water equations and applied for evaluating drainage characteristics at large-sized paddy fields. Manning roughness coefficient was calibrated using the observed inundating depths at drainage tests, and used for validating the model with the results from another drainage test. The simulated results were in good agreement with the observed inundating depths. The result of surface drainage showed that the longer width of the outlet was or the more the number of drainage outlet was, the shorter the drainage time was taken, and the larger the size of the field become, the longer the drainage time was taken, and the field shape had little effect on drainage time. To reduce the drainage time to 24 hours, the outlet is located lower than the elevation of the basin and small drainage ditch is constructed at the field. The results showed that the drainage time was taken short as the small drainage ditch was constructed. The comparison of drainage time as to the size of field constructed small drainage ditch showed the field, 100m ${\times}$ 200 m, can be drained in 24 hours.