• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.37-43
• /
• 2006

Dynamic numerical analysis of geotechnical problems requires a way to simulate the decrease of energy as the domain of interest gets larger. This phenomenon is usually referred to as radiation damping or geometric attenuation and it is distinguished from material damping in which elastic energy is actually dissipated by viscous, hysteretic, or other mechanism. The fact that the domain of analysis in numerical modeling must be chosen, however, causes a need for special attention at the boundary. This observation leads directly to the idea of determining the dynamic response of the interior region from a finite model consisting of the interior region subjected to a boundary condition which ensures that all energy arriving at the boundary is absorbed. This paper presents a simple methodology to simulate transmitting boundaries condition using viscoelastic infinite elements within the recently developed "OpenSees" finite element code. The methodology used here provides that the level of absorption for traveling waves is efficient enough for practical purposes, but unsatisfactory for the case of sharp incident angles. The effectiveness of the infinite elements for the absorption of incident waves at boundaries is evaluated via example analysis.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.116-126
• /
• 2008

In gas hydrate survey, seismic amplitude and frequency characteristics play a very important role in determining whether gas hydrate exists. According to the variation of source frequency and scatterer size, we study seismic amplitude characteristics using elastic modeling applied at staggered grids. Generally speaking, scattering occurs in proportion to the square of source frequency and the scatterer volume, which has an effect on seismic amplitude. The higher source frequency is, the more scattering occurs in gas hydrate bearing zone. Therefore, BSR is hardly observed in high frequencies. On the other side, amplitude blanking zone and BSR is clearly observed in lower frequencies although the resolution is poor as a whole. Seismic reflections traveling through free-gas layer below gas hydrate bearing zone decay so severely a high frequency component that a low frequency term is dominant. Amplitude anomaly of BSR result from high acoustic impedance contrast due to free-gas, which is a very crucial factor to estimate gas hydrate bearing zone. Seismic frequency analysis is carried out using wavelet transform method that frequency component could be decomposed with time variation. In application of wavelet transform to the seismic physical experiments data, we can observe that reflections traveling through air layer, which corresponds to the free-gas layer, decay a high frequency component.

• Economic and Environmental Geology
• /
• v.43 no.6
• /
• pp.603-613
• /
• 2010

The Domi Basin in the South Sea of Korea is located between the Jeju Basin and Ulleung Basins, and is characterized by several sediment sags that are interested to have formed by crustal extension. This paper aims to derive an optimized seismic data processing procedure which helps stratigraphic interpretation of the Domi Basin. In particular, our data processing flow incorporated horizon velocity analysis (HVA) and surface-relative wave equation multiple rejection (SRWEMR) to improve the quality of stack section by enhancing the continuity of reflection events and suppressing peg-leg multiples respectively. As a result of processing procedures in this study, unconformities were recognized in the stack section that defines the early and middle Miocene, Eocene-Oligocene sequences. In addition, the overall quality of the stack section was increased as essential data to investigate the evolution of the basin. The suppression of multiple resulted in the identification of the Cretaceous basement. The data processing scheme evaluated through this study is expected to improve the standardization of processing sequences for seismic data from the Domi and adjacent Sora and north-Sora Basins.

• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.89-97
• /
• 2021

Fault detection in seismic data is well suited to the application of machine learning algorithms. Accordingly, various machine learning techniques are being developed. In recent studies, machine learning models, which utilize synthetic data, are the particular focus when training with deep learning. The use of synthetic training data has many advantages; Securing massive data for training becomes easy and generating exact fault labels is possible with the help of synthetic training data. To interpret real data with the model trained by synthetic data, the synthetic data used for training should be geologically realistic. In this study, we introduce a method to generate realistic synthetic seismic data. Initially, reflectivity models are generated to include realistic fault structures, and then, a one-way wave equation is applied to efficiently generate seismic stack sections. Next, a migration algorithm is used to remove diffraction artifacts and random noise is added to mimic actual field data. A convolutional neural network model based on the U-Net structure is used to verify the generated synthetic data set. From the results of the experiment, we confirm that realistic synthetic data effectively creates a deep learning model that can be applied to field data.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.521-532
• /
• 2007

Elastic wave propagations in the semi-infinite transversely isotropic cylinder under various kinds of longitudinal impact loads are analyzed using the axisymmetric finite element method and Houbolt time-integration scheme. For which the finite element program is newly constructed and verified through the comparison of present numerical results with those by other researchers. E-type glass-epoxy composite cylinders with different fiber volume fractions are adopted and studied in detail with dynamic responses of the isotropic cylinder. Three dimensional wave motions are given in graphic form to show the realistic view of the wave propagation. Nondimensionalized dynamic characteristic variables which relate the size of finite element mesh, the time step, and the wave speed are presented for obtaining accurate and stable numerical results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.629-640
• /
• 2007

In this paper an analytical study is carried out to improve the capacity of absorbing boundary using dashpot, one of the most widely used absorbing boundaries in FEM. Using 2-D harmonic plane wave equation, absorbing boundary condition is modified to maximize its capacity according to the incident angle. Validity of the absorbing boundary conditions which is modified is investigated by adopting the solution of Miller and Pursey. The Miller and Pursey's problem is then numerically simulated using the finite element method. The absorption ratios are calculated by comparing the displacements at the absorbing boundary to those at the free field without the absorbing boundary. The numerical study is carried out through comparison of displacement at the interior region and the boundary of the numerical model.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.15 no.3
• /
• pp.159-166
• /
• 2003

In this study, we develop a finite element model to directly solve the Laplace equation while keeping the same computational efficiency as the models based on the extended mild-slope equation which has been widely used for calculation of wave transformation in shallow water. For this, the computational domain is discretized into finite elements with a single layer in the vertical direction. The velocity potential in the element is then expressed in terms of the potentials at the nodes located at water surface, and the Galerkin method is used to construct the numerical model. A common shape function is adopted in horizontal direction, and the cosine hyperbolic function in vertical direction, which describes the vertical behavior of progressive waves. The model was developed for vertical two-dimensional problems. In order to verify the developed model, it is applied to vertical two-dimensional problems of wave reflection and transmission. It is shown that the present finite element model is comparable to the models based on extended mild-slope equations in both computational efficiency and accuracy.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.2
• /
• pp.32-39
• /
• 1999

In this paper, the characteristics of the ray tracing method (RTM) based on the cylindrical wave are discussed for the high frequency vibration analysis of two-dimensional structures. A ray tube describing the emanating cylindrical wave is used to derive the governing equation for incident reflected, and transmitted ray tubes which satisfies the condition at the coupled boundary. The suggested ray model is applied to panel array structures, and the predicted results for 2-panel, 3-panel, and 4-panel array structures are compared to those by Statistical energy analysis (SEA) and Wave intensity analysis(WIA). More enhanced prediction was obtained compared to the SEA, and similar prediction performance was observed to the WIA. Additionally, the RTM has a novel feature that it can estimate the spatially smoothed distribution of vibration energy and vibration intensity. It is expected that the present RTM can be used as one of the useful tools for the high frequency vibration analysis of two-dimensional coupled structures.

• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.169-177
• /
• 2002

We develop the weighted-averaging finite-element method which uses four kinds of element sets. By constructing global stiffness and mass matrices for four kinds of element sets and then averaging them with weighting coefficients, we obtain a new global stiffness and mass matrix. With the optimal weighting coefficients minimizing grid dispersion and grid anisotropy, we can reduce the number of grid points required per wavelength to 4 for a $1\%$ upper limit of error. We confirm the accuracy of our weighted-averaging finite-element method through accuracy analyses for a homogeneous and a horizontal-layer model. By synthetic data example, we reconfirm that our method is more efficient for simulating a geological model than previous finite-element methods.