We consider a series of questions connected with composition and solution of kinetic equations for chain nuclear fission products into active zone of nuclear toroidal electric generator    . Of course, by spelling of kinetic equations for the reproduction of charged particles and charged splinters of fission this consideration would be accompanied by great share of the convention. The main aim, which here is placed, consists in the description of general regularity of reproduction process and not at all the counting of exact charges quantity up to last electron and proton.
It is important to note, in our view, that on dynamics of growth of charged particles Q number the governing action the growth of number prompt and delayed neutrons N into active zone toroid will be exercising so far as the each new fission act by neutrons will be leading to the new fission of nucleuses and to the formation of another charged particles and charged splinters.
In Sections 2-4, the conclusion of charged kinetic equations for electric generator is kept. Besides the finding equations there are established exponential instability of charged equations solutions for the coefficient of ionites reproduction which over one.
Investigation of charged kinetic equations is going on in Sections 5-9 also when the initial kinetic system is criticized by the influence of small random disturbances. Analysis of stochastic motion properties succeeded come to known scheme with action functional and to finding of corresponding kinetic quasipotential. Initial problem hasn’t analytical solutions in the class of free stable kinetic matrixes. Therefore there is proposed replacement method of the initial problem on the similar problem of kinetics with stable normal matrix.
2. Charge Kinetics of a Nuclear Electric Generator
Scaled demands at wide application of electric energy dictate the necessity to produce nuclear energy by means of fast-neutron electric generator. These nuclear devices are the most favourable due to their high values of the reproduction coefficient.
However, the use of fast-neutron electric generators challenges as there are several demands to maintain the nuclear kinetic process at a fitting level. Above all, it is necessary to keep high average energy of neutrons at the level of several hundred kilovolts   . Here the energy spectrum of neutrons is defined by the properties of combustible fuel, raw materials and construction materials in relation to the inelastic scattering of the neutrons (238U has a big cross-section of inelastic scattering). In addition, to achieve criticality the operation of fast-neutron nuclear devices needs a high concentration of a combustible fuel in relatively small volumes. The latter is caused by low values of the fission cross-section in the area of high energies.
Therefore, our further reasoning is based on the admission that the nuclear device (nuclear electric generator or nuclegen) has a breeder basis.
For simplicity and convenience all charged splinters and particles of both signs (ions, electrons, protons etc.) are called ionites, and their aggregate in quantity Q is ionite gas (ionite gas cord).
As in the case of neutrons, the balance of ionites is characterise by the ionite reproduction coefficient f. The reproduction coefficient f is defined as the ratio of the number of ionites in any generation to the number of ionites in the previous generation.
The average lifetime of one generation of ionites r is supposed to be proportional to the average lifetime of one generation of neutrons :
in which is a proportionality coefficient, is a set of parameters characterizing the electromagnetic influence, air possible residual drag of the medium, dynamic and construction features respectively, all influencing the value r.
For all intents and purposes the dependence (1) may possess not a constant but a rather complex function form with all values included being dependent on time and spatial variables in some way. In work  to analize the simplest charge kinetics value is fixed so that value r coincides with value l and is entirely defined by it. Apparently, this does nothold true what is evidenced by the following the correlation (1).
Unlike the neutron losses (capture without fission, leakage etc.), the losses of ionites, that control their lifetime, are determined by their neutralization degree (i. e. missed charges capture). The neutralization degree depends on the general field , particles velocity etc.
If , nuclear fission with the formation of ionites and neutrons occurs steadily at constant reaction rate ( , k is the coefficient of neutrons reproduction); its active zone is critical with a self-sustained chain reaction. This case implies , elements of inner electromagnetic field including are also close to constant values.
Ignoring delayed ionites the density change of ionite gas q in active zone within one generation equals to . If r is an average life time of one ionite generation, it means that per time unit the density of charged lively particles q changes in times, i.e.
where , with q0 is initial ionite density, is charge period of nuclear electric generator.
3. Linear Theory of Charge Kinetics
As mentioned above delayed neutrons play a crucial role and have a practical application in the controlled process of nuclear chain reaction. The fact that delayed neutrons alon with other β− descendants contribute into ionite gas formation is obvious.
Let γ is delayed ionite part (charged descendants), is instant ionite part. The reproduction coefficient f can be noted as a coefficient consisting of the two parts: , where , . After β−-disintegration decay of splinters the delayed ionites become instant. this means that their effective lifetime equals the sum: , where is delay time (β−-disintegration time), is instant ionite life time.
Hence, average effective lifetime of r ionite generation equals
Since , the charge period of nuclear electric generator is determined by the average time of ionite delay.
The simplest linear theory of charge kinetics needs more thorough consideration. The density of ionite gas q is admitted to be the blend of densities of charged particles and fission splinters of the two types mentioned above, i.e. . Here qm is instant ionite density, qz is delayed ionite density. Then instead of the correlation (2) ignoring the presence of delayed ionites the system of two kinetic equations is possible to be written. These equations describe how the dynamics of ionite gas density (concentration) changes in an nuclear electric generator:
with initial conditions
The right parts of Equations (3) and (4) should be commented. In equation (3) coefficient qm characterizes the density change of instant ionite number in lifetime of one generation. The second summand characterizes the growth of value qm at the cost of qz by the transition of the latter into the category of instant ionites. Equation (4) describes the density change of delayed ionites. In this equation the first summand shows the increase of the concentration of delayed ionites at the cost of formation instant ionites; the second one, on the contrary, shows the decrease of the concentration of the amount of delayed ionites within their lifetime (i.e. within the time of β−-decay). As well seen in the system (3), (4) the density (the concentration of the amount) of delayed ionites qz creates the feedback in relation to the density of instant ionites qm.
Following some simple transformations, the system of Equations (3), (4) is reduced to a homogenous differential equation of the second order with constant coefficients. The equation describes how to regulate kinetics peculiarities of a nuclear electric generator:
in which (or ),
The characteristic equation has two real roots of a clearly determined sign :
Hence, considering the starting conditions (5) the kinetics equations have the following solutions as functions of time:
In solutions (6) for qm and qz the second summand in the right part quicly approaches to zero. Disregarding it, growth dynamics of density (of the number) of ionites qm and qz is supposedly defined by the first summand in the solutions (6). As mentioned above, delayed ionites slow the transient process significantly, and the charge period is almost entirely defined by the average time of delay . This feature is also notable in the correlations (6). Indeed, the following approximations
prove the evaluations
4. Some Preliminary Results
Thus, a rough kinetics analysis of nuclear chain fission in vacuum toroidal shaped electric generator has shown that the achievement of exponential density growth (number) of charged particles moving at high rate and charged splinters may be considered proven. The boost of ionite gas density results in hasty growth of induced of electrical current in the external winding of a nuclear electric generator. It happens according to the electrodynamics peculiarities mentioned above (see the author’s article  , correlations (17)).
In agreement with the principals nuclear electrodynamics the solutions to modified Maxwell’s equations (the work  , formulae (7), (8)) for a general electric field can be found in the form
where -constant positive coefficients, . If the Fourier transform of the distribution of field is known at the initial moment of time
along with the full system of mutually orthogonal vector functions , then a general solution to the problem of unlimited growth of electromagnetic field can be presented using Equation (7) in the type of the following divergent series:
Currently, the possibilities of nuclear power engineering are restricted by chemical and thermal transformations at a significant low level. At the same time physics of nuclear electric engineering today is boosting the amount of neutrons fission in a reactor core as it inevitably results in overheating and thermal explosion.
This paper describes an alternative scheme for vacuum conducting of a chain nuclear reaction in quickly increasing and directional electromagnetic field. Under considered conditions there is nothing to heat except for cylindrical walls of a toroid. However, this problem is solved due to the presence of directional electromagnetic field that screens cascade currents of neutrons and ionite gas inside a nuclear generator.
5. Charge Kinetics under the Action of Small Random Disturbances
Here the problem of stochastic stability under the action of small random disturbances of initial kinetic process is going to be considered     . Stochastic fluctuations in a dynamic system behavior (3), (4) may vary by nature. They may result from the deformation of outer and inner electromagnetic fields, changes of thermal and vacuum regimes, chemical content of active nuclear components of kinetic process, constructive imperfections or operation inertance of control devices. The diversity of reasons causing slight, negligible disturbances to appear can lead to one general catastrophic overall result implying the destruction of the whole kinetic system itself.
The investigation of stochastic stability (instability) of the systems (3), (4) regarding small random disturbances needs revising some basic concepts.
6. Asymptotic Stability and Stochastic Instability
Unperturbed deterministic kinetic system (3), (4) is considered here:
which at the origin of coordinates (at zero point identifiable with a two-dimensional zero) has the stability position: .
The stability position (zero point) is stable (Lyapunov stable) if for any neighborhood there is a neighborhood for which equation solutions are (8) . If, moreover, at , then the stability position is called asymptotically stable.
Then it is assumed that the system (8) is perturbed by small white noise
where is a two-dimensional Wiener process, is a small numerical parameter. The problem of trajectories exit of process from limited region D with stability position is of great interest. The appearance of a small parameter in system (9) can by explained with an assumption that the noise intensity is low compared to deterministic factors affecting the system dynamics.
If matrix A is stable (Hurwitz matrix), i.e. its own numbers and have negative real parts. For this condition it is necessary and satisfactory  that coefficients of characteristical equation fulfiled the following inequalities
It is obvious that the inequalities are held if . As matrix A is considered as stable, the stability position of the unperturbed system (origin of coordinates) is asymptotically stable:
It is worth mentioning that vector function is not norm limited. Consequently solutions to Equation (9) may leave from any neighborhood of the stability position. Even with the presence of asymptotic stability of this stability position  , i.e. trajectories with probability 1 may deviate too far from the stability position. Thus, the system becomes unstable. To make the term “instability” meaningful some concepts are to be introduced using definitions and research schemes from  .
Denote the region is region if state vector , then object (9) functions; if the phase vector leaves D, then the object destroys. The region D is called critical region.
Introduce random value is the time before system (9) destruction, where is its solution with initial condition or in other way τ is the first moment of process exit from D. Probability value can be taken as an instability measure, if time interval of the object functional ability is fixed. Otherwise, mathematical expectation value of can be taken as an instability measure, if the time interval is unknown. The estimations of are given with the aid of action functional.
7. Action Functional and Quasipotential
Denote by a set of continuous functions within time interval with values in . In this space the metric should be considered and functional for absolutely continuous functions should be defined
If function is not absolutely continuous within or the integral (10) is divergent, then we suppose .
For example, in the case of classical mechanical system for which movement along trajectory is a functional motion of this system occurs along the extremals of functional . For a unperturbed dynamic system have
it means that if , then function in is the trajectory of system (11) as is absolutely continuous within and satisfies almost everywhere Equation (11).
Let the D is a limited in (critical) region and is its smooth boundary. If trajectory of system (11) starting at point leaves in a finite time, trajectories of process :
coming from also with probability close to 1 with small ε leave region D within this time. The first exit is more likely to occur near the exit point from D trajectory of system (11).
More precise definition of this result is connected to the approximation that happens while calculating probabilities of events :
with small for functional (10) for time interval . Functional with constant norming is stated as action functional for a set of random processes by analogy with quantum mechanics problems  and described with stochastic Equation (12).
If point is an asymptotically stable stability position of system (11), the quasipotential of system (11) relative to point 0 is function determined by equality
The functional of action is given by Equality (10) with
and the ends of the interval are not fixed. Have: function and is continuous.
The essential property of the quasipotential is given the statement concluding the following. If vector field in system (11) allows for decomposition
with gradient operator . Here potential function is continuously differential in , , , when . This means:
1) quasipotential equals to
for all so that ;
2) the only extremal of functional at a set of functions is defined by equation
It can be shown that for any curve connecting points 0 and the following un equation holds
i.e. and function as the only solution of Equation (15) is the only extremal.
In addition, it is possible to show that asymptotics of average time of process exit (12) from critical region D is determined by limiting equality
where is found quasipotential of dynamic system (11) relative to stability position (point 0) which is assumed as asymptotically stable.
8. Kinetic Quasipotential
The initial kinetic system (8) with a known stable constant matrix of coefficients A is worth mentioning again. Let the vector field permits the decomposition of type (14):
in which ; are non-zero square matrices of dimension 2 to defined. The matrix B is assumed to be a positive definite matrix. Considering the results of the previous section, there appears a task to find matrices B and C, potential , quasipotential and extremal equation of the type (15).
From the correlation (16) there is an equation
and from the condition of orthogonality , there is obtained
in which is a matrix transposed in relation to C.
The equations system (17), (18) leads to a square matrix equation in , which has not analytical solution for an arbitrary value of matrix A (i.e. a solution that does not lead to ). It is worth mentioning it is possible to arrive to the same equation from the Jacobi equation for a quasipotential under the condition of orthogonality of vectors and , in which , :
If A is a normal matrix (that is, ) with a negative definite symmetrical matrix , a solution can be found and take the form
The extremal equation also belongs to the same case when the A being a normal matrix:
To ensure the condition of orthogonality it is necessary to place a enough strict restriction upon matrix A as its normality condition:
The matrix A is not normal in the initial task. Consequently, we arrive to the situation which has a numerical solution in the case to be considered. A more optimistic prognosis is connected to an approximation task, a theoretical analogue of an approximate numerical method in a sense. It deals with a substitution of initial system (8) for a kinetic model which is asymptotically close to it.
9. Approximative Kinetic System
At first, matrix A cannot be a normal one even under the condition with the reproduction coefficient of ionites . As said above, in this case matrix A is stable. Normality of A means its symmetry, i.e. the resulting in the equality of elements :
in which it is designated as .
An important thing is that value can be given the “status” of a small parameter. Indeed, the previous appraisements imply that
Therefore, matrix A is to be approximated by a stable negative definite normal matrix :
where is a small parameter.
The choice of algorithm of approach (19) is to be the following. Let there be coefficient , more precisely in a normal matrix :
The obtained equations are the following
It is obvious that is a small parameter because is a small parameter here.
Thus, initial kinetic system (8) with a “small” disturbing vector for (in a limited time interval is provided for ) of the type
is given approximation by a “normal” stable kinetic system
Further, the scheme described above starts to operate for studying stochastic motion of system (20) disturbed with small white noise
with the aid of quasipotential theory. For system (20) quasipotential relatively to stability position has the type
where is stable negatively defined normal (symmetrical) matrix. The equation for extremals has the form
This allows to draw the following conclusion. Markov trajectories of system (21) with the probability tending to 1 for first leave critical region D near of point , i.e. near of the extremal starting at this point: .
The moment to sum up has come. The main and principal conclusion is contained in following: in this work, the inference of kinetic equations describing the reproduction process for charged particles and charged splinters inside of nuclear electromagnetic generator is conducted. Moreover, the description of solutions behaviour of charged equations in the form of ionite gas motion is realized. Exact analytical questions of the stable stochastic motion of ionite gas under the influence of small white noise on the charge kinetics equations are investigated in detail. We may be confirmed in intention that considered problems allow more effectively to solve the tasks for planning and making of realistic practicable nuclegens in the future.
 Tertychny-Dauri, V.Yu. (2017) Solution of Nuclear Electrodynamics Equa-tions Taking into Account β--Charged Eradiation. Open Access Library Journal, 4, e3794. https://doi.org/10.4236/oalib.1103794