Development for CANDU-6 Moderator Temperature Prediction by Using Porous Media Approach

Author(s):  
Jae Ryong Lee ◽  
Han Young Yoon ◽  
Hyoung Tae Kim ◽  
Jae Jun Jeong

In this study, a thermal hydraulic behavior of the moderator in the CANDU reactor was numerically investigated by using CUPID code. KAERI has been developing a component-scale thermal hydraulics code, CUPID. The aim of the code is multi-dimensional, multi-physics and multi-scale thermal hydraulics analysis. This code adopts a three-dimensional, transient, two-phase and three-field model, and includes physical models and correlations of the interfacial mass, momentum, and energy transfer for the closure. To avoid the complexity to generate computational geometry around the matrix of 440 Calandria tubes, a porous media approach was applied. Flow resistance inside the porous media zone was derived from the empirical correlation of the frictional pressure loss. In order to consider the turbulent jet inflows from the inlet nozzles, the standard k-ε turbulence model was applied. For the grid dependency test, three different grid systems were tested. The moderator test vessel at Stern Laboratories Inc. (SLI) for the validation is a cylinder with a diameter of 2m and a length of 0.2m (a thin “slice” of CANDU-6 Calandria vessel). Since the axial flow is assumed to be invariant, two-dimensional calculation was performed. Vertical profile of the liquid temperature was compared with other calculation results as well as experimental data.

2017 ◽  
Vol 2017 ◽  
pp. 1-9 ◽  
Author(s):  
Zhengzhi Wang ◽  
Chunling Zhu

In view of the rotor icing problems, the influence of centrifugal force on rotor blade icing is investigated. A numerical simulation method of three-dimensional rotor blade icing is presented. Body-fitted grids around the rotor blade are generated using overlapping grid technology and rotor flow field characteristics are obtained by solving N-S equations. According to Eulerian two-phase flow, the droplet trajectories are calculated and droplet impingement characteristics are obtained. The mass and energy conservation equations of ice accretion model are established and a new calculation method of runback water mass based on shear stress and centrifugal force is proposed to simulate water flow and ice shape. The calculation results are compared with available experimental results in order to verify the correctness of the numerical simulation method. The influence of centrifugal force on rotor icing is calculated. The results show that the flow direction and distribution of liquid water on rotor surfaces change under the action of centrifugal force, which lead to the increasing of icing at the stagnation point and the decreasing of icing on both frozen limitations.


2021 ◽  
Vol 19 (2) ◽  
pp. 241
Author(s):  
Ruslan Balokhonov ◽  
Varvara Romanova ◽  
Eugen Schwab ◽  
Aleksandr Zemlianov ◽  
Eugene Evtushenko

A technique for computer simulation of three-dimensional structures of materials with reinforcing particles of complex irregular shapes observed in the experiments is proposed, which assumes scale invariance of the natural mechanical fragmentation. Two-phase structures of metal-matrix composites and coatings of different spatial scales are created, with the particles randomly distributed over the matrix and coating computational domains. Using the titanium carbide reinforcing particle embedded into the aluminum as an example, plastic strain localization and residual stress formation along the matrix-particle interface are numerically investigated during cooling followed by compression or tension of the composite. A detailed analysis is performed to evaluate the residual stress concentration in local regions of bulk tension formed under all-round and uniaxial compression of the composite due to the concave and convex interfacial asperities.


Author(s):  
Chaojian Chen ◽  
Mikhail Kruglyakov ◽  
Alexey Kuvshinov

Summary Most of the existing three-dimensional (3-D) electromagnetic (EM) modeling solvers based on the integral equation (IE) method exploit fast Fourier transform (FFT) to accelerate the matrix-vector multiplications. This in turn requires a laterally-uniform discretization of the modeling domain. However, there is often a need for multi-scale modeling and inversion, for instance, to properly account for the effects of non-uniform distant structures, and at the same time, to accurately model the effects from local anomalies. In such scenarios, the usage of laterally-uniform grids leads to excessive computational loads, both in terms of memory and time. To alleviate this problem, we developed an efficient 3-D EM modeling tool based on a multi-nested IE approach. Within this approach, the IE modeling is first performed at a large domain and on a (laterally-uniform) coarse grid, and then the results are refined in the region of interest by performing modeling at a smaller domain and on a (laterally-uniform) denser grid. At the latter stage, the modeling results obtained at the previous stage are exploited. The lateral uniformity of the grids at each stage allows us to keep using the FFT for the matrix-vector multiplications. An important novelty of the paper is a development of a “rim domain” concept which further improves the performance of the multi-nested IE approach. We verify the developed tool on both idealized and realistic 3-D conductivity models, and demonstrate its efficiency and accuracy.


Author(s):  
Yoshihiko Ishii ◽  
Kazuaki Kitou ◽  
Tomohiko Ikegawa ◽  
Shin Hasegawa ◽  
Hitoshi Ochi

Hitachi utilized three-dimensional transient analysis to design and verify a critical-control mode algorithm of an automatic power regulator (APR). TRACG has a three-dimensional neutron kinetics model based on diffusion theory and a six-equation two-phase flow model. To verify the APR critical-control mode algorithm, an external-neutron-source model that makes possible to simulate a sub-critical initial core, and an APR system model were developed and added on TRACG. The code was verified by comparison of measurements and calculation results of ABWR start-up operation under the critical-control mode. The modified TRACG could simulate neutron count rates of start-up-range neutron monitors (SRNMs), reactor period, control rod operation timing, CR withdrawal length, and time of criticality declaration, well.


2006 ◽  
Vol 28 (3) ◽  
pp. 134-144
Author(s):  
Nguyen The Duc

The paper presents a numerical method to simulate two-phase turbulent cavitating flows in ducts of varying cross-section usually faced in engineering. The method is based on solution of two-phase Reynolds-averaged Navier-Stokes equations of two-phase mixture. The numerical method uses artificial compressibility algorithm extended to unsteady flows with dual-time technique. The discreted method employs an implicit, characteristic-based upwind differencing scheme in the curvilinear grid systems. Numerical simulation of an unsteady three-dimensional two-phase cavitating flow in a duct of varying cross-section with available experiment was performed. The unsteady important characteristics of the unsteady flow can be observed in results of numerical simulation. Comparison of predicted results with experimental data for time-averaged velocity and phase fraction are provided.


Georesursy ◽  
2020 ◽  
Vol 22 (1) ◽  
pp. 4-12
Author(s):  
Timur R. Zakirov ◽  
Maxim G. Khramchenkov

This paper presents the results of numerical simulations of two-phase flows in porous media under capillary forces dominance. For modeling of immiscible two-phase flow, the lattice Boltzmann equations with multi relaxation time operator were applied, and the interface phenomena was described with the color-gradient method. The objective of study is to establish direct links between quantitative characteristics of the flow and invasion events, using high temporal resolution when detecting simulation results. This is one of the few works where Haines jumps (rapid invasion event which occurs at meniscus displacing from narrow pore throat to its wide body) are considered in three-dimensional natural pore space, but the focus is also on the displacement mechanics after jumps. It was revealed the sequence of pore scale events which can be considered as a period of drainage process: rapid invasion event during Haines jump; finish of jump and continuation of uniform invasion in current pore; switching of mobile interfaces and displacement in new region. The detected interface change, along with Haines jump, is another distinctive feature of the capillary forces action. The change of the mobile interfaces is manifested in step-like behavior of the front movement. It was obtained that statistical distributions of pressure drops during Haines jumps obey lognormal law. When investigating the flow rate and surface tension effect on the pressure drop statistics it was revealed that these parameters practically don’t affect on the statistical distribution and influence only on the magnitude of pressure drops and the number of individual Haines jumps.


Sign in / Sign up

Export Citation Format

Share Document