• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제30권11호
• /
• pp.1101-1106
• /
• 2006

To investigate a relation between fine scale eddies and particle dispersion in a near-wall turbulence, direct numerical simulations of turbulent channel flow laden particle are performed for $Re_{\tau}$=180. The motions of 0,8 million particles are calculated for several particle response times ($t_p$) which is the particle response time based on stokes’ friction law. The number density of particles has a tendency to increase with approaching the near-wall regions ($y^+$<20) except for cases of very small and large particle response times (i.e. $t_p$=0.02 and 15). Near the wall, the behavior and distribution of particles are deeply associated with the fine scale eddies, and are dependent on particle response times and a distance from the wall. The Stokes number that causes preferential distribution in turbulence is changed by a distance from the wall. The influential Stokes number based on the Burgers' vortex model is derived by using the time scale of the fine scale eddies. The influential Stokes number is also dependent on a distance from the wall and shows large value in the buffer layer.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제22권7호
• /
• pp.1030-1040
• /
• 1998

One of the main unresolved issues in large-eddy simulation(LES) of wall-bounded turbulent flows is the requirement of high spatial resolution in the near-wall region, especially in the spanwise direction. Such high resolution required in the near-wall region is generally used throughout the computational domain, making simulations of high Reynolds number, complex-geometry flows prohibitive. A grid-embedding strategy using a nonconforming spectral domain-decomposition method is proposed to address this limitation. This method provides an efficient way of clustering grid points in the near-wall region with spectral accuracy. LES of transitional and turbulent channel flow has been performed to evaluate the proposed grid-embedding technique. The computational domain is divided into three subdomains to resolve the near-wall regions in the spanwise direction. Spectral patching collocation methods are used for the grid-embedding and appropriate conditions are suggested for the interface matching. Results of LES using the grid-embedding strategy are promising compared to LES of global spectral method and direct numerical simulation. Overall, the results show that the spectral domain-decomposition grid-embedding technique provides an efficient method for resolving the near-wall region in LES of complex flows of engineering interest, allowing significant savings in the computational CPU and memory.

• Journal of the Korean Society for Marine Environment & Energy
• /
• 제9권2호
• /
• pp.79-84
• /
• 2006

The present study deals with the flow over a flat plate with repeated roughness elements of 2-dimensional rectangular shape, which can be applied into the study on the natural geographical roughness and the turbulent flow on roughened solid surface. Experiment was performed using PIV technique in the circulating water channel. Results showed that the flow over roughness elements was characterized by the high shear flow emanating from top of roughness element and the recirculating region formed at the trough of two roughness elements. In general, the ratio between the spacing and the height of roughness elements plays a crucial role in developing the flow pattern near wall surface.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제34권10호
• /
• pp.925-931
• /
• 2010

The paper describes a Large Eddy Simulation (LES) study of turbulent flow around a cavity. A series of three-dimensional cavities placed in a turbulent boundary layer are simulated at a Reynolds number of $1.0{\times}10^5$ by considering U and h, which represent the velocity at the top and the depth of the cavity, respectively. In order to obtain the appropriate solution for the filtered Navier-Stokes equation for incompressible flow, the computational mesh forms dense close to the wall of the cavity but relatively coarse away from the wall; this helps reduce computation cost and ensure rapid convergence. The Boussinesq hypothesis is employed in the subgrid-scale turbulence model. In order to determine the subgrid-scale turbulent viscosity, the Smagorinsky-Lilly SGS model is applied and the CFL number for time marching is set as 1.0. The results show the flow variations inside cavities of different sizes and shapes.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 한국해양공학회 2003년도 추계학술대회 논문집
• /
• pp.204-212
• /
• 2003

As the first step to investigate the nonlinear interactions between turbulence and marine structures inside a viscous NWT, a LES technique was applied to the turbulent channel flow for $Re_{T}=150$, in this paper. The employed models were 4 types, such as the Smagorinsky model, the Dynamic SGS model, the Structure Function model and the Generalized Normal Stress model. The simulated data in time-series for the LESs were averaged in both time and space and performed statistical analysis. And results of the LESs were compared with those of a DNS developed in the present study and two spectral methods by Yoon et al.(2003) & Kim et al.(1987). It seems to be quite difficult to evaluate their performances to the present problem, but is seen that the accuracy of LESs are still related to the number of grids(or fine grid size).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• 제33권11호
• /
• pp.1267-1273
• /
• 2009

Hot press forming can produce high-strength components by rapidly cooling between closed punch and die after hot forming using quenchable boron steel austenized in a furnace. In the hot press forming process, the cooling rate is influenced by the size, position and arrangement of the cooling channel and the file condition of cooling water in the die. Also, mechanical properties of the final components and operation time are related to cooling rate. Therefore, the design of optimized cooling channel is one of the most important works. In this paper, the effect of position and size of the cooling channel on the cooling rate was investigated by using design of experiment and FE analysis in hot press bending process. Therefore the optimum cooling channel ratio was presented in the HPB.

• Journal of Mechanical Science and Technology
• /
• 제16권5호
• /
• pp.703-709
• /
• 2002

In the present paper, approximate coordinate transformations for simulation of turbulent flows with wall deformation, significantly reducing computational cost with little degradation in numerical accuracy, are presented. The Wavier-Stokes equations are coordinate-transformed with an approximation of Tailor-series truncation. The performance is evaluated by performing numerical simulations of a channel flow at Re$\sub$$\tau$/ = 140 with active wall motions of η$\sub$m/$\^$+/ $\leq$5. The approximate transformations provide flow structures as well as turbulence statistics in good agreement with those from a complete transformation [Phys. Fluids 12, 3301 (2000) ] and allow 25-30% savings in the CPU time as compared to the complete one.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제25권11호
• /
• pp.1467-1475
• /
• 2001

In the present paper, approximate coordinate transformations for simulation of turbulent flows with active wall motions, leading to a significant reduction in the computational cost while maintaining the numerical accuracy, are presented: the Navier-Stokes equations are coordinate-transformed with an approximation of Taylor-series truncation and neglect of some less-significant terms. The performance of the proposed transformations is evaluated in simulation of the channel flow at Re$\sub$$\tau$/=140 with wall deformations of │η$\sub$m/$\^$+/ 5. The approximate transformations provide flow structures as wall as turbulence statistics in good agreement with those from a complete coordinate transformation [Phys. Fluids 12, 3301 (2000)] and allow 25-30% savings in the CPU time as compared to the complete one.

• Korea-Australia Rheology Journal
• /
• 제17권3호
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• Proceedings of the KSME Conference
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.811-814
• /
• 2002

Using a mathematical theory, we show that the optimality condition of a turbulent diffuser with maximum pressure recovery at the exit is zero shear stress along the wall. The optimal diffuser shape is designed through iterative procedures by using the $k-{\varepsilon}-{\nu}^{2}-f$ turbulence model for flow simulation. The Reynolds number based on the bulk mean velocity and the channel height at the diffuser entrance is 18,000. We also perform large eddy simulation to validate the shape design results and investigate the flow characteristics near the zero-skin friction wall. Results from large eddy simulation show that the skin friction is slightly higher than zero but is still very small as compared to that of the flat plate boundary layer flow Although the time-averaged wall shear stress is slightly above zero along the diffuser wall, instantaneous flow reversals occur intermittently. The streamwise mein velocity shows an asymptotic behavior of the half-power-law near the wall where the skin friction is close to zero.