• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.95-103
• /
• 2007

An analytical closed-form solution for wave propagation velocity and damping in saturated porous media is presented in this paper The fully coupled field model with compressible solid Brains and pore water were used to derive this solution. An engineering approach for the analysis of fully saturated porous media was adopted and closed-form solutions for one dimensional wave propagation in a homogeneous domain were derived. The solution is highly versatile in that it considers compression of the solid grains, compression of the pore water, deformation of the porous skeleton, and spatial damping and can be used to compute wavespeeds of first and second kind and damping coefficients in various geologic materials. This solution provides a means of analyzing the influence of material property variations on wavespeed and attenuation. In Part 2 of this work the theoretical solution is incorporated into the numerical code and the code is used in a parametric study on wave propagation velocity and damping.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.315-322
• /
• 2010

The phase-screen and the split-step Fourier migrations, which are implemented in both the frequency-wavenumber and frequency-space domains by using one-way scalar wave equation, allow imaging in laterally heterogeneous media with less computing time and efficiency. The generalized-screen migration employs the series expansion of the exponential, unlike the phase-screen and the split-step Fourier migrations which assume the vertical propagation in frequency-wavenumber domain. In addition, since the generalized-screen migration generalizes the series expansion of the vertical slowness, it can utilize higher-order terms of that series expansion. As a result, the generalized-screen migration has higher accuracy in computing the propagation with wide angles than the phase-screen and split-step Fourier migrations for media with large and rapid lateral velocity variations. In this study, we developed a 2D prestack generalized-screen migration module for imaging a complex subsurface efficiently, which includes various dips and large lateral variations. We compared the generalized-screen propagator with the phase-screen propagator for a constant perturbation model and the SEG/EAGE salt dome model. The generalized-screen propagator was more accurate than the phase-screen propagator in computing the propagation with wide angles. Furthermore, the more the higher-order terms were added for the generalized-screen propagator, the more the accuracy was increased. Finally, we compared the results of the generalizedscreen migration with those of the phase-screen migration for a model which included various dips and large lateral velocity variations and the synthetic data of the SEG/EAGE salt dome model. In the generalized-screen migration section, reflectors were positioned more accurately than in the phase-screen migration section.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.183-188
• /
• 2005

This study targets to investigate thermal wave transfer in a droplet strongly coupled with acoustic pressure and its effects on combustion response. The one-dimensional vaporization model uses SRK equation of state and flash calculation method to obtain more accurate thermophysical properties and compute vapor-liquid equilibrium. Calculations of an n-pentane droplet exposed into a perturbing nitrogen gas is carried out in terms of different ambient gas pressures and wave frequencies. The thermal wave is transferred more effectively at lower frequencies, which results in the decrease in the amplitude of the response.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.70-75
• /
• 2000

In the seismic migration, Kirchhoff and reverse time migration are used in general. In the reverse time migration using wave equation, two-way and one-way wave equation are applied. The approach of one-way wave equation uses approximately computed downward continuation extrapolator, it need tess amounts of calculations and core memory in compared to that of two-way wave equation. In this paper, we applied one-way wave equation to pre-stack reverse time migration. In the frequency-space domain, forward propagation of source wavefield and back propagration of measured wavefield were executed by using monochromatic one-way wave equation, and zero-lag cross correlation of two wavefield resulted in the image of subsurface. We had implemented prestack migration on a massively parallel processors (MPP) CRAYT3E, and knew the algorithm studied here is efficiently applied to the prestck migration due to its suitability for parallelization.

• The Korean Journal of Petroleum Geology
• /
• v.1 no.1 s.1
• /
• pp.47-52
• /
• 1993

A method and results of computations are presented for the 2-D seismic migration process in the frequency-wavenumber domain for the laterally and vertically inhomogeneous medium. In order to take the intrinsic attenuation effect into account in the migration process the complex-valued wave velocity is used in the wavefield extrapolation operator, improving the generalized frequency-wavenumber migration technique. The imaginary part of the complex-valued wave velocity includes the seismic quality factor Q value. In derivation of the solution of the wave equation for the medium of inhomogeneous wave velocity and anelasticity, the inhomogeneous medium is mathematically converted to an equivalent system which consists of a homogeneous medium of averaged slowness and an inhomogeneous distribution of hypothetical wave source. The strength of the hypothetical wave source depends on the deviation of squared slowness from the averaged value of the medium. Results of numerical computation using the technique show more distinct geologic images than those using the convensional generalized frequency-wavenumber migration. Especially, the obscured images due to the wave attenuation by anelasticity are restored to show sharp boundaries of structures. The method will be useful in the imaging of the reflection data obtained in the regions of possible petroleum or natural gas reservoir and of fractured zone.

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

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

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.6
• /
• pp.512-519
• /
• 2013

It is important to understand the vibrating behavior of plate structures for its many engineering applications. In this study, the vibration characteristics of strip plates that have finite width and infinite length are investigated theoretically and numerically. The waveguide finite element(WFE) approach, which is an effective tool for studying waveguide structures, is used in this study. The WFE method requires only a cross-sectional finite element model, and uses theoretical harmonic solutions to assess wave propagation along the longitudinal direction. First, WFE results for a simple strip plate are compared with the theoretical results(i.e., dispersion diagrams and point mobilities) to validate the numerical model. Then, in the numerical analysis, different numbers of longitudinal stiffeners are included in the plate model to investigate the effects of stiffeners in terms of the dispersion curves and mobilities. Finally, the dispersion curves of a stiffened double plate are obtained to examine the characteristics of its wave propagation.

• Journal of the Korean Society for Railway
• /
• v.17 no.1
• /
• pp.7-13
• /
• 2014

Stiffness variations and wave propagation/reflection in railway catenaries are the primary sources of contact loss between a pantograph and a railway contact wire. This paper analyzes which design parameter is more important for 200km/h conventional rail and 300km/h high-speed rail, in order to effectively reduce the contact loss. For the high-speed rail, the wave propagation and reflection in the overhead contact lines are more influential than the stiffness variation over a span. When the high-speed rail needs to speed-up, it is necessary to develop higher strength contact wires in order to increase the wave propagation speed. In addition, the dropper clamp mass should be reduced in order to alleviate the wave reflection. However, it is noted that the increase in the tension to a messenger wire could deteriorate the current collection quality, which contrasts with expectations. For the 200km/h conventional rail, the stiffness variation over a span is more influential than the wave propagation and reflection. Therefore, shortening span length, increasing the tension in the contact wire and optimizing the location of the droppers are recommended for a smoother stiffness variation over the span.

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