• Journal of Korean Society on Water Environment
• /
• v.23 no.1
• /
• pp.38-45
• /
• 2007

A one-dimensional kinematic wave model was used to calculate temporal and spatial changes of the highway runoff. Infiltration into pavement was considered using Darcy's law, as a function of flow depth and pavement hydraulic conductivity ($K_p$). The model equation was calculated using the method of characteristics (MOC), which provided stable solutions for the model equation. 22 storm events monitored in a highway runoff monitoring site in west Los Angeles in the U.S. were used for the model calculation and evaluation. Using three different values of $K_p$ ($5{\times}10^{-6}$, $10^{-5}$, and $2{\times}10^{-5}cm/sec$), total runoff volume and peak flow rate were calculated and then compared with the measured data for each storm event. According to the calculation results, $10^{-5}cm/sec$ was considered a site representative value of $K_p$. The study suggested a one-dimensional method to predict hydrodynamic behavior of highway runoff, which is required for the water quality prediction.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.4
• /
• pp.323-329
• /
• 2013

This study attempts to simulate the structure of a woody leaf netted venation system by using topology optimization techniques. Based on finite element method (FEM) analysis of an incompressible fluid, a topology optimal design is applied to those woody leaf netted venation models. To solve the transverse shear locking problem of a thin plate caused by the Mindlin-Reissner plate model where a leaf netted venation is assumed to be a thin plate, a P1-nonconforming element and selective reduced integration are employed. Topology optimal design is applied to multiple physical domains. Combined with the Darcy-Stokes flow problems and extended to the optimal design of fluid channels, the multiple physical models of the flow system are analyzed and venation patterns of leafs are simulated. The calculated optimal shapes are compared with the natural shapes of woody leaf venation patterns. This interdisciplinary approach may improve our understanding of the leaf venation system.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.1691-1695
• /
• 2010

본 연구에서 AFDA(Approximate Full Distribution Approach)를 사용하여 하수관의 불능확률을 정량적으로 산정할 수 있는 신뢰성 모형이 개발되었다. 여러 도시의 연 최대강우강도(Yearly Maximum Rainfall Intensity)를 이용하여 그 확률분포함수를 분석하였고 우수관(Storm sewer)의 불능확률 산정을 위한 신뢰성 모형에 적용하였다. 연 최대강우강도 자료의 분석결과 우리나라 중부지방의 여러 중소도시에 대한 연 최대강우강도의 확률분포함수는 Gumbel분포와 일치하는 것으로 나타났다. 신뢰성 모형은 불능확률의 신뢰함수를 구하기 위해 하중(Load)을 규정하는 식은 합리식이 사용되었고 용량(Capacity)를 규정하는 식은 Darcy-Weisbach공식과 Manning의 공식이 사용되었다. 이렇게 개발된 신뢰성 모형을 실제 우수관에 적용하여 불능확률을 산정하는 신뢰성 해석을 수행하였다. Y자형 우수관망에서 2개의 관으로 유입하는 각각의 유량이 그 관의 허용유량을 초과할 경우를 불능확률로 가정하였고, 나머지 관의 경우는 두 개의 관으로부터 유입하는 유량과 그 세 번째 관의 매설지역의 우수유입량의 합이 그 관의 허용유량을 초과할 경우를 불능상태(state of system failure)로 간주하여 불능확률을 정량적으로 산정하였다. Darcy-Weisbach공식과 Manning의 공식을 사용한 신뢰성 해석결과를 비교하였으며 우수관 직경의 변화에 따른 불능확률을 산정하였다. 특정한 수치(설계직경)이하일 경우 불능확률이 급격히 증가하는 것으로 나타났다. 따라서 실제 우수관의 유효직경이 설계직경에 항상 가깝도록 불순물을 제거하는 것이 최선의 관리 방법이며 불능확률을 줄이는 최선의 방법일 것이다. 본 연구에서 개발된 신뢰성 모형은 우수관의 운용, 관리, 감독은 물론 설계에 활용이 가능 할 것이다.

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

• Water Engineering Research
• /
• v.2 no.1
• /
• pp.49-61
• /
• 2001

A grid-based KIneMatic wave soil-water EROsion and deposition Model(KIMEROM) that predicts temporal variation and spatial distribution of sediment transport in a watershed was developed. This model uses ASCII-formatted map data supported from the regular gridded map of GRASS (U.S. Army CERL, 1993)-GIS(Geographic Information Systems), and generates the distributed results by ASCII-formatted map data. For hydrologic process, the kinematic wave equation and Darcy equation were used to simulated surface and subsurface flow, respectively (Kim, 1998; Kim et al., 1998). For soil erosion process, the physically-based soil erosion concept by Rose and Hairsine (1988) was used to simulate soil-water erosion and deposition. The model adopts single overland flowpath algorithm and simulates surface and subsurface water depth, and sediment concentration at each grid element for a given time increment. The model was tested to a 162.3 $\textrm{km}^2$ watershed located in the tideland reclaimed ares of South Korea. After the hydrologic calibration for two storm events in 1999, the results of sediment transport were presented for the same storm events. The results of temporal variation and spatial distribution of overland flow and sediment areas are shown using GRASS.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2001.05a
• /
• pp.25-34
• /
• 2001

A grid-based KIneMatic wave soil-water EROsion and deposition Model (KIMEROM) that predicts temporal variation and spatial distribution of sediment transport in a watershed was developed. This model uses ASCII-formatted map data supported from the regular gridded map of GRASS (U.S. Army CERL, 1993)-GIS (Geographic Information Systems), and generates the distributed results by ASCIIl-formatted map data. For hydrologic process, the kinematic wave equation and Darcy equation were used to simulate surface and subsurface flow, respectively (Kim, 1798; Kim et al., 1993). For soil erosion process, the physically-based soil erosion concept by Rose and Hairsine (1988) was used to simulate soil-water erosion and deposition. The model adopts sing1e overland flowpath algorithm and simulates surface and subsurface water depth, and sediment concentration at each grid element (or a given time increment. The model was tested to a 162.3 km$^2$ watershed located in the tideland reclaimed area of South Korea. After the hydrologic calibration for two storm events in 1999, the results of sediment transport were presented for the same storm events. The results of temporal variation and spatial distribution of overland flow and sediment areas are shown using GRASS.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.6
• /
• pp.883-892
• /
• 2002

A numerical study has been carried out to investigate the flow and heat transfer from an aluminum foam heat sink in a confined channel. A uniform heat flux is given at the bottom of the aluminum foam heat sink, which is horizontally placed on the heated surface. The channel walls are assumed to be adiabatic. Cold air is supplied from the top opening of the channel and exhausted to the channel outlet. Comprehensive numerical solutions are acquired to the governing Wavier-Stokes and energy equations, using the Brinkman-Forchheimer extended Darcy model and the local thermal non-equilibrium model f3r the region of porous media. Details of flow and thermal fields are examined over wide ranges of the principal parameters; i.e., the Reynolds number Re, the height of heat sink h/H, porosity $\varepsilon$and pore diameter ratio $R_{H}$.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.06a
• /
• pp.5-15
• /
• 1997

It is the objective of this study that by conducting the serni-solid extrusion using A12024, the effect of various process variables on the quality of extruded product and extrusion force is understood. The results of experiment are compared with those of finite element simulation in order to verify the effectiveness of the developed FE-simulation code. In order to simulate densification in the deformation of serni-solid material, the semi-solid material is assumed to be composed of solid region as porous skeleton following compressible visco-plastic model and liquid region following Darcy's equation for the liquid flow saturated in the interstitial space. Then the flow and deformation of the semi-solid alloy are analyzed by coupling the deformation of the porous skeleton and the flow of the eutectic liquid. It is assumed that initial solid fraction is homogeneous. Yield and plastic potential function presented by Kuhn and constitutive model developed by Gunasekera are used for solid skeleton.

• Journal of Wetlands Research
• /
• v.14 no.1
• /
• pp.75-87
• /
• 2012

Hyporheic zone is a region beneath and alongside a stream, river, or lake bed, where there is mixing of shallow groundwater and surfacewater. Hyporheic exchange controls a variety of physical, biogeochemical and thermal processes, and provides unique ecotones in a aquatic ecosystem. Field and experimental observations, and modeling studies indicate that hyporheic exchange is mainly in response to pressure gradients driven by the geomorphological features of stream beds. In the reach scale of a stream, pool-riffle structures dominate the exchange patterns. Flow over a pool-riffle sequence develops recirculation zones and stagnation points, and this flow structures make irregular pressure gradient which is driving force of the hyporheic exchange. In this study, 3 D hydro-dynamic model solves the Reynolds-averaged Navier-Stokes equations for the surface water and Darcy's Law and the continuity equation for ground water. The two sets of equations are coupled via the pressure distribution along the interface. Simulation results show that recirculation zones and stagnation points in the pool-riffle structures dominantly control the upwelling and downwelling patterns. With decrease of recirculation zones, length of donwelling zone formed in front of riffles is reduced and position of maximum downwelling point moves downward. The numerical simulation could successfully predict the behavior of hyporheic exchange and contribute the field study, river management and restoration.