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

For the large scale computation of turbulent flows around an arbitrarily shaped body, a parallel LES (large eddy simulation) code has been recently developed in which domain decomposition method is adopted. METIS and MPI (message Passing interface) libraries are used for domain partitioning and data communication between processors, respectively. For unsteady computation of the incompressible Wavier-Stokes equation, 4-step splitting finite element algorithm [1] is adopted and Smagorinsky or dynamic LES model can be chosen fur the modeling of small eddies in turbulent flows. For the validation and performance-estimation of the parallel code, a three-dimensional laminar flow generated by natural convection inside a cube has been solved. Then, we have solved the turbulent flow around MIRA (Motor Industry Research Association) model at $Re = 2.6\times10^6$, which is based on the model height and inlet free stream velocity, using 32 processors on IBM SMP cluster and compared with the existing experiment.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.573-578
• /
• 2019

Anisotropic distribution of the turbulent kinetic energy and the near-field excitations are the main causes of the steady state Flow-Induced Vibration (FIV) which could lead to fretting wear damage in vertically arranged supported slender rods. In this article, a combined Computational Fluid Dynamics (CFD) and Computational Structural Mechanic (CSM) approach named two-way Fluid-Structure Interaction (FSI) is used to investigate the modal characteristics of a typical rod's vibration. Performance of an Unsteady Reynolds-Average Navier-Stokes (URANS) and Large Eddy Simulation (LES) turbulence models on asymmetric fluctuations of the flow field are investigated. Using the LES turbulence model, any large deformation damps into a weak oscillation which remains in the system. However, it is challenging to use LES in two-way FSI problems from fluid domain discretization point of view which is investigated in this article as the innovation. It is concluded that the near-wall meshes whiten the viscous sub-layer is of great importance to estimate the Root Mean Square (RMS) of FIV amplitude correctly as a significant fretting wear parameter otherwise it merely computes the frequency of FIV.

• International Journal of Fluid Machinery and Systems
• /
• v.3 no.2
• /
• pp.113-121
• /
• 2010

LES(Large Eddy Simulation) with a cavitation model was performed to calculate an unsteady flow for a mixed flow pump with a closed type impeller. First, the comparison between the numerical and experimental results was done to evaluate a computational accuracy. Second, the torque acting on the blade was calculated by simulation to investigate how the cavitation caused the fluctuation of torque. The absolute pressure around the leading edge on the suction side of blade surface had positive impulsive peaks in both the numerical and experimental results. The simulation showed that those peaks were caused by the cavitaion which contracted and vanished around the leading edge. The absolute pressure was predicted by simulation with -10% error. The absolute pressure around the trailing edge on the suction side of blade surface had no impulsive peaks in both the numerical and experimental results, because the absolute pressure was 100 times higher than the saturated vapor pressure. The simulation results showed that the cavitation was generated around the throat, then contracted and finally vanished. The simulated pump had five throats and cavitation behaviors such as contraction and vanishing around five throats were different from each other. For instance, the cavitations around those five throats were not vanished at the same time. When the cavitation was contracted and finally vanished, the absolute pressure on the blade surface was increased. When the cavitation was contracted around the throat located on the pressure side of blade surface, the pressure became high on the pressure side of blade surface. It caused the 1.4 times higher impulsive peak in the torque than the averaged value. On the other hand, when the cavitation was contracted around the throat located on the suction side of blade surface, the pressure became high on the suction side of blade surface. It caused the 0.4 times lower impulsive peak in the torque than the averaged value. The cavitation around the throat caused the large fluctuation in torque acting on the blade.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.93-99
• /
• 2012

The present study focuses on the case of developing flow with in a channel containing a long array of sinusoidal waves (2a/${\lambda}$=0.1, ${\lambda}$=h, ${\lambda}$ is the wavelength, 2a is the wave height, h is the mean channel depth) at the bottom wall. The Reynolds number defined with channel height, h and the mean velocity, U, is Re=6,700. The channel is sufficiently long such that transition is completed and the flow is fully developed over the downstream half of the channel. For the case of an incoming steady flow with no resolved turbulence, the instantaneous flow fields in the transition region are characterized by the formation of arrays of highly-organized large-scale hairpin vortices whose dimensions scale with that of the roughness elements. The paper explains the mechanism for the formation of these arrays of hairpin vortices and shows these eddies play the primary role in the formation of the large-scale streaks of high and low velocity over the wavy wall region. The presence of resolved turbulence in the incoming flow, reduces the streamwise distance needed for the streaks to develop over the wavy region, but does not affect qualitatively the transition process. In the fully-developed region, isolated and trains of large-scale hairpins play an important role in the dynamics of the streaks over the wavy wall.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.1
• /
• pp.12-19
• /
• 2001

Study of turbulent mixing layers has been a popular subject from the point of view of both practical application and phenomenological importance in engineering field. Turbulent mixing layers can be applied in many fields where rapid transition to turbulence is desirable in order to prevent boundary layer separation or to enhance mixing. The ability to control mixing, structure and growth of the shear flow would obviously have a considerable impact on many engineering applications. In addition to practical applications, free shear flows are one of the simplest flows to understand the fundamental mechanism in the transition process to turbulence. After the discovery of large-scale vortical structure in free shear flows many researchers have investigated the physical mechanism of generation and dissipation processes of the vortical structure. This study investigated the role of the large-scale vortical structures in the turbulent mixing layer using LES(Large-Eddy Simulation). The result shows that the pairing interaction of the vortical structure plays an important role in the growth rate of a mixing layer. It is found that the turbulence quantities depend strongly on the velocity ratio. It is also found that the vorticity in the high-velocity-side can extract energy from the mean flow, while the vorticity in the low-velocity-side lose energy by the viscous dissipation. Finally the results suggest the guideline to obtain the desired flow by control of the velocity ratio.

• Atmosphere
• /
• v.25 no.2
• /
• pp.295-308
• /
• 2015

Based on the CMIP5 historical simulation datasets, we assessed the performance of state-of-the-art climate models in respect to the relationship between interannual variabilities of the North Pacific synoptic eddy (NPSE) and East Asian winter monsoon (EAWM). Observation (ERA-Interim) shows a high negative correlation (-0.73) between the interannual variabilities of East Asian winter monsoon (EAWM) intensity and North Pacific synoptic eddy (NPSE) activity during the period of 1979~2005. Namely, a stronger (weaker) EAWM is related to a weaker (stronger) synoptic eddy activities over the North Pacific. This strong reverse relationship can be well explained by latitudinal distributions of the surface temperature anomalies over East Asian continent, which leads the variation of local baroclinicity and significantly weakens the baroclinic wave activities over the northern latitudes of $40^{\circ}N$. This feature is supported by the distribution of the meridional heat flux (${\overline{{\nu}^{\prime}{\theta}^{\prime}}}$) anomalies, which have negative (positive) values along the latitudes $40{\sim}50^{\circ}N$ for strong(weak) EAWM years. In this study, the historical simulations by 11 CMIP5 climate models (BCC-CSM1.1, CanESM2, GFDL-ESM2G, GFDL-ESM2M, HadGEM2-AO, HadGEM2-CC, IPSL-CM5A-LR, MPI-ESM-LR, MPI-ESM-MR, MRI-CGCM3, and NorESM1-M) are analyzed for DJF of 1979~2005. Correlation coefficient between the two phenomena is -0.59, which is comparable to that of observation. Model-to-model variation in this relationship is relatively large as the range of correlation coefficient is between -0.76 (HadGEM2-CC and HadGEM2-AO) and -0.33 (MRI-CGCM3). But, these reverse relationships are shown in all models without any exception. We found that the multi-model ensemble is qualitatively similar to the observation in reasoning (that is, latitudinal distribution of surface temperature anomalies, variation of local baroclinicity and meridional heat flux by synoptic eddies) of the reverse relationship. However, the uncertainty for weak EAWM is much larger than strong EAWM. In conclusion, we suggest that CMIP5 models as an ensemble have a good performance in the simulation of EAWM, NPSE, and their relationship.

• Journal of the Korean Society for Railway
• /
• v.10 no.2 s.39
• /
• pp.224-230
• /
• 2007

The performance and applicability of FDS code is analyzed for flow simulation in railway tunnel. FDS has been built in NIST(USA) for simulation of fire-driven flow. RANS and DNS's results are compared with FDS's. AJL non-linear ${\kappa}-{\epsilon}$[7,8] model is employed to calculate the turbulent flow for RANS. DNS data by Moser et al.[9] are used to prove the FDS's applicability in the near wall region. Parallel plate is used for simplified model of railway tunnel. Geometrical variables are non-dimensionalized by the height (H) of parallel plate. The length of streamwise direction is 50H and the length of spanwise direction is 5H. Selected Re numbers are 10,667 for turbulent flow and 133 for laminar low. The characteristics of turbulent boundary layer are introduced. AJL model's predictions of turbulent boundary layer are well agreed with DNS data. However, the near wall turbulent boundary layer is not well resolved by FDS code. Slip conditions are imposed on the wall but wall functions based on log-law are not employed by FDS. The heavily dense grid distribution in the near wall region is necessary to get correct flow behavior in this region for FDS.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.3
• /
• pp.156-163
• /
• 2009

In this study, comparative analysis on the back-layer phenomena in the tunnel-fire driven flow is performed using numerical simulation with LES and RANS. FDS(Fire Dynamics Simulator) code is employed to calculate the fire-driven turbulent flow for LES and Smartfire code is used for RANS. Hwang and Wargo's data of scaling tunnel fire experiment are employed to compare with the present numerical simulation. The modeled tunnel is 5.4m(L) ${\times}$ 0.4m(W) ${\times}$ 0.3m(H). Heat Release Rate (HRR) of fire is 3.3kW and ventilation-velocity is 0.33m/s in the main stream. The various grid-distributions are systematically tested with FDS code to analyze the effects of grid size. The LES method with FDS provides an improved back-layer flow behavior in comparison with the RANS (${\kappa}-{\epsilon}$) method by Smartfire. The FDS solvers, however, overpredict the velocity in the center region of flow which is caused by the defects in the tunnel-entrance turbulence strength and in the near-wall turbulent flow in FDS code.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.49-57
• /
• 2013

Analysis of fluid-structure interaction for two nearby underwater vehicles immersed in the sea is quite challenging because simulation of flow around them is very difficult due to the complexity of underwater vehicle shapes. The conventional approach using body-fitted or unstructured grids demands much time in dynamic grid generation, and yields slow convergence of solution. Since an analysis of fluid-structure interaction must be based on accurate simulation results, a more efficient way of simulating flow around underwater vehicles, without sacrificing accuracy, is desirable. An immersed boundary method facilitates implementation of complicated underwater-vehicle shapes on a Cartesian grid system. An LES modeling is also incorporated to resolve turbulent eddies. In this paper, we will demonstrate the effectiveness of the immersed boundary method we adopted, by presenting the simulation results on the flow around a modeled high-speed underwater vehicle interacting with a modeled low-speed one.

• Wind and Structures
• /
• v.4 no.3
• /
• pp.177-196
• /
• 2001

At present the most popular turbulence models used for engineering solutions to flow problems are the $k-{\varepsilon}$ and Reynolds stress models. The shortcoming of these models based on the isotropic eddy viscosity concept and Reynolds averaging in flow fields of the type found in the field of Wind Engineering are well documented. In view of these shortcomings this paper presents the implementation of a non-linear model and its evaluation for flow around a building. Tests were undertaken using the classical bluff body shape, a surface mounted cube, with orientations both normal and skewed at $45^{\circ}$ to the incident wind. Full-scale investigations have been undertaken at the Silsoe Research Institute with a 6 m surface mounted cube and a fetch of roughness height equal to 0.01 m. All tests were originally undertaken for a number of turbulence models including the standard, RNG and MMK $k-{\varepsilon}$ models and the differential stress model. The sensitivity of the CFD results to a number of solver parameters was tested. The accuracy of the turbulence model used was deduced by comparison to the full-scale predicted roof and wake recirculation zone lengths. Mean values of the predicted pressure coefficients were used to further validate the turbulence models. Preliminary comparisons have also been made with available published experimental and large eddy simulation data. Initial investigations suggested that a suitable turbulence model should be able to model the anisotropy of turbulent flow such as the Reynolds stress model whilst maintaining the ease of use and computational stability of the two equations models. Therefore development work concentrated on non-linear quadratic and cubic expansions of the Boussinesq eddy viscosity assumption. Comparisons of these with models based on an isotropic assumption are presented along with comparisons with measured data.