• Journal of the Korean Society for Marine Environment & Energy
• /
• v.20 no.1
• /
• pp.37-44
• /
• 2017

Complexity of multiphase flows due to existence of more than two interface including free-surface in one system, cannot be simulated easily. Since more than two fluids affect to flows and disturb interface, non-linearities such as instabilities can be appeared. Among several instabilities on multiphase flows, one of representative is Rayleigh-taylor instability. In order to examine in importance of density disparity, several cases with numerous Atwood number are set. Moreover, investigation of influence on initial disturbance were also considered. Moving particle simulation (MPS) method, which was employed in this paper, was not widely used for multiphase problem. In this study, by adding new particle interaction models such as self-buoyance correction, surface tension, and boundary condition at interface models, MPS were developed having more strength of physics and robust. By applying newly developed multiphase MPS, considered cases are performed and compared each other. Additionally, though existence of disagreement of magnitude of rising velocity between theoretical values from linear potential theory and that of numerical simulation, agreement of tendency can be proved of similarity of result. the discordance of magnitude can be explained due to non-linear effects on numerical simulation which was not considered in theoretical result.

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

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.591-596
• /
• 2006

This study is base on the fault detection and diagnosis when a crack is happened a structure. A crack in a structure will affect the modal parameters. We are searched a percentage of changes in the natural frequencies according to changes of location and propagation of the crack using the Rayleigh's energy method. This method is presented to identify the presence of a crack and its location. The study is carried out both theoretically and experimentally and the results are presented in this paper. The location of the crack is also moved from the fixed end to the free end along its length. The changes in natural frequencies are observed from theoretically study, due to the presence of the crack at different locations and depths, and the percentage change in frequency values are calculated. These results are confirmed by the experiments. And then, a difference between a cracked beam and uncracked beam observed using the bispectrum as high-order spectrum.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.73-85
• /
• 2006

Array observations of the vertical component of microtremors are frequently conducted to estimate a subsurface layered-earth structure on the assumption that microtremors consist predominantly of the fundamental mode Rayleigh waves. As a useful tool in the data collection, processing and analysis, the spatial autocorrelation (SPAC) method is widely used, which in practice requires a circle array consisting of M circumferential stations and one centre station (called "M-station circle array", where M is the number of stations). The present paper considers the minimum number of stations required for a circle array for efficient data collection in terms of analytical efficacy and field effort. This study first rearranges the theoretical background of the SPAC algorithm, in which the SPAC coefficient for a circle array with M infinite is solely expressed as the Bessel function, $J_0(rk)$ (r is the radius and k the wavenumber). Secondly, the SPAC coefficient including error terms independent of the microtremor energy field for an M-station circle array is analytically derived within a constraint for the wave direction across the array, and is numerically evaluated in respect of these error terms. The main results of the evaluation are: 1) that the 3-station circle array when compared with other 4-, 5-, and 9-station arrays is the most efficient and favourable for observation of microtremors if the SPAC coefficients are used up to a frequency at which the coefficient takes the first minimum value, and 2) that the Nyquist wavenumber is the most influential factor that determines the upper limit of the frequency range up to which the valid SPAC coefficient can be estimated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.347-352
• /
• 2007

Eigenvalue reduction schemes approximate the lower eigenmodes that represent the global behavior of the structures. In the previous study, we proposed a two-level condensation scheme (TLCS) for the construction of a reduced system. In the first step, the selection of candidate elements by energy estimation, Rayleigh quotient, through Ritz vector calculation. In the second step, the primary degrees of freedom are selected by the sequential elimination method from the degrees of freedom connected to the candidate elements in the first step. In the present study, we propose TLCS combined with iterative improved reduced system (IIRS) to increase accuracy of the higher modes in the intermediate range. Also, it is possible to control the accuracy of the eigenvalues and eigenmodes of the reduced system. Finally, numerical examples demonstrate the performance of the proposed method.