Magnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape.
|
|
- Ella Carpenter
- 5 years ago
- Views:
Transcription
1 1 Magnetized ion collection by oblique surfaces including self-consistent drifts: Mach-probes of arbitrary shape I H Hutchinson Plasma Science and Fusion Center and and Engineering Department MIT APS DPP Dallas, Nov 2008
2 Abstract 2 A complete analytic theory for magnetized Mach-probes, when cross-field diffusion is neglected, is presented It is shown that the full self-consistent quasi-neutral fluid drift equations around an ion-collecting probe of arbitrary 3-D shape, in a magnetized isothermal plasma with background parallel and perpendicular flow, can be solved exactly The resulting flux to the probe (per unit area perpendicular to B) is n c s exp( 1 M + M cot θ), where θ is the angle between the surface and B in the plane of background-drift This exponential dependence is in good agreement with prior numerical fits of the diffusive case The (background) perpendicular Mach number, M, is that arising from the sum of ExB and, counter-intuitively, electron (not ion) diamagnetic drifts Fluid displacements in the magnetic presheath are important, and included in this expression, but give rise to small additional terms at some orientations Temperature-gradient diamagnetic drifts can be added, but only approximately: both electron and ion drifts contribute Corrections of order Larmor radius divided by electrode-dimensions are also evaluated They can bias the results for small probes
3 Probes in Tokamaks Used for many important studies of the plasma edge 3 Plasma flows, parallel and perpendicular to B, are important and can presumably be measured Represent generic problem of absorbing object in a flowing magnetized plasma My initial motivation for this longstanding work Layout of scanning probes in Alcator C-Mod
4 Magnetized Plasma Flow to Objects 4 Surface is a sink of ions and electrons Anisotropic collection with crossfield diffusion Ion acceleration along field Electron repulsion Self-consistent potential variation Local sheath structure is multi-layered Ions enter magnetic presheath at v = c s Enter Debye sheath at vk = c s v Magnetic Presheath B k θ Plasma: Quasineutral Drifts Thicknesses: Magnetic Presheath ρ S Solid Surface Debye Sheath Debye Sheath 4λ D
5 External Pre-existing Plasma Flow 5 Flow past the object is some combination of parallel and perpendicular ion drifts v v B Large flows in tokamak scrape-offs Plasma flow Object motion Plasma Flow v Object Space shuttle surface-charging Mach Probe Idea: Measure flow by ion flux density Electrodes facing in different directions should presumably collect different ion flux densities, depending on the external plasma flow θ Collectors B Insulator
6 Practical Mach probes are complicated 6 MacLatchy et al 1992 Gundestrup Peleman et al 2006 (Tokamak) Antoni et al 1996 (RFP) Smick and LaBombard 2008
7 Physics applies to some astrophysical situations 7 Interaction of Jupiter s moons with its magnetosphere is in the same regime: Larmor radius smaller than object (moon) smaller than scale-length
8 What s new about the treatment? 8 Old: Heuristic Diffusion v = D ln n + v Elliptic B v v External n = n, M = M New: Rigorous Convection: v (x) solved Hyperbolic v v v B Instead of invoking heuristic cross-field diffusion, ignore diffusion Treat pure imposed transverse-drift, solve self-consistently, no heuristics More satisfactory physics Rigorous mathematics Analytic solution! Result consistent Dependence on diamagnetic drifts discovered
9 Ion Fluid Drift Equations 9 (nv) = 0 mn(v )v = nze φ p + nze(v B) Uniform B Neglect friction, vicosity Zero diffusion Cross-field drift from (momentum equation) B replaces perpendicular momentum: [ ( Keep = Drift Approx v = φ + 1 ) Drop ] nze m p + Ze (v )v B }{{}}{{} B 2 Electric + Diamagnetic Polarization Drop the polarization drift term in the plasma region (not magnetic presheath) Use simplest possible closure: Uniform temperature T i =const (v ) ln n + v = 0 (v )v + Ze ( m φ + T ) i Ze ln n = 0 Electron parallel momentum φ = (T e /e) ln n e, with uniform T e (for now) Is used to eliminate potential φ by quasi-neutrality (Boltzmann electrons)
10 Hyperbolic system Analyze by Characteristics 10 Define c 2 s = (ZT e + T i )/m and M = v/c s : Mach Number Eliminate φ Then equations become M ln n + M = 0 M M + ln n = 0 which can be rearranged to display explicitly the characteristics (M + )(ln n + M ) = 0 (M )(ln n M ) = 0 Thus the quantities (ln n ± M ) are constant along their corresponding characteristics: dx = (M ± B/B)ds Integration along characteristics is the secret to solving hyperbolic systems
11 Solution by Characteristics 1 11 Assume M simply uniform (in y-direction) Take B along x 2-D problem Analyse the higher-x side of the object Two characteristics pass through any point Object Region P 3 P 4 θ 4 θ P 0 2 P P 1 Positive Characteristic θ m 2 y θ Perturbed Plasma Region M Unperturbed Plasma Region x B Negative Characteristic Negative characteristics start at y in unperturbed plasma: ln n M = ln n M If positive characteristic also starts at y then also: ln n + M = ln n + M Only solution then: n = n, M = M Point (P 0 ) in unperturbed plasma Unperturbed characteristics are straight lines with slope tanθ m = M /(M ± 1) If line from point to lowest part of object is shallower than M /(M + 1), then positive characteristic does not intersect object and point is unperturbed
12 2 The perturbed positive characteristics determine the solution 12 P 4 y M If positive characteristic intersects θ 4 P 2 θ x B object, ln n + M is still constant, but ln n + M Consequently Object Region P 3 θ P 0 2 P 1 Positive Characteristic θ m Perturbed Plasma Region Negative Characteristic Unperturbed Plasma Region M = const, ln n = const along straight +ve characteristics, on which M = M cot θ 1 Boundary condition is unconstrained outflow at object surface: ln n = ln n M + M M is most negative, consistent with overall solution θ is as large as possible This means the characteristic is tangent to surface eg P 1, P 2 Points such as P 3 in concave regions have a non-local tangency (M greater) Angles greater than π/2 are perfectly permissible (eg P 4 )
13 Uniform Velocity Solution* 13 For homogeneous perpendicular velocity M M h, M = M h cot θ 1, n = n exp( 1 M + M h cot θ) throughout the plasma region, where θ is the angle (to B) of the positive characteristic, which is the tangent of greatest possible angle At the plasma boundary with object, the outgoing flux per unit surface area is nc s (M h cos θ M sin θ) = nc s sin θ which is most conveniently written as flux per unit perpendicular area as Γ = nc s = n c s exp( 1 M + M h cot θ) with θ the local surface tangent (in convex regions) This is fine, but who says M is uniform? Generally it is not because of the perturbation to n and φ in the presheath For a 2-D probe (z an ignorable coordinate) the extra drifts are along z: no effect But for a 3-D object what happens? *I H Hutchinson Phys Rev Lett 101, (2008)
14 Now Comes the Magic! Consider a 3-D situation with the full, self-consistent transverse velocity Integrate the electron equation along the field to obtain φ φ (z) = (T e /e) ln(n/n ) Substitute into the drift equation to get transverse velocity in two parts uniform self consistent external v /c s = M = mc 2 s φ + }{{} Ze ln n }{{} M h M 1 Integrating the negative characteristic equation: B c s B 2 = M h + M 1 14 ln n = ln n + M M implies that n = n(m ), so M n Therefore a perpendicular velocity M 1 satisfying M 1 ln n = 0 (and ours does) gives zero contribution to the characteristics: ([M + M h + M 1 ] ± )(ln n ± M ) }{{} original characteristics = ([M + M h ] ± )(ln n ± M ) }{{} modified characteristics We can solve the problem using characterstics that contain just the homogeneous external transverse drift The uniform drift solution just discussed is correct for arbitrary 3-D probe shape = 0
15 Example 3-D Solution: Sphere 15 Contours at 01 spacing of M + 1 = 1 + M + ln(n/n ) Object: unity-radius sphere External perpendicular flow: M h = 024 External parallel flow: M = 02
16 Diamagnetic drift along contour lines 16 Sphere M = 0 M h = 03 M + 1 contours Coincide with ln n contours
17 Example 3-D Solution: Pyramid Probe 17 Contours at 01 spacing of M + 1 = 1 + M + ln(n/n ) Object: C-Mod-like pyramid probe External perpendicular flow: M h = 015 External parallel flow: M = 01
18 n External Diamagnetic Drift Suppose ln n (z) = ẑ/l n Integration along the negative characteristic gives 18 ln n M = ln(n (z )) M where z is the starting z-value of the characteristic The z-motion along the characteristic can t be ignored a priori ln(n/n (z)) = M + M δz/l n where δz z z is not constant But I ll show δz = δz(m ) In which case ln n/n is also function only of M, so write [ ( v /c s = M = φ + T ) ] i Ze ln n + mc2 s Ze ln(n/n ) }{{}}{{} M M h 1 B c s B 2, Write the characteristic equations to take advantage of this form d ( ) ln n ± M = d ( ) ln n/n ± M + Mz /L n = 0 ± ± c s dt c s dt where d c s dt M ± = (M ± 1) ± x + M y y + M z z, when acting on functions of ln n/n, can be d = (M ± 1) ± x + M h y c s dt
19 δz(m k ) analytical solution 19 Define d dm (δz/2l n ) r; so that d δz c s dt ± 2L n = c d s dt M d δz ± dm 2L n = r d c s dt M ± Eliminate M -derivatives from characteristics, eg d c s dt (2M δz + L n ) + M z L n = 0 to get d δz c s dt 2L n = M z 2L n and d δz c s dt + 2L n = M z r 2L n 1 r Then eliminate partial x-derivatives of ln n/n between the positive and negative characteristics, and partial y-derivatives using the z-velocity expression ln M z = ρ n ln n/n s y = ρ s y Substitute for d c s dt δz, and divide through by M z to arrive at the quadratic equation: ± [ ] 1 = M r +1 M r r (M 1) + 4L n M h /ρ s Solve, using the simplifying notation: u 0 2L n M h /ρ s, u M + u 0, to find r = (1/4)[u + 2 ± u 2 + 4u 0 ] Integrate to obtain the z-displacement δz = 2 rdm = 1 [ { u 2 + 4u u u L n u u 0 ln (±u + )} ] u u 2 + 4u 0 /2 u This solution is real only if 1/u o > 1/4, ie the ion diamagnetic drift magnitude M ni must be less than the E B drift magnitude M E
20 δz value proves unimportant I ve proved that the solution exists; that s critical for the presumptions Solution shows the δz values are small: 20 The values are nearly symmetric about M = 1: nearly even parity in M + 1 Perpendicular Mach number is derived from the odd-parity in M + 1 part of flux For practically all parameters, δz/l n is negligible < 2% contribution to M
21 Positive characteristics now curved 21 Schematic shapes of positive characteristics in x-y plane Curvature exaggerated Object Object θ m M y θ m M y 1 + M 1 + M L n positive M ni opposes M E L n negative M ni reinforces M E Notice that the curvature changes the definition of convex Plane surfaces are concave when the characteristics curve towards them We don t actually have to solve for the characteristic trajectories
22 Magnetic Presheath Displacement 22 Summary so far: the flux density from the plasma into the magnetic presheath has been shown, for an arbitrary (convex) 3-D object shape, to be with M h = M E + M ni and Γ = nc s = n (z )c s exp[ 1 (M M h cot θ)] n (z ) = n (z δz) = n (z) exp( δz/l n ) Although I ve solved for it, plasma displacement δz/l n can generally be ignored When n 0 (there s diamagnetic drift) the corresponding displacement in the MPS is not ignorable Must be calculated Consider general surface orientation normal: ˆk sin α = ˆkB/B other directions ˆl B ˆk, ĵ Integrate the ˆk-component of the full momentum equation including inertia to find the total displacement in the ˆl direction through the MPS After significant work, obtain δ l = ρ s [1 sin α] 2 sin α cos α
23 Displacement changes flux density 23 The displacement δ lˆl has a component along the gradient ln n The impact point samples flux entering MPS a distance δ lˆlẑ up the gradient z E 1 End B E E 2 y E B E B This has odd-parity (unlike δz) Effectively substracts a important fraction of M D from measurement x E δ lˆl θ 1 θ 2 δ lˆl E [Also causes potential end-effect] Stalk
24 Temperature Gradients 24 Accommodate the temperature-gradient part of Diamagnetic drifts by assuming T e and T i are functions only of z Two major changes occur First, c s (z) is spatially varying Characteristics become (v ± c s )[ln n/n ± v /c s ] = ±(v /c s )( v z /L c ) v z /L p where L c c s /dc s /dz is scale-length of c s -variation, and L p of ion-pressure Look for solutions that are of the form ln n/n M = g(m ) Reduce to a quadratic in r dg/2dm using same techniques as before (M 1) ( M L c + 1 L p ) ( r 1 + r + (M M + 1) L c + 1 L p ) = 2M h ρ s (1 + 1/2r) For this to be a consistent solution we require w M h L c /ρ s to be constant Solve for r Integrate to find g = 2 r dm
25 Characteristic-variation small 25 Solutions for g = ln n/n M Generally small and/or even-parity provided we avoid 0 < w < 3 [Mathematical Note: All this is necessary because these corrections are not small by ordering They are the same order in ρ s /L as M So we have to evaluate them before we can say if they are important] Incidentally, the geometric effect of g-variation is to make the positive characteristic curved The smallness of g amounts to saying the radius of curvature is bigger than the object
26 New temperature-gradient term looks like electron diamagnetic drift 26 Second T-effect New type of drift term [ ( ) ( Te M = ln( n n ) + φ + T ) i + T i e Ze Ze ln n + mc2 s Ze ln( n n ) }{{}}{{} M M h 1 Combines parallel density difference with perp T e -gradient The M h = ln(n/n )M Te + M Ti + M E + M ni is no longer homogeneous But, with boundary condition, M = M h cot θ 1 gives ln n/n = 1 + g + (M Ti + M E + M ni ) cot θ 1 + M Te cot θ ] B c s B 2 Adding the temperature gradient effects into the MPS analysis results in a similar term but full pressure-gradient effect [Skipping non-trivial algebra on MPS effects]
27 Final Full Flux Density Approximating g = M and writing to first order in M : 27 ln { Γ p nc s } = 1 M + [ (1 + M )M Te + M Di + M E ( ) ] 1 sin α M 1 + sin α D cot θ with M E = E B/c s B 2 ŷ Electric field drift, agreeing with intuition M Di = p B/(c s nzeb 2 )ŷ Total ion diamagnetic drift Intuitive M Te = T e B/(c s eb 2 )ŷ Electron diamagnetic due to T e gradient M D = M Di M De Difference between ion and electron diamag drifts α Angle between B and surface θ Angle in x-y-plane between B and surface The quantity in [ ] is what a transverse Mach probe measures Cf prior diffusive treatment ln { Γ p nc s } = 1 11(M M cot θ) * I H Hutchinson Phys Rev A, 37 (1988) 4358; IHH Principles of Plasma Diagnostics (2002)
28 Discussion 28 Solutions are exact Equations are approximate but good approximations Drift approximation probably better than diffusive for tokamaks Geometry is arbitrary (provided convex) Generality a gift from God! Results agree (within 10%) with standard diffusive M E formula which is based on ad hoc fit to numeric solutions But now given rigorous analytic justification Previously, the dependence on Diamagnetic drifts was simply unknown { } [ ( ) ] Γ p 1 sin α ln = 1 M + (1 + M )M nc Te + M Di + M E M s 1 + sin α D cot θ We now see that there is an extra apparent velocity in the measurements that is equal to M Te This may help to make sense of the Smick/LaBombard observations A more complete formula and fuller understanding of the effect of MPS displacement M D has been obtained The drift approximation is problematic when M D 1 because that requires ρ s /L 1 A local approximation is then impossible [Why solution by expansion impossible]
29 Summary 29 A complete analytic solution of the magnetized plasma drift equations surrounding an arbitrary 3-D absorbing object has been obtained It shows that the (appropriately normalized) flux-density depends only on surface orientation, not position on the object, for convex regions It shows why the (known) Mach probe calibration has the form it does, and what the (previously-unknown) diamagnetic drift dependence is
30 Supplemental 30
31 Eddy size equals or exceeds probe 31 C-Mod scanning probe head (schematic) compared in scale with Example of imaged turbulence eddies in SOL Probe body is comparable in size Probe electrodes are smaller
Oblique ion collection in the drift-approximation: how magnetized Mach-probes really work.
arxiv:0809.1557v2 [physics.plasm-ph] 28 Oct 2008 Oblique ion collection in the drift-approximation: how magnetized Mach-probes really work. I H Hutchinson Plasma Science and Fusion Center and Department
More informationImprovement of a 1-D Fluid Probe Model for Mach Probe Measurements
Improvement of a -D Fluid Probe Model for Mach Probe Measurements P. Peleman, S. Jachmich, M. Van Schoor, G. Van Oost Ghent University Department of Applied Physics, B-9000 Ghent, Belgium Laboratory for
More informationThe Invalidity of a Mach Probe Model
The Invalidity of a Mach Probe Model I. H. Hutchinson Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA, USA Abstract Despite its recent application to the interpretation
More informationFluid 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 informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationHelium-3 transport experiments in the scrape-off layer with the Alcator C-Mod omegatron ion mass spectrometer
PHYSICS OF PLASMAS VOLUME 7, NUMBER 11 NOVEMBER 2000 Helium-3 transport experiments in the scrape-off layer with the Alcator C-Mod omegatron ion mass spectrometer R. Nachtrieb a) Lutron Electronics Co.,
More informationTotal Flow Vector in the C-Mod SOL
Total Flow Vector in the SOL N. Smick, B. LaBombard MIT Plasma Science and Fusion Center APS-DPP Annual Meeting Atlanta, GA November 3, 2009 Motivation and Goals Measurements have revealed high parallel
More informationDrift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer
Drift-Driven and Transport-Driven Plasma Flow Components in the Alcator C-Mod Boundary Layer N. Smick, B. LaBombard MIT Plasma Science and Fusion Center PSI-19 San Diego, CA May 25, 2010 Boundary flows
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationOperational Phase Space of the Edge Plasma in Alcator C-Mod
Operational Phase Space of the Edge Plasma in B. LaBombard, T. Biewer, M. Greenwald, J.W. Hughes B. Lipschultz, N. Smick, J.L. Terry, Team Contributed talk RO.00008 Presented at the 47th Annual Meeting
More informationKinetic solution to the Mach probe problem in transversely flowing strongly magnetized plasmas. L. Patacchini I. H. Hutchinson
PSFC/JA-9-24 Kinetic solution to the Mach probe problem in transversely flowing strongly magnetized plasmas L. Patacchini I. H. Hutchinson October 29 Plasma Science and Fusion Center Massachusetts Institute
More informationBlob sizes and velocities in the Alcator C-Mod scrapeoff
P1-59 Blob sizes and velocities in the Alcator C-Mod scrapeoff layer R. Kube a,b,*, O. E. Garcia a,b, B. LaBombard b, J. L. Terry b, S. J. Zweben c a Department of Physics and Technology, University of
More informationImpact of neutral atoms on plasma turbulence in the tokamak edge region
Impact of neutral atoms on plasma turbulence in the tokamak edge region C. Wersal P. Ricci, F.D. Halpern, R. Jorge, J. Morales, P. Paruta, F. Riva Theory of Fusion Plasmas Joint Varenna-Lausanne International
More informationTURBULENT 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 informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationFlow measurements in the Scrape-Off Layer of Alcator C-Mod using Impurity Plumes
Flow measurements in the Scrape-Off Layer of Alcator C-Mod using Impurity Plumes S. Gangadhara,. Laombard M.I.T. Plasma Science and Fusion Center, 175 Albany St., Cambridge, MA 2139 USA Abstract Accurate
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationSpin Stability of Aysmmetrically Charged Plasma Dust. I.H. Hutchinson. February 2004
PSFC/JA-04-3 Spin Stability of Aysmmetrically Charged Plasma Dust I.H. Hutchinson February 2004 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139 USA This work
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationIntrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas. Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I.
Intrinsic rotation due to non- Maxwellian equilibria in tokamak plasmas Jungpyo (J.P.) Lee (Part 1) Michael Barnes (Part 2) Felix I. Parra MIT Plasma Science & Fusion Center. 1 Outlines Introduction to
More informationConnections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod
Connections between Particle Transport and Turbulence Structures in the Edge and SOL of Alcator C-Mod I. Cziegler J.L. Terry, B. LaBombard, J.W. Hughes MIT - Plasma Science and Fusion Center th 19 Plasma
More informationInnovative Concepts Workshop Austin, Texas February 13-15, 2006
Don Spong Oak Ridge National Laboratory Acknowledgements: Jeff Harris, Hideo Sugama, Shin Nishimura, Andrew Ware, Steve Hirshman, Wayne Houlberg, Jim Lyon Innovative Concepts Workshop Austin, Texas February
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationSpherical probes at ion saturation in ExB fields. L. Patacchini I. H. Hutchinson
PSFC/JA-9-5 Spherical probes at ion saturation in ExB fields L. Patacchini I. H. Hutchinson October 9 Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge MA 39 USA Submitted
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationMAGNETOHYDRODYNAMICS
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 informationMHD 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 informationStudying the Formation of the Pre-Sheath in an Oblique Magnetic Field using a Fluid Model and PIC Simulation
J. Plasma Fusion Res. SERIES, Vol. 8 (2009) Studying the Formation of the Pre-Sheath in an Oblique Magnetic Field using a Fluid Model and PIC Simulation Jernej KOVAI 1, Tomaž Gyergyek 2,1, Milan EREK 1,3
More informationKinematics
Classical Mechanics Kinematics Velocities Idea 1 Choose the appropriate frame of reference. Useful frames: - Bodies are at rest - Projections of velocities vanish - Symmetric motion Change back to the
More informationA note on the plasma sheath and the Bohm Criterion
A note on the plasma sheath and the Bohm Criterion G.D. Severn Dept. of Physics, University of San Diego, San Diego CA 92110 (Dated: April 6, 2006) PACS numbers: 52.27.Aj, 52.27.Cm The word sheath in connection
More informationTokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011
Tokamak Divertor System Concept and the Design for ITER Chris Stoafer April 14, 2011 Presentation Overview Divertor concept and purpose Divertor physics General design considerations Overview of ITER divertor
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More informationTurbulence and flow in the Large Plasma Device
Turbulence and flow in the Large Plasma Device D.A. Schaffner, T.A. Carter, P. Popovich, B. Friedman Dept of Physics, UCLA Gyrokinetics in Laboratory and Astrophysical Plasmas Isaac Newton Institute of
More informationGyrokinetic simulations including the centrifugal force in a strongly rotating tokamak plasma
Gyrokinetic simulations including the centrifugal force in a strongly rotating tokamak plasma F.J. Casson, A.G. Peeters, Y. Camenen, W.A. Hornsby, A.P. Snodin, D. Strintzi, G.Szepesi CCFE Turbsim, July
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationSingle Particle Motion
Single Particle Motion C ontents Uniform E and B E = - guiding centers Definition of guiding center E gravitation Non Uniform B 'grad B' drift, B B Curvature drift Grad -B drift, B B invariance of µ. Magnetic
More informationDirect drive by cyclotron heating can explain spontaneous rotation in tokamaks
Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks J. W. Van Dam and L.-J. Zheng Institute for Fusion Studies University of Texas at Austin 12th US-EU Transport Task Force Annual
More informationAstronomy, Astrophysics, and Cosmology
Astronomy, Astrophysics, and Cosmology Luis A. Anchordoqui Department of Physics and Astronomy Lehman College, City University of New York Lesson VI March 15, 2016 arxiv:0706.1988 L. A. Anchordoqui (CUNY)
More informationDynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit. Transport Physics of the Density Limit
Dynamics of Zonal Shear Collapse in Hydrodynamic Electron Limit Transport Physics of the Density Limit R. Hajjar, P. H. Diamond, M. Malkov This research was supported by the U.S. Department of Energy,
More informationImpact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma
1 Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma D. Galassi, P. Tamain, H. Bufferand, C. Baudoin, G. Ciraolo, N. Fedorczak, Ph.
More informationPhysics of fusion power. Lecture 13 : Diffusion equation / transport
Physics of fusion power Lecture 13 : Diffusion equation / transport Many body problem The plasma has some 10 22 particles. No description is possible that allows for the determination of position and velocity
More informationDiffusion equation, flux, diffusion coefficient scaling. Diffusion in fully ionized plasma vs. weakly ionized plasma. n => Coulomb collision frequency
Last Time Diffusion in plasma: the main reason why we need to control it (i.e. using magnetic field) Diffusion equation, flux, diffusion coefficient scaling o o t nx,t Dn D2 nx,t o D ~ L 2 T Diffusion
More informationComparison of Pellet Injection Measurements with a Pellet Cloud Drift Model on the DIII-D Tokamak
Comparison of Pellet Injection Measurements with a Pellet Cloud Drift Model on the DIII-D Tokamak T.C. Jernigan, L.R. Baylor, S.K. Combs, W.A. Houlberg (Oak Ridge National Laboratory) P.B. Parks (General
More informationShock and Expansion Waves
Chapter For the solution of the Euler equations to represent adequately a given large-reynolds-number flow, we need to consider in general the existence of discontinuity surfaces, across which the fluid
More informationElectron temperature barriers in the RFX-mod experiment
Electron temperature barriers in the RFX-mod experiment A. Scaggion Consorzio RFX, Padova, Italy Tuesday 5 th October 2010 ADVANCED PHYSICS LESSONS 27/09/2010 07/10/2010 IPP GARCHING JOINT EUROPEAN RESEARCH
More informationStudies 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 informationMAGNETIC 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 informationPlasma Science and Fusion Center
Plasma Science and Fusion Center Turbulence and transport studies in ALCATOR C Mod using Phase Contrast Imaging (PCI) Diagnos@cs and Comparison with TRANSP and Nonlinear Global GYRO Miklos Porkolab (in
More informationL-H transitions driven by ion heating in scrape-off layer turbulence (SOLT) model simulations
L-H transitions driven by ion heating in scrape-off layer turbulence (SOLT) model simulations D.A. Russell, D.A. D Ippolito and J.R. Myra Research Corporation, Boulder, CO, USA Presented at the 015 Joint
More informationThe Levitated Dipole Experiment: Towards Fusion Without Tritium
The Levitated Dipole Experiment: Towards Fusion Without Tritium Jay Kesner MIT M.S. Davis, J.E. Ellsworth, D.T. Garnier, M.E. Mauel, P.C. Michael, P.P. Woskov MCP I3.110 Presented at the EPS Meeting, Dublin,
More informationParticle Transport Measurements in the LHD Stochastic Magnetic Boundary Plasma using Mach Probes and Ion Sensitive Probe
Particle Transport Measurements in the LHD Stochastic Magnetic Boundary Plasma using Mach Probes and Ion Sensitive Probe N. Ezumi a*, K. Todoroki a, T. Kobayashi b, K. Sawada c, N. Ohno b, M. Kobayashi
More informationNeoclassical 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 informationEffect of divertor nitrogen seeding on the power exhaust channel width in Alcator C-Mod
Effect of divertor nitrogen seeding on the power exhaust channel width in Alcator C-Mod B. LaBombard, D. Brunner, A.Q. Kuang, W. McCarthy, J.L. Terry and the Alcator Team Presented at the International
More informationAlcator C-Mod. Particle Transport in the Scrape-off Layer and Relationship to Discharge Density Limit in Alcator C-Mod
Alcator C-Mod Particle Transport in the Scrape-off Layer and Relationship to Discharge Density Limit in Alcator C-Mod B. LaBombard, R.L. Boivin, M. Greenwald, J. Hughes, B. Lipschultz, D. Mossessian, C.S.
More informationTransport and drift-driven plasma flow components in the Alcator C-Mod boundary plasma
Transport and drift-driven plasma flow components in the Alcator C-Mod boundary plasma The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationWaves 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 informationFundamentals 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 informationPLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION
PLANAR KINETICS OF A RIGID BODY FORCE AND ACCELERATION I. Moment of Inertia: Since a body has a definite size and shape, an applied nonconcurrent force system may cause the body to both translate and rotate.
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationVerification & Validation: application to the TORPEX basic plasma physics experiment
Verification & Validation: application to the TORPEX basic plasma physics experiment Paolo Ricci F. Avino, A. Bovet, A. Fasoli, I. Furno, S. Jolliet, F. Halpern, J. Loizu, A. Mosetto, F. Riva, C. Theiler,
More informationMark Howell Gonzaga High School, Washington, D.C.
Be Prepared for the Calculus Exam Mark Howell Gonzaga High School, Washington, D.C. Martha Montgomery Fremont City Schools, Fremont, Ohio Practice exam contributors: Benita Albert Oak Ridge High School,
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory Equilibrium theory for plasma discharges of finite length Francis F. Chen and Davide Curreli LTP-6 June, Electrical Engineering Department Los Angeles, California
More informationCross-Field Plasma Transport and Main Chamber Recycling in Diverted Plasmas on Alcator C-Mod
Cross-Field Plasma Transport and Main Chamber Recycling in Diverted Plasmas on Alcator C-Mod B. LaBombard, M. Umansky, R.L. Boivin, J.A. Goetz, J. Hughes, B. Lipschultz, D. Mossessian, C.S. Pitcher, J.L.Terry,
More informationMicrotearing Simulations in the Madison Symmetric Torus
Microtearing Simulations in the Madison Symmetric Torus D. Carmody, P.W. Terry, M.J. Pueschel - University of Wisconsin - Madison dcarmody@wisc.edu APS DPP 22 Overview PPCD discharges in MST have lower
More informationLow 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 informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
More informationEffect of secondary beam neutrals on MSE: theory
Effect of secondary beam neutrals on MSE: theory S. Scott (PPPL) J. Ko, I. Hutchinson (PSFC/MIT) H. Yuh (Nova Photonics) Poster NP8.87 49 th Annual Meeting, DPP-APS Orlando, FL November 27 Abstract A standard
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationNon-perturbative statistical theory of intermittency in ITG drift wave turbulence with zonal flows
Non-perturbative statistical theory of intermittency in ITG drift wave turbulence with zonal flows Johan Anderson 1 and Eun-jin Kim University of Sheffield Department of Applied Mathematics Hicks Building,
More informationPhysics of MR Image Acquisition
Physics of MR Image Acquisition HST-583, Fall 2002 Review: -MRI: Overview - MRI: Spatial Encoding MRI Contrast: Basic sequences - Gradient Echo - Spin Echo - Inversion Recovery : Functional Magnetic Resonance
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationReminders: Show your work! As appropriate, include references on your submitted version. Write legibly!
Phys 782 - Computer Simulation of Plasmas Homework # 4 (Project #1) Due Wednesday, October 22, 2014 Reminders: Show your work! As appropriate, include references on your submitted version. Write legibly!
More informationMeasurements of relativistic emission from runaway electrons in Alcator C-Mod: spectrum, polarization, and spatial structure
Measurements of relativistic emission from runaway electrons in Alcator C-Mod: spectrum, polarization, and spatial structure R. Granetz, R. Mumgaard MIT PSFC APS-DPP New Orleans 2014/10/30 Motivation Modeling
More information196 7 atmospheric oscillations:
196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere
More informationChapter 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 informationFundamentals 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 informationIssues of Perpendicular Conductivity and Electric Fields in Fusion Devices
Issues of Perpendicular Conductivity and Electric Fields in Fusion Devices Michael Tendler, Alfven Laboratory, Royal Institute of Technology, Stockholm, Sweden Plasma Turbulence Turbulence can be regarded
More information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More informationModélisation de sources plasma froid magnétisé
Modélisation de sources plasma froid magnétisé Gerjan Hagelaar Groupe de Recherche Energétique, Plasma & Hors Equilibre (GREPHE) Laboratoire Plasma et Conversion d Énergie (LAPLACE) Université Paul Sabatier,
More informationA neoclassical model for toroidal rotation and the radial electric field in the edge pedestal. W. M. Stacey
A neoclassical model for toroidal rotation and the radial electric field in the edge pedestal W. M. Stacey Fusion Research Center Georgia Institute of Technology Atlanta, GA 30332, USA October, 2003 ABSTRACT
More informationOn the physics of shear flows in 3D geometry
On the physics of shear flows in 3D geometry C. Hidalgo and M.A. Pedrosa Laboratorio Nacional de Fusión, EURATOM-CIEMAT, Madrid, Spain Recent experiments have shown the importance of multi-scale (long-range)
More informationSimple 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 information34.3. Resisted Motion. Introduction. Prerequisites. Learning Outcomes
Resisted Motion 34.3 Introduction This Section returns to the simple models of projectiles considered in Section 34.1. It explores the magnitude of air resistance effects and the effects of including simple
More informationCompound Perpendicular Diffusion of Cosmic Rays and Field Line Random Walk, with Drift
Compound Perpendicular Diffusion of Cosmic Rays and Field Line Random Walk, with Drift G.M. Webb, J. A. le Roux, G. P. Zank, E. Kh. Kaghashvili and G. Li Institute of Geophysics and Planetary Physics,
More informationFlow and dynamo measurements in the HIST double pulsing CHI experiment
Innovative Confinement Concepts (ICC) & US-Japan Compact Torus (CT) Plasma Workshop August 16-19, 211, Seattle, Washington HIST Flow and dynamo measurements in the HIST double pulsing CHI experiment M.
More informationSMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.
SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction
More informationA kinetic neutral atom model for tokamak scrape-off layer tubulence simulations. Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva
A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva CRPP - EPFL SPS Annual Meeting 2014 02.07.2014 CRPP The tokamak
More informationLecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.
Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationModelling 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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationTH/P4-9. T. Takizuka 1), K. Shimizu 1), N. Hayashi 1), M. Hosokawa 2), M. Yagi 3)
1 Two-dimensional Full Particle Simulation of the Flow Patterns in the Scrape-off-layer Plasma for Upper- and Lower- Null Point Divertor Configurations in Tokamaks T. Takizuka 1), K. Shimizu 1), N. Hayashi
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More informationAnalysis of Jeans Instability of Partially-Ionized. Molecular Cloud under Influence of Radiative. Effect and Electron Inertia
Adv. Studies Theor. Phys., Vol. 5, 2011, no. 16, 755-764 Analysis of Jeans Instability of Partially-Ionized Molecular Cloud under Influence of Radiative Effect and Electron Inertia B. K. Dangarh Department
More informationDr. Allen Back. Sep. 10, 2014
Dr. Allen Back Sep. 10, 2014 The chain rule in multivariable calculus is in some ways very simple. But it can lead to extremely intricate sorts of relationships (try thermodynamics in physical chemistry...
More informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
More informationPLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT. Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford
1 PLASMA: WHAT IT IS, HOW TO MAKE IT AND HOW TO HOLD IT Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford 2 Overview Why do we need plasmas? For fusion, among other things
More informationStudy of Enhanced D α H-modes Using the Alcator C-Mod Reflectometer
Study of Enhanced D α H-modes Using the Reflectometer Y. Lin 1, J.H. Irby, E.S. Marmar, R. Nazikian, M. Greenwald, A.E. Hubbard, J. Hughes, I.H. Hutchinson, B. LaBombard, A. Mazurenko, E. Nelson-Melby,
More informationTH/P8-4 Second Ballooning Stability Effect on H-mode Pedestal Scalings
TH/P8-4 Second Ballooning Stability Effect on H-mode Pedestal Scalings T. Onjun 1), A.H. Kritz ), G. Bateman ), A. Pankin ) 1) Sirindhorn International Institute of Technology, Klong Luang, Pathumthani,
More informationCollision Processes. n n The solution is 0 exp x/ mfp
Collision Processes Collisions mediate the transfer of energy and momentum between various species in a plasma, and as we shall see later, allow a treatment of highly ionized plasma as a single conducting
More information