• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.107-115
• /
• 2009

A Computational study was carried out in order to investigate the aerodynamic characteristics of circular cylinder moving near the wavy wall at a low Reynolds number of 50. Lattice Boltzmann method was used to simulate the flow field and immersed boundary method was combined to represent the moving cylinder and wavy wall regardless of the constructed grid in the domain. The aerodynamics characteristics of the cylinder moving near the wavy wall were represented by the comparing the lifting coefficients with various altitudes (H/D) and wave length and amplitudes of wavy wall. It indicated that the twice of increasing-decreasing variations of lifting coefficient are obtained while the cylinder moves near the wavy wall. The first variation is obtained where the cylinder locates near the peak of the wavy wall. Another variation occurs when the distance to the wavy wall becomes longer after passing the peak. It was also classified that three different patterns of relation between the lifting and drag coefficient of the cylinder. However, the classification is limited to the case of the same order of altitude, amplitude and wave length of the wavy wall.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• 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.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1522-1539
• /
• 2002

In the present study, an explicit algebraic stress model is shown to be the exact tensor representation of algebraic stress model by directly solving a set of algebraic equations without resort to tensor representation theory. This repeals the constraints on the Reynolds stress, which are based on the principle of material frame indifference and positive semi-definiteness. An a priori test of the explicit algebraic stress model is carried out by using the DNS database for a fully developed channel flow at Rer = 135. It is confirmed that two-point correlation function between the velocity fluctuation and the Laplacians of the pressure-gradient i s anisotropic and asymmetric in the wall-normal direction. Thus, a novel composite algebraic Reynolds stress model is proposed and applied to the channel flow calculation, which incorporates non-local effect in the algebraic framework to predict near-wall behavior correctly.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.32-38
• /
• 2006

A premixed flame propagating in a tube suffers strong variation in its shape and structure depending on boundary conditions. The effects of thermal boundary conditions and flow fields on flame propagation are numerically investigated. Navier-Stokes equations and species equations are solved with a one-step irreversible global reaction model of methane-air mixture. Finite volume method using an adaptive grid method is applied to investigate the flame structure. In the case of an adiabatic wall, friction force on the wall significantly affected the flame structure while in the case of an isothermal wall, local quenching near the wall dominated flame shapes and propagation. In both cases, variations of flow fields occurred not only in the near field of the flame but also within the flame itself, which affected propagation velocities. This study provides an overview of the characteristics of flames in small tubes at a steady state.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.85-89
• /
• 2003

In a closed square cavity filled with a liquid, a cooled the upper horizontal wall and a heated the lower horizontal wall, the flow isn't generated under the ground-based condition when Rayleigh number is lower than 1700. In such case the flow phenomena near an air bubble under a cooled horizontal wall were investigated. The temperature and the flow fields were studied by using the Thermo-sensitive Liquid-Crystal and the image processing. The qualitative analysis for the temperature and the flow fields were carried out by applying the image processing technique to the original data. Injecting bubble at the center point of upper cooled wall, the symmetry shape of two vortexes near an air bubble was observed. The bubble size increased, the size of velocity and the magnitude of velocity increased. In spite of elapsed time, a pair of two vortexes was the unique and steady-state flow in a square cavity and wasn't induce to the other flow in the surround region.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.397-403
• /
• 2000

In the present study, flow characteristics of turbulent pulsating flow in a square-sectional $180^{\circ}$ curved duct were experimentally investigated. Experimental studies for air flows were conducted to measure axial velocity and wall shear stress distributions and entrance length in a square-sectional $180^{\circ}$ curved duct by using the LDV with the data acquisition and the processing system. The experiment was conducted in seven sections from the inlet (${\phi}=0^{\circ}$) to the outlet (${\phi}=180^{\circ}$) at $30^{\circ}$ intervals of the duct. The results obtained from the experimentation were summarized as follows ; (1) When the ratio of velocity amplitude ($A_1$) was less than one, there was hardly any velocity change in the section except near the wall and any change in axial velocity distributions along the phase. When the ratio of velocity amplitude ($A_1$) was 0.6, the change rate of velocity was slow. (2) Wall shear stress distributions of turbulent pulsating flow were similar to those of turbulent steady flow. The value of the wall shear stress became minimum in the inner wall aid gradually increased toward the outer wall where it became maximum. (3) The entrance length of turbulent pulsating flow reached near the region of bend angle of $90^{\circ}$, like that of turbulent steady flow. The entrance length was changed by the dimensionless angular frequency (${\omega}^+$).

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.

• Journal of Biomedical Engineering Research
• /
• v.13 no.4
• /
• pp.339-352
• /
• 1992

A numerical simulation of the steady and pulsatile flow across the end-to-side anastomosis was performed In order to understand the role of flow dynamics in the preferential bevel opment of distal anastomotic intimal hyperplasla. The finite element technique was employed to solve two-dimensional unsteady pulsatile flow in that region. The results of the steady flow revealed that low shear stresses occur at the proximally occluded host artery and at the recirculation region in the Inner wall just distal to the toe region of the anastomosis. The nor- mal;zed wall shear rate was increased, as was the recirculation zone size in the host artery of the by-pass graft anastomosis, with increased anastomotic junction angle. In order to min imize the size of the low wall shear region which might result in the intimal hyperplasia in the by-pass graft anastomosis, a smaller anastomotic junction angle is recommended. The pulsatile flow simulation revealed flow that regions of low and ascillating mali shear do exist near the anastomosis as In the steady simulation. The shift of stagnation point depends on the pulsation of the flow. As the flow was accelerated at systole, the stagnation point moved downstream, disappered at early diastole and reappeared during late diastole. Low shear stress was also found along both walls of the occluded proximal artery. However, the diastolic flow behavior is quite different from the steady results. The vortex near the occluded artery moved downstream and inwardly during late systole, and disappeared during diastole. Recirculations proximal to the toe and heel regions were significant during diastole. Shear stress oscillation was found along the opposite wall. The results of the present study revealed that tow shear occurs at the proximally occluded host artery aud the recirculation region in the inner wall Just dlstal to the toe region of the anastomosis. The present study suggested that the regions of fluctuated wall shear stress wit flow separation is correlated with the preferential developing regions of anastomosis neointial fibrous hyperplasia.

• Journal of Advanced Marine Engineering and Technology
• /
• v.31 no.2
• /
• pp.166-172
• /
• 2007

The flow around a rotating circular cylinder near a plane wall is investigated by the measurement of the lift acting on the cylinder and by the flow visualization using the hydrogen bubble technique in the circulating water tank. The experimental parameters are the rotating direction of the cylinder. the space ratios $H/D(H/D=0.05{\sim}0.5)$ between cylinder and plane wall and the velocity ratios ${\alpha}({\alpha}=0{\sim}{\pm}2.0)$. In the case of clockwise, the lift on the rotating circular cylinder was increased with the reduction of the space ratios and with the velocity ratios, the upper separation point was more shifted in the rotating direction with them. In the case of anticlockwise, the absolute value of the lift on the rotating circular cylinder was increased with increasing the space ratios and the velocity ratios. the lower separation point was more shifted in the rotating direction with them.