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

난류흐름 거동은 지형이나 수공구조물과 같은 고체 경계면의 변화에 민감하게 반응하며 특징 또한 다양하다. 보나 여수로 등과 같은 단차 구조물을 통과하는 흐름은 구조물의 모서리 같은 흐름 경계면이 급변하는 지점에서는 흐름분리(flow separation)가 발생하는 것이 특징이다. 이러한 흐름분리로 인해 전단층이 발생하며 흐름 재순환(recirculation)이 구조물 하류부에 형성된다. 이 연구에서는 낙차공 형식의 단차 구조물 하류부에서의 흐름 거동을 이해하기 위해 CFD모델링을 통하여 계산된 3차원 유동장을 분석한다. 난류 모의는 하이브리드 LES(large-eddy simulation)/RANS 계산 기법인 IDDES(improved delayed detached-eddy simulation)기법을 적용한다. IDDES의 기본 모형으로는 k-ω SST모형과 Spalart-Allmaras모형을 이용하여 두 모형의 성능을 평가한다. 자유수면의 변동은 VoF(volume of fluid)기법을 이용하여 계산하며, 각 지배방정식은 최소의 수치분산을 유지하면서 수치해의 안정성을 확보할 수 있는 2차 정확도의 유한체적법을 이용하여 이산화하였다. 수치해석 결과는 레이놀즈수 23,400과 후르드수 0.22의 조건에서 기존에 계측된 자료와 비교하여 수치모형의 정확도를 평가하고 하상 단차 하류부에서의 흐름 거동 특성을 분석한다. 계산 결과는 공학적으로 널리 사용되는 RANS 수치모의에서 볼 수 없는 전단층과 난류구조의 동적 거동 특성과 이에 따른 레이놀즈 응력분포의 특성을 설명해준다.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.2
• /
• pp.39-47
• /
• 1997

In the present paper, viscous flow fields around a wigley hull with appendages are analysed to study interactions between the hull and appendages. Navier-Stokes and continuity equations are solved by a finite volume method in a body-fitted coordinate system which conforms three dimensional ship geometries with appendages. A Sub-Grid Scale(SGS) turbulent model is used for a calculation of high Reynolds number flow. Numerical computations has been done for a Wigley hull form at $Rn=1.0{\times}10^6$. The results show that the present approach can predict, at least in qualitative sense, the influence of the appendages upon the flow field around a ship.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.233-240
• /
• 2017

The flowfield around transonic wing-body configuration was simulated using in-house CFD code and compared with the experimental data to understand the influence of several features of CFD(Computational Fluid Dynamics) ; grid dependency, turbulence models, spatial discretization, and viscosity. The wing-body configuration consists of a simple planform RAE Wing 'A' with an RAE 101 airfoil section and an axisymmetric body. The in-house CFD code is a compressible Euler/Navier-Stokes solver based on unstructured grid. For the turbulence model, the $k-{\omega}$ model, the Spalart-Allmaras model, and the $k-{\omega}$ SST model were applied. For the spatial discretization method, the central differencing scheme with Jameson's artificial viscosity and Roe's upwind differencing scheme were applied. The results calculated were generally in good agreement with experimental data. However, it was shown that the pressure distribution and shock-wave position were slightly affected by the turbulence models and the spatial discretization methods. It was known that the turbulent viscous effect should be considered in order to predict the accurate shock wave position.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.1-7
• /
• 2003

The V2F turbulence model, which has shown very good performance in several test cases at low speeds, has been applied to supersonic ramp flow with 20. corner angle at the free stream Mach number of 2.79. The flow is known to manifest strong shock wave/turbulent boundary layer interactions. As a comparative study, low-Reynolds k-ε models are also considered. While the V2F model predicts wall-pressure distribution well, it relatively predicts larger separation bubble and higher skin-friction after the reattachment than the experimental data. Although the ellpticity of f equation is the characteristics of incompressible flows, the converged solutions are acquired in the compressible flow with shock waves. The effect of the realizability constraints used in the model is also examined. In contrast to the result of impinging jet flows, the realizability bounds proposed by Durbin deterioate the overall solutions of the supersonic ramp flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.58-65
• /
• 1994

It is very important to understand characteristics of solution due to the variation of computational grid sizes, especially when turbulence model not incorporating wall-function is used. The present paper performs numerical investigation on the grid dependency of numerical solution for three dimensional turbulent flow field around a ship. In the present study a finite volume method with a modified sub-grid scale turbulence model and a numerically constructed non-orthogonal curvilinear coordinate system capable of conforming complex ship geometries are used. Numerical studies are then performed for a mathematical Wigley hull and the Series 60, $C_B=0.8$ hull forms. The results for various grid sizes are compared with each other and with measured data to show grid dependencies of numerical solutions.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.9-17
• /
• 1997

The performance of several turbulence models in computing an axisymmetric supersonic base flow is investigated. A compressible Navier-Stokes code, which incorporates k-ε, k-ω model and Reynolds stress closure with three kinds of pressure-strain correlation model, has been developed using implicit LU-SGS algorithm with second-order upwind TVD scheme. Numerical computations have been carried out for Herrin and Dutton's base flow. It is observed that the two-equation models give large backward axial velocity approaching to the base and somewhat larger variation of base pressure distribution than the Reynolds stress model. It is also found that the Reynolds stress model with third order pressure-strain model in the anisotropy tensor predicts most accurate mean flow field.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.9
• /
• pp.6305-6314
• /
• 2015

This study presents numerical simulations of mean flow and turbulence structure of an open channel with submerged vegetation. Filtered Navier-Stokes equations are solved using large-eddy simulation (LES). The immersed boundary method (IBM) is employed based on a Cartesian grid. The numerical result is compared with experimental data of Liu et al. (2008) and shows that simulated results coincided reasonably with experimental data within the average error of 10%. Strong vortices are generated at the interface between vegetated and non-vegetated regions with spanwise extent. The generation of turbulence induced by shear at the interface is interfered with wake turbulence, resulting turbulence intensity maximum. Turbulence produced by shear affects the flow in vegetated region and the penetration depth increases with an increase in the submergence ratio. This result can be used to understand sediment transport mechanisms in the vegetated region.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.30-35
• /
• 2004

A numerical study of underexpanded jet and impingement on a wall mounted at various distances from the nozzle exit is presented. The 3-dimensional Wavier-Stokes equations and κ-ω turbulence equations are solved. The grids are constructed as overlapped grid systems to examine the distance effect. The DADI method is applied to obtain steady-state solutions. To avoid numerical instability such as the carbuncle phenomena that sometimes accompany approximate Riemann solver, the HLLE+ scheme is employed for the inviscid flux at the cell interfaces. A goal of this work is to apply a number of two-equation turbulence models based on the w equation to the impinging jet problem.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.517-524
• /
• 2010

This paper covers the numerical studies performed to investigate the characteristics of turbulent wake generated by a submarine, SUBOFF model. A SUBOFF model assumed as an axial-symmetric body was used to generate wake. The numerical simulation was performed by using a commercial S/W, FLUENT, with the same condition as the experiments by Shin et al.(2009). Mainly the cross-sectional distribution of the time-averaged mean wake and turbulent kinetic energy was compared with the experiments. Both results are agreed well with each other in the propeller wake section, but the agreement between both is not so satisfied in the far wake field. It means that more numerous number of grid points and their concentration should be required in that field.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.