• Applied Biological Chemistry
• /
• v.43 no.2
• /
• pp.116-123
• /
• 2000

A numerical model of three-dimensional soil water distribution for drip irrigation management under cropped conditions was developed using Richards equation in Cartesian coordinates. The model accounts for both seasonal and diurnal changes in evaporation and transpiration, and the growth of plant root and the shape of root zone. Solutions were numerically approximated using the Crank-Nicolson implicit finite difference technique on the block-centered grid system and the Gauss-Seidel elimination in tandem. The model was tested under several conditions to allow the flow rates and configurations of drip emitters vary. In general, simulation results agreed well with experimental results and were as follows. The velocity of soil-water flow decreased drastically with distance from the drip source, and the rate of expansion of the wetted zone decreased rapidly during irrigation. The wetting front of wetted zone from a surface drip emitter traveled farther in vertical direction than in horizontal direction. Under this experimental weather condition, water use efficiency of a drip-irrigated apple field was greatest for 4-drip-emitter system buried at 25 cm, resulting from 10% increase in transpiration but 20% reduction in soil evaporation compared to those for surface 1-drip emitter system. Soil moisture retention curve obtained using disk tension infiltrometer showed significant difference from the curve obtained with pressure plate extractor.

• 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.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.20-26
• /
• 2019

Where hydraulic structures are to be installed over the entire width of a river or stream, usually a bed protection structure is to be installed. However, a local scour occurs in which the river bed downstream of the river protection system is eroded due to the influence of the upstream flow characteristics. This local scour is dominant in the flow and turbulence characteristics at the boundary of the flow direction and in the material of the bed materials, and may gradually become dangerous over time. Therefore, in this study, we compared the turbulence patterns in the local scour hole at the downstream of the river bed protection with the results of the analysis of the mobile bed experiment, and compared with the application of OpenFoam, a three dimensional numerical analysis model. The distribution of depth-averaged relative turbulence intensities along the flow direction was analyzed. In addition to this result, the stabilization of scour hole was compared with the bed shear stress and Shields parameter, and the results were compared by changing the initial turbulent flow conditions. From the results, it was confirmed that the maximum depth of generation of the three-stage was dominantly developed by the magnitude of depth-averaged relative turbulence intensity rather than the mean flow velocity. This result also suggests that design, construction or gate control are needed to control the depth-averaged relative turbulence intensities in order to reduce or prevent the local scour faults that may occur in the downstream part of the bed protection.

• Journal of Korean Society of Disaster and Security
• /
• v.14 no.2
• /
• pp.61-75
• /
• 2021

Recently, as the occurrence frequency of sudden floods due to climate change increased and the aging of the existing spillway, it is necessary to establish a plan to utilize an auxiliary spillway to minimize the flood damage of downstream rivers. Most studies have been conducted on the review of flow characteristics according to the operation of auxiliary spillway through the hydraulic experiments and numerical modeling. However, the studies on examination of flood damage in the downstream rivers and the stability of the revetment according to the operation of the auxiliary spillway were relatively insufficient in the literature. In this study, the stability of the revetment on the downstream river according to the outflow conditions of the existing and auxiliary spillway was examined by using 3D numerical model, FLOW-3D. The velocity, water surface elevation and shear stress results of FLOW-3D were compared with the permissible velocity and shear stress of design criteria. It was assumed the sluice gate was fully opened. As a result of numerical simulations of various auxiliary spillway operations during flood season, the single operation of the auxiliary spillway showed the reduction effect of maximum velocity and the water surface elevation compared with the single operation of the existing spillway. The stability of the revetment on downstream was satisfied under the condition of outflow less than 45% of the design flood discharge. However, the potential overtopping damage was confirmed in the case of exceeding the 45% of the design flood discharge. Therefore, the simultaneous operation with the existing spillway was important to ensure the stability on design flood discharge condition. As a result of examining the allocation ratio and the total allowable outflow, the reduction effect of maximum velocity was confirmed on the condition, where the amount of outflow on auxiliary spillway was more than that on existing spillway. It is because the flow of downstream rivers was concentrated in the center due to the outflow of existing spillway. The permissible velocity and shear stress were satisfied under the condition of less than 77% of the design flood discharge with simultaneous operation. It was found that the flood damage of downstream rivers can be minimized by setting the amount allocated to the auxiliary spillway to be larger than the amount allocated to the existing spillway for the total outflow with simultaneous operation condition. However, this study only reviewed the flow characteristics around the revetment according to the outflow of spillway under the full opening of the sluice gate condition. Therefore, the various sluice opening conditions and outflow scenarios will be asked to derive more efficient utilization of the auxiliary spillway in th future.

• Journal of Korean Society of Disaster and Security
• /
• v.14 no.4
• /
• pp.47-60
• /
• 2021

Recently, as the occurrence frequency of sudden floods due to climate variability increased, the damage of aging chuteway slabs of spillway are on the rise. Accordingly, a wide array of field survey, hydraulic experiment and numerical simulation have been conducted to find the cause of damage on chuteway slabs. However, these studies generally reviewed the flow characteristics and distribution of pressure on chuteway slabs. Therefore the derivation of damage on chuteway slabs was relatively insufficient in the literature. In this study, the cavitation erosion and hydraulic jacking were assumed to be the causes of damage on chuteway slabs, and the phenomena were reproduced using 3D numerical models, FLOW-3D and COMSOL Multiphysics. In addition, the cavitation index was calculated and the von Mises stress by uplift pressure distribution was compared with tensile and bending strength of concrete to evaluate the possibility of cavitation erosion and hydraulic jacking. As a result of numerical simulation on cavitation erosion and hydraulic jacking under various flow conditions with complete opening gate, the cavitation index in the downstream of spillway was less than 0.3, and the von Mises stress on concrete was 4.6 to 5.0 MPa. When von Mises stress was compared with tensile and bending strength of concrete, the fatigue failure caused by continuous pressure fluctuation occurred on chuteway slabs. Therefore, the cavitation erosion and hydraulic jacking caused by high speed flow were one of the main causes of damage to the chuteway slabs in spillway. However, this study has limitations in that the various shape conditions of damage(cavity and crack) and flow conditions were not considered and Fluid-Structure Interaction (FSI) was not simulated. If these limitations are supplemented and reviewed, it is expected to derive more efficient utilization of the maintenance plan on spillway in the future.

• Journal of Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1023-1035
• /
• 2015

This study applies 3-D N-S solver based on PBM (Porous Body Model), LED-WASS-3D ver 2.0 to directly analyze non-linear interaction of seawater-freshwater-coastal aquifer in order to simulate the seawater infiltration into coastal aquifer. This numerical simulation is the first trial in Korea, as well as unusual and new numerical analysis abroad. Firstly, to validate the applied numerical model, the validity and effectiveness was verified for the numerical model by comparing and considering it with the result of laboratory experiment for seawater-freshwater interface in coastal aquifer. And then it simulated the seawater infiltration into coastal aquifer considering the changed levels of seawater and groundwater in order to analyze the distribution characteristics of flow field and seawater-freshwater interface of coastal aquifer as the level difference between seawater and groundwater and rate of seawater level (${\Delta}h/h$) increased. In addition, the characteristics of seawater infiltration were analyzed from the vertical salinity in the coastal aquifer by ${\Delta}h/h$, which cannot be obtained from existing non-diffusion numerical models. Finally, it analyzed the effect of ${\Delta}h/h$ on the seawater infiltration distance in coastal aquifer, which was indexed.