• Tunnel and Underground Space
• /
• v.14 no.3
• /
• pp.203-214
• /
• 2004

For total system performance of a potential radwaste repository, a hypothetical site is assumed with feasible boundary conditions. Assuming a coastal repository, the sensitivity of the depth and the location of a repository along with the distance to a joint on groundwater transport pathways is studied. Results from Connectflow analysis could be used as input of the MASCOT-t the probabilistic safety assessment code.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• Tunnel and Underground Space
• /
• v.14 no.3
• /
• pp.229-240
• /
• 2004

Coastal aquifers may serve as major sources fur freshwater. In many coastal aquifers, intrusion of seawater has become one of the major constraints imposed on groundwater utilization. The management of groundwater in coastal acquifers means making decision as to the pumping rate and the spatial distribution of wells. Several numerical techniques for flow and solute transport simulation can provide the means to achieve this goal. As a basic study to predict the intrusion of seawater in coastal phreatic aquifers, the coupled flow and solute transport analysis was conducted by use of the 3-D finite element code, SWICHA. In order to understand how the location and the shape of freshwater-seawater transition zone were affected by the boundary conditions and hydrogeologic variables, parametric study was carried out.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.47-54
• /
• 1993

For the numerical analysis of free oscillation characteristics in a harbor with general boundary and bottom topography, finite element method is applied. The governing Helmholtz equation is transformed into a generalized matrix eigenvalue problem using the standard finite element procedure. A computer code is developed for the numerical evaluation of natural frequencies and free oscillation modes. In the eigensolution process, a shifting strategy is devised for the treatment of numerical singularity. Scaling of coefficient matrix is also found to be effective for the alleviation of numerical ill-conditioning. For the test problems, firstly, analytical and numerical solutions are compared and validity of the code is obtained. Hence the method is successfully applicable for the real-world problems with general geometric boundaries and bottom topography.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2395-2406
• /
• 1992

Laminar natural convection heat transfer from a hot body in a square enclosure has been studied for various center positions of a hot body at Grashof number Gr=1.5$\times$10／sup 5／, Prandtl number Pr=0.71 and dimensionless thermal conductivity K／sub s／／K／sub f／= 14710. In case of vertical cold walls, the natural convection at the dimensionless center position of a hot body, X／sub c／Y／sub c／=0.2, 0.5 shows the most strong and at X／sub c／, Y／sub c／=0.5, 0.8 the most weak. In case of horizontal cold walls, the natural convection at the dimensionless center position of a hot body ; X／sub c／ Y／sub c／=0.5, 0.2 shows the most strong and at X／sub c／, Y／sub c／=0.2, 0.5 the most weak.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.23-30
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Furthermore, numerical simulations for supersonic shock boundary layer interaction with a bleed region were performed and their results were compared with the existing computational results.

• Journal of Korea Society of Industrial Information Systems
• /
• v.4 no.1
• /
• pp.102-109
• /
• 1999

In this paper, Electromagnetic scattering problem by a perfectly conducting strip grating with 2 dielectric layer on a grounded plane by incidence of a electric wave is analyzed by applying the PMM (Point Matching Method) known as a simple procedure. The scattered electromagnetic fields are expanded in a series of Floquet mode functions. The boundary conditions are applied to obtain the unknown field coefficients and the conducting boundary condition is used for the relationship between the tangential electric field and the electric current density on the strip When the incident angle is normal incidence the minimum value of the geometrically normalized reflected power according as relative permittivity is increased it should be noted that the value of the strip width gets moved toward high value. Them most energy by a normal incident wave is scattered in direction of the other angles except normal incident angle.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.3
• /
• pp.203-210
• /
• 2009

If there are some construction works that affect the stability of buried pipelines, the pipelines should be moved to guarantee their safety. In this paper, modeling methods for analyzing the movement of pipelines were sought, and the step-by-step stress estimation method of moving pipelines was developed. Some factors affecting of pipeline response such as the element type, the element size, boundary modeling, and geometric non-linearity were quantitatively investigated. In addition, some conditions in which accuracy and effectiveness can be compromised in the analysis of long pipelines were identified. A neural network was used to estimate the pipeline stress. The inputs to the neural network included step-by-step displacements, and the output was the resulting stress at each movement step. After training the neural network, it can be used to estimate pipeline stresses at some sub-steps that are not included in the training. A Windows-based stress estimation program was developed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.14-14
• /
• 2021

개수로에서는 반드시 자유수면이 있으며, 따라서 물과 공기와의 마찰은 관수로에 비해 상대적으로 매우 작고, 개수로의 전단응력 분포는 관수로와 달리 근본적으로 비대칭이다. 따라서 전단응력은 수로 바닥이나 측면에서만 작용하게 된다. 이러한 평균 전단응력 개념은 흐름에 의해 경계면 구성재료가 이동하는 이동상 수리학에서 유사이송 능력을 해석하는 기준이 되며, 경계면의 전단응력은 힘으로 표시하여 통상 소류력이라 한다. 이러한 복잡한 유체거동은 하천시설물 설계, 시공 및 관리에 있어서 구성재료의 보호능력에 따라 예상하지 못한 조건에서 쉽게 파손될 수 있다. 국내 하천의 경우 한계유속과 한계소류력에 의해 하천설계에 적용되고 있다. 한계 유속의 경우 간단한 수식에 의해 산정될 수 있지만 실제 하천의 보호능력을 대표하기는 힘들기 때문에 한계소류력이 동시에 고려되어야 한다. 한계소류력은 개수로 흐름에서 복잡하게 발생하는 이차류나, 난류 특성에 의해 산정하거나 예측하기는 매우 어렵다. 한계 소류력 뿐만 아니라 하천을 구성하는 재질의 조도계수 역시 균일하지 못하고 매우 예측하기 어렵다. 따라서 본 연구에서는 이러한 복잡한 양상을 나타내는 수리학적 요소에 대해 표준화된 실험수로에서 실험을 통해 평가하고, 체계화된 설계 지침이 되고자 연구를 진행하였다. 본 연구에서는 자연하천과 유사한 조건의 경사를 가지는 경사수로와 경사의 영향에서 자유롭게 평가를 진행하고자 기존 연구를 바탕으로 제작된 소류력 측정장치를 이용하였다. 하천의 설계나 평가에 적용되는 평균 소류력 개념은 복잡한 난류흐름에서 평가지표로써 대표하기 힘들기 때문에 유사 하천환경의 바닥에서 발생하는 소류력을 직접 측정하였다. 본 연구에 사용된 장치는 난류유속 u', v'을 이용하여 Reynolds stress산정하여 Total shear stress를 추정하는 기법을 사용하여 검증하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.518-518
• /
• 2023

자연하천에서 오염물질의 혼합 거동은 비균일한 지형학적 요인으로 인해 매우 복잡한 특성을 나타낸다. 일반적으로 오염물질 거동 모델링에서는 수체에서의 혼합을 Fick의 법칙에 따라 유속에 의한 이송과 난류에 의한 확산으로 계산하고, 국부적인 정체현상 등에 의한 non-Fickian 혼합을 야기하는 하천의 특성을 기하학적 지형 형상으로 구현하여 실제 현상에 근접한 혼합 거동을 재현한다. 하지만 계산의 효율성을 위하여 모델링의 차원을 낮추는 경우, 하천의 지형을 경계조건으로 고려할 수 없게 된다. 특히, 1차원 모델링의 경우 하천의 비균일성을 무시하고 1개의 유선으로 간주하며, 이 경우 non-Fickian 물질이동 해석을 위한 추가적인 현상학적 해석이 필요하다. 지난 50년간, non-Fickian 물질이동 해석을 위한 다양한 현상학적 모형이 제시되어 왔다. 하천을 흐름영역과 정체영역으로 구분하고 두 개의 영역 사이의 물질교환 속도를 모델링하거나, Random walk 개념으로 물질이 이동하는 경우와 이동하지 않는 경우를 확률론적으로 모델링하거나, 물질이 정체되었을 때 다시 빠져나오는 시간을 모델링하는 경우가 그 예이다. 본 연구에서는 선행연구에서 제시한 음함수 형태의 현상학적 모형을 기반으로, 수치적 반복계산 없이 상류 경계에서 임의의 형태의 농도곡선(shape-free breakthrough curve)을 갖는 오염물질운(cloud)이 일정 거리를 유하하며 발생하는 변화를 예측할 수 있는 해를 제시한다. 본 연구의 방법론은 추적법(routing procedure)을 활용한 Fickian 혼합 해석, 전달함수(transfer function) 형태의 정체시간분포 해석, 그리고 라플라스 도메인에서의 해석해 유도를 포함한다. 본 연구에서 제시된 해는 2020년 경상북도 김천시에 위치한 감천의 4.5 km 구간에서 수행한 추적자 실험의 현장 자료를 통해 정확도를 검증하여 타당성을 입증하였다.