• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.169-172
• /
• 2002

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.997-1004
• /
• 2002

The current study describes experimental and computational works on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order 10 control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is copuled with two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indrafi transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is signi ficantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.1
• /
• pp.75-90
• /
• 1995

A mathematical model for the dynamic analysis of a riser system with the inclusion of internal flow and nonlinear effects due to large structural displacements is developed to investigate the effect of internal flow on marine riser dynamics. The riser system accounts fir the nonlinear boundary conditions and includes a steady flow inside the pipe which is modeled as an extensible or inextensible. tubular beam subject to nonlinear three dimensional hydrodynamic loads such as current or wave excitation. Galerkin's finite element approximation and time incremental operator are implemented to derive the matrix equation of equilibrium for the finite element system and the extensibility or inextensibility condition is used to reduce degree of freedom of the system and the required computational time in the case of a nonlinear model. The algorithm is implemented to develop computer programs used in several numerical applications. The investigations of the effect of infernal flow on riser vibration due to current or wave loading are performed according to the change of various parameters such as top tension, internal flow velocity, current velocity, wave period, and so on. It is found that the effect of internal flow can be controlled by the increase of top tension. However, careful consideration has to be given in the design point particularly for the long riser under the harmonic loading such as waves. And it is also found that the consideration of nonlinear effects due to large structural displacements increases the effect of internal flow on riser dynamics.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.4B
• /
• pp.243-251
• /
• 2012

To investigate the characteristics of rip current occurring at Haeundae beach, observations obtained from a buoy and a CCTV were analyzed and numerical experiments were conducted. During observed rip-current events, the CCTV images showed that a couple of wave-trains, which are close to regular waves with slightly different directions, propagated to the beach, and wavelet analyses of data from the buoy showed very narrow-banded spectra with a peak frequency. From the evidences, it was inferred that a known mechanism of generating rip current due to the nodal line area of honeycomb-patterned wave crest was one of the significant factors of rip current occurrences of Haeundae beach. The mechanism has been explained by the following: When two wave-trains with slightly different directions propagate to a beach, wave crests of the incident wave-trains form honeycomb pattern due to nonlinear interaction. The nodal lines of honeycomb pattern are developed in the cross-shore direction. And longshore currents flow toward the nodal line area which has very low wave energy. Consequently their mass flux is expelled through the area toward the sea direction. To confirm the generation, numerical experiments were performed using a nonlinear Boussinesq equation model. In the cases with two incident wave-trains with slightly different directions and with a monochromatic wave propagating over a submerged shoal, it was seen that the honeycomb pattern of wave crests was well developed, and thus rip currents were evolved along the nodal lines.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.2
• /
• pp.923-938
• /
• 2019

Surge motion of top-end platform induced by periodic wave makes marine flexible riser encounter equivalent sheared-oscillatory flow, under which the Vortex-induced Vibration (VIV) response will be more complicated than pure sheared flow or oscillatory flow cases. Based on a time domain force-decomposition model, the VIV response characteristics under sheared-oscillatory flows are investigated numerically in this paper. Firstly, the adopted numerical model is validated well against laboratory experiments under sheared flow and oscillatory flow. Then, 20 sheared-oscillatory flow cases with different oscillation periods and top maximum current velocities are designed and simulated. Under long and short oscillation period cases, the structural response presents several similar features owing to the instantaneous sheared flow profile at each moment, but it also has some different patterns because of the differently varying flow field. Finally, the effects and essential mechanism of oscillation period and top maximum current velocity on VIV response are discussed systematically.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.634-639
• /
• 2001

The current study addresses experimental and computational work of impulse wave discharged from the exit of two kinds of right-angle pipe bends, which are attached to the open end of a simple shock tube. The weak normal shock wave with its magnitude of Mach number from 1.02 to 1.20 is employed to obtain the impulse wave propagating outside the exit of the pipe bends. A Schlieren optical system visualizes the impulse wave discharged from the exit of the pipe bends at an instant. The experimental data of the magnitude of the impulse wave and its propagating directivity are analyzed to characterize the impulse waves discharged from the exit of the pipe bends and compared with those discharged from a straight pipe. Computational results well predict the experimented dynamic behaviors of the impulse wave. The results obtained show that a right-angle miter bend considerably reduces the magnitude of the impulse wave and its directivity toward to the pipe axis, compared with the straight pipe and right-angle smooth bend. It is believed that the right-angle miter bend pipe can playa role of a passive control against the impulse wave.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.3
• /
• pp.139-146
• /
• 2017

In order to understand the efficacy of submerged breakwater constructed for the beach protection, laboratory experiments were carried out by observing the characteristics of flow around a single shore-parallel submerged breakwater. The velocity field near the shoreline was measured by utilizing the LSPIV (Large-Scale Particle Image Velocimetry) technique, and mean surface and wave height distributions were observed around the submerged breakwater, according to various combinations of incident waves and submerged breakwaters. In this experiment, it was found that the mean flow pattern behind the submerged breakwater was determined by the balance among the gradients of mean water surface and excess wave-momentum flux (i.e., radiation stress tensors) which interact with the wave-induced current developed by the gradients on the rear and the side of the submerged breakwater. The divergent and convergent flow patterns behind the submerged breakwater (i.e., accretion and erosion response) of the numerical study of Ranasinghe et al.(2010) were observed in the measured velocity distributions, and their empirical formula mostly agreed with the experimental results. However, for some cases in this experiment, it was difficult to say that the flow pattern was one of them and was agreed with the empirical formula.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.53 no.4
• /
• pp.309-316
• /
• 2017

In this paper, numerical modeling is conducted to analyze the tension of an anchor line by varying the size and drag coefficient of a buoy when the trapnet is influenced by the wave and the current simultaneously. A mass-spring model was used to analyze the behavior of trapnet underwater under the influence of waves and current. In the simulation of numerical model, wave height of 3, 4, 5 and 6 m, a period of 4.4 s, and the flow speed of 0.7 m/s were used for the wave and current condition. The drag coefficients of buoy were 0.8, 0.4 and 0.2, respectively. The size of buoy was 100, 50 and 25% based on the cylindrical buoy ($0.0311m^3$) used for swimming crab trap. The drag coefficient of the trapnet, the main model for numerical analysis, was obtained by a circular water channel experiment using a 6-component load cell. As a result of the simulation, the tension of the anchor line decreased proportional to buoy's drag coefficient and size; the higher the wave height, the greater the decrease rate of the tension. When the buoy drag coefficient and size decreased to one fourth, the tension of the anchor line decreased to a half and the tension of the anchor line was lower than the holding power of the anchor even at 6 m of wave height. Therefore, reducing the buoy drag coefficient and size appropriately reduces the trapnet load from the wave, which also reduces the possibility of trapnet loss.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.2
• /
• pp.183-192
• /
• 2004

New dynamic analysis procedures lot the vertically drilled sea water pile are suggested and demonstrated by the typical design Problem. Pile structure submerged in the sea water as well as forces by the ocean waves and tidal currents are modeled and formulated by finite element method. To obtain wave forces for the finite element equation, Airy's wave theory is tested and selected among others. Lateral lifting forces induced by the vortex shedding of current flow is simply based on the harmonic function with the Strouhal frequency and lifting coefficient. Natural frequencies and frequency responses for the pile are calculated by NASTRAN using the results of the formulation. Dynamic displacement and stress results obtained by these procedures are shown to be applicable to predict the dynamic behaviors of the ocean pile by the wave and lifting forces as a preliminary design analysis.