All stages of investigations have given a great deal of attention to measurement methodologies and techniques, including development of unique velocity, flow rate, pressure, level, temperature and other sensors. Microthermocouples were developed to measure temperature in safety cans with an outer diameter of 0.3 to 0.8 mm operating in a temperature interval of 300 to 1800 °С. Flow meters of different designs were developed to measure flow rates of liquid metals. Methodologies and a technique were developed later for measuring electromagnetically the liquid metal local flow rate (velocity) vectors in channels and fuel rod bundles, and measuring the coolant mixing characteristics in experiments in air with a small fraction of gaseous tracers added in the form of Freon or propane (Rachkov et al. 2018, Sorokin et al. 2021b).

Much attention is given to methods of simulating physically experimental studies on hydrodynamics and heat exchange in liquid metal cooled nuclear power facilities. It has been experimentally shown that it is possible to simulate hydrodynamics of incompressible fluids, including liquid metals in experiments with air, and heat exchange in liquid metals, such as Na, Na-K, Li, Hg, Pb, Pb-Bi, etc., using simulating fluids (Sorokin and Kuzina 2019). The said methodologies have made it possible to undertake a broad range of fundamental and applied experiments.

Extensive studies have been undertaken into hydrodynamics of irregularly shaped channels, including rod bundles, and flow paths of reactor facilities; maximum attention was given to measurement of velocity fields, distribution of tangential stresses, and turbulent characteristics.

The results of experimental studies into hydrodynamic turbulent characteristics in fuel assemblies in air using a Pitot tube have shown (Fig. 1a, b) that there is a local maximum observed in the distribution of tangential stresses along the wetted perimeter of the regular fuel lattice cell at the point with the greatest channel expansion, which can be explained by the secondary vortex impact (Sorokin et al. 2021b).

Figure 1. 

Distributions of tangential stresses on the fuel rod surface and velocity in the channel cross-section in the peripheral region of a model FA with displacers in the peripheral channels with a nominal geometry (a, b) and during a fuel lattice strain in the FA peripheral region (c, d).

In a deformed lattice (Fig. 1c, d), the distribution of tangential stresses along the wetted surfaces is practically symmetrical with respect to the geometrical symmetry axes of the flow section, though anomalies are observed in some of the FA portions, which can be explained by the impact of certain secondary vortexes not only inside the channels, but also at the boundary. The velocity distribution along the normal to the wetted perimeter is described by a universal law if the local tangential stress value is used to calculate the dynamic velocity. There is a substantial intensification of turbulent velocity fluctuations observed in the bundle’s peripheral region as compared with an infinite lattice (Sorokin et al. 2021b).

Most attention has been given to thermohydraulic studies for the most heated and essential component of a reactor plant (the reactor core) affected, in the course of life, by a number of factors, including design, mode, process, radiation and operating factors (Rachkov et al. 2018). It was shown that there was no thermal (contact) resistance at the coolant – heat-exchange surface interface when the concentration of impurities in the coolant did not exceed their solubility at the circulating metal temperatures. In these conditions, heat transfer to liquid metals, such as Na, Na-K, Li, Hg, and Pb-Bi, in tubes is described by a single criterial relationship close to the Lion formula. When the coolant is saturated with impurities, the heat transfer coefficients are expected to decrease by a factor of one and a half to two, which corresponds to experiments by the Energy Institute (ENIN), the Central Boiler and Turbine Institute, and the IPPE (Rachkov et al. 2018).

Experimental and computational studies have shown the need for solving the ‘conjugate’ problem of heat removal from fuel rods with regard for their thermophysical properties. P.A. Ushakov developed a theory of approximate thermal similarity of fuel rods in regular lattices (Sorokin and Kuzina 2019) that has made it possible to simulate fuel rods by multilayer tubes electrically heated from the inside and is used in all investigations for fast reactor cores. Detailed experimental data have been obtained for thermal hydraulics of full-scale core models provided there are fuel rod deflections, asymmetrical displacements and deformations of components, closure of different core parts, and counter flows.

As a result of experimental studies and a computational and theoretical analysis of the mass, pulse and energy exchange among the channels in bundles of smooth and wire wrap finned fuel rods, physically justified methods and programs (ТЕМP, MIF) have been developed for thermohydraulic calculations of reshaped fast neutron reactor core fuel assemblies (Zhukov et al. 1991).

The influence of the fuel rod geometry and materials, and the radiation-induced swelling and creep effects on the FA temperature mode has been investigated, and peculiarities of the core temperature mode formation in the course of operation (life) have been identified for fast neutron reactors. The efficiency of using differently directed wire wraps leading to oppositely directed coolant flows in transverse directions has been shown.

A water flow pattern has been obtained for the flow path of a flat-plate DHS with central supply and side removal of water. It has been found that the liquid flow pattern in the header is defined by the DHS dimensional ratio and design.

The liquid flow in the cylindrical-type DHS flow path is of a complex nature and is defined predominantly by the DHS dimensional ratio and design, the liquid flow pattern, and the hydraulic resistance coefficient for the outlet component flow path. Representative models of liquid flows in the flow path of the above DHS type (Fig. 2) have been obtained with regard for results of various experimental studies (Gabrianovich and Delnov 2016).

Figure 2. 

Typical designs and liquid flow models in axisymmetrical DHSs of a cylindrical type with central supply of water to the header: a, b. Superconstricted DHSs without the central tube extension with and with no distributor respectively; c. constricted DHSs with the central tube extension and a header with a constricted inlet section; d. Constricted DHS with a distributor, the central tube extension and a header with a free inlet section; e, f. Constricted DHSs with a distributor, the central tube extension, and a header with constricted and free inlet sections; 1 – central tube; 2 – tube sheet; 3 – bundle tube; 4 – housing; 5 – stage; 6 – bottom; 7 – header; 8 – side annulus, 9 – distributor.

A water flow pattern was obtained on a water table for the flow path of a flat-plate DHS with side supply and central removal of water. It has been found that the liquid flow pattern in the header is defined by the DHS dimensional ratio and design.

The liquid flow in the flow path of a cylindrical-type DHS is of a complex type and is defined predominantly by the DHS dimensional ratio and design, the liquid flow pattern, and the hydraulic resistance of the lattice. Representative models of the liquid flow in the flow path of the considered DHS type have been obtained with regard for the results of experimental studies (Gabrianovich and Delnov 2016). The most graphic peculiarities of the liquid flow manifest themselves in the inlet, main and outlet portions of the flow path in the DHS under consideration (Fig. 3).

Figure 3. 

Typical designs and liquid flow models in the flow paths of axisymmetrical DHSs of a cylindrical type with side supply of liquid to the header: a, b. constricted DHSs with the displaced tube sheet in a shell and with free and constricted inlet sections respectively; c, d. constricted DHSs with a constricted inlet section, the displaced tube sheet in a shell and inserts of relatively large and small diameters respectively; e, f. superconstricted DHSs with a header constricted by the inlet section, with no lattice displacement in a shell with no and with inserts respectively; 1 – annulus; 2 – header; 3 – shell; 4 – lattice; 5 – bottom; 6 – insert; 7 – housing.

As a result of the studies, the earlier unknown regularity and phenomenon, dealing with nuclear, space, metallurgical and chemical fields of science and technology, have been identified and registered as scientific discoveries.

It has been found that there is an earlier unknown regularity of the liquid distribution at the outlets of the flow paths in distributing header systems that consists in that axisymmetrical regions form as the liquid exits the header, the characteristics of which are defined by the design and process peculiarities of the header system (liquid supply point, motion path, jet parameters, hydraulic resistance, etc.) (Delnov et al. 2019).

It has been found that there is an earlier unknown phenomenon of the hydrodynamic identity occurrence in distributing header systems that consists in the similarity of the hydrodynamic characteristics of the flow paths in axisymmetrical distributing header systems, e.g., nuclear power facilities and heat exchangers with different conditions of supplying and removing the liquid flowing in the system (Delnov 2021).

DHS design differences: a tube sheet and a system of plates are used as the outlet component respectively in cylindrical and flat-plate DHSs:

• axisymmetrical round, cylindrical, conical and annular jets occur in DHSs of a cylindrical type, and jets with a rectangular (square) cross-section are typical of flat-plate DHSs;
• in cylindrical DHSs, one type of jets is transformed upstream into another type, and flat-plate DHSs have the common jet divided into individual parts or individual parts of the jet converge into one jet;
• DHSs with different liquid supply and removal points differ in the sequence for transformation of certain types of jets into others.

Scientific discoveries

• have changed the existing scientific concepts in the field of hydrodynamics of axisymmetrical DHSs;
• have explained the scientific facts and experimental data not rationalized in scientific terms earlier;
• have shown an extremely strong effect of minor changes in the DHS design and dimensional ratio on the DHS liquid flow;
• have made it possible to obtain empirical relationships for determining the hydrodynamic irregularities at the outlet of different DHS designs and types;
• have allowed predicting, prior to experiments and calculations, the DHS designs with the required liquid flow profile at the header outlet;
• contain data on the models of liquid flows in different DHS types;
• have identified the existence of constricted and free DHSs the liquid flow in which differs greatly from that in thoroughly studied classic superconstricted and free DHSs.

Scientific discoveries have been used to justify the flow paths of DHSs in reactors and heat exchangers of nuclear power facilities and to develop and verify the DHS flow fluid dynamics codes.

Experimental studies based on an integrated water model of reactors using a developed theory of approximation simulation, temperature fields and the structure of the nonisothermic coolant movement in the primary circuit components of a fast neutron reactor for forced circulation modes, and changeover to the cooldown mode and emergency cooldown by natural coolant circulation (Opanasenko et al. 2017), demonstrate a substantial and stable temperature stratification of the coolant in the peripheral region of the reactor’s upper (hot) chamber above the side shields, in the cold and discharge chambers, in the elevator baffle, in the reactor vessel cooling system, and at the intermediate and auxiliary heat exchanger outlets in different modes of the heat exchanger operation (Figs 4, 5).

Figure 4. 

Axial upper chamber distributions of the averaged coolant temperature (a) and the intensity of temperature fluctuations (b) obtained during the movement of mobile temperature probes along the upper chamber height in the nominal mode of the plant operation.

Figure 5. 

Averaged temperature field along the upper chamber height in the steady-state natural convection cooldown mode.

High gradients and fluctuations of temperature have been recorded at the interfaces of stratified and recirculating formations. The obtained results can be used for verification of codes and coarse estimation of the reactor plant parameters during recalculation based on similarity criteria.

The Protva and Ugra codes, set up based on a well-developed theory of an anisotropic porous body for the calculation of complex flows in reactors, heat exchangers and steam generators, were used to prove the possibility of using in the BN-800 design the heat exchangers with the same heat-transfer surface as in the BN-600 design. Characteristics of heat exchange, critical heat fluxes and circulation stability have been studied for steam generators of reactor plants with the BN-350, BN-600 and BN-800 reactors and a fundamentally new large-module steam generator for an advanced fast reactor (Grabezhnaya and Mikheev 2015).

Investigations into the liquid metal boiling based on fuel assembly (FA) models have shown three patterns for a two-phase liquid metal flow in fuel rod bundles (bubble, slug and annular dispersed flow patterns) which is limiting with respect to the assembly cooling. It has been shown that long-term cooling of the core is conceptually possible in emergency modes involving boiling of liquid metals. Heat transfer during boiling of liquid metals in fuel rod bundles has been studied, the effects of the fuel rod surface roughness on the development of the boiling process has been investigated, and a diagram of the two-phase liquid metal flow patterns in fuel rod bundles has been plotted [Sorokin et al. 2018, 2019].

The results of computational studies for a system of parallel FAs (Fig. 6), undertaken based on an upgraded version of the SABENA code (Sorokin et al. 2021a) implementing a two-liquid model of a two-phase liquid metal flow in an approximation of equal pressures in a steam or liquid phase, reproduce the variation of temperature, the development of single-phase (bubble, slug) flow patterns, and the liquid metal flow rate fluctuations obtained in experimental studies, and demonstrate antiphase coolant flow rate fluctuations in parallel FAs, the interchannel stability characterized by a major growth in the coolant flow rate fluctuation amplitude in parallel FAs as compared with single FAs, a periodic drop in the FA coolant flow rate practically to zero, and potential FA exposure (departure from nucleate boiling).

Figure 6. 

Comparison of calculated and experimental time-dependent distributions of the fuel simulator surface temperature (a), coolant temperature (b) and coolant flow rate (c) in the left-hand FA during parallel operation of FAs.

To exclude the development of an emergency caused by an operational occurrence with simultaneous loss of power for the reactor coolant pumps and a failure of the reactor scram rods (an ULOF accident), a design solution has been proposed with a ‘sodium plenum’ arranged above the reactor core. Comparing the calculation and experiment results has shown the possibility of heat removal by the boiling coolant in a model FA with a ‘sodium plenum’ during thermal loads of 10% to 15% and the sodium flow rate level of about 5% of the nominal values (Sorokin et al. 2018).

The activities to systematize thermophysical data in the field of alkali cooled fast reactors were undertaken under the auspices of the IPPE’s Thermophysical Data Center in many fields of fast reactor thermal physics: velocity and temperature fields in the core, in the hot chamber, and in heat exchangers and steam generators, hydraulic resistance and heat transfer in channels and fuel rod bundles for a single-phase flow, heat transfer and the diagram of two-phase flow patterns, departure from nucleate boiling in assemblies, and thermophysical properties of coolants and reactor materials. Technical guidelines and reference books have been developed (Kirillov et al. 2010).

• refining methods to calculate local turbulent characteristics of pulse and energy transport for single- and two-phase liquid metal flows in channels and in large volumes taking into account large-scale eddy currents, and the coolant flow stratification effects;
• development and justification of a system of verification tests;
• development of a verified system of codes taking into account the interconnection of nuclear physical, thermohydraulic, physicochemical, thermomechanical, mass-exchange and engineering processes taking place in a nuclear power facility, for justifying its service life with regard for the entire combination of its operating processes and modes;
• evaluation of temperature fluctuations directly in the coolant flow and on the channel walls, and investigation of the effects these fluctuations have on structural strength;
• justification of the core’s temperature mode taking into account random deviations of parameters (geometrical, operating and process parameters, uncertainties of thermophysical properties, calculated constants, etc.);
• analyzing the consequences of potential nonstandard modes (interlocks, emergency heat removal, boiling) and development of measures to prevent these from progressing into a severe accident;
• investigating the dynamics of the sodium boiling region propagation in a real FA with a ‘sodium plenum’ above the core;
• justification of nominal modes that exclude the formation of eddies in the core header and on the sodium surface (gas capture), and passive circulation zones (stratification phenomena, temperature fluctuations);
• estimation of temperature fluctuations directly in the coolant flow and on the channel walls, and investigation of the effects these fluctuations have on structural strength;
• development of a system for detection of anomalies in an individual FA or in a number of FAs prior to the fuel cladding failure;
• justifying, for large-unit SGs, thermohydraulic modes and hydrodynamic stability, and integrated tests of the automatic steam generator protection systems (all scram systems shall be supplied and assembled to support the tests, and the scram system bench shall be put into operation);
• development of the material and design for improving the safety of the large-unit steam generator to provide for the slowdown of the leakage self-development and development processes, and online repair of SGs after the water escape into sodium.