• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1098-1103
• /
• 2004

In this study, the effect of channel flow in the concave surface on local heat transfer characteristics of array jets was investigated experimentally. A TLC method is employed to determine local heat transfer coefficients on the target plate and also flow visualization has been conducted to investigate the behavior of a row of impinging jets and array of impinging jets. Two different array patterns of impinging array jets devices are tested for Reynolds number(Re=10,000). In a row of impinging jets, secondary vortex is strongly maintained by main vortex at nozzle-to-plate distance of H/d=2. Therefore, the Nusselt number slowly decreased at the mid-way region between adjacent jets. In array jets, the local maximum Nusselt number region move further in the downstream direction due to the increase of channel flow velocity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.10 s.241
• /
• pp.1075-1082
• /
• 2005

The characteristics of longitudinal dispersion enhancement by the flow oscillation are numerically studied according to various groove geometries in a 2-D channel in the present study. The length of expanded section l$_{1}$/h$_{1}$ is varied from 0 to 8.75. The oscillating flow condition is given at both side ends, i.e., u = Asin ($2{\pi}ft$) The non-dimensional temperatures at both side ends are set to zero. The bottom and top walls are adiabatic. The local heat sources are located at the middle of the groove wall. In order to solve the governing equations, the SIMPIER algorithm is employed. The present results indicate that maximum longitudinal thermal dispersion can be achieved when the area ratio of the expanded section to the contracted section in the grooved channel becomes 1.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.179-183
• /
• 2006

Surface roughness is present in most of the microfluidic devices due to the microfabrication techniques. This paper presents lattice Boltzmann method (LBM) results for laminar flow in a microchannel with surface roughness. The surface roughness is modeled by an array of rectangular modules placed on top and bottom side of a parallel-plate channel. In this study, LBGK D2Q9 code in lattice Boltzmann Method is used to simulate flow field for low Reynolds number in a micro-channel. The effects of relative surface roughness, roughness distribution, roughness size and the results are presented in the form of the product of friction factor and Reynolds number. Finally, a significant increase in Poiseuille number is detected as the surface roughness is considered, while the effect of roughness on the microflow field depends on the surface roughness.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.118-124
• /
• 2008

The characteristics of heat transfer in a two-dimensional channel obstructed by a square rod is investigated by a turbulence model. The computation is made for the six cases of different rod positions between channel walls. To analyze the wall heat transfer, the heat flux of channel walls is set as a constant value and the $k-{\epsilon}-f_{\mu}$ model is employed. Downstream the square rod, the flow recirculation region appear and they are varied by the rod position. The enhancement of the turbulent heat transfer to the wall is induced by the flow instability using a square rod. The averaged heat transfer rate is maximized at a specific rod position. Finally, the effects of square rod on unsteady flows are scrutinized with the frequency analysis.

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

The commutation errors by the filtering process in the large eddy simulation are considered. It is compared the conventional filter with the inhomogeneous filter that is devised to reduce the commutation errors. The weighting factor of the inhomogeneous filter suggested by Vasilyev is adopted. Also, using the optimizing function that estimates test filter width to eliminate the dissipations in the region excluding the vicinity of the wall, the flow patterns are analyzed. It is evaluated in simulations of the turbulent channel flow at Reynolds number of 1020, based on friction velocity and channel half height. Results show that the commutation errors can be significantly reduced by using the inhomogeneous filter and the optimized test filter width.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.79-82
• /
• 2004

As an alternative numerical method, the lattice Boltzmann method (LBM) is used to simulate a 2-dimensional pressure driven microchannel flow which comes from frequently in MEMS problems. The flow is assumed to be isothermal ideal gas flow. The flow field is calculated with various Knudsen numbers, pressure ratios and aspect ratios of the microchannel. The LBM can show the fundamental characteristics in microchannel flow such as velocity slip and nonlinear pressure drop.

• 한국가시화정보학회:학술대회논문집
• /
• 2007.11a
• /
• pp.84-87
• /
• 2007

Recently, many studies concern on the slip flow and slip length, which allow liquid flow to reduce drag force in microchannel. However, until now not enough investigation is performed experimentally to understand the slip flow in the superhydrophobic microchannel exhibiting riblet structures on vertical wall. Here we investigated and compared the slip flows according to the surface characteristics; hydrophilic, hydrophobic, and superhydrophobic wettabilities. Using the micro-PIV, velocity profiles can be obtained in the glass (hydrophilic), PDMS (hydrophobic), and micro-structured PDMS (superhydrophobic) microchannels. For both PDMS and superhydrophobic PDMS microchannels, we observed the slip effects showing the microscale slip lengths. Due to the micro-riblet, there are two distinctive flow characteristics on the riblet surface and the liquid meniscus in the valleys.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.407-410
• /
• 2006

This study describes a simulation of two-dimensional bubble forming and motion by the Lattice Boltzmann Method with the phase field equation. The free energy model is used to treat the interfacial force and deformation of binary fluids system, drawn into a T-junction the micro channel. A numerical simulation of a binary flow in a cross-junction channel is carried out by using the parallel computation method. The aim in this investigation is to examine the applicability of LBM to numerical analysis of binary fluid separation and motion in the micro channel.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.7
• /
• pp.859-866
• /
• 2003

The Lagrangian dispserion of fluid particles in inhomogeneous turbulence is investigated by a direct numerical simulation of turbulent channel flow. Fluid particle velocity and acceleration along a particle trajectory are computed by employing several interpolation schemes such as linear interpolation, high-order Lagrange polynomial interpolation and the Hermite interpolation schemes. The performances of the schemes are evaluated through comparison of errors in computed particle positions, velocities and accelerations against spectral interpolation. Adopting the four-point Hermite interpolation in the homogeneous directions and Chebyshev polynomials in the wall-normal direction appears to produce most reliable Lagrangian statistics including acceleration correlations with a reasonable amount of computational overhead.