An Introduction to DESC: a Stellarator Equilibrium Solver

Total Page:16

File Type:pdf, Size:1020Kb

An Introduction to DESC: a Stellarator Equilibrium Solver DESC: An Efficient Stellarator Equilibrium Code Daniel Dudt1 Egemen Kolemen1,2 1Princeton University 2Princeton Plasma Physics Laboratory APS DPP November 9, 2020 1 Overview & Motivation DESC1 is a pseudo-spectral stellarator equilibrium solver that: 1. Uses global spectral methods with Fourier & Zernike basis functions • Properly resolves the magnetic axis • Minimizes the system dimensionality • Gives a global solution (no interpolation between flux surfaces) 2. Solves force balance directly in real space (instead of the energy principle) • Avoids numerical issues at rational surfaces • Allows for perturbations to easily search the equilibrium solution space 3. Is written modern Python with high-level structure • Easy to use and extend the code for individual applications • Designed for stellarator optimization: automatic differentiation, GPUs, etc. 4. Currently solves fixed-boundary equilibria with nested flux surfaces APS DPP November 9, 2020 1D. W. Dudt, and E. Kolemen, Phys. Plasmas 27 102513 (2020). 2 “Inverse” Equilibrium Problem • Computation domain is the straight field- line coordinate system (휌, 휗, 휁) • Free variables are the flux surface shapes 푅(휌, 휗, 휁) & 푍(휌, 휗, 휁) • Problem is to find the flux surfaces that satisfy the equilibrium conditions: 푱 × 푩 = 훁푝 훁 × 푩 = 휇0푱 훁 ⋅ 푩 = 0 APS DPP November 9, 2020 3 Fourier-Zernike Basis Set • Discretize flux surfaces with global Fourier- Zernike2 spectral basis functions: 푚 푛 푅 휌, 휗, 휁 = ෍ 푅푙푚푛풵푙 휌, 휗 ℱ (휁) 푚 푛 Z 휌, 휗, 휁 = ෍ 푍푙푚푛풵푙 휌, 휗 ℱ (휁) • Inherently satisfies analytic boundary conditions at the magnetic axis 2 • Number of basis functions scales as 푀 푁Τ2 (about half as many terms as other methods) 푀 = 2 APS DPP November 9, 2020 2F. Zernike, Mon. Not. R. Astron. Soc. 94, 377–384 (1934). 4 Magnetic Field in Flux Coordinates • Assume3 nested flux surfaces: 푩 ⋅ 훁휌 = 0, and Gauss’s law: 훁 ⋅ 푩 = 0 휕휌휓 휕휌푅 푩 = 휄풆휗 + 풆휁 휋 푔 풆휌 = 0 휕휌푍 푩 휌, 휗, 휁 = 푩 푅 휌, 휗, 휁 , 푍 휌, 휗, 휁 , 휄(휌) 휕휗푅 풆휗 = 0 • Using Ampere’s Law: 훁 × 푩 = 휇0푱 휕휗푍 휕 퐵 −휕 퐵 휕 퐵 −휕 퐵 휕 퐵 −휕 퐵 휕휁푅 퐽휌 = 휗 휁 휁 휗, 퐽휗 = 휁 휌 휌 휁, 퐽휁 = 휌 휗 휗 휌 휇0 푔 휇0 푔 휇0 푔 풆휁 = 푅 휕 푍 푱 휌, 휗, 휁 = 푱 푅 휌, 휗, 휁 , 푍 휌, 휗, 휁 , 휄(휌) 휁 APS DPP November 9, 2020 3W. D. D’haeseleer et al., Flux Coordinates and Magnetic Field Structure, (Springer-Verlag, 1991). 5 Force Balance Equations 4 • MHD force balance error : 푭 ≡ 푱 × 푩 − 훁푝 = ퟎ 푅 풙 = 푙푚푛 푍푙푚푛 • Substitute in 푩 and 푱: 푭 = 퐹휌훁휌 + 퐹훽휷 휁 휗 휗 휁 ′ spectral transform 퐹휌 = 푔 퐵 퐽 − 퐵 퐽 − 푝 휌 푅 → 푅 휌, 휗, 휁 퐹훽 = 푔퐽 푙푚푛 푍 → 푍(휌, 휗, 휁) 휷 = 퐵휁훁휗 − 퐵휗훁휁 푙푚푛 훁 ⋅ 푩 = 0 훁 × 푩 = 휇0푱 • Form scalar equations: 푓휌 = 퐹휌 훁휌 푓훽 = 퐹훽 휷 푩 휌, 휗, 휁 푱 휌, 휗, 휁 • An equilibrium is a solution to the system of equations 퐉 × 푩 = 훁푝 풇 풙 ≈ ퟎ, solved at a given set of collocation points 푓휌 휌, 휗, 휁 풇 풙 = 푓훽 휌, 휗, 휁 APS DPP November 9, 2020 4S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553–3568 (1983). 6 Convergence: Heliotron 훽 ≈ 2% dimension ≈ 2 × 105 dimension ≈ 5 × 103 APS DPP November 9, 2020 7 Error: Heliotron 훽 ≈ 2% 푭 훁푝 휌 = 0.5 Δ푅푖 휀푥 = Δ푍푖 2휋 2휋 1 휀퐹 = න න න 푭 2 푔푑휌푑휗푑휁 0 0 0 APS DPP November 9, 2020 8 Equilibrium Perturbations • 1st–order Taylor expansion about an equilibrium solution: 휕풇 휕풇 풇 풙 + Δ풙, 풄 + Δ풄 = 풇 풙, 풄 + Δ풙 + Δ풄 풄 = 휕풙 휕풄 input parameters: 0 0 • pressure profile −1 휕풇 휕풇 • boundary modes Δ풙 = − Δ퐜 • etc. 휕풙 휕풄 • The new equilibrium solution for any perturbation Δ퐜 is trivial to approximate: 풙∗ = 풙 + Δ풙 휕풇 −1 Jacobian matrix 휕풙 • Can be used to find solution branches in parameter space was already computed to solve equilibrium • Has been extended to 2nd–order approximations APS DPP November 9, 2020 9 3D Boundary Perturbation • Perturbing an axisymmetric solution gives an accurate stellarator equilibrium! APS DPP November 9, 2020 10 Quasi-Symmetric Perturbations • Define a measure of quasi-symmetry (no Boozer coordinate transform needed!) 품 풙, 풄 ≡ 훁휓 × 훁퐵 ⋅ 훁 푩 ⋅ 훁퐵 • 1st –order Taylor expansion about an equilibrium QS solution: 휕품 휕품 품 풙 + Δ풙, 풄 + Δ풄 = 품 풙, 풄 + Δ풙 + Δ풄 0 0 휕풙 휕풄 −1 휕품 휕풇 휕풇 휕품 = − + Δ풄 휕풙 휕풙 휕풄 휕풄 • Resulting eigenvalue problem: 푮Δ풄 = ퟎ • Eigenvectors of 푮 corresponding to 휆 = 0 are perturbations that preserve QS APS DPP November 9, 2020 11 Summary DESC is a stellarator equilibrium solver with the following advantages: • Properly resolves the magnetic axis • Minimizes the system dimensionality • Gives a global solution (no interpolation between flux surfaces) • Avoids numerical issues at rational surfaces • Allows for perturbations to easily search the equilibrium solution space • Easy to use and extend the code for individual applications • Designed for stellarator optimization: automatic differentiation, GPUs, etc. APS DPP November 9, 2020 12 Future Development • Improved performance, user interface, documentation • Quasi-symmetry optimization • Ideal MHD stability calculations • Free-boundary equilibria • Magnetic islands & stochastic regions Repository: https://github.com/ddudt/DESC Publication: D. W. Dudt, and E. Kolemen, Phys. Plasmas 27 102513 (2020) APS DPP November 9, 2020 13 References 1D. W. Dudt, and E. Kolemen, Phys. Plasmas 27 102513 (2020). 2F. Zernike, Mon. Not. R. Astron. Soc. 94, 377–384 (1934). 3W. D. D’haeseleer, W. N. G. Hitchon, J. D. Callen, and J. L. Shohet, Flux Coordinates and Magnetic Field Structure, Springer Series in Computational Physics, edited by R. Glowinski, M. Holt, P. Hut, H. B. Keller, J. Killeen, S. A. Orszag, and V. V. Rusanov (Springer-Verlag, 1991). 4S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553–3568 (1983). APS DPP November 9, 2020 14 Bonus Slides APS DPP November 9, 2020 15 PEST1 Flux Coordinates 푍 휁 = 휙 푅 휆 휙 휓 휌 = ൗ 휓푎 휃 푍 푅 휗 = 휃 + 휆 휃, 휙 toroidal coordinates: (푅, 휙, 푍) straight field-line coordinates: (휌, 휗, 휁) APS DPP November 9, 2020 1R. C. Grimm, J. M. Greene, and J. L. Johnson (Academic Press, Inc, New York, 1976) pp. 253–280. 16 Axisymmetric Results: “D-shape” 훽 ≈ 3% 푭 훁푝 휌 = 0.5 APS DPP November 9, 2020 17 Axisymmetric Results: “D-shape” 훽 ≈ 3% • Accuracy metrics: Δ푅푖 휀푥 = Δ푍푖 2휋 2휋 1 휀퐹 = න න න 푭 2 푔푑휌푑휗푑휁 0 0 0 APS DPP November 9, 2020 18 Boundary Condition: Magnetic Axis • An analytic function expanded near the origin of a disc must have a real Fourier series of the form1,2: 푚 2 4 푓 휌, 휗 = ෍ 휌 푎푚,0 + 푎푚,2휌 + 푎푚,4휌 + ⋯ cos 푚휗 푚 푚 2 4 + ෍ 휌 푏푚,0 + 푏푚,2휌 + 푏푚,4휌 + ⋯ sin 푚휗 푚 • The Zernike polynomials inherently satisfy this condition! • Reduces the number of variables by eliminating the unnecessary high- frequency modes near the axis • No additional boundary condition equations need to be solved 1J. P. Boyd and F. Yu, J. Comput. Phys. 230, 1408–1438 (2011). APS DPP November 9, 2020 19 2H. R. Lewis and P. M. Bellan, J. Math. Phys. 31, 2592–2596 (1990). Boundary Condition: Last Closed Flux Surface • Fixed-boundary surface is given as: 푅푏 = 푅푏(휃, 휙), 푍푏 = 푍푏(휃, 휙) • Last closed flux surface is evaluated as: 푅ȁ휌=1 = 푅(휗, 휁), 푍ȁ휌=1 = 푍(휗, 휁) • Introduce 휆(휃, 휙) to convert between coordinates: 휗 = 휃 + 휆(휃, 휙), 휁 = 휙 푅ቚ = ෍ 푅푚푛ℱ(휗, 휁) ⟹ 푅ቚ = ෍ 푅෨푚푛ℱ 휃, 휙 휌=1 푚,푛 휌=1 푚,푛 푍ቚ = ෍ 푍푚푛ℱ(휗, 휁) ⟹ 푍ቚ = ෍ 푍෨푚푛ℱ 휃, 휙 휌=1 푚,푛 휌=1 푚,푛 ෨ 푏 ෨푏 • Boundary condition: σ푙 푅푙푚푛 = 푅푚푛 σ푙 푍푙푚푛 = 푍푚푛 APS DPP November 9, 2020 20 Collocation Nodes • The computational grid is a finite set of discrete points (휌푖, 휗푖, 휁푖) • The force balance errors 푓휌(휌, 휗, 휁) & 푓훽(휌, 휗, 휁) are minimized at these nodes • The equilibrium solution is still valid everywhere, and spectral collocation theory predicts global convergence • Great flexibility in choosing the nodes • Control grid refinement • Avoid rational surfaces APS DPP November 9, 2020 21 Continuation Methods 1. Perturbations to solve for complex equilibria: • vacuum solution → pressure perturbation → finite-β solution • axisymmetric tokamak → boundary perturbation → 3D stellarator 2. Perturbations to optimize for quasi-symmetry: • axisymmetric tokamak → boundary perturbation → QA stellarator • non-QS equilibrium → perturb some inputs → more-QS equilibrium APS DPP November 9, 2020 22 Order of ODE to Solve Order of Derivatives Variables Equations 0 푅, 푍 1 휕푖푅, 휕푖푍 → 푩 훁 ⋅ 푩 = 0 2 휕푖푗푅, 휕푖푗푍 → 푱 푱 × 푩 = 훁푝 3 휕푖푗푘푅, 휕푖푗푘푍 훁 ⋅ 푱 = 0 • The equilibrium equations are a 2nd-order ODE • Rational surface issues arise at the next higher level with 훁 ⋅ 푱 = 0 APS DPP November 9, 2020 23 Equilibrium Example Inputs Axisymmetric Non-Axisymmetric “D-shaped” Tokamak high-β Heliotron 푅푏 = 3.51 − cos 휃 + 0.106 cos 2휃 푅푏 = 10 − cos 휃 − 0.3 cos 휃 − 19휙 푍푏 = 1.47 sin 휃 + 0.16 sin 2휃 푍푏 = sin 휃 − 0.3 sin 휃 − 19휙 휄 = 1 − 0.67휌2 휄 = 1.5휌2 + 0.5 2 2 푝 = 1.65 × 103 1 − 휌2 푝 = 3.4 × 103 1 − 휌2 휓푎 = 1 휓푎 = 1 APS DPP November 9, 2020 24 Perturbation Example Inputs Axisymmetric Non-Axisymmetric 푀 = 6, 푁 = 2 푅푏 = 5 − cos 휃 푅푏 = 5 − cos 휃 − 0.2 cos 휃 − 휙 푍푏 = sin 휃 푍푏 = sin 휃 − 0.2 sin 휃 − 휙 휄 = 1.618 휄 = 1.618 푝 = 0 푝 = 0 휓푎 = 1 휓푎 = 1 APS DPP November 9, 2020 25.
Recommended publications
  • L-Shell and Energy Dependence of Magnetic Mirror Point of Charged
    Soni et al. Earth, Planets and Space (2020) 72:129 https://doi.org/10.1186/s40623-020-01264-5 FULL PAPER Open Access L-shell and energy dependence of magnetic mirror point of charged particles trapped in Earth’s magnetosphere Pankaj K. Soni* , Bharati Kakad and Amar Kakad Abstract In the Earth’s inner magnetosphere, there exist regions like plasmasphere, ring current, and radiation belts, where the population of charged particles trapped along the magnetic feld lines is more. These particles keep performing gyration, bounce and drift motions until they enter the loss cone and get precipitated to the neutral atmosphere. Theoretically, the mirror point latitude of a particle performing bounce motion is decided only by its equatorial pitch angle. This theoretical manifestation is based on the conservation of the frst adiabatic invariant, which assumes that the magnetic feld varies slowly relative to the gyro-period and gyro-radius. However, the efects of gyro-motion can- not be neglected when gyro-period and gyro-radius are large. In such a scenario, the theoretically estimated mirror point latitudes of electrons are likely to be in agreement with the actual trajectories due to their small gyro-radius. Nevertheless, for protons and other heavier charged particles like oxygen, the gyro-radius is relatively large, and the actual latitude of the mirror point may not be the same as estimated from the theory. In this context, we have carried out test particle simulations and found that the L-shell, energy, and gyro-phase of the particles do afect their mirror points. Our simulations demonstrate that the existing theoretical expression sometimes overestimates or underesti- mates the magnetic mirror point latitude depending on the value of L-shell, energy and gyro-phase due to underlying guiding centre approximation.
    [Show full text]
  • Nicholas Christofilos and the Astron Project in America's Fusion Program
    Elisheva Coleman May 4, 2004 Spring Junior Paper Advisor: Professor Mahoney Greek Fire: Nicholas Christofilos and the Astron Project in America’s Fusion Program This paper represents my own work in accordance with University regulations The author thanks the Program in Plasma Science and Technology and the Princeton Plasma Physics Laboratory for their support. Introduction The second largest building on the Lawrence Livermore National Laboratory’s campus today stands essentially abandoned, used as a warehouse for odds and ends. Concrete, starkly rectangular and nondescript, Building 431 was home for over a decade to the Astron machine, the testing device for a controlled fusion reactor scheme devised by a virtually unknown engineer-turned-physicist named Nicholas C. Christofilos. Building 431 was originally constructed in the late 1940s before the Lawrence laboratory even existed, for the Materials Testing Accelerator (MTA), the first experiment performed at the Livermore site.1 By the time the MTA was retired in 1955, the Livermore lab had grown up around it, a huge, nationally funded institution devoted to four projects: magnetic fusion, diagnostic weapon experiments, the design of thermonuclear weapons, and a basic physics program.2 When the MTA shut down, its building was turned over to the lab’s controlled fusion department. A number of fusion experiments were conducted within its walls, but from the early sixties onward Astron predominated, and in 1968 a major extension was added to the building to accommodate a revamped and enlarged Astron accelerator. As did much material within the national lab infrastructure, the building continued to be recycled. After Astron’s termination in 1973 the extension housed the Experimental Test Accelerator (ETA), a prototype for a huge linear induction accelerator, the type of accelerator first developed for Astron.
    [Show full text]
  • Anl-7807 Anl-7807 Survey of Thermonuclear-Reactor
    ANL-7807 ANL-7807 SURVEY OF THERMONUCLEAR-REACTOR PARAMETERS P. J. Persiani, W. C. Lipinski, and A. J. Hatcli U of C-AUA-USAECB ARGONNE NATIONAL LABORATORY, ARGONNE, ILLINOIS Prepared for the U.S. ATOMIC ENERGY COMMISSION under Contract W-31-109-Eng-38 The faciliUes of Argonne National Laboratory are owned by the "-'f •> S'^^f °°^^%'" ment. Under the terms of a contract (W-31-109-Eng-38) between the U. S. ""^^^^J^^'J/ Commission, Argonne Universities Association and The University of Chicago, the ""'J'"=">' employs the staff and operates the Laboratory in accordance with policies and programs formu- lated, approved and reviewed by the Association. MEMBERS OF ARGONNE UNIVERSITIES ASSOCIATION The Ohio State University The University of Arizona Kansas State University Carnegie-Mellon University The University of Kansas Ohio University Case Western Reserve University Loyola University The Pennsylvania State University The University of Chicago Marquette University Purdue University University o£ Cincinnati Michigan State University Saint Louis University Illinois Institute of Technology The University of Michigan Southern Illinois University University of Illinois University of Minnesota The University of Texas at Austin Indiana University University of Missouri Washington University Iowa State University Northwestern University Wayne State University The University of Iowa University of Notre Dame The University of Wisconsin NOTICE This report was prepared as an account of work sponsored by the United States Government. Neither the United States nor the United States Atomic Energy Commission, nor any of their employees, nor any of their contractors, subcontrac­ tors, or their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness or usefulness of any information, apparatus, product or process disclosed, or represents that its use would not infringe privately-owned rights.
    [Show full text]
  • Effect of Quaslconflned Particles and I= 2 Stellarator Fields on the Negative Mass Lnstablllty in a Modified Betatron G
    TABLE II. Coefficients a and b + r of log(ac,) vs log[2/( p 1/p0 +Pel ACKNOWLEDGMENTS p,)J. We thank Jeanette Nelson for a careful reading of the log(ac1 ) =a+ b log(2/( pc/p1 + p 1/p0 ) l when a = 10 manuscript and G. Pollarolo for useful discussions about the NAG FORTRAN LIBRARY. The numerical computations were 0.1 <Pi <Po< 10 carried out on the Vax 11/780 ofINFN Sezione di Torino. M = 100 M=30 M =IO a; b; r 0.44; 0.06; 0.99 0.33; 0.06; 0.98 0.42; 0.04; 0.98 1 where r0 is the Lorentz factor and spans between 1 and 10. A. E. Gill, Phys. Fluids 8, 1428 (1965). The presence of a shear in cylindrical symmetry has not 2 A. Ferrari, B. Trusooni, and L. Zaninetti, Mon. Not. R. Astron. Soc. 196, yet been considered and will presumably produce a cutoff in !OSI (1981). 3H. Cohn, Astrophys. J. 269, 500 (1983). the instabilities where a> 21'Id, with d (in jet radius units) 4 D. G. Payne and H. Cohn, Astrophys. J. 287, 29S (1984). the length characterizing the thickness of the shear. ' A. M. Anile, J.C. Miller, and S. Motta, Phys. Fluids 26, 14SO (1983). Effect of quaslconflned particles and I= 2 stellarator fields on the negative mass lnstablllty in a modified betatron G. Roberts and N. Rostoker Department ofPhy sics, University ofCalifomia, Irvine. California 92717 (Received 5 February 1985; accepted 10 October 1985) A sufficient stability condition for the negative mass instability is derived.
    [Show full text]
  • Publ. Astron. Obs. Belgrade No. 99 (2020), 256 - 257 Invited Lecture
    Publ. Astron. Obs. Belgrade No. 99 (2020), 256 - 257 Invited Lecture CHALLENGES AND PROGRESS ON THE PATH TOWARDS FUSION ELECTRICITY A. J. H. DONNÉ EUROfusion, Garching, Germany E-mail [email protected] Abstract. The European Roadmap to the Realisation of Fusion Energy1 breaks the quest for fusion energy into eight missions: 1. Plasma regimes of operation: Demonstrate plasma scenarios (based on the tokamak configuration) that increase the success margin of ITER and satisfy the requirements of DEMO. 2. Heat-exhaust systems: Demonstrate an integrated approach that can handle the large power leaving ITER and DEMO plasmas. 3. Neutron tolerant materials: Develop materials that withstand the large 14MeV neutron flux for long periods while retaining adequate physical properties. 4. Tritium self-sufficiency: Find an effective technological solution for the breeding blanket that also drives the generators. 5. Implementation of the intrinsic safety features of fusion: Ensure safety is integral to the design of DEMO using the experience gained with ITER. 6. Integrated DEMO design and system development: Bring together the plasma and all the systems coherently, resolving issues by targeted R&D activities 7. Competitive cost of electricity: Ensure the economic potential of fusion by minimising the DEMO capital and lifetime costs and developing long-term technologies to further reduce power plant costs. 8. Stellarator: Bring the stellarator line to maturity to determine the feasibility of a stellarator power plant. Now we are approaching the end of the 8th European Framework Programme (2014- 2020), it is a good moment to look back at the achievements since the establishment of EUROfusion in 2014, while at the same time have a peek into the future, to see which challenges lay ahead of us as well as the strategy to tackle them.
    [Show full text]
  • Observation of Nuclear Fusion Driven by a Pyroelectric Crystal
    letters to nature 1900þ14. I. An interpretive study of BeppoSAX and Ulysses observations. Astrophys. J. 549, electrostatic field of the crystal is used to generate and accelerate 1021–1038 (2001). a deuteron beam (>100 keV and >4 nA), which, upon striking a 10. Gaensler, B. M. et al. Second-epoch VLA observations of SGR 1806220. GRB Circ. Network 2933 (2005). deuterated target, produces a neutron flux over 400 times the 11. Corbel, S. & Eikenberry, S. S. The connection between W31, SGR 1806220, and LBV 1806220: background level. The presence of neutrons from the reaction Distance, extinction, and structure. Astron. Astrophys. 419, 191–201 (2004). D 1 D ! 3He (820 keV) 1 n (2.45 MeV) within the target is 12. Kolpak, M. A., Jackson, J. M., Bania, T. M. & Dickey, J. M. The radial distribution of cold atomic hydrogen in the galaxy. Astrophys. J. 578, 868–876 (2002). confirmed by pulse shape analysis and proton recoil spectro- 13. Corbel, S. et al. The distance of the soft gamma repeater SGR 1806220. Astrophys. J. 478, 624–630 scopy. As further evidence for this fusion reaction, we use a novel (1997). time-of-flight technique to demonstrate the delayed coincidence 14. Hartmann, D. & Burton, W. B. Atlas of Galactic Neutral Hydrogen. Ch. 4, 169 (Cambridge Univ. Press, between the outgoing a-particle and the neutron. Although the Cambridge, 1997). 15. Garwood, R. W. & Dickey, J. M. Cold atomic gas in the inner Galaxy. Astrophys. J. 338, 841–861 (1989). reported fusion is not useful in the power-producing sense, we 16. Figer, D.
    [Show full text]
  • Astron 104 Laboratory #7 Nuclear Fusion and Stars Chapter 12
    Lab #7 Name: Date: Section: Astron 104 Laboratory #7 Nuclear Fusion and Stars Chapter 12 Introduction You are intimately linked to the nuclear fusion that occurs in stars. Life on Earth thrives on sunlight, which originates as energy produced at the center of the Sun as hydrogen is fused to form helium. Even the very atoms and molecules of your body | iron in your blood, calcium in your bones, phosphorous in your DNA | were forged in stars and dispersed to the interstellar medium in spectacular explosions. In this lab, you will explore the creation of heavy elements and release of energy through nuclear fusion in stars. Learning Objectives At the completion of this lab, you should be able to: 1. Describe the components of atomic nuclei and the properties of subatomic particles 2. Describe the force between two protons as a function of their separation distance 3. Describe the temperature required for nuclear fusion of hydrogen and heavier elements 4. Describe how nuclear fusion produces energy 5. Describe why iron is the heaviest element that can be made via nuclear fusion in stars Atomic Nuclei [5 pts each, 15 pts total] Atomic nuclear are made of combinations of sub-microscopic particles called protons and neutrons. The neutron is slight more massive than the proton, but for now we will ignore that difference and assume the neutron and proton each have a mass of 1 unit. Protons have a positive electrical charge (let's call it +1 unit) and neutrons have no electrical charge. Astron 104 Fall 2015 1 Lab #7 1.
    [Show full text]
  • Towards Nuclear Fusion
    TOWARDS NUCLEAR FUSION The outlook for controlled nuclear fusion has many people had apparently been led to believe. become much more complex than it appeared to be - Professor Teller said that while he believed that at least to the general public - when the possibility of thermonuclear energy generation was possible, it was taming thermonuclear reactions for the production of not going to be "quite easy". Professor L.A.Artsim- useful power was optimistically mentioned at the 1955 ovich of the Soviet Union remarked: "We do not Geneva conference on the peaceful uses of atomic wish to be pessimistic in appraising the future of our energy. Indeed, it may now be impossible to re­ work, yet we must not underestimate the difficulties capture the widespread excitement that followed the which will have to be overcome before we learn to prediction about thermonuclear power by the Presi­ master thermonuclear fusion. " Dr. P. C. Thonemann dent of the conference. Dr. Homi J. Bhabha of India, of the United Kingdom thought that an answer to the and the subsequent disclosures that research on con­ question whether electrical power could be generated trolled nuclear fusion was actively under way in some "using the light elements as fuel by themselves" could of the technically advanced countries. be given only in the next decade, and if the answer was "yes" a further ten years would be required to The significance of this possibility was immedi­ answer the question whether such a power source was ately recognized and there were many enthusiastic economically valuable. accounts of what the generation of power from fusion reactions would mean to the world.
    [Show full text]
  • Introduction to Fusion Energy Plasma Physics (K
    BOOK REVIEWS S',".'"" o book, tor review " based on the editor's 'P'"'""' ",,,d'", possible reader E interest and on the availability of the book to the editor. Occasional selections may include II. L><:5.NB; ns books on topics somewhat peripheral to the subject matter ordinarily considered acceptable. 1111. Controlled Nuclear Fusion - Fundamentals of Its experiments and design studies have evidently so consumed Utilization for Energy Supply the time of experienced workers that there is a great short­ age of rigorous but comprehensible review articles that Authors J. Raeder, K. Borass, R. Bunde, W. should serve as the basis for such a text. Danner, R. Klingelhofer, L. Lengyel, In summary, the book by Raeder et al. should serve as F. Leuterer, M. Soli a useful supplementary text for courses on controlled fusion and a useful enough reference to justify its purchase by Publisher John Wiley & Sons, Inc., Somerset, researchers and instructors active in the various fields of New Jersey (1986) tokamak research that it covers. Despite the recent publica­ tion of a number of very good efforts, the definitive, self­ Pages 316 (illustrated) contained, introductory text on fusion reactor design and a more widely useful reference work for tokamak researchers Price $100.00 remain to be written. Clifford E. Singer Reviewer Clifford E. Singer received his PhD at the University of California, Berkeley. He has worked on the theory and This book is a translation of Kontrol/ierte Kernfusion, applied physics ofplasma transport in tokamak experiments written in 1980. Despite the delay in translation, the book and reactors at Princeton Plasma Physics Laboratory (and remains a timely summary of many aspects of tokamak the University ofIllinois) since 1977.
    [Show full text]
  • Insight Manufacturers, Publishers and Suppliers by Product Category
    Manufacturers, Publishers and Suppliers by Product Category 2/15/2021 10/100 Hubs & Switch ASANTE TECHNOLOGIES CHECKPOINT SYSTEMS, INC. DYNEX PRODUCTS HAWKING TECHNOLOGY MILESTONE SYSTEMS A/S ASUS CIENA EATON HEWLETT PACKARD ENTERPRISE 1VISION SOFTWARE ATEN TECHNOLOGY CISCO PRESS EDGECORE HIKVISION DIGITAL TECHNOLOGY CO. LT 3COM ATLAS SOUND CISCO SYSTEMS EDGEWATER NETWORKS INC Hirschmann 4XEM CORP. ATLONA CITRIX EDIMAX HITACHI AB DISTRIBUTING AUDIOCODES, INC. CLEAR CUBE EKTRON HITACHI DATA SYSTEMS ABLENET INC AUDIOVOX CNET TECHNOLOGY EMTEC HOWARD MEDICAL ACCELL AUTOMAP CODE GREEN NETWORKS ENDACE USA HP ACCELLION AUTOMATION INTEGRATED LLC CODI INC ENET COMPONENTS HP INC ACTI CORPORATION AVAGOTECH TECHNOLOGIES COMMAND COMMUNICATIONS ENET SOLUTIONS INC HYPERCOM ADAPTEC AVAYA COMMUNICATION DEVICES INC. ENGENIUS IBM ADC TELECOMMUNICATIONS AVOCENT‐EMERSON COMNET ENTERASYS NETWORKS IMC NETWORKS ADDERTECHNOLOGY AXIOM MEMORY COMPREHENSIVE CABLE EQUINOX SYSTEMS IMS‐DELL ADDON NETWORKS AXIS COMMUNICATIONS COMPU‐CALL, INC ETHERWAN INFOCUS ADDON STORE AZIO CORPORATION COMPUTER EXCHANGE LTD EVGA.COM INGRAM BOOKS ADESSO B & B ELECTRONICS COMPUTERLINKS EXABLAZE INGRAM MICRO ADTRAN B&H PHOTO‐VIDEO COMTROL EXACQ TECHNOLOGIES INC INNOVATIVE ELECTRONIC DESIGNS ADVANTECH AUTOMATION CORP. BASF CONNECTGEAR EXTREME NETWORKS INOGENI ADVANTECH CO LTD BELDEN CONNECTPRO EXTRON INSIGHT AEROHIVE NETWORKS BELKIN COMPONENTS COOLGEAR F5 NETWORKS INSIGNIA ALCATEL BEMATECH CP TECHNOLOGIES FIRESCOPE INTEL ALCATEL LUCENT BENFEI CRADLEPOINT, INC. FORCE10 NETWORKS, INC INTELIX
    [Show full text]
  • Advanced Spheromak Fusion Reactor
    UCRL-JC-124363 Advanced Spheromak Fusion Reactor T. Kenneth Fowler University of California at Berkeley '4 6 E. Bickford Hooper Lawrence Livermore National Laboratory This paper was prepared for submittal to 8th International Conference on Emerging Nuclear Energy Systems Obninsk, Russia June 24-28,1996 June 19,1996 This is a preprint of a paper intended for publication in a journal or proceeding Since changes may be made before publication, this preprint is made available with the understanding that it will not be cited or reproduced.withoutthe permission of the author. DfSTRIBUlION OF TffFS DOCUMENT UNLrMlTEO IS 4? DISCLAIMER This document was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the University of California nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information. apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark. manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or the University of California. and shall not be. used for advertising or product endorsement purposes. DISCLAIMER Portions of this document may be illegible in electronic image products. Images are produced from the best available original document.
    [Show full text]
  • Generation of Astron-Spherator Configuration
    INSTITUTE OF PLASMA PHYSICS NAGOYA UNIVERSITY Generation of Astron-Spherator Configuration Kazumichi NARIHARA, Mitsuru HASEGAWA Yukihiro TOMITA, Tetsuya TSUZUKI Kuninori SATO and Akihiro MOHRI (Received - Dec.27, 1982 ) IPPJ-621 Jan. 1983 RESEARCH REPORT NAGOYA, JAPAN Generation of Astron-Spherator Configuration Kazumichi NARIHARA, Mitsuru HASEGAWA Yukihiro TOMITA, Tetsuya TSUZUKI Kuninori SATO and Akihiro MOHRI (Received - Dec.27, 1982 ) IPPJ-621 Jan. 1983 Further communication about this report is to be sent to the Research Information Center, Institute of Plasma Physics, Nagoya University, Nagoya 464, Japan. Generation of Astron-Spherator Configuration issue ) Kazumiohi NARIHARA, Mitsuru HASEGAWA, Yukihiro TORIITA, Tetsuya TSUZUKI, Kuninori SATO and Akihiro MOHRI Institute of Plasma Physics, Nagoya University Nagoya, 4 64 Synopsis It was experimentally demonstrated that Astron-Spherator configuration is formed by injecting a pulsed relativistic electron bean in a toroidal device SPAC-VI with external toroidal and vertical magnetic fields. A plasma is confined in the extended magnetic region produced by a slender core of current carrying energetic electrons. This configuration continued for 40 ms without fatal instabilities. - 2 - Introduct ion A toroidal magnetic system called 'Astron-Spherator' was first considered for plasma confinement by S.Yoshikawa and N.C.Christfilos so a- to synthesise the favourable features of 2) 3 4) Astron and Spherator * . Spherator is equipped with a conducting hard core carrying toroidal currents. Its magnetic configuration is characterised by deep magnetic well, high magnetic shear and short connection length between the bad and good curvature regions. However, the presence of the hard core,floated and immersed inside the confined plasma, makes it difficult to extend this system to a fusion reactor concept.
    [Show full text]