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

• Wind and Structures
• /
• v.13 no.6
• /
• pp.557-580
• /
• 2010

A new inflow turbulence generation method and a combined dynamic SGS model recently developed by the authors were applied to evaluate the wind effects on 508 m high Taipei 101 Tower. Unlike the majority of the past studies on large eddy simulation (LES) of wind effects on tall buildings, the present numerical simulations were conducted for the full-scale tall building with Reynolds number greater than $10^8$. The inflow turbulent flow field was generated based on the new method called discretizing and synthesizing of random flow generation technique (DSRFG) with a prominent feature that the generated wind velocity fluctuations satisfy any target spectrum and target profiles of turbulence intensity and turbulence integral length scale. The new dynamic SGS model takes both advantages of one-equation SGS model and a dynamic production term without test-filtering operation, which is particular suitable to relative coarse grid situations and high Reynolds number flows. The results of comparative investigations with and without generation of inflow turbulence show that: (1) proper simulation of an inflow turbulent field is essential in accurate evaluation of dynamic wind loads on a tall building and the prescribed inflow turbulence characteristics can be adequately imposed on the inflow boundary by the DSRFG method; (2) the DSRFG can generate a large number of random vortex-like patterns in oncoming flow, leading to good agreements of both mean and dynamic forces with wind tunnel test results; (3) The dynamic mechanism of the adopted SGS model behaves adequately in the present LES and its integration with the DSRFG technique can provide satisfactory predictions of the wind effects on the super-tall building.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.221-225
• /
• 2004

Evaluation of turbulence models is performed for a better prediction of thermal stratification in an upper plenum of a liquid metal reactor by applying them to the experiment conducted at JNC. The turbulence models tested in the present study are the two-layer model, the $\kappa-\omega$ model, the v2-f model and the low-Reynolds number differential stress-flux model. When the algebraic flux model or differential flux model are used for treating the turbulent heat flux, there exist little differences between turbulence models in predicting the temporal variation of temperature. However, the v2-f model and the low-Reynolds number differential stress-flux model better predict the steep gradient o( temperature at the interface of thermal stratification, and only the v2-f model predicts properly the oscillation of temperature. The LES Is needed for a better prediction of the amplitude and frequency of the temperature fluctuation.

• Journal of the Korean Society for Railway
• /
• v.7 no.3
• /
• pp.180-185
• /
• 2004

The present study investigates the effect of fire growth model on fire development characteristics in a carriage. The parallel processing version of FDS code is used to simulate the fire driven flow in a carriage and two types of fire growth model which are flame spread model and t$^2$ model are examined for the same geometrical condition. The heat release rates(HRR) of both model are similar each other until 30 s after ignition, but the flame spread model predicts 5 times higher than those of the t$^2$ fire model during the quasi-steady fire period. Maximum heat release rate in the case of flame spread model reaches about to 12 MW at 100 s after fire ignition. Also, various database of fire properties for combustible materials and more elaborate combustion model considering the flame spreading phenomena are required for better predictions of fire development characteristics using numerical simulation.

• Wind and Structures
• /
• v.13 no.1
• /
• pp.57-82
• /
• 2010

The 486m-long roof of Shenzhen Citizens Centre is one of the world's longest spatial lattice roof structures. A comprehensive numerical study of wind effects on the long-span structure is presented in this paper. The discretizing and synthesizing of random flow generation technique (DSRFG) recently proposed by two of the authors (Huang and Li 2008) was adopted to produce a spatially correlated turbulent inflow field for the simulation study. The distributions and characteristics of wind loads on the roof were numerically evaluated by Computational Fluid Dynamics (CFD) methods, in which Large Eddy Simulation (LES) and Reynolds Averaged Navier-Stokes Equations (RANS) Model were employed. The main objective of this study is to explore a useful approach for estimations of wind effects on complex curved roof by CFD techniques. In parallel with the numerical investigation, simultaneous pressure measurements on the entire roof were made in a boundary layer wind tunnel to determine mean, fluctuating and peak pressure coefficient distributions, and spectra, spatial correlation coefficients and probability characteristics of pressure fluctuations. Numerical results were then compared with these experimentally determined data for validating the numerical methods. The comparative study demonstrated that the LES integrated with the DSRFG technique could provide satisfactory prediction of wind effects on the long-span roof with complex shape, especially on separation zones along leading eaves where the worst negative wind-induced pressures commonly occur. The recommended LES and inflow turbulence generation technique as well as associated numerical treatments are useful for structural engineers to assess wind effects on a long-span roof at its design stage.

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.181.1-181.1
• /
• 2010

기후변화의 주요인이 되는 온실가스 감축을 목표로 화석연료를 대체하기 위한 대체 에너지 개발을 위한 많은 노력이 진행되고 있다. 풍력 에너지와 같은 신재생에너지는 이러한 하나의 해결 수단이 될 수 있으며 풍력 에너지 사업의 활성화를 위해서는 정확한 풍력 정보 제공이 우선이다. 풍력-기상자원지도는 풍력 발전에 유용한 정보 제공을 위한 목적으로 중규모 수치 모델을 이용하여 작성된다. 본 연구에서는 중규모 수치 모델의 정확도 향상을 위한 자료동화 방법으로써 Four-Dimensional Data Assimilation (FDDA) 방법을 이용한다. 풍력-기상자원지도는 공간분해능 1 km 해상도로 개발된다. 풍력-기상자원지도는 1998-2008년까지의 평균적인 상태에 대하여 모의를 하기 위하여 통계적인 방법으로 11년 기간의 평균과 유사한 기간을 선정하였다. 풍력-기상자원지도는 연 평균, 월 평균 풍속과 주 풍향, 주풍향 발생 비율 등의 정보를 제공한다. 우리나라 풍속의 평균 분포는 내륙 산악지역, 남해안, 제주도에서 강풍이 발생하며 주 풍향은 대체로 북서풍이다. 주 풍향의 발생비율은 산악 지역과 남 동해안에서 높아 풍력 발전의 최적지 정보를 제공한다. 1 km 해상도의 모델과 관측의 오차는 서해안 등의 해안지역보다 강원 산악지역에서 오차가 더욱 증가하였다. 이러한 산악 지역의 오차는 복잡한 지형에서는 1km 미만의 수 백 m 해상도 수치모의가 필요함을 지시한다. 따라서 본 연구에서는 WRF-LES 모형을 이용하여 333m 해상도의 기상자원지도를 개발한다. 333m 해상도의 자원지도 영역은 강원도 지역에 대하여 모의되었다. 333m 해상도의 풍속 분포는 1km 해상도의 풍속 분포와 비교하였을 때 풍속의 분포가 보다 세밀하게 표현되었다. 정량적인 검증을 하였을 때 관측소에 따라 차이는 있었으나 1km 해상도에서 과대 모의된 풍속의 분포가 현저히 개선이 되었으며, 시간적인 경향도 잘 일치함을 보였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1661-1665
• /
• 2006

In a natural river, cross sections of a channel vary according to inner or outer parts of meandering. Generally, depth of outer parts is deeper than that of inner parts. This kind of cross section change by meandering can be demonstrated by Beta distribution. The objects of this research is a 3D simulation of primary and secondary flow in the meandering natural channel. FLOW-3D program, a numerical model using CFD technique, and LES method was used for this research. 3D simulations were conducted in the channels having Beta distribution cross sections which have beds of mortar, gravel and vegetation. Two types of water stages and discharge were applied to each channel. In this research, primary flows are located in the outer parts of a top of bend and secondary flows rotate in the bottom on outer parts.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.5
• /
• pp.1-9
• /
• 2017

The cryogenic nitrogen spray of a swirl injector has been numerically investigated using three dimensional LES turbulence model to analyze the dynamic characteristics under supercritical condition. To predict the precise nitrogen properties under supercritical condition, SRK equation of state, Chung's method for viscosity and thermal conductivity and Takahashi's correlation based on Fuller's theory for diffusion coefficient are implemented. The complex flow structures due to interaction between flow field and acoustic field are observed inside and outside the injector under supercritical condition. FFT, POD, and DMD techniques are employed to understand the coherent structures. By implementing the FFT, the dominant frequencies are identified inside and outside the injector. The coherent flow structures related to the dominant frequencies are visualized using the POD and DMD techniques. In addition, the DMD provides the damping coefficient which is related with the instability prediction.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.22 no.3
• /
• pp.181-190
• /
• 2010

In case of constructing submerged breakwaters, the circulation current is occurred around the open inlet because of mean water level difference between front and rear sides of them. The aim of this study is to investigate the flow control structure of new-type submerged breakwater which is able to reduce mean water level at rear side of it. At first, the numerical model (LES-WASS-3D) is validated by comparing with existing experimental data. And then, numerical simulation is carried out to examine wave height, mean water level and mean flow around the newtype submerged breakwater. From the numerical results, it can be pointed out that the new-type submerged breakwater with drainage system reduces the rip current around the open inlet.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.449-452
• /
• 2008

Flow difference between reacting and non-reacting case in a swirl stabilized annular combustor is investigated using 3D Large Eddy Simulation with flamelet turbulent combustion model. The combustor of concern is the LM6000, lean premixed dry low-NOx annular combustor, developed by GEAE. Boundary conditions are based on experimental data. Heat release as a result of combustion put the dilatation of density in primary combustion zone highly increased so that the main swirl stream behind of a swirl cup stretched further downstream than that of non-reacting case. The oval shape of core flow in cross-section to flow direction, which clearly observed in non-reacting case, tends to be circle, and small vorticities in wide range in non-reacting case disappears, but the size of iso-vorticity increase in reacting case.