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

수력발전 사업에 있어 Desander 구조물은 주로 고산지대 수력발전댐의 Run-of-river 형식의 발전방식에서 유사로 인한 터빈의 손상을 방지하기 위한 목적으로 설치된다. Desander의 적정 규모는 터빈의 손상을 일으킬 수 있는 유사 입경에 대해 안정적으로 침전을 시킬 수 있는 폭/길이/깊이로 평가할 수 있으며 상대적으로 Desander의 규모가 크게 설계된 경우 초기 공사비 증가하고 반대로 규모가 작게 설계된 경우 터빈의 교체 주기 단축으로 인한 유지관리비가 증가된다. 현재까지 일반적인 Desander 구조물의 설계 방식은 제거 입경의 침전 속도, 유입유량 및 깊이를 변수로 사용하여 경험식(L. Sudry method, Guicciardis method 및 Rouse method)을 통해 규모를 결정해 왔다. 하지만, 3-D 전산유체해석을 통해 유속 흐름 분석으로 직 간접적 Desander 규모의 적정성을 평가할 수 있는 현 시점에서 경험식으로부터 도출된 결과의 신뢰성과 객관성을 검증할 필요가 있다고 판단된다. 본 연구에서는 노르웨이 NSTU에서 개발한 유사의 이송 및 확산해석 기능이 내장된 범용 소프트웨어인 SSIIM을 이용하였다. SSIIM(Simulation of Sediment movements In water Intakes with Multiblock)은 개수로 흐름 상태에서 유사 이동 및 하상 변동을 분석할 수 있도록 개발된 3-D 해석 프로그램이다. SSIIM은 수치해석 방법으로 유한체적법(Finite Volume Method)를 채택하였으며 Navier-Stokes equations을 통해 유체의 흐름을 해석한다. 입력 자료는 유입 유량($m^3/sec$), 유입 유사량(kg/sec), 유출부 수위 및 해당 Desander Structure grid 자료가 사용되며 해석 결과로 Desander 내 grid 별 유속, 수위, 유사 농도 변화 등을 제공한다. 본 연구에서는 SSIIM을 이용하여 제거 목표 유사 입경의 차집 효율(Trap efficiency)로 Desander의 적정 규모를 평가 할 수 있는 설계법을 제안하며 설계 단계에서 결정되는 최소 제거유사 입자와 차집 효율에 의한 Desander의 적정 규모 평가 분석을 파키스탄 A 프로젝트를 대상으로 수행하였다. 연구 성과로 (1)SSIIM을 통해 해석된 차집 효율을 기초로 Desander의 적정 규모를 계획할 경우 경험적 방식에 비해 설계의 객관성과 신뢰성을 제고할 수 있으며 (2)3-D 수치해석을 통해 grid 별 유사농도를 확인 할 수 있어 Desander 형상과 규모에 대한 평가가 가능하다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.299-304
• /
• 2008

Starting characteristics of the axi-symmetric second throat exhaust diffuser (STED) with zero-secondary flows are numerically investigated. Renolds-Average Navier-Stokes equations with a standard ${\kappa}-{\varepsilon}$ turbulence model incorporated with enhanced wall treatment are solved to simulate the diffusing evolutions of the nozzle plume. Minimum (optimum) starting pressure difference of 20$\sim$25% between 1-D theory and the measured data validated from previous results[5] is also applied to predict the range of an effective diffuser expansion ratio (Ad/At) in this system.

• Tribology and Lubricants
• /
• v.34 no.6
• /
• pp.247-253
• /
• 2018

Temperature rise due to viscous shear of the lubricating oil generates hydrodynamic pressure, even if the lubricating surfaces are parallel. This effect, known as the thermal wedge effect, varies significantly with film-temperature boundary conditions. The bearing conducts a part of the heat generated; hence, the oil temperature varies with the thermal conductivity of the bearing. In this study, we analyze the effect of thermal conductivity on the thermohydrodynamic (THD) lubrication of parallel slider bearings. We numerically analyze the continuity equation, Navier-Stokes equation, energy equation including the temperature-viscosity and temperature-density relations for lubricants, and the heat conduction equation for bearing by creating a 2D model of the micro-bearing using the commercial computational fluid dynamics (CFD) code FLUENT. We then compare the variation in temperature, viscosity, and pressure distributions with the thermal conductivity. The results demonstrate that the thermal conductivity has a significant influence on THD lubrication characteristics of parallel slider bearings. The lower the thermal conductivity, the greater the pressure generation due to the thermal wedge effect resulting in a higher load-carrying capacity and smaller frictional force. The present results can function as the basic data for optimum bearing design; however, the applicability requires further studies on various operating conditions.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.1
• /
• pp.76-84
• /
• 2019

For decades, solitary waves have commonly been used to simulate tsunami conditions in numerical studies. However, the main component of a tsunami waveform acts at completely different spatial and temporal distributions than a solitary waveform. Thus, this study applied a 2-D numerical wave tank that included a non-reflected tsunami generation system based on Navier-Stokes equations (LES-WASS-2D) to directly simulate the run-up of a tsunami-like solitary wave on a slope. First, the waveform and velocity due to the virtual depth factor were applied to the numerical wave tank to generate a tsunami, which made it possible to generate the wide waveform of a tsunami, which was not reproduced with the existing solitary wave approximation theory. Then, to validate the applied numerical model, the validity and effectiveness of the numerical wave tank were verified by comparing the results with the results of a laboratory experiment on a tsunami run-up on a smooth impermeable 1:19.85 slope. Using the numerical results, the run-up characteristics due to a tsunami-like solitary wave on an impermeable slope were also discussed in relation to the volume ratio. The maximum run-up heights increased with the ratio of the tsunami waveform. Therefore, the tsunami run-up is highly likely to be underestimated compared to a real tsunami if the solitary wave of the approximation theory is applied in a tsunami simulation in a coastal region.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Proceedings of the Acoustical Society of Korea Conference
• /
• spring
• /
• pp.469-474
• /
• 2002

In recent years, modularization of engine parts has increased the application of plastic products in air intake systems. Plastic intake manifolds provide many advantages including reduced weight, contracted cost, and lower intake air temperatures. These manifolds, however, have some weakness when compared with customary aluminium intake manifolds, in that they have low sound transmission loss because of their lower material density. This low transmission loss of plastic intake manifolds causes several problems related to flow noise, especially when the throttle is opened quickly. The physical processes, responsible for this flow noise, include turbulent fluid motion and relative motion of the throttle to the airflow. The former is generated by high-speed airflow in the splits between the throttle valve and the inner-surface of the throttle body and surge-tank, which can be categorized into the quadrupole source. The latter induces the unsteady force on the flow, which can be classified into the dipole source. In this paper, the mechanism of noise generation from the turbulence is only investigated as a preliminary study. Stochastic noise source synthesis method is adopted for the analysis of turbulence-induced, i.e. quadrupole noise by throttle at quick opening state. The method consists of three procedures. The first step corresponds to the preliminary time-averaged Navier-Stokes computation with a $k-\varepsilon$ turbulence model providing mean flow field characteristics. The second step is the synthesis of time-dependent turbulent velocity field associated with quadrupole noise sources. The final step is devoted to the determination of acoustic source terms associated with turbulent velocity. For the first step, we used market available analysis tools such as STAR-CD, the trade names of fluid analysis tools available on the market. The steady state flows at three open angle of throttle valve, i.e. 20, 35 and 60 degree, are numerically analyzed. Then, time-dependent turbulent velocity fields are produced by using the stochastic model and the flow analysis results. Using this turbulent velocity field, the turbulence-originated noise sources, i.e. the self-noise and shear-noise sources are synthesized. Based on these numerical results, it is found that the origin of the turbulent flow and noise might be attributed to the process of formulation and the interaction of two vortex lines formed in the downstream of the throttle valve. These vortex lines are produced by the non-uniform splits between the throttle valve and inner cylinder surface. Based on the analysis, we present the low-noise design of the inner geometry of throttle body.

• Journal of Navigation and Port Research
• /
• v.46 no.6
• /
• pp.491-500
• /
• 2022

It has been noted that an accurate estimation of wind loads on offshore structures such as an FLNG (Liquefied Natural Gas Floating P roduction Storage Offloading Units, LNG FPSOs) with a large topside plays an important role in the safety design of hull and mooring system. Therefore, the present study aims to develop a computational model for estimating the wind load acting on an FLNG. In particular, it is the sequel to the previous research by the author. The numerical computation model in the present study was modified based on the previous research. Numerical analysis for estimating wind loads was performed in two conditions for an interval of wind direction (α), 15° over the range of 0° to 360°. One condition is uniform wind speed and the other is the NPD model reflecting the wind speed profile. At first, the effect of sand-grain roughness on the speed profile of the NPD model was studied. Based on the developed NPD model, mesh convergence tests were carried out for 3 wind headings, i.e. head, quartering, and beam. Finally, wind loads on 6-degrees of freedom were numerically estimated and compared by two boundary conditions, uniform speed, and the NPD model. In the present study, a commercial RANS-based viscous solver, STAR-CCM+ (ver. 17.02) was adopted. In summary, wind loads in surge and yaw from the wind speed profile boundary condition were increased by 20.35% and 34.27% at most. Particularly, the interval mean of sway (45° < α <135°, 225° < α < 315°) and roll (60° < α < 135°, 225° < α < 270°) increased by 15.60% and 10.89% against the uniform wind speed (10m/s) boundary condition.