• Journal of Astronomy and Space Sciences
• /
• v.10 no.2
• /
• pp.152-162
• /
• 1993

We developed the S/W which calculate the Planetary and the Moon ephemerides. The ephemeris of the Solar system objects was obtained from a simultaneous numerical integration of the equations of motion for the nine planets and the Moon. The mathematical model includes contributions from (1) point-mass interactions (2) figure effect (3) earth tides (4) the orientations of the Earth and the Moon. The calculated ephemerides are compared with DE200 data produced by JPL (Jet Propulsion Laboratory).

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.180-186
• /
• 2003

We have developed an one-dimensional (ID) inversion program that can invert multiple frequency small-loop EM data from horizontal coplanar (HCP) and vertical coplanar (VCP) configurations. The inverse problem is solved using least-squares method with active constraint balancing (ACB) method and Jacobian matrix is calculated analytically. Tests using synthetic data from simple ID models indicate that conductivity and depth of each layer can be estimated properly when both real and imaginary data are used together.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.109-117
• /
• 2010

To enable high-performance computing (HPC) for applications with large datasets using a $Sony^{(R)}$ $PLAYSTATION^{(R)}$ 3 ($PS3^{TM}$) video game console, we configured a hybrid system consisting of a $Windows^{(R)}$ PC and a $PS3^{TM}$. To validate this system, we implemented the real-time multiplet identifier (RTMI) application, which identifies multiplets of microearthquakes in terms of the similarity of their waveforms. The cross-correlation computation, which is a core algorithm of the RTMI application, was optimised for the $PS3^{TM}$ platform, while the rest of the computation, including data input and output remained on the PC. With this configuration, the core part of the algorithm ran 69 times faster than the original program, accelerating total computation speed more than five times. As a result, the system processed up to 2100 total microseismic events, whereas the original implementation had a limit of 400 events. These results indicate that this system enables high-performance computing for large datasets using the $PS3^{TM}$, as long as data transfer time is negligible compared with computation time.

• Geophysics and Geophysical Exploration
• /
• v.4 no.1
• /
• pp.1-7
• /
• 2001

The application of geophysical survey methods need to be integrated to meet the increasing demands of imaging of the subsurface in the practical application of civil engineering, underground water survey and environmental problems. This paper examines the IP survey which can be surveyed simultaneously with DC resistivity survey. In this study, 3-D IP modeling algorithm was developed. The 3-D IP modeling algorithm was based on 3-D resistivity modeling by finite-element method. The result of 3-D modeling was compared with 2-dimensional modeling result. The result showed that the 3-D modeling algorithm developed in this study was accurate. Finally, the 3-D modeling algorithm developed in this paper will be useful for the study of IP data.

• Journal of the Korean earth science society
• /
• v.24 no.5
• /
• pp.411-419
• /
• 2003

The assessment of environmental impact on NO$_2$ (or TSP) emitted by vehicles is important for local residents from the point of view of their health and environmental protection. In the course of field investigations, correct concentrations are measured and meteorological data are observed for numerical simulation. To determine background concentration for numerical simulation, annual average concentrations of NO$_2$ (or TSP) are estimated using the Puff-Plume model. If the estimated result affects the environment, it must be considered in the environmental conservation activity. To make the process of a estimation of environmental assessment more easily, this system is developed. Moreover, this system was supplied a Graphic User Interface (GH) for the user who calculated the concentration of air pollution exhausted from the traffic on general roads except special roads such as interchanges and entrances to tunnels. This system can offer not only the numerical result but also a graphic display. Even a beginner who is not a professional programmer can calculate the result easily.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.75-80
• /
• 2009

Despite its flexibility to complex geometry, three-dimensional (3D) electromagnetic(EM) modeling schemes using finite element method (FEM) have been faced to practical limitation due to the resulting large system of equations to be solved. An efficient 3D FEM modeling scheme has been developed, which can adopt either direct or iterative solver depending on the problems. The direct solver PARDISO can reduce the computing time remarkably by incorporating parallel computing on multi-core processor systems, which is appropriate for single frequency multi-source configurations. When limited memory, the iterative solver BiCGSTAB(1) can provide fast and stable convergence. Efficient 3D simulations can be performed by choosing an optimum solver depending on the computing environment and the problems to be solved. This modeling includes various types of controlled-sources and can be exploited as an efficient engine for 3D inversion.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.154-160
• /
• 2007

This paper compares the characteristics of three different algorithms for three-dimensional (3D) magnetotelluric (MT) modeling. These methods are developed by Mackie et al. (1994), Sasaki (1999) and Nam et al. (2007). The first and second methods are based on the finite difference method (FDM), while the last one the finite-element method (FEM). MT responses, apparent resistivities and phases, for a COMMEMI 3D-2 model show a good agreement with integral equation solutions and only minor discrepancies are found over the anomalous bodies in the 3D model. The computation time of the two methods based on FDM is short and the static divergence correction contributes to speed up. The FEM modeling scheme is accurate but slow.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.185-190
• /
• 2011

In marine controlled-source electromagnetic (CSEM) modeling, it may be difficult to evaluate primary fields accurately using conventional linear filters because they decay very rapidly with distance. However, since there exists a closed-form solution to the Hankel transform in TM mode for a homogeneous half space, we can assess the accuracy of linear filters for evaluating the Hankel transform. As a result, only nine out of 36 source-receiver pairs show that EM fields decrease linearly in semi-log scale with an increase of source-receiver distance, while EM fields are either 0 or not reduced significantly due to an effect of the air layer. There also exist closed-form solutions for the nine pairs, and the others can be evaluated accurately with a relatively short filter. This paper proposes a method which uses closed-form solutions for TM-mode Hankel transforms and a filter with 61 coefficients for TE-mode ones.

• 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.18 no.3
• /
• pp.133-142
• /
• 2015

Displacement of the Earth's surface due to atmospheric loading has been recognized since a century years ago, and its accurate estimation is required in present day geodesy and surveying, particularly in space geodesy. Atmospheric load deformation in continental region can readily be calculated with the given atmospheric pressure field and the load Green's function, and, in near coastal area, approximate model is used for the calculation. The changes in the Earth's atmospheric circulation and the seasonal variation of atmospheric pressure on two hemispheres of the Earth are the each main causes of variation of the Earth's spin angular velocity and polar motion respectively. Wind and atmospheric pressure do the major role in other periodic and non-periodic perturbations of the positions in the Earth's reference frame and variations in the Earth's spin rotational state. In this reviewing study, the developments of related theories and models are summarized along with brief description of phenomena, and the geodetic perturbing effects of a hypothetical typhoon passing Korea are shown as an example. Finally related existing problems and further necessary studies are discussed in general.