• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1B
• /
• pp.21-27
• /
• 2012

The HLLL scheme, proposed by T. Linde, determines all the wave speeds from the initial states because the middle wave is evaluated by the introduction of a generalized entropy function. The scheme is considered a genuine successor to the original HLL scheme because it is completely separated form the Roe's linearization scheme unlike the HLLE scheme and does not rely on the exact solution unlike the HLLC scheme. In this study, a numerical model was configured by the HLLL scheme with the total energy as a generalized entropy function to solve governing equations, which are the one-dimensional shallow water equations without source terms and with an additional conserved variable relating a concentration. Despite the limitations of the first order solutions, results to three cases with the exact solutions were generally accurate. The HLLL scheme appeared to be superior in comparison with the other HLL-type schemes. In particular, the scheme gave fairly accurate results in capturing the front of wetting and drying. However, it revealed shortcomings of more time-consuming calculations compared to the other schemes.

• Tunnel and Underground Space
• /
• v.17 no.6
• /
• pp.464-474
• /
• 2007

An enhancement in the performance and safety of a high-level waste repository requires a validation of its engineered barrier. An engineering-scale test (named "KENTEX") has been conducted to investigate the thermal-hydro-mechanical behaviors in the engineered barrier of the Korean reference disposal system The validation test started on May 31, 2005 and is still under operation. The experimental data obtained allowed a preliminary and qualitative interpretation of the thermal-hydro-mechanical behaviors in the bentonite blocks. The temperature was higher as it became closer to the heater, while it became lower as it was farther away from the heater. The water content had a higher value in the part close to the hydration surface than that in the heater part. The relative humidity data suggested that a hydration of the bentonite blocks might occur by different drying-wetting processes, depending on their position. The total pressure was continuously increased by the evolution of the saturation front in the bentonite blocks and thereby the swelling pressure. Near the heater region, there was also a significant contribution of the thermal expansion of bentonite and the vapor pressure in the pores of the bentonite blocks.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• International Journal of Automotive Technology
• /
• v.5 no.4
• /
• pp.221-238
• /
• 2004

Wall interactions of direct injection spray were investigated using laser-sheet imaging, shadowgraphy, wetted footprint and phase Doppler interferometry techniques. A narrow-cone high-pressure swirl injector is used to inject iso-octane fuel onto a plate, which has three different impact angles inside a pressurized chamber. Heated air and plate conditions were compared with unheated cases. Injection interval was also varied in the heated case to compare dry- and wet- wall impingement behaviors. High-speed macroscopic Mie-scattering images showed that presence of wall and air temperature has only minor effect on the bulk spray structure and penetration speed for the narrow-cone injector tested. The overall bulk motions of the spray plume and its spatial position at a given time are basically unaffected until a few millimeters before impacting the wall. The surface properties of the impact surface, such as the temperature, the presence of a preexisting liquid film also have a small effect on the amount of wetting or the wetted footprint; however, they have strong influence on what occurs just after impact or after a film is formed. The shadowgraph in particular shows that the plate temperature has a significant effect on vapor phase propagation. Generally, 10-20% faster horizontal vapor phase propagation is observed along the wall at elevated temperature condition. For impingement onto a preexisting film, more splash and evaporation were also observed. Contrary to some preconceptions, there is no significant splashing and droplet rebounding from surfaces that are interposed in the path of the DI gasoline spray, especially for the oblique impact angle cases. There also appears to be a dense spray front consists of large sac spray droplets in the oblique impact angle cases. The bulk of the spray is not impacted on the surface, but rather is deflected by it The microscopic details as depicted by phase Doppler measurements show that the outcome of the droplet impaction events can be significantly influenced. Only droplets at the spray front have high enough Weber numbers for wall impact to wet, splash or rebound. Using the sign of vertical velocity, the time-resolved downward droplets and upward droplets are compared. The Weber number of upward moving droplets, which seldom exceeds unity, also decreases as the impact angle decreases, as the droplets tend to impact less and move along the wall in the deflected spray plume.

• The Journal of Engineering Geology
• /
• v.27 no.1
• /
• pp.21-29
• /
• 2017

Variations in matric suction were compared between field measurements and numerical analysis for a tailings dump slope. We performed an analysis of slope stability using precipitation records measured in the field, selecting the period (72 hours) of highest rainfall intensity during the field monitoring, for which the matric suction and rainfall records measured in the field were analyzed. In addition, we applied the precipitation records of this period to the numerical analysis. SEEP/W and SLOPE/W were used to analyze the seepage flow in the slope due to rainfall and to investigate the slope stability considering the wetting front depth from the ground surface, respectively. The seepage analysis, using SEEP/W, was carried out on the rainfall data obtained in the field. Comparisons between the field monitoring data and simulation results for matric suction show some quantitative difference but similar patterns of temporal variation. According to the results of slope stability analysis using SLOPE/W and the results of seepage analysis, the safety factor of the slope showed a sudden increase at the point of rapid increase in rainfall intensity. It then recovered because subsequent rainfall was scarce. Therefore, the stability of the tailings dump slope can be reasonably estimated if seepage and slope stability analyses, based on precipitation records, have been carried out.

• Korean Journal of Soil Science and Fertilizer
• /
• v.47 no.6
• /
• pp.406-411
• /
• 2014

Serious problems in reclaimed land agriculture are high soil salinity and poor vertical drainage, so desalinization in these soils is very difficult. Also, although desalinization is accomplished in reclaimed top soils, before long, soils are resalinized according to capillary rise of salts from the subsurface soils. To resolve these problems, multi-layered soil columns with subsurface layer of macroporous medium utilizing coal bottom ash (CBA) were constructed and the effects of blocked resalinization of these soils were investigated. In this experiment soil samples were collected from Munpo series (coarse-loamy, nonacid, mixed, mesic, typic Fluvaquents). The soil texture was silt loam and the EC was $33.9dS\;m^{-1}$. As for groundwater seawater was used and groundwater level of 1 cm from the bottom was maintained. The overall rate of capillary rise was $2.38cm\;hr^{-1}$ in soil 60 cm column, $0.25cm\;hr^{-1}$ in topsoil (30 cm) + CBA (5 cm) + subsurface soil (10 cm) column and $0.08cm\;hr^{-1}$ in topsoil (30 cm) + CBA (10 cm) + subsurface soil (10 cm) column. In multi-layered soil columns with CBA 20, 30 cm layer, wetting front due to capillary rise could not be seen in top soil layer. After 70 days capillary rise experiment water soluble Na+ accumulated in top soil of soil columns with CBA 20, 30 cm was diminished by 92.8, 96.5% respectively in comparison with Na+ accumulated in top soil of soil 60 cm column because CBA layer cut off capillary rise of salts from the subsurface soil. From these results we could conclude that the macroporous layer utilizing CBA placed at subsurface layer cut off capillary rise of solutes from subsurface soil, resulting in lowered level of salinity in top soil and this method can be more effective in newly reclaimed saline soil.

• The Korean Journal of Rheology
• /
• v.10 no.4
• /
• pp.185-194
• /
• 1998

Flow-induced voids during resin impregnation and poor fiber wetting give serious effects on the mechanical properties of composites in resin transfer molding process. In order to better understand the characteristics of resin flow and to investigate the mechanism of void formation, flow visualization experiment for the resin impregnation was carried out on plain weaving glass fiber mats using silicon oils with various viscosity values. The permeability and the capillary pressure for the fiber mats of different porosities were obtained by measuring the penetration length of the resin with time and with various injection pressure. At low porosity and low operating pressure, the capillary pressure played a significant role in impregnation process. Video-assisted microscopy was used in taking the magnified photograph of the flow front of the resin to investigate the effect of the capillary pressure on the void formation. The results showed that the voids were formed easily when the capillary pressure was relatively high. No voids were detected above the critical capillary number of 2.75$\times$$10^{-3}, and below the critical number the void content increased exponentially with decrease of the capillary number. The content of void formed was independent of the viscosity of the resin. For a given capillary number, the void content reduced with the lower porosity of the fiber mat.

• The Journal of Engineering Geology
• /
• v.28 no.3
• /
• pp.489-499
• /
• 2018

A flume experiments was used to study the characteristics of the surface displacements and volumetric water contents (VWC) during torrential rain. The surface displacement and VWC of the granite weathered soil were measured for rainfall intensity (100, 200 mm/hr) and initial ground condition (VWC 7, 14, 26%). The test processes were also recorded by video cameras. According to the test results, The shallow failure is classified into three types: retrogressive failure, progressive failure and defined failure. In the case of retrogressive failure and progressive failure, relatively large damage could occur due to the feature that soil is deposited to the bottom of the slope. the shallow failure occurred when the VWC reached a certain value regardless of the initial soil condition. It was found that the shallow failure can be predicted through the increase patton of the VWC under the condition of the ground dry condition (VWC 7%) and the natural condition (VWC 14%). For high rainfall intensity, progressive failure predominated, and rainfall intensity above a certain level did not affect wetting front transition.

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

• The Journal of Engineering Geology
• /
• v.25 no.3
• /
• pp.377-385
• /
• 2015

We constructed a model test apparatus to evaluate the dependence of the saturation velocity (V_s) in soils on rainfall intensity (I_R). The apparatus comprises a soil box, a rainfall simulator, and measuring sensors. The model grounds (60 cm × 50 cm × 15 cm) were formed by Joomunjin standard sand with a relative density of 75%. The rainfall simulator can control the rainfall intensity to reenact the actual rainfall in a soil box. Time Domain Reflectometer (TDR) sensors and tensiometers were installed in the soils to measure changes in the volumetric water content and matric suction due to rainfall infiltration. During the tests, the soil saturation was determined by raising the groundwater table, which was formed at the bottom of the soil box. [Please check that the correct meaning has been maintained.] The wetting front did not form at the ground surface during rainfall because the soil particles were uniform and the coefficient of permeability was relatively high. Our results show that the suction stress of the soils decreased with increasing volumetric water content, and this effect was most pronounced for volumetric water contents of 20%-30%. Based on a regression analysis of the relationship between rainfall intensity and the average saturation velocity, we suggest the following equation for estimating the saturation velocity in soils: V_savg (cm/sec) = 0.068I_R (mm/hr).