• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.3
• /
• pp.529-561
• /
• 2014

The two-phase turbulent flow past an interface-piercing circular cylinder is studied using a high-fidelity orthogonal curvilinear grid solver with a Lagrangian dynamic subgrid-scale model for large-eddy simulation and a coupled level set and volume of fluid method for air-water interface tracking. The simulations cover the sub-critical and critical and post critical regimes of the Reynolds and sub and super-critical Froude numbers in order to investigate the effect of both dimensionless parameters on the flow. Significant changes in flow features near the air-water interface were observed as the Reynolds number was increased from the sub-critical to the critical regime. The interface makes the separation point near the interface much delayed for all Reynolds numbers. The separation region at intermediate depths is remarkably reduced for the critical Reynolds number regime. The deep flow resembles the single-phase turbulent flow past a circular cylinder, but includes the effect of the free-surface and the limited span length for sub-critical Reynolds numbers. At different Froude numbers, the air-water interface exhibits significantly changed structures, including breaking bow waves with splashes and bubbles at high Froude numbers. Instantaneous and mean flow features such as interface structures, vortex shedding, Reynolds stresses, and vorticity transport are also analyzed. The results are compared with reference experimental data available in the literature. The deep flow is also compared with the single-phase turbulent flow past a circular cylinder in the similar ranges of Reynolds numbers. Discussion is provided concerning the limitations of the current simulations and available experimental data along with future research.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.1134-1147
• /
• 1995

POD (proper orthogonal decomposition) is applied to turbulent leading-edge separation bubble to extract coherent structures. A two-dimensional leading-edge separation bubble is simulated by discrete-vortex method, where a time-dependent source forcing is incorporated. Based on the wealth of numerical data, POD is applied in a range of the forcing amplitude ( $A_{o}$ = 0, 0.5, 1.0 and 1.5) and forcing frequency (0 .leq. $f_{F}$H/ $U_{\infty}$.leq. 0.3). It is demonstrated that the structures of POD have noticeable changes with local forcings. In an effort to investigate the mechanism of decreasing reattachment length, dynamic behaviors of the expansion coefficients and contributions of the eigenfunctions of POD are scrutinized. As the forcing amplitude increases, the large-scale vortex structures are formed near the forcing amplitude increases, the large-scale vortex structures are formed near the separation point and the flow structures become more organized and more regular, accompanying with the reduction of reattachment length. By further inverstigation of POD global entropy, it is seen that the reattachment length is closely linked to the degree of organization of the flow structures.es.s.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.59-72
• /
• 1997

A numerical model is constructed to simulate the interactions of oblique shock wave / turbulent boundary layer on a strongly heated wall. The heated wall temperature is two times higher than the adiabatic wall temperature and the shock wave is strong enough to induce boundary layer separation. The numerical diffusion in the finite volume method is reduced by the use of a higher order convection scheme(UMIST scheme) which is a TVD version of QUICK scheme. The turbulence model is Chen-Kim two time scale model. The comparison of the wall pressure distribution with the experimental data ensures the validity of this numerical model. The effect of strong wall heating enlarges the separation region upstream and downstream. In order to eliminate the separation, wall suction is applied at the shock foot position. The bleeding slot width is about same as the upstream boundary layer thickness and suction mass flow is 10% of the flow rate in the upstream boundary layer. The final configuration of the shock reflection pattern and the wall pressure distribution approach to the non-viscous value when wall suction is applied.

• Advances in aircraft and spacecraft science
• /
• v.10 no.4
• /
• pp.303-315
• /
• 2023

The primary goal of this research is to investigate flow separation phenomena using various turbulence models. Also investigated are the effects of free-stream turbulence intensity on the flow over a NACA 0018 airfoil. The flow field around a NACA 0018 airfoil has been numerically simulated using RANS at Reynolds numbers ranging from 100,000 to 200,000 and angles of attack (AoA) ranging from 0° to 18° with various inflow conditions. A parametric study is conducted over a range of chord Reynolds numbers for free-stream turbulence intensities from 0.1 % to 0.5 % to understand the effects of each parameter on the suction side laminar separation bubble. The results showed that increasing the free-stream turbulence intensity reduces the length of the separation bubble formed over the suction side of the airfoil, as well as the flow prediction accuracy of each model. These models were used to compare the modeling accuracy and processing time improvements. The K- SST performs well in this simulation for estimating lift coefficients, with only small deviations at larger angles of attack. However, a stall was not predicted by the transition k-kl-omega. When predicting the location of flow reattachment over the airfoil, the transition k-kl-omega model also made some over-predictions. The Cp plots showed that the model generated results more in line with the experimental findings.

• Journal of Korea Water Resources Association
• /
• v.48 no.4
• /
• pp.281-290
• /
• 2015

Numerical simulations of free surface flow over a broad-crested rectangular weir are conducted by using the volume of fraction (VOF) method and three different turbulence models, the k-${\varepsilon}$, RNG k-${\omega}$ and k-${\omega}$ SST models. The governing equations are solved by a second-order accurate finite volume method and the grid sensitivity study of solutions is carried out. The numerical results are evaluated by comparing the solutions with experimental and numerical results of Kirkgoz et al. (2008) and some non-dimensionalized experimental results obtained by Moss (1972) and Zachoval et al. (2012). The results show that the present numerical model can reasonably reproduce the experimental results, while three turbulent models yield different numerical predictions of two distinct zones of flow separation, the first zone is in front of the upstream edge of the weir and the second is created immediately behind the upstream edge of the weir where the flow is separated to form the separation bubble. The standard k-${\varepsilon}$ model appears to significantly underestimate the size of both separation zones and the k-${\omega}$ SST model slightly over-estimates the first separation zone in front of the weir. The RNG k-${\varepsilon}$ model predicts both separation zones in overall good agreement with the experimental measurement, while the k-${\omega}$ SST model yields the best numerical prediction of separation bubble at the upstream edge of the weir.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.2
• /
• pp.35-41
• /
• 2009

Both the nozzle expansion ratio and the chamber pressure are simultaneously and continuously changed according to pintle movement, resulting in a different internal flow structure and flow separation characteristics. In this paper, the pintle position effect on nozzle flow structure and separation phenomena is analyzed by experimental-aided Computational Fluid Dynamic(CFD). Among the turbulent models for RANS(Reynold Averaged Navier Stokes) in Fluent, Spalart-Allmaras model is better agreement with the nozzle wall pressure distribution attained by cold-flow test than other models. And even if a conical nozzle is used, there is a shock structure similar to cap-shock pattern mainly occurred in contoured or shaped optimized nozzle because of internal shock generated from pintle tip flow separation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.260-269
• /
• 2000

Oil-film flow visualizations and three-dimensional flow measurements have been conducted in the downstream region of a butterfly-type valve used in air-conditioning systems, with the variation of a disk open angle. The flow visualizations in the flow symmetry plane show that there are a pair of counter-rotating separation/recirculation zones as wall as two jet-like near-wall flows. These flow disturbances are strongly depends on the disk open angle. Based on the flow visualization, a qualitative flow model is suggested in the near-field and downstream region of the valve disk. For a small disk open angle, the mean velocities and turbulent intensities have relatively small values in the near-field of the valve disk, but they do not show uniform distributions even in some downstream region. With an increment of the disk open angle, mean velocity variations and turbulent intensities are greatly increased in the immediate downstream region, but uniform distributions are quickly resumed as departing from the valve disk. The mass flow rate remains nearly constant for the disk open angles less than 30 degrees, meanwhile it strongly depends on the disk open angles between 45 and 75 degrees. The pressure loss is found to be about zero for the disk open angles less than 45 degrees, but is substantially increased for those larger than 75 degrees.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.46-53
• /
• 1997

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.306-311
• /
• 2001

Effect of boundary layer thickness on the flow characteristics around a rectangular prism has been investigated by using a PIV(Particle Image Velocimetry) technique. Three different boundary layers(thick, medium and thin)were generated in the Atmospheric Boundary Layer Wind Tunnel at Pusan National University. The thick boundary layer having 670mm thickness was generated by using spires and roughness elements. The medium thickness of boundary layer$(\delta=270mm)$ was the natural turbulent boundary layer at the test section with fully long developing length(18m). The thin boundary layer with 36.5mm thickness was generated by on a smooth panel elevated 70cm from the wind tunnel floor. The Reynolds number based on the free stream velocity and the height of the model was $7.9{\times}10^3$. The mean velocity vector fields and turbulent kinetic energy distribution were measured and compared. The effect of boundary layer thickness is clearly observed not only in the length of separation bubble but also in the reattachment points. The thinner boundary layer thickness, the higher turbulent kinetic energy peak around the model roof. It is strongly recommended that the height ratio between model and approaching boundary layer thickness should be a major parameter.

• Journal of Mechanical Science and Technology
• /
• v.18 no.2
• /
• pp.272-285
• /
• 2004

The losses at off-design points from a compressor cascade occur due to the deviation from a design incidence angle at the inlet of the cascade. The self-noise from the blade cascade at off-design points comes from a separated boundary layer and vortex sheddings. If the incidence angle to the cascade increases, stalling in blades may occur and the noise level increases significantly. This study applied Large-Eddy Simulations (LES) using deductive and deductive dynamic SGS models to low Mach-number, turbulent flow with each incidence angle to the cascade ranging from -40$^{\circ}$ to +20$^{\circ}$ and compared numerical predictions with measured data. It was observed that the oscillating separation bubbles attached to the suction surface do not modify wake flows dynamically for cases of negative incidence angles. However, an incidence angle greater than 8$^{\circ}$ caused a separated vortex near the leading edge to be shed downstream and created stalling. The computed performance parameters such as drag coefficient and total pressure loss coefficient showed good agreement with experimental results. Noise from the cascade of the compressor is summarized as sound generated by a structure interacting with unsteady, turbulent flows. The hybrid method using acoustic analogy was observed to closely predict the measured overall sound powers and directivity patterns at design and off-design points of blade cascade.