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

• Journal of the Korean Society for Nondestructive Testing
• /
• v.20 no.2
• /
• pp.138-149
• /
• 2000

Various modeling techniques for ultrasonic wave propagation and scattering problems in finite solid media are presented. Elastodynamic boundary value problems in inhomogeneous multi-layered plate-like structures are set up for modal analysis of guided wave propagation and numerically solved to obtain dispersion curves which show propagation characteristics of guided waves. As a powerful modeling tool to overcome such numerical difficulties in wave scattering problems as the geometrical complexity and mode conversion, the Boundary Element Method(BEM) is introduced and is combined with the normal mode expansion technique to develop the hybrid BEM, an efficient technique for modeling multi mode conversion of guided wave scattering problems. Time dependent wave forms are obtained through the inverse Fourier transformation of the numerical solutions in the frequency domain. 3D BEM program development is underway to model more practical ultrasonic wave signals. Some encouraging numerical results have recently been obtained in comparison with the analytical solutions for wave propagation in a bar subjected to time harmonic longitudinal excitation. It is expected that the presented modeling techniques for elastic wave propagation and scattering can be applied to establish quantitative nondestructive evaluation techniques in various ways.

• Geophysics and Geophysical Exploration
• /
• v.20 no.2
• /
• pp.72-77
• /
• 2017

We calculated 3D frequency- and Laplace-domain wavefields using time-domain modeling and Fourier transform or Laplace transform. We adopted OpenACC and GPU for an efficient parallel calculation. The OpenACC makes it easy to use GPU accelerators by adding directives in conventional C, C++, and Fortran programming languages. Accordingly, one doesn't have to learn new GPGPU programming languages such as CUDA or OpenCL to use GPU. An OpenACC program allocates GPU memory, transfers data between the host CPU and GPU devices and performs GPU operations automatically or following user-defined directives. We compared performance of 3D wave propagation modeling programs using OpenACC and GPU to that using single-core CPU through numerical tests. Results using a homogeneous model and the SEG/EAGE salt model show that the OpenACC programs are approximately 53 and 30 times faster than those using single-core CPU.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.8
• /
• pp.793-799
• /
• 2015

In this research, we study the propagation of longitudinal and transverse waves through a metal rod including a liquid layer using computational and experimental analyses. The propagation characteristics of longitudinal and transverse waves obtained by the computational and experimental analyses were consistent with the wave propagation theory for both cases, that is, the homogeneous metal rod and the metal rod including a liquid layer. The fluid-structure interaction modeling technique developed for the computational wave propagation analysis in this research can be applied to the more complex structures including solid-liquid interfaces.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.257-263
• /
• 2019

Multilayered structures include lamination by relatively thick plies and thin interlayers. For efficient peridynamic analysis of dynamic fracturing multilayered structures, the interlayer is modeled using ghost peridynamic particles while the ply is formulated via real peridynamics. With the nonlocal ghost interlayer, one may keep the discretization resolution low for the ply. In this study, the characteristics of dynamic wave propagation through the nonlocal ghost interlayer in peridynamic analysis are investigated. It is observed that the interlayer not only binds adjacent plies, but also significantly influences energy transfer between plies, and thereby their deformation and motion. In addition, near a surface or boundary, peridynamic particles do not have a full nonlocal neighborhoods. This causes the effective material properties near the surface to be different from those in the bulk. Surface correction based on neighborhood volumes is employed. The impact of surface correction on wave propagation in multilayered structures is investigated.

• Proceedings of the Acoustical Society of Korea Conference
• /
• 1996.06a
• /
• pp.33-36
• /
• 1996

해양에서 내부파 (internal wave)는 수직적으로 밀도가 크게 다른 두 층의 경계면에서 생성되는데, 일반적으로 경계면은 수온약층 (thermocline)에 해당된다. 본 연구에서는 간단한 단주기(하나의 주기만을 갖는)내부파를 가정하고 이 내부파에 의해서 야기될 수 있는 저주파 (100-400 HZ) 전파손실의 변동 정도를 모델링을 통해 추정하였다. 내부파의 파장은 1Km이고 진폭은 25m이며 수심 50-100m 사이에 존재하는 것으로 가정하였다. 전파손실 계산은 PE (parabolic equation) 기법을 도입한 모델을 이용하였다.수층의 두께는 2000m로 균일하고 퇴적층 및 기반암층의 두께는 각각 500m로 가정하였다. 모델링 결과 단주기 내부파의 1 파장이 진행하는 동안 고정된 수신 깊이별로 20km 거리에서 10 dB 이상의 큰 변동이 야기 될 수 있음을 확인하였다. 이는 내부파 내에서 큰 굴절을 변화로 인한 음파의 산란 때문으로 볼 수 있다. 실제 해양에서는 내부파가 단주기 파동 형태로 나타나기 보다는 여러 주기의 파동이 중첩되어 나타나므로 정확한 내부파 영향을 추정하기 위해서는 보다 실제에 가까운 내부파를 모델링할 필요가 있다.

• Geophysics and Geophysical Exploration
• /
• v.21 no.4
• /
• pp.213-219
• /
• 2018

In this study, we simulated 3D wave propagation modeling using a Xeon Phi x200 processor and compared the parallel computation performance with that using a Xeon CPU. Unlike the 1st generation Xeon Phi coprocessor codenamed Knights Corner, the 2nd generation x200 Xeon Phi processor requires no additional communication between the internal memory and the main memory since it can run an operating system directly. The Xeon Phi x200 processor can run large-scale computation independently, with the large main memory and the high-bandwidth memory. For comparison of parallel computation, we performed the modeling using the MPI (Message Passing Interface) and OpenMP (Open Multi-Processing) libraries. Numerical examples using the SEG/EAGE salt model demonstrated that we can achieve 2.69 to 3.24 times faster modeling performance using the Xeon Phi with a large number of computational cores and high-bandwidth memory compared to that using the 12-core CPU.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.145-152
• /
• 2016

We used graphics processing units for an efficient time-domain 3D wave propagation modeling. Since graphics processing units are designed for massively parallel processes, we need to optimize the calculation and memory management to fully exploit graphics processing units. We focused on the memory management and examined the performance of programs with respect to the memory management methods. We also tested the effects of memory transfer on the performance of the program by varying the order of finite difference equation and the size of velocity models. The results show that the memory transfer takes a larger portion of the running time than that of the finite difference calculation in programs transferring whole 3D wavefield.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.131-135
• /
• 2016

Coarray is a parallel processing technique introduced in the Fortran 2008 standard. Coarray can implement parallel processing using simple syntax. In this research, we examined applicability of Coarray to seismic parallel processing by comparing performance of seismic data processing programs using Coarray and MPI. We compared calculation time using seismic wave propagation modeling and one to one communication time using domain decomposition technique. We also compared performance of parallel reverse-time migration programs using Coarray and MPI. Test results show that the computing speed of Coarray method is similar to that of MPI. On the other hand, MPI has superior communication speed to that of Coarray.