• Journal of Ocean Engineering and Technology
• /
• v.36 no.1
• /
• pp.11-20
• /
• 2022

Generally, tsunamis are generated by the rapid crustal movements of the ocean floor. Other factors of tsunami generation include landslides on coastal and ocean floor slopes, glacier collapses, and meteorite collisions. In this study, two numerical analyses were conducted to examine the formation, propagation, and deformation properties of landslide tsunamis. First, LS-DYNA was adopted to simulate the formation and propagation processes of tsunamis generated by dropping rigid bodies. The generated tsunamis had smaller wave heights and wider waveforms during their propagation, and their waveforms and flow velocities resembled those of theoretical solitary waves after a certain distance. Second, after the formation of the landslide tsunami, a tsunami based on the solitary wave approximation theory was generated in a numerical wave tank (NWT) with a computational domain that considered the stability/steady phase. The comparison of two numerical analysis results over a certain distance indicated that the waveform and flow velocity were approximately equal, and the maximum wave pressures acting on the upright wall also exhibited similar distributions. Therefore, an effective numerical model such as LS-DYNA was necessary to analyze the formation and initial deformations of the landslide tsunami, while an NWT with the wave generation method based on the solitary wave approximation theory was sufficient above a certain distance.

• Nuclear Engineering and Technology
• /
• v.55 no.8
• /
• pp.3088-3101
• /
• 2023

A pressurized vessel-pipe-safety valve (PVPSV) combination is a commonly used configuration in nuclear power plants, and a good numerical model is essential for the system design, sizing and performance optimization. However, owing to the large-scale and cross-scale features, it is still a challenge to build a system level numerical model with both high accuracy and efficiency. To overcome this, a novel system level modeling method which can synthesize the advantages of various models is proposed in this paper. For system modeling, the analytical approach, the method of characteristics (MOC) and the surrogate model approach are respectively adopted to predict the dynamics of the pressure vessel, the connecting pipe and the safety valve, and different models are connected through data interfaces. With this system model, dynamic simulations were carried out and both the stable and the unstable system responses were obtained. For the model verification purpose, the simulation results were compared with those obtained from experiments and full CFD simulations. A good agreement and a better efficiency were obtained, verifying the ability of the model and the feasibility of the modeling method proposed in this paper.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.86-86
• /
• 2020

수리구조물에 관한 기존 연구들은 대부분 기능성과 안정성 측면에서 본체에 작용하는 유체력에 대한 안정성에 주안점을 두고 있다. 수리구조물 상·하류의 수위차에 기인한 기초지반내의 흐름 및 간극수압 변화는 하천 구조물의 안정성을 연구하는데 매우 중요하다. 해양에서는 파랑하중에 의한 과잉간극수압이 액상화를 발생시켜 해안구조물의 안정에 큰 영향을 미치는 것으로 보고되며, 이에 관련 연구들이 활발하게 진행되고 있다. 반면, 하천구조물 주변 지반의 흐름 및 간극수압 뿐 아니라, 액상화에 관한 연구는 아직 미진한 실정이다. 본 연구에서는 수리구조물 주변의 유동 및 와동 현상 뿐 아니라, 수위차에 따른 지반 내부 유동장과 간극수압에 관한 특성을 분석하기 위해 유체-구조물-지반 비선형 상호작용을 고려할 수 있는 수치수조를 새롭게 제안하였다. 그리고 제안하는 수치수조의 타당성 및 유효성을 검증하기 위해 기존 실험값과 비교·검토를 수행하였고, 그 결과는 거의 유사한 경향을 나타내었다. 또한 이 수치수조에 다양한 입사조건(상·하류 수위차)에 적용하여 유체-구조물-지반의 비선형동적상호간섭 해석을 수행하였다. 최종적으로 수치수조에서 측정한 구조물 주변의 유동, 와동, 수위로부터 수리특성을 논의하였다. 게다가 지반내의 흐름과 간극수압을 측정하여 상·하류 수위차가 수리구조물의 안정성에 미치는 영향을 분석할 수 있었다.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.3
• /
• pp.111-121
• /
• 2023

This paper intends to introduce the applicability of HydroQus to a problem of a tanker collision against a semi-submersible type floating offshore wind turbine (FOWT). HydroQus is a plug-in based on potential flow theory that generates interactive hydroforces in a commercial Finite element analysis (FEA) code Abaqus/Explicit. Frequency response analyses were conducted for a 10MW capacity FOWT to obtain hydrostatic and hydrodynamic constants. The tanker was modeled with rigid elements, while elastic-plastic elements were used for the FOWT. Mooring chains were modeled to implement station keeping ability of the FOWT. Two types of fracture models were considered: constant failure strain model and combined failure strain model HC-LN model composed of Hosford-Coulomb (HC) model & localized necking (LN) model. The damage extents were evaluated by hydroforces and failure strain models. The largest equivalent plastic strain observed in the cases where both restoring force and radiation force were considered. Stress triaxiality and damage indicator analysis showed that the application of HC-LN model was suitable. It could be stated that applications of suitable failure strain model and hydrodynamics into the collision simulations were of importance.

• Nuclear Engineering and Technology
• /
• v.55 no.7
• /
• pp.2454-2465
• /
• 2023

Film boiling may lead to burnout of the heating element. Even though burnout does not occur, the heating element is subject to deformation because it is not sufficiently strong to withstand external loads. In particular, the ballooning and rupture of a tube under film boiling are important phenomena in the field of nuclear reactor safety. If the tube-type cladding of nuclear fuel ruptures owing to high internal pressure and thermal load, radioactive materials inside the cladding are released to the coolant. Therefore, predicting the ballooning and rupture is important. This study presents numerical simulations to predict the ballooning behavior and rupture time of a horizontal tube at high internal pressure under saturated film boiling. To do so, a multi-step coupled simulation of conjugated film boiling heat transfer and ballooning using creep model is adopted. The numerical methods and models are validated against experimental values. Two different nonuniform heat flux distributions and four different internal pressures are considered. The three-step simulation is enough to obtain a convergent result. However, the single-step simulation also successfully predicts the rupture time. This is because the film boiling heat transfer characteristics are slightly affected by the tube geometry related to creep ballooning.

• Wind and Structures
• /
• v.37 no.6
• /
• pp.461-471
• /
• 2023

The main objective of this study was to establish design guidelines for three key design variables (spar thickness, spar diameter, and total draft) by examining their impact on the stress distribution and resonant frequency of a 2.5-MW spar-type floating offshore wind turbine substructure under extreme marine conditions, such as during Typhoon Bolaven. The current findings revealed that the substructure experienced maximum stress at wave frequencies of either 0.199 Hz or 0.294 Hz, consistent with previously reported experimental findings. These results indicated that the novel simulation method proposed in this study, which simultaneously combines hydrodynamic diffraction analysis, computational dynamics analysis, and structural analysis, was successfully validated. It also demonstrated that our proposed simulation method precisely quantified the stress distribution of the substructure. The novel findings, which reveal that the maximum stress of the substructure increases with an increase in total draft and a decrease in spar thickness and spar diameter, offer valuable insights for optimizing the design of spar-type floating offshore wind turbine substructures operating in various harsh marine environments.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.3
• /
• pp.11-22
• /
• 2010

An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

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

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.391-401
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and the results are compared to the existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results show that as cross-sectional area decreases with the reduction in downstream pressure, flow rate increases and reaches the maximum when the speed index (mean velocity divided by wave speed) is near the unity at the point of minimum cross-section area, indicating the flow limitation or choking (flow speed equals wave speed) in one-dimensional studies. for further reductions in downstream pressure, flow rate decreases. The flow limitation or choking consist of the main reasons of waterfall effect which occurs in the airways, capillaries of lung, and other veins. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is 2, the area throat is located near the downstream end. As this ratio is increased to 3, the constriction moves to the upstream end of the tube.

• Journal of Korea Water Resources Association
• /
• v.49 no.1
• /
• pp.51-60
• /
• 2016

The purpose of this study is to examine the control function of discharged flow due to the shape and plane arrangement of permeable submerged breakwater. For the discussion on it in detail, 3-dimensional numerical model based on PBM (Porous Body Model), which is able to simulate directly interaction of Fluid Permeable structure Seabed has been used to simulate water discharge in a NWT (Numerical Water Tank). To verify the applicability, LES-WASS-3D is analyzed comparing to the experimental result about propagation characteristics of dam-break wave through a permeable structure. Using the results obtained from numerical simulation, the effects of the shape and plane arrangement of submerged breakwater on reducing velocity and flow induction have been discussed related to the mean flow distribution and vertical distributions of horizontal velocities around ones.