Corresponding author: Valery P. Bereznev (bvp@ibrae.ac.ru)

Academic editor: Yury Kazansky

The importance of calculation of radiation fields inside in-reactor cavities is associated with the necessity to simulate the emergency modes in fast breeder reactors (FBR), as well as reactor states with different coolant levels in special dedicated channels of passive feedback devices in lead-cooled fast reactors (LFR) of BREST type or in sodium cavities in sodium-cooled fast reactors (SFR).

The Last Flight (LF) method (

In addition, it is desirable to implement the cavities calculation methodology within the framework of the approximations used in reactor calculations introducing certain specific features. In particular isotropy of the neutron flux density and the necessity of forced introduction of a “conditional” calculation cell on the boundary surface of the void cavity are assumed in the diffusion approximation. If the LF method is oriented on the connection of the source point with the detection point, then it is necessary to determine in the calculation of neutron field in the cavities the neutrons escaping the surface area of the source and neutrons reaching a certain surface area of the cavity. In order to solve the problem, the authors suggested using the approximate solution presented in the paper.

Thus, an algorithm for calculation of in-reactor cavities using the diffusion approximation was developed and implemented by the authors.

Neutron transport in steady-state conditions in the fast reactor core can be described, for instance, in multi-group diffusion approximation as follows:

where Ф^{g}^{g}^{g}_{r}^{l}^{→}^{g}^{g}^{g}_{f}

The main advantage of the diffusion approximation for practical implementation in computer codes is the high speed with which the solution is obtained and low requirements on the availability of computational resources. Nevertheless, certain applicability limits exist posed by the presence of strongly absorbing medium and pronounced spatial heterogeneity.

In particular, diffusion approximation is not applicable in the presence of gaseous (vacuous) medium. Methods for calculating radiation in such media can be based on the solution of integral neutron transport equation. If the purpose is to define radiation field outside the calculation area, then the last flight method (LF) is widely applied. However, in practice this method is not easily applicable for calculating internal cavities, since uncertainty of the method is high at the distances close to the surface. In such case volume sources used in LF-method can be replaced with surface sources with subsequent numerical integration over the surface separating dense and gaseous media.

Thus, the purpose of the present study is to create combined algorithm with associated calculation module converging the diffusion approximation producing fast solution with semi-analytical solution of integral equation inside in-reactor cavities.

Linear integrodifferential neutron transport equation can be written in the following integral form:

Equation (1) means that neutron flux φ(

We suggest spreading such representation of the solution only on “vacuous” calculation cells filled with gaseous medium where density of the matter is low to such extent that cross-sections of neutron interactions with medium are very small, so that diffusion factor exceeds by several orders of magnitude its values in a typical reactor medium, and not allowing using conventional diffusion approximation in the medium in question.

As it was indicated in (

where integration is performed over the surface separating dense and gaseous media; exp(– Σ^{g}^{2}) is the probability for neutron from ^{g}^{g}

Let us examine transmission of radiation using two-dimensional model of hexagonal calculation cell conventionally applied in the FBR calculations as represented in Figure

Model of transmission of radiation in hexagonal calculation cell with internal hexagonal structure (S is the center point of the length AB).

where integration is performed along the plane ^{g}

Average value of field on the “receiving” plane

In such case the value

is the contribution in the field on the plane

Assuming that field on the plane ^{g}^{g}_{AB}

The obtained integral can be calculated approximately and, in case of gaseous medium (Σ^{g}^{–4} cm^{–1}), Taylor expansion of the integration term can be used with (6) accepting the following form:

where γ = _{B} – x_{a}_{A} + x_{b}_{B}_{b}_{A}_{a}

In the simplest case for

where

Since the problem of radiation transmission in the hexagonal cell (see Fig.

Neutron flux density on the boundary separating the media can be determined from the integrodifferential neutron transport equation written for the boundary surface, i.e. for the zero-measure volume during integration of which all members associated with cell volume will disappear and only the flow members, integral for which will be converted from volume to surface ones, will remain. In such case we obtain equality of flows on the surface separating media

where φ^{g}_{i}_{is}_{si}_{is}_{j}_{j}^{g}_{i}_{i}_{s}^{g}_{s}

Since integration strongly depends on the cell geometry, then let us describe them.

There are cells with external steel cladding and central cells with coolant and some structures (end fittings of fuel rods in reactors of BN type (

Model of calculation cell with drainable internal hexagonal cell and arrangement of calculation nodes.

While the above mentioned homogenous medium fills the channel in the fuel assembly of BN-type reactor in normal operational conditions starting from the cells containing absorber to the top of reactor core and is substituted with gaseous medium in the case of boiling coolant, then in normal operational conditions of BREST-type reactors (

Corresponding model of calculation cells and nodes is represented in Figure

In case of complete coolant drainage from calculation cells of the central channel and their filling with inter gas, diffusion factor in such medium varies from unity to ~100 for positive feedback devices and to ~10^{4} for reactor facilities of BN type. Conventional diffusion approximation does not provide for the obtaining of acceptable result. Auxiliary calculation nodes become determining, the problem of obtaining average diffusion factor on the boundaries of calculation cells is removed and solution of integral neutron transport equation in gaseous medium is ensured.

Possible model of the positive feedback device includes external steel cladding and central steel technological pipe between which coolant is located filling the channel to the reactor core top in normal operational conditions and giving place to gaseous medium (argon) when pressure drops in the lower header of the first reactor cooling loop with forward gas stroke of 70 cm from the reactor core top and lower (Fig.

Model of calculation cell of the positive feedback device of BREST reactor facility and arrangement of calculation nodes.

Model of calculation cells and nodes in the positive feedback device channel is also represented in Figure

Such layout of arrangement of calculation nodes allows using conventional procedure in diffusion approximation with filling intermediate calculation cells with coolant when auxiliary calculation nodes are actually playing supporting role. However, when coolant is drained from auxiliary calculation cells calculation nodes on the surface separating media play determining role in the solution of neutron transport equations in gaseous medium.

Thus, concordance is achieved for the indicated layout between the solution of diffusion equation in dense media and solution of integral equation in gaseous medium.

Specific models of radiation in calculation cells filled with gaseous medium are examined below.

Let us examine possible situations in the propagation of radiation in hexagonal cell with internal hexagonal void.

When hexagonal cell (Fig.

Propagation of radiation from points А, S and B of the hexagonal cell (AS = AB/2).

Values of angles for points

When equation (9) is used for calculations, the initial section _{0}, _{1}, _{2} are the distances from the source to the receiving side, _{0} is the shortest distance; α_{1}, α_{2} are the aspect angles under which receiving sections

Dependence of flight distance on the aspect angle of the receiving section.

Let us examine possible situations with propagation of radiation in the hexagonal cell with gas void located between the external and internal dense media represented in hexagonal form.

First let us examine propagation of radiation from the internal dense hexagonal structure. For surface 7 with flat-to-flat dimension equal to _{c}

Propagation of radiation from the points of section ab on the side 7.

Let us analyze backward radiation from the external hexagonal structure to the internal one through gaseous medium. For this purpose, we will examine propagation of radiation for points of one of the side planes, for instance, the first one, depending on the dimensions of the internal hexagonal cell with flat-to-flat size _{c}_{c}_{c}_{c}_{c}_{c}_{c}_{c}

Propagation of radiation from points of the side plane _{c}

Propagation of radiation from points A, G, P, V and S.

Point

Results of comparison of design-basis drainage of 18 positive feedback devices in the version of design of BREST reactor facility obtained using Monte-Carlo method (

Deviation of efficiency of drainage of the positive feedback device to the height of 70 cm from the reactor core top from the designed value, %.

Diffusion approximation | Solution of integral equation, model in Fig. |
Solution of integral equation, model in Fig. |
---|---|---|

750.0 | –18.1 | –4.5 |

Calculation time for model of BREST-type reactor facility with model of positive feedback device presented in Fig. ^{–6} for neutron flux density of multi-group diffusion approximation.

Description is given of the algorithm for calculation in diffusion approximation of FBR in the presence of calculation cells containing voids. The algorithm for solution of neutron transport equation in gaseous medium is based on the semi-analytical methodology of solution of Peierls integral equation in the assumption that calculation cells with gaseous medium are surrounded by calculation cells in which diffusion approximation is used, which allows keeping high speed of obtaining the solution with high enough accuracy of this solution. No information about earlier application of such algorithm in international calculation practice is known to authors.

Preliminary results of application of the suggested algorithm are of high enough quality with high speed of obtaining the solution. As it was originally expected, direct application of diffusion approximation produces solution with accuracy which is not acceptable.

The developed algorithm for calculation in diffusion approximation of fast breeder reactor with presence of calculation cells with voids can be used for estimation of sodium void reactivity effect (

✩ Russian text published: Izvestiya vuzov. Yadernaya Energetika (ISSN 0204-3327), 2018, n. 2, pp. 67–77.