• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• pp.198-200
• /
• 2002

경제적이며 정확한 3차원 전자탐사 모델링을 위해 위해 Habashy et al. (1993)에 의해 제안된 국소 비선형 근사(localized nonlinear approximation)를 이용하여 전자탐사 모델링 알고리듬을 개발하였다. 전자탐사 수치모델링시 많은 계산시간 및 기억용량을 필요로 하는 Green 텐서 적분을 정확하고 빠르게 계산하기 위해, 단일 미소체를 이용한 공간파수 영역에서의 Green 텐서 적분 알고리듬을 제안하였다. 더욱이 Green 텐서의 송수신 방향 및 상반성을 고려하여 각각의 미소체에 의한 전체 미소체에의 Green 텐서 적분을 한 개의 미소체에 의한 전체 미소체에의 Green 텐서 적분 값으로 구하게 하므로 매우 적은 기억용량 만으로 Green 텐서 적분 행렬을 구성할 수 있어, 역산법에 효과적으로 적용할 수 있다. 이 수치 모델링 알고리듬을 기본으로 하여 평활화 제한을 가한 최소자승 역산 알고리듬을 개발하였다. 이 역산 알고리듬을 지표 전자탐사 및 시추공-지표 전자탐사 등에 적용하여 PC에서도 빠르게 3차원 전자탐사 역산이 수행됨을 보였다.

• Journal of the Korean Geophysical Society
• /
• v.5 no.1
• /
• pp.63-71
• /
• 2002

The integral equation method is a powerful tool for electromagnetic numerical modeling. But the difficulty of this technique is the size of their linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. To overcome this limitation, a lot of approximation and series methods, such as conventional Born, modifed Born and extended Born, were developed. But all the methods need volume integration of Green tensor, which is very time consuming. In electromagnetic theory, Green tensor rapidly decreases as the distance between source and field cell increases. Therefore, the source cell which are far away from the field cell does not make an effect on the electric field of the field cell. Consequently, by ignoring the effect of Green tensor due to far away source cells, computing time for electromagnetic numerical modeling can be reduced dramatically. Comparisons of this new method against a full integral equation, extended Born approximation and series code show that the method is accurate enough much less time consuming.

• Geophysics and Geophysical Exploration
• /
• v.26 no.1
• /
• pp.1-7
• /
• 2023

This study derives the expressions of vector gravity and gravity gradient tensor due to an elliptical cylinder. The vector gravity for an arbitrary three-dimensional (3D) body is obtained by differentiating the gravitational potential, including the triple integral, according to the shape of the body in each axis direction. The vector gravity of the 3D body with axial symmetry is integrated along the axial direction and reduced to a double integral. The complex Green's theorem using complex conjugates subsequently converts the double integral into a one-dimensional (1D) closed-line integral. Finally, the vector gravity due to the elliptical cylinder is derived using 1D numerical integration by parameterizing a boundary of the elliptical cross-section as a closed line. Similarly, the gravity gradient tensor due to the elliptical cylinder is second-order differentiated from the gravitational potential, including the triple integral, and integrated along the vertical axis direction reducing it to a double integral. Consequently, all the components of the gravity gradient tensor due to an elliptical cylinder are derived using complex Green's theorem as used in the case of vector gravity.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.51-56
• /
• 2024

In this paper, the vector gravity and gravity gradient tensor of an elliptical disk are derived. The vector gravity of an elliptical disk is defined by differentiating the gravitational potential due to the elliptical disk expressed by a double integral with respect to each axial direction. The vector gravity defined by the double integral is then transformed into a line integral of a closed curve along the elliptical disk boundary using the complex Green's theorem. Finally, vector gravity due to the elliptical disk is derived by 1D parametric numerical integration along the elliptical disk boundary. The xz, yz, zz components of the gravity gradient tensor due to the elliptical disk are obtained by differentiating the vector gravity with respect to vertical direction. The xx, yy, xy components are derived by differentiating the horizontal components of the vector gravity in the form of a double integral with respect to horizontal directions and then using the complex Green's theorem.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.13-18
• /
• 2000

We describe the concept of the singularity in the integral equation of electromagnetic (EM) modeling in comparison with that in the integral representation of electric fields in EM theory, which would clarify the singular integral problems of the Green's tensor. We have also derived and classified the singular integrals of the Green's tensors in 3-D, 2.5-D and 2-D as well as in the thin sheet integral equations of the EM scattering problem, which have the most important effect on the accuracy of the numerical solution of the problems.

• The Korean Journal of Rheology
• /
• v.8 no.3_4
• /
• pp.215-225
• /
• 1996

평판형의 열가소성 수지를 유리전이온도 이상으로 가열한 다음 압력을 가함으로써 원하는 형상의 제품을 성형하는 열성형공정은 대상 수지가 큰 변형을 일으킬 뿐만 아니라 비선형적 거동을 보이게 된다. 따라서 수지의 변형거동 예측과 최적성형조건의 설정에 많은 어려움과 시행착오를 거치게 되는 바, 열성형 공정의 최적화를 위한 연구의 일환으로 원형 평판위 수지를 대상으로 수지의 부풀림 거동과 이에 따른 두께 분포를 예측할수 있는유한요 소법의 수치모사 알고리듬을 개발하고자 하였다. Piola-Kirchhoff 응력 텐서와 Green 변형 텐서 및 lagrangian 변형 텐서를 사용하여 평판상의 응력-변형에 대한 비선형의 에너지 수 지식을 수립하고 Newton-Raphson 반복수렴법을 이용하여 근사적으로 해석하였으며 수지의 유변학적 구성방정식으로는 neo-Hookean 모델, Mooney-Rivlin 모델 및 Ogden 모델등의 초탄성 모델을 사용하여 그결과를 비교하였다. 수치모사에는 두께가 매우 얇기 때문에 두께 방향의 응력변화를 무시할수 있는 membrane 가정을 도입한 2차원적 해석과 두께 방향의 응력 변화를 고려하는 3차원적 해석을 모두 수행하고 그 차이를 비교하였으며 3차원적 해석 의 경우에는 penalty법을 이용하여 비 압축성을 만족하였다. 일차적으로 내압을 받는 두꺼 운 원통계에 대한 수치모사 해석을 수행하고 완전해와 비교함으로써 개발된 수치모사 알고 리듬의 열성형 공정에의 적용 타당성을 검증하였으며 이를 이용하여 원형 평판의 자유부풀 림거동을 예측한 결과 Treloar 등의 실험결과와 잘 부합함을 확인하였다. 또 간단한 형상의 금형이 있는 경우와 반지름 방향으로의 온도변화에 따른 수지의 변형거동을 해석함으로써 실제 열성형 공정에서 요구되고 있는 성형품의 두께 분포를 균일하게 하기 위한 방안을 제 시하였다.

• Economic and Environmental Geology
• /
• v.27 no.1
• /
• pp.41-47
• /
• 1994

A fast Hankel transform (FHT) algorithm (Anderson, 1982) is applied to numerical evaluation of many Green's tensor integrals encountered in three-dimensional electromagnetic modeling using integral equations. Efficient computation of Hankel transforms is obtained by a combination of related and lagged convolutions which are available in the FHT. We express Green's tensor integrals for a layered half-space, and rewrite those to a form of related functions so that the FHT can be applied in an efficient manner. By use of the FHT, a complete or full matrix of the related Hankel transform can be rapidly and accurately calculated for about the same computation time as would be required for a single direct convolution. Computing time for a five-layer half-space shows that the FHT is about 117 and 4 times faster than conventional direct and multiple lagged convolution methods, respectively.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.64-67
• /
• 2002

• Journal of the Korean Society of Hazard Mitigation
• /
• v.1 no.3 s.3
• /
• pp.123-136
• /
• 2001

The purpose of application of moment tensor inversion method is to determine source parameters, such as, focal mechanism, seismic moment and source depth. This paper presents results of focal mechanism solutions of 14 recent events with magnitudes ranging from M3.3 to M4.8 by using moment tensor inversion method called TDMT_INV. The strike of faults is in the direction of NE-SW and NW-SE with the movement of strike-slip or strike-slip of minor reverse component. The compressional axis of the stress field is predominantly E-W or ENE-WSW except for some faults, which are distributed at Ryongnam Massif and Wonsan, they have a compressional axis of N-S or NNW-SSE.

• The Journal of Engineering Geology
• /
• v.17 no.4
• /
• pp.673-680
• /
• 2007

The seismic source parameters of the Odaesan earthquake on 20 January 2007, including focal depth, focal mechanism, magnitude, and source characteristics, are analysed using seismic moment tensor inversion. The Green's function for different 3 crust models representing the southern Korean Peninsula are used. Final results show that the event, considering 6 seismic moment tensor elements, is caused by the typical strike slip fault with nearly NNE strike. The focal depth is estimated to be about 11km and 6 seismic moment tensor elements with 7.2% CLVD value shows typical double couple seismic source. The consistent characteristics of the strike and epicenter of the event with Odaesan fault imply that Odaesan earthquake is mainly caused by movement of the Odaesan fault.