• The KSFM Journal of Fluid Machinery
• /
• v.9 no.4 s.37
• /
• pp.36-42
• /
• 2006

In the present work, characteristics of the flow in the cage of a steam turbine bypass control valve for thermal power plant are investigated. Experimental measurement for wall static pressure has been carried out to validate numerical solutions. And, the flowfield is analyzed by solving steady three-dimensional Reynolds-averaged Navier-Stokes equations. Shear stress transport (SST) model is used as turbulence closure. The effects of the flow area between stages of the cage on the pressure drop are also found.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.152-157
• /
• 2003

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In this study, 3-dimensional unsteady Incompressible Navier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.2
• /
• pp.313-324
• /
• 1992

Two-dimensional Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system, then mapped into a generalized body-fitted coordinate system. All direction of viscous terms are incorporated and turbulent effects are modeled using the Baldwin-Lomax algebraic model. Equation are discretized using finite difference method on the C-type grids and solved using implicit LU-ADI decomposition scheme. Calculations are made at a VKI turbine cascade flow in a transonic wind-tunnel and compared to experimental data. Present numerical scheme is shown to be in good agreement with the previous experimental results and simulates the two-dimensional viscous flow phenomena.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.142-146
• /
• 2017

지진해일의 수리특성을 분석하기 위하여 유사한 파형특성을 가진 고립파를 많이 이용하고 있다. 그러나 고립파의 근사파형은 실제 지진해일에 비해 상당히 좁은 파형분포를 가지고 있다. 이에 수치모의에서는 기존의 고립파 근사식을 개량하여 고립파형의 지진해일을 수치적으로 생성하고 있다. 본 연구에서는 지진해일의 파형분포에 따른 월파특성을 수치적으로 조사하기 위하여 개량된 고립파 근사식을 2차원 N-S solver에서 적용하였다. 이것에 기초하여 수치파동수조에 직립호안과 그 배후에는 월파수조를 설치하고, 지진해일 월파량을 측정하였다. 수치해석결과로부터 직립호안 주변의 공간파형과 마루 위의 유속분포로부터 파형분포에 따른 월파현상을 분석할 수 있었다. 또한 기존 고립파 근사이론 대비 개량된 고립파의 체적비에 따른 월파량 변화를 정량적으로 조사하였다. 그 결과 지진해일의 체적비가 증가할수록 월파량이 거의 선형적으로 증가하는 경향을 나타내었다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.1
• /
• pp.1-7
• /
• 2003

The aerodynamic optimization method for airfoil design was described in this paper. The Navier-Stokes equations were solved to consider the viscous flow information around an airfoil. The Modified Method of Feasible Direction(MMFD) was used for sensitivity analysis and the polynomial interpolation was used for distance calculation of the minimization. The Message Passing Interface(MPI) library of parallel computation was adopted to reduce the computation time of flow solver by decomposing the entire computational domain into 8 sub-domains and one-to-one allocating 8 processors to 8 sub-domains. The parallel computation was also used to compute the sensitivity analysis by allocating each search direction to each processor. The present optimization reduced the drag of airfoil while the lift is maintained at the tolerable design value.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.55-61
• /
• 2008

The Reynolds-Averaged Navier-Stokes(RANS) analysis of the 3-D steady flow around the NREL Phase VI horizontal axis wind turbine(HAWT) rotor was performed. The CFD analysis results were compared with experimental data at several different wind speeds. The present CFD model shows good agreements with the experiments both at low wind speed which formed well-attache flow mostly on the upper surface of the blade, and at high wind speed which blade surface flow completely separated. However, some discrepancy occurs at the relatively high wind speeds where mixed attached and separated flow formed on the suction surface of the blade. It seems that the discrepancy is related to the onset of stall phenomena and consequently separation prediction capability of the current turbulence model. It is also found that strong span-wise flow occurs in stalled area due to the centrifugal force generated by rotation of the turbine rotor and it prevents abrupt reduction of normal force for higher wind speed than the designed value.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.44-44
• /
• 2015

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.37-45
• /
• 2015

A parametric study on the shapes of bobsleigh bumpers has been performed to reduce the aerodynamic drag. Effects of geometric parameters, such as leading angle of leading bumper, the ratio of minimum width to maximum width of leading bumper, the ratio of leading bumper length to trailing bumper length, trailing angle of trailing bumper, and the ratio of bumper height to installation location of bumper from the bottom of bobsleigh, on the aerodynamic performance of the bobsleigh were estimated using 3-D Reynolds-averaged Navier-Stokes equations. The turbulence was analyzed using the shear stress turbulence model. Reynolds number based on the hydraulic diameter of the external flow channel was in the range of 150,000~1,000,000. Numerical results for drag coefficient were validated compared to experimental data. Ranges of the five geometric parameters were determined according to the rule of Federation Internationale de Bobsleigh et de Tobaganning. The aerodynamic performance of the bobsleigh sled was most sensitive to the leading angle of leading bumper and the ratio of minimum width to maximum width of leading bumper.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1457-1463
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a multi-blade centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k - c turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.