• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.115-119
• /
• 2009

We numerically investigated propagation of various waves in the two-phase flows such as sound wave, shock wave, rarefaction wave, and contact discontinuity in terms of pressure, void fraction, velocity and density of the two phases. The waves have been generated by a hydrodynamic shock tube, a pair of symmetric impulsive expansion, impulsive pressure and impulsive void waves. The six compressible two-fluid two-phase conservation laws with interfacial friction terms have been solved in two fractional steps. The first PDE Operator is solved by the HLL scheme and the second Source Operator by the semi-implicit stiff ODE solver. In the HLL scheme, the fastest wave speeds were estimated by the analytic eigenvalues of an approximate Jacobian matrix. We have discussed how the interfacial friction terms affect the wave structures in the numerical solution.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.8
• /
• pp.1931-1939
• /
• 1993

The wave propagation characteristics of laminated composites subjected to a transverse high-velocity impact of a steel ball is investigated. For this purpose, high-velocity impact experiments were conducted to obtain the strain response histories, and a finite element analysis based on the higher-order shear deformation theory in conjunction with the static contact law is used. Test materials for investigation are glass/epoxy laminated composite materials with $[0^{\circ}/45^{\circ}/0^{\circ}/-45^{\circ}]_{2s}$ and $[90^{\circ}/-45^{\circ}/90^{\circ}-45^{\circ}/90^{\circ}]_{2s}$ stacking sequences. As a result, the strain responses obtained from the experiments represented the wave propagation characteristics in the transversely impact, also the wave propagation velocities obtained from high-velocity impact experiments and wave propagation theory agree well.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.311-314
• /
• 2007

The experimental method of measuring dynamic properties of structures was presented. The method is based on the flexural wave propagation characteristics. Using the method, change in structural dynamic properties due to damage is measured. The crack has much more significant impact on the strain energy than the inertial effects. From this, the sensitivity of the dynamic stiffness on the crack location is estimated by calculating the strain energy. When the wave propagates, the strain and kinetic energies shows cyclic changed over space. The crack that occurred at locations where the wave energy is in the form of the potential energy affected most significantly the wave propagation characteristics. The effects of crack location on the wave propagation were used to determine the crack location.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.21 no.3
• /
• pp.248-257
• /
• 2011

It was reported recently that railpads can be included as a continuous elastic support of the rail and the model was justified from experiments. In general, however, railpads are installed discretely on sleepers with a regular span. The effect of the discrete railpad was not clearly examined so far in such a high frequency range. In this paper, the effect of the railpads in track modelling for high frequencies is investigated by means of the finite element analysis. To do that, the railpads are regarded as 'a continuous elastic support' and 'a discrete elastic support' in this paper. The dispersion relations and decaying features are predicted and compared between the two models up to 80 kHz.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.242-251
• /
• 2022

Automatic differentiation automatically calculates the derivatives of a function using the chain rule once the forward operation of a function is defined. Given the recent development of computing libraries that support automatic differentiation, many researchers have adopted automatic differentiation techniques to solve geophysical inverse problems. We analyzed the advantages, disadvantages, and performances of automatic differentiation techniques using the gradient calculations of seismic full waveform inversion objective functions. The gradients of objective functions can be expressed as multiplications of the derivatives of the model parameters, wavefields, and objective functions using the chain rule. Using numerical examples, we demonstrated the speed of analytic differentiation and the convenience of complex gradient calculations for automatic differentiation. We calculated derivatives of model parameters and objective functions using automatic differentiation and derivatives of wavefields using analytic differentiation.

• Geophysics and Geophysical Exploration
• /
• v.19 no.3
• /
• pp.131-135
• /
• 2016

Coarray is a parallel processing technique introduced in the Fortran 2008 standard. Coarray can implement parallel processing using simple syntax. In this research, we examined applicability of Coarray to seismic parallel processing by comparing performance of seismic data processing programs using Coarray and MPI. We compared calculation time using seismic wave propagation modeling and one to one communication time using domain decomposition technique. We also compared performance of parallel reverse-time migration programs using Coarray and MPI. Test results show that the computing speed of Coarray method is similar to that of MPI. On the other hand, MPI has superior communication speed to that of Coarray.

• Geophysics and Geophysical Exploration
• /
• v.21 no.4
• /
• pp.213-219
• /
• 2018

In this study, we simulated 3D wave propagation modeling using a Xeon Phi x200 processor and compared the parallel computation performance with that using a Xeon CPU. Unlike the 1st generation Xeon Phi coprocessor codenamed Knights Corner, the 2nd generation x200 Xeon Phi processor requires no additional communication between the internal memory and the main memory since it can run an operating system directly. The Xeon Phi x200 processor can run large-scale computation independently, with the large main memory and the high-bandwidth memory. For comparison of parallel computation, we performed the modeling using the MPI (Message Passing Interface) and OpenMP (Open Multi-Processing) libraries. Numerical examples using the SEG/EAGE salt model demonstrated that we can achieve 2.69 to 3.24 times faster modeling performance using the Xeon Phi with a large number of computational cores and high-bandwidth memory compared to that using the 12-core CPU.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.10a
• /
• pp.267-273
• /
• 2004

This paper targets to define controlling factors of pressure-coupled combustion response and estimate their effects on droplet evaporation process. Dynamic characteristics of hydrocarbon propellant vaporization perturbed by acoustic pressure are numerically simulated and analyzed. 1-D droplet model including phase equilibrium between two phases is applied and acoustic wave is expressed by harmonic function. Effects of various design factors and acoustic pressure on combustion response are investigated with parametric studies. Results show that driving frequency of acoustic perturbation and ambient pressure have important roles in determining magnitude and phase of combustion response. On the other hand, other parameters such as gas temperature, initial droplet size and temperature, and amplitude of acoustic wave cause only minor changes to magnitude of combustion response. Resultant changes in phase of heat of vaporization and thermal wave in droplet highly influence magnitude and phase of combustion response.

• The Journal of the Acoustical Society of Korea
• /
• v.30 no.7
• /
• pp.361-367
• /
• 2011

Propagating waves (flexural, longitudinal and shear waves) travelling with constant amplitudes and evanescent waves decaying exponentially are generated on a cylindrical shell. Evanescent waves are generally generated in the vicinity of an vibration excitation point and near ends of the shell. But the evanescent waves can generates strong axial vibration at the ends of the cylindrical shell. The strong end axial vibration due to those evanescent waves has been observed in an author's previous paper dealing with measurements of the in-plane axial vibration of a finite cylindrical shell. In this paper the strong end axial vibration due to the evanescent waves has been theoretically analyzed. In order to analyze the vibration of the cylindrical shell, wave propagation approach has been implemented. Comparison between theoretical and experimental results for the axial vibration of the shell showed that the strong evanescent wave can be generated due to mode conversion (conversion from flexural wave to evanescent wave) at the ends of cylindrical shell. It also showed that the evanescent wave can generate the strong axial vibration near the ends of the cylindrical shell and that it can have effect even on 1/3 of the total length of the shell.