• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.4
• /
• pp.33-38
• /
• 2003

The Physics of the Plume-induced shock and separation Particularly at a high Plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with and without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG $\kappa$-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The control methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated The computational results show the main effect of the porous extension on plume-afterbody interactions is to restrain the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect rut plume interference.

• Journal of Ocean Engineering and Technology
• /
• v.27 no.4
• /
• pp.90-97
• /
• 2013

Oil booms are one of the most widely used types of equipment for the protection of coastal areas against oil spills. In some situations, however, there are several types of oil leaks from the oil boom. Important factors regarding these phenomena include the surrounding ocean environment, such as waves, the density and viscosity of oil, the length of the oil boom skirt, etc. To estimate the performance of the oil boom, it is necessary to predict the behavior of the spilled oil and oil boom. In the present study, the prediction of oil boom performance in waves was carried out using the Pusan-National-University-modified Moving Particle Semi-implicit (PNU-MPS) method, which is an improved version of the original MPS proposed by Koshizuka and Oka (1996). The governing equations, which consist of continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators in the governing equations are replaced by the particle interaction models based on a kernel function. The simulation results were validated through a comparison with the results of Violeau et al. (2007)..

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.127-138
• /
• 1999

In this study steady state solutions of cavity flows driven by two moving walls are studied. The north and east walls of the cavity are movable where as the remaining two walls are fixed in space. Numerical experiments for three different driving schemes for moving walls are done at two different Reynolds numbers of Re=40 and 400. The first scheme is to accelerate north and east walls simultaneously. In the second one, the north wall is started first and the east wall is accelerated later. In the third one the east wall starts first. It is usually expected that all these three cases yield the same steady state solution after sufficiently long time. However, present numerical experiments show that such a usual belief is valid only when the Reynolds number is low enough (Re=40). At higher Reynolds number (Re=400), the flow develops to three different steady states depending on the history of the boundary condition change.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.5
• /
• pp.476-483
• /
• 2007

One of the main function of Ieodo Ocean Research Station is to service the information about the weather and fishing grounds condition which are collected through calibrating convection flow and ocean current around the station. However, due to the influence of the station's structure below sea level, it is difficult to obtain the exact flow data. Therefore, it is required to research on the effect of the structure and the method to evaluate and revise the observed data. In this paper, as a basic study, it deals with the algorithm that simulate the interaction between ocean current and the station structure, followed by discussions about the way to applicate the algorithm. Through 3-dimensional computational fluid dynamics analyses (using Navier-Stokes equations with K-turbulence model), the influence of the station and submerged rocks are quantitatively evaluated, and we would suggest methods how to obtain accurate flow information from the measured rough data.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.1-8
• /
• 2005

For study on the unsteady wall interference effect, flows around a circular cylinder in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results showed that the unsteady pressure gradient over the cylinder is enhanced by the wall interference, and as a result the fluctuations of lift and drag are augmented. The drag is further increased because of the lower base pressure. The vortex shedding frequency is also increased by the wall interference. The pressure on the test section wall shows the harmonics having the shedding frequency contained in the wall effect.

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

The effects of Triangle Groove Strips (TGS) on Vortex-induced Vibration (VIV) suppression of marine riser are numerically investigated using Computational Fluid Dynamics (CFD) method. The range of Reynolds number in simulations is 4.0 × 10⁴ < Re < 1.2 × 10⁵. The two-dimensional unsteady Reynolds-Averaged Navier-Stokes (RANS) equations and Shear Stress Transport (SST) k-ω turbulence model are used to calculate the flow around marine riser. The Newmark-β method is employed for evaluating the structure dynamics of marine riser. The effect of the height ratio (ε) of TGS on VIV suppression is evaluated. The amplitude responses, frequency responses, vortex patterns and the flow around the structures are discussed in detail. With the increase of the height ratio of TGS, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When ε=0.04, the suppression effect of TGS is the best. Compared with the VIV responses of smooth marine riser, the amplitude ratio is reduced by 38.9%, the peak of the lift coefficient is reduced by 69% and the peak of the drag coefficient is reduced by 40% when Re=6.0 × 10⁴. With the increase of Reynolds number, the suppression effect of TGS on VIV suppression is improved firstly and then weakened. When the Reynolds number is 7.0 × 10⁴, the amplitude ratio can be reduced by 40.1%. As to the large-amplitude vibration cases, the TGS show nice suppression effect on VIV.