Overview of FRC-related modeling (July 2014-present)

Size: px
Start display at page:

Download "Overview of FRC-related modeling (July 2014-present)"

Transcription

1 Overview of FRC-related modeling (July 2014-present) Artan Qerushi AFRL-UCLA Basic Research Collaboration Workshop January 20th, 2015 AFTC PA Release# 15009, 16 Jan 2015 Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

2 Outline 1 FRC configuration Illustration of FRC configuration FRC formation with RMF 2 2D (r θ) model of RMF-formed FRCs Original publications of RMF-formed FRCs Model equations and their numerical solution 3 Collision Radiative model 0D model equations Test calculations 4 Magnetized plasma closure for electron-ion-neutral mixture Magnetized plasma closure for electron-ion-neutral mixture Applied field modules 5 Conclusions Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

3 Illustration of FRC configuration θ pinch formed FRC. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

4 FRC formation with RMF (Rotating Magnetic Field) RMF-formed FRC. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

5 Original publications of modeling RMF-formed FRCs W. N. Hugrass and R. C. Grimm, A numerical study of the generation of an azimuthal current ina plasma cylinder using a transverse rotating magnetic field, J. Plasma Physics 26, (1981). R. D. Milroy, A numerical study of rotating magnetic fields as a current drive for field reversed configurations, Phys. Plasmas 6, (1999). R. D. Milroy, A magnetohydrodynamic model of rotating magnetic field current drive in a field-reversed configuration, Phys. Plasmas 7, (2000). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

6 Illustration of 2D (r θ) RMF model Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

7 2D (r θ) model equations and their numerical solution Coupled equations for B z (r, θ) and A z (r, θ) A z t B z t = η µ 0 2 A z + = η µ 0 2 B z + 1 n e eµ 0 r 1 n e eµ 0 r [( ) ( Az Bz r θ [( Az θ ) ) r 2 A z ( ) ( )] Az Bz θ r ( ) ] Az r θ 2 A z Solve numerically by expanding B z and A z in Fourier series and deriving 1D equations for the Fourier coefficients (coupled). Advance the equations for the Fourier coefficients in time using an ODE solver (method of lines). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

8 0D model equations dn e dt dn k dt dn z dt dt e dt = q i n i k=1,z = n e (n k 1 S k 1 n k (S k + α k ) + n k+1 α k+1 ), k = 1, 2,, Z 1 = n e (n Z 1 S Z 1 n Z α Z ), = n k L k (T e ) k=0,z All coefficients S k, α k, L k are functions of T e. NIST (FLYCHK code) provides these coefficients for all the periodic table in the temperature range 0.5 [ev] to 100 [kev]. Have implemented in C++ the coefficients for H, He, O, Ar, Kr and Xe. Little work needed to add other elements if that is needed. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

9 Test calculation for Oxygen R. A. Hulse, Numerical studies of impurities in fusion plasmas, Nuclear Technol/Fusion 3, (1983). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

10 Test calculation for Oxygen (time-dependent radiation loss) R. A. Hulse, Numerical studies of impurities in fusion plasmas, Nuclear Technol/Fusion 3, (1983). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

11 Magnetized plasma closure (Introduction) (ρ j v j ) t = Πj + k R jk ɛ e t = q j Πj: v j + k Q jk indicate terms not related to collisions, R jk represents transfer of momentum between the species j and k due to collisions, q j represents flux of heat due to the temperature gradient T j and Ohmic heating, Π represents the off-diagonal part of the pressure tensor, Q jk represents collisional energy transfer between the species j and k. This term is proportional to the temperature difference T k T j. The purpose of a fluid closure is to provide explicit expressions for R jk, q j, Π and Q jk in terms of n j, v j, T j and the magnetic field strength. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

12 Magnetized plasma closure (literature) N. A. Bobrova et. al, Magnetohydrodynamic two-temperature equations for multicomponent plasma, Phys. Plasmas 12, (2005). E. T. Meier and U. Shumlak, A general nonlinear fluid model for reacting plasma-neutral mixture, Phys. Plasmas 19, (2012). E. M. Epperlein and M. G. Heines, Plasma transport coefficients in a magnetic field by direct numerical simulation of the Fokker-Planck equation, Phys. Fluids 29, (1986). A. Decoster, Fluid equations and transport coefficients in plasmas, in Modeling Collisions, Elsevier Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

13 Magnetized plasma closure (results) E. M. Epperlein and M. G. Haines, Plasma transport coefficients in a magnetic field by direct numerical solution of the Fokker-Planck equation, Phys. Fluids 29, (1986). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

14 Conclusions We have implemented in stand-alone C++ code The CR data for H, He, O, Ar, Kr and Xe. Magnetized plasma closure for electron-ion-neutral mixture. Applied field modules for the FRC experiment setup including the DC coils and the RMF antenna. The 2D r θ model is currently implemented in modern fortran. A C++ version which uses the SUNDIALS ODE solver suite is being written to be included in our software framework (SMURF). Work on multi-dimensional fluid models using the Finite Element method and unstructured grids is ongoing and will be reported in the near future. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

15 Extra Slides. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

16 Assumptions of 2D RMF model Frequency condition for current drive ω ci < ω < ω ce. Infinitely long plasma cyclinder lying in a uniform magnetic field. All quantities are assumed independent of axial location z. The ions form a uniformly distributed neutralizing background of fixed, massive positive charges. The plasma resistivity, η, is taken to be a scalar quantity which is constant in time and uniform on space. In particular, η is assumed to be of the form η = m eν ei n e e 2 Electron inertia is neglected; that is it is assumed that ω ω ce, ν ei. The displacement current is neglected. This implies that only systems for which ωr p /c, where r p is the plasma radius, are considered. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

17 FRC formation with RMF (Rotating Magnetic Field) Illustration of 2D (r θ) RMF formation FRC. Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

18 Test calculation for Iron R. A. Hulse, Numerical studies of impurities in fusion plasmas, Nuclear Technol/Fusion 3, (1983). Artan Qerushi (AFRL) FRC modeling overview (July 2014-present) January 20th, / 18

Fluid equations, magnetohydrodynamics

Fluid equations, magnetohydrodynamics Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics

More information

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University

More information

MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION

MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends

More information

Simple examples of MHD equilibria

Simple examples of MHD equilibria Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will

More information

Progress in the Plasma Science and Innovation Center

Progress in the Plasma Science and Innovation Center 1 THP/2-02 Progress in the Plasma Science and Innovation Center U. Shumlak 1, C. Akcay 1, A. H. Glasser 1, C. J. Hansen 1, E. D. Held 2, T. R. Jarboe 1, J.-Y. Ji 2, C. Kim 1, W. Lowrie 1, V. S. Lukin 4,

More information

The Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is

The Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is and Stability of Field-Reversed Equilibrium with Toroidal Field Configurations Atomics General Box 85608, San Diego, California 92186-5608 P.O. APS Annual APS Meeting of the Division of Plasma Physics

More information

Heating and current drive: Radio Frequency

Heating and current drive: Radio Frequency Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)

More information

R. Clark, D. Humbert, K. Sheikh Nuclear Data Section

R. Clark, D. Humbert, K. Sheikh Nuclear Data Section Calculation of Atomic Data for Plasma Modeling: Introduction and Atomic Structure Part 1 R. Clark, D. Humbert, K. Sheikh Nuclear Data Section Overview Plasmas in fusion research Data needs for plasma modeling

More information

Self-Organization of Plasmas with Flows

Self-Organization of Plasmas with Flows Self-Organization of Plasmas with Flows ICNSP 2003/ 9/10 Graduate School of Frontier Sciences,, National Institute for Fusion Science R. NUMATA, Z. YOSHIDA, T. HAYASHI ICNSP 2003/ 9/10 p.1/14 Abstract

More information

Weibel instability and filamentary structures of a relativistic electron beam in plasma

Weibel instability and filamentary structures of a relativistic electron beam in plasma Mitglied der Helmholtz-Gemeinschaft 7 th Direct Drive and Fast Ignition Workshop Prague, 3-6 May, 009 Weibel instability and filamentary structures of a relativistic electron beam in plasma Anupam Karmakar

More information

Monte Carlo Collisions in Particle in Cell simulations

Monte Carlo Collisions in Particle in Cell simulations Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature

More information

Waves in plasma. Denis Gialis

Waves in plasma. Denis Gialis Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.

More information

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models 0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,

More information

A half submerged metal sphere (UIC comprehensive

A half submerged metal sphere (UIC comprehensive Problem 1. exam) A half submerged metal sphere (UIC comprehensive A very light neutral hollow metal spherical shell of mass m and radius a is slightly submerged by a distance b a below the surface of a

More information

Neoclassical transport

Neoclassical transport Neoclassical transport Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 28 th January 2013 Dr Ben Dudson Magnetic Confinement Fusion (1 of 19) Last time Toroidal devices

More information

THE SHEARED-FLOW STABILIZED Z-PINCH

THE SHEARED-FLOW STABILIZED Z-PINCH THE SHEARED-FLOW STABILIZED Z-PINCH U. Shumlak* J. Chadney R.P. Golingo D.J. Den Hartog M.C. Hughes S.D. Knecht W. Lowrie V.S. Lukin B.A. Nelson R.J. Oberto J.L. Rohrbach M.P. Ross and G.V. Vogman *University

More information

Multi-fluid Simulation Models for Inductively Coupled Plasma Sources

Multi-fluid Simulation Models for Inductively Coupled Plasma Sources Multi-fluid Simulation Models for Inductively Coupled Plasma Sources Madhusudhan Kundrapu, Seth A. Veitzer, Peter H. Stoltz, Kristian R.C. Beckwith Tech-X Corporation, Boulder, CO, USA and Jonathan Smith

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem .65, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Formulation of the Stability Problem Hierarchy of Formulations of the MHD Stability Problem for Arbitrary 3-D Systems. Linearized equations

More information

Progress in Vlasov-Fokker- Planck simulations of laserplasma

Progress in Vlasov-Fokker- Planck simulations of laserplasma Progress in Vlasov-Fokker- Planck simulations of laserplasma interactions C. P. Ridgers, M. W. Sherlock, R. J. Kingham, A.Thomas, R. Evans Imperial College London Outline Part 1 simulations of long-pulse

More information

NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE

NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Collaborators:

More information

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,

More information

Overview of Accelerated Simulation Methods for Plasma Kinetics

Overview of Accelerated Simulation Methods for Plasma Kinetics Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math

More information

Modeling of ELM Dynamics for ITER

Modeling of ELM Dynamics for ITER Modeling of ELM Dynamics for ITER A.Y. PANKIN 1, G. BATEMAN 1, D.P. BRENNAN 2, A.H. KRITZ 1, S. KRUGER 3, P.B. SNYDER 4 and the NIMROD team 1 Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015

More information

0 Magnetically Confined Plasma

0 Magnetically Confined Plasma 0 Magnetically Confined Plasma 0.1 Particle Motion in Prescribed Fields The equation of motion for species s (= e, i) is written as d v ( s m s dt = q s E + vs B). The motion in a constant magnetic field

More information

Fluid Equations for Rarefied Gases

Fluid Equations for Rarefied Gases 1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel

More information

APPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics

APPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)

More information

The Virial Theorem, MHD Equilibria, and Force-Free Fields

The Virial Theorem, MHD Equilibria, and Force-Free Fields The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely

More information

MAGNETOHYDRODYNAMICS

MAGNETOHYDRODYNAMICS Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study

More information

TURBULENT TRANSPORT THEORY

TURBULENT TRANSPORT THEORY ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical

More information

Physics of High Pressure Helicon Plasma and Effect of Wavenumber Spectrum

Physics of High Pressure Helicon Plasma and Effect of Wavenumber Spectrum (1)Title Physics of High Pressure Helicon Plasma and Effect of Wavenumber Spectrum Interdisciplinary Graduate School of Engineering Sciences, Kyushu Univeristy, Japan Shunjiro SHINOHARA Scientific Center

More information

AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code

AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code AMSC 663 Project Proposal: Upgrade to the GSP Gyrokinetic Code George Wilkie (gwilkie@umd.edu) Supervisor: William Dorland (bdorland@umd.edu) October 11, 2011 Abstract Simulations of turbulent plasma in

More information

A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model

A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model A Study of 3-Dimensional Plasma Configurations using the Two-Fluid Plasma Model B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program University of Washington IEEE International Conference

More information

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX

More information

Module I: Electromagnetic waves

Module I: Electromagnetic waves Module I: Electromagnetic waves Lecture 9: EM radiation Amol Dighe Outline 1 Electric and magnetic fields: radiation components 2 Energy carried by radiation 3 Radiation from antennas Coming up... 1 Electric

More information

PLASMA-NEUTRAL MODELING IN NIMROD. Uri Shumlak*, Sina Taheri*, Jacob King** *University of Washington **Tech-X Corporation April 2016

PLASMA-NEUTRAL MODELING IN NIMROD. Uri Shumlak*, Sina Taheri*, Jacob King** *University of Washington **Tech-X Corporation April 2016 PLASMA-NEUTRAL MODELING IN NIMROD Uri Shumlak*, Sina Taheri*, Jacob King** *University of Washington **Tech-X Corporation April 2016 Plasma-Neutral Model Physical Model is derived by E. Meier and U. Shumlak*

More information

On existence of resistive magnetohydrodynamic equilibria

On existence of resistive magnetohydrodynamic equilibria arxiv:physics/0503077v1 [physics.plasm-ph] 9 Mar 2005 On existence of resistive magnetohydrodynamic equilibria H. Tasso, G. N. Throumoulopoulos Max-Planck-Institut für Plasmaphysik Euratom Association

More information

Integrated Modeling of Fast Ignition Experiments

Integrated Modeling of Fast Ignition Experiments Integrated Modeling of Fast Ignition Experiments Presented to: 9th International Fast Ignition Workshop Cambridge, MA November 3-5, 2006 R. P. J. Town AX-Division Lawrence Livermore National Laboratory

More information

Generating of fusion plasma neutron source with AFSI for Serpent MC neutronics computing Serpent UGM 2015 Knoxville, TN,

Generating of fusion plasma neutron source with AFSI for Serpent MC neutronics computing Serpent UGM 2015 Knoxville, TN, Generating of fusion plasma neutron source with AFSI for Serpent MC neutronics computing Serpent UGM 2015 Knoxville, TN, 14.10.2015 Paula Sirén VTT Technical Research Centre of Finland, P.O Box 1000, 02044

More information

EUV spectra from the NIST EBIT

EUV spectra from the NIST EBIT EUV spectra from the NIST EBIT D. Kilbane and G. O Sullivan Atomic and Molecular Plasma Physics group, UCD, Ireland J. D. Gillaspy, Yu. Ralchenko and J. Reader National Institute of Standards and Technology,

More information

Estimating the plasma flow in a recombining plasma from

Estimating the plasma flow in a recombining plasma from Paper P3-38 Estimating the plasma flow in a recombining plasma from the H α emission U. Wenzel a, M. Goto b a Max-Planck-Institut für Plasmaphysik (IPP) b National Institute for Fusion Science, Toki 509-5292,

More information

Low-field helicon discharges

Low-field helicon discharges Plasma Phys. Control. Fusion 39 (1997) A411 A420. Printed in the UK PII: S0741-3335(97)80958-X Low-field helicon discharges F F Chen, X Jiang, J D Evans, G Tynan and D Arnush University of California,

More information

A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere

A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere P. Song Center for Atmospheric Research University of Massachusetts Lowell V. M. Vasyliūnas Max-Planck-Institut

More information

Multi-Scale Simulations for Fast Ignition. H. Nagatomo, T. Johzaki, A. Sunahara, and K. Mima Institute of Laser Engineering g Osaka University

Multi-Scale Simulations for Fast Ignition. H. Nagatomo, T. Johzaki, A. Sunahara, and K. Mima Institute of Laser Engineering g Osaka University Multi-Scale Simulations for Fast Ignition and Related Laser Plasma Physics H. Nagatomo, T. Johzaki, A. Sunahara, and K. Mima Institute of Laser Engineering g Osaka University Y. Sentoku University of Nevada

More information

Theory and Simulation of Neoclassical Transport Processes, with Local Trapping

Theory and Simulation of Neoclassical Transport Processes, with Local Trapping Theory and Simulation of Neoclassical Transport Processes, with Local Trapping Daniel H. E. Dubin Department of Physics, University of California at San Diego, La Jolla, CA USA 92093-0319 Abstract. Neoclassical

More information

Lectures on basic plasma physics: Introduction

Lectures on basic plasma physics: Introduction Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters

More information

THE TOROIDAL HELICON EXPERIMENT AT IPR

THE TOROIDAL HELICON EXPERIMENT AT IPR THE TOROIDAL HELICON EXPERIMENT AT IPR Manash Kumar Paul Supervisor Prof. Dhiraj Bora D.D.G, ITER Acknowledgement Dr. S.K.P.Tripathi UCLA Institute for Plasma Research Gandhinagar, Gujarat (India) MOTIVATION

More information

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学, RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation

More information

Collisions and transport phenomena

Collisions and transport phenomena Collisions and transport phenomena Collisions in partly and fully ionized plasmas Typical collision parameters Conductivity and transport coefficients Conductivity tensor Formation of the ionosphere and

More information

Fundamentals of Plasma Physics

Fundamentals of Plasma Physics Fundamentals of Plasma Physics Definition of Plasma: A gas with an ionized fraction (n i + + e ). Depending on density, E and B fields, there can be many regimes. Collisions and the Mean Free Path (mfp)

More information

Propagation of Radio Frequency Waves Through Density Filaments

Propagation of Radio Frequency Waves Through Density Filaments PSFC/JA-15-13 Propagation of Radio Frequency Waves Through Density Filaments A. K. Ram and K. Hizanidis a May 015 a National Technical University of Athens (part of HELLAS) School of Electrical and Computer

More information

INTRODUCTION TO FLUID MECHANICS June 27, 2013

INTRODUCTION TO FLUID MECHANICS June 27, 2013 INTRODUCTION TO FLUID MECHANICS June 27, 2013 PROBLEM 3 (1 hour) A perfect liquid of constant density ρ and constant viscosity µ fills the space between two infinite parallel walls separated by a distance

More information

Plasma Simulation in STAR-CCM+ -- towards a modern software tool

Plasma Simulation in STAR-CCM+ -- towards a modern software tool Plasma Simulation in STAR-CCM+ -- towards a modern software tool 1 Introduction Electrical arcs and other forms of industrial plasmas have a large number of technical applications, including circuit breakers,

More information

Feature-level Compensation & Control

Feature-level Compensation & Control Feature-level Compensation & Control 2 Plasma Eray Aydil, UCSB, Mike Lieberman, UCB and David Graves UCB Workshop November 19, 2003 Berkeley, CA 3 Feature Profile Evolution Simulation Eray S. Aydil University

More information

The importance of including XMHD physics in HED codes

The importance of including XMHD physics in HED codes The importance of including XMHD physics in HED codes Charles E. Seyler, Laboratory of Plasma Studies, School of Electrical and Computer Engineering, Cornell University Collaborators: Nat Hamlin (Cornell)

More information

Predicting the Rotation Profile in ITER

Predicting the Rotation Profile in ITER Predicting the Rotation Profile in ITER by C. Chrystal1 in collaboration with B. A. Grierson2, S. R. Haskey2, A. C. Sontag3, M. W. Shafer3, F. M. Poli2, and J. S. degrassie1 1General Atomics 2Princeton

More information

Jet Stability: A computational survey

Jet Stability: A computational survey Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,

More information

14. Energy transport.

14. Energy transport. Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation

More information

PHYSICS OF HOT DENSE PLASMAS

PHYSICS OF HOT DENSE PLASMAS Chapter 6 PHYSICS OF HOT DENSE PLASMAS 10 26 10 24 Solar Center Electron density (e/cm 3 ) 10 22 10 20 10 18 10 16 10 14 10 12 High pressure arcs Chromosphere Discharge plasmas Solar interior Nd (nω) laserproduced

More information

Sheared Flow Stabilization in the Z-Pinch

Sheared Flow Stabilization in the Z-Pinch 1 IF/P7-28 Sheared Flow Stabilization in the Z-Pinch U. Shumlak, C.S. Adams, J.M. Blakely, B.J. Chan, R.P. Golingo, S.D. Knecht, B.A. Nelson, R.J. Oberto, M.R. Sybouts, and G.V. Vogman Aerospace & Energetics

More information

A cylinder in a magnetic field (Jackson)

A cylinder in a magnetic field (Jackson) Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles

More information

Real Plasma with n, T ~ p Equilibrium: p = j B

Real Plasma with n, T ~ p Equilibrium: p = j B Real Plasma with n, T ~ p Equilibrium: p = j B B lines must lie in isobaric surfaces. Since B = 0, only possible if isobaric surfaces are topological tori. Magnetic field lines must form nested tori. Equilibrium

More information

August 2013 Qualifying Exam. Part II

August 2013 Qualifying Exam. Part II August 2013 Qualifying Exam Part II Mathematical tables are allowed. Formula sheets are provided. Calculators are allowed. Please clearly mark the problems you have solved and want to be graded. Do only

More information

Plasma Spectroscopy Inferences from Line Emission

Plasma Spectroscopy Inferences from Line Emission Plasma Spectroscopy Inferences from Line Emission Ø From line λ, can determine element, ionization state, and energy levels involved Ø From line shape, can determine bulk and thermal velocity and often

More information

Plasma instabilities. Dr Ben Dudson, University of York 1 / 37

Plasma instabilities. Dr Ben Dudson, University of York 1 / 37 Plasma instabilities Dr Ben Dudson, University of York 1 / 37 Previously... Plasma configurations and equilibrium Linear machines, and Stellarators Ideal MHD and the Grad-Shafranov equation Collisional

More information

Fluid Dynamics and Balance Equations for Reacting Flows

Fluid Dynamics and Balance Equations for Reacting Flows Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and

More information

Studies of Lower Hybrid Range of Frequencies Actuators in the ARC Device

Studies of Lower Hybrid Range of Frequencies Actuators in the ARC Device Studies of Lower Hybrid Range of Frequencies Actuators in the ARC Device P. T. Bonoli, Y. Lin. S. Shiraiwa, G. M. Wallace, J. C. Wright, and S. J. Wukitch MIT PSFC, Cambridge, MA 02139 59th Annual Meeting

More information

Fokker-Planck Simulation Study of Hot-Tail Effect on Runaway Electron Generation in ITER Disruptions )

Fokker-Planck Simulation Study of Hot-Tail Effect on Runaway Electron Generation in ITER Disruptions ) Fokker-Planck Simulation Study of Hot-Tail Effect on Runaway Electron Generation in ITER Disruptions ) Hideo NUGA, Akinobu MATSUYAMA 1), Masatoshi YAGI 1) and Atsushi FUKUYAMA Kyoto University, Kyoto-Daigaku-Katsura,

More information

MHD Linear Stability Analysis Using a Full Wave Code

MHD Linear Stability Analysis Using a Full Wave Code US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto

More information

Low Temperature Plasma Technology Laboratory

Low Temperature Plasma Technology Laboratory Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California

More information

Plasma instability during ITBs formation with pellet injection in tokamak

Plasma instability during ITBs formation with pellet injection in tokamak Plasma instability during ITBs formation with pellet injection in tokamak P. Klaywittaphat 1, B. Chatthong 2, T. Onjun. R. Picha 3, J. Promping 3 1 Faculty of Engineering, Thaksin University, Phatthalung,

More information

Predictive Engineering and Computational Sciences. Detailed and simplified kinetic mechanisms for high enthalpy flows

Predictive Engineering and Computational Sciences. Detailed and simplified kinetic mechanisms for high enthalpy flows Marco Panesi Non-equilibriumst ionization January 31 st, 211 1 / 22 PECOS Predictive Engineering and Computational Sciences Detailed and simplified kinetic mechanisms for high enthalpy flows Application

More information

Theoretical Analysis of Formation and Sustainment Methods for Compact Toroids

Theoretical Analysis of Formation and Sustainment Methods for Compact Toroids Brazilian Journal of Physics, vol. 34, no. 4B, Decemer, 24 1621 Theoretical Analysis of Formation and Sustainment Methods for Compact Toroids Ricardo Farengo, Agustín F. Lifschitz, Hugo E. Ferrari, Seastián

More information

Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect )

Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Yusuke SHIMIZU, Takaaki FUJITA, Atsushi OKAMOTO, Hideki ARIMOTO, Nobuhiko HAYASHI 1), Kazuo HOSHINO 2), Tomohide

More information

Linear and Nonlinear Dust Acoustic Waves, Shocks and Stationary Structures in DC-Glow-Discharge Dusty Plasma Experiments.

Linear and Nonlinear Dust Acoustic Waves, Shocks and Stationary Structures in DC-Glow-Discharge Dusty Plasma Experiments. 53rd Annual Meeting of the APS Division of Plasma Physics BI2.00005 Monday November 14, 2011 Linear and Nonlinear Dust Acoustic Waves, Shocks and Stationary Structures in DC-Glow-Discharge Dusty Plasma

More information

RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE

RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE Abhay K. Ram Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139. U.S.A. Joan Decker

More information

Equilibrium and transport in Tokamaks

Equilibrium and transport in Tokamaks Equilibrium and transport in Tokamaks Jacques Blum Laboratoire J.-A. Dieudonné, Université de Nice Sophia-Antipolis Parc Valrose 06108 Nice Cedex 02, France jblum@unice.fr 08 septembre 2008 Jacques Blum

More information

Fluid Equations for Rarefied Gases

Fluid Equations for Rarefied Gases 1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel

More information

26. Non-linear effects in plasma

26. Non-linear effects in plasma Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless

More information

Modelling of plasma edge turbulence with neutrals

Modelling of plasma edge turbulence with neutrals Modelling of plasma edge turbulence with neutrals Ben Dudson 1 1 York Plasma Institute, Department of Physics, University of York, Heslington, York YO1 5DD, UK 7 th IAEA TM on Plasma Instabilities 4-6

More information

Fundamentals of Plasma Physics Transport in weakly ionized plasmas

Fundamentals of Plasma Physics Transport in weakly ionized plasmas Fundamentals of Plasma Physics Transport in weakly ionized plasmas APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Luís L Alves (based on Vasco Guerra s original slides) 1 As perguntas

More information

Effect of Ion Orbit Loss on Rotation and the Radial Electric Field in the DIII-D Tokamak

Effect of Ion Orbit Loss on Rotation and the Radial Electric Field in the DIII-D Tokamak Effect of Ion Orbit Loss on Rotation and the Radial Electric Field in the DIII-D Tokamak by T.M. Wilks 1 with W.M. Stacey 1 and T.E. Evans 2 1 Georgia Institute of Technology 2 General Atomics Presented

More information

Non-ohmic ignition scenarios in Ignitor

Non-ohmic ignition scenarios in Ignitor Non-ohmic ignition scenarios in Ignitor Augusta Airoldi IFP, EURATOM-ENEA-CNR Association, Milano, Italy Francesca Bombarda, Giovanna Cenacchi Ignitor Group, ENEA, Italy Bruno Coppi MIT, USA DPP1 APS Meeting

More information

MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS*

MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS* MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS* Kapil Rajaraman** and Mark J. Kushner*** **Department of Physics ***Department of Electrical and

More information

IMPURITY ANALYSIS AND MODELING OF DIII-D RADIATIVE MANTLE DISCHARGES

IMPURITY ANALYSIS AND MODELING OF DIII-D RADIATIVE MANTLE DISCHARGES IMPURITY ANALYSIS AND MODELING OF DIII-D RADIATIVE MANTLE DISCHARGES J. Mandrekas, W.M. Stacey Georgia Institute of Technology M. Murakami, M.R. Wade ORNL G. L. Jackson General Atomics Presented at the

More information

Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating

Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating Tests of Profile Stiffness Using Modulated Electron Cyclotron Heating by T.C. Luce in collaboration with J.C. DeBoo, C.C. Petty, J. Pino, J.M. Nelson, and J.C.M. dehaas Presented at 9th EU-US Transport

More information

Electromagnetics in COMSOL Multiphysics is extended by add-on Modules

Electromagnetics in COMSOL Multiphysics is extended by add-on Modules AC/DC Module Electromagnetics in COMSOL Multiphysics is extended by add-on Modules 1) Start Here 2) Add Modules based upon your needs 3) Additional Modules extend the physics you can address 4) Interface

More information

Modelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves

Modelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves Modelling of low-temperature plasmas: kinetic and transport mechanisms L.L. Alves llalves@tecnico.ulisboa.pt Instituto de Plasmas e Fusão Nuclear Instituto Superior Técnico, Universidade de Lisboa Lisboa,

More information

Mathematical Modelling of a Magnetized Target Fusion Reactor

Mathematical Modelling of a Magnetized Target Fusion Reactor Mathematical Modelling of a Magnetized Target Fusion Reactor Michael PIC Assistant Adjunct Professor, Mathematics UCLA April 20, 2016 Context A Canadian research company General Fusion (2002-Present, Michel

More information

INFLUENCE OF THERMODIFFUSIVE PARTICLE TRANSPORT ON THERMOMAGNETIC CONVECTION IN MAGNETIC FLUIDS

INFLUENCE OF THERMODIFFUSIVE PARTICLE TRANSPORT ON THERMOMAGNETIC CONVECTION IN MAGNETIC FLUIDS MAGNETOHYDRODYNAMICS Vol. 49 (2013), No. 3-4, pp. 473 478 INFLUENCE OF THERMODIFFUSIVE PARTICLE TRANSPORT ON THERMOMAGNETIC CONVECTION IN MAGNETIC FLUIDS TU Dresden, Chair of Magnetofluiddynamics, Measuring

More information

Application Of Faraday s Law

Application Of Faraday s Law Application Of Faraday s Law Dr Miguel Cavero September 2, 2014 Application Of Faraday s Law September 2, 2014 1 / 23 The PHYS120 Exam will be divided into three sections as follows: Section A: Short Questions

More information

Application of Rarefied Flow & Plasma Simulation Software

Application of Rarefied Flow & Plasma Simulation Software 2016/5/18 Application of Rarefied Flow & Plasma Simulation Software Yokohama City in Japan Profile of Wave Front Co., Ltd. Name : Wave Front Co., Ltd. Incorporation : March 1990 Head Office : Yokohama

More information

Nautilus: Robust, positivity and divergence preserving code for multi-fluid, multi-species electromagnetics and plasma physics applications

Nautilus: Robust, positivity and divergence preserving code for multi-fluid, multi-species electromagnetics and plasma physics applications Nautilus: Robust, positivity and divergence preserving code for multi-fluid, multi-species electromagnetics and plasma physics applications Ammar H. Hakim and John Loverich I. Background In this paper

More information

Chapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current

Chapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with

More information

12. MHD Approximation.

12. MHD Approximation. Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal

More information

Critical edge gradients and flows with reversed magnetic field in Alcator C-Mod

Critical edge gradients and flows with reversed magnetic field in Alcator C-Mod Critical edge gradients and flows with reversed magnetic field in B. LaBombard, N. Smick, M. Greenwald, J.W. Hughes B. Lipschultz, K. Marr, J.L. Terry, Team Contributed talk JO1.00004 Presented at the

More information

Numerical Modeling of Radiative Kinetic Plasmas

Numerical Modeling of Radiative Kinetic Plasmas 2014 US-Japan JIFT Workshop on Progress in kinetic plasma simulations Oct.31-Nov.1, 2014, Salon F, New Orleans Marriott, New Orleans, LA, U.S.A Numerical Modeling of Radiative Kinetic Plasmas T. Johzaki

More information

2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson

2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson 2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves

More information

Experimental evaluation of nonlinear collision effect on the beam slowing-down process

Experimental evaluation of nonlinear collision effect on the beam slowing-down process P-2 Experimental evaluation of nonlinear collision effect on the beam slowing-down process H. Nuga R. Seki,2 S. Kamio M. Osakabe,2 M. Yokoyama,2 M. Isobe,2 K. Ogawa,2 National Institute for Fusion Science,

More information

Alex M Imai, Y. Ohta and A. Itoh Department of Nuclear Engineering, Kyoto University

Alex M Imai, Y. Ohta and A. Itoh Department of Nuclear Engineering, Kyoto University Alex M Imai, Y. Ohta and A. Itoh Department of Nuclear Engineering, Kyoto University Joint IAEA-NFRI Technical Meeting on Data Evaluation for Atomic, Molecular and Plasma-Material Interaction Processes

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information