• 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.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

• 연구논문집
• /
• s.33
• /
• pp.17-25
• /
• 2003

In this work, thermal design of a PCR chip for LOC is systematically conducted. From the numerical simulation of a PCR chip based on the finite volume method, how to control the average temperature of a PCR chip and the temperature difference between the denaturation zone and the annealing zone is presented. The average temperature is shown to be controlled by adjusting heat input and a cooler as well as a heater is shown to be necessary to obtain three individual temperature zones for polymerase chain reaction. To reduce the time required, a heat sink for the cooler is not included in the calculation domain for the PCR chip and heat sink design is conducted separately by using a compact modeling method, the porous medium approach.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.4
• /
• pp.565-574
• /
• 2000

A new ground penetrating radar imaging method for the estimation of buried artificial structures location and their approximate shapes in dispersive lossy ground is investigated. Fundamental idea is based on estimating delayed time and amplitude retrieval coefficients from scattered signals by buried scatterers. Using absolute value integration of each scanning site not only improve the accuracy of measured scattered signal, but also offers convenient ways to extract the image of buried structures. Multi-term Debye model was employed to describe a dispersive and lossy ground medium. We used the finite difference time domain method to discretize the wave equation in continuous form into the machine suitable form. This imaging method uses a new wave path tracing technique in time domain, which is helpful to identify the exact position of buried structures against the ground surface fluctuations.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.7
• /
• pp.1065-1073
• /
• 1999

This paper analyzed the ground bounce caused by the power plane resonance in the multilayered printed circuit board(PCB) using the Haar-wavelet-based Multiresolution Time-Domain (MRTD). In conventional Finite-Difference Time-Domain(FDTD), the highly fine vertical cell is needed to represent the distance between $V_{cc}$ plane and ground plane since the two planes are very close. Therefore the time step $\Deltat$ must be very small to satisfy the stability condition. As a result, a large number of iterations are needed to obtain the response in wanted time. For this problem, this paper showed that the computation time can be reduced by application of the MRTD method. The results obtained by the MRTD agree very well with those by FDTD method and analytic solutions. In conclusion, this paper proved the efficiency and accuracy of MRTD method for analyzing the EMI/EMC problems in PCB.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.5
• /
• pp.678-686
• /
• 1998

This paper describes an anisotropic perfectly matched layer (APML) for mesh truncation of the Finite Difference Time Domain(FDTD) method. The APML method can classified by a split type and an unsplit type, in case of the split type be made up 12 equations or 8 equations largely, and in case of the unsplit type be made of 6 equations. Therefore the latter is more simple as compare with the former. For the APML method presented in this paper is the latter, is directly derived from the time domain equations of Maxwell and extend to the three dimensional problem for the method suggested by Chen. Especially, in the edge and corner parts except the planes, the APML method effectively treated as compound with the Chen's method and Gedney's method newly. The results of the numerical method in this paper show the radiation patterns and the time responses of electromagnetic fields of the wire antennas according to wavelengths and the APML results are compared with Mur's first order absorbing boundary condition and Kraus's analytical results. Eventually, Mur's first order condition have many errors at the edge and corner. On the other hand, in comparison with Kraus's analytical results, it is quite good agreement, and the validity of present method is confirmed.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.12 s.105
• /
• pp.1378-1388
• /
• 2005

Spatial control of sound is essential to deliver better sound to the listener's position in space. As it can be experienced in many listening environments. the quality of sound can not be manifested over every Position in a hall. This motivates us to control sound in a region we select. The primary focus of the developed method has to do with the brightness and contrast of acoustic image in space. In particular, the acoustic brightness control seeks a way to increase loudness of sound over a chosen area, and the contrast control aims to enhance loudness difference between two neighboring regions. This enables us to make two different kinds of zone - the zone of quiet and the zone of loud sound - at the same time. The other perspective of this study is on the direction of sound. It is shown that we can control the direction of perceived sound source by focusing acoustic energy in wavenumber domain. To begin with, the proposed approaches are formulated for pure-tone case. Then the control methods are extended to a more general case, where the excitation signal has broadband spectrum. In order to control the broadband signal in time domain, an inverse filter design problem is defined and solved in frequency domain. Numerical and experimental results obtained in various conditions certainly validate that the acoustic brightness, acoustic contrast, direction of wave front can be manipulated for some finite region in space and time.

• Proceedings of the IEEK Conference
• /
• 2002.07a
• /
• pp.137-140
• /
• 2002

In this paper, a new algorithm for two-dimensional (2-D) finite-difference time-domain (FDTD) method is presented. By this new method, the maximum time step size can be increased over the Courant-Friedrich-Levy (CFL) condition restraint. This new algorithm is adapted from an Alternating-Direction Implicit FDTD (ADI-FDTD) method. However, unlike the ADI-FDTD algorithm. the alternation is performed with respect to the blocks of fields rather than with respect to each respective coordinate direction. Moreover. this method can be efficiently simulated with parallel computation. and it is more efficient than the conventional FDTD method in terms of CPU time. Numerical formulations are shown and simulation results are presented to demonstrate the effectiveness and efficiency of our proposed method.

• Economic and Environmental Geology
• /
• v.33 no.1
• /
• pp.41-50
• /
• 2000

Geo-radar survey which has the advantage of high-resolution and relatively fast survey has been widely used for engineering and environmental problems. Three-dimensional effects have to be considered in the interpretation of geo-radar for high-resolution. However, there exists a trouble on the analysis of the three dimensional effects. To solve this problem an efficient three dimension numerical modeling algorithm is needed. Numerical radar modeling in three dimensional case requires large memory and long calculating time. In this paper, a finite difference method time domain solution to Maxwell's equations for simulating electromagnetic wave propagation in three dimensional media was developed to make economic algorithm which requires smaller memory and shorter calculating time. And in using boundary condition Liao absorption boundary. The numerical result of cross-hole radar survey for tunnel is compared with real data. The two results are well matched. To prove application to three dimensional analysis, the results with variation of tunnel's incident angle to survey cross-section and the result when the tunnel is parallel to the cross-section were examined. This algorithm is useful in various geo-radar survey and can give basic data to develop dat processing and inversion program.

• Journal of the Institute of Electronics Engineers of Korea TC
• /
• v.40 no.5
• /
• pp.181-186
• /
• 2003

This paper presents a numerical dispersion relation for the two-dimensional finite-difference time-domain method based on the alternating-direction implicit time-marching scheme(2-D ADI-FDTD), which method has the potential to considerably reduce tile number of time iterations especially in case where the fine spatial lattice relative to the wavelength is used to resolve fine geometrical features. The proposed analytical relation for 2-D ADI-FDTD is compared with those relations in the Previous works. Through numerical tests, the dispersion equation of this work was shown as correct one for 2-D ADI-FDTD.