• Geophysics and Geophysical Exploration
• /
• v.7 no.4
• /
• pp.234-244
• /
• 2004

In order to calculate 3-dimensional wavefield using finite-element method in frequency domain, we must factor so huge sparse impedance matrix. Because of difficulties of handling of this huge impedance matrix, 3-dimensional wave equation modeling is conducted mainly in time domain. In this study, we simulate the 3-D wavefield using finite-element method in Laplace domain by combining high-performance parallel finite-element solver and SWEET (Suppressed Wave Equation Estimation of Traveltime) algorithm which can calculate the traveltime and the amplitude. To verify this combination, we applied it to the SEG/EAGE 3D salt model in serial and parallel computing environments.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.197-202
• /
• 1995

비정상(unsteady) 압축성(compressible) 유동에 의한 공력음향(aeroacoustics)을 모사하여 공력소음원을 해석하기 위해서는 고차(high order)의 정확도와 높은 해상도(resolution)를 가지며, 상대적으로 계산시간을 많이 필요로 하지 않는 외재적(explicit) 유한차분법이 필수적으로 요구된다. 이것은 주어진 차분방식과 격자계로써 공간과 시간상에 존재하는 미소크기의 파동성분들을 충분히 구현하여야 만족할 만한 수치해를 얻을 수 있기 때문이다. 본 연구에서는, 그러한 유한차분법 중 최근에 관심의 대상이 되고있는 삼각(tridiagonal)또는 오각(pentadiagonal) 집적유한차분법(compact finite difference scheme)이 최대의 해상도를 갖도록 하는 수학적인 방법을 개발하고, 이 방법으로써 새롭게 집적유한차분법을 최적화하였다. 개발된 최적화 방법은, 푸리에 해석법(Fourier analysis)을 통하여 파동수(wavenumber) 영역에서 수학적으로 계산된 위상오차(phase error)를 최소화하는 것이며, 이러한 개념과 방법은 본 연구에서 처음으로 집적유한차분법에 적용되었다. 여러가지 절단정확도(truncation order)에 대해서 최적화 된 집적유한차분법들이 실제 공간과 시간상에서 보여주는 정확도와 오차특성을 알아보기 위하여, 이 방법들을 1차원 선형파동방정식에 적용하였고, 이 결과를 통하여 가장 정확하고 효과적인 절단정확도의 집적유한차분법을 선별하였다. 특히, 오각(pentadiagonal)법에 비해 더욱 효율적인 6차 삼각(tridiagonal)법을 1차원 Euler방정식에 적용하여, 비선형 파동에 대한 모사를 수행할 수 있었다.

• Korean Journal of Optics and Photonics
• /
• v.6 no.2
• /
• pp.156-164
• /
• 1995

Based on the scalar wave equation of optical fibers, the dispersion characteristics of arbitrarily profiled fibers were analyzed. We used the I-D FEM employing quadratic interpolation fucntions to solve the scalar wave equation. We simulated the DC optical fibers as objects, and searched for the refractive index distribution to minimize the total dispersion. In DC fibers, we found the design parameters for which the total dispersion was almost zero at $\lambda=1.3{\mu}m and 1.55{\mu}m$ simultaneously. We also found the design parameters where the dispersion was flattened, less than 1.0 ps/km.nm for$\lambda=1.4~1.7{\mu}m$1. and the dispersion was as low as 0.65 ps/km.nm at $\lambda=1.55{\mu}m$..

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

• Journal of the Korea Convergence Society
• /
• v.9 no.6
• /
• pp.175-182
• /
• 2018

In this study, the theoretical method and the finite element analysis were designed in parallel to fabricate basic research data on the production of tool horn for cutting machine with ultrasonic vibration energy. In order to perform high-performance ultrasonic cutting, it is necessary to vibrate only with longitudinal vibration instead of transverse vibration. In order to efficiently transmit the mechanical vibration energy, the maximum amplitude should be generated at the output portion. Therefore, the tool horn must be designed so that the excitation frequency of the oscillator and the natural frequency of the tool horn are the same. In order to design the resonance of the tool horn, there are a theoretical approach using the one-dimensional wave equation and a method of reflecting the finite element analysis result to the design model. In this study, the approximate dimensions of the tool horn are first determined through the one- Based on the results of the finite element analysis, the optimal model was selected and reflected in the final shape of the tool horn. We will use this information as the basic data of actual tool horn for cutting, and will compare the production and experimental data with the contents of this research.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.1
• /
• pp.1-8
• /
• 2005

For the simulation of tsunami propagation an active dispersion-correction two-dimensional finite element model has been developed based on a shallow-water wave equation. This model employs an arbitrary triangular mesh and an explicit time integration scheme. However, the physical dispersion effects as included in the Boussinesq equations can be taken into account in the computation. The validity of the dispersion-correction scheme developed in this study is verified through the comparison of numerical solutions calculated using the new scheme with analytical ones considering dispersion effect of waves. As a result, the present model is shown to be considerably accurate.

• The Korean Journal of Rheology
• /
• v.10 no.3
• /
• pp.160-164
• /
• 1998

본 연구에서 Doi-Edwards 점탄성 조성방정식의 Hadamard 안정성 분석을 행하였 다. Hadamard 안정성은 방정식의 탄성 성질과 연관되는 특성으로 파장이 짧고 진동수가 큰 파동에 의한 외란 하에서 식의 안정성을 의미한다. 먼저 안정성을 위한 일반 3차원 조건을 수립하고 단순한 1차원과 2차원 외란하에서 필요조건을 구하였다. Doi-Edwards 이론을 따 르는 물질의 단순전단유동을 고려함에 의하여 순간 전단변형률이 1.8786을 넘어설 때 파장 이 짧고 진동수가 큰 외란에 의하여 불안정성이 나타남이 증명되었다. 이 안정성의 임계치 는 실제 고분자공정 뿐 아니라 실험실에서도 쉽게 도달할수 있는 값으로 이와 같은 불안정 유동은 mi-crophase separation과 같은 물리적 현상과는 관련이 있다는 증거가 없으므로 조 성방정식 자체가 지니는 수학적 모순점에 기인한 것이라 할수 있다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.859-871
• /
• 2014

In order to analyze the density current due to salt and temperature difference, this study develops new numerical model (LES-WASS-3D ver. 2.0) by introducing state equation for salt and temperature and 3D advection-diffusion equation to existing 3D numerical wave tank (LES-WASS-3D ver. 1.0). To verify the applicability, the newly-developed numerical model is analyzed comparing to the experimental result of existing numerical model. In the result, it well implement the behavior and vertical salt concentration of advected and diffused seawater as well as flow velocity and temperature of the discharged warm water. This confirms the validity and effectiveness of the developed numerical model.