• Water for future
• /
• v.28 no.5
• /
• pp.129-140
• /
• 1995

In this study we derived steady-state relationships between surface saturation area and subsurface outflow, and between surface saturation area and soil moisture storage through numerical experiments with Richards equation on a hillslope. Numerical experiments analyzed the sensitivity of topographic and soil hydraulic properties on steady-state relationships between surface saturation area and subsurface outflow. And the power law for the extent of surface saturation area was determined as a function of subsurface outflow or soil moisture storage.

• Korean Journal of Soil Science and Fertilizer
• /
• v.36 no.4
• /
• pp.181-192
• /
• 2003

We developed a mathematical simulation model to portray the vertical distribution of soil water from the measured weather data and the known soil hydraulic properties, and then compared simulation results with the periodically measured soil water profiles obtained on Jungdong sandy loam to verify the model, In this model, we solved potential-based Richards' equation by the implicit finite difference method superimposed on the predictor-corrector scheme. We presumed that: soil hydraulic properties are homogeneous; soil water flows isothermally; hysteresis is not considered; no vapor flows; no heat transfers into the soil profiles; and water added to soil surface is distributed along the soil profile following partial displacement principle. The input data were broadly classified into two groups: (1) daily weather data such as rainfall, maximum and minimum air temperatures, relative humidity and solar radiation and (2) soil hydraulic data to approximate unsaturated hydraulic conductivity and water retention. Each hydraulic polynomial function approximated using the Chebyshev polynomial and least square difference technique in tandem showed a fairly good fit of the given set of data. Vertical distribution of soil water as approximations to the Richards' equation subject to changing surface and phreatic boundaries was solved numerically during 53 days with a comparatively large time increment, and this pattern agreed well with field neutron scattering data, except for the surface 0.1 m slab.

• Economic and Environmental Geology
• /
• v.42 no.6
• /
• pp.561-574
• /
• 2009

Unsaturated/saturated groundwater flow and solute transport model, VSFRT2D(Variable Saturated Flow and Reactive Transport model) was developed considering effects of pumping, irrigation, and denitrification. VSFRT2D employed Richards equation as governing equation for groundwater flow and previously existing unsaturated models modified by including computational procedure of evapotranspiration at surface using Thornthwaite method when precipitation doesn't occur. Bioremediation processes based on monod kinetics are described using four nonlinear contaminant transport equations and three nonlinear microbes transport equations. The developed model was applied to field data in Hongsung area contaminated with nitrate. In order to identify the effect of precipitation, pumping, evapotranspiration, irrigation, fertilizer application, and various bioremediations on groundwater flow and contaminant transport, individual processes were separated and simulated. Then all results obtained from the individual processes are compared with each other. The simulation results show that bioremediation had a negligible effect on nitrate concentration change. However, pumping for irrigation, precipitation, and nitrogen fertilizer application showed profound influences on nitrate concentration change.

• Water for future
• /
• v.28 no.4
• /
• pp.215-222
• /
• 1995

The water retention function which has the hysteresis phenomena is required to analyze the Richards equation which is a governing equation of the unsaturated flow, and its hysteresis phenomena has influence upon the characteristics of the unsaturated flow. The accuracy of the published hysteresis models is compared by using experimental data of the water retention function. The apparatus to experiment the hysteresis phenomena on the soil is developed, and experimental data for the main wetting process and the main drying process of the water retention function are obtained. The parameters of the van Genuchten equation are calibrated by using experimentally obtained data. As a result of the comparison of the selected hysteresis models which simulate the main drying curve from the main wetting curve, the Model I-1(Mualem) overestimates and the Model II-1(Mualem) underestimates but the Model III-2(Park and Sonu) similarly estimates the experimental data of the main drying curve.

• Journal of Forest and Environmental Science
• /
• v.30 no.2
• /
• pp.201-207
• /
• 2014

This study was conducted to estimate the best-fit nonlinear height-DBH growth models for Larix kaempferi in Gangwon and North Gyeongsang province in South Korea. Exponential, Modified Logistic, Chapman-Richards, and Weibull function were used for estimating height-DBH models. To evaluate the selected models, $R^2$, RMSE, MD, MAD, and residual plots were performed in each model. Also, the coefficients and patterns in models of the previous studies were compared with those in this study. The result showed that Weibull equation was found to be the best-fit model with $R^2$=0.9837, RMSE=2.6133, MD=0.0089, and MAD=2.0896. All model parameters in our study had similar values to those in the previous models except for asymptotic parameter a. Our research result showed that Gangwon and North Gyeongsang province were superior to other provinces with regard to height growth for Larix kaempferi.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1137-1143
• /
• 1994

Through the numerical analysis of the capillary pressure-based Richards equation with and without the effect of the capillary hysteresis under the boundary condition having an antecedent rainfall. the moving tendency of the wetting front, the redistribution of the moisture content, infiltration rate, cumulative infiltration etc, were computed. The effect of the capillary hysteresis cannot be neglected in analyzing an unsaturated flow, and the more accurate results may be obtained by the consideration of the hysteresis effect. If the effect of the hysteresis cannot be considered, the analysis by the use of the main wetting curve may give more reliable result than that of the main drying curve.

• Water for future
• /
• v.29 no.2
• /
• pp.143-151
• /
• 1996

The dynamic responses in the subsurface outflow, surface saturation area, soil moisture storage are established by numerical experiments with Richards equation. In addition to this, the dynamical relationship between surface saturation area and subsurface outflow, and between surface saturation area and soil moisture storage are also determined by varying the hillslope shape, soil type, and boundary condition. The simulation results indicate that the dynamical relationships between surface saturation area and subsurface outflow, and between surface saturation area and soil moisture storage are approximated by the steadystate relationships. And the dynamic responses of subsurface outflow and surface saturation area are characterized by the double peaks although the rainfall pattern is asimple pulse input. As a result of numerical simulation, the physical mechanism for the occurrence of the double peaks is explained using the concept of variable source area.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.3
• /
• pp.125-131
• /
• 1997

In order to analyze the feasibility of the drain layer construction method, which is one of the beach protection methods, a hybrid model is constructed by combining the wave model and the seepage flow model. The used wave model is the analytic solution given by Shuto (1972). and the seepage flow model is used by Richards equation which governs the saturated-unsaturated flow in the porous media. It is concluded by the sensitivity analysis of the hybrid model that the most sensitive parameter in the flow field within the beach is the saturated hydraulic conductivity. The developed hybrid model will be efficiently used in the analysis of the parameter when the drain layers are constructed in the beach, if the field datas are obtained more.

• Journal of Korean Society of Forest Science
• /
• v.108 no.3
• /
• pp.364-371
• /
• 2019

This study was conducted to provide the basic information for a reasonable forest management plan and sustainable forest management by developing a dominant tree height growth model using diameter at breast height (DBH) and site index curves for Pinus densiflora and Chamaecyparis obtusa growing in Jeolla-do. The altitude, slope, orientation, soil type, height and DBH of a dominant tree, and the ages of trees were measured for 3055 Pinus densiflora trees (611 plots) and 3345 Chamaecyparis obtusa trees (699 plots), and these data were used to develop a customized afforestation map. In the dominant tree height growth model, the relationship to DBH was used in the Petterson, Michailow, and log equations. Also, a dominant tree height growth model in relationship to age used the Chapman-Richards, Schumacher, and Gompertz equations. The Petterson equation, which has a lower mean square error, was used to model dominant tree height growth in relationship to DBH. In the model of dominant tree height growth in relationship to age, three kinds of equations were considered to have little statistical difference. Therefore, the Chapman-Richards equation was chosen for modeling on the national level. Thirtyyears was used as the base age, which is an important factor for estimating the site index curves. In the results, a more varied range of site index family curves with 6-18 was developed for Pinus densiflora, and with 6-22 for Chamaecyparis obtusa. As the new site index curves indicated influences on growth of Pinus densiflora and Chamaecyparis obtusa, a reasonable forest management plan will be possible in the future for Jeolla-do.

• Wind and Structures
• /
• v.16 no.1
• /
• pp.117-136
• /
• 2013

A field study of wind-induced internal pressures in a flexible and porous industrial warehouse with a single dominant opening, of various sizes for a range of moderate wind speeds and directions, is reported in this paper. Comparatively weak resonance of internal pressure for oblique windward opening situations, and hardly discernible at other wind directions, is attributed to the inherent leakage and flexibility in the envelope of the building in addition to the moderate wind speeds encountered during the tests. The measured internal pressures agree well with the theoretical predictions obtained by numerically simulating the analytical model of internal pressure for a porous and flexible building with a dominant opening. Ratios of the RMS and peak internal to opening external pressures obtained in the study are presented in a non-dimensional format along with other published full scale measurements and compared with the non-dimensional design equation proposed in recent literature.