Un problème inverse: l identification du profil de courant dans un Tokamak
|
|
- Curtis Fitzgerald
- 5 years ago
- Views:
Transcription
1 Un problème inverse: l identification du profil de courant dans un Tokamak Blaise Faugeras en collaboration avec J. Blum et C. Boulbe Université de Nice Sophia Antipolis Laboratoire J.-A. Dieudonné Nice, France INRIA, Mars 2009 BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
2 Outline 1 Introduction 2 Mathematical modelling of axisymmetric equilibrium 3 The inverse equilibrium problem 4 Numerical method 5 Conclusion BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
3 JET Tokamak BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
4 JET : vacuum vessel and plasma BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
5 Introduction Equilibrium of a plasma : a free boundary problem Equilibrium equation inside the plasma, in an axisymmetric configuration : Grad-Shafranov equation Right-hand side of this equation is a non-linear source : the toroidal component of the plasma current density j φ Aim of this work Perform the real-time identification of this non-linearity from experimental measurements BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
6 Mathematical modelling of the equilibrium Momentum equation : ρ( u t + u. u) + p = j B At the slow resistive diffusion time scale ρ( u t + u. u) p Equilibrium equations p = j B (Conservation of momentum).b = 0 (Conservation of B) B = µj (Ampere s law) + axisymmetric assumption => Grad-Shafranov equation BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
7 Magnetic surfaces p B = 0 and p j = 0 => Field lines and current lines are on isobaric surfaces (p = cst) = magnetic surfaces BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
8 Axisymmetric configuration B = 0, conservation of B Cylindrical coordinates (e r, e φ, e z ) B independent of φ B = (B r e r + B z e z ) + B φ e φ = B p + B φ Poloidal flux : ψ(r, z) = 1 Bds = 2π D r 0 B z rdr B p = 1 r ψ = 1 r ( ψ e φ) B. ψ = 0. ψ = cst on each magnetic surface. p = p(ψ) BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
9 B = µ 0 j, Ampere s law j = 0 j = j p + j φ. j p = 1 r ( ( f ) e φ ) µ 0 j. p = 0. f p = 0. f = f (ψ) diamagnetic function B φ = f r e φ j φ = ( ψ)e φ with. = r ( 1 µ 0 r. r ) + z ( 1 µ 0 r. z ) BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
10 Grad-Shafranov equation j B = p, equilibrium equation (j p + j φ e φ ) (B p + B φ e φ ) = 1 µ 0 r B 1 φ f + j φ ψ = p. r p = p (ψ) ψ and f = f (ψ) ψ Grad-Shafranov equation In the plasma ψ = rp (ψ) + 1 µ 0 r (ff )(ψ) In the vacuum ψ = 0 BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
11 Tokamak BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
12 Definition of the free plasma boundary Two cases : magnetic separatrix : hyperbolic line with an X-point (left) outermost flux line inside the limiter (right) BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
13 Equation for ψ(r, z) inside the vacuum vessel where ψ = λ[ r A( ψ) + R 0 R 0 r B( ψ)]χ Ωp ψ = h on Ω R 0 major radius of the torus p ( ψ) = λ 1 R 0 A( ψ) 1 (ff )( µ ψ) = λr 0 B( ψ), 0 ψ max ψ ψ Ω = ψ b max ψ [0, 1] in Ω p Ω χ Ωp is the characteristic function of Ω p in Ω BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
14 The inverse equilibrium problem : experimental measurements magnetic measurements on the vacuum vessel ψ(m i ) = g i and 1 r ψ n (N j) = h j on Ω, I p = interferometry and polarimetry on several chords n e (ψ) ψ n e (ψ)dl = α m, C m C m r n dl = β m Ω p j φ dx kinetic pressure motional Stark effect p(r, 0) = p d (r) f j (B r (M j ), B z (M j ), B φ (M j )) = γ j BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
15 BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
16 Statement of the inverse problem State equation ψ = λ[ r A( R ψ) + R 0 0 r B( ψ)]1 Ωp ψ = g in Ω on Ω Least square minimization with J(A, B, n e ) = J 0 + J 1 + J 2 + J ɛ J 0 = k w k 2(1 r J 1 = k w k 2( J 2 = k w 2 k ( ψ n (N k) h k ) 2 C k n e ( ψ) ψ r n dl α k) 2 C k n e ( ψ)dl β k ) 2 BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
17 Statement of the inverse problem (2) Tikhonov regularization 1 J ɛ = ɛ A (A (x)) 2 dx 0 1 +ɛ B (B (x)) 2 dx 0 +ɛ ne 1 0 (n e (x)) 2 dx Find A 0, B 0, n e0 such that : J(A 0, B 0, n e0 ) = inf J(A, B, n e ) BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
18 Numerical identification Finite element resolution Find ψ H 1 with ψ = g on Ω such that v H0 1 1, Ω µ 0 r ψ vdx = λ[ r A( ψ) + R 0 Ω p R 0 r B( ψ)]vdx with Fixed point A(x) = i a if i (x), B(ψ) = i b if i (x), u = (a i, b i ) Kψ = Y (ψ)u + g K modified stiffness matrix, u coefficients of A and B, g Dirichlet BC Direct solver : (ψ n, u) ψ n+1 ψ n+1 = K 1 [Y (ψ n )u + g] BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
19 Numerical identification (2) Least-square minimization J(u) = C(ψ)ψ d 2 + u T Au d : experimental measurements A : regularization terms Approximation J(u) = C(ψ n )ψ d 2 + u T Au, with ψ = K 1 [Y (ψ n )u + g] J(u) = C(ψ n )K 1 Y (ψ n )u + C(ψ n )K 1 g d 2 + u T Au = E n u F n 2 + u T Au Normal equation. Inverse solver : ψ n u (E nt E n + A)u = E nt F n BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
20 EQUINOX One equilibrium reconstruction First guess (ψ 0, u 0, λ 0, Ω 0 p) Direct and Inverse problem fixed-point iterations : ne identification (normal equation, ψ n ) Update λ n+1 = I p / Ω n p[ r R 0 A n ( ψ n ) + R 0 r Bn ( ψ n )]dx and scale DOFs (m = max a i, u 1 m u) Inverse solver : ψ n u n+1 Direct solver : (ψ n, u n+1 ) ψ n+1, Ω n+1 p Stopping condition ψn+1 ψ n ψ n < ɛ BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
21 EQUINOX Real-time Quasi-static approach. Maximum number of fixed-point iterations Cheap operations : K = LU and K 1 computed and stored once for all Normal equation : 10 basis functions small linear system Tikhonov regularization parameters unchanged Expensive operations : ( ) Cmag Update of matrices C = and Y (ψ) C polar (ψ) Computation of products C(ψ)K 1 and C(ψ)K 1 Y (ψ) BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
22 Noise free twin experiment. Magnetics only. Identified A and B BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
23 Definitions Current density averaged over magnetic surfaces r 0 < j(r, ψ) r The average < A > is defined as < A >= Safety factor q For one field line q = φ 2π. >= λa( ψ) + λ r 2 0 < r 2 > B( ψ) C ψ q( ψ) = 1 2π Adl / B p C ψ C ψ B φ rb p dl dl B p BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
24 Noise free twin experiment. Magnetics only. Mean current density and safety factor BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
25 Noisy twin experiment. Magnetics only. Identified A and B BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
26 Noisy twin experiment. Magnetics only. Mean current density and safety factor BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
27 Noisy twin experiment. Magnetics and polarimetry. Identified A and B BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
28 Noisy twin experiment. Magnetics and polarimetry. Mean current density and safety factor BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
29 Numerical Results : Tore Supra and JET characteristics ToreSupra Finite element mesh JET Number of triangles Number of nodes functions A and B Basis type Bspline Bspline Number of basis func. 8 8 Computation time (1.80GHz) One equilibrium 20 ms 60 ms Real-time requirement : 100 ms BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
30 Tore Supra. Magnetics and polarimetry. BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
31 Jet Magnetics only. BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
32 Jet Magnetics and polarimetry. BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
33 Conclusion EQUINOX : Real-time identification of the current density and equilibrium reconstruction Makes possible future real-time control of the current profile BF (Université de Nice) Identification équilibre plasma INRIA, Mars / 33
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 informationAssimilation of polarimetry Stokes vector measurements in tokamak free-boundary equilibrium reconstruction
Assimilation of polarimetry Stokes vector measurements in tokamak free-boundary equilibrium reconstruction Blaise Faugeras CASTOR Team-Project INRIA and Laboratoire J.A. Dieudonné CNRS UMR 7351 Université
More informationarxiv: v3 [math.na] 17 Mar 2011
arxiv:99.4474v3 [math.na] 7 Mar Reconstruction of the equilibrium of the plasma in a Tokamak and identification of the current density profile in real time J. Blum a, C. Boulbe a, B. Faugeras a a Laboratoire
More informationIdeal MHD Equilibria
CapSel Equil - 01 Ideal MHD Equilibria keppens@rijnh.nl steady state ( t = 0) smoothly varying solutions to MHD equations solutions without discontinuities conservative or non-conservative formulation
More informationModel based optimization and estimation of the field map during the breakdown phase in the ITER tokamak
Model based optimization and estimation of the field map during the breakdown phase in the ITER tokamak Roberto Ambrosino 1 Gianmaria De Tommasi 2 Massimiliano Mattei 3 Alfredo Pironti 2 1 CREATE, Università
More informationMagnetic Diagnostics Basics
Abstract We report measurement of the equilibrium plasma current profiles in the Levitated Dipole Experiment (LDX) that exhibit a peak beta in excess of 10 percent. The beta of an LDX plasma is calculated
More informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More information0 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 informationFinite Volume for Fusion Simulations
Finite Volume for Fusion Simulations Elise Estibals, Hervé Guillard, Afeintou Sangam To cite this version: Elise Estibals, Hervé Guillard, Afeintou Sangam. Finite Volume for Fusion Simulations. Jorek Meeting
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
More informationThe 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 informationNIMEQ: MHD Equilibrium Solver for NIMROD
NIMEQ: MHD Equilibrium Solver for NIMOD E.C.Howell, C..Sovinec University of Wisconsin-Madison 5 th Annual Meeting of Division of Plasma Physics Dallas, Texas, Nov. 17-Nov. 1,8 1 Abstract A Grad-Shafranov
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 informationMagnetic Field Models and their Application in Optimal Magnetic Divertor Design
Magnetic Field Models and their Application in Optimal Magnetic Divertor Design Maarten Blommaert, Martine Baelmans, Holger Heumann, Yannick Marandet, Hugo Bufferand, N. R. Gauger, D. Reiter To cite this
More informationReal 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 informationDIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH
DIAGNOSTICS FOR ADVANCED TOKAMAK RESEARCH by K.H. Burrell Presented at High Temperature Plasma Diagnostics 2 Conference Tucson, Arizona June 19 22, 2 134 /KHB/wj ROLE OF DIAGNOSTICS IN ADVANCED TOKAMAK
More informationGenerating 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 informationPlasma formation in MAST by using the double null merging technique
1 Plasma formation in MAST by using the double null merging technique P. Micozzi 1, F. Alladio 1, P. Costa 1, A. Mancuso 1, A. Sykes 2, G. Cunningham 2, M. Gryaznevich 2, J. Hicks 2, M. Hood 2, G. McArdle
More informationEquilibrium reconstruction improvement via Kalman-filter-based vessel current estimation at DIII-D
Fusion Engineering and Design 82 (2007) 1144 1152 Equilibrium reconstruction improvement via Kalman-filter-based vessel current estimation at DIII-D Y. Ou a, M.L. Walker b, E. Schuster a,, J.R. Ferron
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
More informationComputations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations
Computations of Vector Potential and Toroidal Flux and Applications to Stellarator Simulations NIMROD Team Meeting Torrin Bechtel April 30, 2017 Outline 1 Project Goals and Progress 2 Vector Potential
More informationEffect of ideal kink instabilities on particle redistribution
Effect of ideal kink instabilities on particle redistribution H. E. Ferrari1,2,R. Farengo1, P. L. Garcia-Martinez2, M.-C. Firpo3, A. F. Lifschitz4 1 Comisión Nacional de Energía Atómica, Centro Atomico
More informationPOWER DENSITY ABSORPTION PROFILE IN TOKAMAK PLASMA WITH ICRH
Dedicated to Professor Oliviu Gherman s 80 th Anniversary POWER DENSITY ABSORPTION PROFILE IN TOKAMAK PLASMA WITH ICRH N. POMETESCU Association EURATOM-MECTS Romania, University of Craiova, Faculty of
More informationParallel Kelvin-Helmholtz instability in edge plasma
Journal of Physics: Conference Series OPEN ACCESS Parallel Kelvin-Helmholtz instability in edge plasma To cite this article: H Guillard et al 2014 J. Phys.: Conf. Ser. 561 012009 View the article online
More informationThe 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 informationIntegrated equilibrium reconstruction and MHD stability analysis of tokamak plasmas in the EU-IM platform
EUROFUSION WPCD-PR(16) 15379 R Coelho et al. Integrated equilibrium reconstruction and MHD stability analysis of tokamak plasmas in the EU-IM platform Preprint of Paper to be submitted for publication
More informationIs the Troyon limit a beta limit?
Is the Troyon limit a beta limit? Pierre-Alexandre Gourdain 1 1 Extreme State Physics Laboratory, Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA The plasma beta,
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationPhysic-based Preconditioning and B-Splines finite elements method for Tokamak MHD
Physic-based Preconditioning and B-Splines finite elements method for Tokamak MHD E. Franck 1, M. Gaja 2, M. Mazza 2, A. Ratnani 2, S. Serra Capizzano 3, E. Sonnendrücker 2 ECCOMAS Congress 2016, 5-10
More informationReconstruction of the Pressure Profile of LDX High Beta Plasma
Reconstruction of the Pressure Profile of LDX High Beta Plasma M. Mauel, I. Karim, A. Boxer, J. Ellsworth, D. Garnier, A. Hansen, J. Kesner, and E. Ortiz Dept. of Applied Physics and Applied Mathematics,
More informationDesign of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system
Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system J.M. Ané 1, V. Grandgirard, F. Albajar 1, J.Johner 1 1Euratom-CEA Association, Cadarache, France Euratom-EPFL
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 informationRole of the Electron Temperature in the Current Decay during Disruption in JT-60U )
Role of the Electron Temperature in the Current Decay during Disruption in JT-60U ) Yoshihide SHIBATA, Akihiko ISAYAMA, Go MATSUNAGA, Yasunori KAWANO, Seiji MIYAMOTO 1), Victor LUKASH 2), Rustam KHAYRUTDINOV
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 informationConfiguration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift )
Configuration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift ) Shoichi OKAMURA 1,2) 1) National Institute for Fusion Science, Toki 509-5292, Japan 2) Department of Fusion Science,
More informationMHD 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 informationHeating and Confinement Study of Globus-M Low Aspect Ratio Plasma
EX/P5- Heating and Confinement Study of Globus-M Low Aspect Ratio Plasma N.V. Sakharov ), V.V. Dyachenko ), B.B. Ayushin ), A.V. Bogomolov ), F.V. Chernyshev ), V.K. Gusev ), S.A. Khitrov ), N.A. Khromov
More informationIon plateau transport near the tokamak magnetic axis. Abstract
Ion plateau transport near the tokamak magnetic axis K.C. Shaing and R.D. Hazeltine Institute for Fusion Studies, The University of Texas at Austin, Austin, Texas 78712 (December 2, 1997) Abstract Conventional
More informationTokamak Fusion Basics and the MHD Equations
MHD Simulations for Fusion Applications Lecture 1 Tokamak Fusion Basics and the MHD Equations Stephen C. Jardin Princeton Plasma Physics Laboratory CEMRACS 1 Marseille, France July 19, 21 1 Fusion Powers
More informationGeneralized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration
Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration D.A. Kaltsas and G.N. Throumoulopoulos Department of Physics, University of Ioannina, GR 451 10 Ioannina,
More informationFull-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies
Full-wave Electromagnetic Field Simulations in the Lower Hybrid Range of Frequencies P.T. Bonoli, J.C. Wright, M. Porkolab, PSFC, MIT M. Brambilla, IPP, Garching, Germany E. D Azevedo, ORNL, Oak Ridge,
More informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
More informationProspects of Nuclear Fusion Energy Research in Lebanon and the Middle-East
Prospects of Nuclear Fusion Energy Research in Lebanon and the Middle-East Ghassan Antar Physics Department American University of Beirut http://www.aub.edu.lb/physics/lpfd Outline 1. Introduction and
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
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 informationDKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron Drift Kinetic Equation
DKE: a fast numerical solver for the 3-D relativistic bounce-averaged electron Drift Kinetic Equation Joan Decker and Y Peysson Association EURATOM-CEA sur la Fusion CEA-Cadarache, F-38 Saint Paul-lez-Durance,
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationCalculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions
Calculation of alpha particle redistribution in sawteeth using experimentally reconstructed displacement eigenfunctions R. Farengo, H. E. Ferrari,2, M.-C. Firpo 3, P. L. Garcia-Martinez 2,3, A. F. Lifschitz
More informationModeling 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 informationToroidal Multipolar Expansion for Fast L-Mode Plasma Boundary Reconstruction in EAST
Plasma Science and Technology, Vol.3, No.3, Jun. 2 Toroidal Multipolar Expansion for Fast L-Mode Plasma Boundary Reconstruction in EAST GUO Yong ), XIAO Bingjia ), LUO Zhengping ) Institute of Plasma Physics,
More informationMHD equilibrium calculations for stellarators. Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE)
MHD equilibrium calculations for stellarators Antoine Cerfon, MIT PSFC with F. Parra, J. Freidberg (MIT NSE) March 20, 2012 MAGNETIC FIELD LINE HAMILTONIAN Consider a general toroidal coordinate system
More informationIntegrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport
1 Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport N. Hayashi, T. Takizuka, T. Ozeki, N. Aiba, N. Oyama Japan Atomic Energy Agency, Naka, Ibaraki-ken, 311-0193 Japan
More informationInternal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak
Journal of Physics: Conference Series Internal transport barriers in some Hamiltonian systems modeling the magnetic lines dynamics in tokamak To cite this article: D Constantinescu et al 2005 J. Phys.:
More information(Motivation) Reactor tokamaks have to run without disruptions
Abstract A database has been developed to study the evolution, the nonlinear effects on equilibria, and the disruptivity of locked and quasi-stationary modes with poloidal and toroidal mode numbers m=2
More informationDT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift )
DT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift ) Tsuguhiro WATANABE National Institute for Fusion Science, 322-6 Oroshi-cho, Toki 509-5292, Japan (Received
More informationImplicit kinetic relaxation schemes. Application to the plasma physic
Implicit kinetic relaxation schemes. Application to the plasma physic D. Coulette 5, E. Franck 12, P. Helluy 12, C. Courtes 2, L. Navoret 2, L. Mendoza 2, F. Drui 2 ABPDE II, Lille, August 2018 1 Inria
More informationOverview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields
Overview of edge modeling efforts for advanced divertor configurations in NSTX-U with magnetic perturbation fields H. Frerichs, O. Schmitz, I. Waters, G. P. Canal, T. E. Evans, Y. Feng and V. Soukhanovskii
More informationAbstract. Department of Physics. University of Wisconsin UCLA UCLA
Abstract Internal magnetic field fluctuations and equilibrium poloidal magnetic field have been measured in the MST reversal field pinch by a 11 chord far-infrared polarimeterinterferometer system with
More informationDetermination of q Profiles in JET by Consistency of Motional Stark Effect and MHD Mode Localization
Determination of q Profiles in JET by Consistency of Motional Stark Effect and MHD Mode Localization R. De Angelis 1), F. Orsitto1), M. Baruzzo 2), P. Buratti 1), B. Alper 3), L. Barrera 4), A. Botrugno
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationMassachusetts Institute of Technology 22.68J/2.64J Superconducting Magnets. February 27, Lecture #4 Magnetic Forces and Stresses
Massachusetts Institute of Technology.68J/.64J Superconducting Magnets February 7, 003 Lecture #4 Magnetic Forces and Stresses 1 Forces For a solenoid, energy stored in the magnetic field acts equivalent
More informationNumerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations
PHYSICS OF PLASMAS 12, 072513 2005 Numerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations J. N. Talmadge and S. P. Gerhardt a HSX Plasma Laboratory, University
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationBunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (
Title The Finite Beta Effects on the Toro Tokamak Plasma Author(s) Bunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (2013 Issue Date 2013-02 URL http://hdl.handle.net/2433/173038
More informationCurrent Holes and other Structures in Motional Stark Effect Measurements
Max-Planck-Institut für Plasmaphysik Garching ASDEX Upgrade Current Holes and other Structures in Motional Stark Effect Measurements Dissertation von Doris Merkl Technische Universität München ii Technische
More informationNon-Solenoidal Plasma Startup in
Non-Solenoidal Plasma Startup in the A.C. Sontag for the Pegasus Research Team A.C. Sontag, 5th APS-DPP, Nov. 2, 28 1 Point-Source DC Helicity Injection Provides Viable Non-Solenoidal Startup Technique
More informationInfluence of Beta, Shape and Rotation on the H-mode Pedestal Height
Influence of Beta, Shape and Rotation on the H-mode Pedestal Height by A.W. Leonard with R.J. Groebner, T.H. Osborne, and P.B. Snyder Presented at Forty-Ninth APS Meeting of the Division of Plasma Physics
More informationAPPENDIX 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 informationGuiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates
Contrib. Plasma Phys. 48, No. -3, 4 8 (8) / DOI./ctpp.839 Guiding Center Orbit Studies in a Tokamak Edge Geometry Employing Boozer and Cartesian Coordinates Y. Nishimura,Y.Xiao,and Z. Lin Department of
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationEdge Momentum Transport by Neutrals
1 TH/P3-18 Edge Momentum Transport by Neutrals J.T. Omotani 1, S.L. Newton 1,2, I. Pusztai 1 and T. Fülöp 1 1 Department of Physics, Chalmers University of Technology, 41296 Gothenburg, Sweden 2 CCFE,
More informationProgress on Optimization Results for Steady-State Plasma Solutions
Progress on Optimization Results for Steady-State Plasma Solutions A.D. Turnbull, H. St John, R. Buttery, M. Choi, L.L Lao, S. Smith General Atomics ARIES Team Meeting Bethesda Md September 26, 27 2012
More informationHighlights from (3D) Modeling of Tokamak Disruptions
Highlights from (3D) Modeling of Tokamak Disruptions Presented by V.A. Izzo With major contributions from S.E. Kruger, H.R. Strauss, R. Paccagnella, MHD Control Workshop 2010 Madison, WI ..onset of rapidly
More informationForce-Free Magnetosphere of an Accreting Kerr Black Hole
Force-Free Magnetosphere of an Accreting Kerr Black Hole Dmitri A. Uzdensky Princeton University, Princeton, NJ 08540, USA I consider a stationary axisymmetric force-free degenerate magnetosphere of a
More informationPerformance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK
Performance limits Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 24) Previously... In the last few
More informationSpring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014)
Spring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014) The first problem in this assignment is a paper-and-pencil exercise
More informationElectrode Biasing Experiment in the Large Helical Device
1 EXC/P8-07 Electrode Biasing Experiment in the Large Helical Device S. Kitajima 1), H. Takahashi 2), K. Ishii 1), J. Sato 1), T. Ambo 1), M. Kanno 1), A. Okamoto 1), M. Sasao 1), S. Inagaki 3), M. Takayama
More informationComputations with Discontinuous Basis Functions
Computations with Discontinuous Basis Functions Carl Sovinec University of Wisconsin-Madison NIMROD Team Meeting November 12, 2011 Salt Lake City, Utah Motivation The objective of this work is to make
More informationPlasma models for the design of the ITER PCS
Plasma models for the design of the ITER PCS G. De Tommasi 1,2 on behalf of the CREATE team 1 Consorzio CREATE, Naples, Italy 2 Department of Electrical Engineering and Information Technology, University
More informationIntegration of Fokker Planck calculation in full wave FEM simulation of LH waves
Integration of Fokker Planck calculation in full wave FEM simulation of LH waves O. Meneghini S. Shiraiwa R. Parker 51 st DPP APS, Atlanta November 4, 29 L H E A F * Work supported by USDOE awards DE-FC2-99ER54512
More informationNIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU. General Atomics. Nimrod Project Review Meeting July 21 22, 1997
NIMROD FROM THE CUSTOMER S PERSPECTIVE MING CHU General Atomics Nimrod Project Review Meeting July 21 22, 1997 Work supported by the U.S. Department of Energy under Grant DE-FG03-95ER54309 and Contract
More informationUnderstanding Edge Harmonic Oscillation Physics Using NIMROD
Understanding Edge Harmonic Oscillation Physics Using NIMROD J. King With contributions from S. Kruger & A. Pankin (Tech-X); K. Burrell, A. Garofalo, R. Groebner & P. Snyder (General Atomics) Work supported
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 informationSAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD
1 SAROD Conference, Hyderabad, december 2005 CONTINUOUS MESH ADAPTATION MODELS FOR CFD Bruno Koobus, 1,2 Laurent Hascoët, 1 Frédéric Alauzet, 3 Adrien Loseille, 3 Youssef Mesri, 1 Alain Dervieux 1 1 INRIA
More informationResponse of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Alfven Frequency
Response of a Resistive and Rotating Tokamak to External Magnetic Perturbations Below the Alfven Freuency by M.S. Chu In collaboration with L.L. Lao, M.J. Schaffer, T.E. Evans E.J. Strait (General Atomics)
More informationHeat Flux Management via Advanced Magnetic Divertor Configurations and Divertor Detachment.
Heat Flux Management via Advanced Magnetic Divertor Configurations and Divertor Detachment E. Kolemen a, S.L. Allen b, B.D. Bray c, M.E. Fenstermacher b, D.A. Humphreys c, A.W. Hyatt c, C.J. Lasnier b,
More informationThe RFP: Plasma Confinement with a Reversed Twist
The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed
More information3D Reconstruction of Plasma Equilibria using Magnetic Diagnostics on the Compact Toroidal Hybrid. Benjamin Adam Stevenson
3D Reconstruction of Plasma Equilibria using Magnetic Diagnostics on the Compact Toroidal Hybrid by Benjamin Adam Stevenson A dissertation submitted to the Graduate Faculty of Auburn University in partial
More informationFormation 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 informationToroidal confinement of non-neutral plasma. Martin Droba
Toroidal confinement of non-neutral plasma Martin Droba Contents Experiments with toroidal non-neutral plasma Magnetic surfaces CNT and IAP-high current ring Conclusion 2. Experiments with toroidal non-neutral
More informationGA A22863 PLASMA PRESSURE AND FLOWS DURING DIVERTOR DETACHMENT
GA A22863 PLASMA PRESSURE AND FLOWS DURING DIVERTOR DETACHMENT by M.J. SCHAFFER, J.A. BOEDO, N.H. BROOKS, R.C. ISLER, and R.A. MOYER AUGUST 1998 DISCLAIMER This report was prepared as an account of work
More informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationNautilus: 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 informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationDisruption mitigation in ITER
Disruption mitigation in ITER S.Putvinski, R.Pitts, M.Sugihara, L.Zakharov Page 1 Outline Introduction Massive gas injection for mitigation of thermal loads Forces on VV and in vessel components Suppression
More informationToroidal flow stablization of disruptive high tokamaks
PHYSICS OF PLASMAS VOLUME 9, NUMBER 6 JUNE 2002 Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland 20742-3511 Received 4 February 2002; accepted
More informationMechanisms of intrinsic toroidal rotation tested against ASDEX Upgrade observations
Mechanisms of intrinsic toroidal rotation tested against ASDEX Upgrade observations William A. Hornsby C. Angioni, E. Fable, P. Manas, R. McDermott, Z.X. Lu, S. Grosshauser 2, A. G. Peeters 2 and the ASDEX
More information