• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.27 no.11
• /
• pp.119-125
• /
• 1990

This paper describes fiber optic voltage sensor using EOM-BSO (Electro-Optic Modulator-Bismuth Silicon Oxcide). Transceiver has an electical/optical converter and an optical/electrical converter which consist of light emitting diode, PIN-PD, and electronic circuits. Multimode fiber cable of $100/140{\mu}m$ core/clad diameter is used for connecting the transceiver to fiber cable and fiber optic voltage sensor. Before our experiments, by applying the Maxwell equations and wave equations, We derive matrix equation on wave propagation in the BSO single crystal. And also we derive optimal equation on intensity modulation arising through an analyzer. According to experi-mental results, fiber optic voltage sensor has maximum $2.5{\%}$ error within the applied AC voltage of 800V. As the applied voltage increases, saturation values of voltage sensor also increase. This phenomenon is caused by optical rotatory power of BSO single crystal. And temperature dependence of sensitivity for fiber optical rotatory power of BSO single crystal. And temperature dependence of sensitivity for fiber optic voltage sensor in the temperature range from$-20^{\circ}C\to\60^{\circ}C$ are measured within ${\pm}0.6{\%}$. And frequency characteristics of the voltage sensor has good frequency characteristics from DC to 100kHz.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.220-226
• /
• 2011

In the elastic wave equations, both horizontal and vertical displacements are defined. Since we can measure both the horizontal and vertical displacements in field acquisition, these displacements compose a displacement vector. In this study, we propose a frequency-domain elastic waveform inversion technique taking advantage of the magnitudes of displacement vectors to define objective function. When we apply this displacement-vector objective function to the frequency-domain waveform inversion, the inversion process naturally incorporates the back-propagation algorithm. Through the inversion examples with the Marmousi model and the SEG/EAGE salt model, we could note that the RMS error of the solution obtained by our algorithm decreased more stably than that of the conventional method. Particularly, the density of the Marmousi model and the low-velocity sub-salt zone of the SEG/EAGE salt model were successfully recovered. Since the gradient direction obtained from the proposed objective function is numerically unstable, we need additional study to stabilize the gradient direction. In order to perform the waveform inversion using the displacementvector objective function, it is necessary to acquire multi-component data. Hence, more rigorous study should be continued for the multi-component land acquisition or OBC (Ocean Bottom Cable) multi-component survey.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Journal of Navigation and Port Research
• /
• v.29 no.10 s.106
• /
• pp.853-858
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation have been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use this method for calculations in the case of random waves because long and stable computation (e.g for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand if the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, can be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes' 5th order wave propagation and a random wave propagation.

• Geophysics and Geophysical Exploration
• /
• v.18 no.2
• /
• pp.74-84
• /
• 2015

To yield physically meaningful images through elastic reverse-time migration, the wavefield separation which extracts P- and S-waves from reconstructed vector wavefields by using elastic wave equation is prerequisite. For expanding the application of the elastic reverse-time migration to anisotropic media, not only the anisotropic modelling algorithm but also the anisotropic wavefield separation is essential. The anisotropic wavefield separation which uses pseudo-derivative filters determined according to vertical velocities and anisotropic parameters of elastic media differs from the Helmholtz decomposition which is conventionally used for the isotropic wavefield separation. Since applying these pseudo-derivative filter consumes high computational costs, we have developed the efficient anisotropic wavefield separation algorithm which has capability of parallel computing by using GPUs (Graphic Processing Units). In addition, the highly efficient anisotropic elastic reverse-time migration algorithm using MPI (Message-Passing Interface) and incorporating the developed anisotropic wavefield separation algorithm with GPUs has been developed. To verify the efficiency and the validity of the developed anisotropic elastic reverse-time migration algorithm, a VTI elastic model based on Marmousi-II was built. A synthetic multicomponent seismic data set was created using this VTI elastic model. The computational speed of migration was dramatically enhanced by using GPUs and MPI and the accuracy of image was also improved because of the adoption of the anisotropic wavefield separation.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2005.10a
• /
• pp.153-159
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation rave been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use VOF method for calculations in the case of random waves because long and stable computation ( e.g. for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand of the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, am be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes'5th order wave propagation and a random wave propagation.

• 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 Computational Structural Engineering Institute of Korea
• /
• v.28 no.2
• /
• pp.123-130
• /
• 2015

This paper presents a new formulation for transient simulations of microwave propagation in heterogeneous unbounded domains. In particular, perfectly-matched-layers(PMLs) are introduced to allow for wave absorption at artificial boundaries used to truncate the infinite extent of the physical domains. The development of the electromagnetic PML targets the application to engineering mechanics problems such as structural health monitoring and inverse medium problems. To formulate the PML for plane electromagnetic waves, a complex coordinate transformation is introduced to Maxwell's equations in the frequency-domain. Then the PML-endowed partial differential equations(PDEs) for transient electromagnetic waves are recovered by the application of the inverse Fourier transform to the frequency-domain equations. A mixed finite element method is employed to solve the time-domain PDEs for electric and magnetic fields in the PML-truncated domain. Numerical results are presented for plane microwaves propagating through concrete structures, and the accuracy of solutions is investigated by a series of error analyses.

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