• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회(2)
• /
• pp.1437-1439
• /
• 2005

• Coupled systems mechanics
• /
• 제7권1호
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회(1)
• /
• pp.757-759
• /
• 2005

• 한국수산과학회지
• /
• 제36권2호
• /
• pp.178-186
• /
• 2003

Seasonal variation of phytoplankton was investigated with surface mixed layer ecosystem model in the East Sea. The model consisted of four compartments (phytoplankton, zooplankton, nutrient, detritus) forced by mixed layer depths, photosynthetically available radiation and nutrient concentrations. From model results we estimated entrainment rate $2.5-4.0\;m{\cdot}day^{-1}$ to reproduce the two annual blooms, and reproduced seasonal variation of phytoplankton at southern and northern regions by the difference of surface winter mixed layer depth (MLD) using the entrainment rate value $3.0\;m{\cdot}day^{-1}$. The spring blooms in the southern and northern regions closely related to deepening of a winter surface MLD. In the southern region where MLD was shallow and phytoplankton spring bloom occurs one month in advance to the northern region where MLD was deep. The amount of light increases within the MLD during the onset of stratification and water temperature increases faster in spring in the southern region than the northern region. Decrease of phytoplankton was mainly affected by zooplankton grazing in the southern region and by nutrient exhaustion in the northern region. The fall bloom in the two regions was caused by the nutrient availability and entrainment on the phytoplankton.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2003년도 총회 및 춘계학술발표회
• /
• pp.163-167
• /
• 2003

Based on the transient finite difference solution of Richards' equation, an infiltration model is developed to analyze temporal variation of precipitation recharge in the region of monsoon climates. Simulation results obtained by using time series data of 20-year daily precipitation and pan evaporation indicate that a linear relationship between the annual precipitation and the annual recharge holds for the soils under the monsoon climates with varying degrees of the correlation coefficient depending on the soil types. A sensitivity analysis reveals that the water table depth has little effects on the recharge for the sandy soil, whereas, for the loamy and silty soils, rise of the water table at shallow depths causes increase of evaporation by approximately 100㎜/yr and a corresponding decrease in recharge. A series of simulations for two-layered soils illustrate that the amount of recharge is dominantly determined by the soil properties of the upper layer, although the temporal variation of recharge is affected by both layers.

• Water Engineering Research
• /
• 제6권1호
• /
• pp.17-30
• /
• 2005

The bed evolution of the stretch of the River Rhine between km-812.5 and km-821.5 is characterised by general bed degradation as a result of the river training works and dredging activities of the last two centuries. The degradation of the river bed affects the water levels, and so the navigation conditions. To combat the erosion of the river bed with the aim to keep up the shipping traffic and to avoid the ecological system damages due to water level reductions, sand-gravel-mixtures were added to the river (so called artificial grain feeding activities). This paper presents the results of an application of a graded sediment transport model in order to study morpholodynamical characteristics due to artificial grain feeding activities in the river stretch. The finite element code TELEMAC2D was used for flow calculation by solving the 2D shallow water equation on non-structured grids. The sediment transport module SISYPHE has been developed for graded sediment transport using a multiple layer model. The needs to apply such graded sediment transport approaches to study morphological processes in the domain are discussed. The calculations have been carried out for the case of middle water flow and different size-fraction distributions. The results show that the grain feeding process could be well simulated by the model.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2008년도 공동학술대회
• /
• pp.47-50
• /
• 2008

방조제에 의해 조성된 담수호 바닥 퇴적물의 특성 파악을 위하여, 스트리머 케이블을 이용한 수상 쌍극자배열 전기비저항탐사를 수행하였다. 담수호에서 얻어진 전기비저항 탐사 자료는 전극전개수가 작은 경우에는 낮은 겉보기 비저항 값을 보이며, 전극전개수가 증가함에 따라 증가하는 양상을 보인다. 따라서 고정된 전류전극에 대하여 전극 전개수가 증가할 경우 겉보기 비저항값이 크게 감소하는 자료는 역산자료에서 제외하는 기법을 적용하는 것이 필수적인데, 하부층에 대한 상부층의 비저항이 0.6$\sim$0.8 사이의 값을 적용하는 경우 역산 결과가 양호한 것으로 나타났다. 또한 바닥 퇴적물을 2층과 3층으로 설정한 후 역산 결과 담수호 물의 비저항을 고려한 역산 결과가 하부 퇴적물의 비저항 분포를 파악하는데 효과적인 것으로 나타났다.

• 한국해안해양공학회지
• /
• 제15권3호
• /
• pp.159-166
• /
• 2003

본 연구에서는 천해 파랑 계산에 널리 사용되어지고 있는 확장형 완경사 방정식과 계산 효율은 같게 유지하면서 Laplace방정식을 직접 풀 수 있는 유한 요소 모형에 대해서 연구하였다. 기존의 확장형 완경사 방정식을 사용하는 경우와 같은 계산효율을 유지하기 위하여 파동장을 수심방향으로 1층인 유한요소로 나누고, 요소내의 포텐셜을 수면에 위치한 절점에 대한 포텐셜만으로 표시하도록 한 후, Galerkin 기법을 적용하여 수치모형을 구성하였다. 요소 내 수평방향에 대해서는 통상의 보간함수를 채택하였으며, 수심방향에 대해서는 진행파의 수심방향 거동인 함수를 사용하여 보간함수를 구성하였다. 모형의 개발은 우선 연직 2차원 문제를 대상으로 하였다. 개발된 모형의 검증을 위하여 연직 2차원에서의 파랑 반사 및 전달문제에 적용한 결과, 개발된 유한 요소 모형은 계산상의 효율면에서나 해의 정확도 면에서 기존의 확장형 완경사 방정식에 기초한 모형과 같은 수준을 보임을 확인하였다.

• 한국작물학회:학술대회논문집
• /
• 한국작물학회 2017년도 9th Asian Crop Science Association conference
• /
• pp.315-315
• /
• 2017

In Japan, more than 80 % of soybean growing area is converted fields and excess water is one of the major problems in soybean production. For example, recent study (Yoshifuji et al., 2016) suggested that in the fields of shallow groundwater level (GWL) (< 1m depth), rising GWL even in a short period (e.g. 1 day) causes inhibition of soybean growth. Thus it becomes more and more important to predict GWL and soil moisture in detail. In addition to conventional surface drainage and underdrain, FOEAS (Farm Oriented Enhancing Aquatic System), which is expected to control GWL in fields adequately, has been developed recently. In this study we attempted to predict GWL and soil moisture condition at the converted field with FOEAS in Biwa lake reclamation area, Shiga prefecture, near the center of the main island of Japan. Two dimensional HYDRUS model (Simuinek et al., 1999) based on common Richards' equation, was used for the calculation of soil water movement. The calculation domain was considered to be 10 and 5 meter in horizontal and vertical direction, respectively, with two layers, i.e. 20cm-thick of plowed layer and underlying subsoil layer. The center of main underdrain (10 cm in diameter) was assumed to be 5 meter from the both ends of the domain and 10-60cm depth from the surface in accordance with the field experiment. The hydraulic parameters of the soil was estimated with the digital soil map in "Soil information web viewer" and Agricultural soil-profile physical properties database, Japan (SolphyJ) (Kato and Nishimura, 2016). Hourly rainfall depth and daily potential evapo-transpiration rate data were given as the upper boundary condition (B.C.). For the bottom B.C., constant upward flux, which meant the inflow flux to the field from outside, was given. Seepage face condition was employed for the surrounding of the underdrain. Initial condition was employed as GWL=60cm. Then we compared the simulated and observed results of volumetric water content at depth of 15cm and GWL. While the model described the variation of GWL well, it tended to overestimate the soil moisture through the growing period. Judging from the field condition, and observed data of soil moisture and GWL, consideration of soil structure (e.g. cracks and clods) in determination of soil hydraulic parameters at the plowed layer may improve the simulation results of soil moisture.

• 지질공학
• /
• 제27권3호
• /
• pp.305-312
• /
• 2017

The most effective way to distinguish subsurface interfaces that produce various geophysical responses is through the integration of multiple geophysical methods, with each method detecting both a complementary and unique set of distinct physical properties relating to the subsurface. In this study, shallow seismic reflection (SSR) and ground penetrating radar (GPR) surveys were conducted at the Cheongju-Gadeok site of the Korea National Groundwater Monitoring Network to map the water table, which was measured at 12 m depth during the geophysical surveys. The water table proved to be a good target reflector in both datasets, as the abrupt transition from the overlying unsaturated weathered rock to the underlying saturated weathered rock yielded large acoustic impedance and dielectric constant contrasts. The two datasets were depth converted and integrated into a single section, with the SSR and GPR surveys conducted to ensure subsurface imaging at approximately the same wavelength. The GPR data provided detailed information on the upper ~15 m of the section, whereas the SSR data imaged structures at depths of 10-45 m. The integrated section thus captured the full depth coverage of the sandy clay, water table, weathered rock, soft rock, and hard rock structures, which correlated well with local drillcore and water table observations. Incorporation of these two geophysical datasets yielded a synthetic section that resembled a simplified aquifer model, with the best-fitting seismic velocity, dielectric constant, and porosity of the saturated weathered layer being $v_{seismic}=1000m/s$, ${\varepsilon}_r=16$, and ${\phi}=0.32$, respectively.