• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.1B
• /
• pp.1-6
• /
• 2006

Amount of coastal ground water discharge(CGD) of surficial aquifer via coastlines of Nakdong River watershed, Seomjin River watershed, Youngsan River watershed and Keum River watershed is estimated. Compared to other major hydrological components, such as evapo-transpiration and river discharge, CGD is not so large in the amount. However, it is important since coastal ground water can be developed relatively free of environmental impact on downstream area and since most of coastal areas currently suffers water shortage. Regional groundwater investigation data and assessment based on Darcy's law are used for estimating coastal groundwater discharge. In this work the amount of CGD across the coastlines of the four rivers is estimated as 1.8 billion cubic meter per year and that is about 2.3 percentage of total amount of annual precipitation. Nakdong River watershed is most appropriate region in view of developing groundwater.

• Journal of the Korea Society of Computer and Information
• /
• v.29 no.6
• /
• pp.23-29
• /
• 2024

In this paper, we propose a fast GPU-based method for representing pore flow, absorption, emission, and diffusion effects represented by cloth-liquid interactions using smoothed particle hydrodynamics (SPH), a particle-based fluid solver: 1) a unified framework for GPU-based representation of various physical effects represented by cloth-liquid interactions; 2) a method for efficiently calculating the saturation of a node based on SPH and transferring it to the surrounding porous particles; 3) a method for improving the stability based on Darcy's law to reliably calculate the direction of fluid absorption and release; 4) a method for controlling the amount of fluid absorbed by the porous particles according to the direction of flow; and finally, 5) a method for releasing the SPH particles without exceeding their maximum mass. The main advantage of the proposed method is that all computations are computed and run on the GPU, allowing us to quickly model porous materials, porous flows, absorption, reflection, diffusion, etc. represented by the interaction of cloth and fluid.

• Journal of Korea Water Resources Association
• /
• v.45 no.1
• /
• pp.39-52
• /
• 2012

Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed. In general, these interactions are considered to be one of the most difficult areas of the discipline, particularly for the modeler who intends simulate the dynamic relations between these two major domains of the hydrological cycle. In essence, one major complexity is the spatial and temporal variations in the dynamically interacting system behavior. The proper simulation of these variations requires the need for providing an appropriate coupling mechanism between the surface and subsurface components of the system. In this study, an approach for modelling surface-subsurface flow and transport in a fully intergrated way is presented. The model uses the 2-dimensional diffusion wave equation for sheet surface water flow, and the Boussinesq equation with the Darcy's law and Dupuit-Forchheimer's assumption for variably saturated subsurface water flow. The coupled system of equations governing surface and subsurface flows is discretized using the finite volume method with central differencing in space and the Crank-Nicolson method in time. The interactions between surface and subsurface flows are considered mass balance based on the continuity conditions of pressure head and exchange flux. The major module consists of four sub-module (SUBFA, SFA, IA and NS module) is developed.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.6
• /
• pp.27-44
• /
• 2015

An analytical solution of dynamic responses for seabed in finite and infinite thicknesses including shallow has been developed under flow and partial standing wave with arbitrary reflection ration coexisting field at a constant water depth condition. In the analytical solution, a field was simply transited to a coexisting field of progressive wave and flow when reflection ratio was 0 and to a coexisting field of fully standing wave and flow when reflection ratio was 1. Based on the Biot's consolidation theory, the seabed was assumed as a porous elastic media with the assumptions that pore fluid is compressible and Darcy law governs the flow. The developed analytical solution was compared with the existing results and was verified. Using the analytical solution the deformation, pore pressure, effective and shear stresses were examined under various given values of reflection ratio, flow velocity, incident wave's period and seabed thickness. From this study, it was confirmed that the dynamic response of seabed was quite different depending on consideration of flow, which causes changing period and length of incident and reflection waves. It was also confirmed that dynamic response significantly depends on the magnitude of reflection ratio.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.2
• /
• pp.118-134
• /
• 2015

An analytical solution of dynamic responses for seabed in shallow, finite and infinite thicknesses has been developed under flow and standing wave coexisting field at a constant water depth condition. To do this, based on the Biot's consolidation theory, the seabed is assumed as a porous elastic media with the assumptions that pore fluid is compressible and Darcy law governs the flow. The developed analytical solution is compared with the previous results and is verified. Using the analytical solution the deformation, pore pressure, effective and shear stresses of seabed are examined under various given values of flow velocity, incident wave period and seabed thickness. From this study, it is confirmed that the seabed response is quite different depending on consideration of flow, which causes changing period and length of incident and reflection waves.

• Journal of Korea Water Resources Association
• /
• v.37 no.11
• /
• pp.939-947
• /
• 2004

Amount of ground water discharging through coastline(CGD) is small compared to the amount of river discharges. However, CGD is important for gound water development in coastal areas and for contamination of coastal waters by nutrients transported by ground water. In the previous study, a method was proposed to estimate CGD from simple data such as ground water levels at observation wells and hydraulic conductivities. In this study, the method is enhanced by taking into account complexity of coastlines. A new method is proposed to estimate average hydraulic gradient, length of coastline, and average thickness of freshwater The new method is verified against numerical solutions. It is shown that the method is able to estimate CGD from complex coastlines more accurately than the previous method.

• Transactions of the Korean hydrogen and new energy society
• /
• v.30 no.6
• /
• pp.629-634
• /
• 2019

Dicyclopentadiene is a low viscosity resin which forms a poly-dicyclopentadiene rapidly through ring opening metathesis polymerization (ROMP). This poly-dicyclopentadiene has outstanding properties of low-temperature, water and impact resistances. Due to these advantages, military and offshore structures try to apply the DCPD composites by using liquid composite molding process. In this study, 14%, 38% volume fraction fiber glass strand mat reinforced p-DCPD composites processed by structural reaction injection molding (S-RIM) which has resin-catalsyt mixing head and glass fiber preform in the mold. Additionally, S-RIM numerical analysis was conducted to predict the process time depending on fiber volume fraction and mold temperature. The process time is shorter when it has the lower fiber volume fraction or the higher mold temperature. At higher mold temperature, it is necessary to set the maximum mold temperature considering the resin curing time.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.979-990
• /
• 1989

Natural convection heat transfer has been investigated numerically in the vertical annulus filled withsaturated porous material for the aspect ratio less than unity. The inner wall of the annulus is exposed to constant heat flux condition and the outer wall is cooled to keep isothermal condition. The upper and the lower horizontal wall are assumed to be insulated. Under conditions ranging 50 .leq. Ra .leq. 10000, 1 .leq. RD .leq. 12, the characteristics of flow and heat transfer have been investigated. The results show that average Nusselt numbers increase when the radius ratio increases and the multicellular flows are not detected under the present conditions. Isothermal lines are plotted within the porous media. Temperatures of the inner wall with constant heat flux conditions and the local heat flux rate of the cooled outer wall with constant temperature are also obtained.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.7
• /
• pp.75-89
• /
• 2013

This paper describes a GIS-based geohydrologic methodology, called YSGWF (YonSei GroundWater Flow) for predicting the rainfall infiltration-groundwater flow of slopes. This physical-based model was developed by the combination of modified Green-Ampt model that considers the unsaturated soil parameters and GIS-based raster model using Darcy's law that reflects the groundwater flow. In the model, raster data are used to simulate the three dimensional inclination of bedrock surface as actual topographic data, and the groundwater flow is governed by the slope. Also, soil profile is ideally subdivided into three zones, i.e., the wetting band zone, partially saturated zone, and fully saturated zone. In the wetting band and partially saturated zones the vertical infiltration of water (rainfall) from surface into ground is modeled. When the infiltrated water recharges into the fully saturated zone, the horizontal flow of groundwater is introduced. A comparison between the numerical calculation and real landslide data shows a reasonable agreement, which indicate that the model can be used to simulate real rainfall infiltration-groundwater flow.