Modeling of ELM Dynamics for ITER


 Lilian Skinner
 1 years ago
 Views:
Transcription
1 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 University of Tulsa, Tulsa, Oklahoma 3 TechX, Boulder, CO General Atomics, San Diego, CA ITPA meeting October, Chengdu, China
2 Outline Goal: Integrated modeling of pedestal and ELMs Constrains of numerical modeling Elements of reduced model for ELMs Triggering conditions for ELMs Equlibria generated with TOQ and TEQ codes Ideal MHD stability analysis with ELITE, DCON, and BALOO codes Ballooning/peeling marginal stability condition parameterized for use as ELM trigger within integrated simulations Integrated modeling simulation of ELMy Hmode plasma with ASTRA code Simulation of ELM crash with the MHD code NIMROD Nonideal effects in NIMROD code Resistivity, viscosity, anisotropic thermal transport, flow shear, twofluid and FLR effects Linear results obtained using NIMROD code Compared with corresponding results from other MHD codes Nonlinear ELM evolution computed using NIMROD code Mode coupling, filaments and explosive growth observed Effect of flow shear on ELM evolution
3 Numerical Constraints [ ] [ # of dimensions Mesh points in each dimension ] Q 7 εp = 3 10 [ time step < sec> ] [# of runs per year][ Problem time < sec> ] Algorithms Constraints Q  codealgorithm requirements P  peak hardware performance ε  hardware efficiency Tflop meshpoint timestep Tflop sec Algorithm performance and problem requirements with available cycles Requirements: Required problem time: at least 1 sec, at least 500 runs per year Conclusions: 3D (i.e., fluid) calculations for times of ~ 10 msec (single ELM cycle) within reach Longer times require next generation computers (or better algorithms) Higher dimensional (kinetic) long time calculations unrealistic Integrated effects must come through low dimensionality closures (transport modeling) Assumption: time step= sec Mesh points in each dimensions= cases/year Q= TFlop/meshpoint/timestep
4 ELM Model for Integrated Simulations A long range goal is to develop a model for ELMs for use in predictive whole device modeling code To simulate ITER plasmas on macro time scales Elements needed for integrated model for ELMs Model for triggering conditions for ELM crashes Model for ELM width Model for particle, current and heat losses Important effects to include Flow shear effect Effect of neutrals Effect of dust particles Effect of NTMs and sawteeth on pedestal and ELMs Effect of Xpoint and vertical asymmetry FLR and twofluid effects Effect of scrapeoflayer
5 ITER Equilibrium TOQ code inverse equilibrium solver Fast and can be used in a batch mode with the ideal MHD stability codes to generate of ELM stability diagram. Equilibrium extended into scrapeoff region with supplementary code VACUUM, called from MHD stability code DCON Nonideal MHD NIMROD code requires an even more robust equilibrium that includes the scrapeoff region TEQ code direct equilibrium solver Extracted from CORSICA code as NTCC module Can be used both for prescribed boundary and free boundary equilibria Parameterized pressure and current density profiles generated in a similar manner as with TOQ code Freeboundary TEQ equilibrium solver applied to generate new equilibria, including scrapeoff region, for use with nonideal MHD NIMROD code
6 ITER PeelingBallooning Stability Map Ideal MHD stability codes, DCON, ELITE, and BALOO, used to produce peeling/ballooning stability map in the pedestal region Stability boundary parameterized for ELM threshold condition Stability boundaries for the ITER normalized pressure gradient and bootstrap current are comparable to the values obtained for the DIIID high triangularity case described by M. Murakami et al., Nucl. Fusion 40, 1257 (2000) Parameters for the ITER reference case: a = 2.0 m, R = 6.2 m, δ = 0.49, κ = 1.85, B 0 = 5.3 T, I = 15.0 MA, n 0 = m 3, n ped = m 3, T 0 = 20.0 kev, T ped = 4.0 kev 100 j (A/cm 2 ) 50 ITER Peeling Stable 100 j (A/cm 2 ) 50 DIIID High Triangularity Peeling Stable 0 0 Ballooning 2 α Ballooning μ0 p V V α= π ψ ψ 2 π R α 2 4 6
7 Integrated Approach to Edge Modeling Vision for integrated (whole device) modeling with edge Currently, integrated modeling codes include: 1. Transport modules are called directly from code: GLF23 and MMM models for anomalous transport NCLASS for neoclassical transport 2. Ideal MHD stability codes used to parametrize peelingballooning triggering conditions for ELM crashes BALOO, DCON, ELITE, and MISHKA In the future, integrated modeling codes will include: 3. Nonlinear NIMROD MHD code results for ELM crash Including effects of resistively, viscosity, anisotropic heat flux, FLR and twofluid effects 4. SOLPS (B2/EIRINE) code results for effects of neutrals
8 Integration of Ideal MHD Stability Analysis Within Transport Simulations Results of ASTRA simulations of DIIID using the new parameterized ideal MHD stability model ELM frequency increases with heating power and triangularity Pedestal and ELMs form spontaneously after auxiliary heating power turned on A C B
9 ELM Modeling with NIMROD Code The NIMROD code numerically advances the resistive MHD equations in 3D geometry B = E + κ divb B t E= V B+ηJ μ 0 J = B n + ( nv) = D n t V ρ + V V = J B p Π t Highorder finite element representation of the poloidal plane: Accuracy for MHD and transport anisotropy at realistic parameters: S>10 6, χ /χ perp >10 9 n T T p n ( ) ˆ ˆ + V = V + χ χ + χ T + Q γ 1 t bb I visc Flexible spatial representation Temporal advance with semiimplicit and implicit methods Multiple timescale physics from ideal MHD (ms) to transport (ms) time scales
10 Nonlinear ELM Modeling with NIMROD Code Electron temperature Velocity field
11 DIIID Equilibria for Stability Analysis Sequence of DIIIDlike equilibria Mode structure changes with pedestal height The linear growth rates as a function of mode number for different equilibria as computed by NIMROD
12 ELM Modeling for ITER with NIMROD Code Contour plots of vector and scalar fields computed with nonlinear NIMROD simulation of an ELM crash for ITER (at 400 time steps)
13 Nonlinear ELM Modeling for ITER with NIMROD Code Time dynamics of kinetic energies for different mode numbers for ITER This nonlinear NIMROD simulation is limited to 21 toroidal modes The n=21 mode is being driven by the unstable modes with lower mode number More modes will have to be included in the computation to demonstrate convergence of the nonlinear mode coupling effects Accurate account for twofluid and FLR effects might also stabilize modes with high mode numbers
14 NIMROD Code: Twofluid and FLR terms In addition to previously added nonideal effects, recent advances have enabled twofluid treatment Hall and diamagnetic terms in Ohm s law More complete stress tensor in momentum equation V ρ + V V = p + J B Πvisc Π Π i t 1 E= V{ B + ( J B p ) e Π e + η { J Ideal MHD ne Resistive MHD Twofluid effects (Hall + diamagnetic) + Equations for p and p, + closures i e gv Π gv = p 4Ω [( b W) ( I + 3bb) + transpose] W = V + V T 2 3 I V
15 Extended MHD Effects Important Linearly The drift velocity is in the poloidal direction, which can be stabilizing to the edge mode, especially nonlinearly, is a critically important effect Linear Growth Rate B ω * ψ Hall and Gyroviscous Effects Included Single Fluid Stable Poloidal Drift Frequency c ω *e,i = q e,i nb k (B p ) 2 θ e,i θ Computationally, the twofluid model with drift effects requires more temporal resolution than the MHD model. Nonlinear effect on ELMs not yet explored. Complete picture includes both toroidal and poloidal flows nonlinearly
16 Effect of Flow Shear on ELM Evolution Flow destabilizes linear spectrum especially at high n in case shown which can possibly change qualitative linear picture Ideal MHD high n modes stabilized Nonlinear energy spectrum similar to that without flow although nonlinear evolution strongly affected by flow Toroidal Flow Shear Only B V No Flow resonance Flow ψ Linear
17 With Flow Shear the Mode Structure is Limited in Radial Propagation Without flow, the high temperature filaments propagate to the wall With flow, the perturbation is dragged by the fluid, significant localization in long term evolution Toroidal Flow Shear Only B V T e ψ Without flow With flow Filaments Sheared Lower amplitude poloidal fluctuation and reduced radial gradients with flow shear. Possible healing. Toroidal flow shear must now be combined with drift effects for completeness
18 Summary Progress in development of reduced model for ELM for integrated modeling transport simulation is reported Triggering conditions for ELMs in ITER and DIIID studied with ideal MHD stability codes BALOO, DCON and ELITE ELM crash studied with nonideal MHD code NIMROD Linear and nonlinear stages of ELM crash investigated Nonideal MHD effects included in NIMROD Additional effects include resistivity (S>10 6 ), viscosity, particle transport, parallel and perpendicular thermal transport (χ /χ >10 9 ) Hall and diamagnetic terms added to Ohm s law Stress tensor in momentum equation Linear growth rate is enhanced by flow shear Growth rate strongly affected by resistive SOL region Filaments are observed during nonlinear evolution of ELM crash Filaments are dragged and twisted by flow shear before they extend all the way to the wall Results being used to construct reduced models for whole device for integrated modeling simulations
19 Extra Slides
20 Transport model for ELMs in ASTRA code Theorybased model for pedestal and ELMs at the edge of Hmode plasmas Flow and magnetic shear stabilization of anomalous transport χ = χ F ( s, ω ) +χ F ( s, ω ) +χ +χ i i i i i ITG ITG ExB RB RB ExB KB i,neo χ = χ F ( s, ω ) +χ F ( s, ω ) +χ e e e e e TEM TEM ExB RB RB ExB KB ETG e,neo j where are shear suppression functions F F l () 1 αlj j s l ( s, ω ExB) = (2) 2 1+ αlj ( ωexbτlj) () 1 ( 2 ) ( j ) lj lj l, where +χ + χ ( ITG / TEM, RB) l = j = ( ions, electrons) α and α are adjustable constants, χ is anomalous thermal diffusivity caused by () l mode turbulence, ωexb is flow shear ExB rate, and τlj are 2 ( j ) turbulence correlation times, which are estimated as τ = L χ. lj l l Pedestal width and critical pressure gradients are parameters that responsible for the pressure at the top of the pedestal Pedestal width is result from the model of shear suppression of anomalous transport Critical pressure gradient can be computed with ideal MHS stability codes
ELMs and Constraints on the HMode Pedestal:
ELMs and Constraints on the HMode Pedestal: A Model Based on PeelingBallooning Modes P.B. Snyder, 1 H.R. Wilson, 2 J.R. Ferron, 1 L.L. Lao, 1 A.W. Leonard, 1 D. Mossessian, 3 M. Murakami, 4 T.H. Osborne,
More informationGA A25229 FINAL REPORT ON THE OFES ELM MILESTONE FOR FY05
GA A25229 ELM MILESTONE FOR FY05 by D.P. BRENNAN, E.D. HELD, S.E. KRUGER, A.Y. PANKIN, D.D. SCHNACK, AND C.R. SOVINEC APRIL 2006 DISCLAIMER This report was prepared as an account of work sponsored by an
More informationMHD Analysis of the Tokamak Edge Pedestal in the Low Collisionality Regime Thoughts on the Physics of ELMfree QH and RMP Discharges
MHD Analysis of the Tokamak Edge Pedestal in the Low Collisionality Regime Thoughts on the Physics of ELMfree QH and RMP Discharges P.B. Snyder 1 Contributions from: H.R. Wilson 2, D.P. Brennan 1, K.H.
More informationTH/P84 Second Ballooning Stability Effect on Hmode Pedestal Scalings
TH/P84 Second Ballooning Stability Effect on Hmode Pedestal Scalings T. Onjun 1), A.H. Kritz ), G. Bateman ), A. Pankin ) 1) Sirindhorn International Institute of Technology, Klong Luang, Pathumthani,
More informationEdge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas
Edge Rotational Shear Requirements for the Edge Harmonic Oscillation in DIII D Quiescent H mode Plasmas by T.M. Wilks 1 with A. Garofalo 2, K.H. Burrell 2, Xi. Chen 2, P.H. Diamond 3, Z.B. Guo 3, X. Xu
More informationUnderstanding Edge Harmonic Oscillation Physics Using NIMROD
Understanding Edge Harmonic Oscillation Physics Using NIMROD J. King With contributions from S. Kruger & A. Pankin (TechX); K. Burrell, A. Garofalo, R. Groebner & P. Snyder (General Atomics) Work supported
More informationCharacteristics of the Hmode H and Extrapolation to ITER
Characteristics of the Hmode H Pedestal and Extrapolation to ITER The Hmode Pedestal Study Group of the International Tokamak Physics Activity presented by T.Osborne 19th IAEA Fusion Energy Conference
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, Ibarakiken, 3110193 Japan
More informationEFFECT OF EDGE NEUTRAL SOUCE PROFILE ON HMODE PEDESTAL HEIGHT AND ELM SIZE
EFFECT OF EDGE NEUTRAL SOUCE PROFILE ON HMODE PEDESTAL HEIGHT AND ELM SIZE T.H. Osborne 1, P.B. Snyder 1, R.J. Groebner 1, A.W. Leonard 1, M.E. Fenstermacher 2, and the DIIID Group 47 th Annual Meeting
More informationInfluence of Beta, Shape and Rotation on the Hmode Pedestal Height
Influence of Beta, Shape and Rotation on the Hmode Pedestal Height by A.W. Leonard with R.J. Groebner, T.H. Osborne, and P.B. Snyder Presented at FortyNinth APS Meeting of the Division of Plasma Physics
More informationExponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1
Exponential Growth and Filamentary Structure of Nonlinear Ballooning Instability 1 Ping Zhu in collaboration with C. C. Hegna and C. R. Sovinec University of WisconsinMadison Sherwood Conference Denver,
More informationSelfconsistent modeling of ITER with BALDUR integrated predictive modeling code
Selfconsistent modeling of ITER with BALDUR integrated predictive modeling code Thawatchai Onjun Sirindhorn International Institute of Technology, Thammasat University, Klong Luang, Pathumthani, 12121,
More informationProgress in characterization of the Hmode pedestal
Journal of Physics: Conference Series Progress in characterization of the Hmode pedestal To cite this article: A W Leonard 2008 J. Phys.: Conf. Ser. 123 012001 View the article online for updates and
More informationGA A26891 A FIRST PRINCIPLES PREDICTIVE MODEL OF THE PEDESTAL HEIGHT AND WIDTH: DEVELOPMENT, TESTING, AND ITER OPTIMIZATION WITH THE EPED MODEL
GA A26891 A FIRST PRINCIPLES PREDICTIVE MODEL OF THE PEDESTAL HEIGHT AND WIDTH: DEVELOPMENT, TESTING, AND ITER OPTIMIZATION WITH THE EPED MODEL by P.B. SNYDER, R.J. GROEBNER, J.W. HUGHES, T.H. OSBORNE,
More informationInnovative Concepts Workshop Austin, Texas February 1315, 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 informationComparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry
Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry D. D. Schnack University of Wisconsin Madison Jianhua Cheng, S. E. Parker University of Colorado
More informationLH transitions driven by ion heating in scrapeoff layer turbulence (SOLT) model simulations
LH transitions driven by ion heating in scrapeoff 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 informationSawteeth in Tokamaks and their relation to other TwoFluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other TwoFluid 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 informationTheory, Simulation and Modelling for the FY2011 Pedestal Milestone
Theory, Simulation and Modelling for the FY2011 Pedestal Milestone Phil Snyder for the ECC General Atomics, San Diego CA, USA ECC Conference Call August 21, 2009 Draft 2011 Theory Milestone The performance
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 3110102, Japan 1) University
More informationControl of Neoclassical tearing mode (NTM) in advanced scenarios
FIRST CHENGDU THEORY FESTIVAL Control of Neoclassical tearing mode (NTM) in advanced scenarios ZhengXiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline
More informationGA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIIID
GA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIIID by C.C. PETTY, P.A. POLITZER, R.J. JAYAKUMAR, T.C. LUCE, M.R. WADE, M.E. AUSTIN, D.P. BRENNAN, T.A. CASPER, M.S. CHU, J.C. DeBOO, E.J. DOYLE,
More informationGA A25592 STABILITY AND DYNAMICS OF THE EDGE PEDESTAL IN THE LOW COLLISIONALITY REGIME: PHYSICS MECHANISMS FOR STEADYSTATE ELMFREE OPERATION
GA A25592 STABILITY AND DYNAMICS OF THE EDGE PEDESTAL IN THE LOW COLLISIONALITY REGIME: PHYSICS MECHANISMS FOR STEADYSTATE ELMFREE OPERATION by P.B. SNYDER, K.H. BURRELL, H.R. WILSON, M.S. CHU, M.E.
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 informationGA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIIID
GA A27857 IMPACT OF PLASMA RESPONSE ON RMP ELM SUPPRESSION IN DIIID by A. WINGEN, N.M. FERRARO, M.W. SHAFER, E.A. UNTERBERG, T.E. EVANS, D.L. HILLIS, and P.B. SNYDER JULY 2014 DISCLAIMER This report was
More informationComputational Modeling of Fully Ionized Magnetized Plasmas Using the Fluid Approximation
Computational Modeling of Fully Ionized Magnetized Plasmas Using the Fluid Approximation Dalton D. Schnack Center for Energy and Space Science Center for Magnetic SelfOrganization in Space and Laboratory
More informationW.A. HOULBERG Oak Ridge National Lab., Oak Ridge, TN USA. M.C. ZARNSTORFF Princeton Plasma Plasma Physics Lab., Princeton, NJ USA
INTRINSICALLY STEADY STATE TOKAMAKS K.C. SHAING, A.Y. AYDEMIR, R.D. HAZELTINE Institute for Fusion Studies, The University of Texas at Austin, Austin TX 78712 USA W.A. HOULBERG Oak Ridge National Lab.,
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, EURATOMCIEMAT, Madrid, Spain Recent experiments have shown the importance of multiscale (longrange)
More informationPossible Experimental Tests of Pedestal Width/ Structure, Motivated by Theory/Modeling
Possible Experimental Tests of Pedestal Width/ Structure, Motivated by Theory/Modeling R.J. Groebner PERSISTENT SURVEILLANCE FOR PIPELINE PROTECTION AND THREAT INTERDICTION Presented at DIIID Pedestal
More informationINTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS
INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS A. Ishizawa and N. Nakajima National Institute for Fusion Science F. L. Waelbroeck, R. Fitzpatrick, W. Horton Institute for Fusion Studies, University
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft Xray. 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), XZ. Tang 2) 1) Massachusetts Institute of Technology,
More informationTRANSPORT PROGRAM CMOD 5 YEAR REVIEW MAY, 2003 PRESENTED BY MARTIN GREENWALD MIT PLASMA SCIENCE & FUSION CENTER
TRANSPORT PROGRAM CMod CMOD 5 YEAR REVIEW MAY, 2003 PRESENTED BY MARTIN GREENWALD MIT PLASMA SCIENCE & FUSION CENTER CMOD  OPPORTUNITIES AND CHALLENGES Prediction and control are the ultimate goals
More informationHybrid KineticMHD simulations with NIMROD
in NIMROD simulations with NIMROD Charlson C. Kim 1 Dylan P. Brennan 2 Yasushi Todo 3 and the NIMROD Team 1 University of Washington, Seattle 2 University of Tulsa 3 NIFS, TokiJapan December 2&3, 2011
More informationDevelopment and Validation of a Predictive Model for the Pedestal Height (EPED1)
Development and Validation of a Predictive Model for the Pedestal Height (EPED1) P.B. Snyder 1 with R.J. Groebner 1, A.W. Leonard 1, T.H. Osborne 1, M. Beurskens 3, L.D. Horton 4, A.E. Hubbard 5, J.W.
More informationDEPENDENCE OF THE HMODE PEDESTAL STRUCTURE ON ASPECT RATIO
21 st IAEA Fusion Energy Conference Chengdu, China Oct. 1621, 2006 DEPENDENCE OF THE HMODE PEDESTAL STRUCTURE ON ASPECT RATIO R. Maingi 1, A. Kirk 2, T. Osborne 3, P. Snyder 3, S. Saarelma 2, R. Scannell
More informationITER operation. Ben Dudson. 14 th March Department of Physics, University of York, Heslington, York YO10 5DD, UK
ITER operation Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 14 th March 2014 Ben Dudson Magnetic Confinement Fusion (1 of 18) ITER Some key statistics for ITER are:
More informationCurrentdriven instabilities
Currentdriven instabilities Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 21 st February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 23) Previously In the last
More informationDriftDriven and TransportDriven Plasma Flow Components in the Alcator CMod Boundary Layer
DriftDriven and TransportDriven Plasma Flow Components in the Alcator CMod Boundary Layer N. Smick, B. LaBombard MIT Plasma Science and Fusion Center PSI19 San Diego, CA May 25, 2010 Boundary flows
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade MaxPlanckInstitut 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 informationNonlinear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations.
Nonlinear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations. M. Becoulet 1, F. Orain 1, G.T.A. Huijsmans 2, P. Maget 1, N. Mellet 1, G. DifPradalier 1, G. Latu 1, C. Passeron
More informationNew bootstrap current formula valid for steep edge pedestal, and its implication to pedestal stability
1 TH/P412 New bootstrap current formula valid for steep edge pedestal, and its implication to pedestal stability C.S. Chang 1,2, Sehoon Koh 2,*, T. Osborne 3, R. Maingi 4, J. Menard 1, S. Ku 1, Scott
More informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center  University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
More informationPedestal Stability and Transport on the Alcator CMod Tokamak: Experiments in Support of Developing Predictive Capability
1 EX/P415 Pedestal Stability and Transport on the Alcator CMod Tokamak: Experiments in Support of Developing Predictive Capability J.W. Hughes 1, P.B. Snyder 2, X. Xu 3, J.R. Walk 1, E.M. Davis 1, R.M.
More informationProgress in Modeling of ARIES ACT Plasma
Progress in Modeling of ARIES ACT Plasma And the ARIES Team A.D. Turnbull, R. Buttery, M. Choi, L.L Lao, S. Smith, General Atomics H. St John, G. Staebler C. Kessel Princeton Plasma Physics Laboratory
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 USEU Transport Task Force Annual
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 informationIntermediate Nonlinear Development of a Linetied gmode
Intermediate Nonlinear Development of a Linetied gmode Ping Zhu University of WisconsinMadison In collaboration with C. C. Hegna and C. R. Sovinec (UWMadison) A. Bhattacharjee and K. Germaschewski
More informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (singlefluid)
More informationBounceaveraged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas
Bounceaveraged gyrokinetic simulations of trapped electron turbulence in elongated tokamak plasmas Lei Qi a, Jaemin Kwon a, T. S. Hahm a,b and Sumin Yi a a National Fusion Research Institute (NFRI), Daejeon,
More informationQTYUIOP LOCAL ANALYSIS OF CONFINEMENT AND TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL.
LOCAL ANALYSIS OF CONFINEMENT AND TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR Presented by D.P. SCHISSEL for the DIII D Team* Presented to 16th IAEA Fusion Conference
More informationParticleincell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas
Institute of Physics Publishing Journal of Physics: Conference Series 16 (25 16 24 doi:1.188/17426596/16/1/2 SciDAC 25 Particleincell simulations of electron transport from plasma turbulence: recent
More informationShear Flow Generation in Stellarators  Configurational Variations
Shear Flow Generation in Stellarators  Configurational Variations D. A. Spong 1), A. S. Ware 2), S. P. Hirshman 1), J. H. Harris 1), L. A. Berry 1) 1) Oak Ridge National Laboratory, Oak Ridge, Tennessee
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, selfconsistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More informationComputational modeling of fully ionized magnetized plasmas using the fluid approximation a
PHYSICS OF PLASMAS 13, 058103 2006 Computational modeling of fully ionized magnetized plasmas using the fluid approximation a D. D. Schnack b Center for Energy Space Science, Science Applications International
More informationImpact of Toroidal Flow on ITB HMode Plasma Performance in Fusion Tokamak
Impact of oroidal Flow on I HMode Plasma Performance in Fusion okamak oonyarit Chatthong 1,*, hawatchai Onjun 1, Roppon Picha and Nopporn Poolyarat 3 1 School of Manufacturing Systems and Mechanical Engineering,
More informationTwo Fluid Dynamo and EdgeResonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and EdgeResonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of WisconsinMadison, Madison,
More informationOverview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations
Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations Lecture by: P.H. Diamond Notes by: C.J. Lee March 19, 2014 Abstract Toroidal rotation is a key part of the design of ITER
More informationPredicting 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 informationExtended MHD simulation of RayleighTaylor/KelvinHelmholtz instability
Extended MHD simulation of RayleighTaylor/KelvinHelmholtz instability R. Goto (a), T. Hatori (a), H. Miura (a, b), A. Ito (a, b) and M. Sato (b) (a) The Graduate University for Advanced Studies (Sokendai)
More informationNonlinear MHD Simulations of Edge Localized Modes in ASDEX Upgrade. Matthias Hölzl, Isabel Krebs, Karl Lackner, Sibylle Günter
Nonlinear MHD Simulations of Edge Localized Modes in ASDEX Upgrade Matthias Hölzl, Isabel Krebs, Karl Lackner, Sibylle Günter Matthias Hölzl Nonlinear ELM Simulations DPG Spring Meeting, Jena, 02/2013
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 informationStationary, High Bootstrap Fraction Plasmas in DIIID Without Inductive Current Control
Stationary, High Bootstrap Fraction Plasmas in DIIID Without Inductive Current Control P. A. Politzer, 1 A. W. Hyatt, 1 T. C. Luce, 1 F. W. Perkins, 4 R. Prater, 1 A. D. Turnbull, 1 D. P. Brennan, 5 J.
More informationarxiv: v2 [physics.plasmph] 12 Sep 2015
Dynamics of Ion Temperature Gradient Turbulence and Transport with a Static Magnetic Island arxiv:156.7438v2 [physics.plasmph] 12 Sep 15 O. Izacard, 1, C. Holland, 2 S. D. James, 3 and D. P. Brennan 4
More informationGyrokinetic Transport Driven by Energetic Particle Modes
Gyrokinetic Transport Driven by Energetic Particle Modes by Eric Bass (General Atomics) Collaborators: Ron Waltz, Ming Chu GSEP Workshop General Atomics August 10, 2009 Outline I. Background Alfvén (TAE/EPM)
More informationExponential Growth of Nonlinear Ballooning Instability. Abstract
Exponential Growth of Nonlinear Ballooning Instability P. Zhu, C. C. Hegna, and C. R. Sovinec Center for Plasma Theory and Computation University of WisconsinMadison Madison, WI 53706, USA Abstract Recent
More informationCharacterization of Edge Stability and Ohmic Hmode in the PEGASUS Toroidal Experiment
Characterization of Edge Stability and Ohmic Hmode in the PEGASUS Toroidal Experiment M.W. Bongard, J.L. Barr, M.G. Burke, R.J. Fonck, E.T. Hinson, J.M. Perry, A.J. Redd, D.J. Schlossberg, K.E. Thome
More informationMHD Linear Stability Analysis Using a Full Wave Code
USJapan 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 informationTH/P614 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a)
1 TH/P614 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 Chang, C.S., 1 Ku, S., 2 Adams M., 3 D Azevedo, G., 4 Chen, Y., 5 Cummings,
More informationOperational Phase Space of the Edge Plasma in Alcator CMod
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 informationValidation of Theoretical Models of Intrinsic Torque in DIIID and Projection to ITER by Dimensionless Scaling
Validation of Theoretical Models of Intrinsic Torque in DIIID and Projection to ITER by Dimensionless Scaling by B.A. Grierson1, C. Chrystal2, W.X. Wang1, J.A. Boedo3, J.S. degrassie2, W.M. Solomon2,
More informationITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model
1 THC/33 ITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model J.E. Kinsey, G.M. Staebler, J. Candy, and R.E. Waltz General Atomics, P.O. Box 8608, San Diego, California
More informationR B. Here the first term represents
TH/61Rb Nonlinear MHD modelling of Edge Localized Modes and their interaction with Resonant Magnetic Perturbations in rotating plasmas. M.Bécoulet 1, F.Orain 1, J. Morales 1, X. Garbet 1, G. DifPradalier
More informationHybrid KineticMHD simulations with NIMROD
simulations with NIMROD 1 Yasushi Todo 2, Dylan P. Brennan 3, KwangIl You 4, JaeChun Seol 4 and the NIMROD Team 1 University of Washington, Seattle 2 NIFS, TokiJapan 3 University of Tulsa 4 NFRI, DaejeonKorea
More informationA mechanism for magnetic field stochastization and energy release during an edge pedestal collapse
A mechanism for magnetic field stochastization and energy release during an edge pedestal collapse S. S. Kim, Hogun Jhang, T. Rhee, G. Y. Park, R. Singh National Fusion Research Institute, Korea Acknowledgements:
More informationEdge Impurity Dynamics During an ELM Cycle in DIII D
Edge Impurity Dynamics During an ELM Cycle in by M.R. Wade 1 in collaboration with K.H. Burrell, A.W. Leonard, T.H. Osborne, P.B. Snyder, J.T. Hogan, 1 and D. Coster 3 1 Oak Ridge National Laboratory General
More informationTHE DIII D PROGRAM THREEYEAR PLAN
THE PROGRAM THREEYEAR PLAN by T.S. Taylor Presented to Program Advisory Committee Meeting January 2 21, 2 3 /TST/wj PURPOSE OF TALK Show that the program plan is appropriate to meet the goals and is wellaligned
More informationShear flows at the tokamak edge and their interaction with edgelocalized modes a
PHYSICS OF PLASMAS 14, 056118 2007 Shear flows at the tokamak edge and their interaction with edgelocalized modes a A. Y. Aydemir b Institute for Fusion Studies, The University of Texas at Austin, Austin,
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 DEFG0395ER54309 and Contract
More informationUCLA. Broadband Magnetic and Density Fluctuation Evolution Prior to First ELM in DIIID Edge Pedestal. Presented by G. Wang a, In collaboration with
Broadband Magnetic and Density Fluctuation Evolution Prior to First ELM in DIIID Edge Pedestal Presented by G. Wang a, In collaboration with W.A. Peebles a, P.B. Snyder b, T.L. Rhodes a, E.J. Doyle a,
More informationGTC Simulation of Turbulence and Transport in Tokamak Plasmas
GTC Simulation of Turbulence and Transport in Tokamak Plasmas Z. Lin University it of California, i Irvine, CA 92697, USA and GPSTTBP Team Supported by SciDAC GPSTTBP, GSEP & CPES Motivation Firstprinciples
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 informationGA A26565 LIMITS TO HMODE PEDESTAL PRESSURE GRADIENT IN DIIID
GA A26565 LIMITS TO HMODE PEDESTAL PRESSURE GRADIENT IN DIIID by R.J. GROEBNER, P.B. SNYDER, T.H. OSBORNE, A.W. LEONARD, T.L. RHODES, L. ZENG, E.A. UNTERBERG, Z. YAN, G.R. McKEE, C.J. LASNIER, J.A. BOEDO,
More informationEnergeticIonDriven MHD Instab. & Transport: Simulation Methods, V&V and Predictions
EnergeticIonDriven MHD Instab. & Transport: Simulation Methods, V&V and Predictions 7th APTWG Intl. Conference 58 June 2017 Nagoya Univ., Nagoya, Japan Andreas Bierwage, Yasushi Todo 14.1MeV 10 kev
More informationArbiTER studies of filamentary structures in the SOL of spherical tokamaks
ArbiTER studies of filamentary structures in the SOL of spherical tokamaks D. A. Baver, J. R. Myra, Research Corporation F. Scotti, Lawrence Livermore National Laboratory S. J. Zweben, Princeton Plasma
More informationResistive Wall Mode Control in DIIID
Resistive Wall Mode Control in DIIID by Andrea M. Garofalo 1 for G.L. Jackson 2, R.J. La Haye 2, M. Okabayashi 3, H. Reimerdes 1, E.J. Strait 2, R.J. Groebner 2, Y. In 4, M.J. Lanctot 1, G.A. Navratil
More informationCharacterization of neoclassical tearing modes in highperformance I mode plasmas with ICRF mode conversion flow drive on Alcator CMod
1 EX/P422 Characterization of neoclassical tearing modes in highperformance I mode plasmas with ICRF mode conversion flow drive on Alcator CMod Y. Lin, R.S. Granetz, A.E. Hubbard, M.L. Reinke, J.E.
More informationModelling plasma scenarios for MASTUpgrade
Modelling plasma scenarios for MASTUpgrade Neutral beam requirements, sensitivity studies and stability D. Keeling R. Akers, I. Chapman, G. Cunningham, H. Meyer, S. Pinches, S. Saarelma, O. Zolotukhin
More informationPossible Pedestal Transport Theory Models And Modeling Tests
Possible Pedestal Transport Theory Models And Modeling Tests J.D. Callen, University of Wisconsin, Madison, WI 537061609 DIIID Pedestal Transport Workshop, General Atomics, San Diego, CA, February 17,
More informationCMod Transport Program
CMod Transport Program PAC 2006 Presented by Martin Greenwald MIT Plasma Science & Fusion Center 1/26/2006 Introduction Programmatic Focus Transport is a broad topic so where do we focus? Where CMod
More informationInvestigation of Intrinsic Rotation Dependencies in Alcator CMod
Investigation of Intrinsic Rotation Dependencies in Alcator CMod D. Kwak, A. E. White, J. E. Rice, N. T. Howard, C. Gao, M. L. Reinke, M. Greenwald, C. Angioni, R. M. McDermott, and the CMod and ASDEX
More informationGA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELMSUPPRESSED REGIMES, AND ITER PREDICTIONS
GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELMSUPPRESSED REGIMES, AND ITER PREDICTIONS by P.B. SNYDER, T.H. OSBORNE, M.N.A. BEURSKENS, K.H. BURRELL, R.J. GROEBNER, J.W. HUGHES, R. MAINGI,
More informationCMod Core Transport Program. Presented by Martin Greenwald CMod PAC Feb. 68, 2008 MIT Plasma Science & Fusion Center
CMod Core Transport Program Presented by Martin Greenwald CMod PAC Feb. 68, 2008 MIT Plasma Science & Fusion Center Practical Motivations for Transport Research Overall plasma behavior must be robustly
More informationarxiv: v1 [physics.plasmph] 11 Mar 2016
1 Effect of magnetic perturbations on the 3D MHD selforganization of shaped tokamak plasmas arxiv:1603.03572v1 [physics.plasmph] 11 Mar 2016 D. Bonfiglio 1, S. Cappello 1, M. Veranda 1, L. Chacón 2 and
More informationInterlinkage of transports and its bridging mechanism
Interlinkage of transports and its bridging mechanism Katsumi Ida National Institute for Fusion Science 17 th International Toki Conference 1519 October 27, Toki OUTLINE 1 Introduction 2 particle pinch
More informationIntegrated Simulation of ELM Energy Loss and Cycle in Improved Hmode Plasmas
1 Integrated Simulation of ELM Energy Loss and Cycle in Improved Hmode Plasmas N. Hayashi 1), T. Takizuka 1), N. Aiba 1), N. Oyama 1), T. Ozeki 1), S. Wiesen 2), V. Parail 3) 1) Japan Atomic Energy Agency,
More informationCoarsegraining the electron distribution in turbulence simulations of tokamak plasmas
Coarsegraining the electron distribution in turbulence simulations of tokamak plasmas Yang Chen and Scott E. Parker University of Colorado at Boulder Gregory Rewoldt Princeton Plasma Physics Laboratory
More informationExperimental test of the neoclassical theory of poloidal rotation
Experimental test of the neoclassical theory of poloidal rotation Presented by Wayne Solomon with contributions from K.H. Burrell, R. Andre, L.R. Baylor, R. Budny, P. Gohil, R.J. Groebner, C.T. Holcomb,
More informationThe EPED Pedestal Model: Extensions, Application to ELMSuppressed Regimes, and ITER Predictions
The EPED Pedestal Model: Extensions, Application to ELMSuppressed Regimes, and ITER Predictions P.B. Snyder 1, T.H. Osborne 1, M.N.A. Beurskens 2, K.H. Burrell 1, R.J. Groebner 1, J.W. Hughes 3, R. Maingi
More informationFineScale Zonal Flow Suppression of Electron Temperature Gradient Turbulence
FineScale Zonal Flow Suppression of Electron Temperature Gradient Turbulence S.E. Parker, J.J. Kohut, Y. Chen, Z. Lin, F.L. Hinton and W.W. Lee Center for Integrated Plasma Studies, University of Colorado,
More informationUnderstanding and Predicting Profile Structure and Parametric Scaling of Intrinsic Rotation. Abstract
Understanding and Predicting Profile Structure and Parametric Scaling of Intrinsic Rotation W. X. Wang, B. A. Grierson, S. Ethier, J. Chen, and E. Startsev Plasma Physics Laboratory, Princeton University,
More information