• Geomechanics and Engineering
• /
• v.19 no.1
• /
• pp.79-91
• /
• 2019

Renewed interest in the long-term pile foundations has been driven by the increase in offshore wind turbine installation to generate renewable energy. A monopile subjected to repetitive loads experiences an evolution of displacements, pile rotation, and stress redistribution along the embedded portion of the pile. However, it is not fully understood how the embedded pile interacts with the surrounding soil elements based on different pile geometries. This study investigates the long-term soil response around offshore monopiles using finite element method. The semi-empirical numerical approach is adopted to account for the fundamental features of volumetric strain (terminal void ratio) and shear strain (shakedown and ratcheting), the strain accumulation rate, and stress obliquity. The model is tested with different strain boundary conditions and stress obliquity by relaxing four model parameters. The parametric study includes pile diameter, embedded length, and moment arm distance from the surface. Numerical results indicate that different pile geometries produce a distinct evolution of lateral displacement and stress. In particular, the repetitive lateral load increases the global lateral load resistance. Further analysis provides insight into the propagation of the shear localization from the pile tip to the ground surface.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.685-692
• /
• 2003

A numerical simulation of flame propagation in a micro combustor was carried out. Combustor has a sub -millimeter depth cylindrical internal volume and axisymmetric one-dimensional was used to simplify the geometry. Semi-empirical heat transfer model was used to account for the heat loss to the walls during the flame propagation. A detailed chemical kinetics model of $H_2/Air$ with 10 species and 16 reaction steps was used to calculate the combustion. An operator-splitting PISO scheme that is non-iterative, time-dependent, and implicit was used to solve the system of transport equations. The computation was validated for adiabatic flame propagation and showed good agreement with existing results of adiabatic flame propagation. A full simulation including the heat loss model was carried out and results were compared with measurements made at corresponding test conditions. The heat loss that adds its significance at smaller value of combust or height obviously affected the flame propagation speed as final temperature of the burnt gas inside the combustor. Also, the distribution of gas properties such as temperature and species concentration showed wide variation inside the combustor, which affected the evaluation of total work available of the gases.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.2
• /
• pp.14-21
• /
• 2015

Design lifetime of a wind turbine is required to be at least 20 years. The most important step to ensure the deign is to evaluate the loads on the wind turbine as accurately as possible. In this study, extreme design load of a offshore wind turbine using Garrad Hassan (GH) Bladed and National Renewable Energy Laboratory (NREL) FAST codes are calculated considering structural dynamic loads. These wind turbine aeroelastic analysis codes are high efficiency for the rapid numerical analysis scheme. But, these codes are mainly based on the mathematical and semi-empirical theories such as unsteady blade element momentum (UBEM) theory, generalized dynamic wake (GDW), dynamic inflow model, dynamic stall model, and tower influence model. Thus, advanced CFD-dynamic coupling method is also applied to conduct cross verification with FAST and GH Bladed codes. If the unsteady characteristics of wind condition are strong, such as extreme design wind condition, it is possible to occur the error in analysis results. The NREL 5 MW offshore wind turbine model as a benchmark case is practically considered for the comparison of calculated designed loads. Computational analyses for typical design load conditions such as normal turbulence model (NTM), normal wind profile (NWP), extreme operation gust (EOG), and extreme direction change (EDC) have been conducted and those results are quantitatively compared with each other. It is importantly shown that there are somewhat differences as maximum amount of 18% among numerical tools depending on the design load cases.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.