• Geomechanics and Engineering
• /
• v.15 no.6
• /
• pp.1161-1171
• /
• 2018

The interaction between particles and fluid was investigated by IB-SEM numerical method which is a combination of combing the spectral/hp element method and the rigid immersed boundary method. The accuracy of this numerical method was verified based on the computed results with the traditional body-fitted mesh in numerical simulation of the flow through the cylinder. Then the governing equations of particles motion and contact in fluid are constructed. The movement of the particles and the interaction between the fluid and the particles are investigated. This method avoided the problem of low computational efficiency and error caused by the re-division of the grid when the solids moved. Finally, the movement simulation of multi particles in the fluid was carried out, which can provide a completely new numerical simulation method.

• Wind and Structures
• /
• v.9 no.3
• /
• pp.231-243
• /
• 2006

An improvement to the spectral representation algorithm for the simulation of wind velocity fields on large scale structures is proposed in this paper. The method proposed by Deodatis (1996) serves as the basis of the improved algorithm. Firstly, an interpolation approximation is introduced to simplify the computation of the lower triangular matrix with the Cholesky decomposition of the cross-spectral density (CSD) matrix, since each element of the triangular matrix varies continuously with the wind spectra frequency. Fast Fourier Transform (FFT) technique is used to further enhance the efficiency of computation. Secondly, as an alternative spectral representation, the vectors of the triangular matrix in the Deodatis formula are replaced using an appropriate number of eigenvectors with the spectral decomposition of the CSD matrix. Lastly, a turbulent wind velocity field through a vertical plane on a long-span bridge (span-wise) is simulated to illustrate the proposed schemes. It is noted that the proposed schemes require less computer memory and are more efficiently simulated than that obtained using the existing traditional method. Furthermore, the reliability of the interpolation approximation in the simulation of wind velocity field is confirmed.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.137-139
• /
• 2007

To enhance capability on discerning local and regional seismic phases, such as, Pn, Pg, Sn, Rg, etc, within the crust, 2-D numerical forward modeling will be applied to the data obtained from local seismic stations by simulating almost all waves including not only body wave but also surface wave generated without having to explicitly include them under consideration of Q factor. In this study, after getting rid of instrumental response by deconvolution, pseudo-spectral method instead of relying on typical numerical methods, such as, FEM(Finite Element Method) and FDM(Finite Difference Method), will be implemented for 2-D numerical forward modeling by considering velocities of P-wave and S-wave, density, and Q factors. Ultimately, the Power of reaching the enhanced capability on discerning local and regional seismic phases will make it easier for us to identify the seismic source, whether it is originated from man-made explosion or pure earthquake.

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

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.587-614
• /
• 2014

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.

• Structural Engineering and Mechanics
• /
• v.14 no.5
• /
• pp.505-520
• /
• 2002

Based on the modal coordinates of the structure, a finite-element and CQC (complete quadratic combination) method for analyzing the coupled buffeting response of long-span bridges is presented. The formulation of nodal equivalent aerodynamic buffeting forces is derived based on a reasonable assumption. The power spectral density and variance of nodal displacements and elemental internal forces of the bridge structure are computed using the finite-element method and the random vibration theory. The method presented is very efficient and can consider the arbitrary spectrum and spatial coherence of natural winds and the multimode and intermode effects on the buffeting responses of bridge structures. A coupled buffeting analysis of the Jiangyin Yangtse River Suspension Bridge with 1385 in main span is performed as an example. The results analyzed show that the multimode and intermode effects on the buffeting response of the bridge deck are quite remarkable.

• Journal of Power System Engineering
• /
• v.15 no.6
• /
• pp.27-34
• /
• 2011

The increasing capabilities of the computers enable us to utilize various numerical schemes for the time-domain simulations concerned with 3-dimensional free-surface wave problems. There are still difficulties to solve such kind of problems, however. That's because long time simulations with large computational domain are needed in time-domain analysis. So, we need faster and more efficient numerical schemes to get the solutions practically for these problems. In this paper, a high-order spectral/boundary-element method is used for the numerical investigation of physics involved in wave-body interaction. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. To get the robust study in these topics, various numerical tests are performed and compared with others' works.

• Structural Engineering and Mechanics
• /
• v.55 no.5
• /
• pp.1055-1086
• /
• 2015

The dynamic response of two-dimensional unbounded domain on the rigid bedrock in the time domain is numerically obtained. It is realized by the modified scaled boundary finite element method (SBFEM) in which the original scaling center is replaced by a scaling line. The formulation bases on expanding dynamic stiffness by using the continued fraction approach. The solution converges rapidly over the whole time range along with the order of the continued fraction increases. In addition, the method is suitable for large scale systems. The numerical method is employed which is a combination of the time domain SBFEM for far field and the finite element method used for near field. By using the continued fraction solution and introducing auxiliary variables, the equation of motion of unbounded domain is built. Applying the spectral shifting technique, the virtual modes of motion equation are eliminated. Standard procedure in structural dynamic is directly applicable for time domain problem. Since the coefficient matrixes of equation are banded and symmetric, the equation can be solved efficiently by using the direct time domain integration method. Numerical examples demonstrate the increased robustness, accuracy and superiority of the proposed method. The suitability of proposed method for time domain simulations of complex systems is also demonstrated.

• Journal of the Optical Society of Korea
• /
• v.15 no.3
• /
• pp.216-221
• /
• 2011

A diffraction grating is a highly symmetric optical element with a physical structure that is invariant under translational spatial movements. The translational symmetry is reflected in the fields that are diffracted from the grating. Here, we introduce a plane-parallel mirror pair onto the grating, which translates the fields through double reflections, and we describe a method of exploiting the symmetry to enhance the spectral resolution of a diffraction grating beyond the limit that is set by the number of grooves. The mirror pair creates another virtual grating beside the original one, effectively doubling the number of grooves. Addition of more mirror pairs can further increase the effective number of grooves despite the increased complexity and difficulty of experimental implementation. We experimentally demonstrate the spectral linewidth reduction by a factor of four in a neon fluorescence spectrum. Even though the geometrical restriction on the mirror deployment limits our method to a certain range of the whole spectrum, as a practical application example, a bulky spectrometer that is nearly empty inside can be made compact without sacrificing the resolution.

• Journal of Ocean Engineering and Technology
• /
• v.8 no.1
• /
• pp.23-32
• /
• 1994

A numerical procedure is described for simultaneously predicting the motion and structural responses of tension leg platforms (TLPs) in multi-directional irregular waves. The developed numerical approach is based on a combination of a three dimensional source distribution method, the finite element method for structurally treating the space frame elements and a spectral analysis technique of directional waves. The spectral description for the linear responses of a structure in the frequency domain is sufficient to completely define the responses. This is because both the wave inputs and the responses are stationary Gaussian ran dom process of which the statistical properties in the amplitude domain are well known. The hydrodynamic interactions among TLP members, such as columns and pontoons, are included in the motion and structural analysis. The effect of wave directionality has been pointed out on the first order motion, tether forces and structural responses of a TLP in multi-directional irregular waves.