• International Journal of Aeronautical and Space Sciences
• /
• 제16권4호
• /
• pp.560-570
• /
• 2015

We present a simulation of a hybrid Reynolds-averaged Navier Stokes / Large Eddy Simulation (RANS/LES) based on detached eddy simulation (DES) for a Burrows and Kurkov supersonic planar mixing experiment. The preliminary simulation results are checked in order to validate the numerical computing capability of the current code. Mesh refinement studies are performed to identify the minimum grid size required to accurately capture the flow physics. A detailed investigation of the turbulence/chemistry interaction is carried out for a nine species 19-step hydrogen-air reaction mechanism. In contrast to the instantaneous value, the simulated time-averaged result inside the reactive shear layer underpredicts the maximum rise in $H_2O$ concentration and total temperature relative to the experimental data. The reason for the discrepancy is described in detail. Combustion parameters such as OH mass fraction, flame index, scalar dissipation rate, and mixture fraction are analyzed in order to study the flame structure.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2023년도 학술발표회
• /
• pp.130-130
• /
• 2023

제트류는 다양한 크기와 운동량의 에디가 복잡하게 혼합되어 이루어져 있으며, 이를 정확하게 모델링하고 이해하기 위해서는 제트류의 다양한 특성들을 잘 반영하여 연구를 수행해야 한다. 다양한 연구 수행 방법 중 수치해석 방법은 상대적으로 공간 및 시간적 비용이 적게 들어서 널리 사용되고 있다. 이러한 수치해석 방법에는 DNS(Direct Numerical Simulation), LES(Large Eddy Simulation), RANS(Reynolds Averaged Navier Stokes) 등이 있으며, 그중 LES는 난류 모델링을 사용하는 RANS 방법에 비해 더욱 정확한 흐름 모델링을 제공하는 장점이 있다. 이러한 LES는 대규모 에디는 직접 해석하면서, 일정 크기 이하의 에디는 모델링을 사용해 해석하는 것이 특징이다. 하지만, LES를 사용하기 위해서는 적절한 그리드 크기를 결정하는 것이 중요하며, 이는 모델의 정확성과 연산 비용에 큰 영향을 미친다. 하지만, 여전히 적절한 그리드 크기를 결정하는 것은 어려운 문제이다. 이러한 LES 모델링을 사용할 때 적절한 그리드 크기를 결정하기 위해서는 정확한 시간 평균 속도 변동을 연구하는 것이 앞서 선행되어야 한다. 따라서, 본 연구에서는 기계학습 기반 접근 방식을 사용하여 난류 제트 내 시간 평균 속도 변동을 예측하는 연구를 진행하였다. 즉, 난류 제트 역학을 이해하는 데 중요한 파라미터인 시간 평균 유속을 이용하여 시간 평균 속도 변동을 예측하는 데 초점을 맞추었다. 모델의 성능은 평균 제곱 오차와 R-제곱 등 다양한 지표를 사용하여 평가되었다.

• 한국수자원학회논문집
• /
• 제45권6호
• /
• pp.545-555
• /
• 2012

컴퓨터 기술의 발달과 더불어 수치해석을 이용한 파랑변형에 대한 연구는 꾸준히 발전하고 있으며 점점 중요한 역할을 수행하고 있다. 하지만 수치모형을 이용한 연구에는 다양한 문제점이 발생할 우려가 있는데, 그 중 가장 빈번하게 발생하는 문제 중의 하나가 파랑의 조파지점에서 발생하는 수치수조내로의 재반사 문제이다. 재반사를 막기 위한 방법으로는 내부조파 기법을 이용하는 것이 일반적이다. Navier-Stokes 방정식 모형에서는 질량 원천항을 이용한 내부조파 기법을 주로 사용해 왔으나, 기존의 연구는 대부분 연직 2차원 수치모형을 이용한 연구에 국한되어 있었다. 그러나 3차원 수치모형을 이용한 연구가 점차 활발해지면서 3차원 Navier-Stokes 방정식 모형의 내부조파 기법에 대한 필요성이 증대되고 있다. 최근 RANS(Reynolds averaged Navier-Stokes) 방정식 모형에서 Boussinesq 방정식의 운동량 원천항을 활용하여 파랑을 내부조파하는 기법이 발표되어 3차원 공간에서 경사지게 입사하는 파랑을 성공적으로 재현하였다. 본 연구에서는 LES(large eddy simulation) 기반의 3차원 Navier-Stokes 방정식 수치모형에 운동량 원천항을 이용한 내부조파 기법을 적용하여 목표파랑을 조파하고 해석해와 비교하여 이를 검증하였다.

• 대한기계학회논문집B
• /
• 제41권6호
• /
• pp.373-380
• /
• 2017

본 연구에서는 레이놀즈 평균 나비어-스톡스(Reynolds-averaged Navier-Stokes, RANS) 방정식을 적용하여 스러스터(thruster) 상호작용에 대한 점성 유동 CFD 해석을 수행하였다. 수치해석은 상용 프로그램인 STAR-CCM+를 이용하였으며, 프로펠러와 No.19a 덕트로 스러스터 모델을 구성하였다. 프로펠러의 회전에 의한 동적인 움직임의 해석을 위해 슬라이딩 격자(sliding mesh)기법을 사용하였으며, 격자구성은 다중영역으로 구분하여 육면체 격자(hexahedral element)로 구성하였고, 계산 시간의 경제성을 고려하여 비등각(non-conformal) 격자를 이용하였다. 스러스터 배치를 그대로 유지한 상태에서 상류방향 스러스터의 방위각(azimuth angle)과 축 기울기에 따라 스러스터 상호작용이 크게 변화하였다. 이러한 결과를 통해 축 기울기와 방위각이 추진성능에 중요한 영향을 미치는 것을 확인하였다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제6권2호
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• 한국유체기계학회 논문집
• /
• 제10권6호
• /
• pp.17-23
• /
• 2007

In the present work, flow analysis has been performed in the steam turbine bypass control valve (single-path type) for two different cases i.e., case with steam only and case with both steam and water. The numerical analysis is performed by solving three-dimensional Reynolds-averaged Navier-Stokes (RANS) equations. The shear stress transport (SST) model and $k-{\varepsilon}$ model are used to each different case as turbulence closure. Symmetry condition is applied at the mid plane of the valve while adiabatic condition is used at the outer wall of the cage. Grid independency test is performed to find the optimal number of grid points. The pressure and temperature distributions on the outer wall of the cage are analyzed. The mass flow rate at maximum plug opening condition is compared with the designed mass flow rate. The numerical analysis of multiphase mixing flow(liquid and vapor) is also performed to inspect liquid-vapor volume fraction of bypass valve. The result of volume fraction is useful to estimate both the safety and confidence of valve design.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.104-111
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, the boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Using the newly modified code, numerical simulations were performed and compared with other computational results as well as the experimental data for the supersonic flows over an oblique shock with a bleed region.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.1377-1382
• /
• 2004

Flow fields in a half ducted propeller fan have been investigated by three-dimensional Reynolds-averaged Navier-Stokes (RANS) simulations and a vortex core identification technique. The simulation at the design operating condition shows that the tip vortex onset point is located at 30 percent tip chord of the suction surface on the blade tip. There is no interaction between the tip vortex and the adjacent blade, so that the tip vortex smoothly convects to the rotor exit. However, the high vorticity in the tip vortex causes the wake and the tip leakage flow to be twined around the tip vortex and to interact with the pressure surface of the adjacent blade. This flow behavior corresponds well with experimental results by Laser Doppler Velocimetry. On the contrary, the simulation at the low-flowrate operating condition shows that the tip vortex onset point is located at the 60 percent tip chord of the suction surface. In contrast to the design operating condition, the tip vortex grows almost tangential direction, and impinges directly on the pressure surface of the adjacent blade.

• 한국유체기계학회 논문집
• /
• 제14권1호
• /
• pp.66-76
• /
• 2011

The heat transfer and pressure loss characteristics of a rotating two-pass channel with a guide vane in the turning region have been studied using three-dimensional Reynolds-averaged Navier-Stokes (RANS) analysis, and the shape of the guide vane has been optimized using surrogate modeling optimization technique. For the optimization, thickness, location and angle of the guide vanes have been selected as design variables. The objective function has been defined as a linear combination of the heat transfer and the friction loss related terms with a weighting factor. Latin hypercube sampling has been applied to determine the design points as design of experiments. A weighted-average surrogate model, PBA has been used as the surrogate model. The guide vane in the turning region does not influence the heat transfer in the first passage upstream of the turning region, but enhances largely the heat transfer in the turning region and the second passage. In an example of the optimization, the objective function has been increased by 13.6%.

• Asian Journal of Atmospheric Environment
• /
• 제11권4호
• /
• pp.270-282
• /
• 2017

To evaluate the deposition amount on a ground surface, mesoscale numerical models coupled with atmospheric chemistry are widely used for larger horizontal domains ranging from a few to several hundreds of kilometers; however, these models are rarely applied to high-resolution simulations. In this study, the performance of a dry and wet deposition model is investigated to estimate the amount of deposition via computational fluid dynamics (CFD) models with high grid resolution. Reynolds-averaged Navier-Stokes (RANS) simulations are implemented for a cone and a two-dimensional ridge to estimate the dry deposition rate, and a constant deposition velocity is used to obtain the dry deposition flux. The results show that the dry deposition rate of RANS generally corresponds to that observed in wind-tunnel experiments. For the wet deposition model, the transport equation of a new scalar concentration scavenged by rain droplets is developed and used instead of the scalar concentration scavenged by raindrops falling to the ground surface just below the scavenging point, which is normally used in mesoscale numerical models. A sensitivity analysis of the proposed wet deposition procedure is implemented. The result indicates the applicability of RANS for high-resolution grids considering the effect of terrains on the wet deposition.