• 대한기계학회논문집B
• /
• 제21권7호
• /
• pp.911-918
• /
• 1997

The turbulent viscous wake flows behind a single airfoil, two-dimensional stationary blade row and three-dimensional rotating blade row were calculated, and the numerical results were compared with experimental ones. The numerical technique was based on the SIMPLE algorithm using three turbulent closure models, standard k-.epsilon. model(WFM), low Reynolds number k-.epsilon. model(LRN) and Reynolds stress model (RSM). In the case of a single airfoil, WFM, LRN and RSM presented fairly good velocity distributions in the wake compared with experimental data. In the case of the stationary blade row, LRN and RSM presented better results than WFM for wake velocity distribution, and especially LRN showed best results among these three turbulent models. In the case of the rotating blade row, WFM and LRN showed fairly good agreement with experimental data of the three-dimensional velocity component distributions in the range from hub to mid span region. LRN was also superior to WFM in accuracy of prediction for the wake velocity distribution as same with the cases of a airfoil and the stationary blade row.

• 대한기계학회논문집B
• /
• 제20권1호
• /
• pp.304-320
• /
• 1996

A multiple production .epsilon. equation model was developed in the low Reynolds number $\kappa$-$\varepsilon$ model with the aids of DNS data. We derived the model theoretically and avoided the use of empirical correlations as much as possible in order for the model to have generality in the prediction of complex turbulent flow. Unavoidable model constants were, however, optimized with the aids of DNS data. All the production and dissipation models in the $\varepsilon$ equation were modified with damping functions to satisfy the wall limiting behavior. A new $f_{\mu}$ function, turbulent diffusion and pressure diffusion model for the k and .epsilon. equations were also proposed to satisfy the wall limiting behavior. By, computational investigation on the plane channel flows, we found that the multiple production model for .epsilon. equation could improve the near wall turbulence behavior compared with the standard production model without the complicated empirical modification. Satisfication of the wall limiting conditions for each turbulence model term was found to be most important for the accurate prediction of near wall turbulence behaviors.

• 대한기계학회논문집
• /
• 제18권10호
• /
• pp.2738-2750
• /
• 1994

The Effects of radiative heat transfer on turbulent flow in a partitioned enclosure is studied numerically. The enclosure is partially divided by a thin, poorly conducting vertical divider projecting from the ceiling of the enclosure. The low Reynolds number $k-{\epsilon}$ model is adopted to calculate the turbulent flow field. The solutions to the radiative transfer equations are obtained by the discrete ordinates method(DOM). This method is based on control volume method and is compatible with the SIMPLER algorithm used to solve the momentum and energy equations. The effects of optical thickness and Planck number on the flow, temperature fields and heat transfer rates are investigated for a moderate Rayleigh number($=10^9$). The changes in buoyant flow fields and temperature distributions due to the variation of baffle length are also analyzed. From the predictions, radiant heat exchange between the baffle and the sidewalls strongly influences the temperature distribution in the baffle and its vicinity and total heat transfer increases as the optical thickness and the baffle length decrease. It is possible to neglect the radiative heat transfer effect when Planck number is over one.

• 대한기계학회논문집B
• /
• 제22권6호
• /
• pp.757-770
• /
• 1998

A Navier-Stokes code with a low Reynolds number k-.epsilon. turbulence model was tested to investigate its predictability for the unsteady transitional boundary layer flow due to rotor-stator interaction. A preliminary calculation with three different numbers of time steps 300, 600, and 1000 for a rotor wake passing period was carried out to see the effects of time steps on the unsteady flow and pressure fields due to rotor-stator interaction. Numerical solutions showed that unsteady pressure was much more sensitive to the number of time steps and over 600 time steps should be used to get a numerical solution independent of the number of time steps for a rotor wake passing period. The original low Reynolds number k-.epsilon. turbulence model showed very poor prediction of the unsteady transitional boundary layer flow due to rotor-stator interaction. This was due to the excessive production of turbulent kinetic energy near the leading edge. A modification suggested by Launder was incorporated and the modified model captured well the wake induced transitional strip. Present solutions also showed improved prediction over previous Euler/boundary layer solution in terms of the onset of unsteady transition and its extent.

• 대한기계학회논문집B
• /
• 제21권6호
• /
• pp.735-746
• /
• 1997

A nonlinear low-Reynolds-number k-.epsilon. model of Park and Sung was extended to predict the flows over a step with inclined wall, where a boundary layer flow without separation and a separated and reattaching flow coexist. For a better prediction of the flows, a slight modification was made on the function of the wall damping( $f_{\mu}$) and the model constant ( $C_{{\epsilon}1}$) in the .epsilon.-equation. The model performance was validated by comparing the model predictions with the experiment. It was shown that the flows over a step with inclined wall are simulated successfully with the present model.ent model.

• 대한기계학회논문집B
• /
• 제22권6호
• /
• pp.771-787
• /
• 1998

A Navier-Stokes code with a modified low Reynolds number k-.epsilon. turbulence model was used to study the unsteady transitional boundary layer flow due to rotor-stator interaction. The modification, proposed by Launder, to improve prediction of stagnation flows was incorporated to the low Reynolds number k-.epsilon. turbulence model by Fan-Lakshminarayana-Barnett. Numerical solution is shown to capture well the calmed laminar flow as well as the wake induced transitional strip due to rotor-stator interaction and shows improvement, in terms of onset of transition and its length, over previous Euler/boundary layer solution. The turbulent kinetic energy shows local maximum along the upstream rotor wake in the wake induced transitional strip and this characteristics is observed untill the end of transition. The wake induced strip also shown apparent even in the laminar sublayer as the upstream rotor wake penetrates inside the boundary layer.

• 대한기계학회논문집B
• /
• 제20권12호
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• 설비공학논문집
• /
• 제14권12호
• /
• pp.1056-1062
• /
• 2002

This paper reports the characteristics of the three dimensional turbulent flow in the 180 degree bends with decreasing cross-sectional area by numerical method. Calculated pressure and velocity, Reynolds stress distributions are compared to the experimental data. Turbulence model employed are low Reynolds number k-epsilon model and algebraic stress model. The results show that the main vortex generated from the inlet part of the bend maintained to outlet of the bend because of the contraction of cross-sectional area. The rate of increase of turbulent kinetic energy through the bend are lower than that of mean flow. Secondary flow strength of the flow is lower about 60% than that of square duct flow.

• 대한기계학회논문집
• /
• 제19권12호
• /
• pp.3329-3343
• /
• 1995

The interaction of mixed convection and surface radiation in a printed circuit board(PCB) is investigated numerically. The electronic equipment is modeled by a two-dimensional channel with three hot blocks. In order to calculate the turbulent flow characteristics, the low Reynolds number k-.epsilon. model which is proposed by Launder and Sharma is applied. The S-4 approximation is used to solve the radiative transfer equation. The effects of the Reynolds number and geometric configuration variation of PCB on the flow and heat transfer characteristics are analyzed. As the results of this study, it is found that the thermal boundary layer occured at adiabatic wall in case with thermal radiation included, and the effect of radiation is also found to be insignificant for high Reynolds numbers. It is found, as well, that the heat transfer increases as the Reynolds number and block space increase and the channel height decreases and the heat transfer of vertical channel is greater than that of horizontal channel.

• 한국항해항만학회지
• /
• 제41권2호
• /
• pp.33-38
• /
• 2017

A marine flapped rudder is designed to improve the effective lift generated by the rudder; this also improves the maneuverability of the ship. The flap is a high lift device installed at the trailing edge of the rudder to augment lift. In this paper, the characteristics of a thick flapped rudder are analyzed at a low Reynolds number with various ratios of flap chord length to total chord length and various aspect ratios, based on the computational fluid dynamics technique. The performance of the rudder with respect to lift, drag, and center of pressure are investigated, and the efficient ratio of flap chord length to total chord length and improved aspect ratio are determined. Ed: highlight - or 'superior'. As a case study, the flow on the flapped rudder of an NACA0021 section shape in free stream condition is simulated. The standard k-epsilon turbulence model is used to model the flow around the flapped rudder. The results indicate that the efficient ratio of the flap chord length to total chord length and aspect ratio are 0.3 and 1.4, respectively.