• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.3
• /
• pp.369-375
• /
• 2006

Large eddy simulation has been applied to simulate turbulent flow around a cubic obstacle mounted on a channel surface for a Reynolds number of 40000(based on the incoming bulk velocity and the obstacle height) using a Smagorinsky model and a Lagrangian dynamic model. In order to develop the LES to the practical engineering application, the effect of upwind scheme, turbulent sub-grid scale model were investigated. The computed velocities. turbulence quantifies, separation and reattachment length were evaluated by compared with the previous experimental results.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.8
• /
• pp.849-854
• /
• 2005

The laminar flame concept in turbulent reacting flow is considered applicable to many practical combustion systems For turbulent premixed combustion under widely used flamelet concept, the flame surface is described as an infinitely thin propagating surface that such a Propagating front can be represented as a level contour of a continuous function G. In this study, for the Purpose of validating the LES of G-equation combustion model. LES of turbulent Premixed combustion with dynamic SGS model of G-equation in turbulent channel flow are carried out A constant density assumption is used. The Predicted flame propagating speed is goof agreement with the DNS result of G. Bruneaux et al.

• Journal of Power System Engineering
• /
• v.4 no.3
• /
• pp.40-47
• /
• 2000

The turbulent flow with wake, reattachment and recirculation flow is very important from the viewpoint of engineering. But that is still difficult because of especially the unsteady problems which are related with the vehicle dynamics and the aerodynamics noise. This paper evaluate LES that can analyze about all fluid flow region including the laminar, transition and turbulent. So we compare the results of LES with those of PIV measurement and Reynolds averaging models. In conclusion, LES predicts flow behavior better than Reynolds averaging models.

• Korean Journal of Remote Sensing
• /
• v.36 no.6_3
• /
• pp.1643-1652
• /
• 2020

To investigate the characteristics of detailed flows in a building-congested district, we coupled a computation fluid dynamics (CFD) model to the local data assimilation and prediction system (LDAPS), a current operational numerical weather prediction model of the Korea Meteorological Administration. For realistic numerical simulations, we used the meteorological variables such as wind speeds and directions and potential temperatures predicted by LDAPS as the initial and boundary conditions of the CFD model. We trilinearly interpolated the horizontal wind components of LDAPS to provide the initial and boudnary wind velocities to the CFD model. The trilinearly interpolated potential temperatures of LDAPS is converted to temperatures at each grid point of the CFD model. We linearly interpolated the horizontal wind components of LDAPS to provide the initial and boundary wind velocities to the CFD model. The linearly interpolated potential temperatures of LDAPS are converted to temperatures at each grid point of the CFD model. We validated the simulated wind speeds and directions against those measured at the PKNU-SONIC station. The LDAPS-CFD model reproduced similar wind directions and wind speeds measured at the PKNU-SONIC station. At 07 LST on 22 June 2020, the inflow was east-north-easterly. Flow distortion by buildings resulted in the east-south-easterly at the PKNU-SONIC station, which was the similar wind direction to the measured one. At 19 LST when the inflow was southeasterly, the LDAPS-CFD model simulated southeasterly (similar to the measured wind direction) at the PKNU-SONIC station.

• 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.

• Journal of computational fluids engineering
• /
• v.20 no.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.

• Wind and Structures
• /
• v.23 no.4
• /
• pp.313-350
• /
• 2016

Taipei 101 Tower, which has 101 stories with height of 508 m, is located in Taipei where typhoons and earthquakes commonly occur. It is currently the second tallest building in the world. Therefore, the dynamic performance of the super-tall building under strong wind actions requires particular attentions. In this study, Large Eddy Simulation (LES) integrated with a new inflow turbulence generator and a new sub-grid scale (SGS) model was conducted to simulate the wind loads on the super-tall building. Three-dimensional finite element model of Taipei 101 Tower was established and used to evaluate the wind-induced responses of the high-rise structure based on the simulated wind forces. The numerical results were found to be consistent with those measured from a vibration monitoring system installed in the building. Furthermore, the equivalent static wind loads on the building, which were computed by the time-domain and frequency-domain analysis, respectively, were in satisfactory agreement with available wind tunnel testing results. It has been demonstrated through the validation studies that the numerical framework presented in this paper, including the recommended SGS model, the inflow turbulence generation technique and associated numerical treatments, is a useful tool for evaluation of the wind loads and wind-induced responses of tall buildings.

• Journal of Korea Water Resources Association
• /
• v.55 no.6
• /
• pp.393-404
• /
• 2022

Reliable terrestrial rainfall observations from satellites at finer spatial resolution are essential for urban hydrological and microscale agricultural demands. Although various traditional "top-down" approach-based satellite rainfall products were widely used, they are limited in spatial resolution. This study aims to assess the potential of a novel "bottom-up" approach for rainfall estimation, the parameterized SM2RAIN model, applied to the C-band SAR Sentinel-1 satellite data (SM2RAIN-S1), to generate high-spatial-resolution terrestrial rainfall estimates (0.01° grid/6-day) over Central South Korea. Its performance was evaluated for both spatial and temporal variability using the respective rainfall data from a conventional reanalysis product and rain gauge network for a 1-year period over two different sub-regions in Central South Korea-the mixed forest-dominated, middle sub-region and cropland-dominated, west coast sub-region. Evaluation results indicated that the SM2RAIN-S1 product can capture general rainfall patterns in Central South Korea, and hold potential for high-spatial-resolution rainfall measurement over the local scale with different land covers, while less biased rainfall estimates against rain gauge observations were provided. Moreover, the SM2RAIN-S1 rainfall product was better in mixed forests considering the Pearson's correlation coefficient (R = 0.69), implying the suitability of 6-day SM2RAIN-S1 data in capturing the temporal dynamics of soil moisture and rainfall in mixed forests. However, in terms of RMSE and Bias, better performance was obtained with the SM2RAIN-S1 rainfall product over croplands rather than mixed forests, indicating that larger errors induced by high evapotranspiration losses (especially in mixed forests) need to be included in further improvement of the SM2RAIN.