• Journal of Power System Engineering
• /
• v.18 no.3
• /
• pp.79-86
• /
• 2014

Mixed convective vortex flow in the three-dimensional rectangular channel filled with high viscous fluid(Pr=909) is investigated computationally under various operating conditions. The Reynolds number is varied from 0 to $5{\times}10^{-1}$, the Rayleigh number from $10^3$ to $5{\times}10^4$. The three-dimensional governing equations are discretized using the finite volume method. The effects of Reynolds number and Rayleigh number are presented and discussed. From a parametric study, it is found that vortex flow pattern of mixed convection in rectangular channels can be classified into three flow patterns basically, but the new vortex flow structures containing wave rolls are found, which are affected by Rayleigh number and Reynolds number. From this results, we can draw a flow regime map to delineate various vortex flow patterns in the high viscosity fluid mixed convective flow.

• Journal of Environmental Science International
• /
• v.1 no.2
• /
• pp.105.2-117
• /
• 1992

The formation mechanism of the vortex streets in the lee of the mountain Is Investigated by a three-dimensional numerical model. The model is based upon the hydrostatic Boussinesq equations in which the vertical turbulent momentum flux is estimated by a turbulence parameterization scheme, but the horizontal viscosity is assumed to be constant. The results show that Karman vortex streets can form even without surface friction in a constant ambient flow with uniform stratification. The vortex formation is related to breaking of the mountain wave, which depends on the Froude number (Fr). In the case of a three-dimensional bell-shaped mountain, the wave breaking occurs when Fr is less than about 0.8, while a barman vortex forms when Fr is less than about 0.22. Vortex formation also depends on Reynolds number, which is estimated from the horizontal diffusivity. The vortex formation can be explained by the wave saturation theory given by Lindzen (1981) with some modification. Simulations in this study show that in the case of Karman vortex formation the momentum flux in the lower level is much larger than the saturated momentum flux, whereas it is almost equal to the saturated momentum at the upper levels as expected from the saturation theory. As a result, large flux divergence is produced in the lower layer, the mean flow is decelerated behind the mountain, and the horizontal wind shear forms between unmodified ambient wind. The momentum exchange between the mean flow and the mountain wave is produced by the turbulence within a breaking wave. From the result, well developed vortices like Karman vortex can be formed. . The results of the momentum budget calculated by the hydrostatic model are almost the same as nonhydrostatic results as long as horizontal scale of the mountain is 10 km. A well developed barman vortex similar to the hydrostatic one was simulated in the nonhydrostatic case. Therefore, we conclude that the hydrostatic assumption is adequate to investigate the origin of the Km8n vortex from the viewpoint of wave breaking.

• Environmental Sciences Bulletin of The Korean Environmental Sciences Society
• /
• v.1 no.2
• /
• pp.105-117
• /
• 1997

The formation mechanism of the vortex streets in the lee of the mountain is investigated by a three-dimensional numerical model. The model is based upon the hydrostatic Boussinesq equations in which the vertical turbulent momentum flux is estimated by a turbulence parameterization scheme, but the horizontal viscosity is assumed to be constant. The results show that Karman vortex streets can form even without surface friction in a constant ambient flow with uniform stratification. The vortex formation is related to breaking of the mountain wave, which depends on the Froude number (Fr). In the case of a three-dimensional bell-shaped mountain, the wave breaking occurs when Fr is less than about 0.8, while a Karman vortex forms when Fr is less than about 0.22. Vortex formation also depends on Reynolds number, which is estimated from the horizontal diffusivity. The vortex formation can be explained by the wave saturation theory given by Lindzen (1981) with some modification. Simulations in this study show that in the case of Karman vortex formation the momentum flux in the lower level is much larger than the saturated momentum flux whereas it is almost equal to the saturated momentum at the upper levels as expected from the saturation theory. As a result, large flux divergence is produced in the lower layer, the mean flow is decelerated behind the mountain, and the horizontal wind shear forms between unmodified ambient wind. The momentum exchange between the mean flow and the mountain wave is produced by the turbulence within a breaking wave. From the result, well developed vortices like Karman vortex can be formed. The results of the momentum budget calculated by the hydrostatic model are almost the same as nonhydrostatic results as long as horizontal scale of the mountain is 10 km. A well developed Karman vortex similar to the hydrostatic one was simulated in the nonhydrostatic case. Therefore, we conclude that the hydrostatic assumption is adequate to investigate the origin of the Karman vortex from the viewpoint of wave breaking.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.775-786
• /
• 1992

Studies are made of the turbulent separation bubble in a two-dimensional semi-infinite blunt plate aligned to a uniform free stream when the oncoming free stream contains a pulsating component. The discrete-vortex method is applied to simulate this flow situations because this approach is effective to represent the unsteady motions of turbulent shear layer and the effect of viscosity near the solid surface. The two key external paramenters in the free stream, i.e., the amplitude of pulsation, A, and the frequency parameter St[=fH/ $U_{1}$], are dealt with in the present numerical computations, A particular frequency gives a minimum reattachment which is related to the drag reduction and the most effective frequency is dependent on the most amplified shedding frequency. The turbulent flow structure is scrutinized. A comparison between the unperturbed flow and the perturbed at the particular frequency of the minimum reattachment length of the separation bubble suggests that the large-scale structure is associated with the shedding frequency and the flow instabilities.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.20-28
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic mode1. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach arc used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vortex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the k-w SST turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2162-2167
• /
• 2003

It is well-known that high anisotropic characteristic of turbulent flow field is dominant inside tip leakage vortex. This anisotropic nature of turbulence invalidates the use of the conventional isotropic eddy viscosity turbulence model based on the Boussinesq assumption. In this study, to check whether an anisotropic turbulence model is superior to the isotropic ones or not, the results obtained from steady-state Reynolds averaged Navier-Stokes simulations based on the RNG ${\kappa}-{\varepsilon}$ and the Reynolds stress model in two test cases, such as a linear compressor cascade and a forward-swept axial-flow fan, are compared with experimental data. Through the comparative study of turbulence models, it is clearly shown that the Reynolds stress model, which can express the production term and body-force term induced by system rotation without any modeling, should be used to predict the complex tip leakage flow, including the locus of tip leakage vortex center, quantitatively.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1806-1813
• /
• 2003

The stirred tank reactor is one of the most commonly used devices in industry for achieving mixing and reaction. Here we report on results obtained from the large eddy simulations of flow inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius. Stereoscopic PIV measurements (Hill et $al.^{(1)}$) along with the theoretical model of the impeller-induced flow (Yoon et $al.^{(2)}$) were used in defining the impeller-induced flow as superposition of circumferential, jet and tip vortex pair components. Large eddy simulation of flow in a stirred tank was carried out for the three different Reynolds numbers of 4000, 16000 and 64000. The effect of different Reynolds numbers is well observed in both instantaneous and time averaged flow fields. The instantaneous and mean vortex structures are identified by plotting an isosurfaces of swirling strength for all Reynolds numbers. The Reynolds number dependency of the nondimeansional eddy viscosity, resolve scale and subgrid scale dissipations is clearly shown in this study.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.7
• /
• pp.1201-1208
• /
• 2022

In general, to maximize the supply and efficiency of floating offshore wind power generation energy, the motion caused by wave attenuation of the substructure must be reduced. According to previous studies, the motion response was reduced due to the vortex viscosity generated by the damping plate installed in the lower structure among the waves. In this study, a 5 MW semi-submersible OC5 platform and two platforms with attenuation plates were designed, and free decay experiments and numerical calculations were performed to confirm the effect of reducing motion due to vortex viscosity. As a result of the model test, when the heave free decay tests were conducted at drop heights of 30 mm, 40 mm, and 50 mm, compared with the OC5 platform, the platform with two types of damping plates attached had relatively improved motion damping performance. In the model test and numerical calculation results, the damping plate models, KSNU Plate 1 and KSNU Plate 2, were 1.1 times and 1.3 times lower than OC5, respectively, and the KSNU Plate 2 platform showed about two times better damping performance than OC5. This study shows that the area of the damping plate and the vortex viscosity are closely related to the damping rate of the heave motion.

• Particle and aerosol research
• /
• v.11 no.2
• /
• pp.47-55
• /
• 2015

This study compares the results of collection efficiency of difference gas temperature in cyclone dust collector. The previous researcher's experiment results were used to confirm the reliability of CFD(Computational Fluid Dynamics) model. Based on this verified CFD model, we extended the analysis on the cyclone dust collectors. In CFD study, we used RNG k-epsilon model for analysis of turbulence flow, fluid is air, the velocity at inlet is 10 m/s, the temperature of air is $20^{\circ}C$, $100^{\circ}C$, $200^{\circ}C$, $300^{\circ}C$, $600^{\circ}C$ and $1000^{\circ}C$. As the temperature decreases, the average velocity of outer vortex and collection efficiency is increased, showed the highest collection efficiency at $20^{\circ}C$. It can be inferred smooth flow in cyclone dust collector is difficult because air viscosity increases as temperature increases. The power required at $1000^{\circ}C$ is almost 18 times greater than that of $20^{\circ}C$ to get the similar collection efficiency.