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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.8
• /
• pp.1088-1094
• /
• 2002

A channel flow with a high Reynolds number but coarse grids is numerically studied to investigate the prediction possibility of its turbulence which is three-dimensional and time-dependent. In the present paper, a Reynolds-Averaged Navier-Stokes (RANS) model, a Large Eddy Simulation (LES) and a Navier-Stokes equation with no model are tested with a new approach of hybrid RANS/LES, which reduces to RANS model in the boundary layers and at separation, and to Smagorinsky-like LES downstream of separation, and then compared with each other. It is found that the simulations of hybrid RANS/LES method sustain turbulence like those of LES and with no model, and the results are stable and fairly accurate. This indicates strongly that gradual improvements could lead to a simple, stable, and accurate approach to predict turbulence phenomena of wall-bounded flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.6
• /
• pp.753-760
• /
• 1999

Experiments were conducted to investigate the effect of the initial turbulent intensity on the flow and heat transfer characteristics for a two-dimensional impinging jet. A square rod was installed at the nozzle exit to increase initial turbulent intensity. A hot wire probe and thermochromic liquid crystal technique were used to measure the turbulent intensity and the surface temperature. All measurements were made over a range of nozzle-to-plate distance from 1 to 10 at Re=20,000. When the rod is not installed, the maximum stagnation point Nusselt number is occurred at H/B=9. A higher initial turbulent intensity enhanced the heat transfer on the surface. A correlation between stagnation point Nusselt number and turbulent intensity are presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.128-133
• /
• 2002

The heat transfer characteristics of a turbulent flow in a ribbed two-dimensional channel have been investigated numerically. The fully elliptic governing equations, coupled with a four-equation turbulence model, $\kappa-\omega-\bar{t^2}-\epsilon_t$, are solved by a finite volume method of SIMPLE type. Calculations have been carried out for three rib cross-sections : square, triangular, and semicircular, with various rib pitches and Reynolds numbers. The procedure appears to be satisfactory as the results for the square rib compare favorably with available experimental data and earlier calculation. The optimal rib pitch that yields the maximum heat transfer has been identified. It is also found that the square rib is most effective in enhancing the heat transfer. The semicircular rib, on the other hand, incurs the least amount of pressure drop but the improvement in heat transfer is substantially lower.

• Proceedings of the SAREK Conference
• /
• 2008.11a
• /
• pp.55-60
• /
• 2008

Turbulent flow characteristics on the gap of two parallel channels are investigated using LES(large eddy simulation) approach. Two parallel channels have the same cross-section area and are connected by the narrow channel named the gap. Turbulent flow near the gap makes the flow pulsation along the streamwise direction of two channels. The flow condition is the Reynolds number of $2.5{\times}10^{-5}$. We compared the predicted results with the previous experimental results and presented the axial mean velocity, turbulent intensities, Reynolds shear stresses and turbulent kinetic energy.

• Aerospace Engineering and Technology
• /
• v.2 no.2
• /
• pp.142-150
• /
• 2003

Turbulent flows and related heat transfer in a square heated duct is investigated by a turbulence model and a large eddy simulation. The cooling channel of calorimeter is modeled to the square duct. The nonlinear k-ε-fμ model of Park et al. [3] is slightly modified and their explicit heat flux model is employed. The Reynolds number is varied in the range 4000≤Reｂ≤20000. The heat transfer is closely linked to the secondary flows which driven by the turbulent motion. Its magnitude is 1~3% of the mean streamwise velocity. The relation of Nu~Re0.8Pr0.34 is validated by comparing with the predicted Nu of k-ε-fμ model. Finally, the coherent structures and thermal fluctuations are scrutinized.

• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.5
• /
• pp.638-644
• /
• 2010

Numerical predictions of a fully developed turbulent flow through a square duct ($30mm{\times}30mm$) with twisted tape inserts and with twisted tape plus interrupted ribs are respectively conducted to investigate regionally averaged heat transfer and flow patterns. A rib height-to-channel hydraulic diameter(e/$D_h$) of 0.067 and a lengthto-hydraulic diameter(L/$D_h$) of 30 are considered at Reynolds number ranging 8,900 to 29,000. The interrupted ribs are axially arranged on the bottom wall. The twisted tape is 0.1 mm thick carbon steel sheet with diameter of 28 mm, length of 900 mm, and 2.5 turns. Each wall of the square channel is composed of isolated aluminum sections. Two heating conditions are investigated for test channels with twisted tape inserts and rib turbulators: (1) electric heat uniformly applied to four side walls of the square duct, and (2) electric heat uniformly applied to two opposite walls of the square channel. The results show that uneven surface heating enhances the heat transfer coefficient over uniform heating conditions, and significant improvements can be achieved with twisted tape inserts plus interrupted ribs.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.35-40
• /
• 2010

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number, 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorrector algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.56-59
• /
• 2009

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number is 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorretor algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.