• Nuclear Engineering and Technology
• /
• 제55권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
• /
• 제30권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
• /
• 제29권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.

• 동력기계공학회지
• /
• 제4권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.

• 대한원격탐사학회지
• /
• 제36권6_3호
• /
• pp.1643-1652
• /
• 2020

본 연구에서는 국지예보시스템(LDAPS)과 전산유체역학(CFD) 모델을 접합하여, 부산 중구 광복동에 소재한 건물 밀집 지역의 상세 흐름을 조사하였다. 도시 지역 기상에 대한 보다 현실적인 수치 모의를 위해 기상청이 현업으로 운용하는 LDAPS 결과로부터 기상 요소(풍향, 풍속, 온위)를 추출하여 CFD 모델의 초기·경계장으로 사용하였다. LDAPS로부터 추출한 기상 요소(U, V)는 공간 선형 보간법을 이용하여 CFD 모델의 수평·연직 격자에 맞게 내삽하였고, 공간 선형 보간법을 이용해 내삽한 온위는 CFD 모델의 각 격자점에서 온도로 변환된다. 본 연구에서는 건물 옥상에서 측정한 3차원 초음파 풍향·풍속계에서 측정한 풍속과 풍향을 이용하여 수치 모의 결과를 검증하였다. 대상 기간(2020년 6월 22일) 동안 측정 지점(PKNU-SONIC)에서는 바람이 약하게 나타났고, 새벽 시간에는 북동풍과 북서풍이 불었으며, 주간에는 주로 남동풍이 불었다. LDAPS-CFD 접합 모델은 측정 풍향과 풍속을 유사하게 수치 모의하였다. 07시에는 동북동쪽에서 유입되는 흐름이 지형과 건물에 의해 변화하였고, PKNU-SONIC 지점에서 수치 모의된 풍향(동남동풍)은 측정 풍향과 유사하다. 19시에는 남동쪽에서 유입되는 흐름이 PKNU-SONIC 지점까지 유입되어 측정 풍향과 유사하게 모의하였다.

• Nuclear Engineering and Technology
• /
• 제55권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.

• 한국전산유체공학회지
• /
• 제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.

• Wind and Structures
• /
• 제23권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.

• 한국수자원학회논문집
• /
• 제55권6호
• /
• pp.393-404
• /
• 2022

위성에서 보다 미세한 공간 분해능으로 신뢰할 수 있는 지상 강우 관측은 도시 수문학적 및 미시적 농업 수요에 필수적이다. 전통적으로 "톱다운" 접근 방식 기반 위성 강우 산출물이 널리 사용되고 있지만 공간 분해능에 한계가 있다. 본 연구는 C-밴드 SAR Sentinel-1 위성 데이터(SM2RAIN-S1)에 적용되는 매개 변수화된 SM2RAIN 모델인 강우 추정을 위한 새로운 "상향식" 접근 방식의 가능성을 평가하여 중부지방에 대한 높은 공간 분해능 지상 강우 추정치(0.01° 그리드/6일)를 생성하는 것을 목표로 한다. 그것의 성능은 중부지방 두 개의 다른 하위 지역, 즉 혼합 산림 중심, 중간 하위 지역, 그리고 경작 중심, 서해안 하위 지역의 1년 기간 동안 기존의 재분석 프로덕트와 우량계 네트워크의 각각의 강우 데이터를 사용하여 공간 및 시간적 가변성에 대해 평가되었다. 평가결과에 따르면 SM2RAIN-S1 프로덕트는 중부지방의 일반적인 강우 패턴을 포착할 수 있고, 서로 다른 토지 피복으로 지역 규모에서 공간 분해능 강우량 측정 가능성을 보유할 수 있으며, 강우량 관측치에 대한 편중된 강우량 추정치가 제공되었다. 또한 SM2RAIN-S1 강우량은 피어슨의 상관 계수(R = 0.69)를 고려할 때 혼합림에서 더 우수했으며, 이는 혼합림에서 토양 수분과 강우의 시간 역학을 포착하는 데 6일 SM2RAIN-S1 데이터의 적합성을 암시했다. 그러나, RMSE와 바이어스 측면에서, 혼합림보다는 경작지의 SM2RAIN-S1 강우 생성물에서 더 나은 성능을 얻었으며, 이는 높은 증발증산 손실(특히 혼합림)에 의해 유도된 더 큰 오류를 SM2RAIN의 추가 개선에 포함해야 한다는 것을 나타낸다.