• Journal of the Korean Society of Hazard Mitigation
• /
• v.6 no.2 s.21
• /
• pp.61-68
• /
• 2006

The earthquake parameters are known to be responsible for determination of the initial free surface profiles of tsunamis. This study investigates the effects of earthquake parameters to variation of initial free surface profiles of tsunamis triggered by an impulsive undersea earthquake. The target event is the 1983 Central East Sea Tsunami recorded as the most devastating tsunami in Korea during last several decades. Among the earthquake parameters, the strike angle may play a most significant role in determining the initial free surface.

• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.28-31
• /
• 2011

An attempt was made to compute the free surface deformation due to the impact of a water droplet. The Cauchy Poisson, i.e. the initial value problem, was solved with the kinematic and dynamic free surface boundary conditions linearized. The zero order Hankel transformation and Laplace transform were applied to the related equations. The initial condition for the free surface profile was derived from a captured video image. The effect of the surface tension was not significant with the water mass used in this investigation. The computed and observed free surface deformations were compared.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.9 no.3
• /
• pp.135-142
• /
• 2009

The initial free surface profile and energy of a tsunami are closely related to the wave heights of a tsunami in nearshore and can be determined by using the earthquake parameters. Along the Western coast of Japan, the possibility of a tsunami triggering by undersea earthquakes is very high. Many seismologists have attempted to predict the parameters of earthquakes that could occur in these regions, but it is difficult to accurately predict them. As such, several case studies have been conducted involving behaviors of an unexpected tsunami that occurred in this region. If a relationship between the earthquake parameters and the wave heights of a tsunami is found, it would be easier to examine the effects of the tsunami. In this study, several virtual tsunami events have been simulated, and the wave heights of the tsunami are computed by varying the earthquake parameters to examine the relationship between the earthquake parameters and the tsunami wave heights. Numerical simulations have been conducted in virtual topography.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.112-119
• /
• 2002

In this paper, an Extended Finite Element Method is proposed by adding discontinuity and singularity enrichment functions to the standard FEM approximation. In this method, the singularity and the discontinuity of the crack are efficiently modeled by using initial regular mesh without refining mesh near the crack tip, so that it enables express the asymptotic stress field near crack tip and crack surface successfully. The developed method was verified by evaluating crack tip stress profile and stress intensity factor of mode Ⅰ/mode Ⅱ fracture problems and the results showed the effectiveness and robustness for fracture problem.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.3
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.4
• /
• pp.403-413
• /
• 1985

Hot-wire measurements of second and third-order mean products of velocity fluctuations have been made in the separated, reattached, and redeveloping boundary layer behind a vertical fence. Mean velocity, wall static pressure distributions have also been measured in the whole flow field. Upstream of the reattachment point, the separated shear layer developes as a free mixing layer, but the gradient of the maximum slope thickness, turbulent intensities and the Reynolds shear stress are higher than that of the mixing layer due to initial streamline curvature and the effects of highly turbulent recirculating flow region. In the reattachment region, Reynolds shear stress and triple products near the surface is far more rapid than the decrease of the shear stress; that is the presence of the solid wall has a marked effect on the apparent gradient diffusivity of intensity or shear stress and throws doubts upon the usefulness of the simple gradient diffusivity model in this region.

• Journal of Korea Water Resources Association
• /
• v.54 no.11
• /
• pp.943-953
• /
• 2021

The propagation of impact wave induced by landslide and debris flow occurred on the slope of lake, reservoir and bays is a three-dimensional natural phenomenon associated with strong interaction of debris flow and water flow in complex geometrical environments. We carried out 3D numerical modeling of such impact wave in a bay using a multiphase turbulence flow model and a rheology model for non-Newtonian debris flow. Numerical results are compared with previous experimental result to evaluate the performance of present numerical approach. The results underscore that the reasonable predictions of both thickness and speed of debris flow head penetrating below the water surface are crucial to accurately reproduce the maximum peak height and free surface profiles of impact wave. Two predictions computed using different initial debris flow thicknesses become different from the instant when the peaks of impact waves fall due to the gravity. Numerical modeling using relatively thick initial debris flow thickness appears to well reproduce the water surface profile of impact wave propagating across the bay as well as wave run-up on the opposite slope. The results show that the maximum run-up height on the opposite slope is not sensitive to the initial thickness of debris flows of same total volume. Meanwhile, appropriate rheology model for debris flow consisting of inviscid particle only should be employed to more accurately reproduce the debris flow propagating along the channel bottom.

• Journal of The Korean Astronomical Society
• /
• v.36 no.3
• /
• pp.89-95
• /
• 2003

Observations of dark matter dominated dwarf and low surface brightness disk galaxies favor density profiles with a flat-density core, while cold dark matter (CDM) N-body simulations form halos with central cusps, instead. This apparent discrepancy has motivated a re-examination of the microscopic nature of the dark matter in order to explain the observed halo profiles, including the suggestion that CDM has a non-gravitational self-interaction. We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. The SIDM particles scatter each other elastically, which results in an effective thermal conductivity that heats the halo core and flattens its density profile. These similarity solutions are relevant to galactic and cluster halo formation in the CDM model. We assume that the local density maximum which serves as the progenitor of the halo has an initial mass profile ${\delta}M / M {\propto} M^{-{\epsilon}$, as in the familiar secondary infall model. If $\epsilon$ = 1/6, SIDM halos will evolve self-similarly, with a cold, supersonic infall which is terminated by a strong accretion shock. Different solutions arise for different values of the dimensionless collisionality parameter, $Q {\equiv}{\sigma}p_br_s$, where $\sigma$ is the SIDM particle scattering cross section per unit mass, $p_b$ is the cosmic mean density, and $r_s$ is the shock radius. For all these solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of $r_s$. We find two different regimes for these solutions: 1) for $Q < Q_{th}({\simeq} 7.35{\times} 10^{-4}$), the core density decreases and core size increases as Q increases; 2) for $Q > Q_{th}$, the core density increases and core size decreases as Q increases. Our similarity solutions are in good agreement with previous results of N-body simulation of SIDM halos, which correspond to the low-Q regime, for which SIDM halo profiles match the observed galactic rotation curves if $Q {\~} [8.4 {\times}10^{-4} - 4.9 {\times} 10^{-2}]Q_{th}$, or ${\sigma}{\~} [0.56 - 5.6] cm^2g{-1}$. These similarity solutions also show that, as $Q {\to}{\infty}$, the central density acquires a singular profile, in agreement with some earlier simulation results which approximated the effects of SIDM collisionality by considering an ordinary fluid without conductivity, i.e. the limit of mean free path ${\lambda}_{mfp}{\to} 0$. The intermediate regime where $Q {\~} [18.6 - 231]Q_{th}$ or ${\sigma}{\~} [1.2{\times}10^4 - 2.7{\times}10^4] cm^2g{-1}$, for which we find flat-density cores comparable to those of the low-Q solutions preferred to make SIDM halos match halo observations, has not previously been identified. Further study of this regime is warranted.