SOVIET ATOMIC ENERGY VOLUME 18, NUMBER 1

Document Type: 
Collection: 
Document Number (FOIA) /ESDN (CREST): 
CIA-RDP10-02196R000700010001-6
Release Decision: 
RIFPUB
Original Classification: 
K
Document Page Count: 
118
Document Creation Date: 
January 4, 2017
Document Release Date: 
September 24, 2013
Sequence Number: 
1
Case Number: 
Publication Date: 
January 1, 1965
Content Type: 
REPORT
File: 
AttachmentSize
PDF icon CIA-RDP10-02196R000700010001-6.pdf7.09 MB
Body: 
Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 4t) Volume 18, Number 1 January, 1965 SOVIET ATOMIC ENERGY ATOMHAFI 3HEPrl1fl (ATOMNAYA iNERGIYA) TRANSLATED FROM RUSSIAN CONSULTANTS BUREAU Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 NEW BOOKS I CONSULTANTS BUREAO/T PLENUM PRESS REVIEWS OF PLASMA PHYSICS A systematic, multi-volume review of the present status of plasma theory, serving.both as an introduction for students and for re- searchers entering the field, and _a's a convenient, 'authoritative, up-to-date presentation. of current knowledge for workers in plasma physics. This continuing series, prepared by interna- tionally known Soviet experts in specific, fields, is under the editorship of Academician M. A. Leontovich, of the Xurchatovin stitute of Atomic,Energy. Each volume contains a number of integrated tutorial reviews, covering in depth' and in breadth specific aspects of ithe theory of the given field of plasma physics. In many cases, new material is presented. Transfaled by Herbert Lashinsky, University of Maryland. Volume 1: A comprehensive introduction ?to "classical" Flama physics, contains authoritative papers on: Orbit Theory, by D. V. Sivukhin; Collisions in Fully Ionized Plasma, by B.A. Trubnikov; Plasma Transport- Phenomena, by, S. I. Braginskii; and Plasma Thermodynamics, by A. A. Vedenov. Much of the material in the first two papers is presented here for the first time. Although the theoretical analyses are quite advanced, the experimental aspects of the subject are kept firpily in view throughout. This is especially .true of ? the article on transport phenomena; in which the kinetic approach is developed in derallel with qualita- tive physical descriptions of transport phenomena, including some of the less familiar "transverse" thermal transport effects plasma in magnetic .fields. Many physical examples' and'ap- plications of the theory are given. _ 336 pages ' ($12.50 Volume 2: Contains four review papers concerned primarily with the problem of /plasma confinement: Magnetic Field .Ge,, ometries; Plasma Equilibrium in Magnetic Fields; liydromagnetic Plasma Stability; and Motion of Charged Particles in Electro- ' manetic Fields.. ' Approx. ,300 pages $12.50 Volume 3: Devoted to plasma waves, includes: Electromagnetic Waves in a Plasma; Oscillations of an Inhomogeneous Plasma Introduction to the Theory of a Weakly Turbulent Plasma; and Symmetric Magnetohydromagnetic Flow and'Helical Waves in a Circular Plasma Cylinder. ' Approx. 300 pages ? $12.50 Volume 4: Contains, three papers: Hydrodynamic Descriptions of a Highly Rarefield Plasma; Collective Phedomena in Collision- less Plasma Shock Waves, and Coulomb Collisions in Fully loniied Plasma. The latter paper contains new material on the relevance of this topic in mirror machines, provided 'by the author (D. B. Sivukhin) for the English edition. Approx. 290 pages '$12.50 Volume 5: To be published late in 1965 in the Soviet Union is to comprise a comprehensive review of radiation phenomena in plasma. An English translation will be published promptly. LEBEDEV PHYSICS SERIES Complete English translations of the Proceedings'("Trudy") of the 'ferried Lebedev Physics Institute of the USSR Academy of Sciences published as Consultants Bureiu Special Research Re'- pits. Series edited by D. V. Skobel'tsyn. The fi6t three 'volumes in this distinguished series of ten Special Research Reports in translation, beginning with Volume 25, are: ? ' ? OPTICAL METHODS,' - OF INVESTIGATING SOLID BODIES "Trudy" Volume 25 Includes papers by N. D. Zhevandrpv (on polarized luminescence', of crystals), V. P. Cheremisinov (on vibrational spectra and struc- ture Of oxides in the Crystalline and glassy states) and L. A. Vainshtein on calculation of cross sections for excitation of .atoms and ions 'by electron impact). ' 194 pages 1965 '$22.50 COSMIC RAYS "Trudy" Volume 26 Contains an account of experimental investigations into nuclear , and 'electromagnetic interactions at high and ultra-high energies. Theoretical 'articles on problems of The passage' of high-energy ' ? ,electrons and photons through matter are included, as well as ,some results of the investigations into cosmic rays by means of artificel earth satellites. ? 2,62 pages 1965 . ? $27.50 RESEARCH'IN MOLECULAR SPECTROSCOPY - "Trudy" Volume,27 Devoted to spectroscopic' investigations into matter in various [states of aggregation by' the methods of Raman scattering 'and infrared absorption. A special section is devoted to the metho: dological problem of Correcting-'measured quantities forInstrit: mental errors. Great ,attention is' paid to the investigation 'of -- the temperature dependence of sdch spectroscopic quantilie:i as , the profile, width, and integrated intensity of Raman lines and absorption bands. 214 pages 1965 $22.50 CONSULTANTS 'BUREAU! PLENUM PRESS 227 ly.,17th St., New York, N.,Y: 10011 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 ATOMNAYA ikERGIYA EDITORIAL BOARD A. I. Alikhanov A. A. Bochvar N. A. DollezhaP K. E. Erglis V. S. Fursov I. N. Golovin V. F. Kalinin N. A. Kolokol'tsov (Assistant Editor) A. IC Krasin I. F. Kvartskhava A. V. Lebedinskii A. I. Leipunskii M. G. Meshcheryakov M. D. Millionshchikov (Editor-in-Chief) I. I. Novikov V. B. Shevchenko A. P. Vinogradov N. A. Vlasov (Assistant Editor) M. V. Yakutovich A. P. Zefirov SOVIET ATOMIC ENERGY A translation of ATOMNAYA iNERGIYA A publication of the Academy of Sciences of the USSR @ 1966 CONSULTANTS BUREAU ENTERPRISES, INC. 227 West 17th Street, New York, N.Y. 10011 Volume 18, Number 1 January, 1965 CONTENTS PA ENG. I GE RUSS. The Seventieth Birthday of Academician A. L. Mints 1 3 High-Frequency Oscillations Excited on Interaction of an Electron Beam with Plasma ?A. K. Berezin, Ya. B. Fainberg, L. I. Bolotin, and G. P. Berezina Interaction of Plasmoids with an Electromagnetic Wave?V. I. Veksler, I. R. Gekker, E. Ya. Gol'ts, G. A. Delone, B. P. Kononov, 0. V. Kudrevatova, G. S. Luk'yanchikov, M. S. Rabinovich, M. M. Savchenko, K. A. Sarksyan, K. F. Sergeichev, V. A. Silin, and L. E. Tsopp 12 14 Determining the Perturbations of the Parameters in the Magnetic and Accelerating Systems of an Electron Synchrotron on the Basis of an Analysis of Information Regarding the Beam?I. P. Karabekov 17 18 Phase Stability of a System of Particles in Self-Regulated Accelerators? E. A. Zhil'kov and A. N. Lebedev 22 22 Measurement of the Photoneutron Yield from Thick Copper and Water Targets and Determination of the Excitation Function of the (y, n) Reaction for 016 and Cu63 by Means of the Belen'kii-Tamm Equilibrium Photon Spectrum?I. A. Grishaev, D. I. Sikora, V. A. Shkoda-Ul'yanov, and B. I. Shramenko 29 28 Transient Processes and the Measurement of Reactivity of a Reactor Containing Beryllium ?S. S. Lomakin and Yu. A. Nechaev 35 33 ' The Crystal Hydrate UF4.4/3 H2O?Yu. V. Gagarinskii, E. I. Khanaev, N. P. Galkin, L. A. Anan'eva and S. P. Gabuda 43 40 The Relative Volatility of Solutions of HTO in H2O?Ya. D. Zel'venskii, V. A. Shalygin, V. S. Tatarinskii, and D. A. Nikolaev 49 46 Determination of y-Ray and Neutron Absorbed Dose in Polymers?F. A. Makhlis and I. M. Kolpakov 52 48 Increasing the Depth of Prospecting for Concealed Uranium Ore Bodies by Means of the Primary Aureole?S. V. Grigoryan 57 52 LETTERS TO THE EDITOR Phase Stability of Particle Blobs in Accelerators with Automatic Control?E. A. Zhil'kov . . 62 58 Simple Method for Measuring the Frequency of Free Transverse Oscillations in Cyclotrons ?S. A. Kheifets and S. K. Esin 65 60 Nomograms for Determining the Potential Barrier's Height and for the Breit?Wigner Formula ?G. N. Potetyunko 67 61 Annual Subscription: $ 95 Single Issue: $30 Single Article: $15 All rights reserved. No article contained herein may be reproduced for any purpose what- ? soever without permission of the publisher. Permission may be obtained from Consultants Bureau Enterprises, Inc., 227 West 17th Street. New York City, United States of America. Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 CONTENTS (continued) P AGE EN6. I RUSS. Simple Unsteady-State Kinetic Equation?V. G. Morogov and S. A. Kholin 70 62 Angular Distribution of y-Quanta in U233, U235, and PO" Fission on Thermal Neutrons G. A: Petrov, D. M. Kaminker, G. V. Varskii, and L. A. Popeko 72 64 Determination of the Absolute Yield of the 74-keV U239 and 87-keV Th 233 y-Lines N. Yurova and A. V. Bushilev 75 65 Spatial Distribution of Neutrons with Energies of 3 and 15 MeV in Beryllium?S. P. Belov, V. A. Minn, Yu. A. Kazanskii, and S. G. Tsypin 78 67 Reduction of the Capture y-Radiation from the Reactor's Structural Materials by Screening Them with Boron-Containing Screens?B. P, Gromov, D. V. Pankratov, M. A. Solodyankin, and M. M4 Sokolov 80 69 Dependence of the Density of RadiatiOn Damage to the Reactor Vessel on the Composition of the Ferro-Aqueous Thermal Shield?K. K4 Popkov and S. M. Rubanov 83 70 Antifriction Characteristics of Neutron-Irradiated Steel?E. A. Markovskii and M. M. Krasnoshchekov 85 72 Method of Measuring Radioactive Preparations and Checking Stability?V. M. Malykhin . . . 87 73 The Role of Thermal Peaks in the FOrmation of Defects?L. G. Gurvich and N. S. Bespalova 91 76 SCIENCE AND ENGINEERING NEWS Agreement on Collaboration in Desalinization Efforts 94 78 International Conference on the Quantiim Theory of Systems Having Many Degrees of Freedom?P. S. Isaev - 96 79 All-Union Conference on Nuclear Meteorology?S. G. Malakhov and I. V. Yagodovskii . . . 99 80 Application of Methods of Nuclear Geophysics in Ore Prospecting, Exploration, and Development?S. I. Savosin and V. I. Sinitsyn 100 81 Radioactive Chlorine-36 in Monitoring the Production and Processing Of Hexachloran ?G. M4 Strongin and M. N. Kulikova ifiBLIOGRAPHY 104 84 New Books. 106 8S The Ibissian date "Podpisano k pechati" of this issue was 12/19/1964 . This is equivalent to "approved for printing." Publication did not occur prior to this date, but must be assigned to have taken place reasonably soon thereafter. Publisher Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 THE SEVENTIETH BIRTHDAY OF ACADEMICIAN A. L. Mints Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 3-4, January, 1965 Aleksandr L'vovich Minch January 8 marked the passage of 70 years since the birth of Hero of Socialist Labor, Academician Aleksandr L'vovich Mints? the Direction of the Radio Engineering Institute, Academy of Sciences of the USSR Alexandr L'vovich Mints is widely known in our country and abroad as a prominent scientist and engineer, talented or- ganizer, and director of the development and installation of high- power broadcasting stations and gigantic charged particle ac- celerators. In the Civil War, Aleksandr L'vovich was an active fighter for the victory of Soviet power, taking part in the battles in the Caucasus, on the Polish and Crimean fronts as a member of the legendary first mounted army, in the capacity of the Radio Divi- sion. After the end of the Civil War, A. L. Mints continued serving for some time in the Red Army, occupying a position of command in the advanced communications War College in the Scientific Research Communications Institute, and then, around 1924, he turned to the development and installation of high- power radio station in the Soviet Union, the installation of which to some extent or other was not connected with the name of Alexandr L'vovich. At that time, he was responsible for the de- sign of new high-power radio stations, and was in charge of construction. A. L. Mints developed and built dozens of radio stations and radio centers in the Soviet Union, starting with the 20 kilowatt ASPopov radio station (1925- 1927), and ending with the superhigh-power 1200-kilowatt radio station (1941-1943). This period clearly exhibited one of the basic characteristics of A. L. Mints ?the effort to bring his scientific developments to practical applica- non in the shortest possible time. The scientific work of A. L. Mints embraces literally all fields of high-power radio equipment: the theory and development of methods for designing radio telephone modulation systems, the development, investigation, and installation of new antenna systems for long and short-wave radio stations, the de- velopment of new vacuum tube devices, and much, much else. In 1930-1932, A. L. Mints proposed and brought to practical realization a completely new unit system, which has been used up to the present time in building all the large radio stations in the USSR as well as abroad. To a considerable extent, as the result of A. L. Mints' work, Soviet high-power radio work occupied a leading place in the world even in those years. The advanced ideas which A. L. Mints carried out in high-power radio stations exerted an enormous effect on the development of Soviet radio engineering, as well as the radio industry. Building radio stations and radio centers revealed the brilliant organizational abilities of A. L. Mints, who showed in practice that the time required to build and adjust large radio engineering installations may be greatly reduced by a correct and exact organization of the work. The creative and organizational abilities of Alexandr L'vovich were brilliantly revealed in the period which marked the .beginning of construction of the first large charged particle accelerators in the USSR In 1946, a group of specialists headed by A. L. Mints, M. G. Meshcheryakov, and D. V. Efremov took over the construction of the first high-power proton accelerator in the USSR ? the 500-MeV synchrocyclotron in Dubna (later its energy was raised 1 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 to 680 MeV). The accelerator was built in an unusually short time, and began to operate in 1949. Then, likewise in Dubna, a group from the radio engineering laboratory of the Academy of Sciences of the USSR, together with other institutes, engaged in the development and construction of what was at that time the largest proton accelera- tor the 10-BeV synchrophasotrOn, which was started up in 1957. Building acceleratorS required the solution of very difficult radio engineering problems: In the synchrocyclotron, the frequency of a high-power oscillator had to in- crease by a factor of two during a very short oscillation cycle (3 ft-Bee), While in the synchrophasotrori, the frequency had to be changed by as much as a factor of 10 according to a definite law, and at the same time the value of the frequency had to be maintained with an accuracy of 3 ? 10-4 at all times. These Very complicated problems were solved, not only successfully and elegantly from the engineering point of view, but with a high degree of reliability. In later years, the range of A. L. Mints' scientific interests continued to increase: In addition to radio electron- ic systems for the ring-type stiff focusing 7-BeV proton accelerator (started up in 1961 at the Institute of Theoretical and Experimental Physics) and the 70-BeV one (installed in the Institute of High-Energy Physics of the GKAt), Alexandr' L'vovich directed the development and building of the 30-MeV linear electron accelerator at the I. V. Kurchatov Atomic Energy Institute, the 20-MeV and 100-MeV linear proton accelerators, and other machines. Since 1962, development of a project for building a gigantic ring-type proton accelerator for energies of 1000 BeV has been going on under A. L. Mints' direction. This project includes new ideas among them, ideas for the cybernetization of the acceleration processes so-that the accelerator has been given the name cybernetic. At the same time a Model of this accelerator is being built for an energy of 1 BeV. .Since 1958, at the suggestion of-academician I. V. Kurchatov, a group from the Radio Engineering Institute of the Academy of Sciences,USSR, headed by A. L. Mints, has Started work on high-frequency plasma equipMent. In 1959, work began on the investigation of the ionosphere, of Space near the earth, and Of interplanetary plasma. The years that Academician A. L. Mints spent On building accelerator's were also the years of the ,founding and development of the Radio Engineering Institute of the Academy of Sciences, USSR, which grew from a small laboratory into a large scientific research center. A. L. Mints was the organizer of this center, and is its per- manent director. To the comparatively young group at the institute, A: L. Mints communicates his character- istic Scientific boldness, persistence in the achievement of a goal, and clearly defined organization, as Well as organic Unity in profound theoretical analysis of the most complex questions of present day radio engineering and experimental physics and the practical realization of Scientific ideas. A. L. Mints has educated many talented young scientists, known for their work, among them not a few State and Lenin prize winners, candidates and doctors of science. During the 45 years of his scientific, engineering, and social activity, A. L. Mints published more than 100 scientific and engineering science papers, books and pamphlets, made more than 30 inventions in various fields of radio engineering and electronics, and has presented more than 70 papers, popular scientific, and publicizing ar- ticles. A. L. Mints carries on a large amount of social work, and is a member of the Scientific Councils of several Institutes, a member of the Bureau of the Division of General and Applied Physics, Academy of Sciences of the USSR, as well as an honorary member of the ASPOpov Scientific and Technological Radio Engineering and Elec- trical Communications Society. The Scientific and Engineering accomplishments of A. L. Mints are highly esteemed by the Communist Party of the Soviet Union, and by the Soviet Government. He has been awarded three orders of Lenin, two Orders of The Red Banner of Labor, the Order of the Red Star, and medals. He has received the high distinction of Hero of Socialist Labor, and has been distinguished by two state prizes of the first grade, and by the Lenin prize, The Presidi- um of the Academy of Sciences of the USSR has awarded A. L. Mints the ASPopov Gold Medal. The high devotion to principle, ability as an organizer, profound feeling of responsibility, and great self criticism, scientific boldness and ability to act, harmoniously blended in Alexandr t'vovich with a fine and spritual relation to his fellow workers, have given him great and well deserved authority. Academician A. L. Mints is a splendid example of a prominent Soviet Scientist, Engineer and Organizer. 2 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 HIGH-FREQUENCY OSCILLATIONS EXCITED ON INTERACTION OF AN ELECTRON BEAM WITH PLASMA UDC 533.9 A. K. Berezin, Ya. B. Fainberg, L. I. Bolotin, and G. P. Berezina Translated from Atomnaya tnergiya, Vol. 18, No. 1, pp. 5-14, January, 1965 Original article submitted June 4, 1964 Some results of experiments on the observation and study of oscillations excited in a beam and plasma as a result of their interaction are presented. The experiments were made under condi- tions in which wo < cuil, where wo is the electron Langmuir frequency of the plasma, and WH is the electron cyclotron frequency. The conditions of excitation of the waves in the plasma were determined, together with their frequency spectra, phase velocities, and gain factors, the inten- sity of the electric field, and the absolute values and spectral distribution of the power of the os- cillations excited. The experimental results for the frequencies, gain factors, and phase ve- locities of the oscillations excited in the plasma are in satisfactory agreement with calculated data. As shown earlier in [1-3], on passing a heavy-current pulsed electron beam (current 5 to 8.5 A, energy 15 keV, pulse length 3.6 ?sec) through a plasma in a longitudinal magnetic field (intensity 400 to 1320 Oe), under certain conditions it loses a considerable part of its initial energy (10 to 25%). In these experiments the plasma is created by the beam itself. It follows from theoretical considerations [4, 5] that this energy must pass into excitation of oscillations in the beam and plasma, and also into heating them. As we know [4-6], in the case described in the present paper, forward waves (normal dispersion) should be ex- cited. For large electron beam densities, frequencies associated with the electronic Langmuir oscillations of the beam must be excited together with Langmuir electronic frequencies of the plasma. The measurements were made under standing and traveling wave conditions for current 5 and 8.5 A in a longitudinal magnetic field of intensity 720 to 1320 G. In the first case, a metal reflector plate (current collector) was placed at the end of the chamber, and in the second an adiabatic absorbing load (graphite cone 12 cm long). Study of Oscillations Excited in the Plasma The block diagram of the apparatus is shown in Fig. 1. Let us examine some of the components in more detail. For recording the Hgo-, Hp -, Hz-components of the hf field, the oscillations are detected by a loop, and for recording the Ez component by a half-wave dipole situated in the plasma and oriented in the appropriate direc- tion, connected via a movable coaxial cable to a power-calibrated resonance wavemeter (band half-width ?1.5Mc). Changing the wavemeter over to frequency and recording the power variation of the signal detected, the fre- quency spectrum of the oscillations excited in the plasma could be obtained for a fixed position of the loop or di- pole inside the chamber. By moving the coaxial cable bearing the loop or dipole in the direction of motion of the electron beam and recording the variation in the hf power of the oscillations excited in the plasma for given fre- quency, the intensity distribution of the corresponding components of these hf oscillations could be obtained. In order to eliminate vibration, the coaxial cable bearing the loop or dipole was moved inside a guide (copper tube 4 x 6 mm, with a side slot along the whole length). The guide was fixed inside the chamber near to the edge, parallel to the axis of the system. The movable cable came out of the side slot in the guide and ended in the loop (area ?1 cm2) or half-wave dipole. The loops and dipole were placed at a distance of 10 mm from the beam bound- 3 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Beam at e ray oscine ra h. Fig. 1. Block diagram of apparatus for studying oscillations excited in the plasma under standing and traveling wave conditions: 1) Guide; 2) loop or dipole; 3) chamber; 4) fixed dipole; 5) current collector; 6) movable coaxial cable; 7) leak; 8) synchronous motor. ary. The coaxial cable passed through to the outside by a way of a special vacuum seal, and was connected via a calibrated attenuator to a power-calibrated resonance wavemeter. The signal from the output of the wavemeter fell simultaneously on to an oscillograph and a wide band amplifier, then to a cathode follower with integrating circuit, and then to an automatic electronic potentiometer EPP-09. The envelope of the hf oscillations was ob- served on the oscillograph screen and its amplitude measured (to an accuracy of 5%). Calibration of the wave- meter for given frequency (in an assigned frequency range) consisted of the following. Oscillations of known power from a standard signal generator operating under continuous conditions fell on the wavemeter with crystal detector, the steady voltage at the output being measured. The generator power was determined with a low-power meter (IMM-6) to an accuracy of 7 to 10% Such calibration was carried out before and after measuring the hf power of the oscillations excited. Only measurements for which the calibration did not alter during the experiment were taken into account. On measuring the power in this way, the total error was ?17 to 20%, the error associated with the calibra- tion of the wavemeter being ? 7 to 10%, and that in determining the attenuator factors being ?10% The magnetic field strength of the corresponding components of the hf oscillations were calculated from the formula fr7los P p ? (OS (1) where P is the hf power Collected by the loop in W, p is the input resistance of the wavemeter (p = 7512), S is the area of the loop in cm2 (S = 1 cm2); and u.) is the oscillation frequency in cps. The Most intense oscillations in the plasma were excited in the range 825 to 835 Mc (half-width 50 to 70 Mc) in the pressure range 4 to 7 ? 10" mm Hg, An example of such an oscillogram taken for a frequency of 825 Mc appears in Fig. 2. From measurements of the hf field distribution in standing wave conditions, the phase velocities were deter- mined; measurements under traveling wave conditions gave the gain coefficient of the excited oscillations. The relative phase velocity 3 was calculated as follows: where X and X0 are the wavelengths in the plasma and free space respectively. from the relation 1 in P2 V = 2 (z2?zi) Pi (z2?z1 In -112 .115 The gain factor (2) was calculated (3) Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Fig. 2. Variation of the amplitude of the hf oscilla- tions with time. Oscillation frequency 825 Mc; volt- age pulse length 3.6 ?sec; scale 1 ?sec. 100 SO 70 50 40 30 .111 20 W, rel. units 10 7 3 2 0 3 6' 9 72 15 1, cm a 100 80 70 50 40 30 20 10 7 3 18 21 24 where H1 and H2 are the field intensities of the propagat- ing wave at points z1 and z2, while P1 and P2 are the cor- responding hf powers. Graphs of the distribution of the Hp -components of the electromagnetic hf field along the direction of motion of the beam (z axis) at frequency 835 Mc under standing wave conditions are presented for various mag- netic field strengths (1320, 960, and 720 G) in [2], in Fig. 7a,b, and c. In these graphs, the axis of ordinates represents the hf power proportional to lizp (in relative units) collected by a loop suitably oriented and situated at a given point along z. The axis of abscissas represents the distance (in cm) reckoned from the end of the region of interaction. Figures 3a and b show distributions of the Hp -com- ponent of the hf field in the plasma along the direction of motion of the beam for both standing and traveling waves. Analogous graphs were obtained for the Ha-, and Ez-components of the hf field in the plasma. The slight peaks on the curve of Fig. 3a may be explained by the fact that, owing to the imperfect load- ing, together with absorption, there is a partial reflection of the wave (cone length 12 cm, i.e., approximately equal to the wavelength in the system at frequency 835 Mc). Measurements show that on an average the wave- length of the oscillations in the plasma is ?10 to 12 cm, i.e., 3 0.3 (oscillation frequency 835 Mc). Below we give the gain factors as a function of the 27 longitudinal magnetic field intensity (frequency 835 Mc): = IIIIIIIIIMINNIMII MIII !Ilii II I ii 1 Sill r or= 16 IAN in IN 11 Wail , 111 3 6' 9 12 15 18 21 24, 27 1, cm Fig. 3. Graphs showing the spatial distribution of the H -component of the hf field at frequency 835 Mc. The electron beam has current 5 A and en- ergy 15 keV; the pressure of the working gas is 6 ? 10-4 mm Hg, and the magnetic field intensity 1320 G; a) traveling waves; b) standing waves. The beam moves from right to left. Magnetic field intensity, G Gain factor, cm 1320 0.21 960 0.23 720 0.26 On changing the current and velocity of the beam, the gain factors also change, increasing with increasing beam current and decreasing beam velocity. Thus, in the present experiments, when an elec- tron beam interacts with a plasma, slow electromagnetic waves are excited in the latter, and these rise in inten- sity along the direction of motion of the beam. The mag- netic field intensity of these oscillations at the end of the chamber, calculated from formula (1), reaches 0.3 to 0.4 G. Variation of the longitudinal magnetic field in- tensity from 1320 to 720 G hardly affects the frequency of the oscillations, but alters their intensity substantially: as the magnetic field intensity rises that of the oscilla- tions falls. The frequency spectra were measured in the fre- quency range 400 to 3200 Mc at the end of the region of 5 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 f ,Mc 1400 1200 7000 800 500 400 7 I I 2 1320 Oe '4 ' I ? I I I I 10x10-1 f , Mc 1400 1200 1000 800 600 400 2 960 Oe J, rel. units J, rel. units 250- 250 200 200 150 - 150 _ 100 1011 tjO 50 0 i Ii 2 4 6 8 10x10-1 2 ' !III :z 8 10x10 Pressure, mm Hg a J, rel. units 250 200 150 100 50 I I IL 0 4 6 8 7010-4 Pressure, mm Hg Fig. 4. a) Excited frequencies II and 12 and spectral half-widths, and b) maximum intensity of ex- cited oscillations as functions of air pressure in the plasma chamber for a 5-A, 15-keV electron beam for various longitudinal magnetic field intensities. 2 4 8 8 10 x10 1- interaction of the beam with the plasma. For pressures below 8 ? 10-4 mm Hg, the frequency spectra have two (sometimes three) sharp maxima. For the high pressure region (8 ? 10-4 to 2 ? 10-3 mm Hg) there is a character- istic plateau in place of the maxima. From the data obtained, we constructed graphs on which the values of the excited frequencies, the spectral half-widths, and the maximum intensity of the excited oscillations are given as functions of air pressure in the plasma chamber, current (5 and 8.5 A), and longitudinal magnetic field intensity (720, 960, and 1320 G). The beam energy was 15 keV. The graphs for the 5 A current are shown in Fig. 4a and b. The graphs for the 8.5 A current are of the same form. The frequency ranges between the intensity maxima are shaded in this figure. The frequency of one of the intensity maxima (11) depends on the air pressure in the plasma chamber, increasing as this rises. On the other hand, the frequency of the second maximum (12) is independent of pressure, remaining almost constant; further, 11 > 12. We should note that f2 S20/27r, where 00 is the electronic Langmuir frequency of the beam. As seen from the graphs given in Fig. 4, the spectral half-widths of the excited frequencies rise on increasing the air pressure in the plasma chamber. This may evidently be explained by the fact that these half-widths are directly proportional to the frequency of collisions in the plasma, which increases on raising the pressure of the working gas. With increasing velocity of the electron beam, the frequency fj rises. Figure 4b shows the maximum intensity of the oscillations at frequency f as a function of the air pressure in the plasma chamber. From these graphs we may draw the following conclusions; 1) the intensity of the excited oscillations reaches a maximum at pressures 4 to 7 . 10-4 mm Hg; 2) on increasing or decreasing the pressure rela- 6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 ? Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 1 Current, A a A*, MC 520/2n, MC 5 0.53 450 560 8.5 0.44 600 730 tive to the optimum value the intensity of the excited oscillations falls; 3) on lowering the longitudinal mag- netic field intensity from 1320 to 720 G the intensity of the oscillations rises. Measurements also show that on raising the velocity of the electron beam the intensity of the excited oscillations falls. Similar laws were ob- served for experiments in both standing and traveling wave conditions. Since in our case the plasma density N is proportional to the air pressure in the chamber, we constructed graphs relating f to the plasma density for fixed current and energy of the electron beam. Such graphs appear in Fig. 5a and b for currents 5 and 8.5 A respectively. The values of plasma density as a function of air pressure in the chamber were taken from Fig. 3 of [2]. The plasma density was measured by two methods: by means of a hf interferometer (X 3.2 cm), and by means of a method based on the use of the dispersion properties of plasm-a waveguides. 14 12 10 \ 8 4 2 to fo 1 2 N, cm-3 a Cl U 2 /o 2 3 4 5 6 N, cm-3 Fig. 5. Variation of the square of the excited fre- quencies in the plasma (fl) as a function of plasma density for an electron beam of energy 15 keV and current (a) 5 and (b) 8.5 A: ED 1320 G; 0) 960 G; A) 720 G. mov As we know, there is a cutoff band in our wave- guide between the limits coo and (.0H [for the case WH > wo (see Fig. 9, curve ?o?)]. Hence after deter- mining the frequency corresponding to the beginning of the cutoff we may find the plasma density [2]. We may therefore suppose that the squares of the ex- cited frequencies (A) vary in direct proportion to the electron density of the plasma. The excited frequencies rise on increasing the current in the electron beam. The equation of the straight lines shown in Fig. 5a and b may be written in the form a21(2) (A*)2, (4) where a is the coefficient of proportionality between f and f, A* is the frequency of the excited oscilla- tions in the absence of plasma (fo = 0), and f = w0/27r. From the calculations made we obtain the values given in the table (remember that no is the electronic Lang- muire frequency of the beam). As seen from the table, the values Ai' = 450 Mc and A2* = 600 Mc are in agreement with the values of the Langmuir frequency for the electrons in the beam for the above current values. If we neglect a small correction in the investigated frequency region As to fi, the frequency of the excited oscillations fi is 0.53 f for a 5 A current and 0.44 f for an 8.5 A current. We also studied oscillations with frequencies in the 2400 Mc range. These oscillations were only ob- served under standing wave conditions. The measure- ments showed that these hf oscillations in the plasma had only Hz- and Hp -components. The intensity dis- tribution of the oscillations along the z axis at fre- quency 2405 Mc for current 5 A and energy 15 keV with magnetic field intensity 1320 G is shown in Fig. 6. As seen from the graph, the intensity variation of the oscillations along the z axis at frequency 2405 Mc dif- fers from the distributions found earlier., The radiation consists of narrow lines (two or three lines with a half- 7 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 I, rel. units 160 150 140 130 120 110 100 90 80 70 60 50 48 30 20 ? 10 %.? 0 ; 2 3 4 6. 6' Ju_ I 8 3 10 11 12 13 14 1516' 12L,cm Fig. 6. Distribution of the Hp -component of the hf field in the plasma along the axis of the system at frequency 2405 Mc. Beam energy 15 keV, current 5 A, longitudinal magnetic field intensity 1320 G. Direction of beam motion in the oscillogram right to left. width of the order of 3 to 5 Mc). The maximum field intensity of these oscillations is 0.05 to 0.15 G. The fre- quency of the oscillations is almost independent of the pressure of the working gas in the chamber and the magnetic field intensity, but depends greatly on the current and velocity of the electron beam, rising as these increase. Thus we may conclude that the excitation of os- cillations in the frequency region 2400 Mc depends on the boundary conditions, while in the range 825 to 835 Mc it is independent of these. The most distinctive feature of the oscillations at frequency 825 Mc is the fall in their intensity on lowering the magnetic field strength. Analogous measurements were made using argon and hydrogen as working gas instead of air. The de- pendence of Hp on z for air and argon was much the same if the working gas pressure in the chamber was roughly the same and other experimental conditions were identical. In order to obtain the same kind of os- cillations after admitting hydrogen to the chamber, its pressure had to be raised to 5 ? 10-3 mm Hg (i.e., to some six times the-value for air). These results maybe explained by noticing that the specific ionization co- efficients &(V) in air and argon are roughly the same, while for hydrogen the value is six or seven times smaller [7]. Study of Oscillations Excited in the Electron Beam After Passing through the Plasma The block diagram of the apparatus for studying these oscillations is shown in Fig. 7. The electron beam (not modulated at the input), after interacting with the plasma, passed through a helical junction, by means of which hf power in the beam was selected. 8 tEl rEl0tEl 00 o ZICZ1 I1,21/NNYMMMY11111 IVIVIVAIIIINI7A7AVAVNIAMVAltim G') 0000000 Pumping Cs4 Attenuator , i!meri CZ0 ShU-1 0 0 scillo- graph I Pumping Integrating CirCuit To oscillograph Fig. 7. Block diagram of the apparatus for studying oscillations excited in the beam: 1) Tube for creating pressure drop; 2) chamber; 3) inner helix of helical junction; 4) glass tube; 5) outer helix; 6) .current collector; 7) solenoid for creating longitudinal magnetic field; 8) leak. Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 It 3 c.) k 2 The helical junction was made in the following way. In- side a vacuum-sealed glass tube of internal diameter 18 mm was placed a helix of 1.5 mm diameter molybdenum wire. The pitch of the helix was chosen so that the phase velocity of a wave propagated along the helix should be equal to the velocity k- 7 of the electron beam passing inside the same helix. On the out- side, the helix was surrounded by a second, shorter helix with turns in the opposite direction. This latter was connected to a coaxial cable, the wave resistance of which was 75 O. At one end, the inner helix was connected to an electrode at zero ("ground") potential. The ends of the glass tube were covered with a layer of aquadag. Thus electrons falling on the ends of the glass tube and the inner helix spiral flowed to ground. 'The length of the helical junction was -15 cm. The standing wave voltage coefficients (SWVC) of the input into these helical junc- tiorfs-clid not exceed 2.0 in the frequency range 650 to 5000 Mc. These measurements were made on systems without an electron beam. Beyond the helical junction was placed a collector, and the electron beam fell into this. The experiments were mainly carried out with a 5-A electron beam of 15-keV energy at a pressure of 6 ? 10-4 mm Hg, and a longitudinal magnetic field intensity of 1320 G. Measurements showed that the spectra of the frequencies excited in the beam were analogous to the earlier described spectra of oscillations excited in the plasma under the same conditions. The power of the oscillations was determined by a method similar to that described earlier, using fixed attenuators and a calibrated wavemeter. As the measurements showed, the absolute power of the oscillations at frequency 825 Mc was of the order of 150 to 200 W, and at frequency 2400 Mc some 1 to 2 kW. The total power of the oscillations radiated in the 825-Mc frequency range was 3 to 4 kW pulsed. Thus, the total power of the radiation recorded in the beam was 6 to 8 kW pulsed. Estimates show that, for the above beam and plasma parameters, the electric field strength of the excited wave at the end of the beam-plasma interaction region (f - 825 Mc) equals 0.5 to 0.6 kV/cm. The properties of the oscillations excited in the beam in the frequency region -.825 Mc are analogous to those of the oscilla- tions excited in the plasma and described earlier. The time at which the oscillations excited in the beam at frequencies in the range 825 to 1100 and -2400 Mc appeared was also determined. This time was determined from oscillograms similar to those shown in Fig. 2. The time was reckoned from the beginning of the voltage pulse applied to the cathode of the electron gun up to the maximum amplitude of the high frequency envelope. A graph showing the variation of the time of appearance of the oscillations with the air pressure is shown in Fig. 8. As seen from the graph, oscillations at frequency 2400 Mc are excited in the beam earlier than oscillations in the frequency range 825 to 1100 Mc. With increasing pres- sure, the time of excitation of the oscillations in the beam diminishes. 0 104 2 6 810-4 2 It 6 8 10-3 2 p, mm Hg Fig. 8. Time of appearance of oscillations excited in the beam as a function of the pres- sure of the working gas; 1) Frequency 825 Mc; b) frequency -2400 Mc. Discussion of Results Obtained and Comparison with Theoretical Data Let us compare the experimental results obtained with theory. As indicated in [3-6], the interaction of an electron beam with plasma situated in a longitudinal magnetic field may result in the development of instabilities in both beam and plasma owing to the Cherenkov-Vavilov and anomalous Doppler effects. In order to determine the possible spectrum of excited frequencies and gain factors we used a dispersion equation describing the interac- tion of a bounded electron beam with a bounded plasma situated in a longitudinal magnetic field, obtained by M. F. Gorbatenko [611 In the case in which the radius of the plasma column b equals the radius of the beam a, the dispersion equation has the form where 11 1 + + it- 9 N'Y2 [ U.- X- ?(y ?x)" x2 ( 81 n y?x)2J b2coti SY,2, k3v N' = 8= , 0 = (")2 , y = , x = ? A.pv2 wo (00 (00 'All the succeeding calculations were also made by M. F. Gorbatenko. Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 (4) 9 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Li-311 Y Wo -2-u 2 -1 0 1 2 u ?? 4 Fig. 9. Solution of dispersion equations; ?0-0? for plasma, ? for beam with plasma (values of parameters indicated in text). (here k3 = w /v, v = beam velocity, vo = constant lying within the limits 2.4 to 3.8). Dispersion relations for our beam and plasma parameters are given in Fig. 9, where u2 = 8, N' = 2, e = 0.125, which corresponds to N 2 ? 1010 cm-3, H = 1200 G, current = 5 A, v 6 ? 109 cm/sec, and b = 2.5 cm. It follows from an analysis of the dispersion equation that in our case the regions of excitation of the oscilla- tions are situated close to the frequencies coo and WH. In the wo region the excitation may be caused by the Cherenkov-Vavilov effect (region 1) and anomalous Doppler effect (region 3) on the forward wave. In the WH re- gion excitation may be caused by the Cherenkov-Vavilov effect (region 2) and the anomalous Doppler effect (re- gion 4) on the backward wave. The growth increments (imaginary part of x) calculated for these regions in our case have the following values; for region 1 the maximum growth increment equals 0.26, for regicn 2 it is 0.12, for region 3 it is 0.07, and for region 4 it is 0.06. Thus the growth increment has the largest value in region 1, in which the oscillations are excited by the Cherenkov-Vavilov effect on the forward wave (normal dispersion). The frequency fi of these oscillations is close to the electronic Langmuir frequency of the plasma f, but not coincident with this (fi ?0.6 f 0) owing to the bounded state of the beam and plasma and the influence of the strong magnetic field. The experimentally measured frequency of the oscillations most amplified is fi ? 0.53f for a current of 5 A (see Fig. 5a and table), in very good agreement with calculated data. The frequency of the most amplified oscilla- tions was also measured for an 8.5 A current; this equals 0.44 f 0 (see Fig. 5b and table). It follows from the calculations that in our case slow electromagnetic waves with a wavelength of about 7.2 cm in the plasma (frequency 835 Mc) must be excited. The measured wavelength in the plasma is 10 to 12 cm (13 0.3) (see Fig. 3a). Thus slow waves with a phase velocity approximately equal to the beam velocity (vo v) are in fact excited in the present experiments. The calculated value of the gain factors for the conditions ruling in the experiment equals 0.32 cm-1; the ex- perimental values lie between 0.21 and 0.26 cm-1 (see above). Thus the gain factors are also in agreement with calculated values. 10 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 As indicated in [8, 9], on interaction between an electron beam and plasma a convective (transport) instabil- ity should be developed. The graphs of Fig. 3a and b confirm the validity of this principle. During the measurements we failed to observe any marked excitation of frequencies in the neighborhood of the electron cyclotron frequency WH. This corresponds to theory, since in our case the increments in this region of fre- quencies equals 0.12, which is considerably less than that for frequencies in the neighborhood of coo (0.26). Hence the oscillations excited in the frequency range ?800 to 1100 Mc are longitudinal waves in the beam and plasma caused by the Cherenkov-Vavilov effect (forward waves). The experimental results agree with the main theoretical principles. Together with oscillations in the frequency range 800 to 1100Mc, for certain conditions oscillations of fre- quency 2400 Mc are also excited. The frequency of these oscillations is practically independent of the plasma density. In time, they develop earlier than the oscillations in the frequency range 800 to 1100 Mc (see Fig. 8). Their frequency does not depend on the longitudinal magnetic field strength, but depends on the velocity of the electron beam. The intensity of the oscillations in the region of 2400 Mc, in contrast to the plasma oscillations, rises on increasing the longitudinal magnetic field strength. It should be mentioned that the excitation of oscillations at frequencies of 2400 Mc only takes place in the case in which a reflecting metal surface is placed at the end of the region in which beam and plasma interact. All this strongly suggests that such oscillations are caused by the electron beam itself, not depending on the param- eters of the plasma. It may be assumed that in the present case parametric excitation of oscillations takes place. Additional measurements must nevertheless be made in order to decide the question completely. It follows from the measurements made in [1-3] that, for current 5 A, energy 15 keV, and pressure 4 to 6. 10-4 mm Hg, the electron beam loses up to 18% of its initial energy as a result of interaction with the plasma (-13 kW pulsed). This energy goes into excitation of oscillations in the beam and plasma, and also into "heating" these. As measurements show, the power of the excited oscillations is 8 kW pulsed, i.e., at least 60% of the total beam energy loss. We must note, however, that in the experiments carried out the coupling of the beam with the helical junc- tion, and hence the power of the excited oscillations was not optimal, so that it is quite possible that only part of the power in the hf oscillations was taken off from the electron beam during its passage through the helical junc- tion, and hence the power of the excited oscillations may be in fact higher. Some part of the energy goes into heating the beam and plasma. The authors express their thanks to M. F. Gorbatenko for carrying out the calculations. LITERATURE CITED 1. A. K. Berezin, et al., Atomnaya energiya, 14, 249 (1963). 2. A. K. Berezin, et al., Collection; Plasma Physics and Problems of Controlled Thermonuclear Synthesis [in Russian] (Kiev, Izd. AN UkrSSR, 1963), Vol. 3, p. 125. 3. A. K. Berezin, et al., Atomnaya energiya, 18, No. 3 (1965). 4. ya. B. Fainberg, Dissertation [in Russian] (FIAN SSSR, 1960); Atomnaya energiya, 11, 313 (1961). 5. V. D. Shapiro and V. I. Shevchenko, ZhtTF, 42, 1515 (1962); V. D. Shapiro, ZhtTF, 44, 613 (1963); Dis- sertation [in Russian] (0IYaI, Dubna, 1963). 6. M. F. Gorbatenko, Zh. tekhn. fiz., 33, 173, 1070 (1963); Dissertation [in Russian] (Khar'kov Gos. Univ., 1964). 7. A. Engel' and M. Shteenbek, Physics and Technology of the Electrical Discharge in Gases [in Russian] (Moscow- Leningrad, ONTI, 1935), Part 1. 8. P. Sturrock, Phys. Rev., 112, 1488 (1958). 9. Ya. B. Fainberg, V. I. Kurilko, and V. D. Shapiro, Zh. tekhn. fiz., 31, 633 (1961). 11 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 INTERACTION OF PLASMOIDS WITH AN ELECTROMAGNETIC WAVE (UDC 621.384.623) V. I. Veksler, I. R. Gekker, . Ya. Gol'ts, G. A. Delone, B. P. Kononov, 0. V. Kudrevatova, G. S. Luk'yanchikov, M. S. Rabinovich, M. M. Savchenko, K. A. Sarksyan, K. F. Sergeichev, V. A. Silin, and L. E. Tsopp Translated from Atomnaya biergiya, Vol. 18, No. 1, pp. 14-18, January, 1965 Original article submitted April 22. 1964 Some preliminary results of a study of plasma acceleration in circular waveguides are presented. The study was made in the 10-cm range on systems with H01 and Hli type waves and various plasma injectors. Plasmoids with initial particle concentration 1012 cm-3 and above were in- jected with an initial velocity of 5 ? 106 cm/sec from a spark source, or formed directly on the waveguide axis by means of a plasma source with pressure drop, the working vacuum in the accelerator being 10-7 to 10-6 mm Hg. Diagnostics were conducted by means of electric probes uhf methods, and an electrostatic particle energy analyzer. The plasma was contained by auxili- ary external magnetic fields of various configurations. ,Accelerated ions of energy greater than 10 keV were obtained. A proposal was made in 1954 [1] to use the coherent interaction of charged particles with an electromagnetic wave for accelerating plasma. Subsequent theoretical studies showed that, on fulfilment of certain conditions, prolonged acceleration of plasmoids could be achieved, their form and dimensions remaining unchanged [2, 3]. In the present investigation, the initial stage of acceleration and the behavior of the plasma in the wave- guide in the presence of a strong electromagnetic wave were studied experimentally. No attempt was made in these experiments to secure stable acceleration of the plasmoids or reach high energies of accelerated particles. Two systems were prepared; these differed in the type of accelerating wave used and the method of plasma injec- tion. External magnetic fields were used to reduce the diffusion of plasma to the walls of the waveguide. Plasmoids were created in the waveguide at a vacuum of 10-7 to 10-6 mm Hg, both with a spark injector [4] and a pressure- drop plasma source [5]. The investigations thus made enabled us to obtain the first experimental data on the radia- tion acceleration of plasma. Acceleration of Plasma in an H01 Wave The choice of the H01 wave was determined by the following considerations; 1. The field configuration of the H01 wave is such that it permits metal rods of small cross section to be in- troduced along the radius and along the axis without seriously distorting the field pattern. This enables us to place the plasma gun and probes inside the accelerating waveguide. 2. The H01 wave is a wave of the lower type and quite easily excited. 3. The electric field configuration of the H01 wave eliminates the occurrence of uhf discharges to the wave- guide walls. In order to achieve radiation acceleration of plasma in an H01 wave, the system shown schematically in Fig. 1 was constructed. Power from the uhf generator passes into the waveguide tract through a ferrite valve into a King- type wave transformer. Here the H01 type wave in the rectangular waveguide is transformed into an H01 type wave in the circular waveguide. The accelerating waveguide is formed by a stainless steel tube with wall thickness 1 mm. This kind of tube ensures low damping of the uhf waves (around 0.2 dB), and screens the external pulsed magnetic field very little. Between the accelerating waveguide and the transformer lies a vacuum-tight uhf window. The 12 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Fig. 1. H01 wave accelerator (schematic): 1) uhf gen- erator; 2) ferrite valve; 3) wave transformer 4) accelerating waveguide; 5) spark plsma injector; 6) magnetic system; 7) pumping section; 8) uhf absorbing load; 9) screened electric probe; 10) diagnostic window; 11) vacuum-tight uhf window; 12) section with two de- tector heads; 13) high-vacuum line. Fig. 2. Screened electric probe ion current oscillograms: 1) uhf power switched off; 2) uhf power switched on; 3) calibrating signal, period 2 ?sec. ??????\ 3 Fig. 3. Signals from an uhf probe placed at the window of the accelerating waveguide (time base 23 ?sec): 1) Envelope of uhf pulse without the plasma; 2) the same on discharging the plasma gun; 3) signal from gun. 2 Fig. 4. Signals from two detector heads displaced by a distance of 10 cm (time base 17 ?sec): 1) First head; 2) second head. accelerating waveguide is terminated by an absorbing uhf load coinciding with the pumping section. Windows for plasma diagnostics are set in the accelerating wave- guide tube. The 10-cm wave uhf generator operates on a single pulse 8 ?sec long. The mean power flow density through the cross section of the accelerating wave- guide is not greater than 8 ? 103 W/cm2. The standing wave coefficient (voltage) of the whole waveguide sys- tem of the accelerator (without plasma) is not worse than 1.3. A plasma injector of the spark type is introduced along the radius into the center of the waveguide and has a symmetrical "counterweight," reducing distortion of the fundamental wave. The maximum current through the gun in the aperiodic condition is 1.8 kA for a duration of 0.3 ?sec. The spark gun creates a plasma containing around 50% ions of atomic hy- drogen. The total number of ions is 1018 to 1018. In order to reduce the passage of plasma to the walls of the accelerating waveguide, both a longi- tudinal magnetic field and also a magnetic field created by a system of straight conductors parallel to the waveguide axis, fed by currents moving in al- ternate directions, are employed. During the motion of the plasma in the accelerating waveguide (20 ?sec), the magnitude of the magnetic fields remains prac- tically unaltered. The accelerated ions are recorded by means of screened electric probes [6] introduced from the end of the waveguide. Figure 2 shows a superposition of two oscillograms giving the probe ion currents for uhf power "on" and "off"; the longitudinal magnetic field here equalled 300 G, and the gun to probe distance was 45 cm. We see from Fig. 2 that switching on the accelerating field leads to a shift in the velocity dis- tribution of particles in the bunch in the direction of higher velocities. It was not possible to determine the total number of accelerated particles using the probe method, since this can only separate out particles with velocities exceeding the maximum particle velocity in the ab- sence of the uhf accelerating field. Estimates showed that there were not less than 1012 accelerated par- ticles, the maximum velocity of some of the ions ex- ceeding 108 cm/sec. It should be noted that these ex- 1)erimental values agree satisfactorily with calcula- tions in order of magnitude, if we assume that the momentum of the uhf wave is transferred completely to the plasma. Varying the number of particles created by the plasma injector by an order of magnitude (by varying the conditions of the spark gun) does not great- ly affect the number of accelerated particles. 13 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 1 2 Gas Fig. 5. Principle of H11 wave accelerator; 1) uhf generator; 2) ferrite valve; 3) detector head; 4) wave transformer (H ri---03-.14)1); 5) vacuum-tight uhf window; 6) high vacuum pumping line; 7) electric probe; 8) absorbing uhf-load; 9) electron multiplier; 10) electrostatic (127?) par- ticle energy analyzer; 11) magnetic system; 12) accelerating waveguide; 13) pressure-drop plasma source. Fig. 6. Oscillograms (time base 6.4 psec); 1) uhf envelope; 2) signal from electric probe. I, rel. units 1,0 0.6 0,4 0,2 0 2 4 6 8 10 12 14 Energy, keV Fig. 7. Energy spectrum of accelerated ions; 1) uhf wave power P = Pmax ; 2) P CI. 5 Pmax? Operation with the magnetic field created by the system of straight conductors required the introduction of a glass tube inside the metal accelerating waveguide, since, as experiments proved, the plasma losses during motion in a metal tube with this field were large (two orders in density lost in a distance of 1 cm). This is connected with the fact that the polarization of the plasma arising from the presence of the transverse magnetic field is taken by the metal walls, and polarization currents flow constantly in the plasma, retarding its motion. The introduction of a glass tube leads to the insulation of the plasma from the metal walls. The results obtained on accelerating plasma by an uhf wave in this field were analogous to those described earlier. 14 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 By means of an uhf probe (detector head) set at the diagnostic windows of the accelerating waveguide, it was observed (Fig. 3) that, after discharge of the plasma gun, the waveguide practically cut off, that is a "block" was set up for the uhf-wave independently of the external magnetic field. This was confirmed by measurements made with the help of two detector heads [7] in respect of uhf power reflected from the plasma. The presence of the reflected wave appeared in the form of modulation of signals sent from the detector heads (Fig. 4). Analysis of the experimental data showed considerable reflection of the uhf wave from the bunch, reaching some 900/o. Acceleration of Plasma in an H11 Wave The intensity polarized H11 wave is a wave of the lower type in a circular waveguide. A possible failing of the H11 wave is the considerable electric field strength at the waveguide wall, facilitating the development of uhf breakdown. On the waveguide axis, in a condition well removed from critical, the H11 wave takes on a form close to a plane wave, and in certain cases rthis may make for easier comparison between experimental and computed data. Figure 5 shows the principle of a system for accelerating plasma in an H11 wave. The uhf generator is the same as that in the arrangement of Fig. 1. Excitation of the linearly polarized Hn wave is effected by smooth transition from a rectangular waveguide. The accelerating waveguide also consists of a stainless steel tube 1 mm thick. The magnetic field is formed by a system of six straight conductors with currents flowing in alternate di- rections. A pressure-drop source is used to form the plasma. The principle of this source is described in [5]. A plasma pinch with particle concentration above 1012 cm-3 is formed in a discharbe between electrodes A and K (see Fig. 5) in a longitudinal magnetic field. The pressure of the hydrogen admitted at, the anode A and cathode K equals 104 mm Hg. The pressure in the accelerating waveguide is 10-6 mm Hg. The pressure difference is maintained by means of four diffusion pumps. After the plasma pinch has been formed, the A-K voltage is switched off and the magnetic field of the straight conductors is switched on (the time of growth of-the field is 25 ?sec). After the magnetic field has reached a value of the order of 103 G (at the stops) and the plasma is squeezed away from the walls, the uhf field is switched on. A displacement of the standing waves takes place in the waveguide as the plasma moves. This is recorded by a detector head (see Fig. 5). The observed modulation of the uhf oscillation envelope is analogous to that shown in Fig. 4. The accelerated particles were recorded by an electric probe screened from the uhf field (asymmetrical double probe); this was introduced on the side of the uhf-load and could be moved along the waveguide axis. Figure 6 shows an oscillogram of the uhf oscillation envelope and signal from the electric probe situated at a distance of 30 cm from the point of injection. In the absence of the uhf wave, no signals are observed from the probe. From the delay in the arrival of the plasma we may estimate its mean translational velocity, which in or- der of magnitude is 107 cm/sec. In order to measure the energy of the accelerated ions, an electrostatic analyzer rotating the particles through 127? was used together with an electron multiplier as ion detector. The analyzer enables us to measure ion energies within the range 1 to 100 keV; it is placed at the end of the accelerating waveguide, 70 cm from the injector. Figure 7 shows the detector current as a function of the energy of the recorded ions. With increasing uhf power the energy spectrum shifts towards higher energies. On increasing the sensitivity of the system, ions with en- ergies up to 50 keV were recorded. Conclusions As a result of the experiments described, radiation acceleration of plasma was established. Accelerated ions were obtained in both systems, independently of the type of wave in the waveguide and the form of the plasma in- jectors. The energy of the accelerated ions increased on raising the uhf power. The total number of accelerated particles is of the order of 1012, which corresponds to the momentum of the uhf wave transferred to the plasma. The maximum energy of the accelerated particles is up to 50 keV. Regarding the small value of the mean energy obtained in the experiments, there is a number of explanations. 74: First, our plasma injectors give a very large number of particles (up to 1013 to 1016) which, for a given uhf field power, cannot all be markedly accelerated. The background of unaccelerated particles blocks the waveguide. Secondly, almost all the hf power is reflected from the plasma. Thus the acceleration region is small (a few cm), and the length of the accelerating waveguide (30 to 100 cm). Thirdly, a wave of the H11 type evidently cannot 15 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 exert a stabilizing influence on the size and direction of motion of the bunch, as would be expected from theory, owing to the large dimensions of the bunch [2]. The results obtained for the acceleration of plasma in waves of the H01 and H11 types are similar to one another (a wave of the H01 type should not stabilize the dimensions of the bunch in thedirection of acceleration). This is also explained by the action of the background particles in screening the uhf wave. Higher energies can evidently be obtained not only by increasing the uhf power but also by forming com- pact bunches at the beginning of the acceleration. The use of inhomogeneous fields [3] for stabilizing the transverse dimensions of the bunches justified itself completely. We had hardly any plasma losses to the walls of the waveguide on using quadrupole or sextupole,mag- , netic fields. Up to the present, however, we have not succeeded in stabilizing the longitudinal dimensions by means of an inhomogeneous field [3]. This is explained by the low initial velocity of the bunch and its small conductivity. It would appear, that, in the initial stages of acceleration at least, we should use two types of wave simultaneously in the circular waveguide, as was proposed earlier [2]. Thus the preliminary experiments carried out offer a basis for developing a fairly clear program to increase the energy and number of accelerated particles. LITERATURE CITED 1. V. I. Veksler, CERN, Symposium, Geneva (1956), Vol. 1, p. 80; Atomnaya energiya, 2, 427 (1957). 2. M. L. Levin, M. S. Rabinovich, and G. A. Askar'yan, Proc. Internat. Conf., CERN (1959). 3. G. A. Askar'yan, et al., Nucl. Fusion, Suppl., Part 2, 797 (1962). 4. W. Bostic, Phys. Rev., 104, 292 (1956). 5. B. P. Kononov and K. A. Sarksyan, Zh. tekhn. fiz., 31, 1294 (1961). 6. K. D. Sinernikov, et al., Plasma Physics and Problems of Controlled Thermonuclear Synthesis [in Russian] (Kiev, Izd. AN Ukr.SSR, 1962), p. 102. 7. A. Bloch, F. Fisher, and G. Hunt, Proc, IEE, 100, 93 (1953). 16 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 DETERMINING THE PERTURBATIONS OF THE PARAMETERS IN THE MAGNETIC AND ACCELERATING SYSTEMS OF AN ELECTRON SYNCHROTRON ON THE BASIS OF AN ANALYSIS OF INFORMATION REGARDING THE BEAM (UDC 621.384.60) I. P. Karabekov Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 18-22, January, 1965 Original article submitted May 13, 1964 A method of selecting and analyzing information regarding the center of gravity of the beam in the cross section of an accelerator vacuum chamber is proposed; on the basis of this, the deviation of certain fundamental parameters of the electron synchrotron may be obtained. The various components making up an accelerator can only be constructed with finite accuracy. Thus the parameters of the magnetic and accelerating systems also can only hold to a certain accuracy around their cal- culated values. The deviation of these parameters by more than the permissible values leads to a considerable loss of particles. Hence the question of measuring the parameters of the magnetic and accelerating systems from beam data and also stabilizing them (by introducing special corrections) demands special attention in setting up modern accelerators. At the present time, methods of controlling the frequency of the accelerating voltage from information regarding the radial position of the beam [1, 2], methods of measuring and stabilizing the frequency of betatron oscillations [3, 4], and a method of correcting the space harmonics of the magnetic field from beam data [3, 5] have been or are being developed. Also being developed are methods of correcting the magnetic field before injection in order to make the start-up of the first rotation of the beam automatic [6]. Different arrangements of the mutual disposition of the sensory beam elements and controlling devices are considered, as well as the number of controllable sections of orbit required [7]. The distortion of the beam trajectory depends on many factors associated with the deviation of various ac- celerator parameters from their nominal values. For effective correction of the beam trajectory and reduction of particles loss during acceleration, it is important to determine the distortion arising correctly, i.e., to separate out the perturbations corresponding to the various parameters. The possibility of separating the perturbations due to different parameters on the basis of an analysis of beam information and of determining the operations required to correct these perturbations enables us to make the tolerances laid on various components of the accelerator less stringent, so making construction cheaper and more rapid. In the present paper, methods of handling beam information and determining therefrom the perturbation of certain leading parameters of the magnetic and accelerating systems are discussed. Determining the Value of the Field Perturbation of the Magnetic Sections of the Ac- celerator, and Corrective Operations for Compensating These Let us suppose that we know the coordinates of the trajectory of the center of gravity of the beam and its derivative with respect to azimuth at the ends of the magnetic section at points 01 and 02. Let us 'assume that our magnet has no field deviation; then the motion of the particles within the magnet in question must be described by the differential equation (1) 2np where 1 is the length of the period of the magnetic system; p the radius of curvature in the magnets, n(0) the field falloff index. Hence 17 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Col% CV*, r' =C0F' cv*/, (2) where F and F are the Floquet function and its complex conjugate, Co arid Co* are complex conjugate conStants, and 6 Varies by 2it Within the lirnits of the gradient period. Solving the algebrait system of equations (2) with respect to Co at points 61 and 02; we obtain r1F1' ? riFt ?2i c 02_ r2Fr ?2i ' (3) (4) in which C01 is identically equal to CO2 for AH /H0 = 0. The existence of some perturbation, however, destroys this identity; AC =CO2 ? C = r1PI ?2i (5) At the same time, the quantity AC, as we know already [8], is associated with the value Of the perturbation(AH /H0)8 via the expression 02 ( 's 2 I (. AC = ) 2i01 \ (6) Let us measure ri, r2, ri, and 4 and determine AC from (5). Expression (6) connects the quantity AC with the per- turbation (Ai-1/ tio)8 in the unit; this connectian, however, is clearly not unequiVocal. What is important in prac- tice is not the exact measurement of the function (AI-I /1-10)e , but the determination of such equivalent Corrective perturbations as will restore the condition AC = 0, i.e., make the atilt perturbation-free. (The deviation of the beam inside the unit must not be tOo great.) Since AC, generally speaking, is a complex quantity, two parameters are required in order to compensate it. Let us divide the unit into two parts (these may, for exampit, be focusing and defocusing parts). Then expression (6) may alas be written in the form On 02 AC = '\ 2 1 ( AHI r ) F* -F K 92-t ) 2ip Ho /eq. 61 (7) where (Alli /H)N is the required equivalent perturbation created by the first correcting cbil,andK is the ratio be- tween the value8 of the perturbations created by the first and second correcting coils. By due choice of the signs arid Magnitudes of (A1-11/Ho)e,4 and K, the right-hand sides of formulas (5) and (`7) may be made equal. EqUating the arguments of the complex expressions iii (5) and (7), we obtain where on On, r ? 1 IM dO e 1 F* a'0 1 //i1(Z) Re (Z) 01 81 K 02 Z Ini() [ F* aqj [ F* (101 Re (Z) Cid Olt Correspondingly the quantity (Mi. /1-ideci is determined by the expression 18 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 (8) 1 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 i AH, _ r2F:' ?riF;.?riFf' -F-riFt \, Ho )eq ? ?N. 62 , ( ?12 a - Y -= { F* dO+K S F* (10} 271 Q el en (9) where en is the coordinate of the middle of the block. For the simple case in which AH/ Ho = const over the whole length of the block, K = 1 and (Nil/ Ho)eq = AH /Ho. Formulas (8) and (9) do not depend on normalization and the initial phase of the Floquet function. Still more im- portant, the results of determining (6111/H0)eq .and K do not depend on the number of passages of the beam through the magnet over which the quantities r and r' trieasured by the signal electrodes are averaged. This enables us to increase the precision of determining (Alii /H0)eq and K for a given precision of measuring the momentary values of r and r'. In order to determine the sensitivity of the method, calculations were made for, the Erevan annular electron accelerator EKU with the FOFDOD magnetic structure; the radius of curvature of the magnets was -25 m, and the mean ring radius -30 m. Calculations showed that use of the values of r and r' measured at the ends of the block yielded a reliable value of AH /Ho = 0.5% for the perturbation. Thus, the occurrence of a deviation AH /Ho = 0.50/0, constant over the length of the bloc, leads to a change of Ar2 = 1.21 mm and Ari = 1.18 mm/rad in the coordinate at the exit from the block. The coordinates of the beam axis are conveniently measured at the ends of the mag- netic block, since the derivatives with respect to azimuth do not alter in the free spaces. For greater reliability in the determination of the derivatives we may then use information from sensing devices situated at the entrance and exit of two neighboring blocks. The fair distance between these points ensures good precision in determining r'. Thus for the EKU accelerator for the length of the free space between the blocks AO = 0.86 rad and .l-1/ H0 = 0.50/0, Ar'6,0 = 1.00 mm. The value of the derivatives may be made more exact by additional measurements in the free spaces. The method enables us to determine and compensate perturbations in the magnetic system caused by im- precise installation of the blocks and chance displacements. In the case in which the beam-sensing elements are fixed in the vacuum chamber, such perturbations will be sensed as phase shifts of the Floquet function. In order to distinguish the case of loss of beam inside the magnetic block from that of its, emergence with zero parameters, we must use the total signal from the sensing elements situated at the end of the block. If the over-all signal at the block entrance is zero, then a special signal proportional to the maximum possible displacement of the center of gravity of the beam at the block exit must be applied to the analyzing device. From this signal and the measured values of 1.1 and r1 we may determine the values of (AH1 /H0)eci and K. For successive correction of the blocks, however, this case in not very probable. Thus, for the EKU, loss of the beam inside the block in the case of zero entrance parameters (r1 = 0, r1 = 0) means a field deviation of AH /Ho r:-.1 300/0 in the block. Determining the Magnitude of the Frequency Deviation of Betatron Oscillations Suppose we know for certain that in one of the magnetic blocks the magnitude of the perturbation is constant along the whole length. If r1, ri and r2, 1.2 are the measured values of the coordinates and the derivatives of the beam center of gravity at the ends of this magnetic block, then the value of K I determined from Eq. (8) is in- dependent of the value of AH /Ho. Let us suppose that the frequency of the betatron oscillations varies by Av. Here we obtain a new distribution of the Floquet function with respect to azimuth, leading to an additional dis- tortion of the beam trajectory. If the values of r1 and r1 remain the same as for Ay = 0, the r2 and ti acquire cer- tain increments Ar2 and A*. If the new values of r2 and r2 are put into expression (8), the value of K will differ from unity. Let us examine the variation with Ay of the phase of the expression Z d-r;Fi*. (10) The new values of the Floquet function and its derivative for a change of Lv in the frequency of the betatron os- cillations will to a first approximation be given by the expressions; FejAv? F iAv0F, cp'= iAvFeim'e?F'eim'o F' (1+ iAv0)? iAvF. 19 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 From the measured values of r1 and ri and the new values of the Floquet function, we determine 72 and Making some fairly simple mathematical transformations and rejecting terms containing Av2, we obtain r2= rzfr ?r2fo ?Re {K (riri)+ L (rolAv, (12) r2= r2 fr -1- r2fo ? Re {N (riri)+ E (n)} Av , (13) where r2fr and r2fr, are the values of the coordinate of the center of gravity and its derivative at the point 02 for Av =0 and AH /Ho = 0 ("free"), and r2f0 and r2f0, are the "forced" solutions at the point 02 for A v = 0; K(ri, ri) and N(ri, ri) are functions determined by the values of r1 and ri, and L(II) and E(II) are functions determined by the magnitude of the perturbation in the magnetic block. These functions have the form K (r jr) = ? F2 (02? 01) ? riF:F 21, (14) N (riri) = ?4-{[F2--F; (02 ? 01)] [riF ,*" ? ? riF:F;} , (15) Here 82 92 L (n) = F211F*0 d0 ?02 c F* do) , Eli ea 02 E (n) = -2- [ (02F; ?F2) F* de ? F; F*0 de . ?2t,IQ Ho ? Substituting expressions (12) and (13) into (10), rejecting terms containing Av2, and allowing for the condition r2fr we obtain (16) (17) (18) (19) Z = r2 fo F:' ? AvRe {K (r iri)+ L (n)} F:' ?r fo F; + Av Re {IV (r tri)+ E (II)} F. (20) Let us form the ratio of the imaginary part of expression (20) to the real part. Expanding the resultant expression in series in powers of Av and restricting consideration to the term containing the first power of Av, we obtain; Im ()= r2 fo ImFt" ?ri fo InzF Av T Re (Z) r2f0 Rely ?r'2 fo ReF An expression for T(ririlI) can easily be found by making some simple transformations. In expression (21) 02 m F* T2 fo I mFt' ?r 2' fo IniF; r2 fo ReFr fo 02 ,"="- const. Re F* del (21) (22) 'Thus for the t KU accelerator this quantity equals ?0.195. The function T(ri,rilI) is a linear function of the prin- cipal coordinates, their derivatives, and the value of the perturbation. Hence the value averaged over many rota- tions of the beam is = T(II), i.e., constitutes a certain constant determined by the magnitude of the perturbation specially introduced into the block. Substituting expression (21) into (8) and averaging over many turns, we determine the deviation of K from unity. This deviation may, for a known value of T(11), be graduated in units of A v. 20 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Determining the Amplitude and Frequency Coherent Phase Oscillations In electron synchrotrons calculated for a particle energy of the order of several BeV, owing to the high mul- tiplicity of the high frequency voltage the direct method of measuring the. amplitude and frequency of coherent phase oscillations presents considerable difficulty. This information may be obtained far more simply by measuring the coherent deviation of the particle energy from the equilibrium value arising for phase oscillations. It is known [8] that for Es >> mc2 da) AE dt E, (1), (23) where q is the multiplicity of the frequency of the accelerating voltage, c4)0 is the frequency of rotation of the par- ticles in the orbit, and cc is the logarithmic derivative of the length of the orbit with respect to momentum. Let us suppose that the momentary phase of the accelerating voltage for which the center of gravity of the beam passes through the middle of the accelerating gap under the influence of an external perturbation, is given by the expression 0= Os+ Acp cos Qt, (24) where cps is the value of the equilibrium phase, Acp the amplitude of the phase oscillation, and Q the oscillation frequency. On this assumption, expression (23) will have the form AE ? ? AyS2 sin SU Tooa (25) Let us measure the coordinates and their derivatives with respect to azimuth at the ends of one magnetic block and determine the value of the total perturbation arising as a result of deviations of the magnetic field values in the block and the mean particle energy in the beam: 02 * All 7.21T ?r;F* ?riF,*: +r/?* =(?-)2 1 ?c (0) dO ( 1 2 AE F* dO. 23t 'CO . Ho 23t j es el. 0.1 The deviation of the magnetic field in the blocks during the acceleration cycle changes considerably more slowly than the coherent deviation of the particle energy, owing to the phase oscillations. Hence the duly separated high frequency component of expression (24) determines the amplitude and frequency of the oscillation of (AE/E)t. The frequency of the variation of (AE/E0)t may be measured by many well known methods. The amplitude of the phase oscillations will be connected with the maximum deviation of the beam energy by the expression qcooa AE = Q max. Ao In conclusion, the author considers it his pleasant duty to thank S. M. Rubchinskii, t. L. Burshtein, A. A. Vasil'ev, Yu. F. Orlov, S. A. Kheifets, and V. M. Kharitonov for discussion of the work, and M. A. Garzoyan for help in carrying it out. LITERATURE CITED 1. A. A. Vasil'ev, A. A. Kuz'min, and Yu. S. Ivanov, Pribory i tekhnika eksperimenta, No. 4, 111 (1962). 2. 'k . L. Burshtein, Yu. S. Ivanov, and A. A. Kuz'min, Pribory i tekhnika eksperimenta, No. 4, 102 (1962). 3. t. L. Burshtein, et al., Atomnaya energiya, 12, No. 2 (1962). 4. A. A. Vasil'ev, A. A. Kuz'min, and V. A. Uvarov, Pribory i tekhnika eksperimenta, No. 4, 134 (1962). 5. A. A. Vasil'ev, Transactions of the International Conference on Accelerators, Dubna, 1963 [in Russian] (Moscow, Atomizdat, 1964), p. 871. 6. A. A. Vasil'ev, Dokl. AN SSSR, 148, 577 (1963). 7. A. I. Dzergachard, V. A. Karpov, See [5], p. 867. 8. S. A. Khaifets, The Electron Synchrotron [in Russian] (Erevan, Izd. AN Arm. SSR, 1963). Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 (26) 21 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 PHASE STABILITY OF A SYSTEM OF PARTICLES IN SELF-REGULATED ACCELERATORS (UDC 621.384.60) f. A. Zhil'kov and A. N. Lebedev Translated from Atomnaya tnergiya, Vol. 18, No. 1, pp. 22-28, January, 1965 Original article submitted February 13, 1964 The problem of the phase stability of a sytem of particles in a cyclic accelerator with frequency autocorrection of the accelerating field along the beam is considered. The stability of nonlinear synchrotron oscillations is investigated by the kinetic-equation method for arbitrary characteristics of the self-regulating system. Some general stability criteria are obtained. One of the most promising routes in the development of cyclic accelerators is the wide application of beam autocontrol systems, which peimit a considerable reduction to be made in the allowance for different parameters of the magnet and high-frequency systems [1, 2]. In practice the most important and at present most widely used system of this type is the system of automatic frequency trimming of the accelerating field by radial (or phase) shift of the accelerated particles, thus permitting considerable increase in the stability of coherent synchrotron os- cillations relative to various perturbations. An investigation, which was carried out in linear approximation [3], showed that the frequency of coherent synchrotron oscillations, i.e., oscillations of the beam center, is increased by a factor of (1 + K)112 where K is the amplification factor in the self-regulation cycle with respect to the radial shift of the beam. In addition, it was proved possible to achieve strong damping of the coherent oscillations. The oscillations around the beam center, discussed in [3] in linear approximation, are not subject to the effect of self-adjustment and have the normal frequency of synchrotron oscillations. Naturally, the nonlinear nature of the phase motion should alter this picture, leading to coupling of the beam-center oscillations and the free oscilla- tions around it. The representation of the nature of this coupling can be obtained in this case if the first moments of the distribution are considered as in [4] and, more rigorously, in the first and second sections of this paper. However, for the subsequent statement of the problem the discussion should center on the stability of a system having a large, almost infinite number of degrees of freedom corresponding to a large number of accelerated par- ticles. It is shown below that this investigation can be carried out by the kinetic-equation method with self- consistent interaction. It was found possible to obtain a very general characteristic equation, useful for a wide class of regulatory systems and giving in principle an answer to the problem .as posed. The paper has been limited to the investigation of certain general properties of this equation which are in- dependent of the choice of the feedback circuit parameters. The calculation of an actual system, requiring com- putation of a number of purely radiotechnological considerations, will be carried out in a separate report. 1. Statement of the Problem In order to describe the synchrotron oscillations of an individual particle, we shall use the canonical conjugate variables 23520 dE eV sin cps (E) ' Es where w(E)is the rotational frequency of a particle with energy E; V and ca are the amplitude and phase, respec- tively, of the accelerating potential; 00 is the frequency of the linear synchrotron oscillations in the absence of self-regulation. We shall neglect the explicit time dependence of the parameters 00, V, etc. 22 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 If we choose T = Qt as an independent variable, then the synchrotron oscillation equation can be written in the form uF (q)) ?-= cos q)? cos cps Sin cp.= u+ A (T), (1.2) where A(T) is the correcting adjustment to the frequency of the accelerating field, associated with the energy of the self-regulation system and expressed in units of Q. According to Eq. (1.2) the kinetic equation which describes the change of the distribution function has the form f(u, T) [-I-q, T) = O. (1.3) The quantity A can be related to the various characteristics of the beam as a function of the following (or tracking) system parameters: to the coordinates of its center of gravity, dimensions, etc. In the general case, if we denote by K(u, , w) the reaction of the system with frequency u.) to a shift of a point-like beam with respect to u, and yo, then in Fourier presentation the quantity A(w) can be described in the form of a linear function of the dis- tribution function A (6)) =K (u, cp, o) f (a, cp, co) du dcp. (1.4) Equations (1.3) and (1.4) represent a closed integral-differential system describing the behavior of the entire buildup of particles in an accelerator with self-regulation. We shall proceed further from the fact that this system permits a stationary (equilibrium) solution: f o(u, co), Ao = const. We note that the constant component of the cor- rection signal Ao, associated with the equilibrium distribution, can be set equal to zero without limiting the gen- erality, since this leads simply to an insignificant indeterminacy of the equilibrium energy. It is not difficult to see that an extremely wide range of equilibrium distributions is possible, on which a unique requirement is im- posed ? the constancy of f 0 along the phase trajectory [see Eq. (3.2)]. In the general case, the center of the equi- librium phase cos as a consequence of the nonlinearity of the synchrotron oscillations, i.e., in conseque ce of the deviation of the potential well ? f Fcico from parabolic. This shift is found to be greater, the larger the equilibrium phase of the beam. In accordance with the assumptions made above concerning the absence in the equilibrium state of a correction signal Ao, it should be assumed that the center of the equilibrium distribution with respect to energy coincides with Es, i.e.; uo = 0. 2. Method of Moments First of all we shall consider the case of a beam of quite small dimensions, since by this example it is easier to trace the physical significance of the relationships obtained below. Multiplying the kinetic equation (1.3) by ( ? 0)Juk and integrating with respect to the entire phase space, we obtain a system of equations for the set of distribution moments Mi,k = (99 ? vo)juk (the bar denotes the average): II; j-1, /AM5,h kF (9)) =0. (2.1) In accordance with the assumption made above concerning the smallness of the beam, the distribution moments should be reduced rapidly with increase of their number. Having expanded the quantity F(yo) in a Taylor series around the point cpo, Eq. (2.1) leads conveniently to the form j, Pi 41? /A (t) M_1, co F (1) (W m? 1! j+17 h-1-= 0. 1=0 (2.2) We note that the system of equations (2.2) is, generally speaking, nonlinear since the correction signal A is con- nected with the moments by the relationship (1.4). We now use the fact that the phase dimensions of the beam 23 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 A > 1) the "rigidity" of the oscillations of the center of gravity of the beam can be increased strongly and even damping of them can be originated (for example, for Re K2 < 0). The natural oscillation frequency of the second moments, which in this approximation remain free, is equal to twice the frequency of the synchrotron oscillations in the ab- sence of self-regulation. By taking into account terms ?A4, the corrections to the natural frequency of the first two moments can be found (terms ?A3 make it possible to determine only the partial frequencies a the third moments, equal of course to ?1 and ?3). These calculations are quite cumbersome since they are associated with the calculation of a deter- minant of the 10th order. We shall mention only the most important result ? the correction of the first approxima- tion to the real partial frequency w = 2 of the second moments, since small corrections to the complex frequency (2.6) cannot change the nature of the motion qualitatively. If we assume that 0./ = 2 + a(1 a I 0) and buildup (Im a < 0) of the oscillations are also possible at this frequency. The con- dition of stability for I K1,2 >> 1 is, obviously, the inequality Kt Re ?K2 1 Po?(Atillo)Pie(131-p?)1+ (A:1/140) Pie(Pi?P?)t (13) In order to select the cycle for the reactivity measurement with respect to the positive period, calculations were carried out of values of 6 for various reactivities and holding times of the system at a constant power level. The results of the calculations are given in Figs. 1 and 2. The deviation of the measured period from asymptotic was determined experimentally for a reactivity of 0.10/0 in the critical beryllium oxide assembly being investigated. As the asymptotic period in these measurements, the steady state period was used as measured over 200 sec after a positive reactivity surge (for the stated time the value of 6, as shown by the calculations, does not exceed 0.1% The results of?the experiment and of the calculation are compared in Fig. 3. Measurement of Reactivity The calculated and experimental data obtained concerning the magnitude of the deviation of the measured period from asymptotic permits selection of the reactivity measurement cycle to be made. 38 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 In the experiments on the critical assembly, .the results of which are given in the present paper, the holding time at constant power was taken equal to 1000 sec, A shorter holding time, although advantageous from the point of view of reducing 6, is inadvisable since prior to measurement of the positive period it is necessary to determine the critical state of the assembly. The positive period was measured over 100,200 sec after the reactivity range. The measurement results were analyzed by the "inhour" equation.. The cycle taken for the measurement of the period and the equipment used permitted the reactivity to be measured to an error not exceeding ?1.55. A section of the control rods was,cali- .. brated by the method mentioned. The data obtained are shown in Fig, 4, In addition to the positive reactivities, the negative reactivities were measured by an integral method [7]. It has been shown that after a Laplace transform of Eq. (1) and (2) expression (5) is obtained which, for Ulieff ,lXicio/keffRiTio and for p = 0 is transformed to the form where whence N (0) ? -Q keff ' co N (0) = n (t) dt, to no xi lkicio yipi N (0) k p?n 2,,? ? eff o If the reactor has previously operated for sufficiently long at stationary power, so that the condition ci0Xil/k5in0 = 1 is fulfilled for all the delayed neutron groups, then -p is determined from the expression __ Pole \-1 ViN N(0)1t-eii ? In contrast from the "inhourn equation, Eq, (14) depends to a great extent on the quantity )Li. The presence of photorieurrons in the system increases the significance of this fact. However, the difficulties can be overcome, to a considerable extent, if the photoneutron precursors do not reach saturation and if the quantity cio be deter- mined, i.e? ci at the instant of time to, when the source is eliminated from the solution of the equation (14) (14') dci (t) kpi = n (t) Xici (t). dt 1 For the case of constant value of k, the solution has the form to c0 (t)----= (143i11) e?mo n (t) ekt Cit. 0 After substituting Eq. (16) in Eq. (14), the equation for the reactivity being measured assumes the form Q to k E viPiS n(t) (t-td dl (15) (16) (17) The effect of long lived groups of photoneutrons on the magnitude of the reactivity can almost be eliminated by change of to, and the time of the experiment can be shortened significantly. 39 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 1.0 0 4 RIO Delay time, sec a ZOO Relative contribution reactivity 40 Delay time, sec Fig. 5. Undersaturation factor of delayed photoneutron groups (a) and relative contribution to reactivity injected by the photoneutron groups (b) as a function of holding time of the neutron source in a sub- critical assembly for p = -1 ? 10-3 (figures near the curves are the number of groups). Preliminary results of the use of the integral method in the absence of saturation of photoneutron precursors were obtained on an analog computer, by which the behavior of the reactor was reproduced with six delayed neutron groups and four photoneutron groups. By choice of scale factors it was presumed that as a result of saturation of the precursors Un = Ucn and then nokikeff E (v/k)cucio/und ?Q = n (t) dt to (17') The complex yi i3i/ Xi of each group of delayed neutrons enters into this formula with its undersaturation factor .13 = U /Uno, depending on the holding time of the source in the recator. cio The data obtained (Fig. 5a and 5b) permit conclusions to be drawn concerning the magnitudes of the under- saturation factor and the relative contribution to the reactivity of each photoneutron group as a function of the holding time of the source in the reactor, and to distinguish two possibilities for measuring the reactivity by the integral method; by using the undersaturation factor [Eq. (17)] or by choosing an operating cycle which will permit the photoneutron groups to be neglected, within the limits of the stated accuracy. We have used the first possibility. A section of the control rod was calibrated by the integral method discussed in the critical assembly being in- vestigated. The data obtained are given in Fig. 4, where both methods of reactivity measurement are compared. Determination of the Efficiency of the Delayed Neutrons and Photoneutrons In describing transient processes or in measuring reactivity, it is necessary to know the value of y or 13eff = y i3, the effective fraction of delayed neutrons and photoneutrons. For the critical assembly being investigated the quantity y was measured by the method of substitution of fuel by absorber and by determining the change of reactivity originated by this substitution. It follows from perturbation theory [8, 9] that the change of reactivity of the system as a result of changing t he quantity of fuel and absorber is expressed as L\keffI beff = r dv VC?r fa ,u +1 fTvTF TFT ) dv vF*SE/Fdu-F6X/FITT] ? SE,F*F du? 6M,T.F TF1,1 , n, - tir u where F is the neutron flux; F* is an adjoint function; the suffix T denotes a thermal group. 40 (18) Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 If the absorption in the absorber inserted is equal to the absorption in the spent fuel at the site of location of the absorber, then (SEc = 0 and there will be no terms in expression (18) which will take absorption into account. If the system is divided into N parts and this substitution of fuel by absorber is performed each part, then we can write 1=1 S vF*SEfF du dv + . . . + S S vF*15EfF du dv eff ti v N u k eff 5 5 vF*.EfF du dv V u =1 (19) (for brevity of writing, the thermal group is included here under the integral sign). The value of Akeff/keff = ky/keffEi Sipo/(po + Xi) is determined in the experiments to measure the reac- tivity, if y is assumed constant for all groups of delayed neutrons and photoneutrons. Then, = N PiPo 1=1 keff )t 1 (20) If the absorption cross section in the absorber inserted, Ea, does not match the absorption cross section of the fueld, c E c, than we obtain from Eq. (18) Akeff S (I! F*F du dv keff vu vF*F.Ef du dv U (21) and the correction to unity should be obtained by a numerical method. Expression (21) was used to determine the experimental error associated with the inaccuracy in Ete and 4, since rs (Av+ Ica) F*F du dv i) U L Ay = 1 k ek ff 1=-1. iPo 5 S vF*F2/ du dv v Po+ki u (22) Experiments to determine y were carried out on a critical beryllium oxide assembly with fuel elements on a Teflon-4 base [10]. Specially prepared plates of Teflon-4 with a boron filler were used as the absorbing elements. The ex- periment consisted in the substitution of the fuel elements by absorbers and in measuring the change of reactivity caused by this substitution. Usually, one or two fuel elements were substituted during a single procedure. As a re- sult of this, in order to reduce the number of substitution and measurements the symmetry used for the assembly was represented by a rectangular prism with square cross section. In order to ascertain the inaccuracy in the quan- tity y, the same experiments were carried out for two cases when the absorption in the absorbing element was not equal to the absorption in the fuel element. It was found for the assembly being studied that y = 1.15; the experi- mental error was 4%. For comparison with the experiment, calculations were undertaken of the quantity y. The difference between the prompt and delayed neutrons, as is well known, is due to their different initial energy. The energy.of the de- layed neutrons from U235 fission products was used in accordance with [2, 11, 12]. The energy of the photoneutrons was determined from consideration of the Be9(y, n)Be8 reaction and from analysis of the y -radiation which forms the photoneutrons (see table). In the formation of a Be% nucleus in the ground state as a result of the reaction Be9(y, n)Be8, the energy of the photoneutron of the i-th group was determined from the expression E.= ?1_67 MeV, where EY i is the average energy of the gamma radiation for the i-th photoneutron group; 1.67 MeV is the Be9(y,n)Be8 reaction threshold. The average energy with which photoneutrons are formed was found to be 670 keV. 41 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Since in the experiment carried out the value of y obtained was relative to all groups of delayed neutrons and photoneutrons, the calculated value of y was determined in accordance with [13] from the expression y eX2 AT, where AT is the difference between the increase of fission neutrons and delayed neutrons, and photoneutrons, whose spectrum was obtained. For the assembly investigated x2A T = 0.13 and y = 1.14, Thus, consideration of the effect of the Be9(y, n)Be8 and Be9(n, 2n)Be8 reactions on the kinetics of systems con- taining beryllium has shown the necessity for taking these reactions into account for reactivity measurements and in choosing the measurement cycle. Analysis of the data obtained shows that by choosing the holding time of the reactor at a constant power level to be about 1000 sec, by measuring the period over 200 sec after the reactivity surge, the deviation from asymptotic of the period to be measured does not exceed 0.1% for periods within the inter- val 10-100 sec. In measuring negative reactivities the taking into account, as suggested, of changing the amount of photo- neutron precursors by the undersaturation factor makes it possible to use the integral method, which gives ? as com- parison with positive period measurements show ? a satisfactory result. Determination of the efficiency of the delayed neutrons and photoneutrons by the method of substituting fuel by absorber and comparison of the experimental and calculated data showed the good applicability of the numerical formula for the case with photoneutrons. The special significance of the experiment in substituting fuel by absorber consists in the fact that it makes - possible the determination of the constant in front of the summation in the expression yk/keff 6ipo Apo + Xi) in the "inhour" formula, by eliminating the effect of existing inaccuracies in the delayed neutrons, lphotoneutron and neutron parameters of the Be9(n, 2n)Be8 on the result of the reactivity measurement. It is desirable always to carry out a normalized experiment of this nature prior to the precision experiments associated with reactivity measure- ments. In conclusion, the authors tender sincere thanks to N. N. Ponornarev-Stepnoi for valuable advice and interest in the work, and to Ya. V. Shevlov for discussion of the results. LITERATURE CITED 1. G. Keepin, et al., Phys. Rev., 101, 1044 (1957). 2. G. Keepin, et al., J. Nucl. Energy, 6, No. 1/2 (1957). 3. G. Keepin, Nucleonics, 2, 151 (1962). 4. A. K. Krasin, et al., In the Book: Proceedings of the Second International Conference on the Peaceful Uses of Atomic Energy, Geneva, 1958 [in Russian] (Moscow, Atomizdat, 1959), Vol. 2, p. 39. 5. P. Benoist, et al., Proceeding of the 2nd Intern. Conf, United Nations, Geneva (1958), Vol. 12, p. 89. 6. S. Bernstein, et al., J. Appl. Phys., 27, 18 (1956). 7. U. Hojgn, Nucl. Sci. and Engng., 8, 518 (1960). 8. L. N. Usachev, Reactor Construction and the Theory of Reactors. Report of the Soviet Delegation at the International Conference on the Peaceful Uses of Atomic Energy [in Russian] (Moscow, Izd-vo AN SSSR, 1955), p. 251. 9. G. I. Marchuk, Numerical Methods of Assessing Nuclear Reactors [in Russian] (Moscow, Atomizdat, 1958), p. 205. 10. N. N. Ponomarev-Stepnoi, S. S. Lomakin, and Yu. G. Degal'tsev, Atomnaya energiya, 15, 259 (1963). 11. R. Batchelor and H. McHuder, J. Nucl. Energy, 3, 7 (1956). 12. D. Hughes, et al., Phys. Rev., 73, 111 (1948). 13. E. Cross and J. Marable, Nucl. Sci. and Engng., 7, 281 (1960). 42 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 THE CRYSTAL HYDRATE UF4 ? 4/3 H20 (UDC 546.791.4) Yu. V. Gagarinskii, E. I. Khanaev, N. P. Galkin, L. A. Anan'eva, and S. P. Gabuda Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 40-45, January, 1965 Original article submitted July 24, 1964 X-ray diffraction, refractometry, infrared spectroscopy, nuclear magnetic resonance and thermo- graphy are used to study a new hydrated form of uranium tetrafluoride, UF4 ? 4/31120, and its de- hydration products. This is a hitherto unknown monoclinic crystal hydrate of UF4. Its water of crystallization is attached by a hydrogen bond to the fluorine. The water molecules can be divided into three groups of different bond strength, corresponding to three peaks in the absorption bands of the valence and deformation vibrations of the 0 ?H bonds. The substance dehydrates in two stages. Its crystalline class is retained down to the composition UF4 ? H20. On further dehydra- tion down to 0.5H20 it undergoes a phase transformation accompanied by change of structure. The crystalline lattice thus formed is very close to that of the cubic hydrate. [1] and [2] describe two crystal hydrates of uranium tetrafluoride: UF4 ? 2.5H20, which is orthorhombic, and UF4 ? nH20 (0.5 < n < 2), which is cubic. We have found another, previously unknown, hydrate, whose properties are described below. It is a bright grass-green, contrasting with the deep green cubic and pale green orthorhombic hydrates. Its water content, after washing with alcohol and ether and drying in air, corresponds approximately to the composition UF4 ? 1.4H20. By the calorimetric method [3] it was established that the formula is UF4 ? 4/3H20. We investigated this hydrate by x-ray diffraction, optical methods (refractometry and infrared spectroscopy), thermography and nuclear magnetic resonance (NMR). Under the microscope it appears as long thin greenish crys- tals with spherical concretions. It is homogeneous: the crystals are optically anisotropic and have normal and oblique extinctions, showing that they belong to the monoclinic system. Table 1 gives the refractive indices n of a specimen of composition UF4 ? 1.4H20 (measured by immersion, error of order ?0.005), and the density(deter- mined by pycnometry in toluene with error ?0.01 g/cm3). For comparison, the same data, taken from [4], are given for the cubic and orthorhombic hydrates and the monoclinic anhydrous form of UF4. Table 1 shows that the new crystal hydrate differs markedly in crystallographic properties from the cubic and orthorhombic hydrates. The refractive indices are close to those of anhydrous UF4 (only nmin differs appreciably). These optical data show clearly that this actually is a new crystal hydrate of uranium tetrafluoride. X-ray diffraction photographs of the substance were taken with a URS-50I apparatus in filtered Cu radiation. For comparison, photographs were also taken for the stable monoclinic form of anhydrous UF4, the higher, ortho- rhombic, crystal hydrate UF4 ? 2.5H20 and the lower, cubic, hydrate UF4 ? 1.5 H20. They show that the substance under investigation does not resemble any of the previously known forms of UF4 and is not a mixture of them. This confirms the conclusion, drawn from the optical data, that it is a new substance. The x-ray diffraction spectra are shown in Fig. la. To determine the nature of the water bonds in the new hydrate, we investigated its infrared absorption spectra and nuclear proton magnetic resonance spectra. The method of measuring the former spectrum is given in [4], that for the latter in [5]. The infrared absorption spectrum of monoclinic UF4 ? 1.4H20 is given in Fig. 2 (curve 1), the positions of the absorption maxima in Table 2. The absorption band corresponding to valence vibration of the 0 ?H bond has three maxima at 2950, 3365 and 3480 cm-1, that corresponding to deformation vibration of the 0?H bond has less clearly visible maxima at 1565, 1625, and 1645 cm-1. The presence of absorption in this region shows that the compound is of the crystal hydrate type, the water being present as molecules. The absorption band corre- 43 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 44 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 ,130 020 ,110 ,100 ,90 ,80 ,I0 ,60 ,50 ,130 ,120 ,t10 ,I00 190 ,40 ,_30 ,20 , to 180 170 ,60 _,50 ,40 ,30 20 130 100 ,90 ,86 ,70 ,60 28 Fig. 1. X-ray diffraction spectra of UF4 ? 4/31-120, and the products of its dehydration and subsequent rehydration. a) Original crystal hydrate; b) UF4. H20; 8 0F4 ? 0.70H20; d) UF4 ? 0.45H20; e) UF4 ? 0.20 H20; f) UF4 ? 1.40H20, obtained by hydrating specimen of composition UF4 ? 0.20H20. 0 ,60 140 ,30 ,j0 00 0 (D (I) (I) (D -To -To (D -n (D (T) (D o.) 0 0 -0 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 1. Refraction Constants of the Various Forms of Uranium Tetrafluoride Composition Crystal class nmax nmed nmin Ti d425,g/cm3 f2D,cm3 UF4.1.4H20 Monoclinic 1.596 1.586 1.567 1.582 5.79 19.77 UF4.1.5H20 Cubic ? 1.523 ? 1.523 5.08 20.51 UF4.2.5H20 Orthorhombic. ? . 1.545 1.537 1.529 1.537 4.70 23.86 UF4 Monoclinic 1.594 1.584 1.549 1.576 6.68 15.58 TABLE 2. Maxima of Absorption Bands (in cm-1) in Infrared Spectra of Specimens Obtained by Dehydrating Crystal Hydrates of Monoclinic Uranium Tetrafluoride Curve 14o. in Fig. 2 Composition Dehydration temperature,?K Deformation vibrations of 0?H bond, cm-1 Valence vibrations of 0?H bond,cm-1 1 2 3 4 5 UF4 ? 1.401126 U F4 ? 1 00H20 U F4 ? 0 . 301120 U F4.0.101120 ' UFO . 401120 150 200 250 Obtained by hy- drating I UF4 ? 0.101-120 1565; 1625; 1645 1590; 1645; 1675 1615; 1640; 1675 1650 1680 2950; 3365; 3480 2935; 3365; 3480 3465 3470 3400 sponding to valence vibration of the water molecule is markedly broadened and displaced towards longer wave- lengths. It follows that in the new hydrate the water is also attached by a hydrogen bond. The presence of three maxima in the absorption bands corresponding to valence and deformation vibration of the 0?H bond leads to the supposition that in UF4 ? 1.4H20 the water molecules can be subdivided into three groups according to bond strength. Figure 3 shows the NMR spectrum of the crystal hydrate, taken at 290?K, and the derivative spectra at 90 and 290?K. The second moment of the NMR lines is 27.5 ? 0.3 0e2 at 90?K and 23.6 ? 0.9 0e2 at 290?K. At 290? the relative intensity of the central peak, which determines the fraction of mobile water molecules, is about 50/0. The shape at 90? of the new hydrate's NMR spectrum can be explained by the presence of a hydrogen bond be- tween the water molecules and a fluorine atom, OH?F. In this case we should have a three-spin system whose spec- trum should be a triplet (or unresolved triplet degenerating into a singlet). At 90?K the second moment of the NMR spectrum is less for the monoclinic hydrate than for the other hy- drates [6]. This means that, provided the distance between the protons in the water molecules remains constant, the mean H?F distance in this hydrate is relatively greater tharrin the other hydrates of UF4. Figure 4 shows the thermogram for dehydration of the new-hydrate in a vacuum, recorded by means of an PFK-59 Kurnakov pyrometer. The graph shows that the dehydration takes place in two stages. As in the cases of the other two hydrate, an unstable form of anhydrous UF4 is produced and on further heating this undergoes an ir- reversible transformation to the stable form. This process is accompanied by heat evolution, as shown by the exo- thermic peak at 330?K. X-ray and optical investigations were made of the products formed by dehydrating UF4 ? 1.4H20 in a vacuum at various temperatures. The results are given in Tables 2 and 3. As seen from Table 3, during vacuum dehydra- tion the crystal class of the original substance remains unchanged at least until the composition UF4 ? H20 is reached; only its refractive index changes slightly. The x-ray diffraction spectra show that the structure remains unchanged. On reaching the composition UF4 ? 0.45H20 and after further dehydration, the optical properties and x-ray diffraction patterns change appreciably. If a specimen of composition UF4 ? 0.20H20 is washed with alcohol and ether, dried in air for 24 h and then rehydrated, the resulting substance has a very slight anisotropy; all the strong and medium lines in its x-ray diffraction spectrum correspond to those for the cubic hydrate. Additional, very weak lines betray the presence of a small quantity of another phase (which apparently causes the slight anisotropy). Figure 2 gives the infrared absorption spectra of similar specimens obtained by dehydrating the monoclinic hydrate in various degrees. The positions of the absorption maxima are given in Table 2. This table shows that, 46 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 3. Optical Properties of Substances Obtained by Dehyrating Monoclinic Crystal Hydrate of UF4 Composition Dehydration temperature, ?K N Remarks UF4?1.40H20 1.582 Anisotropic spherulites UF4?1.00H20 150 1.570?0.008 The same UF4. O. 45H20 175 1.52?0.01 Almost isotropic UF4?0.20H20 250 1.504?0.006 Weakly anisotropic UF4?1?10H20 Obtained by hydrating 1.545?0.003 Weakly anisotropic UF4?0.20H20 e-60 0-40 .g.2? 4 400 600 1400 1800 2800 3200 3600 4000 Fig. 2. Infrared spectra of monoclinic lower crys- tal hydrate of uranium tetrafluoride and of speci- mens obtained by dehydrating it. 1-5) Specimen numbers in Table 2. 210 -14-12-10-8-6-4-2 0 2 4 6 8 10 12 14 Fig. 3. a) Derivative NMR spectra of monoclinic crystal hydrate of uranium 184 tetrafluoride of composition UF4?1.40H20 at 90 and 290?K; b) [integrated] spec- trum at 290?K. ) spectrum at 290?K; ----) broad component; - .-.-.-.) narrow component corresponding to highly mobile water. Vertical scale chosen so that total area is equal to unity. in the frequency region of valence deformation vibrations of the 0?H bond, the absorption bands of the specimen of composition UF4 ? 1.00H20 retain the structures of the corresponding bands for the original crystal hydrate. On transition to composition UF4 ? 0.30H20 the spectrum of valence vibrations of 0?H becomes simpler. The maxima at the long-wave end disappear, the spectrum being left with one band at 3470 very similar to the absorption band for valence vibrations of 0?H for the cubic hydrate of identical composition [4]. Fig. 4. Thermogram for dehydration of monoclinic lower hydrate UF4 ? 1.40H20. Heating rate 4 deg/min. 47 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 UF4 1.401120 obtained by hydrating UF4 ? 0.10H20 has an absorption even more like that of the cubic hydrate of identical composition. The maxima of its infrared absorption spectra are at the same positions as those of the cubic hydrate. The above results thus show that the original hydrate undergoes a phase transition between the composition UF4 ? 1.00H20 and UF4 ? 0.5H20, accompanied by a structural change; dehydration products containing less than 0.5H20 per UF4 molecule cannot be rehydrated to the original UF4 ? 4/3H20. In conclusion, the authors would like to thank S. S. Batsanov for the refractometric work, recording the in- frared spectra and discussing the results, and L. A. Khripin for recording the thermograms. LITERATURE CITED 1. W. Zachariasen, Acta crystallogr., 2, 388 (1949). 2. J. Dawson, R. D'Eye, and A. Truswell, J. Chem. Soc., November 3922 (1954). 3. Yu. V. Gagarinskii and V. P. Mashireev, Zh. neorgan. khim., 4, 1246 (1959). 4. S. S. Batsanov and Yu. V. Gagarinskii, Zh. struktur. khim., 4, 387 (1963). 5. S. P. Gabuda, et al., Zh. stxuktur. khim., 5, 303 (1964). 6. Yu. V. Gagarinskii, S. P. Gabuda, and G. M. Mikhailov, ibid., p. 383. 48 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 ? ? Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 THE RELATIVE VOLATILITY OF SOLUTIONS OF HTO IN H2O (UDC 546.23:536.432.1) Ya. D. Zel'venskii, V. A. Shalygin, V. S. Tatarinskii, and D. A. Nikolaev Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 46-48, January, 1965 Original article submitted January 27, 1964 A method of simple distillation, with calculations performed according to Rayleigh's equation, was used to determine the separation factor a of solutions of HTO in H2O in the 38-100?C tem- perature range. The results obtained are expressed by the equation logna = 38.80/T ?0.0935. The data of our study are close to the separation factor values calculated on the basis of earlier measurements of the vapor pressure of T20 (see Atomnaya Energiya, Vol. 8, No. 5, p. 420, 1960). The data available in the literature regarding the values of the relative volatility (separation factor) a of so- lutions of tritium water are unreliable and contradictory. In [1] it was stated that the boiling point of hydrogen- tritium water, HTO, is lower than that of ordinary water, H20, i.e., that the former is more volatile. This con- clusion was contradicted by the data obtained by Price [2], who used the ratio of the activities of the liquid and gaseous phases to determine the separation factor of HTO and H2O solutions in the 25-80?C temperature range. According to Price, for example, at 55?C the separation factor of the HTO-H20 system is 1.13. Considerably lower values of a were found in [3] by a method using a fractionating column: a = 1.051- 1.053 at 70?C and a = 1.036 at 100?C. At the two temperatures investigated, according to [3], (aT ?1) /(aD ?1) = 1.37 ? 0.02, where aD is the separation factor of the HDO ?H2O system. Theoretically, the ratio of the enrichment factors, (aT ?1)/ (aD ?1) should be 4/3[4]. In 1960 M. Popov and F. Tazetdinov published the results of their measurements of the vapor pressure of triti- um water, T20, in the 12-95?C temperature range [5]. Using water specimens with concentrations of 83.40/o and 98.12% T20, the authors of [5] conducted their measurements by a static method. The resulting values of vapor pressure, converted to the values for 100% T20, are expressed by the generalized equation igp = 7.9957 1654.9 t ? 222 ? Popov and Tazetdinov computed the vapor pressure of HTO on the assumption that it is equal to the geometric mean of the vapor pressures of T20 and H20. At 25?C, according to [5], a 'PH20/PHTO = 1.097, instead of the value of 1.295 found in [2]. There is less disagreement between the data of [5] and [3]. In order to check and improve the accuracy of such very divergent data on the value of the separation factors of solutions of tritium water, the study described below was carried out. Experimental Procedure In order to determine the separation factor, we used the method of simple distillation, with calculations per- formed according to Rayleigh's equation. This method is relatively simple and, at the same time, sufficiently rigorous and free from unverified assumptions, which makes it possible to obtain accurate results [6]. The experi- mental technique consists in evaporating a large amount of the initial solution, under isobaric equilibrium condi- tions, until a small residual amount is left. For dilute solutions of the less volatile component, the separation factor is calculated by means of the equa- tion [6]: Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 49 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Fig. 1, Schematic of equipment used for determining separation factor: 1) Still; 2) shell of air thermo- stat; 3) heater; 4) autotransformer; 5) stirrer; 6) magnet; 7) motor; 8) water condenser; 9) drying column with zeolite; 10) manometer; 11) manostat; 12) vacuum pump; 13) condensate collector. 1g (Wo/WE) a =-- lg (Wo/WR)+ Ig (Ao/Als) (1) where Wo/Wf is the ratio of the amount of the initial solution to the amount of the final solution, and Ao/Af is the ratio of specific activities. A schematic of the equipment used for the distillation is shown in Fig. 1. The still 1 was placed in an air thermostat 2, whose temperature was kept 10-20?C higher than the boiling point of water in order to prevent con- densation of vapor on the walls of the still. The still was equipped with a magnetically driven stirrer. In accord- ance with the results of the preliminary tests, the rotation rate was kept between 80 and 100 rpm. The steam, being liquefied in be condenser 8, flowed into the collector 9. A predetermined constant pressure was maintained in the apparatus by means of a vacuum pump 12 and an automatic manostat 11. The distillation was carried out at a slow rate, in order not to disturb the vapor-liquid equilibrium. For a separation factor a which is close to unity, it is desirable to increase the ratio Wo/Wf if we wish to im prove the accuracy of the results. However, if the volume of the still is made excessively large, this will upset the conditions on which Eq. (1) is based (according to our data, the residual liquid at the end of the distillation should be not less than 1% of the original amount). We therefore carried out the distillation in two stages, with a volume reduction factor of Wo/Wf = 50-60 in each stage. To achieve this, after the first stage all of the residual liquid was poured into a smaller still, in which the second stage of the distillation was carried out. The amount of water before and after the distillation was determined by weighing. The activity of the triti- um water before and after the distillation was measured by a scintillation method. The scintillator used was a so- lution of paraterphenyl in dioxane containing an admixture of the preparation POPOP (5 g of paraterphenyl and 0.1 g of POPOP per liter of dioxane). To make the measurements, on the basis of preliminary, tests, we mixed 0.1 ml of the water to be analyzed w th 9.9 ml of a scintillating solution of the indicated composition. The resulting sOlution was poured into a glass cuvette sealed with Canada balsam to the face of the FEU-19M photomultiplier. The latter was combined with a "Volna" instrument complex. When we measured the back- ground, we poured a mixture consisting of 9.9 ml of the scintillating solution and 0.1 ml of ordinary distilled water into the same cuvette. The activity was measured with an accuracy of lob. The specific activity of the tritium water used was about 0.5 mCi/ml. 50 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Results of the Measurements factor determined 11111111111=1111111111 The separation of dilute solutions of tritium water was 1111M1111?111111111 at pressures of 50-750 mm Hg, which corresponds to varying the boiling point from 38 to 100?C. The results obtained are given below. 1111111MMIMIIII MINE1111111111N1111 7.15 11111111111?11111111111M Pressure, mm Hg 50 100 300 750 111111UMMIIIM Boiling point, ?C 38.1 52.0 76.8 100.0 M111111?111111111111110 Volume reduction in first 61.0 59.2 56.4 62.2 1113111111111111111= stage of distillation, , Wo/WI 11111111111111=1011111 410 RIRIUi Volume reduction in 52.8 61.6 52.3 52.5 second stage of distilla- tion, Wi/Wf 1?1111111111MINII Total volume reduction, 3220 3646 2940 3260 11111?111110111111111 1.0 1111111111?11111172111111111 Wo/Wf Activity, Ao, of original water, pulses per 100 sec 58,600 58,720 58,740 58,790 111111111MMIN. 11111111111111MMIk Activity, Af, of remainder 95,400 105,330 80,670 72,760 after distillation Separation factor a cal- 1079 1064 1041 1027 11111111111111111111M Fig. tion data Present 0 20 411 60 80 t.?C culated according to 2. Comparison of separa-Eq. (1) factors according to the Separation factor acalc of various studies: 0) calculated according study; A) [2]; -F) [3]; to Eq. (2). 1080 1065 1041 1027 The variation of the separation factor as function of temperature can be represented by the equation log10 a = A/T - B. Accordingly, the resulting ex- perimental data can be expressed by the general equation g a = 38.80 0.0935. As can be seen from the tables, the deviation of acak from the experimental values of a is no more than 0.01. (2) Figure 2 shows the values of the separation factors in comparison with the data obtained in other studies. As can be seen from the figure, the results of the present study are in satisfactory agreement with the data of [5], which were obtained by a different method. There is also comparatively little difference between the two values of a ob- tained in [3] by a third method. Price's data [2] are considerably higher and must be considered erroneous. The assertion of Libby and Cornog [1] that the vapor pressure of HTO is higher than that of H20 is also incorrect. LITERATURE CITED 1. W. Libby and R. Cornog, Phys. Review, 60, 171 (1941). 2. A. Price, Nature, 181, 262 (1958). 3. P. Avinur and A. Nir, Nature, 188, 652 (1960). 4. I. Bigeleisen, Tritium in the Physical and Biological Sciences, Vienna, IAEA (1962), Vol. 1, p. 161. 5. M. M. Popov and F. M. Tazetdinov, Atomnaya Energiya, 8, 420 (1960). 6. Ya. D. Zel'venskii and V. A. Shalygin, Zh. prikl. khimii, 31, 1501 (1957). 51 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 DETERMINATION OF y -RAY AND NEUTRON ABSORBED DOSE IN POLYMERS (UDC 639.121.7) F. A. Makhlis and I. M. Kolpakov Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 48-52, January, 1965 Original article submitted January 9, 1964 Together with results, a method is given for calculating the 7-ray and neutron absorbed dose in several polymers which is based on the elemental composition of the materials and the energy of the radiation. In radiochemical studies and in the radiation treatment of various materials and biological subjects, the basic physical problem is the determination of the amount of energy which is absorbed by the irradiated medium, i.e., the determination of absorbed dose. Strictly speaking, this criterion, although it is fundamental, need not be unique. Clearly, radiochemistry lacks a concept resembling the coefficient of relative biological effectiveness (RBE) in radiation biology. Actually, for the same absorbed energy, the changes in the properties of a polymer under neutron irradiation can differ from the corresponding changes under 7-irradiation [1]. However, the required quan- titative data are lacking, and we are forced to proceed solely on the basis of the magnitude of the absorbed energy, on which the final effect of the irradiation process depends to a considerable extent, other conditions being equal. Gamma Irradiation With the exception of calorimetry, there are no known experimental methods by which one might determine the magnitude of the absorbed dose for each specific material under 7 irradiation. However, the use of the calori- metric method is made difficult by the complex configuration of the absorbed dose field which is produced by ac- tual irradiation equipment. Consequently, one ordinarily resorts to indirect methods for determining the absorbed dose, using various liquid dosimetric systems of air ionization. A number of important considerations connected with the determination of absorbed dose from the results of radiation dose measurements were discussed in [2]. Knowing the magnitude of the absorbed dose D1 in a definite volume occupied by one material, one can ob- tain the corresponding magnitude of the absorbed dose D2 for another material, other conditions being equal, from the relation D (1-ta/Ot n (N/Q)2 (1) where (p a /p)1 and (p a /p)2 are the true mass absorption coefficients for y radiation in the first and second ma- terials. If the radiation dose in air is known, then the conversion to absorbed dose in the irradiated material is ac- complished on the basis of formula (1). Indeed, if it is assumed that the energy needed for the formation of an ion pair in air is 34 eV, the energy equivalent of the roentgen in air with fulfillment of electron equilibrium is 87.7 erg/g. From that, we obtain DM == 0.877 (ILta/Q)m- [rad/R], (11aMair (2) where the subscript "M" indicates the irradiated material. The true mass absorption coefficient for a multicomponent system (pa /p)syst consisting of n elements is defined as (Ra/Q)syst = ti?gi (11.70i, where gi is the fractional weight of the j-th element in the system. 52 (3) Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 1. Values of Dm for Different Elements, rad/r . '2' y -Energy, MeV 0.1 0.2 0.5 1.0 2.0 5.0 10.0 H 1.54 1.73 1.74 1.74 1.74 1.60 1.36 C 0.811 0.870 0.877 0.877 0 877 0.864 0.834 N 0.844 0.873 0.877 0.877 0.877 0.876-0.866 0 0.877 0.883 0.877 0.877 0.877 0.880 0.894 F 0.905 0.845 0,832 0.830 0.835 0.850 0.880 Na 1.072 0.870 0.837 0.840 0.843 0.877 0,832 Mg 1.255 0.912 0.864 0.852 0.870 0.910 0.990 Al 1.395 0.899 0.845 0.842 0.852 0.900 0.991 Si 1.70 0.96 0.875 0.867 0.876 0.945 1,06 S 2.28 1.01 0.884 0.870 0.877 0,968 1.100 Cl 2.55 1,02 0.850 0.832 0.840 0.940 1.11 K 3.37 1.11 0.864 0.850 0.859 0,970 1.55 Ca 4.13 120 0.894 0.870 0.877 1.01 1.213 Fe 8.34 1.60 0.864 0.819 0.820 1.0 1.261 In order to simplify the computations, values of Dm are given in Table 1 for elements frequently en- countered in radiochemistry and radiobiology. These values were obtained by using the values of ga/p in [3]. Table 2 gives values of Dm for water and a number of polymers. Neutron Irradiation For reasons which are given in [4], radiation treat- ment of materials with neutrons is possible under very limited conditions, but the calculation of absorbed dose is necessary in order to estimate the radiation resistance of materials intended for operation in nuclear reactors. In the calculations given below for the absorbed dose in polymer irradiation with thermal and fast neutrons, the assumption was made that the neutrons undergo only one collision during penetration of the irradiated material. The validity of such an assumption is verified by the following considerations. 1. The thickness of polymer products (for example, rubber products) which are irradiated usually does not ex- ceed the neutron mean free path. Thus, the thermal neutron ranges in water and polyethylene are 0.28 and 0.24 cm, respectively; for 1 MeV neutrons, they are 2.45 and 2.3 cm, and 11.0 and 7.7 cm for 10 MeV neutrons. 2. The value of the absorbed energy we obtained was close to the maximum absorbed energy for the thickness of the irradiated object. This followed from a comparison of the experimentally determined distribution of neutron absorbed depth dose in a paraffin phantom [5] with computed values which were obtained with the "first collision" approximation. The deviation between the corresponding values did not exceed 40/0 for neutron energies of 0.24- 1.0 MeV. Thermal Neutron Irradiation For the elements which are present in a majority of the materials used in radiochemistry, the following basic reactions are characteristic of interactions with thermal neutrons: N14 (n, p) C14, a = 1.76 barn, Q = 0.623 MeV, K = 0.996, 0" (n, a) C", a -= 0.5 barn, Q = 1,72 MeV, K = 0.00039, C135 (n, p) S35, a = 0.17 barn, Q = 0.62 MeV, K = 0.7543. Here, a is the cross section for the respective reaction, cm2; Q is the energy of the emitted particles, MeV; K is the isotopic fraction in a natural mixture of isotopes. The contribution from the 017(n, a)C14 reaction is negligibly small, because of the insignificant amount of 017 in a natural mixture of isotopes. One can assume that the energy transferred by the products of the other two reactions is absorbed in the immediate vicinity of the location of the neutron-nucleus interaction. Then the absorbed dose, normalized to unit neutron flux, is Du= 1.6?10-8NiKiQuaii (rad/(neut.cm-2.)]. Here, Ni is the number of atoms of element j per cm3; the subscript i identifies the\ reaction under consideration, and the subscript j, the material under consideration. For the N14(n, p)C14 reaction, Do is 7.5 ? 10-10 rad/neut-cm-2 and for the C135(n, p)S35 reaction, it is 0.239 ? 10-1? rad/neut-cm-2. In the emission of protons or a-particles, one can neglect the energy of recoil nuclei. (4) 53 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 2. Values of Dm for Water and Several Polymers, rad/R Polymer y-Energy, MeV 0.1 0.5 0.8 i 1.5 2 5 10 Polyethylene, SKEP, butyl rubber 0.91 1.0 1.00 1.0 1,0 1.0 0.97 0.91 Natural rubber 0.90 0.99 0.98 0.98 0.99 0.99 0.96 0.90 SKS, SKB 0.89 0.97 0.97 0.97 0.97 0.97 0.94 0.90 Nitrile rubber 0.88 0,96 0.96 0.96 0.96 0.96 0.93 0.88 Silicone rubber 1.2 0.95 0:94 0.94 0.94 0.94 0.95 0.97 Fluoresin, "Viton" type 0.89 0.86 0.76 0.76 0.76 0.86 0.87 0.87 Iiy2alon 1.38 0.95 0.95 0.94 0.94 0.94 0.95 0.96 Teflon 0.88 0.84 0,84 0.84 0.84 0.84 0.85 0.86 Water 0.96 0.97 0.97 0.97 0.97 0,97 0,96 0.94 TABLE 3. Values of the Absorbed Dose, Normalized to Unit Neutron Flux, for Water and Several Polymers, rad/neut-cm 10-10-2. Material Neutron energy, MeV thermal 0.1 0.5 0.8 1 2 5 10 Polyethylene, SKEP, butyl rubber 0.155 9.0 21.9 28.5 31.6 42.7 61.2 75.5 Natural rubber 0.136 8.1 20.6 25.8 28.6 38.8 55.5 .69.0 Nitrile rubber 0.67 6.74 17.2 21.5 23.9 32.5 46.7 58,7 Silicon rubber 0.124 5.15 13.3 16.6 19.1 24.9 35.8 44,4 Fluoresin, "Viton" type 0.091 2.19 4.35 5.5 6.0 8.8 12.9 18.0 Fluoresin, Kel-F type 6,10 0.77 2.67 3.43 3,74 6.04 8.53 12.6 lIy2alon 4.05 5.9 15.0 18.8 20.9 28,2 40.8 50.7 Tenon ' 0.019 1.03 1.39 1.84 1.93 7.3 5.1 9.0 Water 0.139 6.8 17.8 21.6 27.0 32.4 45.5 56.8 TABLE 4. Values ofDij Normalized to Unit Fast Neutron Flux, rad/neut-cm-2 ? 10-1? Element Neutron energy, MeV 0.1 0.2 0.5 1.0 2.0 5.0 10.0 II 60.5 92 148 205 278 390 455 C 0.512 0.96 1.94 2.93 3.95 7.41 13.1 N -- 0.64 1.03 2.00 2.86 6,42 12.0, 0 0.23 0.47 1.67 5.0 9.07 4.33 8.15 F 1.07 0.39 0.65 1.09 3.10 4.46 7.72 Na 0.11 0.49 0.35 1.20 2.01 3.77 7.55 Mg 0.25 0.35 0.77 0.83 1.91 3.24 1.86 Al 0.15 0.22 0.43 0.64 1.67 2,71 4.19 Si 0.023 0.50 0.35 1.04 1.08 2,75 4.12 S 0.13 0.18 1.19 0.47 1,04 2.07 1.6 Cl 0:025 0.064 0.16 0.33 0.90 1.99 3.05 Ca 0.005 0.022 0.11 0.23 0.46 1.15 2,30 Fe 0.23 0.05 0.09 0,13 0.28 0,70 1.01 The energy released in radiative capture can also be found with formula (4), the only difference being that Qij is understood to be the binding energy released with the addition of a neutron to the nucleus [6]. Then the y -ray energy which is absorbed by a thin layer of irradiated material of thickness r can be approximated by the relation 54 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 = Duliar, (5) where ma is the y -quantum linear absorption coefficient for the material under consideration. If a single y -quantum is emitted in neutron capture by the hydrogen nucleus and, with some approximation, by the carbon nucleus, then one can assume with a reasonable amount of accuracy that three y -quanta with energies each are emitted in the case of the materials of interest to us. The contribution of nuclear recoil energy and induced activity to the absorbed dose is negligibly small in this case [7]. Results from a computation of the magnitude of the absorbed dose, normalized to unit thermal neutron flux, are given in Table 3 for water and several polymers. Fast Neutron Irradiation The energy which is transferred to a polymer as the result of elastic scattering of fast neutrons is determined from the relation Di; =1.94-1016E E [rad/neut-cm-2], (AJ-1-1)2 (6) where o ij is the neutron elastic scattering cross section for nuclei of element j, cm2;g ? is the proportion of element j in the mixture; Aj is the atomic weight of element j; E is the neutron energy, MeV.l In hydrogenous materials, the principal process which is responsible for the transfer of energy to matter is scattering by hydrogen. The fraction of the energy which is transferred by this process is 85-90% of the total absorbed dose [8]. To facilitate computations using formula (6), values of Dij (rad/neut-cm-2) are given in Table 4 for ele- ments often encountered in radiochemistry (only elastic scattering of neutrons was considered in the calculations). The relative role of inelastic neutron scattering as compared with elastic scattering increases with increasing energy and with the atomic number of the irradiated material. Of the light elements which make up polymers, only fluorine has an inelastic scattering cross section of 1.25-0.5 barns for neutron energies of 0.5-1.0 MeV; further- more, y-quanta with energies of 0.109 and 0.197 MeV are emitted. At neutron energies exceeding 1 MeV, inelastic scattering is also possible in aluminum, magnesium, iron, and other metals found in resins. The energy of the ac- companying y -quanta is approximately 1 MeV. Therefore, the contribution of y -radiation from inelastic scatter- ing to the absorbed dose can compete with the contribution from elastic scattering only for neutron energies above 0.5 MeV and only if the proportion of metal in the composition of the irradiated material is sufficiently great. The fast neutron radiative capture cross section is hundreds of times smaller than the scattering cross section, and only a small fraction of the energy released is absorbed in light elements in the case of thin, irradiated objects; consequently, one can neglect the contribution to the absorbed dose from y -rays in (n, y) reactions. A similar conclusion is reached with regard to secondary charged particles, recoil nuclei, and induced activity. The results of absorbed dose calculations, normalized to unit fast neutron flux, are given in Table 3 for water and several polymers. With reactor irradiation, the contributions of fast and thermal neutrons to the total absorbed dose depends on the composition of the polymer. For the majority of polymers, the absorbed dose from fast neutrons exceeds by far the corresponding value for thermal neutrons. Polymers which contain a considerable quantity of fluorine are an exception. Calculations which were carried out for the irradiation of polyethylene and Kel-F in the core of the VVR-M thermal reactor [9] showed that the thermal neutron absorbed dose in polyethylene did not exceed 1%, and for Kel-F was 501o, of the fast and intermediate neutron absorbed dose. For irradiation in the vicinity of the beryllium reflector, the analogous values were 3% for polyethylene and 100% for Kel-F. In conclusion, we wish to thank E. N. Smagin for valuable discussions. LITERATURE CITED 1. A. Charlesby, Nuclear Radiations and Polymers [Russian translation] (Moscow, Izd-vo inostr. lit., 1962). 2. Yu. V. Sivintsev, Atomnaya energiya, 9, 39 (1960). 55 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 3. R. Berger, Rad. Res., 15, 1 (1961). 4. F. A. Makhlis and A. Kh. Breger, Kauchuk i resina, No. 8, 18 (1964). 5.. A. M. Kogan, et al., Atomnaya energ i y a , 7, 351 (1959). 6. B. Price, C. Horton, and K. Spinney, Nuclear Radiation Shielding [Russian translation] (Moscow, Izd-vo inostr. lit., 1959). 7. M. L Shal'nov, Neutron Tissue Dose [in Russian] (Moscow, Atomizdat, 1960). 8. R. Eger, Dosimetry and Radiation Shielding [in Russian] (Moscow, Gosatomizdat, 1961). 9. V. V. Goncharov, et al., Proceedings of the Second International Conference on the Peaceful Use of Atomic Energy [in Russian] (Dokl. soy. uchenykh, Moscow, Atomizdat, 1959), Vol. 2, p. 293. 56 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 INCREASING THE DEPTH OF PROSPECTING FOR CONCEALED URANIUM ORE BODIES BY MEANS OF THE PRIMARY AUREOLE (UDC 550.8/546.791) S. V. Grigoryan Translated from Atomnaya Energiya, Vol. 18, No. 1, pp. 52-57, January, 1965 Original article submitted February 14, 1964 A study of the primary uranium aureoles in one hydrothermal deposit showed that, by using sys- tematic chemical analysis (sampling analysis of the mobile uranium content), the primary aureole can be more fully revealed, and the depth of aureole prospecting for hidden uranium ores thus increased. It is shown that this method of analysis is most effective for elements whose aureoles are of small extent and intensity, owing to high background content in the country rocks, low ore content, etc. Recent investigations have shown that uranium ores of hydrothermal origin are surrounded by aureoles of urani- um and other elements; these can be used in prospecting for concealed uranium ore bodies [1, 2]. The main ele- ments used as indicators of uranium ores are those forming the most extensive aureoles. A study of the distribution characteristics of elements in a number of hydrothermal uranium deposits showed that, in any cross section, as one moves laterally away from the ore, the concentrations of uranium and its ac- cessory minerals decrease in geometric progression [2, 3]. The distribution of any element around the ore can thus be represented by a graph (Fig. 1). Since the concentration decreases in geometric progression, the distribution is represented by a straight line in semilog coordinates. In practice, aureoles are mapped with respect to the least anomalous content of an element, just as in prospec- ting for deposits the limits of commerical ores are determined from the boundary content of the required component [2]. In Fig. 1 the half-width of the aureole corresponds to the intercept cut off by the distribution graph on the horizontal line representing the least anomalous content. In actual fact the aureole will be considerably wider; with increasing distance from the ore body, the element concentrations introduced during ore formation decrease, not down to the least anomalous content, but to the content in the interstitial solutions which saturated the rocks when the ore and aureole were formed. This is due to the fact that, when the primary aureoles were formed, the element contents were equalized between the ore-bearing solutions and the interstitial solutions saturating the country rock. The aureoles formed by this equalization will also have the maximum possible size. Let us call them "true," and those mapped from the least anomalous content "apparent." The difference between these aureoles will increase with the least anomalous content of any element. The vertical extent of the apparent aureoles is much less than that of the true ones. The relation between the dimensions of the apparent and true aureoles was studied in detail on one hydrothermal uranium deposit occurring in granites. The deposit appeared as numerous fine tarry fluorite veinlets, forming stockwork ore bodies of various sizes, elongated in a northwesterly direction. The dip of the ore bodies is steep. The boundaries of the commerical ores are known from the results of assays. In mapping the aureoles of the elements, the least anomalous uranium content was determined from the geo- chemical background value and standard deviation in the rocks [4]. The method of mapping is given in more de- tail in [2]. 57 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Least anomalous content (apparent aureole) 71i Geochemical background /0 a.) 8 : ? 6 0 U 4- ? 3 Content of element in interstitial solutions (true aureole) Cs I I 2 ? 4,1 Distance from ore deposit, m Ore body 444,. ? 1%. Fig. 1. Distribution of a chemical element around an ore deposit. To determine the vertical extents of the apparent and true aureoles, we drew graphs of the linear productivityl of the uranium aureoles in a vertical cross section through the ore body (Fig. 2). The linear productivity of the aureole at a given horizon was found by multiplying the mean uranium content by the width of the aureole. The abscissae in Fig. 2 represent the linear productivity calculated from the aureole of mobile uranium [2, 4]. The or- dinates represent the distance of the sampling horizons from the ore body, plotted on a linear scale. It is seen that the linear productivity of the uranium aureole is represented by a nearly straight line. To determine the vertical extent of the uranium aureole above the ore body, this line must be produced to meet the line corresponding to the least anomalous linear productivity. (This productivity corresponds to the aureole of minimum breadth and with mean content equal to the least anomalous content.) The apparent minimum breadth is equal to the sampling inter- val (5 m). However, we took it as 10 m, since the aureole must be mapped by at least two samplings. Conversely, samples in which the high concentration of elements is caused by random factors (e.g., background variation) may be mistaken for an aureole. For the ore deposit described, the least anomalous uranium content was taken as 0.00161, which exceeds the geochemical background (GB) of the surrounding granites, 0.0007%, by twice the standard devia- tion (GB + 2a). The minimum anomalous productivity will thus be 0.0016 ? 10 = 0.016m %. The extent of the uranium aureole above the ore body was found by extrapolation to be 180m (Fig.2). That of the true aureole above the same ore body was determined in a similar manner. It is equal to the distance from the ore body to the point of intersection of the straight line representing the productivity with a vertical line corre- sponding to the linear productivity (given by the product of the uranium content in the interstitial solutions and the minimum aureole breadth, 10 m). Graphical determination shows that, in this cross section, the true aureole ex- tends about 500 m above the ore body, i.e., more than 2.5 times the extent of the apparent aureole. From the above it is clear that the depth of an aureole survey can be increased by decreasing the least anom- alous productivity of the aureole (preferably by reducing the least anomalous content, since decreasing the sampling interval leads to a considerable rise in working expenses). To choose an appropriate method of systematic sample analysis (the essential point of the problem), we studied the forms of uranium occurrence within the aureole. The uranium content of an aureole accumulates from the primary distribution in the country rock (the Clarke content) and from uranium introduced by the hydrothermal solutions and forming the aureole. Many investigators, having studied the occurrence of primary uranium in magmatic rocks, remark that part of the uranium can be extracted from the rocks by relatively weak solvents, but the whole is extracted only by breaking down the crystalline structures of the rock-forming and accesory minerals. Tauson [6] distinguishes be- tween two main forms of uranium in granitoids: non-isomorphic and isomorphic. He assigns the name "non-iso- iThe idea of linear productivity of a diffuse aureole is due to A. P. Solovov [5]. 58 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 TABLE 1. Degree of Uranium Extraction by Na2CO3 from Pitchblende and Sooty Uraninite Sample No. Uranium mineral Total U content, oh Extraction into so- lution, oh A-3 Pitchblende 3.6.10-4 75 A-12 * 9.10-4 89 A-5 # 8.10-4 75 C-89 * 0.2 90 C-1 # 0.31 74 A-9 Uraninites 3.6?I0-4 100 s400 .0 z jog morphic" to that part of the uranium which can be extracted by boiling with 507o ammonium carbonate solution containing hydrogen peroxide. In these con- ditions the crystalline lattice of the rock-forming and accessory minerals is practically undisturbed. Tauson considers that, with this method of analysis, the urani- um from separate minerals (such as oxides) also passes into solution. Experiments have shown that granitoid rocks have a negligible content of non-isomorphic uranium. Our investigations showed that uranium, intro- duced into the country rocks during the formation of an ore body, occurs mainly as the individual minerals pitchblende and the sooty uraninite group of oxides [2]. In order to bring. the latter into solution without destroying the crystal lattices of the rock-forming and accessory minerals, weighed samples were boiled for 30 min with 20h soda solution containing hydrogen per- oxide as oxidizer. o To determine the degree of uranium extraction 0 0 from pitchblende and oxides, an artificial mixture of I I . . "8 - . these minerals with pure quartz sand was treated as 0 E ? Zoo 02 Surface ?E o s o . . described. The total uranium content of the samples was determined by melting weighed portions with soda. O 2 Horizon I ..1). II % The results showed that pitchblende and uraninite are satisfactorily dissolved by soda solution (see Table 1). w foo "8 ...47 It should be noted that Na2CO3 solution ex- Horizon III -8 % N tracts only a tiny part of the primary, so-called mobile Horizon IV .ct Yi VI % . uranium from the country granites outside the limits of . the aureole [6-8]. The least anomalous content of - ? . Horizon V ... ' o ? ?r .0 0 r- mobile uranium for the granites of the deposit de- 4 HorrizonVI -o. ; scribed is equal to 2 ? 10-407, which also exceeds the o geochemical background (0.000040h) by twice the stand- cool o.00t aot of to to ard deviation (GB + 20). The (linear) least anomalous ' productivity for mobile uranilum is 0.002 m ?A As shown by the graphical determinations, the extent of Fig. 2. Extent of uranium aureoles above an ore body. the mobile uranium aureole above the ore body is I, II, III) Linearized changes of productivity of urani- more than 11/2 times that revealed by measuring the urn aureole in a vertical direction. total U content (see Fig. 2). In other words, the above method of systematic sample analysis considerably increases the effective depth of prospecting for concealed ore bodies by means of uranium aureoles. This conclusion, based on graphical construction, is confirmed by mapping the aureoles in a vertical Cross section through the ore body. Figure 3 shows the uranium aureoles around a concealed ore body. The graph of the linear productivity in this cross section is given in Fig. 2. The aureole is mapped twice, according to sample assays of mobile uranium and of total U content. To reduce the number of analyses, the total U content was determined only for the samples taken at the surface and first horizon (see Fig. 3). The mobile uranium aureole is the wider and longer of the two. The mobile U content of all the surface samples was anomalous, whereas that of total U in the same samples was less than the minimum anomalous value ? i.e., no total-U aureole is observed at the surface of this cross section. At the first horizon the width of the U aureole is small. It follows from Fig. 3 that, if the total U aureole were used in prospecting for this ore body, it would not be located. A Horizon II 1/4 1/4 ? HorizonVII ?..? ? . . Linear productivity of uranium aureol, m% 59 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 a Scale 20 0 20 'Cm 1=1 ZNN1 4[= 5=I 6LJ 7E-1 8 Fig. 3. Primary uranium aureoles; a) Aureole of total U content; b) aureole of mobile urani- um. 1) Granites; 2) ore body; 3) faults; 4) workings and boreholes; 5) sampling intervals. Uranium contents, 0/0; 6) 0.002-0.0015; 7) 0.0015-0.005; 8) >0.005. a Scale 10 0 10 20 30 40.m Wimi==imamichnowni CSO 2R 3 MOE F-1 VW 6M 7 Fig. 4. Primary uranium aureoles; a) Of total U content; b) of mobile U. 1) Ore body; 2) working and boreholes; 3) sampling intervals. Uranium contents, 010: 4) 0.0002-0.0005; 5) 0.0005-0.0015; 6) 0.0015-0.005; 7) > 0.005. Figure 4 shows the uranium aureoles around two contiguous concealed uranium ore bodies. The latter are revealed at the surface by only two samples, if the least anomalous total U contents are used. For mobile U, how- ever, all the surface samples fall within the aureole. Furthermore, at the surface directly above the ore bodies an anomalous field of mobile uranium is found, its width exceeding that of the aureole of total U content (see Fig. 4). 60 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 The method of analysis for mobile uranium is easier than that for the total content: It facilitates the bulk sample analysis which is so necessary in geochemical prospecting. We can expect this systematic analysis method to be most effective for elements with small apparent aureoles of low intensity. Figure 2 shows that the use of sys- tematic analysis increases the extension of the apparent maximum-productivity aureole of uranium by a factor of 1.2 (curve III), while the extension of the aureole with minimum productivity is increased nearly ten times (curve II). The following conclusions may be drawn. 1. The use of systematic chemical analysis (determination of the mobile uranium content) permits fuller ex- posure of primary uranium aureoles and increased depth of prospecting for concealed uranium ores by means of pri- mary aureoles. 2.- The use of systematic chemical analysis will be most effective for elements with apparent aureoles which are, owing to high background concentration in the country rocks, low concentration in the ores, etc., of small dimen- sions and intensities. 3. By means of systematic chemical analysis it will be possible to find primary aureoles for a number of new elements which cannot be discovered by total content determinations, as their background levels are high and, con- sequently, so are their least anomalous contents. LITERATURE CITED 1. A D. Kablukov and G. L Vertepov, Geologiya rudnykh mestorozhdenii, No. 2, 20 (1960). 2. E. M. Yanishevskii, et aL, Endogenous Diffusion Aureoles of Certain Hydrothermal Deposits [in Russian] (Moscow, Gosgeoltekhizdat, 1963). 3. H. Morris and T. Lovoring, Econ. Geol., 47, 7, 685 (1952). 4. S. V. Grigoryan, Razvedka i okhrana nedr, No. 2, 26 (1964). 5. A. P. Solovov, Razvedka i okhrana nedr, No. 4, 48 (1952). 6. L. V. Tauson, Geokhimiya, No. 3, 9 (1956). 7. L. S. Evseeva and A. I. Perel'man, Geochemistry of Uranium in the Supergene Zone [in Russian] (Moscow, Gosatomizdat, 1962). 8. A. S. Saukov, Radioactive Elements of the Earth [in Russian] (Moscow, Gosatomizdat, 1961). 61 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 LETTERS TO THE EDITOR PHASE STABILITY OF PARTICLE BLOBS IN ACCELERATORS WITH AUTOMATIC CONTROL (UDC 621.384.60) A. Zhil'kov Translated from Atomnaya Energiya, Vol: 18, No. 1, pp. 58-59, January, 1965 Original article submitted May 8, 1964 The characteristic equation for a system of particles performing nonlinear synchrotron oscillations in an ac- celerator with automatic correction of the accelerating field's frequency (4)0= qws with respect to the beam was derived in [1]; +03 d/0 K _hk2Th 1=2ni, deQ(e) E co?kg2 (e) ? de ?Co (1) Here, u) is the natural oscillation frequency in 11(0) units, where 0(o) is the frequency of linear synchrotron oscilla- tions; e and Q(e) are the energy and the frequency of nonlinear synchrotron oscillations; fo(e) is the equilibrium distribution function of the blob; K(w) is the system's reaction to the shift of the point blob with respect to u and 9 at the frequency w; 9k and K_k are the coefficients of the expansions of 9(c, T) and K(e, T) as periodic func- tions of the variable T ,)u (e,) Wi in Fourier series at the section Es ( u 2n Q (0) ?dE Q(e) 12(e)) ' ? eV sin Ws .1 cos (E) Certain general stability criteria for small-size blobs were found in [1] under the assumption that the automatic control system corrected the position of the blob's center of gravity. Here, we shall investigate the stability of small blobs in the presence of a system of automatic control with respect to second moments and also the stability of higher moments. Consider the stability of a small blob [1] in the case where feedback with respect to the blob dimensions is provided [2], i.e., when not only the center of gravity, but also deviations of the blob dimensions with respect to vo their equilibrium values 4,(1= = e;W0 = 0 are corrected. The automatic control system is described by the equation K (w) = K (0)) u K2 (C)) (P+ K3 (0))('Y?i'0) + K4 (a)) (17 ?TM+ K5 MI 4 7, (2) where Ki(w) are the transfer constants of the feedback circuits with respect to the corresponding blob characteristics, and 4, = ? 90)2 , V = u2 and W = (9 ? 90)u are the characteristics of the blob dimensions. By substituting Eq. (2) in (1) and using the dependence of the phase 9 on the time T [3] where (T) cos QT ctg cps e L 8 ctg cps cos 2QT 8 Y28 2 m 6 96 1 5 Q 1 --8- e (1+ ?3 ctg2ps) , (1 + ctg2 Ts) cos 3QT + (3) we find the characteristic equation for determining the natural oscillation frequency of the blob; 62 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 a 1 ';Tr?max -411.6max. Fig. 1. Regions of stability (I) and instability (II) of a blob in the plane of the convex variable El. a) Linearly falling equilibrium distribution; b) linearly rising equilibrium distribution. correcting don cycle the center of gravity. This makes it possible by choosing suitable transfer constants K3 K while 5 s -=???:(J. (1+ vtg2ip,3) , (5 1 + -3- ctg2 cps ) (3 + 2iK2-1C1) 4i ctg cps [ et% (Ps (K2+ 2iKi)? + K3 iK5 e2 de j =C dfo 3 de (02_4+e ( ?5 ctgq),) 3 (4) ' (5) ce=2 (6) If we consider that the distribution is uniform [1], we obtain the solution of (4) in the following form: co-2 ?1 (1 +1 ctg2 cps) i ctg cps-8 [ etg6 (Ps (K2+2iKi)?K4+ K3 i K5 3+2iK2?Ki --e_cmax where 2 co-2'M The results obtained in [1], where the parameter g denotes the expression defined by (5), remain valid for a triangular distribution. Consequently, the introduction of feedback with respect to the blob dimensions produces the same effect as the consideration of a slight nonlinearity in to secure the blob's stability throughout the entire accelera- 4, and K5. We shall analyze the stability of a blob described by the three first distribution moments. In - (1 5 + ? ctg2q) {8 + 3iK2-1Ci?ie [K2 (0.94 +0.12 ctg2 TO+ ilf (0 8+0.1 ctg2 cps)? 3 8} ?3.2 ctg cps (K3?K4-1.5iK5)]} this approximation ( 0 3 (8) (0=3. (9) (10) 3 c t g2 Ts) [4. (1+ ctg2 cps) (K2+ 3iKi)?ctg cps (21(4?K3+ 1.5iK5) For a uniform distribution, we find from expression (1) is2 [-h-(1 + ctg2 cps) (K2+ 3iKi) + cg cps kK3 (1 6tg2 21C4-1.5iK5 .co (1+ ?5 ctg2cps + cps) 4 3 8 8+3iK2-1Ci If we use a triangular distribution, the characteristic equation assumes the following form: 1 y 1 max 4(9...2) Y-= 5 9emax (1 + ctg2 cps ) ? The logarithm branch is chosen in such a manner that the logarithm value for Im y > 0 lies on the first sheet. The figure shows the stability region for the blob, plotted in the plane of the convex variable g for the linearly falling (a) and the linearly rising (b) distributions. 63 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6, The stability criteria for blobs characterized by higher moments are established in a similar manner. It should be mentioned that it is rather difficult to consider the stability of higher moments by using the method of moments [1]. We shall formulate a general statement concerning the stability of higher moments for a small blob. The char- acteristic equation for the k-th moment is given by it--2=gtfit. (1+-35-- ctg cps) , where gk_2 is a (k -2)-power polynomial with respect to &max, and j dfo eh de 4(k2_2) h 3 de e?x ' 5 k2emax (1?-3-ctg2cPs) (12) For a uniform distribution (natural damping is absent), the correction for the frequency of the k-th moment ak that is secured by the automatic control system is proportional to the (k -1)th moment of the equilibrium distribu- tion function [al, ehm-atx = 810 (e) de] . By calculating Jk for any other equilibrium distribution function, we readily see that the order of natural damping for the k-th moment is ekm-a2x. Therefore, the effect con- nected only with the automatic control system, which is a small quantity of a higher order with respect to &max, is less important. Thus, the stability of higher moments is basically determined by natural damping, i.e., by the form of the blob's equilibrium distribution function. For small blobs, the moments rapidly decrease with an increase in their numbers. The stability (instability) regions of higher moments given by Eq. (11) became narrower. In prac- tice, the stability criterion for small blobs depends on the first, the second, and, perhaps, the third, moments. Higher moments apparently have no significance. At high frequencies co >> 1, the stability of small blobs is determined by the criterion for the stability of the center of gravity. The above results naturally cannot be applied to blobs whose dimensions cannot be considered as being small, since expansions of the type (3), cannot be used in this case. The author acknowledges his deep gratitude to A. N. Lebedev for his guidance in the work and to A. A. Kolomenskii for the discussion. LITERATURE CITED 1. E. A. Zhil'kov and A. N. Lebedev, Present issue, p. 22. 2. H. G. Hereward, Proc. of Intern. Conf. on High Energy Accel. Brookhaven (1961), p. 236. 3. L. D. Landau and E. M. Lifshits, Mechanics [in Russian] (Moscow, Fizmatgiz, 1958), Paragraph 29. 4. E. A. Zhil'kov, Pribory i tekhnika jksperimenta, No. 1 (1965). 64 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 SIMPLE METHOD FOR MEASURING THE FREQUENCY OF FREE TRANSVERSE OSCILLATIONS IN CYCLOTRONS (UDC 621.384.612) S. A. Khe,ifets and. S. K. Esin Translated from Atomnaya nergiya, Vol. 18, No. 1, p. 60, January, 1965 Original article submitted January 29, 1964 The usual method I'm measuring the frequency of transverse oscillations of particles in cyclotrons consists in exciting coherent betatron oscillations and measuring the frequency of the voltage induced at the signal electrodes [1, 2]. The energy scattering of particles in the circulating beam leads to a rapid disappearance of the coherent part of oscillations. If it is necessary to measure the frequency fairly often (in comparison with the damping time), repeated pulse excitations of oscillations may result in a considerable increase in the beam's aperture. We shall demonstrate that the frequencies of transverse oscillations can also be measured by introducing local perturbations of the equilibrium orbit which are constant in time. We shall assume that a certain perturbation exists at the azi- muth so. In this case, the right-hand side appears in the equation describing transverse oscillations; where d2x + (02 (s) x =llso (s), ds2 f 0, 0 ;:: Shield i :?..4.c 10 20 30 40 Fig. 1. Density of sources of capture y -rays in steel. qVo ) Without a boron screen; ql ) with a boron screen. V 5 0 5 10 15 pt Fig. 2. Dependence of the reduction factor for cap- ture y -radiation secured by the boron screen on the "optical" thickness pt (in mean free path units) in the shield. N Experimental result given in [2]. reducing the strength of sources of capture y rays in the steel mass, a boron-containing screen also influences the form of their angular distribution at the steel surface. In this case, the angular distribution assumes a predominantly cosine character. It is well known that sources with equal strength, but different angular distributions (all other con- ditions being equal), produce different dose rates. This is due to the fact that, along the path of attenuation of y - rays in the shield, isotropic sources are transformed into anisotropic sources due to intensive filtration, so that the reduction factor for the dose of capture y radiation (the blocking factor) constitutes a function of the "optical" thickness pt (measured in mean free path units) between the sources and the detector. For the confirmation of what has been said above, we calculated the spatial and energy distributions in steel screens and the reactor vessel by using the 18-group meth- od in the P2-approximation by means of an electronic com- puter for reactors with and without a boron-containing screen. Figure 1 shows the spatial dependence of the density of capture y -rays, calculated for the entire thickness of the steel mass with and without a boron screen. By using these sources of capture y -rays, we can readily calculate the dose rate produced by these sources for different "optical" thicknesses pt (measured in mean free path units). The ratio of the dose rate produced by a source without the boron screen for the thickness pt to the dose rate produced by the same source with a boron screen characterizes the 'advantage provided by the boron screen. This factor is given by S0 q7(r) _.:_vy-sou (r') dr' 0 4rc I r 12 dV 4;0/.12 dV Figure 2 shows the reduction factor for the dose produced by capture y -radiation in dependence on the "optical" thickness pt between the source and the detector. For pt 0.5, the calculated value of the above factor is in satis- factory agreement with the experimental value given in [2, 3]. If a boron screen is provided in the reactor structure, the screen reduces the dose by a factor of 13(pt) for the thickness pt. This factor can be useful in selecting shield variants, where the dose at the point pt in question is cal- culated without considering a boron screen and is theh reduced by'a factor Of 3 if a boron screen is to be provided. It should be noted that the factor 8 depends on the composition and thickness of the steel screens and the re- actor vessel as well as on the spec rum of incident neutrons. 81 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 The authors are deeply grateful to S. G. Tsypin and Yu. A. Kazanskii for their interest in the work and critical remarks. LITERATURE CITED 1. G. I. Marchuk, Methods for Calculating Nuclear Reactors [in Russian] (Moscow, Gosatomizdat, 1961). 2. A. G. Bakov, et al., Atomnaya Energiya, 13, No. 7 (1962). 3. D. L. Broder, et al., Atomnaya nergiya, 8, 49 (1960). 82 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 Declassified and Approved For Release 2013/09/24: CIA-RDP10-02196R000700010001-6 DEPENDENCE OF THE DENSITY OF RADIATION DAMAGE TO THE REACTOR VESSEL ON THE COMPOSITION OF THE FERRO-AQUEOUS THERMAL SHIELD (UDC 621.039.553.3) K. K. Popkov and S. M. Rubanov Translated from Atomnaya ?Energiya, Vol. 18, No. 1, pp. 70-71, January, 1965 Original article submitted January 15, 1964 The present article is concerned with an investigation of ferro-aqueous mixtures with different compositions for thermal shields from the point of view of the effect of the mixture composition on the radiation damage to the reactor vessel caused by neutrons. For this purpose, we considered a primary-shield composition consisting of 25cm Fig. 1. Geometry of the compo- sitions. 1) Reactor core; 2) ferro-