• 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.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.1
• /
• pp.37-44
• /
• 2016

In this study, a parametric study of a centrifugal pump with double volute has been performed numerically using three-dimensional Reynolds-averaged Navier-Stokes equations. The shear stress transport model was selected as turbulence closure through turbulence model test. The finite volume method and unstructured grid system were used for the numerical analysis. The optimal grid system in the computational domain was determined through a grid dependency test. The expansion coefficient, circumferential and radial starting positions and length of divider were selected as the geometric parameters to be tested. And, the hydraulic efficiency and the radial thrust coefficient were considered as performance parameters. It was found that the radial thrust and hydrualic efficiency are more sensitive to the expansion angle and circumferential starting position of the divider than the other geometrical parameters.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.410-415
• /
• 2001

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

• The KSFM Journal of Fluid Machinery
• /
• v.5 no.3 s.16
• /
• pp.40-45
• /
• 2002

Since the previous cut-and-try design algorithm requires much cost and time, the automated design technique with the CFD and optimum design algorithm has recently been concerned. In this work, the Navier-Stokes equation was solved to gain more detailed viscous flow information of cascade with rotor-stator interactions. The H-grid embedded by O-grid was generated to obtain more accurate solution by eliminating the branch cut of H-grid near airfoil surface. To handle the relative motion of the rotor to the stationary stator, the sliding multiblock method was applied and the cubic-spline interpolation was used on the block interface boundary. To validate present procedure, the time-averaged aerodynamic loads were compared with experimeatal data. A good agreement was obtained. The Modified Method of Feasible Direction (MMFD) was used to carry out the sensitivity analysis of the change of aerodynamic performance by the changes of the cascade geometry. The present optimization of the cascade gave a dramatic reduction of the drag while the lift maintains at the value within the user-specified tolerance.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.71-80
• /
• 2010

The violent free-surface motions and the corresponding impact loads are numerically simulated by using the Moving Particle Semi-implicit (MPS) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flows. In the original MPS method, there were several shortcoming including non-optimal source term, gradient and collision models, and search of free-surface particles, which led to less-accurate fluid motions and non-physical pressure fluctuations. In the present study, how those defects can be remedied is illustrated by step-by-step improvements in respective processes of the revised MPS method. The improvement of each step is explained and numerically demonstrated. The numerical results are also compared with the experimental results of Martin and Moyce (1952) for dam-breaking problem. The current numerical results for violent free-surface motions and impact pressures are in good agreement with their experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.1
• /
• pp.49-56
• /
• 1997

Heat-transfer enhancement is seeked through modifications of fin surface. Real life plate-fin heat exchangers have complex three-dimensional geometries. Fins can have arrays of dimples and are attached to rows of penetrating tubes. To isolate the effect of surface modification, we model the real flow by a two-dimensional channel flow with a dimple on one side. The flow is analysed by solving the incompressible Navier-Stokes equation by a finite volume method on a generalized boundary-fitted coordinate. Results show a trapped vortex inside the dimple for all cases computed. Local maximum of Nusselt number occurs near the downstream end of the dimple, due to such a vortex. Location of the vortex does not change with respect to the wall temperature change, but moved downstream when Reynolds number increases. This, together with the results that in all cases vortex core is somewhat downstream of the dimple center, suggests that the mean flow above continuously feeds the kinetic energy to the recirculating flow. Heat transfer enhancement and pressure losses are studied through analysing the relevant dimensionless parameters like, Nusselt number and friction factor. In all cases computed, dimpled channel flow experiences less pressure loss than two-dimensional Poiseuille flow.