• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.90-98
• /
• 1998

본 논문은 2차원 해면효과의 수치계산 결과를 정리하였다. 지면으로부터의 높이변화에 따른 점성유동장을 계산하기 위하여 지배방정식으로는 비압축성 RANS방정식을, 시간에 대하여 서는 음해법으로 프로그램을 구성하였다. 압력항은 가상압축성을 도입 4차 수치확산항을 추가하는 것에 의해 계산하였으며, 높은 레이놀즈수에서의 효과적인 계산을 위해 Baldwin-Lomax 난류 모델을 도입하였다. 해면효과가 없는 무한유중에서의 NACA-0012단면 계산결과를 실험데이터와 비교하는 것에 의해 프로그램의 타당성을 확인하였다. NACA-6409와 두께비 $4.6\%$의 날개에 대하여 해면효과를 고려한 계산을 수행하였다. 높이의 변화에 따라 계산된 무차원계수, 압력 및 속도분포는 해면효과의 특성을 잘 보여주고 있다.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.704-708
• /
• 2000

A new efficient numerical method for computing unsteady, incompressible flows, DRANS (Decoupled Reynolds-Averaged Navier-Stokes), is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used fer both the diffusion and convection terms. is adopted. Based on decomposition method, the velocity-turbulent quantity decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully-implicit time advancement scheme. Since the iterative procedures for the momentum, ${\kappa}\;and\;{\varepsilon}$ equations are not required, the components decouplings bring fourth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to turbulent boundary layer with local forcing.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.3
• /
• pp.227-239
• /
• 2017

The accuracy of turbulence models for the Reynolds-Averaged Navier-Stokes (RANS) equations in rough-wall flows is evaluated using data from large-eddy simulations (LES) of boundary layers with favourable and adverse pressure gradients. Some features of the flow (such as flow reversal in the roughness sublayer) cannot be captured accurately by any model, due to the fundamental model formulation. In mild pressure gradients most RANS models are sufficiently accurate for engineering applications, but if strong favourable or adverse pressure gradients are applied (especially those leading to separation) the model performance rapidly degrades.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.2094-2097
• /
• 2007

본 연구에서는 사석경사제의 전면부에 불투과성 수중방파제를 설치하여, 불규칙파의 반사에 대한 수치모의와 수리모형실험을 수행하였다. 수치해석 모형에는 Reynolds Averaged Navier-stokes (RANS) 방정식을 지배방정식으로 사용하였고, 난류해석을 위하여 $\kappa-\varepsilon$ 방정식을 사용하였다. 자유수면변위의 정확한 예측을 위하여 VOF 기법을 적용하였다. 수리모형실험은 한양대학교 수리실험실의 조파수로에서 수행되었다. 본 연구에서는 입사파랑으로 불규칙파를 조파하였으며, Bretsch neider-Mistuyasu 스펙트럼을 목표스펙트럼으로 하여 재현하였다. 반사율의 산정에서 입사파와 반사파를 분리하기 위하여 수치모형실험에서는 3점법을 사용하였고, 수리모형실험에서는 2점법을 사용하였다. 수치모의를 통하여 예측된 반사율과 수리모형실험에서 관측한 결과는 서로 잘 일치하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1287-1292
• /
• 2007

Recently, the frequency of unexpecting heavy rains has been increased due to abnormal climate and extreme rainfall. There was a limit to analyze 1D or 2D stream flow that was applied simple momentum equation and fixed energy conservation. Therefore, hydrodynamics flow analysis in rivers has been needed 3D numerical analysis for correct stream flow interpretation. In this study, CFD model on FLOW-3D was applied to stream flow analysis, which solves three dimenson RANS(Reynolds Averaged Navier-Stokes Equation) control equation to find out physical behavior and the effect of hydraulic structures. Numerical simulation accomplished those results was compared by using turbulence models such as ${\kappa}-{\varepsilon}$, RNG ${\kappa}-{\varepsilon}$ and LES. Those numerical analysis results have been illustrated by the turbulence energy effects, velocity of flow distributions, water level pressure distributions and eddy flows around the piers at Jangwall bridge in urbarn stream.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.6
• /
• pp.427-435
• /
• 2015

This paper provides the numerical results of a simulation of the flow around a propeller working beneath the free surface. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes (URANS) equations, where the wave-making problem is solved using a volume-of-fluid (VOF) method. The numerical analysis focuses on the propeller wake structure affected by the free surface, where we consider another free surface boundary condition that treats the free surface as a rigid wall surface. The propeller wake under the effect of these two free surface conditions shows a reduction in the magnitude of the longitudinal and vertical flow velocities, and its vortical structures strongly interact with the free surface. The thrust and torque coefficient under the free surface effect decrease about 3.7% and 3.1%, respectively. Finally, the present numerical results show a reasonable agreement with the available experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.60-60
• /
• 2021

Investigating Backward-Facing Step(BFS) flow is important in that it is a representative case for separation flows in various engineering flow systems. There have been a wide range of experimental, theoretical, and numerical studies to investigate the flow characteristics over BFS, such as flow separation, reattachment length and recirculation zone. However, most of such previous studies were concentrated only on the perpendicular step angle. In this study, several numerical investigations on the flow pattern over BFS with various step angles (10° ~ 90°) and expansion ratios (1.48, 2 and 3.27) under different Reynolds numbers (5000 ~ 64000) were carried out, mainly focused on the reattachment length. The numerical simulations were performed using an open source 3D CFD software, OpenFOAM, in which the velocity profiles and turbulence intensities are calculated by RANS (Reynolds Averaged Navier-Stokes equation) and 3D LES (Large Eddy Simulation) turbulence models. Overall, it shows a good agreement between simulations and the experimental data by Ruck and Makiola (1993). In comparison with the results obtained from RANS and 3D LES, it was shown that 3D LES model can capture much better and more details on the velocity profiles, turbulence intensities, and reattachment length behind the step for relatively low Reynolds number(Re < 11000) cases. However, the simulation results by both of RANS and 3D LES showed very good agreement with the experimental data for the high Reynolds number cases(Re > 11000). For Re > 11000, the reattachment length is no longer dependent on the Reynolds number, and it tends to be nearly constant for the step angles larger than 30°.) Based on the calibrated and validated numerical simulations, several additional numerical simulations were also conducted with higher Reynolds number and another expansion ratio which were not considered in the experiments by Ruck and Makiola (1993).

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.78-84
• /
• 2016

In this paper, turbulence models for considering roughness in the open source code(OpenFOAM) was investigated. Wall function in the RANS(Reynolds-averaged Navier - Stokes) turbulence model was modified considering roughness on the flat plate by using roughness function. Correlation between the first layer height in the CFD model and roughness height of the plate was observed, and the most proper roughness function, and the first layer height from the plate wall in the CFD analysis was suggested in this paper.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.3
• /
• pp.229-236
• /
• 2016

Exergetic analysis was introduced in optimization of a rotating equilateral triangular internal cooling channel with staggered square ribs to maximize the net exergy gain. The objective function was defined as the net exergy gain considering the exergy gain by heat transfer and exergy losses by friction and heat transfer process. The flow field and heat transfer in the channel were analysed using three-dimensional Reynolds-averaged Navier-Stokes equations under the uniform temperature condition. Shear stress transport turbulence model has been selected as a turbulence closure through the turbulence model test. Computational results for the area-averaged Nusselt number were validated compared to the experimental data. Three design variables, i.e., the angle of rib, the rib pitch-to-hydraulic diameter ratio and the rib width-to-hydraulic diameter ratio, were selected for the optimization. The optimization was performed at Reynolds number, 20,000. Twenty-two design points were selected by Latin hypercube sampling, and the values of the objective function were evaluated by the RANS analysis at these points. Through optimization, the objective function value was improved by 22.6% compared to that of the reference geometry. Effects of the Reynolds number, rotation number, and buoyancy parameter on the heat transfer performance of the optimum design were also discussed.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.157-164
• /
• 2015

The present work performs three-dimensional flow calculations based on Reynolds Averaged Navier-Stokes (RANS) and Delayed Detached Eddy Simulation (DDES) to investigate the flow field of a transonic rotor (NASA Rotor 37) at near-stall condition. It is found that the DES approach is likely to predict well the complex flow characteristics such as secondary vortex or turbulent flow phenomenon than RANS approach, which is useful to describe the flow mechanism of a transonic compressor. Especially, the DES results show improvement of predicting the flow field in the wake region and the model captures reasonably well separated regions compared to the RANS model. Besides, it is discovered that the three-dimensional vortical flows after the vortex breakdown from the rotor tip region are widely distributed and its vortex structures are clearly present. Near the rotor leading edge, a part of the tip leakage flows in DES solution spill over into next passage of the blade owing to the separation vortex flow and the backflow is clearly seen around the trailing edge of rotor tip. Furthermore, the DES solution shows strong turbulent eddies especially in the rotor hub, rotor tip section and the downstream of rotor trailing edge compared to the RANS solution.