• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.46-47
• /
• 2003

Direct numerical simulations are performed for a turbulent flow subjected to a sudden change in pressure gradient. The calculations are started from a fully-developed turbulent channel flow at $Re_{\tau}=180$. The pressure gradient of the channel flow is then changed abruptly. The responses of the turbulence quantities (e.g., turbulence intensities, Reynolds shear stress, and vorticity fluctuations) and the near-wall turbulence structure to the pressure gradient change are investigated. It is found that there are two different relaxations: a fast relaxation at the early stage and a slow one at the later stage. The early response of the velocity fluctuations shows an anisotropic response of the near-wall turbulence.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권1호
• /
• pp.37-50
• /
• 2016

A large eddy simulation (LES) of a turbulent channel flow is performed by using the third order low-storage Runge-Kutta method in time and second order finite difference formulation in space with staggered grid at a Reynolds number, $Re_{\tau}=590$ based on the channel half width, ${\delta}$ and wall shear velocity, $u_{\tau}$. To reduce the calculation cost of LES, algebraic wall model (AWM) is applied to approximate the near-wall region. The computation is performed in a domain of $2{\pi}{\delta}{\times}2{\delta}{\times}{\pi}{\delta}$ with $32{\times}20{\times}32$ grid points. Standard Smagorinsky model is used for subgrid-scale (SGS) modeling. Essential turbulence statistics of the flow field are computed and compared with Direct Numerical Simulation (DNS) data and LES data using no wall model. Agreements as well as discrepancies are discussed. The flow structures in the computed flow field have also been discussed and compared with LES data using no wall model.

• Journal of Ship and Ocean Technology
• /
• 제9권4호
• /
• pp.1-10
• /
• 2005

The effect of buoyancy orientation on turbulent channel flow has been investigated using DNS (direct numerical simulation). Grashof number is kept at 9.6 $\times 10^{5}$ while changing the orientation of the buoyancy vector to be parallel or perpendicular to the channel walls. Four study cases can be distinguished during this research namely; streamwise, wall-normal unstable stratification, wall-normal stable stratification and spanwise oriented buoyancy. The driving mean pressure gradient used in all cases is adjusted to keep mass flow rate constant while friction Reynolds number is around 150. At this Grashof number, the skin friction shows decrement in the unstable and stable stratification and increment in the other two cases. Analyses of the changes of flow structure for the four cases are presented highlighting on the mean quantities and second order statistics.

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

• Nuclear Engineering and Technology
• /
• 제52권9호
• /
• pp.1945-1954
• /
• 2020

Subchannel code is one of the effective simulation tools for thermal-hydraulic analysis in nuclear reactor core. In order to reduce the computational cost and improve the calculation efficiency, empirical correlation of turbulent mixing coefficient is employed to calculate the lateral mixing velocity between adjacent subchannels. However, correlations utilized currently are often fitted from data achieved in central channel of fuel assembly, which would simply neglect the wall effects. In this paper, the CFD approach based on spectral element method is employed to predict turbulent mixing phenomena through gaps in 3 × 3 bare tight lattice rod bundle and investigate the flow pulsation through gaps in different positions. Re = 5000,10000,20500 and P/D = 1.03 and 1.06 have been covered in the simulation cases. With a well verified mesh, lateral velocities at gap center between corner channel and wall channel (W-Co), wall channel and wall channel (W-W), wall channel and center channel (W-C) as well as center channel and center channel (C-C) are collected and compared with each other. The obvious turbulent mixing distributions are presented in the different channels of rod bundle. The peak frequency values at W-Co channel could have about 40%-50% reduction comparing with the C-C channel value and the turbulent mixing coefficient β could decrease around 25%. corrections for β should be performed in subchannel code at wall channel and corner channel for a reasonable prediction result. A preliminary analysis on fluctuation at channel gap has also performed. Eddy cascade should be considered carefully in detailed analysis for fluctuating in rod bundle.

• 대한기계학회논문집B
• /
• 제29권3호
• /
• pp.314-321
• /
• 2005

Discrete Wavelet Transform (DWT) has been applied to the Direct Numerical Simulation (DNS) data of turbulent channel flow. DWT splits the turbulent flow into two orthogonal parts, one corresponding to coherent structures and the other to incoherent background flow. The coherent structure is extracted from not vorticity field but velocity's since the channel flow is not isoropic. By comparing DWT's result of channel flow with that of isotropic flow, it is shown that coherent structure maintains the properties of original channel flow. The velocity field of coherent structures can be represented by few wavelet modes and that these modes are sufficient to reproduce the velocity probability density function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of DWT, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1813-1818
• /
• 2004

Discrete Wavelet Transform (DWT) has been applied to the Direct Numerical Simulation (DNS) data of turbulent channel flow. DWT splits the turbulent flow into two orthogonal parts, one corresponding to coherent structures and the other to incoherent background flow. The coherent structure is extracted from not vorticity field but velocity's since the channel flow is not isotropic. By comparing DWT's result of channel flow with that of isotropic flow, it is shown that coherent structure maintains the properties of original channel flow. The velocity field of coherent structures can be represented by few wavelet modes and that these modes are sufficient to reproduce the velocity probability distribution function (PDF) and the energy spectrum over the entire inertial range. The remaining incoherent background flow is homogeneous, has small amplitude, and is uncorrelated. These results are compared with those obtained for the same compression rate using large eddy simulation (LES) filtering. In contrast to the incoherent background flow of DWT, the LES subgrid scales have a much larger amplitude and are correlated, which makes their statistical modeling more difficult.

• 한국전산유체공학회지
• /
• 제20권3호
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집B
• /
• pp.461-466
• /
• 2000

Direct numerical simulations (DNS) is carried out to study fully-developed turbulent concentric annular pipe flow with two radius ratios at $Re_{Dh}\;=\;8900$. In case of $R_1/R_2\;=\;0.5$, the present result for the mean flow is in good agreement with the previous experimental data. Because of the transverse curvature effects, the distributions of mean flow and turbulent intensities are asymmetric in contrast to those of other fully-developed flows (channel and pipe flow). From the distributions of skewness of radial velocity fluctuations, it co be identified that all of the characteristics of channel, pipe and turbulent flow on a cylinder in axial flow can be appeared in concentric annular pipe flow.

• 대한기계학회논문집B
• /
• 제25권8호
• /
• pp.1060-1067
• /
• 2001

A large-scale bursting event has been analyzed in a turbulent channel flow using a data obtained from a direct numerical simulation (DNS). Large-scale, plume-like structures have been frequently observed in many experimental results, but the origin of those structures is far from being fully understood. It is believed that those large scale events occur occasionally but contribute significantly to the generation of Reynolds shear stress in the outer layer. This paper attempts to give detailed examples of those large-scale motions observed in a turbulent channel flow at relatively low Reynolds number.