Fukushima plutonium effect and blow-up regimes in neutron-multiplying media

It is shown that the capture and fission cross-sections of 238U and 239Pu increase with temperature within 1000-3000 K range, in contrast to those of 235U, that under certain conditions may lead to the so-called blow-up modes, stimulating the anomalous neutron flux and nuclear fuel temperature growth. Some features of the blow-up regimes in neutron-multiplying media are discussed.


Introduction
It is known that after the loss of coolant at three nuclear reactors during Fukushima nuclear accident its nuclear fuel melted. It means that the temperature in the active zone at some moments reached the melting point of uranium-oxide fuel 1 , i.e. ∼3000 • C.
Surprisingly enough, scientific literature today contains absolutely no either experimental or even theoretically calculated data on behavior of the 238 U and 239 P u capture and fission cross-sections depending on temperature at least in 1000-3000 • C range. At the same time there are serious reasons to believe that the cross-section values of the mentioned elements increase with temperature. At least we may point, for example, to the qualitative estimates by Ukraintsev [1], Obninsk Institute of Atomic Energetics (Russia), that confirm the possibility for 239 P u cross-sections growth in 300-1500 • C range.
Obviously, such anomalous temperature dependence of 238 U and 239 P u capture and fission cross-sections may change the neutron and heat kinetics of nuclear reactors drastically. This is also true for the perspective new generation fast reactors (uranium-plutonium of Feoktistov [2] and thorium-uranium of Teller [3] type), that we classify as fast TWR reactors. Hence it is very important to know the anomalous temperature behavior of 238 U and 239 P u capture and fission cross-sections, and furthermore it becomes critically important to know their influence on the heat transfer kinetics, since it may become a reason of the positive feedback 2 (PF) with neutron kinetics leading to undesirable solution stability loss (the nuclear burning wave) as well as to a trivial reactor runaway with a subsequent nontrivial catastrophe.
A special case of PF is a non-linear PF, which leads to the system evolution in socalled blow-up mode [4][5][6][7][8][9], or in other words, in such a dynamic mode when one or several modeled values (e.g. temperature and neutron flux) grows to infinity at a finite time. In reality, a phase transition is observed instead of the infinite values in this case, and this can in its turn become a first stage or a precursor of the future technogenic disaster.
Investigation of the temperature dependence of 238 U and 239 P u capture and fission cross-sections in 300-3000 • C range, and correspondingly, the heat transfer kinetics and its influence on neutron kinetics in TWR, is the main goal of the present paper.
2 Temperature blow-up regimes in neutron-multiplying media Heat transfer equation for uranium-plutonium fissile medium is: ρ ( r, T, t) ·c ( r, T, t) · ∂T ( r, t) ∂t = = ℵ ( r, T, t) · ∆T ( r, t) + ∇ℵ ( r, T, t) · ∇T (r, t) + q f T ( r, T, t) , where the effective medium density is ρ i are tabulated values, N i ( r, T, t) are the concentrations of the medium components, while the effective specific heat capacity (accounting for the medium components heat capacity values c i ) and fissile material heat conductivity coefficient (accounting for the medium components heat conductivity coefficients ℵ i (T )) respectively are: Here q f T ( r, T, t) is the heat source density produced by nuclear fissions N i of fissile metal components which vary in time.
Theoretical temperature dependence of heat capacity c(T ) for metals is known: at low temperatures c(t) ∼ T 3 , and at high temperatures c(T ) → const, where the constant value (const ≈ 6 Cal/(mol · deg)) is determined by Dulong-Petit law. At the same time it is known that thermal expansion coefficient is small for metals, therefore the specific heat capacity at constant volume c v almost equals to the specific heat capacity at constant pressure c p . On the other hand, the theoretical dependence of heat conductivity ℵ i (T ) at high temperature of "fissile" metals is not known, while it is experimentally determined that the heat conductivity coefficient ℵ(T ) of fissile medium is a non-linear function of temperature (e.g. see [10], where heat conductivity coefficient is given for α-uranium 238 and for metallic plutonium 239, and also [11]).
Further for solving thermal conductivity equation we used the following initial and boundary conditions: where j n is the normal (to the fissile medium boundary) heat flux density component, ℵ(T ) is the thermal conductivity coefficient, is the fissile medium boundary, T 0 is the temperature of the medium adjacent to the active zone.
Obviously, if the cross-sections of some fissile nuclides increase with temperature, then due to exothermic nature of the nuclei fission reaction, the significantly non-linear kinetics of mother and daughter nuclides in the nuclear reactor will immediately result in an autocatalytic increase of generated heat, just like in the autocatalytic processes of exothermic chemical reactions. The heat source density q f T ( r, Φ, T, t) which characterizes the amount of generated heat in this case will be: where is the total neutron flux density; Φ ( r, E, T, t) is the density of a neutron flux with energy E; Q f i is the mean heat released in one fission event of the i-th nuclide; is the fission cross-section of the i-th nuclide averaged over the neutron spectrum; is the probability density function of the neutron energy distribution; σ i f (E, T ) is the microscopic fission cross-section of the i-th nuclide, which is known to depend on the neutron energy and fissile medium temperature (Doppler effect [12]); N i ( r, T, t) is the density of the i-th nuclide nuclei.
As follows from (6), in order to build a density function of the heat source, it is necessary to solve a problem related to the construction of a theoretical dependence of the cross-sections σ i f ( r, T, t) averaged over the neutron spectrum on the temperature of reactor fuel (fissile medium). As is known, the impact of the nuclei thermal motion in the medium comes to a broadening and lowering of the resonances. By optical analogy this phenomenon is usually referred to as Doppler effect [12]. Since the resonance levels in the low energy region are observed for heavy nuclei only, the Doppler effect is noticeable only during the interaction of neutrons with such nuclei. And the higher is the environment temperature, the stronger is the effect.
Therefore a program was developed using Microsoft Fortran Power Station 4.0 (MFPS 4.0) that allows at the first stage to calculate the cross-sections of the resonance neutron reactions depending on neutron energy taking into account the Doppler effect. The crosssections dependence on neutron energy for reactor nuclides from ENDF/B-VII database [13] corresponding to 300K environment temperature were taken as the input data for the calculations. For example, the results for radioactive neutron capture cross-sections dependence on neutron energy for 235 U are given in Fig. 1 for different temperatures of the fissile medium in 300K-3000K temperature range. Using this program, the dependence of scattering, fission and radioactive neutron capture cross-sections for essential reactor fuel nuclides 235 92 U, 238 92 U, 239 92 U and 239 94 Pu were obtained for different temperatures in 300K-3000K range.
At the second stage a program was developed to calculate the dependence of the cross-sections σ i f ( r, T, t) averaged over the neutron spectrum for main reactor nuclides and for main reactor nuclear neutron reactions for the specified temperatures. The averaging of neutron cross-sections for the Maxwell distribution was performed using the following expression: where E lim is the the upper limit of neutrons thermalization, while for the procedure of neutron cross-sections averaging over Fermi spectrum the following expression was used: During further calculations in our programs we used the results obtained at the first stage i.e. the dependence of reaction cross-sections on neutron energy and medium temperature (Doppler effect). The neutron spectrum was specified in a combined way -below the limit of thermalization E lim the neutron spectrum was described by Maxwell spectrum Φ M (E n ); above E lim but below E F (upper limit for Fermi neutron energy spectrum) the neutron spectrum was described by Fermi spectrum Φ F (E) for a moderating medium with absorption; above E F , but below maximal neutron energy E max n the spectrum was described by 239 P u fission spectrum [14,15]. Here the neutron gas temperature for Maxwell distribution was given by (7), described in [12]. According to this approach [12], the drawbacks of standard slowing-down theory for thermalization area may be formally reduced if a variable ξ(x) = ξ(1 − 2/z) is introduced instead of the average logarithmic energy loss ξ, which is almost independent of neutron energy (as is known, for environment consisting of nuclei with A > 10 the statement ξ ≈ 2/A is true). Here z = E n /kT , E n is the neutron energy, T is the environment temperature. Then within such a framework the following expression may be used for the temperature of the neutron gas in Maxwell spectrum of thermal neutrons 1 : where T 0 is the fuel medium temperature, Σ a (kT 0 ) is the absorption cross-section for energy kT 0 , η = 1.8 is the dimensionless constant, ξ is averaged over the whole energy interval of Maxwell spectrum ξ(z) at kT = 1 eV . Fermi neutron spectrum for a moderating medium with absorption (we considered carbon as a moderator and 238 U, 239 U and 239 Pu as absorbers) was set in the form [12,16]: where S is the total volume neutron generation rate, ξ = i ξ i Σ i S /Σ S , ξ i is the average logarithmic decrement of energy loss, Σ i S is the macroscopic scattering cross-section of the i-th nuclide, Σ t = i Σ i S + Σ i a is the total macroscopic cross-section of the fissile material, S is the total macroscopic scattering cross-section of the fissile material, Σ a is the macroscopic absorption cross-section, E F is the upper neutron energy for Fermi spectrum.
The upper limit of neutron thermalization E lim in our calculation was considered as a free parameter, setting the neutron fluxes of Maxwell and Fermi spectra at a common energy limit E lim equal: The high energy neutron spectrum part (E > E F ) was defined by fission spectrum [16][17][18] in our calculations. Therefore for the total volume neutron generation rate S in the expression for the Fermi spectrum (8) the following expression may be written: where E max n is the maximum energy of neutron fission spectrum (usually taken as E max n ≈ 10 M eV ), E F is the neutron energy, below which the moderating neutrons spectrum is described as Fermi spectrum (usually taken as E F ≈ 0.2 M eV );P ( r, E, T, t) is the probability for the neutron not to leave the boundaries of the fissile medium, which depends on the fissile material geometry and conditions at its border as well (e.g. presence of a reflector).
The obtained calculation results show that the cross-sections averaged over the spectrum may increase (Fig. 2 for 239 P u and Fig. 4 for 238 U ) as well as decrease (Fig. 3 for 235 U ) with fissile medium temperature increase. As follows from the obtained results, the arbitrariness in selection of the limit energy for joining Maxwell and Fermi spectra does not alter the character of these dependences evolution significantly. This can be justified by the fact that 239 Pu resonance region starts from significantly lower energies than that of 235 U and with fuel temperature increase, the neutron gas temperature increases producing Maxwell's neutron distribution maximum shift to higher neutron energies. In other words, the neutron gas spectrum hardening, when more neutrons fit into resonance area of 239 Pu, is the cause of the averaged cross-sections growth.
For 235 U this process in not as significant because its resonance region is located at higher energies. As a result, the 235 U neutron gas spectrum hardening related to fuel temperature increase (in the considered interval) does not result in a significant increase a) b) Figure 2: Temperature dependences for the fission cross-section (a) and radioactive capture cross-section (b) for 239 P u, averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra joining energy and η = 1.8 (see (7)). a) b) Figure 3: Temperature dependences for the fission cross-section (a) and radioactive capture cross-section (b) for 235 U , averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra joining energy and η = 1.8 (see (7)).
a) b) Figure 4: Temperature dependences for the fission cross-section (a) and radioactive capture cross-section (b) for 238 92 U, averaged over the combined Maxwell and Fermi spectra depending on Maxwell and Fermi spectra joining energy and η = 1.8 (see (7)). of the number of neutrons fitting into the resonance region. Therefore according to the known dispersion relations for 235 U giving the neutron reactions cross-sections behaviour depending on their energy E n for non-resonance areas, we observe a dependence for the averaged cross-sections σ nb ∼ 1/ √ E n . The data on the averaged fission and capture cross-sections of 238 U presented at Fig. 4 show that the averaged fission cross-section for 238 U is almost insensitive to the neutron spectrum hardening caused by the fuel temperature increase, because of the high fission threshold ∼1 MeV (see fig. 4a). On the other hand, they confirm the dependence of the capture cross-section on temperature, because its resonance region is located as low as for 239 Pu. Obviously, fuel enrichment with 235 U essentially makes no difference in this case, because the averaged cross-sections for 235 U, as described above, behave in a standard way.
And finally we performed a computer estimate of the heat source density dependence q f T ( r, Φ, T, t) (6) on temperature for different compositions of uranium-plutonium fissile medium with a constant neutron flux density, presented at fig. 5. We used the dependences presented above at Fig.2-4 for these calculations. Let us note that our preliminary calculations were made not taking into account the change in composition and density of the fissile uranium-plutonium medium that is a direct consequence of the constant neutron flux assumption.
The necessity of such assumption is caused by the following. It is obvious that for a reasonable description of the neutron source density q f T ( r, Φ, T, t) (6) dependence on temperature, a system of three equations must be solved. Two of them correspond to the neutron kinetics equation (flux and fluence) and to the system of equations for the parental and child nuclides nuclear density kinetics (e.g. see [19,20]), while the third equation corresponds to a heat transfer equation of (1) type. However, some serious difficulties arise at this point, associated with the limited computational capabilities available. And here is why.
One of the principal physical peculiarities of TWR is the fact [21] that fluctuation residuals of plutonium (or 233 U in Th-U cycle) over its critical concentration burn out for the time comparable with the reactor lifetime of a neutron τ n (x, t) (not considering delayed neutrons), or at least comparable with the reactor period 1 T (x, t) (considering delayed neutrons). Meanwhile, the new plutonium (or 233 U in Th-U cycle) is formed in a few days (or a month) and not at once. It means [21] that numerical calculation must be performed with a temporal step around 10 −6 -10 −7 for the case of not taking into account the delayed neutrons and ∼10 −1 -10 0 otherwise. At first glance, taking into account the delayed neutrons, according to [21], really "saves the day", however it is not always true. If the heat transfer equation contains a significantly non-linear source, then in the case of a blow-up mode, the temperature at some conditions may grow extremely fast and in 10-20 steps (with time step 10 −6 -10 −7 s) reaches the critical amplitude that may lead to (as a minimum) a solution stability loss or (as a maximum) blow-up bifurcation of the phase state, almost unnoticeable with a rough time step.
Here we should also mention that modern scientific literature has absolutely no theoretical or experimental data on heat capacity c p and heat conductivity ℵ of the fissile material for the temperatures over 1500K, which makes the further model calculations quite problematic.
According to these remarks, and considering the goal and format of this paper, we didn't aim at finding the exact solution of some specific system of three joined equations described above. However, we found it important to illustrate -at the qualitative levelthe consequences of the possible blow-up modes in case of non-linear heat source existence in the heat transfer equation. As it was described above, we performed the estimate computer calculations of the heat source density q f T ( r, Φ, T, t) (6) dependence on temperature in 300-1400K range for selected compositions of uranium-plutonium fissile medium at a constant neutron flux (Fig. 5).
The obtained dependences for the heat source density q f T ( r, Φ, T, t) were successfully approximated by a power function of temperature with an exponent of 4 (Fig. 5). In other words, we obtained a heat transfer equation with a significantly nonlinear heat source in the following form: where δ > 1 in the case of a non-linear heat conductivity dependence on temperature [4][5][6][7][8]. The latter means that the solutions of the heat transfer equation (1) describe the so-called Kurdyumov blow-up modes [4][5][6][7][8][9], i.e. such dynamic modes when one of the modeled values (e.g. temperature) turns into infinity for a finite time interval. As noted before, in reality instead of reaching infinite values, a corresponding phase transition is observed (a final phase of the parabolic temperature growth), that requires a separate model and is a basis for an entirely new problem. Mathematical modeling of blow-up modes was performed mainly using Mathematica 5.2-6.0, Maple 10, Matlab 7.0, utilizing multiprocessor calculations for effective application. Runge-Kutta method of 8-9 th order and the numerical methods of lines [22] were also applied to the calculations. The numerical error estimate did not exceed 0.01%. The coordinate and temporal step were variable and chosen by the program in order to fit the given error at every calculation step.
Below we give the solutions for the heat transfer equation (1) with nonlinear exponential heat source (11) in uranium-plutonium fissile medium for the boundary and initial parameters corresponding to those of the technical reactors. The calculations were done for a cube of fissile material with different spatial size, boundary and initial temperature values. Since the temperature dependences of the heat source densities were obtained without account for changing composition and density of the uranium-plutonium fissile medium, different blow-up modes can take place (HS-mode, S-mode, LS-mode) depending on the ratio between the exponents of the temperature dependences of thermal conductivity and heat source according to [4][5][6][7][8][9]. Therefore we considered cases for 1 st , 2 nd and 4 th temperature order sources. Here the power of the source also varied by varying the proportionality factor in (11) (const = 1.00J/(cm 3 ·s·K) for 1 st temperature order source; 0.10J/(cm 3 · s · K 2 ), 0.15 J/(cm 3 · s · K 2 ) and 1.00 J/(cm 3 · s · K 2 ) for 2 nd temperature order source; 1.00 J/(cm 3 · s · K 4 ) for 4 th temperature order source). While calculating the the heat capacity c p (Fig. 6a) and thermal conductivity ℵ (Fig. 6b) dependences on the fissile medium temperature in 300-1400K range, the specified parameters were given by the analytic expressions, obtained by approximation of the experimental data for 238 U based on the polynomial progression: And finally a solution of the heat transfer equation (1) was obtained for the constant thermal conductivity (27.5 W/(m · K)) and heat capacity (11.5 J/(K · mol)) values, presented in Fig. 7a, as well as the solutions of the heat transfer equation considering their temperature dependences (Fig. 7b-d).
Preliminary results point directly to a possibility of the local melting of the uraniumplutonium fissile medium, having melting temperature almost identical to that of 238 U, that equals 1400K (Fig. 6a-d). Moreover, these regions of the local melting are not the areas of the so-called thermal spikes [23], and probably are the anomalous areas of the uranium surface melting observed by Laptev and Ershler [24] that were also mentioned in [25]. More detailed analysis of the probable temperature scenario associated with the blow-up modes is discussed below.

Physical peculiarities of the blow-up regimes in neutron-multiplying media
Earlier we noted the fact that due to the coolant loss at nuclear reactors during Fukushima nuclear accident the fuel was melted, or in other words, temperature inside the active zone at some moment reached the melting temperature of the uranium-oxide fuel, i.e. ∼3000 • C.
On the other hand, we already know that the coolant loss may become a cause of the nonlinear heat source formation inside the nuclear fuel and thus become a cause of the   (1) solution for 3D case (crystal sizes 0.001×0.001×0.001 mm; initial and boundary temperatures equal to 100K): a) The source is proportional to the 4 th order of temperature; const =1.00 J/(cm 3 · s · K 4 ), heat capacity and thermal conductivity are constant and equal to 11.5 J/(K · mol) and 27.5 W/(m · K) respectively; b) the source is proportional to the 4 th order of temperature; const =1.00 J/(cm 3 · s · K 4 ); c) The source is proportional to the 2 nd order of temperature; const =1.00 J/(cm 3 · s · K 2 ); d) the source is proportional to the 2 nd order of temperature; const = 0.10 J/(cm 3 · s · K 2 ). Note: in cases b)d) the heat capacity and thermal conductivity were determined by (12) and (13) respectively. temperature and neutron flux blow-up mode occurrence. A natural question arises as to whether it is possible to use such blow-up mode (in terms of temperature and neutron flux) for the initiation of certain controlled physical conditions, under which the nuclear burning wave would regularly "experience" the so-called controlled blow-up regime. It is quite difficult to answer this question definitely, because such fast process has some physical vaguenesses, any of which can become experimentally unsurmountable for the its controlling.
Nevertheless such process is very elegant and beautiful from the physical point of view and therefore requires more detailed phenomenological description. Let us try to make it in short.
As one can see from the plots of the capture and fission cross-sections evolution for 239 Pu (Fig. 2), the blow-up mode may develop actively at ∼1000-2000K (depending on the real value of Fermi and Maxwell spectra joining boundary), but it returns to almost the initial cross-sections values at temperatures over 2500-3000K. If we turn on some effective heat sink at that point, the fuel may return to its initial temperature state. However, while the blow-up mode develops, the fast neutrons already penetrate to the adjacent fuel areas, where the new fissile material starts accumulating and so on (see cycles (1) and (2)). After some time the similar blow-up mode will develop in this adjacent area and everything starts over again. In other words, such hysteresis blow-up mode, closely time-conjugated to a heat takeoff procedure, will appear against the background of a stationary nuclear burning wave in a form of the periodic impulse bursts.
In order to demonstrate the marvelous power of such process, we investigated the heat transfer equation with non-linear exponential heat source in uranium-plutonium fissile medium with the boundary and initial parameters emulating the heat takeoff process. In other words, we investigated the blow-up modes in the fast Feoktistov-type uraniumplutonium reactor [2] where the temperature was deliberately fixed at 6000K inside and outside the boundary. This temperature is defined by the following important question: "Is it possible to obtain a solution, i.e. the spacio-temporal temperature distribution not in a form of a δ-function at some local spatial area, but as some kind of a stationary and limited by amplitude solitary wave under such conditions (6000K), which emulate the time-conjugated heat takeoff (see Fig. 2)?" As shown below, such suggestion proved productive.
Below we present some calculation characteristics and parameters. During these calculations we used the following expression for the heat conductivity coefficient: which was obtained using Wiedemann-Franz law and the data on electric conductivity of metals at temperature 6000K [26]. Specific heat capacity at constant pressure was set to c p ≈ 6 cal/(mol · deg) according to Dulong and Petit law.
The fissile uranium-plutonium medium was modeled as a cube of the size 10.0×10.0×10.0 m (Fig.8) during the calculations. Here we used the 2nd order temperature dependence for the heat source (see (11)).
And finally Fig. 8a-d present a set of solutions of the heat transfer equation (1) with nonlinear exponential heat source (11) in uranium-plutonium fissile medium with boundary and initial conditions emulating such process of heat takeoff that initial and boundary temperatures remain constant and equal to 6000K.
It is important to note here, that the solution set presented at Fig. 8 demonstrates the temporal evolution of the solution to its "stationary" state quite clearly. This is achieved using the so-called "magnifying glass" approach when the solutions of the same problem are deliberately investigated at the different timescales. For example, Fig. 8a shows the solution at the time scale t ∈ [0, 10 −6 s], while Fig. 8b describes the spatial solution of the problem (temperature field) for t = 10 −6 s. The Fig. 8c-d present the solution (the spatial temperature distribution) at t = 0.5 s and t = 50 s. As one can see, the solution (Fig. 8d) is completely identical to the previous (Fig.8c), i.e. to the distribution established in the medium in 0.5 seconds, which allows us to make a conclusion about the temperature field stability, starting from some moment. It is interesting that the established temperature field creates the conditions enough for thermonuclear synthesis reaction, i.e. reaching 10 8 K, and the time of such temperature field existence is not less than 50 s. These conditions are highly favorable for a stable thermonuclear burning given the necessary nuclei concentration to enter the thermonuclear synthesis reaction.
That said, one should remember that the results of the current chapter are only demonstrative, because their accuracy is very relative and requires careful investigations involving the necessary computational resources. However, qualitative peculiarities of these solutions should attract the researchers' attention to the nontrivial properties of the blow-up modes, at least, with respect to the obvious problem of internal TWR safety violation.

Conclusion
Below we give short conclusions stimulated by significant problems, that can be formulated in the following form.
1. The consequences of the anomalous 238 U and 239 Pu cross-sections behavior with temperature increase. It is shown that the capture and fission cross-sections of 238 U and 239 Pu manifest a monotonous growth in 1000-3000K range. Obviously, such anomalous temperature dependence of 238 U and 239 Pu cross-sections changes the neutron and heat kinetics drastically in the nuclear reactors, and in TWR in particular. It becomes crucial to know their influence on kinetics of heat transfer because it may become the cause of a positive feedback with neutron kinetics, which may lead not only to undesirable loss of solution stability (the nuclear burning wave), but also to a trivial reactor runaway with a subsequent nontrivial disaster.
2. Blow-up modes and the problem of nuclear burning wave stability. One of the causes of a possible fuel temperature growth may lie, for instance, in a deliberate or spontaneous coolant loss, analogous to what happened during the Fukushima nuclear accident. As shown above, the coolant loss may become a cause of the nonlinear heat source formation in the nuclear fuel, and consequently emergence of the mode with the temperature and neutron flux blow-up. In our opinion, the preliminary investigation of the heat transfer equation with nonlinear heat source points to an extremely important phenomenon of the anomalous of the temperature and neutron flux blow-up modes behavior. This result poses a natural nontrivial problem of fundamental nuclear burning wave stability and, orrespondingly, of a physically reasonable application of Lyapunov method to this problem, which is a basis for the motion stability theory, and thus a reliable basis for justification of the Lyapunov functional minimum existence. It is shown that some variants of solution stability loss are caused by anomalous nuclear fuel temperature evolution. They can be not only the cause of TWR internal safety loss, but can lead to a new stable mode when nuclear burning wave would periodically "experience" the so-called controlled blow-up regime through a bifurcation of states (which is very important!). At the same time, it is noted that such fast (blow-up regime) process has a number of physical vaguenesses, any of which may happen to be experimentally insurmountable for the control of such process.
3. On-line remote neutrino diagnostics of intra-reactor processes. Due to the fact that a high-power TWR or a nuclear fuel transmutation reactor are the projects of the single-load and fuel burn-up with a subsequent burial of the reactor apparatus, there is an obvious necessity for remote system of neutrino monitoring of the neutron-nuclear burning wave in normal operation mode and control of the neutron kinetics in emergency situation. The calculation of the spatio-temporal distribution of the isotope composition in the TWR active zone in the framework of the inverse problem of neutrino diagnostics of intra-reactor processes is presented in [27][28][29] in detail.