• 동력기계공학회지
• /
• 제18권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.

• 한국환경과학회지
• /
• 제1권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
• /
• 제1권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.

• 대한기계학회논문집
• /
• 제16권4호
• /
• pp.775-786
• /
• 1992

본 연구에서는 DVM의 이론적 배경과 수치계산에 대해 자세히 다루었다. 음향 교란이 가졌을 때의 재부착에 대한 수치해석결과는 Kiya등의 실험결과와 비교하 였으며 만족할만한 일치를 보여주었다. 결과 및 고찰에서는 주로 음향교란이 있는 경우의 재부착길이를 최소화하는 주파수와 교란이 없는 유동의 난류구조해석을 평균속 도 및 압력과 그의 섭동치, 그리고 파워 스펙트럼과 상관계수등을 통해 자세히 비교검 토하였다.

• 한국항공우주학회지
• /
• 제32권8호
• /
• pp.20-28
• /
• 2004

경사충격파와 와류간의 상호작용에 관한 수치적 해석은 경사충격파에 기인한 와류붕괴 의 발생에 관한 이론적 모델 및 실험치와 비교하기 위해 수행되었다. 본 연구의 해석을 위해 마하수에 기초한 Roe의 기법 (RoeM) 과 2-방정식 난류 모델을 이용하여 3차원 난류 유동장을 효과적으로 계산하였다. 자유류 마하수 2.49에서 와유세기, 축방향속도 성분, 충격파 세기의 변화에 따른 상호작용의 영향을 연구하였다. 2-방정식 난류모델인 k-wSST 난류모델과 적절히 모델링된 와류모델을 이용하여 수치연구를 수행함으로써 실험적으로 관찰된 여러 유동장을 정확히 모사할 수 있었다. 경사충격파에 기인한 와류붕괴의 발생에 관한 기준을 결정하기 위한 본 연구결과는 이론적 모델과 실험치와 일치된 결과를 보여주고 있음을 확인하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• 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년 춘계학술대회논문집
• /
• 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).

• 해양환경안전학회지
• /
• 제28권7호
• /
• pp.1201-1208
• /
• 2022

일반적으로 부유식 해상풍력발전 에너지의 수급성과 효율을 극대화하기 위해서는 하부구조물의 파랑 감쇠로 인한 운동을 저감시키는 것이 중요하다. 선행 연구들에 따르면 파도 중 하부구조물에 설치된 감쇠판에 의해서 발생한 와류점성으로 인해 운동 응답이 감소되는 것으로 나타났다. 본 연구에서는 5MW급 반잠수식 OC5 플랫폼과 감쇠판이 부착된 두가지 플랫폼을 설계하고, 와류점성으로 인한 운동저감효과를 확인하기 위해 자유감쇠실험과 수치계산을 수행하였다. 모형시험 결과로 낙하 높이를 30 mm, 40 mm, 50 mm에서의 상하 자유감쇠실험을 수행하였을 때 OC5 플랫폼 대비 두 가지의 형태의 감쇠판이 부착된 플랫폼이 상대적으로 운동감쇠성능이 향상되었다. 모형시험과 수치계산 결과에서 형상화한 감쇠판 모델(KSNU Plate 1, KSNU Plate 2)들이 각각 OC5 대비 상하운동 진폭이 1.1배, 1.3배 각각 감소했으며, KSNU Plate 2 플랫폼은 KSNU Plate 1 플랫폼보다 OC5 대비 약 2배 감쇠성능이 좋아진 것으로 나타났다. 본 연구에서는 감쇠판의 면적과 와류점성이 상하동요의 감쇠율과 밀접한 관련을 보여준다.

• 한국입자에어로졸학회지
• /
• 제11권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.