• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.402-407
• /
• 2006

본 연구는 곡선좌표계에서 유한차분기법(finite difference method)을 이용하여 2차원 흐름이 모의가능한 수치모형을 개발하는 것이다. 기존의 연구는 대부분 직교좌표계(cartesian coordinate system)에서의 격자망을 대상으로 개발되고 적용되었기 때문에 불규칙한 흐름의 경계 및 형상을 올바로 표현하기 어려웠다. 유한요소법이나 유한체적법같은 수치모의기법들이 개발되어 비구조격자체계를 구성하고 자연현상에 가까운 경계 표현할 수 있도록 개발되었다. 하지만 위의 기법들은 질량과 운동량과 같은 물리량을 보존하기 위해서 매우 조밀한 격자체계를 가져야만 한다. 이에 본 연구에서는 기존의 문제점들을 해결하기 위하여 곡선좌표계(curvilinear coordinate system)를 이용하여 지배방정식을 표현하고 2차원 흐름을 모의할 수 있는 모형을 구축한다. 수치모형은 leap-frog기법과 1차 정확도의 풍상차분기법(upwind scheme)을 사용하여 구성하였다. 본 연구에서 개발된 모형을 사각수조 및 만곡수로흐름에 적용하여 모의결과를 해석해 및 실험관측값과 비교하였다. 이로부터 본 수치모형이 해석해 및 실측치와 잘 일치하고 있음을 알 수 있었다.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05b
• /
• pp.707-712
• /
• 1996

중대사고시 원자로 압력용기내 또는 원자로 공동(cavity) 내에서의 노심용융물은 주입되는 물로 인하여 물과 접촉하는 표면이 냉각되면서 피막층(crust)이 형성된다. 이러한 피막층의 형성은 노심용융물과 냉각수 사이의 열전달 현상에 영향을 미치며 중대사고 발생시 사고 진행에 중요한 역할을 한다. 본 연구에서는 이러한 용융물의 피막층 형성의 해석모델을 수립하기 위해 전이현상과 전도와 대류를 포함하는 2차원 열전달과 상변화를 수반하는 문제를 포함하는 운동량방정식과 에너지방정식을 2차원으로 구성하였으며 에너지방정식은 엔탈피의 함수로 나타내었다. 그리고 이러한 2차원 지배방정식을 해석하기 위해 유한차분법 및 SIMPIER 알고리즘을 이용하였다. 비교대상으로는 한국원자력연구소에서 수행한 냉각수의 비등과 기체주입 효과가 고려되지 않은 실험을 대상으로 하였다. 계산결과 용융물의 피막층은 파동(wave) 형태로 형성되었으며 일정시간이 경과하면 변화가 없는 안정한 상태가 되었다. 용융물 내에서의 온도분포는 액체상태일 경우에는 하부가열면과 상변화가 일어나는 경계면부근을 제외하고는 거의 일정한 온도분포를 나타내고 있으며 용융물이 고화된 피막층에서는 급격한 온도변화를 보여주고 있다.

• Economic and Environmental Geology
• /
• v.29 no.3
• /
• pp.345-355
• /
• 1996

The effectiveness of various electrode configurations in horizontal mappings and 1-D inversions of vertical sounding data for delineating 2-D structures was studied. Apparent resistivity values of three point, dipole-dipole, Wenner, and Schlumberger mappings were simulated for such structures as vertical dyke, tabular prism, buried vertical fault, ramp and complex structure by finite difference method (FDM) and they were compared with each other. Also 2-D cross sections for three structures obtained by interpolation of 1-D inverted sounding data in terms of three layers were compared for Schlumberger and Wenner arrays. On these cross sections, horizontal and vertical resistivity interfaces of the 2-D structures are revealed relatively clearly. Apparent resistivity curves of Schlumberger mapping show vertical resistivity discontinuities very well. On the whole, Schlumberger array is superior to the other arrays in electric sounding as well as mapping. This study clearly indicates that interpretations of 2-D structures based on 1-D inversion are possible.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.4
• /
• pp.839-849
• /
• 1990

본 연구에서는 위와 같은 열 점소성 해석의 여러 수치해석 기법을 첫째, 변형 과 온도 해석을 연계 방식과 비연계 방식으로 계산한 후 결과를 검토하고, 둘째, 온도 해석을 유한 요소법과 유한 차분법으로 계산한 후 각각의 장단점과 효율적 방법을 검 토하는데 목적이 있다.또한 이 결과로 계속적인 3차원 열 점소성 해석 연구수행의 기초를 마련하고자 한다.

• Geophysics and Geophysical Exploration
• /
• v.22 no.3
• /
• pp.99-106
• /
• 2019

Resistivity method has been used for the dam safety inspection and, for the convenience of fieldwork, two-dimensional (2D) resistivity data has been usually measured along the dam crest. However, since the dam has three-dimensional (3D) structure, 2D resistivity survey along the dam crest violates 2D assumption and 3D effects caused by 3D topography and material properties in the dam distort the inversion result of 2D resistivity data acquired along the dam crest. Furthermore, it is really hard to evaluate the 3D structure of the dam and 3D leakage pathway using 2D resistivity survey because 2D resistivity survey can provide only 2D resistivity section beneath the survey line. In this study, 3D resistivity survey was conducted at a dam in Korea. By comparing the results from 3D and 2D resistivity surveys, merit and demerits of 3D survey were investigated. Finally, it was confirmed that 3D survey can provide more accurate information about the dam status and 3D leakage pathway compared to the 2D survey. Therefore the 3D resistivity survey should be actively expanded for more accurate dam safety inspection even though more time and expense are required.

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

본 연구에서 개발한 하천흐름해석모형 HDM-2Di는 2차원 추적모형에 필요한 격자생성기(RAMS-G) 및 GUI (RAMS+)와 연계하여 요소망 생성기능 및 요소망 처리기능을 통해 실제 지형을 반영할 수 있는 전처리 과정을 수행할 수 있다. 장구간 하천에서 신속하고 정확하게 모의결과를 도출할 수 있도록 입출력 체계를 구조화하다. 또한 하천의 복잡한 지형과 자연 하천구조의 동역학적인 흐름환경을 효과적으로 재현할 수 있으며 흐름해석모의 결과인 2차원 평면상의 유속과 수심을 2차원 유해화학물질 추적모형(CTM-2D-TX)의 입력자료로 활용할 수 있도록 개발하였다. HDM-2Di 모형의 성능 검증을 위하여 정상류 및 준부정류, 부정류 조건 흐름해석 성능 테스트를 진행하였으며, 실제 자연하천의 물리적 구조를 재현한 사행수로 추적자 실험 결과와 모형 결과를 비교분석하였다. 또한 자연하천 적용을 위하여 구미보와 칠곡보 구간을 대상으로 ADCP 실측 결과와 상용모형인 Nays-2DH 모의 결과와의 비교를 통해 HDM-2Di 흐름해석모형의 정확성 및 적용성을 평가하였다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.41-44
• /
• 2003

In spite of many researches made on the supersonic jets until now, detailed three-dimensional structures of supersonic jets are not well hewn. In the current study, the detailed structures of three-dimensional supersonic jets are numerically investigated using a CFD method. The total variation diminishing (TVD) scheme is used to solve the unsteady, three-dimensional, compressible Euler equations. Computational results are visualized to investigate the major features of supersonic jets. The three-dimensional computation results show that the structures of the supersonic jets are significantly different from those of the two-dimensional or axisymmetric supersonic jets.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.815-821
• /
• 2009

하천의 2차원 흐름 및 하상변동, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있으며, 국내의 경우 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 널리 사용되고 있다. 하천을 유한 요소화 하는 전처리 과정은 전체 해석 과정을 자동화 하는데 있어 필수적인 요소이며, 주로 삼각 요소망 또는 사각 요소망을 이용하여 해석을 수행하게 된다. 삼각 요소망의 경우 상대적으로 자동화하기 쉬운 반면 사각 요소망의 생성은 절점 생성 자체가 삼각 요소망 보다 더 많은 기하학적 제한 요소를 가지고 있기 때문에 상대적으로 완성도 높은 알고리즘을 구현하기가 어렵다 할 수 있다. 이에 따라 본 연구에서는 2차원 상에서 사각 요소망(quadrilateral elements)을 생성할 수 있는 Paving method를 중심으로 한 요소망 생성 알고리즘에 대해 고찰하고, 국내 최초의 범용 수치해석 모형인 RAMS(River Analysis and Modeling System)에 적용하였다. Paving method는 1990년에 Blacker and Stephenson에 의해 제안되었으며, Sandia National Laboratories에 의해 완성되었다. Paving Method는 advancing front style의 요소망을 생성하게 되고, 바깥쪽에서 안쪽으로 element layer를 생성하면서 채워나간다. 본 연구에서는 기존의 요소망 생성 프로세스에서 element 삽입 전의 검증 기능을 강화한 새로운 버전의 paving method를 적용하엿다.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.6
• /
• pp.38-46
• /
• 2011

In order to simulate the ignition-gas injection in the interior ballistics, a two-dimensional analytic code for two-phase flows has been developed. The Eulerian-Lagrangian approach and the low-dissipation simple high-resolution upwind scheme(LSHUS) have been adopted in the numerical code for the propellant combustion of the gun propelling charges. The ghost-cell extrapolation method has been used for the moving boundary in the chamber with the projectile movement. The calculation results of the developed code have been compared and verified through those of the dimensionless IBHVG2 code and the previous one-dimensional code. In comparison with the two-phase flows according to the drag models, the numerical analysis of the muzzle velocity has been affected by the drag model.

• Journal of the Korean earth science society
• /
• v.30 no.6
• /
• pp.699-708
• /
• 2009

Recently as the interest in the development of domestic ore deposits has increased, we can easily find some studies on exploration geophysics-based ore-deposit survey in literature. Based on the fact that mineralized zone are generally more conductive than surrounding media, electrical resistivity survey among several geophysical surveys has been applied to investigate metallic ore deposits. Most of them are grounded on 2-D survey. However, 2-D inversion may lead to some misinterpretation for 3-D geological structures. In this study, we investigate the feasibility of the 3-D electrical resistivity survey to 3-D vein-type ore deposits. We first simulate 2-D dipole-dipole survey data for survey lines normal to the strike and 3-D pole-pole survey data, and then perform 3-D inversion. For 3-D ore-body structures, we assume a width-varying dyke, a wedge-shaped, and a fault model. The 3-D inversion results are compared to 2-D inversion results. By comparing 3-D inversion results for 2-D dipole-dipole survey data to 3-D inversion results for 3-D pole-pole survey data, we could note that the 2-D dipole-dipole survey data yield better inversion results than the 3-D pole-pole data, which is due to the main characteristic of the pole-pole array. From these results, we are convinced that if we have certain information on the direction of the strike, it would be desirable to apply 2-D dipole-diple survey for the survey lines normal to the strike. However, in most cases, we do not have any information on the direction of the strike, because we already developed the ore deposit with the outcrops and the remaining ore deposits are buried under the surface. In that case, performing 3-D pole-pole electrical resistivity survey would be a reasonable choice to obtain more accurate interpretation on ore body structure in spite of low resolution of pole-pole array.