Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
安部 諭; 柴本 泰照
Annals of Nuclear Energy, 202, p.110461_1 - 110461_16, 2024/07
During a severe accident in a nuclear containment vessel, jets released from the primary system exhibit complex thermohydraulic behavior due to buoyancy effects and impingement on internal obstacles such as inner walls and floors. Thus, the obstacle-influenced jets are of interest in recent research activities. This paper describes an experimental investigation of the behavior of jets passing through a grid-type obstacle. The flow field was acquired by a particle image velocimetry system. The experiment captured the jet fragmentation by the grid-type obstacle and their recoupling. The mean velocity field obtained by postprocessing indicates a "Rectifying effect," with the axial velocity increasing at the center and the magnitude of the radial velocity decreasing. The meandering flow was suppressed due to this effect. In the near grid-obstacle region, the axial turbulence intensity was relatively large at the edge of each fragmented region due to shear stress. Moreover, the spatial distribution of the radial turbulence fluctuation became more complex. Further investigation is required to clarify the budget of the transport equation for turbulence fluctuation. The experimental data shown in this paper is useful for computational fluid dynamics validation.
廣瀬 意育; 安部 諭; 石垣 将宏*; 柴本 泰照; 日引 俊*
Progress in Nuclear Energy, 169, p.105085_1 - 105085_13, 2024/04
被引用回数：0Immersed boundary methods (IBMs) have been developed as complementary methods for computational fluid dynamics (CFD). They allow a flow simulation in a mock-up model that includes complex-shaped inner structures and/or boundaries with a non-body conformal mesh. Such a model might force us to create a complicated body-fitted mesh with a high cost in the conventional CFD (CCFD) approach. We focus on the Brinkman penalization (BP) method and its extended version, which we call here the extended Brinkman penalization method (EBP), among the different types of IBMs, aiming to apply them to the phenomena that occur during severe accidents in a nuclear reactor containment vessel and explore the possibility that the methods can partially replace the CCFD. In this paper, as a preliminary step to validate the applicability of these methods, we measure the jet flow rectified by a grating-type structure used for the validation of numerical techniques and apply them to simulate the behavior of an upward jet rectified by a horizontally placed grating-type structure modeled as an immersed body. This type of structure is generally used in reactor buildings, and it is crucial to evaluate their influence on gaseous flows because the behaviors of hydrogen produced during severe accidents may be influenced by them. The structure is selected as our subject because it has moderate complexity, enabling us to examine the effects of the IBMs and compare them with CCFD. We investigate whether these methods can reproduce a result of corresponding CCFD in which the grating is modeled as body-conformal mesh and show that the former can produce the latter with equivalent accuracy. All these results are also compared with the experimental data on the flow velocity distributions downstream of the grating measured using particle image velocimetry.
相馬 秀; 石垣 将宏*; 安部 諭; 柴本 泰照
Nuclear Engineering and Design, 416, p.112754_1 - 112754_18, 2024/01
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)The wall function (WF) enables analyzing condensation flow in a nuclear reactor containment vessel with reasonable computational costs. However, conventional wall treatments rely on the logarithmic laws for velocity, temperature, and concentration, limiting applicability. In this paper, we applied the analytical wall function approach to the condensation flow analysis of steam/air mixtures. This approach features the analytical integration of transport equations considering the buoyancy, the material property change, and the convective terms. We conducted CFD analysis with the analytical wall function models for the forced, mixed, and natural convection and confirmed good prediction, especially when the log law does not hold.
相馬 秀; 石垣 将宏*; 安部 諭; 柴本 泰照
Nuclear Engineering and Technology, 10 Pages, 2024/00
When analyzing containment thermal-hydraulics, computational fluid dynamics (CFD) is a powerful tool because multi-dimensional and local analysis is required for some accident scenarios. According to the previous study, neglecting steam bulk condensation in the CFD analysis leads to a significant error in boundary layer profiles. Validating the condensation model requires the experimental data near the condensing surface, however, available boundary layer data is quite limited. It is also important to confirm whether the heat and mass transfer analogy (HMTA) is still valid in the presence of bulk condensation. In this study, the boundary layer measurements on the vertical condensing surface in the presence of air were performed with the rectangular channel facility WINCS, which was designed to measure the velocity, temperature, and concentration boundary layers. We set the laminar flow condition and varied the Richardson number (1.0-23) and the steam volume fraction (0.35-0.57). The experimental results were used to validate CFD analysis and HMTA models. For the former, we implemented a bulk condensation model assuming local thermal equilibrium into the CFD code and confirmed its validity. For the latter, we validated the HMTA-based correlations, confirming that the mixed convection correlation reasonably predicted the sum of wall and bulk condensation rates.
岡垣 百合亜; 柴本 泰照; 和田 裕貴; 安部 諭; 日引 俊詞*
Journal of Nuclear Science and Technology, 60(8), p.955 - 968, 2023/08
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)Pool scrubbing is an important filtering process that prevents radioactive aerosols from entering the environment in the event of severe accidents in a nuclear reactor. In this process of transporting aerosol particles using bubbles, bubble hydrodynamics plays a crucial role in modeling pool scrubbing and significantly affects particle removal in a bubble. The pool scrubbing code based on Lumped Parameter (LP) approach includes the particle removal model, and its hydrodynamic parameters are determined based on simple assumptions. We aim to apply the three-dimensional Computer Fluid Dynamics (CFD) approach to understand the detailed bubble interaction. This study validated the applicability of the CFD simulation to bubble hydrodynamics at the flow transition from a globule to a swarm region, which is critical in the stand-alone pool scrubbing code-SPARC-90. Two types of solvers based on the Volume Of Fluid (VOF) and the Simple Coupled Volume Of Fluid with Level Set (S-CLSVOF) methods were used to capture the gas-liquid interface in the CFD simulation. We used the experimental data for validation. As a result, the VOF and S-CLSVOF methods accurately predicted the bubble size and void fraction distributions. In addition, we confirmed that the bubble rise velocity of the S-CLSVOF method almost agreed with the experimental results.
安部 諭; 柴本 泰照
Nuclear Engineering and Technology, 55(5), p.1742 - 1756, 2023/05
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)The hydrogen behavior in a nuclear containment vessel is a significant issue when discussing the potential of hydrogen combustion during a severe accident. After the Fukushima-Daiichi accident in Japan, we have investigated in-depth the hydrogen transport mechanisms by utilizing experimental and numerical approaches. Computational fluid dynamics is a powerful tool for better understanding the transport behavior of gas mixtures, including hydrogen. This paper describes a large-eddy simulation of gas mixing driven by a high-buoyancy flow. We focused on the interaction behavior of heat and mass transfers driven by the horizontal high-buoyant flow during density stratification. For validation, the experimental data of the Containment InteGral effects Measurement Apparatus (CIGMA) facility were used. With a high-power heater for the gas-injection line in the CIGMA facility, a high temperature flow of approximately 390
C was injected into the test vessel. By using the CIGMA facility, we can extend the experimental data to the high temperature region. The phenomenological discussion in this paper help understand the heat and mass transfer induced by the high-buoyancy flow in the containment vessel during a severe accident.
安部 諭; 岡垣 百合亜
Nuclear Engineering and Design, 404, p.112165_1 - 112165_14, 2023/04
被引用回数：0 パーセンタイル：0.01(Nuclear Science & Technology)Pressurized Thermal Shock (PTS) is induced potentially by the rapid cooling of the cold-leg and downcomer wall in the primary system of a Pressurized Water Reactor (PWR) due to the initiation of Emergency Core Cooling System (ECCS). Thus, fluids mixing in a horizontal cold-leg and downcomer should be predicted accurately; however, turbulence production and damping often hinders this prediction due to the presence of the density gradients. Hence, the Fifth International Benchmark Exercise, the cold-leg mixing Computational Fluid Dynamics (CFD) Benchmark, was conducted under the support of OECD/NEA. The experiment was designed for visualization of the mixing phenomena of two liquids with different densities. The heavy liquid was a simulant of cold water from ECCS, in a horizontal leg and downcomer. We used the Large-eddy Simulation (LES) to investigate the time fluctuation behaviors of velocity and liquid concentration. The CFD simulation was performed with two turbulence models and three different numerical meshes. We investigated the characteristics of the appearance frequency of the heavy liquid concentration with the statistical method. Based on our findings, we propose further experiments and numerical investigations to understand the fluid mixing phenomena related to PTS.
Hamdani, A.; 安部 諭; 石垣 将宏; 柴本 泰照; 与能本 泰介
Progress in Nuclear Energy, 153, p.104415_1 - 104415_16, 2022/11
被引用回数：3 パーセンタイル：68.71(Nuclear Science & Technology)This paper describes the computational fluid dynamics (CFD) analysis and validation works from the previous experimental study on the natural convection driven by outer surface cooling in the presence of density stratification consisting of air and helium (as a mimic gas of hydrogen). The experiment was conducted in the Containment InteGral effects Measurement Apparatus (CIGMA) facility at Japan Atomic Energy Agency (JAEA). The numerical simulation was carried out to analyze the detailed effect of the cooling region on the erosion of the helium stratification layer. The temporal and spatial evolution of the helium concentration and the gas temperature inside the containment vessel was predicted and validated against the experimental data. In addition, two stratification behaviors that depend on the cooling location were presented and discussed. The CFD simulation confirmed that an upper head cooling caused two counter-rotating vortexes in the helium-rich zone. Meanwhile, the upper half body cooling caused two counter-rotating vortexes in the helium-poor zone. These findings are important to understand the mechanism of the density stratification process driven by natural convection in the containment vessel.
石垣 将宏*; 廣瀬 意育; 安部 諭; 永井 亨*; 渡辺 正*
Fluids (Internet), 7(7), p.237_1 - 237_18, 2022/07
To estimate thermal flow in a nuclear reactor during an accident, it is important to improve the accuracy of computational fluid dynamics simulation. Temperature and flow velocity are not homogeneous and have large variations in a reactor containment vessel because of its very large volume. In addition, Kelm et al (2016) pointed out that the influence of variations of initial and boundary conditions was important. Therefore, it is necessary to set the initial and boundary conditions taking into account the variations of these physical quantities. However, it is a difficult subject to set such complicated initial and boundary conditions. Then, we can obtain realistic initial and boundary conditions by the data assimilation technique, and we can improve the accuracy of the simulation result. In this study, we applied the data assimilation by local ensemble transform Kalman filter (Hunt et al., 2007) to the simulation of natural convection behavior in density stratification, and we performed a twin model experiment. We succeeded in the estimation of the flow fields and improving the simulation accuracy by the data assimilation, even if we applied the boundary condition with error for the true condition.
石垣 将宏*; 安部 諭; Hamdani, A.; 廣瀬 意育
Annals of Nuclear Energy, 168, p.108867_1 - 108867_20, 2022/04
被引用回数：4 パーセンタイル：68.71(Nuclear Science & Technology)It is essential to improve computational fluid dynamics (CFD) analysis accuracy to estimate thermal flow in a containment vessel during a severe accident. Previous studies pointed out the importance of the influence of initial and boundary conditions on the CFD analysis. The purpose of this study is to evaluate the influence of initial and boundary conditions by numerical analysis of natural convection experiments in a large containment vessel test facility CIGMA(Containment InteGral effects Measurement Apparatus). A density stratification layer was initially formed in the vessel using helium and air, and external cooling of the vessel surface-induced natural convection. In this study, we carried out numerical simulations of the density stratification erosion driven by the natural convection using the RANS (Reynolds averaged Navier-Stokes) model. As a result, the temperature boundary condition of the small internal structure in the vessel had a significant influence on the fluid temperature distribution in the vessel. The erosion velocity of the density stratification layer changed depending on the initial gas concentration distribution. Then, appropriate settings of the temperature and gas concentration conditions are necessary for accurate analysis.
相馬 秀; 安部 諭; 柴本 泰照; 石垣 将宏*
Proceedings of 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19) (Internet), 13 Pages, 2022/03
The experimental data of boundary layer profiles are necessary to validate condensation models applied in numerical simulation of CFD codes and also to develop wall treatment models for heat and mass transfer in the presence of significant buoyancy and suction. The available data for velocity, temperature, and concentration boundary layer, however, is quite limited. In this study, we present experimental results of measuring boundary layer profiles by our experimental facility WINCS (WINd tunnel for Condensation of Steam and air mixture). WINCS is a once-through type of wind tunnel having a 1.5m-long measuring part of rectangular duct. The velocity and temperature profiles were obtained with laser doppler velocimetry and thermocouple, respectively. The temperature data was then used to calculate the steam concentration boundary layer by assuming local thermodynamic equilibrium. The laminar boundary layer profiles were obtained in the present condition. The dropwise condensation and fimwise condensation are available in this apparatus. We also conducted numerical simulations with CFD codes and compared the experimental and numerical results of boundary layer profiles and heat fluxes. The wall condensation model based on Stefan flow and bulk condensation model were used in the numerical analysis. The comparison, in general, shows good agreement between the experimental and numerical results.
安部 諭; 小尾 善男*; 佐藤 聡; 岡垣 百合亜; 柴本 泰照
Proceedings of 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19) (Internet), 13 Pages, 2022/03
A modeling of heat transfers with boiling transition (BT) and that after the occurrence of BT, called post-BT, is one of the key technical issues to estimate the duration of surface dryout and the peak cladding temperature during DBA (Design Basis Accident) and BDBA (Beyond Design Basis Accident) in light water reactors. Recently, CFD (Computational Fluid Dynamics) has emerged as a powerful tool for representing the heat transfer mechanism. Our main purpose is to obtain in-depth physical insight into the BT and post-BT phenomena by combining experiment and CFD simulation research. This paper introduces our developing activity for an integrated three-field two-fluid CFD methodology based on the Eulerian-Eulerian approach toward the accurate prediction of the dryout behavior from annular-mist to mist flow regimes. We implemented following interaction terms and functions into OpenFOAM ver.7, an opensource code developed by OpenFOAM foundation, as (1) Interaction terms between liquid film and droplets due to deposition and entrainment, (2) Interaction terms on the liquid interface between the liquid film and the gas phase on friction and heat conduction, (3) Heat transfer from the heated wall to the liquid film, (4) Dryout occurrence judgement and the switching function on the boundary condition. The dryout occurrence judgment is based on a correlation on critical film thickness, which is originally applied into the MARS (Multi-dimensional Analysis of Reactor Safety) code. A trial calculation with the developed solver called two Phase Three Field Euler Foam was performed to check the solver operation. The CFD could simulate temperature increase behavior due to the dryout occurrence, whereas here were still challenges in reproducing the transition from mist flow to single-phase vapor flow.
安部 諭; Hamdani, A.; 石垣 将宏*; 柴本 泰照
Annals of Nuclear Energy, 166, p.108791_1 - 108791_18, 2022/02
被引用回数：5 パーセンタイル：56.94(Nuclear Science & Technology)This paper describes an experimental investigation of natural convection driven by outer surface cooling in the presence of density stratification consisting of an air-helium gas mixture (as mimic gas of hydrogen) in an enclosed vessel. The unique cooling system of the Containment InteGral effects Measurement Apparatus (whose test vessel is a cylinder with 2.5-m diameter and 11-m height) is used, and findings reveal that the cooling location relative to the stratification plays an important role in determining the interaction behavior of the heat and mass transfer in the enclosed vessel. When the cooling region is narrower than the stratification thickness, the density-stratified region expands to the lower part while decreasing in concentration (stratification dissolution). When the cooling region is wider than the stratification thickness, the stratification is gradually eroded from the bottom with decreasing layer thickness (stratification breakup). This knowledge is useful for understanding the interaction behavior of heat and mass transfer during severe accidents in nuclear power plants.
安部 諭; 岡垣 百合亜; 佐藤 聡; 柴本 泰照
Annals of Nuclear Energy, 159, p.108321_1 - 108321_12, 2021/09
被引用回数：3 パーセンタイル：45.99(Nuclear Science & Technology)Heat transfer in the rod bundle is augmented by the mixing vanes on the spacer grid. We conducted a computational fluid dynamics (CFD) simulation with three isotropic turbulence models - standard 
-epsilon, realizable -epsilon, and SST -omega models - to investigate the relationship between heat transfer and turbulence behaviors downstream a simulant spacer (with four vanes) in a single tube under single-phase flow conditions. Quantitatively, the predicted heat transfer coefficient (HTC), secondary flow intensity, and turbulence intensity with the SST -omega model displayed a better agreement (than the other isotropic models) with the correlation based on previous studies. Furthermore, the turbulence production was localized in the near-spacer region (z/D 
10, where D is the inner diameter), which corresponds with the HTC argumentation region. These results indicate that examining the turbulence production when discussing the HTC augmentation downstream the spacer is essential.
安部 諭; Studer, E.*; 石垣 将宏; 柴本 泰照; 与能本 泰介
Nuclear Engineering and Design, 368, p.110785_1 - 110785_14, 2020/11
被引用回数：10 パーセンタイル：75.92(Nuclear Science & Technology)The hydrogen behavior in a nuclear containment vessel is one of the significant issues raised when discussing the potential of hydrogen combustion during a severe accident. Computational Fluid Dynamics (CFD) is a powerful tool for better understanding the turbulence transport behavior of a gas mixture, including hydrogen. Reynolds-averaged Navier-Stokes (RANS) is a practical-use approach for simulating the averaged gaseous behavior in a large and complicated geometry, such as a nuclear containment vessel; however, some improvements are required. We implemented the dynamic modeling for 
based on the previous studies into the OpenFOAM ver 2.3.1 package. The experimental data obtained by using a small scale test apparatus at Japan Atomic Energy Agency (JAEA) was used to validate the RANS methodology. Moreover, Large-Eddy Simulation (LES) was performed to phenomenologically discuss the interaction behavior. The comparison study indicated that the turbulence production ratio by shear stress and buoyancy force predicted by the RANS with the dynamic modeling for was a better agreement with the LES result, and the gradual decay of the turbulence fluctuation in the stratification was predicted accurately. The time transient of the helium molar fraction in the case with the dynamic modeling was very closed to the VIMES experimental data. The improvement on the RANS accuracy was produced by the accurate prediction of the turbulent mixing region, which was explained with the turbulent helium mass flux in the interaction region. Moreover, the parametric study on the jet velocity indicates the good performance of the RANS with the dynamic modeling for on the slower erosive process. This study concludes that the dynamic modeling for is a useful and practical approach to improve the prediction accuracy.
安部 諭; 岡垣 百合亜; 石垣 将宏; 柴本 泰照
Proceedings of OECD/NEA Workshop on Virtual CFD4NRS-8; Computational Fluid Dynamics for Nuclear Reactor Safety (Internet), 11 Pages, 2020/11
The fifth international benchmark exercise (IBE-5), the cold-leg mixing CFD benchmark, was conducted under the support of OECD/NEA. The experiment for IBE-5 was designed to visualize the mixing phenomena of two liquids with different density in a horizontal leg (as a simulant of the cold-leg) and downcomer. This paper shows our CFD result on the open test condition in IBE-5. We selected the Large-eddy simulation (LES) solving the filtered equation of flow and concentration fields. Regarding the eddy-viscosity to model the turbulence flux of the momentum at sub-grid scale (SGS), Wall-adapting locale eddy-viscosity (WALE) model, a modified version from the Smagorinsky model, was applied. The experimental geometry was resolved with three different numerical mesh systems. The CFD analysis predicted the laminar-like flow behavior in the horizontal leg. Due to the large density difference between the two liquids, the turbulence production was suppressed strongly, and the velocity fluctuation in the horizontal leg became very slow and small. In contrast, the strong turbulence mixing in the downcomer was predicted. The plume from the horizontal leg entrained with the surroundings and spread circumferentially in the downcomer. The comparison with the TAMU experimental data reveals the good performance of the WALE model. In addition, we discuss the appearance characteristics of the high concentration of the heavy liquid in the downcomer in the LES. The Probability Density Function (PDF) and Cumulative Distribution Function (CDF) are derived based on the predicted time-series of the heavy liquid concentration. The PDF around the mean concentration in the case with the low mesh resolution is larger than that predicted by the higher resolution due to the excessive homogenization of the heavy fluid concentration. This study reveals the importance to note the required mesh resolution to predict the appearance event of the high concentration.
岡垣 百合亜; 柴本 泰照; 安部 諭
Proceedings of OECD/NEA Workshop on Virtual CFD4NRS-8; Computational Fluid Dynamics for Nuclear Reactor Safety (Internet), 12 Pages, 2020/11
A bubbly flow with a single injection orifice is numerically analyzed for pool scrubbing phenomena using different computational fluid dynamics (CFD) methodologies. The calculation covers the total regime of pool scrubbing from air injection to bubble swarm through the transition region. Such two-phase flow behaviors strongly affect particle removal in a bubble. The experimental bubbles are known to be oblate spherical and exhibit secondary motion, including path instability and shape oscillations. Moreover, bubbles in a swarm are subject to coalescence and breakup. While these may well affect bubble internal heat/mass transfer and particle capture, no established way is available for considering such influences in practical calculations. Pool scrubbing code SPARC-90 uses an oblate spherical bubble model but assumes a steady, rectilinear bubble rise without secondary motion. The 3-D CFD has the potential to capture the bubble interaction in the swarm region in detail. In the present study, the experiment by Abe et al. (Nuclear Engineering and Design 337, 2018) was referred for the calculation, and their data were used to validate if the CFD simulation can predict the flow transition accurately. Two types of solvers based on the volume of fluid (VOF) method and the simple coupled volume of fluid with level set (S-CLSVOF) method are used for the interface capture. The two solvers were validated by comparing with the experimental results. As a result, the void fraction profiles along the vertical central axis were in good agreement with the experimental data, regardless of the solvers, and those along horizontal lines in a central plane slightly improved with the S-CLSVOF method by the more accurate calculation of the surface tension.
石垣 将宏; 安部 諭; 柴本 泰照; 与能本 泰介
Nuclear Engineering and Design, 367, p.110790_1 - 110790_15, 2020/10
被引用回数：14 パーセンタイル：84.32(Nuclear Science & Technology)シビアアクシデント時の格納容器(CV)内の流体や構造物を冷却するための有効なアクシデントマネジメントとして、CVの外面冷却が期待されている。一方、以下のような可能性も考えられる。第一に、シビアアクシデント時に水-ジルコニウム反応により水素ガスが発生し、外表面冷却により水素を含む非凝縮性ガスが蓄積し、密度成層が形成される可能性がある。第二に、非凝縮性ガスの蓄積は熱伝達を低下させ、CVの冷却を阻害する可能性がある。これらの課題については、これまで多くの研究が行われてきた。しかし、外表面冷却によって生じる密度成層挙動や成層崩壊の条件に着目した体系的な検討は十分に行われていない。また、水素の蓄積による伝熱劣化を定量的に評価していない。そこで、実験設備CIGMA(Containment InteGral effects Measurement Apparatus)を構築し、CIGMA設備を用いて容器外面冷却時の格納容器熱流動挙動の実験的研究を行った。さらに、安定な密度成層が維持できる条件を考慮することで、自然対流が密度成層化挙動に与える影響を議論した。
廣瀬 意育; 石垣 将宏; 安部 諭; 柴本 泰照
Proceedings of International Topical Meeting on Advances in Thermal Hydraulics (ATH '20) (Internet), p.757 - 767, 2020/10
This paper describes an application of an immersed boundary method (IBM) to internal structure in a containment vessel as a means to reduce the modeling cost for accident analysis. A brinkman penalization method, that is one of the IBMs, is selected and implemented in OpenFOAM code. The calculation is performed for a grating-type structure in a rectangular vessel and an upward jet flow impinging on the grating is analyzed. The IBM is an appealing approach for solving flow passing complex geometry, whereas a very complicated body fitted mesh with fine resolution might be required in a conventional CFD approach. The results were compared with experimental data of flow velocity distribution through the grating measured by Particle Image Velocimetry technique. The flow immediately downstream of the grating is characterized by multimodal profile with intervals determined by a lattice width of the grating. This multimodal flow merges into one with unimodal shape in the downstream far from the grating. The present analysis reproduced the distinctive flow structure observed in the experiment. The model predicted a serrated profile consistent with the experimental results and this could be reproduced even with a relatively regular computational mesh compared to the traditional method.
Hamdani, A.; 安部 諭; 石垣 将宏; 柴本 泰照; 与能本 泰介
Energies (Internet), 13(14), p.3652_1 - 3652_22, 2020/07
被引用回数：6 パーセンタイル：18.77(Energy & Fuels)In the case of a severe accident, natural convection plays an important role in the atmosphere mixing of nuclear reactor containments. In this case, the natural convection might not in the steady-state condition. Hence, instead of steady-state simulation, the transient simulation should be performed to understand natural convection in the accident scenario within a nuclear reactor containment. The present study, therefore, was aimed at the transient 3-D numerical simulations of natural convection of air around a cylindrical containment with unsteady thermal boundary conditions (BCs) at the vessel wall. For this purpose, the experiment series was done in the CIGMA facility at Japan Atomic Energy Agency (JAEA). A 3-D model was simulated with OpenFOAM, applying the unsteady Reynolds-averaged Navier-Stokes equations (URANS) model. Different turbulence models were studied, such as the standard k-
, standard k-
, k- shear stress transport (SST), and low-Reynolds-k- Launder-Sharma. The results of the four turbulence models were compared versus the results of experimental data.
