• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.295-304
• /
• 2012

This study proposed a two-dimensional horizontal numerical model based on the nonlinear shallow water wave equations to simulate tsunami propagation and coastal inundation. We numerically investigated the possible impacts of tsunami caused by the triple interlocked Tokai, Tonankai and Nankai Earthquakes on the Jeju coastal areas, using the proposed model. The simultaneous Tokai, Tonankai and Nankai Earthquakes were created a virtual tsunami model of an M9.0 earthquake. In numerical analysis, a grid nesting method for the local grid refinement in shallow coastal regions was employed to sufficiently reproduce the shoaling effects. The numerical model was carefully validated through comparisons with the data collected during the tsunami events by 2011 East Japan Earthquake and 1983 central East Sea Earthquake (Nihonkai Chubu Earthquake). Tsunami propagation triggered by the combined Tokai, Tonanakai and Nankai, Earthquakes was simulated for 10 hours to sufficiently consider the effects of tsunami in the coastal areas of Jeju Island. The numerical results revealed that water level fluctuation in tsunami propagation is greatly influenced by water-depth change, refraction, diffraction and reflection. In addition, the maximum tsunami height numerically estimated in the coastal areas of Jeju Island was about 1.6 m at Sagye port.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.150-168
• /
• 1991

A practical procedure for the ultimate compressive strength-based safety and reliability assessment of the double skin upper deck structure is described. The external compressive stress acting on the upper deck structure which is due to the still water and wave-induced sagging moment is approximately estimated by using the existing rule of classification society. The ultimate compressive stress of double skin structure under the action of sagging moment is analyzed by using idealized structural unit method. Here an idealized plate element subjected to uniaxial load is formulated by idealizing the nonlinear behaviour of the actual element taking account of the initial imperfections in the form of initial deflection and welding residual stress. The interaction effect between the local and global failure in the structure is also taken into consideration. The accuracy of the present method is verified comparing with the present solution and the existing numerical and experimental results for unit member and welded box columns. The safety of the structure is evaluated using the concept of conventional central safety factor and the reliability assessment is made by using Cornel's MVFOSM method. The present procedure is then applied to upper deck structure of double skin product oil carrier. The influence of the initial imperfections and the yield stress of the material on the safety and reliability of the structure is investigated.

• Sleep Medicine and Psychophysiology
• /
• v.14 no.1
• /
• pp.42-48
• /
• 2007

Introduction: Detrended fluctuation analysis (DFA) is used as a way of studying nonlinearity of EEG. In this study, DFA is applied on sleep EEG of normal subjects to look into its nonlinearity in terms of EEG channels and sleep stages. Method: Twelve healthy young subjects (age：$23.8{\pm}2.5$ years old, male：female=7：5) have undergone nocturnal polysomnography (nPSG). EEG from nPSG was classified in terms of its channels and sleep stages and was analyzed by DFA. Scaling exponents (SEs) yielded by DFA were compared using linear mixed model analysis. Results: Scaling exponents (SEs) of sleep EEG were distributed around 1 showing long term temporal correlation and self-similarity. SE of C3 channel was bigger than that of O1 channel. As sleep stage progressed from stage 1 to slow wave sleep, SE increased accordingly. SE of stage REM sleep did not show significant difference when compared with that of stage 1 sleep. Conclusion: SEs of Normal sleep EEG showed nonlinear characteristic with scale-free fluctuation, long-range temporal correlation, self-similarity and self-organized criticality. SE from DFA differentiated sleep stages and EEG channels. It can be a useful tool in the research with sleep EEG.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.281-288
• /
• 2014

As a gas explosion is the most fatal accident in shipbuilding and offshore plant industries, all safety critical elements on the topside of offshore platforms should retain their integrity against blast pressure. Even though many efforts have been devoted to develop blast-resistant design methods in the offshore engineering field, there still remain several issues needed to be carefully investigated. From a procedure for calculation of explosion design pressure, impulse of a design pressure model having completely positive side only is determined by the absolute area of each obtained transient pressure response through the CFD analysis. The negative pressure phase in a general gas explosion, however, is often quite considerable unlike gaseous detonation or TNT explosion. The main objective of this study is to thoroughly examine the effect of the negative pressure phase on structural behavior. A blast wall for specific FPSO topside is selected to analyze structural response under the blast pressure. Because the blast wall is considered an essential structure for blast-resistant design. Pressure time history data were obtained by explosion simulations using FLACS, and the nonlinear transient finite element analyses were performed using LS-DYNA.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.741-749
• /
• 1999

This study is to develop a method of measuring the soil water concentration by using TDR, which is based on the relationships between the bulk soil electrical conductivity of soil and the reflected wave of TDR. The proposed monitoring method is combined with two important relationships. One is that between the bulk soil electrical conductivity and the solute concentration, which is known to be linear at a constant volumetric soil water content and the other is that between the relative bulk soil electrical conductivity and the water content at a constant concentration. Some formulas have been proposed to solve the second relationship, but a new formula and the critical water content are proposed to improve the accuracy of measurement. This proposed formula estimates the relative bulk soil electrical conductivity for water contents which is divided to two regions, linear and nonlinear, by the critical water content. As the result of the comparison with other formulas, the proposed formula is proved to be superior to other formulas and to be an available method to apply to the unsaturated transient solute transport.

• Korean Journal of Optics and Photonics
• /
• v.30 no.4
• /
• pp.154-158
• /
• 2019

We have used a mid-IR (mid-infrared) continuous-wave (cw) optical parametric oscillator (OPO), developed previously and described in Ref. 12, to build a performance-evaluation setup for a mid-IR spectrometer. The used CW OPO had a wavelength tuning range of $ 2.5-3.6{\mu}m$ using a pump laser with a wavelength of 1064 nm and a fan-out MgO-doped periodically poled lithium niobate (MgO:PPLN) nonlinear crystal in a concentric cavity design. The OPO was combined with a near-IR integrating sphere and a Fourier-transform IR optical spectrum analyzer to build a performance-evaluation setup for mid-IR spectrometers. We applied this performance-evaluation setup to evaluating a mid-IR spectrometer developed domestically, and demonstrated the capability of evaluating the performance, such as spectral resolution, signal-to-noise ratio, spectral stray light, and so on, based on this setup.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.5
• /
• pp.320-333
• /
• 2019

In this study, a physics-based 3D morphology model for the estimation of an erosion rate of nourished beach is newly proposed. As a hydrodynamic module, IHFOAM toolbox having its roots on the OpenFoam is used. On the other hand, the morphology model comprised a transport equation for suspended sediment, and Exner type equation derived from the viewpoint of sediment budget with the bed load being taken to accounted. In doing so, the incipient motion of sediment is determined based on the Shields Diagram, while the bottom suspended sediment concentration, the bed load transport rate is figured out using the bottom shearing stress directly calculated from the numerically simulated flow field rather than the conventional quadratic law and frictional coefficient. In order to verify the proposed morphology model, we numerically simulate the nonlinear shoaling, breaking over the uniform beach of 1/m slope, and its ensuing morphology change. Numerical results show that the partially skewed, and asymmetric bottom shearing stresses can be successfully simulated. It was shown that sediments suspended and eroded at the foreshore by wave breaking are gradually drifted toward a shore and accumulated in the process of up-rush, which eventually leads to the formation of swash bar. It is also worth mentioning that the breaker bar formed by the sediments dragged by the back-wash flow which commences at the pinnacle of up-rush as the back-wash flow gets weakened due to the increased depth was successfully duplicated in the numerical simulation.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.5
• /
• pp.221-228
• /
• 2008

To estimate array shape with reference sources in SONAR systems, nearfield signal modeling is required for the reference sources near a towed array. Array shape estimation method based on the nearfield signal modeling generally exploits the spatial covariance matrix of the received reference sources. Among those method, nearfield eigenvector method uses the eigenvector corresponding to the maximum eigenvalue as a steering vector of the reference source. In this paper, we propose a simplified subspace fitting method based on the nearfield signal modeling with spherical wave modeling. Furthermore, we analyze performance of the array shape estimation methods based on the nearfield signal modeling for various environments. The results of the numerical experiments indicate that the simplified subspace fitting method and the nearfield eigenvector method with single reference source shows almost similar performance. Furthermore, the simplified subspace fitting method with 2 reference sources consistently estimates the shape of the array regardless of the incident angle of the reference sources, whereas the nearfield eigenvector method cannot apply for the case of 2 reference sources.

• Atmosphere
• /
• v.23 no.3
• /
• pp.293-305
• /
• 2013

Nonhydrostatic effects on convectively forced mesoscale flows in two dimensions are numerically investigated using a nondimensional model. An elevated heating that represents convective heating due to deep cumulus convection is specified in a uniform basic flow with constant stability, and numerical experiments are performed with different values of the nonlinearity factor and nonhydrostaticity factor. The simulation result in a linear system is first compared to the analytic solution. The simulated vertical velocity field is very similar to the analytic one, confirming the high accuracy of nondimensional model's solutions. When the nonhydrostaticity factor is small, alternating regions of upward and downward motion above the heating top appear. On the other hand, when the nonhydrostaticity factor is relatively large, alternating updraft and downdraft cells appear downwind of the main updraft region. These features according to the nonhydrostaticity factor appear in both linear and nonlinear flow systems. The location of the maximum vertical velocity in the main updraft region differs depending on the degrees of nonlinearity and nonhydrostaticity. Using the Taylor-Goldstein equation in a linear, steady-state, invscid system, it is analyzed that evanescent waves exist for a given nonhydrostaticity factor. The critical wavelength of an evanescent wave is given by ${\lambda}_c=2{\pi}{\beta}$, where ${\beta}$ is the nonhydrostaticity factor. Waves whose wavelengths are smaller than the critical wavelength become evanescent. The alternating updraft and downdraft cells are formed by the superposition of evanescent waves and horizontally propagating parts of propagating waves. Simulation results show that the horizontal length of the updraft and downdraft cells is the half of the critical wavelength (${\pi}{\beta}$) in a linear flow system and larger than ${\pi}{\beta}$ in a weakly nonlinear flow system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.306-318
• /
• 2018

In order to evaluate the seismic capacity of massive vertical type breakwaters which have intensively been deployed along the coast of South Korea over the last two decades, we carry out the preliminary numerical simulation against the PoHang, GyeongJu, Hachinohe 1, Hachinohe 2, Ofunato, and artificial seismic waves based on the measured time series of ground acceleration. Numerical result shows that significant sliding can be resulted in once non-negligible portion of seismic energy is shifted toward the longer period during its propagation process toward the ground surface in a form of shear wave. It is well known that during these propagation process, shear waves due to the seismic activity would be amplified, and non-negligible portion of seismic energy be shifted toward the longer period. Among these, the shift of seismic energy toward the longer period is induced by the viscosity and internal friction intrinsic in the soil. On the other hand, the amplification of shear waves can be attributed to the fact that the shear modulus is getting smaller toward the ground surface following the descending effective stress toward the ground surface. And the weakened intensity of soil as the number of attacking shear waves are accumulated can also contribute these phenomenon (Das, 1993). In this rationale, we constitute the numerical model using the model by Hardin and Drnevich (1972) for the weakened shear modulus as shear waves go on, and shear wave equation, in the numerical integration of which $Newmark-{\beta}$ method and Modified Newton-Raphson method are evoked to take nonlinear stress-strain relationship into account. It is shown that the numerical model proposed in this study could duplicate the well known features of seismic shear waves such as that a great deal of probability mass is shifted toward the larger amplitude and longer period when shear waves propagate toward the ground surface.