• Economic and Environmental Geology
• /
• v.29 no.3
• /
• pp.325-333
• /
• 1996

numerical model is developed to estimate gas flow in the landfill site. Darcy's law, the mass conservation law, and the ideal gas state equation are combined to compose the governing equation for the steady-state and transient-state gas flows. The finite element method (FEM) is used as the numerical solution scheme. Two-dimensional radial symmetric triangular ring element is used to discretize the simulation domain. The steady state model developed in this study is compared with AIRFLOW that is a commercial model developed by Hydrologic Inc. Mass balance test is performed on the transient gas flow simulation. The developed model is applied to analyze the gas extraction experiment performed by Daewoo Institute of Construction Technology at the Nanjido landfill in 1993. The developed model was registered at Korea Computer Program Protection Foundation.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.55 no.6
• /
• pp.145-154
• /
• 2013

Recently, damage occurrence by wind erosion has been increasing in society. In times past, such problems only took place in desert area ; however, in recent years, the wind erosion problem is spreading out to agricultural land. Wind erosion in agricultural land can cause loss of loam soils, the disturbance of the photosynthesis of the crop fields and serious economic losses. To overcome the mentioned problems, installation of windbreak fence can be recommended which function as disturbing strong wind and wind erosion. However, there is still no proper guideline to install the windbreak fence and the installation used to rely on the intuition of the workers due to the lack of related studies. Therefore, this study measured the aerodynamic resistance of screens of the windbreak fence using the apparatus for testing screens. The apparatus for testing screens was designed to measure pressure loss around the screen. Measured pressure loss by wall friction compensated for pressure loss to calculate the aerodynamic resistance of screens. The result of pressure loss by regression analysis derived the aerodynamic coefficient of Darcy-Forchheimer equation and power law equation. The aerodynamic resistance was constant regardless of the overlapped shape when the screen was overlapped into several layers. Increasing the number of layers of the screen, internal resistance increased significantly more, and pressure loss caused by the screen also increased linearly when the wind speed was certain conditions, but permeability had no tendency. In the future, the results of this study will be applied to the computational fluid dynamics simulation. The simulation models will be also validated in advance by wind tunnel experiments. It will provide standard of a design for constructing windbreak fence.

• Biomaterials and Biomechanics in Bioengineering
• /
• v.2 no.2
• /
• pp.111-125
• /
• 2015

For engineers, generating a mesh in porous media (PMs) sometimes represents a smaller computational load than generating realistic stent geometries with computer fluid dynamics (CFD). For this reason, PMs have recently become attractive to mimic flow-diverter stents (FDs), which are used to treat intracranial aneurysms. PMs function by introducing a hydraulic resistance using Darcy's law; therefore, the pressure drop may be computed by test sections parallel and perpendicular to the main flow direction. However, in previous studies, the pressure drop parallel to the flow may have depended on the width of the gap between the stent and the wall of the test section. Furthermore, the influence of parameters such as the test section geometry and the distance over which the pressure drops was not clear. Given these problems, computing the pressure drop parallel to the flow becomes extremely difficult. The aim of the present study is to resolve this lack of information for stent modeling using PM and to compute the pressure drop using several methods to estimate the influence of the relevant parameters. To determine the pressure drop as a function of distance, an FD was placed parallel and perpendicular to the flow in test sections with rectangular geometries. The inclined angle method was employed to extrapolate the flow patterns in the parallel direction. A similar approach was applied with a cylindrical geometry to estimate loss due to pipe friction. Additionally, the pressure drops were computed by using CFD. To determine if the balance of pressure drops (parallel vs perpendicular) affects flow patterns, we calculated the flow patterns for an ideal aneurysm using PMs with various ratios of parallel pressure drop to perpendicular pressure drop. The results show that pressure drop in the parallel direction depends on test section. The PM thickness and the ratio of parallel permeability to perpendicular permeability affect the flow pattern in an ideal aneurysm. Based on the permeability ratio and the flow patterns, the pressure drop in the parallel direction can be determined.

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

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

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.4
• /
• pp.265-273
• /
• 2002

In this paper, wave absorbing characteristics of a horizontal submerged punching plate are investigated throughout the calculation and the experiment. The punching plate with the array of circular holes can force the flow to separate and to form eddies of high vorticity and cause significant energy loss. As an analytic tool, the linear water wave theory and the eigenfunction expansion method is applied. Darcy's law that the normal velocity of the fluid passing through the punching plate is linearly proportional to the pressure difference between two sides of the punching plate is assumed. The proportional constant called the porous coefficient is deeply dependent to the porosity. To obtain the relationship between the porosity and the porous coefficient the systematic model test for the punching plates with 6 different porosities is conducted at 2-dimensional wave tank. It is found that the porous coefficient is linearly proportional to the porosity(b=57.63P-0.9717). It is also noted that the optimal porosity value is near P=0.1 and the optimal range of submergence depth is $d/h\\leq0.2$ within entire frequency range.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.21 no.1
• /
• pp.30-38
• /
• 2009

In recent years, there's been strong demand for coastal structures that have a permeability that serves water affinity and disaster prevention from wave attack. The aim of this study is to examine the wave transformation, including wave run-up that propagates over the coastal structures with a steep slope. A numerical model based on the nonlinear shallow water equation, together with the unsteady nonlinear Darcy law for fluid motion in permeable underlayer and laboratory measurements was carried out in terms of the free surface elevations and fluid particle velocities for the cases of regular and irregular waves over 1:5 impermeable and permeable slopes. The numerical results were used to evaluate the application and limitations of the PBREAK numerical model. The numerical model could predict the cross-shore variation of the wave profile reasonably, but showed less accurate results in the breaking zone that the mass and momentum influx is exchanged the most. Except near the wave crest, the computed depth averaged velocities could represent the measured profile below the trough level fairly well.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.1
• /
• pp.63-69
• /
• 2011

In this paper, the wave forces on a cylindrical aquaculture fish cage, which consists of the porous mesh with the uniform porosity, are analyzed using matched eigenfunction expansion method. The boundary condition on the porous net is derived based on the Darcy's law, which implies that the velocity of the fluid passing through the net is linearly proportional to the pressure difference between two sides of the net. The wave forces and wave responses are investigated by changing the porous parameter of porous net as well as the submerged position (floating type, bottom-mounted type) of an aquaculture fish cage. It is found that the wave forces on a bottom-mounted type are largely decreased compared with that on a floating type. Also, the porosity of the netting structure plays an important role in reducing the wave forces and the wave elevation in the vicinity of an aquaculture fish cage.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.2B
• /
• pp.139-147
• /
• 2012

This study is derived relationships of the resistance coefficients of Darcy-Weisbach and Manning for flow resistance and the dimensionless velocity using many field measurements for 1,875 rivers consist of sand 179, gravel 992, cobble 651 and boulder 53 channels in natural rivers, respectively. The relationships of power law forms are developed as a function of flow discharge, friction slope, and relative submergence by the regression and the semi-empirical method. The measurements distribution of Manning resistance coefficients by the Box-Whisker Plots show the values which ranges from 0.004~0.151 for sand, 0.008~0.250 for gravel, 0.015~0.327 for cobble, 0.023~0.444 for boulder in natural rivers, respectively. Relationships of these semi-empirical and resistance coefficients will be useful to give information in hydraulic engineering.