• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.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.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.130-130
• /
• 2023

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

• Journal of Korea Water Resources Association
• /
• v.45 no.6
• /
• pp.545-555
• /
• 2012

A numerical modeling has become increasingly popular and more important to the study of water waves with a rapid advancement of computer technology. However, different types of problems are induced during simulating wave motion. One of the key problems is re-reflection to a computation domain at the incident boundary. The internal wave generating-absorbing boundary conditions have been commonly used in numerical wave models to prevent re-reflection. For the Navier-Stokes equations model, the internal wave maker using a mass source function of the continuity equation has been used to generate various types of waves. Nonetheless, almost every numerical experiment is performed in two dimensions and only a few tests have been expanded to three dimensions. More recently, a momentum source function of the Boussinesq equations is applied to generate essentially directional waves in the three dimensional Navier-Stokes equations model. In this study, the internal wave maker using a momentum source function is employed to generate targeted linear waves in the three-dimensional LES model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.6
• /
• pp.373-380
• /
• 2017

In this study, we conducted computational fluid dynamics (CFD) simulations for the unsteady hydrodynamic interaction of multiple thrusters by solving Reynolds averaged Navier-Stokes equations. A commercial CFD software, STAR-CCM+ was used for all simulations by employing a ducted thruster model with combination of a propeller and No. 19a duct. A sliding mesh technique was used to treat dynamic motion of propeller rotation and non-conformal hexahedral grid system was considered. Four different combinations in tilting and azimuth angles of the thrusters were considered to investigate the effects on the propulsion performance. We could find that thruster-hull and thruster-thruster interactions has significant effect on propulsion performance and further study will be required for the optimal configurations with the best tilting and relative azimuth angle between thrusters.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.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.

• The KSFM Journal of Fluid Machinery
• /
• v.10 no.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.04a
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• 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.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.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
• /
• v.11 no.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.