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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.323-334
• /
• 2008

Present study examined the functionality of the solitary wave (tsunami) control of the two-rowed porous submerged breakwater by numerical experiments, using a numerical wave tank which is based on the Navier-Stokes equation to explain fluid fields and uses a Volume of Fluid (VOF) method to capture the free water surface. Solitary wave was generated by the internal wave source installed within the computational zone in the numerical wave tank and its wave transformations by structure were compared with those in the previous study. Comparisons with the precious numerical results showed a good agreement. Based on these results, several tow-dimensional numerical modeling investigations of the water fields, including wave transformations, reflection, transmission and energy flux, by the one- and two-rowed permeable submerged breakwater under solitary waves were performed. Even if, it is a research of the limited scope, in case of two-rowed permeable submerged breakwater with $h_0/h=0.925$ ($h_0$ is height of submerged breakwater and h is water depth), the wave height damping in range of $l/L_{eff}>0.4$($L_{eff}$ is effective distance of solitary wave) can reach nearly 60% of the incident wave height. In addition, it is found that reflection coefficient increases nearly 47% and transmission coefficient decreases nearly 18% than one-rowed one. The numerical results revealed that the tow-rowed submerged breakwater can control the incident solitary wave economically and more efficiently than the one-rowed one.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• The Journal of the Acoustical Society of Korea
• /
• v.33 no.3
• /
• pp.191-199
• /
• 2014

Wave reflection and transmission characteristics in waveguides are an important issue in many engineering applications. A combined spectral element and finite element (SE/FE) method is used to investigate the effects of local non-uniformities but limited at relatively low frequencies because the SE is formulated by using a beam theory. For higher frequency applications, a method named a combined spectral super element and finite element (SSE/FE) method was presented recently, replacing spectral elements with spectral super elements. This SSE/FE approach requires a long computing time due to the coupling of SSE and FE matrices. If a local non-uniformity has a uniform cross-section along its short length, the FE part could be further replaced by SSE, which improves performance of the combined SSE/FE method in terms of the modeling effort and computing time. In this paper SSEs are combined to investigate the characteristics of waves propagating along waveguides possessing geometric non-uniformities. Two models are regarded: a rail with a local defect and a periodically ribbed plate. In the case of the rail example, firstly, the results predicted by a combined SSE/FE method are compared with those from the combined SSEs in order to justify that the combined SSEs work properly. Then the SSEs are applied to a ribbed plate which has periodically repeated non-uniformities along its length. For the ribbed plate, the propagation characteristics are investigated in terms of the propagation constant.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.28 no.5
• /
• pp.677-697
• /
• 1995

The analysis of long- period sea level variations with tidal record data around Korea, Japan, and Russia shows that about half of the variations are due to atmospheric influences. The sea level variation by water movements is the largest in the coasts along the Tsushima Current, and becomes smaller in the distant areas. It suggests that the sea level varications are related with the Tsushima Current. The effect of sea level variations to ocean circulation has been studied with a numerical model allowing barotropic sea level fluctuations, like the result with GCM (Semtner) model by Pang et al.(1993), the present model also shows that waters basically flow along isobaths over the last China Sea after geostyophic adjustment around Taiwan. However, barotropic sea level fluctuation makes the basic circulation in the Yellow Sea, which waters flow into the central Yellow Sea and out along the west coast of the Korean Peninsula. Besides this, barotropic sea level fluctuation makes long period waves over the shelf area as the Kuroshio varies. By the waves, the basic circulation in the Yellow Sea is disturbed, so that the flow pattern of oppositely flowing into the Yellow Sea along the west roast of the Korean Peninsula appears. In the Yellow Sea circulation, it seems that northwest winds strengthen the basic circulat ion In winter, and southeast winds strengthen the disturbed circulation in summer. Another point appeared by the long period wave is that the Tsushima Current possibly originates in different areas. There have been two opposing argues on the area in which the Tsushima Current originates the southwest sea of Kyushu Island and the adjacent sea of Taiwan. Through this study, we found that both of them seem to be important areas for the origin of the Tsushima Current, and one of them is possibly strengthened by long period waves. The long period waves given by the variation of the Kuroshio Current in the adjacent sea of Taiwan propagate to the Korea Strait as forced waves. The wave continuously propagates to the last Sea through the eastern channel, but reflects in the western channel due to bottom topography. The reflected waves propagate southwestward along the last China Sea as free waves and determine the sea level variations with forced waves.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.126-131
• /
• 2001

충격파-와동 간섭은 다양한 비선형 물리 현상을 포함한다. 그 중에서도, 충격파 분리와 충격파-충격파 간섭에 대해서는 잘 알려져 있는 반면, 미끄럼 층에서 반사되는 파동이나 관통 충격파 등에 대해서는 지금까지 적절한 연구가 이루어지지 못했다. 저자들은 복잡한 충격파-와동 간섭이 시간에 따라 몇 개의 국지적인 소모델 (submodels)로 귀착됨을 발견하였다. 이는 와동에 의한 충격파의 분리, 충격파-충격파 간섭, 충격파-미끄럼층 간섭, 와동 중심의 충격파통과, 그리고 충격파-소와동 간섭이다. 이러한 5 개의 소모델은 탐구 범위 내에서 충격파-와동 간섭의 전체 구조를 구성한다.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1998.09a
• /
• pp.134-137
• /
• 1998

항내의 파랑장을 구하는 것은 항만설계상 중요하다. 요즈음에는 많은 경비와 노력을 요하는 수리모형실험대신에 다양한 수치해석이 행하여 지고있다. 그러나 회절 및 다중반사, 굴절 등 복잡한 파의 변형을 표현할 수 있는 수치해석 방법으로 여러 가지가 제안되어 있지만 아직 충분한 수준은 아니라고 판단된다. 등수심에 대해서는 이미 확립되어져 있는 Green함수법을 변수심장에 응용하기 위해서 변수심장에 대한 점원파(source waves)의 기본해(unit solution)를 구하는 것이 본 연구의 목적이다. (중략)

• Journal of Industrial Technology
• /
• v.20 no.A
• /
• pp.303-309
• /
• 2000

Seismic reflection method is applied to detect shallow location of limestone in Chechon area. The data using hammer source is compared with that of weight drop. Small size hammer and weight-drop are used as energy source and 100Hz geophones are used for data aquisition. Data processing is conducted utilizing the available processing technique of "Geobit", which is seismic data processing software developed by KIGAM. The result of above data processing, the velocity of topsoil layer is 1,250m/sec. The velocity of this area is higher than other area because loading trucks pass this area and make this layer compact. And in limestone area, hammer is proposed to energy source instead of weight drop because the energy propagates the layer very well.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.6
• /
• pp.23-28
• /
• 2011

In this paper, a new generation method for ship induced waves is proposed using the fully non-linear 3-D numerical model with non-reflected wave generation system (LES-WASS-3D). A ship induced wave generated by the newly proposed method is examined in comparison with that obtained by an empirical formula. It is then shown that there is a good agreement in free surface the elevation between them. As a result, it is revealed that a ship induced wave in a 3-D numerical wave field can be simulated well using LES-WASS-3D.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.505-510
• /
• 2002

The attenuation of waves transmitted through non-conservative joints that are shown in many paractical structures, is affected by the impedance and the orientation of the joint. In this paper, the joints between plate structures are assumed to be modeled as linear spring-dashpot systems and the transmission and reflection of vibration energy in the medium to high frequency ranges are investigated. The calculated power transmission and reflection coefficients are applied to the PFA method for the prediction of energy density and intensity in structures.