• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.105-105
• /
• 2019

보나 여수로와 같은 수공구조물의 주변에서 발생하는 흐름 거동은 구조물 모서리에서 발생하는 흐름분리(flow separation)와 이에 따른 전단층(shear layer)과 재순환(recirculation) 흐름 영역의 발달 그리고 분리된 흐름의 재부착(reattachment)이 특징이다. 공학적으로 난류의 해석에 있어서 이러한 흐름 거동들을 정확하게 예측하는 것은 수공구조물 설계에 있어서 중요하다. 이 연구에서는 흐름 분리와 재순환 영역의 발달 그리고 흐름 재부착을 포함하는 후방계단(backward-facing step) 흐름을 155,000의 레이놀즈수 조건에서 하이브리드 RANS/LES 모델을 적용하여 해석결과를 평가한다. 하이브리드 모델로는 벽에 인접한 격자의 해상도에 상대적으로 민감하지 않은 SST(shear-stress transport) 난류 모델을 이용하는 DES(detached-eddy simulation) 기법을 적용하였다. 계단 높이가 h인 계산영역은 흐름방향 길이가 34h, 높이는 계단 상류와 하류에서 각각 1h와 2h 그리고 폭은 $2{\pi}$이다. 계단은 상류단으로부터 10h 하류부 지점에 위치한다. 경계조건으로 상부와 하부 벽면에 대해서는 비활조건을 적용한다. 상류부 수로에서 완전 발달한 흐름을 재현하기 위해서 유입경계조건은 유입부 하류 $2{\pi}h$ 지점에서 계산된 유속과 난류량을 매핑(mapping)기법을 이용하여 반복적으로 적용한다. 총 3.1백만개와 7.3백만개의 셀로 계산영역을 구현한 두 개의 계산격자 그리고 약 3.1백만개의 셀을 이용했지만 벽면 근처에서의 격자 구성을 다른 방식으로 설정한 두 가지 격자를 이용하여 격자 해상도가 DES 수치해석 결과에 미치는 영향을 분석하였다. 수치해석결과는 본 연구에서 상류단 조건으로 적용한 매핑기법이 대상 수로에서 완전 발달한 흐름을 잘 재현함을 보여주며, 합리적인 DES 해석 결과를 얻기 위해서는 벽에 수직한 방향으로 적절한 격자의 해상도와 분포가 필요함을 보여준다.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.3
• /
• pp.45-52
• /
• 2003

Detailed flowfield resulting from the secondary sonic gas injection into a divergent section of supersonic conical nozzle has been numerically investigated. The three-dimensional flowfield associated with the bow-shock/boundary-layer interaction inside the nozzle has been solved by Reynolds-averaged Navier-Stokes equations with an algebraic and $\kappa$-$\varepsilon$ turbulence model. The numerical results have been compared with the experimental results for the identical flow conditions, and it is shown that the comparison is satisfactory Effects of different injection pressures of the secondary jet on the shock/boundary-layer interactions and the overall flow structure inside the nozzle have been investigated. The vortex structures behind the shock interaction and wall pressure variations have also been studied.

• Journal of Korea Water Resources Association
• /
• v.55 no.2
• /
• pp.147-157
• /
• 2022

The backward-facing step (BFS) is a benchmark geometry for analyzing flow separation occurred at the edge and resulting development of shear layer and recirculation zone that are occupied by turbulent flow. It is important to accurately reproduce and analyze the mean flow and turbulence statistics of such flows to design physically stable and performance assurance structure. We carried out 3D RANS computations with widely used, two representative turbulence models, k-ω SST and RNG k-ε, to reproduce BFS flow at the Reynolds number of 23,000 and the Froude number of 0.22. The performance of RANS computations is evaluated by comparing numerical results with an experimental measurement. Both RANS computations with two turbulence models appear to reasonably well reproduce mean flow in the shear layer and recirculation zone, while RNG k-ε computation results in about 5% larger velocity between the outer edge of boundary layer and the free surface above the recirculation zone than k-ω SST computation and experiment. Both turbulence models underestimate the shear stress distribution experimentally observed just downstream of the sharp edge of BFS, while shear stresses computed in the boundary layer downstream of reattachment point are agree reasonably well with experimental measurement. RNG k-ε modeling reproduces better shear stress distribution along the bottom boundary layer, but overestimates shear shear stress in the approaching boundary layer and above the bottom boundary layer downstream of the BFS.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.112-119
• /
• 2002

Three dimensional numerical heat transfer analysis was carried out against periodically arranged fuel injectors of the liquid rocket engine. A finite volume method based on SIMPLE algorithm was adapted which gave a good agreement with the published results of the heat transfer problem of a backward facing step. The Nusselt number and pressure drop increased as the distance between the injector elements decreased. When the Reynolds number increased, the Nusselt number increased but nondimensionalized pressure drop decreased slightly.

• Journal of Advanced Marine Engineering and Technology
• /
• v.18 no.4
• /
• pp.94-101
• /
• 1994

실험 데이터는 급확대비 3:1 팽창의 시험관에서의 실험결과를 나타내고 있으며, 실험에 이용된 동작유체로써는 공기가 사용되었다. 입구관에서 레이놀즈수는 60,000으로부터 120,000까지 변하게 하였고, 스월강도는 0으로부터 16까지 변화되게 하였다. 균일한 열 플럭스 경계조건이 사용되었는데, 그 결과 관벽온도 및 체적온도는 24$^{\circ}C$로부터 71$^{\circ}C$까지에 걸쳐 나타났다. 플롯상에 국소 Nusselt수는 최대 열전달점에서 정점을 이루는 모습을 보여 주고 있다. 스월강도가 0으로부터 최대값으로 증가 되었을때, 최고 Nusselt수의 위치는 시험관에서 4로부터 1스텝 하이트로 변경되는 것이 조사되었다. 이러한 최대 Nusselt수의 상류부 이동은 완전 발달된 유동에서의 값보다 2.2배에서 8.8배나 많은 그의 크기를 증가시킨다고 할 수 있다.

• Tribology and Lubricants
• /
• v.37 no.1
• /
• pp.14-24
• /
• 2021

This study shows the effect of the thermal boundary condition around the tilting pad journal bearing on the static and dynamic characteristics of the bearing through a high-precision numerical model. In many cases, it is very difficult to predict or measure the exact thermal boundary conditions around bearings at the operating site of a turbomachine, not even in a laboratory. The purpose of this study is not to predict the thermal boundary conditions around the bearing, but to find out how the performance of the bearing changes under different thermal boundary conditions. Lubricating oil, bearing pads and shafts were modeled in three dimensions using the finite element method, and the heat transfer between these three elements and the resulting thermal deformation were considered. The Generalized Reynolds equation and three-dimensional energy equation that can take into account the viscosity change in the direction of the film thickness are connected and analyzed by the relationship between viscosity and temperature. The numerical model was written in in-house code using MATLAB, and a parallel processing algorithm was used to improve the analysis speed. Constant temperature and convection temperature conditions are used as the thermal boundary conditions. Notably, the conditions around the bearing pad, rather than the temperature boundary conditions around the shaft, have a greater influence on the performance changes of the bearing.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.28-40
• /
• 2002

Turbulent flows around two-dimensional wing sections in ground effect are analysed by incompressible RANS equations and a finite difference method. The Baldwin-Lomax algebraic turbulence model is used to simulate high Reynolds number flows. The main purpose of this study is to clarify the two-dimensional ground effect and its flow characteristics due to different ground boundary conditions, i.e., moving and fixed bottom boundary. As a first step, to validate the present numerical code, the computational result of Clark-Y(t/C 11.7%) is compared with published numerical results and experimental data. Then, NACA4412 section in ground effect is calculated for various ground clearances with two bottom boundary conditions. According to the computational results, the difference in the lift and moment simulated with the two bottom boundary conditions is negligible, but the drag force simulated by the fixed bottom is to some extent smaller than that by the moving bottom. Therefore, it can be concluded that the drag force measured in a wind tunnel with the fixed bottom could be smaller than that with the moving bottom.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.941-951
• /
• 2013

A three dimensional numerical simulation of laminar flow and heat transfer in fully developed curved pipe flow has been performed to study the effects of Dean number and pipe curvature on the flow and temperature fields under the thermal boundary condition of axially uniform wall heat flux. The Reynolds number under consideration ranges from 100 to 4000, and the Prandtl number is 0.71. The curvature ratios are 0.01, 0.025, 0.05 and 0.1. The axial velocity and temperature profiles and the local Nusselt number obtained from the present study are in good agreement with the previous numerical and experimental results currently available. To show the effects of pipe curvature on the flow and heat transfer, the resistance coefficients and heat transfer coefficients are computed and compared with the results of the previous theoretical and experimental studies. The averaged Nusselt number is correlated with Dean and Prandtl numbers. Furthermore, the critical Reynolds number for transition to turbulent flow is observed to depend upon the curvature ratio.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.566-570
• /
• 2012

본 연구에서는 비선형 k-$\varepsilon$ 모형을 이용하여 직사각형 개수로에서 평균흐름과 난류구조를 모의하였다. 표준 k-$\varepsilon$ 난류모형은 난류의 등방성을 가정하여 국부적 평형상태에서 계산하기 때문에 유선에 따른 레이놀즈 응력의 변형이 큰 경우나 이방성이 강한 경우 이를 계산하지 못한다. 이를 보완하기 위하여 제시된 것이 비선형 k-$\varepsilon$ 난류모형이다. 본 연구에서는 표준 k-$\varepsilon$ 모형과 비선형 k-$\varepsilon$ 모형에 의한 모의결과를 비교하였다. 난류모형을 검증하기 위하여 직사각형 개수로에 흐름을 완전 발달된 등류로 가정하여 해석하였다. 지배방정식을 해석하기 위해 Patankar와 Spalding (1972)이 제시한 SIMPLER 알고리즘을 사용하였고 유한체적법을 이용하여 이산화하고 엇갈린 격자체계를 사용하여 계산에서 발생하는 과도한 진동을 줄였다. 또한 차분기법은 Patankar (1980)가 제시한 Power-law 기법을 채택하였으며 경계조건으로 2층 벽법칙 모형과 Hossain과 Rodi (1993)의 모형을 이용하였다. 두 모형의 적용성을 검증하기 위하여 실측자료를 이용하여 비교하였고 그 결과 비선형 k-$\varepsilon$ 모형이 표준 k-$\varepsilon$ 모형에 비해 좀 더 실측지에 가깝게 모의하는 것을 볼 수 있었다.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.