• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.8
• /
• pp.2069-2078
• /
• 1993

The buoyancy-driven turbulent thermal convection is predicted using an anisotropic hybrid turbulence model, which is incorporated with a low Reynolds k-.epsilon. turbulence model and an anisotropic buoyant part of algebraic stress model(ASM). The numerical predictions are compared with the Davidson's model,(1) the full ASM and the experimental results of Cheesewright et al.(2) All the models are shown to predict good agreements with the experiments for the averaged turbulence quantities. It is found that the effect of an anisotropic part on the Reynolds stress and the turbulent heat fluxes is substantial. In this study, the present hybrid model gives a fairly reasonable prediction in terms of the computational accuracy, convergence and stability. The contribution of an anisotropic buoyant part to turbulent heat fluxes are also scrutinized over the range of Rayleigh numbers $(4.79{\times}10^{10}{\le}Ra{\le}7.46{\times}10^{10}).$

• Nuclear Engineering and Technology
• /
• v.54 no.11
• /
• pp.4359-4372
• /
• 2022

Turbulent mixing in buoyant flows is an essential mechanism involved in many scenarios related to nuclear safety in nuclear power plants. Comprehensive understanding and accurate predictions of turbulent buoyant flows in the reactor are of crucial importance, due to the function of mitigating the potential detrimental consequences during postulated accidents. The present study uses URANS methodology to investigate the buoyancy-influenced flows in the reactor pressure vessel under the main steam line break accident scenarios. With a particular focus on the influence of turbulent heat flux closure models, various combinations of two turbulence models and three turbulent heat flux models are utilized for the numerical simulations of three ROCOM tests which have different characteristic features in terms of the flow rate and fluid density difference between loops. The simulation results are compared with experimental measurements of the so-called mixing scalar in the downcomer and at the core inlet. The study shows that the anisotropic turbulent heat flux models are able to improve the accuracy of the predictions under conditions of strong buoyancy whilst in the weak buoyancy case, a major role is played by the selected turbulence models with essentially a negligible influence of the turbulent heat flux closure models.

• Nuclear Engineering and Technology
• /
• v.36 no.6
• /
• pp.559-570
• /
• 2004

A computational fluid dynamics model for predicting moderator circulation inside the Canada deuterium uranium (CANDU) reactor vessel has been developed to estimate the local subcooling of the moderator in the vicinity of the calandria tubes. The buoyancy effect induced by the internal heating is accounted for by the Boussinesq approximation. The standard $k-{\varepsilon}$ turbulence model with logarithmic wall treatment is applied to predict the turbulent jet flows from the inlet nozzles. The matrix of the calandria tubes in the core region is simplified to a porous media in which the anisotropic hydraulic impedance is modeled using an empirical correlation of pressure loss. The governing equations are solved by DFX-4.4, a commercial CFD code developed by AEA technology. The resultant flow patterns of the constant-z slices containing the inlet nozzles and the outlet port are "mined-type", as observed in the former 2-dimensional experimental investigations. With 103% full power for conservatism, the maximum temperature of the moderator is $82.9^{\circ}C$ at the top of the core region. Considering the hydrostatic pressure change, the minimum subcooling is $24.8^{\circ}C$.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.162-171
• /
• 1993

This paper deals with the anisotropy of the mushy region during solidification process of a binary mixture. A theoretical model which specifies a permeability tensor in terms of pricipal values is proposed. Also, the governing equations are modified into convenient forms for the numerical analysis with the existing algorithm. Some test computations are performed for soeidification of aqueous ammonium chloride solution contained in a square cavity. Results show that not only the present model is capable of resolving fundamental characteristics of the tranport phenomena, but also the anisotropy significantly affects the interdendritic flow structure, i.e., double-diffusive convection and macrosegregation patterns.