• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.9
• /
• pp.1051-1057
• /
• 2004

A numerical optimization procedure for the shape of three-dimensional channel with angled ribs mounted on one of the walls to enhance turbulent heat transfer is presented. The response surface method is used as an optimization technique with Reynolds-averaged Wavier-Stokes analyses of flow and heat transfer. SST turbulence model is used as a turbulence closure. Computational results for local heat transfer rate show reasonable agreements with experimental data. The pitch-to-height ratio of the rib and rib height-to-channel height ratio are set to be 9.0 and 0.1, respectively, and width-to-rib height ratio and attack angle of the rib are chosen as design variables. The objective function is defined as a linear combination of heat-transfer and friction-loss related terms with weighting factor. Full-factorial experimental design method is used to determine the data points. Optimum shapes of the channel have been obtained in the range from 0.0 to 0.1 of weighting factor.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.104-110
• /
• 2008

Detached Eddy Simulation (DES) is applied to an axisymmetric base flow at supersonic mainstream. DES is a hybrid approach to modeling turbulence that combines the best features of the Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) approaches. In the Reynolds-averaged mode, the model is currently based on either the Spalart-Allmaras (S-A) turbulence model. In the large eddy simulation mode, it is based on the Smagorinski subgrid scale model. Accurate predictions of the base flowfield and base pressure are successfully achieved by using the DES methodology with less computational cost than that of pure LES and monotone integrated large-eddy simulation (MILES) approaches. The DES accurately resolves the physics of unsteady turbulent motions, such as shear layer rollup, large-eddy motions in the downstream region, small-eddy motions inside the recirculating region. Comparison of the results shows that it is necessary to resolve approaching boundary layers and free shear-layer velocity profiles from the base edge correctly for the accurate prediction of base flows. The consideration of an empirical constant CDES for a compressible flow analysis may suggest that the optimal value of empirical constant CDES may be larger in the flows with strong compressibility than in incompressible flows.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.3
• /
• pp.29-36
• /
• 2002

A numerical analysis based on two-dimensional and three-dimensional incompressible Wavier-Stokes equations has been carried out for double-circular-arc compressor cascades and the results are compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE algorithm adopt pressure weighted interpolation method for non-staggered grid and hybrid scheme for the convective terms. Turbulence modeling is very important for prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. Considering computation times, $\kappa$-$\varepsilon$ turbulence model with wall function is used.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.29-34
• /
• 1999

To develop the efficient numerical optimization method for the design of an airfoil, an evaluation of various methods coupled with two-dimensional Naviev-Stokes analysis is presented. Simplex method and Hook-Jeeves method we used as direct search methods, and steepest descent method, conjugate gradient method and DFP method are used as indirect search methods and are tested to determine the search direction. To determine the moving distance, the golden section method and cubic interpolation method are tested. The finite volume method is used to discretize two-dimensional Navier-Stokes equations, and SIMPLEC algorithm is used for a velocity-pressure correction method. For the optimal design of two-dimensional airfoil, maximum thickness, maximum ordinate of camber line and chordwise position of maximum ordinate are chosen as design variables, and the ratio of drag coefficient to lift coefficient is selected as an objective function. From the results, it is found that conjugate gradient method and cubic interpolation method are the most efficient for the determination of search direction and the moving distance, respectively.

• Nuclear Engineering and Technology
• /
• v.16 no.1
• /
• pp.21-28
• /
• 1984

The reduction of the truncation error including numerical diffusion, has been one of the most important tasks in the development of numerical methods. The stream line method is used to cancel cross numerical diffusion and some of the non-diffusion type truncation error. The two-step stream line method which is the combination of the stream line method and finite difference methods is developed in this work for the solution of the govern ing equations of incompressible buoyant turbulent flow. This method is compared with the finite difference method. The predictions of both classes of numerical methods are compared with experimental findings. Truncation error analysis also has been performed in order to the compare truncation error of the stream line method with that of finite difference methods.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2615-2620
• /
• 2007

This paper describes on pressure drop in a circular pipe of refuse collecting system. The flow characteristics inside the circular pipe are analyzed by three-dimensional Navier-Stokes analysis. In numerical analysis, an organic waste is modeled using the data obtained by site survey. Pressure drop obtained by numerical simulation is compared to the value obtained by experimental measurements for the two kinds of pipe; straight and bended type. The pressure drop obtained by numerical simulation has a good agreement with that of experiments. It is noted that the accurate prediction of pressure drop in the waste pipe is very important to determine the performance of turbo blower used in making a suction pressure in the waste pipe. Especially, the pressure drop for an organic waste is analyzed according to the mass flow rate of waste.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.49-55
• /
• 2011

Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${\alpha}/{\lambda}$) of wavy wall is 0.05. The computational domain is $2{\lambda}{\times}1.05{\lambda}{\times}{\lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.06a
• /
• pp.179.2-179.2
• /
• 2010

In this study, we have constructed the method of design about H-Darrieus wind turbine, a kind of VAWT(vertical axis wind turbine). The NACA 0012 airfoil is chosen for the blade, and DMS(double multiple streamtube) theory is used for the analysis. The flow field is computed with numerical solution of rotating Navier-Stokes equations. From the result of experimental data of power coefficient curves, the validity of the present research is checked. Through the non-dimensional parameter analysis for the wind turbine design, we estimated the efficiency of wind turbine with the resultant Cp's, with which an efficient design of VAWT is achieved, and aerodynamic characteristics are presented systematically.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.3
• /
• pp.183-188
• /
• 2017

This paper provides numerical results for the estimation of the efficiency of KRISO energy saving devices in the design stage. A finite volume method is used to solve Reynolds averaged Navier-Stokes (RANS) equations, where the SST k-$\omega$ model is selected for turbulence closure. The propeller rotating motion is determined using a rigid body motion (RBM) scheme, which is called a sliding mesh technique. The numerical analysis focuses on predicting the power reduction by the designed KRISO devices (K-DUCT) under a self-propulsion condition. The present numerical results show good agreement with the available experimental data. Finally, it is concluded that CFD can be a useful method, along with model tests, for assessing the performance of energy saving devices for propulsion efficiency improvement.