• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.5
• /
• pp.1435-1442
• /
• 2014

This study presented numerical simulations of smooth-bed flows in the rectangular open-channel using the source code by OpenFOAM. For the analysis of the turbulent flow, Large Eddy Simulations were carried out and the dynamic sub-grid scale model proposed by Germano et al. (1991) is used to model the residual stress term. In order to analyze the coherent structure, the uw quadrant method proposed by Lu and Willmarth (1973) is used and the contribution rate and the fraction time of the instantaneous Reynolds stress are obtained in the Reynolds stress. The results by the present study are analyzed and compared with data from previous laboratory studies and direct numerical simulations. It is found that the contribution rate of the ejection events is larger than that of sweep events over the buffer layer in the open-channel flow over the smooth bed, however, the frequency of the sweep event is higher than that of the ejection events.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.236-245
• /
• 2021

In this study, direct numerical and large eddy simulations of transitional flows around studs were conducted to investigate the effectiveness of turbulence stimulators at very low speeds for the minimum propulsion power condition of four knots. For simplicity, the studs were assumed to be installed on a flat plate, while the wake was observed up to 0.23 m downstream behind the second stud. For applicability to a model ship, we also studied the flow characteristics behind the first and second studs installed on a curved plate, which was designed to describe the geometry of a bulbous bow. A laminar-to-turbulent transition was observed in the wake at Re_D ≥ 921 (U_∞≥0.290 m/s), and the wall shear stress at Re_D = 1162 (U_∞ = 0.366 m/s) in the second wake was similar to that of the fully developed turbulent boundary layer after a laminar-to-turbulent transition in the first wake. At Re_D = 581 (U_∞ = 0.183 m/s), no turbulence was stimulated in the wake behind the first and second studs on the flat plate, while a cluster of vortical structures was observed in the first wake over the curved plate. However, a cluster of vortical structures was revealed to be generated by the reattachment process of the separated shear layer, which was disturbed by the first stud rather than directly initiated by the first stud. It was quite different from a typical process of transition, which was observed at relatively high Re_D that the spanwise scope of the turbulent vortical structures expanded gradually as it went downstream.

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.1594-1600
• /
• 2011

Recently developed FDS5 CFD code has employed a near-wall flow treatment method which is Werner-Wengle wall law provided by NIST(National Institute of Standards and Technology). In this study, the wall law has been verified against DNS(Direct Numerical Simulation) data in the parallel plate. The $y^+$ was kept above 11 to fulfill the near-wall flow requirement in the grid generation. The total grid was $32{\times}32{\times}32$. The boundary condition for inlet and outlet was periodic condition and for both side, symmetric condition was used. The fully developed turbulent flow was generated and Re = 10,700. The simulated results were compared with DNS data. RANS results were also used for verification.

• Wind and Structures
• /
• v.8 no.6
• /
• pp.427-442
• /
• 2005

Numerical and wind tunnel simulations of pollutant dispersion around rectangular obstacles with five aspect ratios have been conducted in order to identify the effects of flow patterns induced by buildings on plume dispersion in the near wake of buildings. An emission from a low source located upwind of obstacles was used in this simulation. The local flow patterns and concentrations around a cubical obstacle were initially investigated using three RANS turbulence models, (the standard $k-{\varepsilon}$, Shear Stress Transport (SST), Reynolds-Stress RSM turbulence model) and also using Large-eddy simulation (LES). The computed concentrations were compared with those measured in the wind tunnel. Among the three turbulence models, the SST model offered the best performance and thus was used in further investigations. The results show, for normal aspect ratios of width to height, that concentrations in the near wake are appreciably affected because of plume capture by the horseshoe vortex and convection by the vertical vortex pairs. These effects are less important for high aspect ratios. Vertical vortex pairs present a strong ability to exchange mass vertically and acts efficiently to reduce ground-level concentrations in the near wake.

• The Transactions of the Korea Information Processing Society
• /
• v.7 no.1
• /
• pp.57-64
• /
• 2000

We describe a parallel implementation for the large-eddy simulation(LES) of stratified and rotating turbulence based on MPI. The parallelization strategy is specified by eliminating the tridiagonal solver with explicit method and by domain decompositions for solving the poisson equation. In this simulation we have run on CRAY-T3E under the message passing platform MPI with a various domain decomposition and the scalability of this parallel code of LES are also presented. The result shows that we can gain up to 16 times faster speed up on 64 processors with xyz-directional domain decomposition and scalable up to $128{\times}128{\times}$ which processing time is almost similar to that of $40{\times}40{\times}40$ on a single processor machine with a sequential code.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3874-3897
• /
• 2023

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High-Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next-generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds-averaged Navier-Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

• The Journal of the Acoustical Society of Korea
• /
• v.36 no.5
• /
• pp.314-320
• /
• 2017

A pressure relief valve is generally used to prevent piping systems from being broken due to high pressure gas flows. However, the sudden pressure drop caused by the pressure relief valve produces high acoustic energy which propagates in the form of compressible acoustic waves in the pipe and sometimes causes severe vibration of the pipe structure, thereby resulting in its failure. In this study, internal aerodynamic noise due to valve flow is estimated for a simple contraction-expansion pipe by combining the LES (Large-Eddy Simulation) technique with the wavenumber-frequency analysis, which allows the decomposition of fluctuating pressure into incompressible hydrodynamic pressure and compressible acoustic pressure. In order to increase the convergence, the steady Reynolds-Averaged Navier-Stokes equations are numerically solved. And then, for the unsteady flow analysis with high accuracy, the unsteady LES is performed with the steady result as the initial value. The wavenumber-frequency analysis is finally performed using the unsteady flow simulation results. The wavenumber-frequency analysis is shown to separate the compressible pressure fluctuation in the flow field from the incompressible one. This result can provide the accurate information for the source causing so-called acoustic-induced-vibration of a piping system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of the Korean Society for Railway
• /
• v.17 no.5
• /
• pp.328-335
• /
• 2014

Numerical prediction methods were applied to investigate the turbulent air currents and air-conditioning efficiency in an underground subway station, and the results compared to experimental data. The Shin-gumho Station($8^{th}$ floor underground and 43.6m in depth) in Seoul was selected for the analysis. The entire station was covered for simulation and the ventilation mode was ordinary. The ventilation diffusers were modeled as 95 square shapes of $0.6m{\times}0.6m$ in the lobby and as 222 square shapes in the platform. Cooling air of $47,316m^3/h$ was supplied and the returned air of $33,980m^3/h$ is exhausted in the lobby and the cooling air of $33,968m^3/h$ is supplied and the returned air of $76,190m^3/h$ was exhausted in the platform which is the same as the experimental data. The cases of the screen-door-closed and open were respectively investigated. A total of 7.5million grids were generated and the whole domain divided into 22 blocks for MPI efficiency of calculation. Large eddy simulation (LES) was applied to solve the momentum and energy equation.

• Wind and Structures
• /
• v.28 no.5
• /
• pp.299-313
• /
• 2019

In order to examine the effects of different wind deflectors on the wind load distribution characteristics of extra-large cooling towers, a comparative study of the distribution characteristics of wind pressures on the surface of three large cooling towers with typical wind deflectors and one tower without wind deflector was conducted using wind tunnel tests. These characteristics include aerodynamic parameters such as mean wind pressures, fluctuating wind pressures, peak factors, correlation coefficients, extreme wind pressures, drag coefficients and vorticity distribution. Then distribution regularities of different wind deflectors on global and local wind pressure of extra-large cooling towers was extracted, and finally the fitting formula of extreme wind pressure of the cooling towers with different wind deflectors was provided. The results showed that the large eddy simulation (LES) method used in this article could be used to accurately simulate wind loads of such extra-large cooling towers. The three typical wind deflectors could effectively reduce the average wind pressure of the negative pressure extreme regions in the central part of the tower, and were also effective in reducing the root of the variance of the fluctuating wind pressure in the upper-middle part of the windward side of the tower, with the curved air deflector showing particularly. All the different wind deflectors effectively reduced the wind pressure extremes of the middle and lower regions of the windward side of the tower and of the negative pressure extremes region, with the best effect occurring in the curved wind deflector. After the wind deflectors were installed the drag coefficient values of each layer of the middle and lower parts of the tower were significantly higher than that without wind deflector, but the effect on the drag coefficients of layers above the throat was weak. The peak factors for the windward side, the side and leeward side of the extra-large cooling towers with different wind deflectors were set as 3.29, 3.41 and 3.50, respectively.