• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.178-184
• /
• 2002

For economy and convenience, seismic refraction survey is widely used in surveying for large civil engineering work. The purpose of this study is to obtain the numerical responses of various models using Kirchhoff migration, and to analyze its application to the real data processing. Synthetic traveltime curve was calculated by vidale's algorithm, and various models such as 2 or 3 layer model and irregular topography model are tested to simulate the response of real structure. In order to compare the effect of initial velocity model, true velocity models, inversion results by tomography, smooth velocity models are used as an initial guess. The responses of model data show that the algorithm of this study is more sensitive to initial velocity model than the reflection survey, so choosing a suitable initial velocity model will be the most important thing in real data processing.

• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.129-136
• /
• 2017

3D seismic data processing methods such as full waveform inversion or reverse-time migration require 3D wave propagation modeling and heavy calculations. We compared efficiency and accuracy of a Xeon Phi coprocessor to those of a high-end server CPU using 3D frequency-domain wave propagation modeling. We adopted the OpenMP parallel programming to the time-domain finite difference algorithm by considering the characteristics of the Xeon Phi coprocessors. We applied the Fourier transform using a running-integration to obtain the frequency-domain wavefield. A numerical test on frequency-domain wavefield modeling was performed using the 3D SEG/EAGE salt velocity model. Consequently, we could obtain an accurate frequency-domain wavefield and attain a 1.44x speedup using the Xeon Phi coprocessor compared to the CPU.

• Proceedings of the KSRS Conference
• /
• 2008.10a
• /
• pp.293-296
• /
• 2008

Land surface Temperature (LST) is a very useful surface parameter for the wide range of applications, such as agriculture, numerical and climate modelling community. Whereas operational observation of LST is far from the needs of application community in the spatial Itemporal resolution and accuracy. So, we developed split-window type LST retrieval algorithm to estimate the LST from MTSAT-IR data. The coefficients of split-window algorithm were obtained by means of a statistical regression analysis from the radiative transfer simulations using MODTRAN 4 for wide range of atmospheric profiles, satellite zenith angle and lapse rate conditions including the surface inversions. The sensitivity analysis showed that the LST algorithm reproduces the LST with a reasonable quality. However, the LST algorithm overestimates and underestimates for the strong surface inversion and superadiabatic conditions especially for the warm temperature, respectively. And the performance of LST algorithms is superior when satellite zenith angle is small. The accuracy of the retrieved LST has been evaluated with the Moderate Resolution Imaging Spectroradiometer (MODIS) LST data. The validation results showed that the correlation coefficients and RMSE are about 0.83${\sim}$0.98 and 1.38${\sim}$4.06, respectively. And the quality of LST is significantly better during night and winter time than during day and summer. The validation results showed that the LST retrieval algorithm could be used for the operational retrieval of LST from MTSAT-IR and COMS(Communication, Ocean and Meteorological Satellite) data with some modifications.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.

• Atmosphere
• /
• v.16 no.2
• /
• pp.67-83
• /
• 2006

The statistical analysis for the springtime windstorm in Korea shows that Yeongdong region has the highest occurrence frequency during recent 10 years. The objective of this study is to find possible mechanisms for the downslope windstorm formation in the Yeongdong region by using a mesoscale numerical model, WRF. Dynamical process, wave breaking (hereafter WB), is qualitatively investigated as the candidate mechanism for a windstorm event occurred in 5 April, 2005. WB is developed in upper troposphere downstream, since stable air is lifted by the Taebaek mountain. This process can cause and maintain the severe downslope windstorm by drawing the upper flow down to the surface. And the intensified downslope wind leads the hydraulic jump (hereafter HJ) in downstream region. Froude numbers at Chuncheon (upslope side), Seorak Mountain (crest), Yangyang (lee side), and the East Sea (distant downstream position) are estimated by about 0.4, 1.0, 1.6, and 0.6, respectively. This result implies that the accelerated and supercritical (Fr>1) flow adjusts to the ambient subcritical (Fr<1) conditions in the turbulent HJ. In addition, we find the formation of upstream inversion near top level of the mountain cause the intensification of HJ. Experiments to examine the orographic effect on the mechanisms suggest that the magnitudes of WB and HJ are larger in the experiment of higher topography, but there is no significant difference of windstorm magnitude among the experiments. Another important result from these sensitivity experiments is that the intensity of downslope windstorm strongly depends on the magnitude of upper (2~4 km) wind in upstream side.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.45 no.4
• /
• pp.1-8
• /
• 2008

In orthogonal frequency division multiplexing (OFDM) system the time selectivity of wireless channel introduces intercarrier interference (ICI), which degrades system performance in proportion to Doppler frequency. To mitigate the ICI effect, we can generally employ a classical zero-forcing (ZF) equalizer. However, the ZF scheme requires an inverse of a large matrix, which results in prohibitively high computational complexity. In this paper, we propose a low complexity ZF equalization scheme for suppressing the ICI caused by highly time-varying channels in OFDM systems. From the fact that the ICI on a subcarrier is mainly caused by several neighboring subcarriers, the proposed scheme exploits a numerical approximation for matrix inversion based on Neumann's Series (truncated second order). To further improve performance, the partial ICI cancellation technique is also used with reduced complexity. Complexity analysis and simulation results show that the proposed scheme provides the advantage of reducing computational complexity significantly, while achieving almost the same performance as that of the classical ZF a roach.

• Journal of the Korean Geophysical Society
• /
• v.3 no.4
• /
• pp.251-260
• /
• 2000

Applicabilities of two numerical methods, the 2-dimensional and the 3-dimensional method, are evaluated to inverse test results obtained from the underwater SASW(Spectral -Analysis-of-Surface-Waves) method. As a result of this study, it has been found that the 2-dimensional method can supposed to be applicable for the cases where stiffness of soil layer increases gradually with depth, and the stiffness is relatively low. For the other cases, however, it has been concluded that the 3-dimensional method needs to be applied to determine realistic theoretical dispersion curves. An example is also shown that in situ soil profile underwater is estimated from experimental dispersion curves using the 3-dimensional method. As a results, it can be concluded that the underwater SASW method can be effectively applied to explore the underwater soil condition.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.59-69
• /
• 2001

Five kinds of modified electrode arrays were proposed to overcome the weak points of the commonly used arrays using dipole and/or pole in two-dimensional resistivity surveys. The modified pole-pole array was suggested to overcome the inefficiency caused by distant earthing in pole-pole array. Four kinds of modified arrays using dipole were designed to enhance the signal-to-noise ratio of the conventional dipole-dipole and pole-dipole arrays through boosting up the measured potential difference. In the numerical experiments using the two-dimensional modeling and inversion, the effects of the ambient electrical noise and the resolving power were examined and the results showed the validity of the modified arrays proposed in this study.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.18 no.4
• /
• pp.1042-1058
• /
• 2024

Terahertz (THz) communication is becoming a key technology for future 6G wireless networks because of its ultra-wide band. However, the implementation of THz communication systems confronts formidable challenges, notably beam splitting effects and high computational complexity associated with them. Our primary objective is to design a hybrid precoder that minimizes the Euclidean distance from the fully digital precoder. The analog precoding part adopts the delay-phase alternating minimization (DP-AltMin) algorithm, which divides the analog precoder into phase shifters and time delayers. This effectively addresses the beam splitting effects within THz communication by incorporating time delays. The traditional digital precoding solution, however, needs matrix inversion in THz massive multiple-input multiple-output (MIMO) communication systems, resulting in significant computational complexity and complicating the design of the analog precoder. To address this issue, we exploit the characteristics of THz massive MIMO communication systems and construct the digital precoder as a product of scale factors and semi-unitary matrices. We utilize Schatten norm and Hölder's inequality to create semi-unitary matrices after initializing the scale factors depending on the power allocation. Finally, the analog precoder and digital precoder are alternately optimized to obtain the ultimate hybrid precoding scheme. Extensive numerical simulations have demonstrated that our proposed algorithm outperforms existing methods in mitigating the beam splitting issue, improving system performance, and exhibiting lower complexity. Furthermore, our approach exhibits a more favorable alignment with practical application requirements, underlying its practicality and efficiency.