• The Journal of the Acoustical Society of Korea
• /
• v.15 no.2
• /
• pp.72-78
• /
• 1996

This thesis deals with a method to solve a two-and-one-half-dimensional ($2\frac12$ D) problem, which means that the ocean environment is two-dimensional whereas the source is fully three-dimensionally propagating, including three-dimensional refraction phenomena and three-dimensional back-scattering, using two-dimensional two-way parabolic equation method combined with Fourier synthesis. Two dimensional two-way parabolic equation method uses Galerkin's method for depth and Crank-Nicolson method and alternating direction for range and provides a solution available to range-dependent problem with wave-field back-scattered from discontinuous interface. Since wavenumber, k, is the function of depth and vertical or horizontal range, we can reduce a dimension of three-dimensional Helmholtz equation by Fourier transforming in the range direction. Thus transformed two-dimensional Helmholtz equation is solved through two-way parabolic equation method. Finally, we can have the $2\frac12$ D solution by inverse Fourier transformation of the spectral solution gained from in the last step. Numerical simulation has been carried out for a canonical ocean environment with stair-step bottom in order to test its accuracy using the present analysis. With this spectral parabolic equation method, we have examined three-dimensional acoustic propagation properties in a specified site in the Korean Straits.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.6
• /
• pp.100-105
• /
• 2013

This paper presents a sea surface reconstruction method that uses measured radar images by applying filtering techniques and identifying wave characteristics of the surrounding the Ieodo ocean research station using WaveFinder (X-band wave measurement radar), which is installed in the station. In addition, the results obtained from real radar images are used to verify the reconstructed sea surface. WaveFinder is a marine system that was developed to measure wave information in real time. The WaveFinder installed in the station could acquire sequent images for the sea surface at constant time intervals to obtain real time information (Wave height, mean wave period, wave directionality, etc.) for the wave by getting a three-dimensional spectrum by applying an FFT algorithm to the acquired sequent images and wave dispersion relation. In particular, we found the wave height using the SNR (Signal to noise ratio) of the acquired images. The wave information measured by WaveFinder could be verified by comparing and analyzing the results measured using the wave measurement instrument (Sea level monitor) in the station. Additionally, the wave field around the station could be reconstructed through the three-dimensional spectrum and the inverse FFT filtering from the analyzed results for the measured radar images. We verified the applicability of the sea surface reconstruction method by comparing the measured and simulated sea surfaces.

• 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 Institute of Electronics Engineers of Korea TC
• /
• v.39 no.4
• /
• pp.201-208
• /
• 2002

In this paper, we present a new formulation using time-domain electric field integral equation (TD-EFIE) to obtain transient scattering response from arbitrarily shaped three-dimensional conducting bodies. The time derivative of the magnetic vector potential is approximated with a central finite difference and the scalar potential is time averaged by dividing it into two terms. This approach with an implicit method using central difference results in accurate and more stable transient scattering responses from conducting objects. Detailed mathematical steps are included and several numerical results are presented and compared with the inverse discrete Fourier transform (IDFT) of the frequency-domain solution.

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

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.4
• /
• pp.427-435
• /
• 2008

In this paper, a stable marching-on in time(MOT) method with a time domain combined field integral equation(CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time domain electric field integral equation(EFIE) with the magnetic field integral equation(MFIE). The time derivatives in the EFIE and MFIE are approximated using a central finite difference scheme and other terms are averaged over time. This time domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. Numerical results with the proposed MOT scheme are presented and compared with those obtained from the conventional method and the inverse discrete Fourier transform(IDFT) of the frequency domain CFIE solution.