• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.214-214
• /
• 2003

During a model test of Hutton TLP, a "ringing" response was first observed about 20 years ago. This phenomenon is a resonant build up over the time of wave period and this burst-like motion can cause the extreme load on the TLP tether. It is often detected in the large and steep irregular waves but the generation mechanism leading to the "ringing" is not yet well understood. According to the research since then, the higher order harmonic components may account for the "ringing" on the floating offshore structures. The main purpose of the present research is, thus, to measure the higher harmonic forces exerted on a vertical truncated circular column and to compare them with available data. A vertical truncated cylinder with a diameter of 3.5inch and a draft of 10.5inch is used as a test structure, which is a scaled model of ISSC TLP column. The cylinder is installed at a distance of 45ft from the wave maker in order to avoid parasitic waves created in the wave flap. Attached to the upper part of the cylinder are two force gages to measure the horizontal (surge) and vertical (heave) forces on the cylinder. The incoming waves are Stokes waves with a slope ranging from 0.06 to 0.24. The forces and waves are measured for 60 seconds with a sampling rate of 50 Hz. Among the recorded data, the first 10 waves are excluded because of transient behavior of the waves and the next The horizontal and vertical forces are analyzed up to 5th order harmonics. The horizontal forces are then compared to the values from the theoretical model called "FNV model". In addition, force transfer functions are also investigated. Major findings in this research are below. 1) The first order forces measured are slightly larger than the theoretical values of "FNV model" 2) The "FNV model" considerably overpredicts the second order forces. 3) The larger the amplitude and more extreme the wave slope, the smaller the predictions are compared to the experimental. 4) The higher harmonic forces are significantly smaller than the first harmonic force for all wave parameters. 5) The normalized forces vs. waves slopes are almost constant in the lower harmonics but vary a lot in the higher harmonics. 6) The trend of forces is more nonlinear in the horizontal forces than in the vertical forces as the wave slope increases. 7) The part of the results above is also observed by other researchers and confirmed again through the present work.

• Journal of Navigation and Port Research
• /
• v.29 no.10 s.106
• /
• pp.853-858
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation have been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use this method for calculations in the case of random waves because long and stable computation (e.g for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand if the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, can be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes' 5th order wave propagation and a random wave propagation.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2005.10a
• /
• pp.153-159
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation rave been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use VOF method for calculations in the case of random waves because long and stable computation ( e.g. for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand of the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, am be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes'5th order wave propagation and a random wave propagation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1511-1521
• /
• 2002

High-speed moist air or steam flow has long been of important subject in engineering and industrial applications. Of many complicated gas dynamics problems involved in moist air flows, the most challenging task is to understand the nonequilibrium condensation phenomenon when the moist air rapidly expands through a flow device. Many theoretical and experimental studies using supersonic wind tunnels have devoted to the understanding of the nonequilibrium condensation flow physics so far. However, the nonequilibrium condensation can be also generated in the subsonic flows induced by the unsteady expansion waves in shock tube. The major flow physics of the nonequilibrium condensation in this application may be different from those obtained in the supersonic wind tunnels. In the current study, the nonequilibrium condensation phenomenon caused by the unsteady expansion waves in a shock tube is analyzed by using the two-dimensional, unsteady, Navier-Stokes equations, which are fully coupled with a droplet growth equation. The third-order TVD MUSCL scheme is applied to solve the governing equation systems. The computational results are compared with the previous experimental data. The time-dependent behavior of nonequilibrium condensation of moist air in shock tube is investigated in details. The results show that the major characteristics of the nonequilibrium condensation phenomenon in shock tube are very different from those in the supersonic wind tunnels.

• 한국연소학회:학술대회논문집
• /
• 2000.12a
• /
• pp.83-91
• /
• 2000

A numerical study was carried out to investigate the 'unstart' process of thermally-choked combustion in model scramjet engines. The combustion mechanism of supersonic combustor will be compared with the experimental results obtained from the T3 free-piston shock tunnel at ANU (Australian National University) and the high enthalpy supersonic wind tunnel at UT (University of Tokyo). For the numerical simulation of supersonic combustion. multi-species Navier-Stokes equations were considered. and detailed chemistry reaction mechanism of $H_2$-Air were adopted. The governing equations were solved by Roe's FDS method and LU-SGS method with MUSCL scheme. In this study. it is found that the thermal choking process could result from excessive heat release due to combustion. In detail, sufficient heat release could be generated at local region of very high temperature increased by reflection of shock waves or vortex sheets. Accordingly the flow of downstream of the combustor fell to subsonic field propagated upstream along the combustor. Sometimes the subsonic flow field propagated into isolator could generate precombustion shock waves in the isolator.

• Journal of computational fluids engineering
• /
• v.7 no.4
• /
• pp.35-41
• /
• 2002

The projectile afterbodies for zero-lift drag reduction has been analyzed using the Navier-Stokes equations with the κ-εturbidence model. The numerical method of a second order upwind scheme has been used on an unstructured adaptive grid system. Base drag reduction methods that have been found effective on axisymmetric bodies are boattailing, base bleed, base combustion, locked vortex afterbodies and multistep afterbodies. In this paper, turbulence flow and pressure charateristics have been studied for geometries of multistep afterbodies. The important geometrical and flow parameters relevant to the design of such afterbodies have been identified by step number, length and height. The flow over multistep aftoerbodies or base have many kinds of compressible flow characteristics including expansion waves at the trailing edge, recompression waves, separation and recirculating flow in the base region, shear flow and wake flow. The numerical results have been compared and analyzed with the experimental data. The flow characteristics have been clearly shown.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.62-66
• /
• 2006

High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation, reattachment, shock waves and expansion waves. The general cavity flow phenomena includes the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions. The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio (L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyzed and compared with the results of Rossiter's Eq.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.525-532
• /
• 2010

In this paper, wave breakers which occur in two dimensional coasts are simulated using a SPH(Smoothed Particle Hydrodynamics) method which represents the movement of fluidic physical volume with particles. As continuative fluid is approximated to the particles, the simulations are performed using fully Lagrangian method without any grid system. Two-dimensional Navier-Stokes equations and continuity equation are used for the numerical simulations. To generate incident waves, a piston type wavemaker is employed. The accuracy of the wave which is numerically generated by the wavemaker is verified by comparing with analytical results. The computations are carried out with various wave heights and slopes. The wave patterns generated through the numerical simulations are compared with several existing experimental and computational results. Agreement between the experimental data and the computation results is comparatively good. Also, the breaker depth index and the breaker height index from the present calculations are compared with the existing experimental results, and the tendency is very similar.