• 지구물리와물리탐사
• /
• 제9권3호
• /
• pp.241-249
• /
• 2006

한국의 지질 환경은 암석 분포가 매우 다양하고 복잡한 구조 활동의 영향을 받아 지하매질의 이방성 특성이 국부적으로 심하게 변화한다. 기존의 이방성 주시 모델링의 경우 지질 모델을 2차원으로 단순화시킴으로써 이러한 복잡한 지질 환경을 제대로 고려할 수 없었다. 또한 약 이방성 가정을 사용하여 실제로 나타날 수 있는 지하 매질의 심각한 이방성 영향을 주시 모델링에서 고려할 수 없었다. 이에 이 연구에서는 보다 실제적이고 복잡한 3차원 횡등방성 매질(transversely isotropic media)에서 q-P파의 초동 주시 양상을 모사할 수 있는 주시 모델링 알고리듬을 개발하였다. 이 알고리듬에서는 2차원 비선형 주시 내삽(2D nonlinear traveltime interpolation) 기법과 주시의 3차원 격자 채움법(mapping)을 이용한 직접 전파법(direct calculation)을 통해 급격한 물성의 변화에도 주시 계산이 가능하도록 하였다. 또한, 최소 주시 계산과정에서 수치 미분을 통한 최대 경사법(steepest descent method)을 사용하여 약 이방성 가정을 극복하였다. 개발된 알고리듬은 해석해와 비교하여 그 타당성을 검증하였고 3차원 2층구조에 대한 주시 계산을 수행하여 물성이 급격히 변화하는 모델에 대해서도 안정적으로 주시 계산이 이루어짐을 확인하였다. 이 연구에서 개발한 3차원 주시 모델링 알고리듬은 향후 구조보정이나 토모그래피 알고리듬 개발에 사용될 수 있을 것으로 기대한다.

• 한국전산유체공학회지
• /
• 제20권2호
• /
• pp.46-51
• /
• 2015

In this paper, we have developed numerical model for particulated flow through narrow slit using Eulerian-Granular method. Commercial software (FLUENT) was utilized as simulation tool and main focus was to identify the effect from various numerical options for modeling of solid particles as continuos phase in granular flow. Gidaspow model was chosen as basic model for solid viscosity and drag model. And lun-et-al model was used as solid pressure and radial distribution model, respectively. Several other model options in FLUENT were tested considering the cross related effect. Mass flow rate of the particulate through the slit was measured to compare. Due to the high volume density of the stacked particulates above the slit, effect from various numerical options were not significant. The numerical results from basic model were also compared with experimental results and showed very good agreement. The effects from the characteristics of particles such as diameter, angle of internal friction, and collision coefficient were also analyzed for future design of velocity resistance layer in solar thermal absorber. Angle of internal friction was found to be the dominat variable for the particle mass flow rate considerably. More defined 3D model along with energy equation for complete solar thermal absorber design is currently underway.

• 대한토목학회논문집
• /
• 제33권3호
• /
• pp.965-974
• /
• 2013

2011년 7월 27일 우면산 유역에서 일련의 토석류가 발생하여 많은 인명과 재산 피해를 야기했다. 이러한 토석류 피해를 저감하기 위해서는 토석류가 어떻게 이동하고 퇴적되어 하류부에 피해를 주는지를 예측할 수 있는 기술이 필요하다. 이 연구에서는 우면산 토석류에 대해서 FLO-2D 수치모형을 적용하여, 모형의 성능을 평가하고, 이들 토석류의 유동 매개변수를 분석하였다. 우면산 북측사면에서 발생한 토석류의 현장 관측자료와 FLO-2D 모의 결과 분석을 통해서 역산된 우면산 토석류의 대표 항복응력은 1,022 Pa 그리고 대표 동점성은 652 $Pa{\cdot}s$ 인 것으로 나타났다. 선정된 매개변수를 이용하여 우면산 토석류의 퇴적영역을 산정한 결과, 래미안 유역의 경우 63 - 85% 정확도의 피해유역을 재현됨으로서 FLO-2D 모형은 피해유역 산정에 적합한 것으로 나타났다. 하지만, FLO-2D 모형이 침식물질과 우수의 연행작용을 고려할 수 없는 고유의 한계 때문에 연행작용이 현저한 래미안 및 신동아 유역 토석류의 유하속도는 과소평가하는 것으로 나타났다. 유하부에서 침식에 의한 연행작용이 활발하지 않은 형촌마을 유역 토석류는 수로와 저수지 상부 부근에서의 퇴적을 우수한 정확도로 모의되었다. 민감도 분석 결과 유사체적농도와 바닥조도계수 모두 토석류의 흐름두께와 유속 모두에 상당한 영향을 미치는 것으로 나났으므로, 이들 매개변수의 선정에 신중해야 함을 보여주었다.

• Structural Engineering and Mechanics
• /
• 제33권5호
• /
• pp.621-633
• /
• 2009

In this paper we compare the solution of the one-dimensional beam model and the numerical solution of the two-dimensional linearized elasticity problem for rectangular domain of the beam-like form. We first derive the beam model starting from the two-dimensional linearized elasticity, the same way it is derived from the three-dimensional linearized elasticity. Then we present the numerical solution of the two-dimensional problem by finite element method. As expected the difference of two approximations becomes smaller as the thickness of the beam tends to zero. We then analyze the applicability of the one-dimensional model and verify the main properties of the beam modeling for thin beams.

• 한국환경과학회지
• /
• 제20권9호
• /
• pp.1151-1163
• /
• 2011

The bed change analysis near the opening gate of a dam or weir to release deposited sediments have been conducted mostly using the numerical models. However, the use of unverified input parameters in the numerical model is able to produce the different results with natural and real conditions. Also, the bed changes near the opening gate of a dam or weir calculated with a numerical model could be varied depending on the geometry extent included the downstream area with supercritical flow in the model. In addition, the different time steps could provide different results in the bed change calculation, even though other conditions such as input parameters, geometries, and total simulation time were same. Therefore, in this study, hydraulic experiments were performed to validate the eddy viscosity coefficient which is the one of important input parameters in the RMA2 model and relevant to variation of simulation results. The bed changes were calculated using the SED2D model based on flow results calculated in the RMA2 model with the verified and selected eddy viscosity coefficient and also compared with experimental results. The bed changes near the opening gate were underestimated in the numerical model comparing with experimental results except only the numerical case without the modeling section of sediment release pipe and downstream area where the supercritical flow was produced. For the simulation of minimum time steps, different shapes of scour hole were produced in numerical and physical modeling.

• Structural Engineering and Mechanics
• /
• 제48권2호
• /
• pp.151-171
• /
• 2013

The current paper presents the numerical blind prediction of nonlinear seismic response of two full-scale, three dimensional, one-story reinforced concrete structures subjected to bidirectional earthquake simulations on shaking table. Simulations were carried out at the laboratories of LNEC (Laboratorio Nacional de Engenharia Civil) in Lisbon, Portugal. The study was motivated by participation in the blind prediction contest of shaking table tests, organized by the challenge committee of the 15th World Conference on Earthquake Engineering. The test specimens, geometrically identical, designed for low and high ductility levels, were subjected to subsequent earthquake motions of increasing intensity. Three dimensional nonlinear analytical models were implemented and subjected to the input base motions. Reasonably accurate reproduction of the measured displacement response was obtained through appropriate modeling. The goodness of fit between analytical and measured results depended on the details of the analytical models.

• Computers and Concrete
• /
• 제18권3호
• /
• pp.319-336
• /
• 2016

In this study, both two- and three-dimensional (2D and 3D) finite-volume-based models were developed to analyze the heat transfer mechanisms through the porous structures of cellular concretes under steady-state heat transfer conditions and to investigate the differences between the 2D and 3D modeling results. The 2D and 3D reconstructed pore networks were generated from the microstructural information measured by 3D images captured by X-ray computerized tomography (X-CT). The computed effective thermal conductivities based on the 2D and 3D calculations performed on the reconstructed porous structures were found to be nearly identical to those evaluated from the 2D cross-sectional images and the 3D X-CT images, respectively. In addition, the 3D computed effective thermal conductivity was found to agree better with the measured values, in comparison with the 2D reconstruction and real cross-sectional images. Finally, the thermal conductivities computed for different reconstructed porous 3D structures of cellular concretes were compared with those obtained from 2D computations performed on 2D reconstructed structures. This comparison revealed the differences between 2D and 3D image-based modeling. A correlation was thus derived between the results of the 3D and 2D models.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회A
• /
• pp.222-227
• /
• 2008

This study is focused on an integrated numerical modeling enabling one to investigate the dynamic behavior and failure of 2-D textile composite and 3-D orthogonal woven composite structures weakened by micro-cracks and subjected to an impact load. The integrated numerical modeling is based on: I) determination of governing equations via a three-level hierarchy: micro-mechanical unit cell analysis, layer-wise analysis accounting for transverse strains and stresses, and structural analysis based on anisotropic plate layers, II) development of an efficient computational approach enabling one to perform transient response analyses of 2-D plain woven and 3-D orthogonal woven composite structures featuring the matrix cracking and exposed to time-dependent loads, III) determination of the structural characteristics of the textile-layered composites and their degraded features under various geometrical yarn shapes, and finally, IV) assessment of the implications of stiffness degradation on dynamic response to impact loads.

• 한국터널지하공간학회 논문집
• /
• 제19권5호
• /
• pp.799-812
• /
• 2017

쉴드 TBM 터널에 대한 2차원 수치해석 시 기존 연구에서는 뒤채움재 및 실제적인 내공변위 분포를 적절히 고려하지 않았다. 따라서 본 연구에서는 이를 보완하고자 gap parameter 적용 시 굴착면에서 발생하는 실제적인 내공변위분포를 고려하고, 적정한 뒤채움재 물성치를 적용한 모델링 방법을 제시하였다. 이를 위해 단일층 풍화토 지반을 대상으로 gap parameter와 토피고가 굴착면의 지반손실량과 지표손실량으로 구성되는 체적손실량에 미치는 영향을 분석하였다. 그 결과 대부분의 지표침하가 굴착 직후에 발생하였으며, 토피고 및 gap parameter가 증가할수록 수치해석을 통해 얻은 지표침하곡선의 영향범위가 제안된 이론식보다 넓게 나타났다. 따라서 2차원 수치해석 시 쉴드 TBM 터널의 거동을 정확히 평가하기 위해서는 하중분담률을 이용하여 gap parameter를 적용하고, 뒤채움재의 물성치도 적절히 고려해야 한다고 판단된다.

• 지질공학
• /
• 제26권1호
• /
• pp.1-13
• /
• 2016

Groundwater flow paths and groundwater ages at a radioactive waste repository located in a coastal area of South Korea were evaluated using the hydrochemical and hydrogeological characteristics of groundwater, surface water, rain water, and seawater, as well as by numerical modeling. The average groundwater travel time in the top layer of the model, evaluated by numerical modeling and groundwater age (34 years), approximately corresponds to the groundwater age obtained by chlorofluorocarbon (CFC)-12 analysis (26-34 years). The data suggest that the groundwater in wells in the study area originated up-gradient at distances of 140-230 m. Results of CFC analyses, along with seasonal variations in the δ¹⁸O and δD values of groundwater and the relationships between ²²²Rn concentrations and δ¹⁸O values and between ²²²Rn concentrations and δD values, indicate that groundwater recharge occurs in the summer rainy season and discharge occurs in the winter dry season. Additionally, a linear relationship between dissolved SiO₂ concentrations and groundwater ages indicates that natural mineralization is affected by the dilution of groundwater recharge in the rainy summer season.