• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.193-198
• /
• 2003

A gas turbine blade has an internal cooling passage equipped with ribs, which can be modeled as a ribbed channel. We have studied a flow inside a ribbed channel using large eddy simulaton (LES) with a dynamic subgrid-scale model. The simulation results are compared with the experimental ones. The turbulence intensity and local heat transfer near the rib have not been well captured by the conventional Reynolds averaged Navier-Stokes simulation (RANS). However, these variables obtained by the present LES agree well with those from experiments. From the instantaneous velocity and temperature fields, we explain the mechanisms responsible for the local peaks in the heat transfer distribution along the channel wall. We have also investigated the effect of rotation on the flow and heat transfer in the ribbed channel.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

• 한국전산유체공학회지
• /
• 제16권2호
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

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

• 대한기계학회논문집B
• /
• 제27권9호
• /
• pp.1290-1298
• /
• 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 PlY 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 non-dimensional eddy viscosity, resolved scale and subgrid scale dissipations is clearly shown in this study.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.819-822
• /
• 2002

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.

• 대한기계학회논문집B
• /
• 제30권5호
• /
• pp.438-445
• /
• 2006

We performed the simulation of the unsteady three dimensional flow over a square cylinder in a wind tunnel in moderate Reynolds number range, $100{\sim}2500$ by using LBM. SGS model was applied for the turbulent flow. Frist of all we compared LBM(Lattice Boltzmann Method) solution of Poiseuille flow applied Farout and bounce back boundary conditions with the analytical and FOAM solutions to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured grids and prescribed uniform velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback and wind tunnel boundary conditions were applied to the cylinder walls and the boundaries of calculation domain respectively. The maximum Strouhal number of the vortex shedding is 0.2025 at Re = 750. and the number maintains the constant value of 0.18 when Re>1000. We also predicted that the critical reynolds number of the turbulent flow is in the range of $250{\sim}500$.

• 대한토목학회논문집
• /
• 제34권5호
• /
• pp.1435-1442
• /
• 2014

본 연구에서는 OpenFOAM에서 제공하는 소스코드를 이용하여 매끄러운 하상의 직사각형 개수로 흐름에 대해 수치모의를 수행하였다. 난류 해석을 위해 LES를 수행하였는데, 비등방성 잔여 응력항을 모델링하기 위해서 Germano et al. (1991)이 제시한 Dynamic Subgrid-scale 모형을 이용하였다. 조직구조를 분석하기 위하여 Lu and Willmarth (1973)가 제시한 uw 사분면기법을 이용하여 순간레이놀즈 응력이 레이놀즈 응력에 미치는 영향을 기여율과 시간비로 나누어 분석하였다. LES 모의 결과를 토대로 기존 실험 및 DNS 모의 결과와 비교하고 분석하였다. 매끈한 하상을 가진 개수로 흐름에서 완충층 이후의 구간에서 분출현상이 쓸기현상에 비해 레이놀즈 응력의 양의 생성에 기여하는 바가 크지만, 분출현상에 비해 쓸기현상의 발생빈도가 큰 것으로 확인되었다.