• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.56-59
• /
• 2010

포화된 다공질 매체의 수치해석은 일반적으로 혼합유한요소방법(Mixed Finite Element Method)이 쓰인다. 이 혼합유한요소 방법은 고체변형과 유체의 이동을 동시에 고려하게 되는데 고체의 변형이 거의 없이 유체만 이동할 경우나 고체와 유체의 변형이 없이 간극수압만 존재할 경우에는 요소잠김현상(Element Locking)이 발생하여 혼합유한요소방법으로 해석하기에는 수치적으로 불안정해 진다. 본 논문에서는 이러한 수치적 불안정성을 해결한 스태거드 방법(Park and Tak 2010)을 소개하고 수치적 효율성을 위해 MPI(Message-Passing Interface) 라이브러리를 이용한 병렬해석 기법이 적용된다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.369-378
• /
• 2010

The mixed finite element analysis is the most widely used method for saturated porous media. Generally, in this method, direct method and iterative method are proposed to obtain unknown variable, however, the iterative method is recommended because the method provide numerical stability and accuracy under the material properties for solid and fluid are different. In this paper, we introduce staggered method which has strong numerical stability, and FETI(Finite Element Tearing and Interconnecting) which is one of decomposition methods are applied into the method in order to obtain numerical efficiency. In which, Lagrange Multipliers and conjugated gradient method to solve decomposed domain are proposed, and then, the proposed method is verified numerical efficiency by point to point MPI(Message Passing Interface) library.

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

본 연구에서는 수리모형실험을 통해 해안지역 대수층에서 해수침투가 일어나지 않는 지하수의 최적 양수량을 산정하였다. 해안 대수층 및 양수정을 구현하기 위해 가로 140 cm, 세로 70 cm, 두께 10 cm 의 직사각형 Sandtank 모형과 5 mm 원형관을 이용하여 수행하였다. 해수의 염도, 해안가의 지표경사 등을 기준으로 설정하였으며, 다공질 매체의 입자크기에 따른 수리전도도 2가지를 경계조건으로 하여 각각의 양수정 위치에 따라 최대 담수량을 모의하였으며, 이 때, 하나의 과잉 담수량에 따른 최적 염수량을 산정하여 염수침투 제어 효과를 검증하였다. 다공질매체의 수리전도도가 높을수록 최대 담수량이 증가되며 이에 따라 염수 침투 제어를 위한 최적염수량이 증가하였다. 또한 양수정의 위치가 염수조에 가까울수록 최대담수량 및 최적염수량이 감소되었다. 이러한 모의를 통하여 설정된 각각의 양수정 위치에서 최적의 염수량으로 염수침투를 제어함으로써 지속시간 이용가능한 최대담수량을 나타내었다.

• Journal of the KSME
• /
• v.16 no.1
• /
• pp.49-53
• /
• 1976

Groundwater flow in a two-layer aquifer is studied theoretically and experimentally. Nonlinear differential equations are independently linearized for the upper and the lower layer, and approximate solutions are obtained. Fine sands and coarse sands are used in experiments. For steady flow, theoretical solutions agrees excellently with experimental results For unsteady flow, however, agreement is not so good, because the experimental conditions do not match with theoretical assumptions.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 1997.11a
• /
• pp.447-450
• /
• 1997

Temperature distribution measurements in the mockup apparatus of reactor vessel were performed to determine the effective thermal conductivity of Al powder porous media where stainless steel tubes were installed with different geometry. The temperature distributions at four separated sections with different arrangements of porous media have different slopes according to the geometrical configuration. From the measured temperature distribution, effective thermal conductivity have been derived using the least square fitting method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1997.05a
• /
• pp.475-479
• /
• 1997

다공질매체에서의 포화-비포화 흐름 거동을 파악하기 위하여는 시간에 따른 함수량의 변화과정을 정확하고 빠르게 측정하여야 한다. 본 연구는 실험실에서 함수량을 측정하는 방법의 하나로서 TDR(Time Domain Reflectometer)을 사용하는 실험 방법에 관한 연구이다. TDR 이란 전기신호의 전도특성 이용하여 토양내 설치된 탐침(probe)의 전기신호 반향시간을 측정하는 기기로서 이 반향시간과 토양의 유전율상수(dielectric constant)의 관계를 이용하여 함수량을 추정할 수 있다. 본 연구에서는 TDR의 원리설명과 이를 이용한 함수량측정방법을 제시하고 있다.

• The Journal of Engineering Geology
• /
• v.3 no.2
• /
• pp.203-214
• /
• 1993

Underground Injection Control regulations for Class I injection wells require that the vertical dufusion be considered as a mechanism for transport of contaninants in evaluating containment Due to the low permeability in the confining aquitard, the movement of contaminants over the long term is controlled by the molecular diffusion. The movement can be predicted, using the Fick' S second law of diffusion. The diffusion coefficient in Fick' s law has been determined experimentally in this study. Instantaneous injection of $^{51}Cr$ was used to trace the distribution of Cr(VI) in soil plugs and an analytical solution was applied to calculate the diffusion coefficients. This study shows the effect of environmental factors, such as temperature, chloride concentration, applied amount of $^{51}Cr$, and bulk density of injection formations on diffusion of Cr(VI)

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.20-25
• /
• 2008

The finite element analysis for porous media is severe job because constituents have different physical peoperties, and element's continuity and stability should be considered. Thus, we propose the new mixed finite element method in order to overcome the problems. In this method, multi time step, remeshing step, and sub iteration step are introduced. The multi time step and remeshing step make it possible to satisfy a stability and an accuracy during sub iteration in which global time is determined. Finally, the proposed method is compared with the ABAQUS(2007) software and exact solution(Schiffman 1967) through two dimensional consolidation model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.235-242
• /
• 2003

The solids and the fluids in porous media have a relative velocity to each other. Due to physically and chemically different material properties and their relative velocity, the behavior of saturated porous media is extremely complicated. Thus, in order to describe and clarify the deformation behavior of saturated porous media, constitutive models for deformation of porous media coupling several effects such as flow of the fluids or thermodynanical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian elements, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of the solids and the fluids. In this work, governing equations of porous media based on ALE description are obtained from governing equations in frame of updated Lagrangian description. Then, weak forms of these equations are derived using arbitrary weighting functions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.381-388
• /
• 2006

A porous medium is composed of solids, fluids, and gas which have different physical and chemical properties. In addition, these constituents have a relative velocity between each other. So far, in order to analyze porous media using finite element method, Lagrangian or Eulerian method has been used. However, the numerical analyses for porous media have a defect that the methods do not describe the movements of constituents. In this paper, numerical analysis for unsaturated porous media was performed in frame of ALE method which has advantages of Lagrangian and Eulerian. Namely, the Lagrangian description was used in solid phase, and the Eulerian description was used in fluid or gas phase in a porous medium Then the relationship between each other was controlled by the convective term in ALE method. Finally, the numerical results of ALE were compared with tile results of Lagrangian analysis.