Patent Number: 048760579
Section: claims

1. Process for the control of a nuclear reactor having a core containing fissile material assemblies, the fission of the nuclei of said materials being brought about by interactions with neutrons and producing in turn neutrons, said nuclear reactor also having a means for regulating the neutron flux in the reactor core and the power given off by said reactor core, said regulating means comprising at least one assembly of control rods displaceable in said reactor core, said process being characterized in that: (A) iteratively there is a determination of the neutron fluxes .phi.(j,g) for each zone or mesh j 1.ltoreq.j.ltoreq.J of a group of meshes corresponding to the volume of the reactor core and for each velocity group g 1.ltoreq.g.ltoreq.G of a plurality of velocity groups for the neutrons; and the numbers of neutrons or sources S(j) emitted by each mesh j, 1.ltoreq.j.ltoreq.J, per volume and time unit, the determination of said values consisting of repeating the following sequence of operations until said values converge: (B) the regulating means is controlled as a function of the neutron fluxes .phi.(j,g) and powers P(j).  effective macroscopic diffusion section .SIGMA.s(j,g g') of the neutrons of the velocity group g, 1.ltoreq.g.ltoreq.G in the velocity group g', g'.noteq.g,  effective macroscopic absorption section .SIGMA.a(j,g) of the neutrons of velocity group g, 1.ltoreq.g.ltoreq.G,  effective macroscopic fission section .SIGMA.f(j,g) of the neutrons of velocity group g, 1.ltoreq.g.ltoreq.G.  effective macroscopic diffusion sections .SIGMA.s.sup.0 (j,g g'), whose values are the possible values for said core,  an effective macroscopic absorption section .SIGMA.a.sup.0 (j,g), whose value is a possible value for said core,  an effective macroscopic fission section .SIGMA.f.sup.0 (j,g) of value equal to zero. 2. Process according to claim 1, wherein the distribution of neutrons in the core is defined by all the following parameters for each j, 1.ltoreq.j.ltoreq.J: 3. Process according to claim 1, wherein the predetermined interaction probabilities of the neutrons in the core can be defined for each mesh j, 1.ltoreq.j.ltoreq.J and each velocity group g, 1.ltoreq.g.ltoreq.G by: 4. Process according to claim , wherein for each mesh j 1.ltoreq.j.ltoreq.J all the second neutron flux components .phi..sup.1 (j,g), 1.ltoreq.g.ltoreq.G is determined as the solution of the system with G linear equations: ##EQU21## 5. Process according to claim 1, wherein the first component .phi..sup.0 (j,g) are calculated in accordance with the relations: ##EQU22## in which v(j) is the volume of the mesh j. 6. Process according to claim 1, wherein the new sources NS(j) are calculated in accordance with the relations: ##EQU23## in which .nu. is the average number of new neutrons per volume unit in the core. 7. Process according to claim 3, wherein G=2, the neutron fluxes .phi.(j,1) and .phi.(j,2) for each mesh j, 1.ltoreq.j.ltoreq.J being calculated according to the relations: ##EQU24## 8. Process according to claim 7, wherein the new sources NS(j) are calculated according to the equation: EQU NS(j)=.nu...SIGMA.f(j,1)..phi.(j,1)+.nu...SIGMA.f(j,2)..phi.(j,2) 9. Process according to claim 3, wherein for each mesh j 1.ltoreq.j.ltoreq.J, all the second neutron flux components .phi..sup.1 (j,g).sup.(m) 1.ltoreq.g.ltoreq.G is determined iteratively, the values .phi..sup.1 (j,g).sup.(n,m) obtained at the nth iteration being calculated as the solution of the system with G linear equations: ##EQU25## in which m is a time index, the neutron fluxes being calculated at instants t=t.sub.0 +m..DELTA.t, n is the iteration index for the calculation of the neutron fluxes at each instant t, .phi.(j,g).sup.(m-1) is the neutron flux obtained at the time iteration m-1 and v.sub.g is the average velocity of the neutrons of the velocity group g. 10. Process according to claim 1, wherein the first neutron flux components .phi.(j,g) are calculated as a function of predetermined coupling matrixes [.psi.g.sup.PREC ] for the delayed neutrons, whereof each element .psi.g.sup.PREC (j,k) represents the mean flux in the mesh j, for the velocity group g associated with a source uniformly distributed in mesh k and emitting one neutron per second in accordance with a delayed neutron emission spectrum .chi..sub.g.sup.PREC. 11. Process according to claim 10, wherein the first neutron flux components are calculated in accordance with the equation: ##EQU26## in which m is a time index, the neutron fluxes being calculated at instants t=t.sub.0 +m..DELTA.t, n is the iteration index for the calculation of the neutron fluxes at each instant t and S.sub.k.sup.PREC is the source, in mesh k, linked with delayed neutron precursors. 12. Process according to claim 11, wherein the sources S.sub.k.sup.(n,m) and S.sub.k.sup.PREC(n,m) are calculated in accordance with the equations: ##EQU27## in which .beta. is the total number of delayed neutrons per fission neutron and ##EQU28## in which .lambda.i is the radioactive decay constant of the precursor and Ci(j).sup.(n,m) is the concentration of the precursor i in mesh j at instant t.sub.0 +m..DELTA.t and at iteration n. 13. Process according to claim 12, wherein the new sources are calculated in accordance with equation: ##EQU29## 14. Process according to claim 1, wherein the neutron fluxes .phi.(j,g) and sources S(j) being iteratively calculated, said neutron fluxes are initialized to the value zero and said sources to a constant non-zero value. 15. Process according to claim 1, wherein the neutron fluxes .phi.(j,g) and sources S(j) are iteratively calculated at an instant t, the said neutron fluxes and said sources being initialized with the values obtained by said process at a preceding instant t-1.