Status of TRANSP and PTRANSP Presented at APS-DPP

Status of TRANSP and PTRANSP Presented at APS-DPP 2009, Nov. 2 --6, 2009 *Supported by U. S. DOE Contract No. DE-AC 02 -09 CH 11466 Nov. 2, 2009 D. Mc. Cune 1

Abstract – GP 8. 00051 D. Mc. Cune, R. Andre, E. Feibush, M. Gorelenkova, K. Indireshkumar, C. Ludescher-Furth, L. Randerson, PPPL*, G. Bateman, A. Kritz, Lehigh University – This poster describes the status of TRANSP and PTRANSP code development and run production operations. Production rates continue to climb as new users and tokamaks are added; statistics will be shown, including utilization of the recently added TRANSP MPI capability. New code features include greatly expanded flexibility in specification of transport models for density, temperature, and angular momentum profile prediction in PTRANSP runs. A new replay option enables retrieval of sources from TRANSP analysis runs for PTRANSP validation. TRANSP / PTRANSP free boundary MHD modeling options have been enhanced. New RF modeling capabilities, such as the ability to model ICRF with full toroidal mode spectrum, have been added. NUBEAM upgrades include improvements to deposition atomic physics and loss orbit distribution capture. Fusion Grid post-processing of TRANSP results has been improved using extraction of SWIM-Sci. DAC Plasma State files from TRANSP archives. Opportunities for PTRANSP Support of Sci. DAC and FSP will be presented. *PPPL work performed under auspices of DOE contract DE-AC 02 -76 CH 03073 Nov. 2, 2009 D. Mc. Cune 2

TRANSP: Vision Statement Provide a comprehensive end -to-end modeling capability for magnetic confinement fusion energy experiments of today and tomorrow. Nov. 2, 2009 D. Mc. Cune 3

Traditional TRANSP: Overview Experiments (Asdex-U, C-Mod, DIII-D, ITER, JET, KSTAR, MAST, NSTX) Preliminary data Analysis and Preparation (largely automated) Diagnostic Hardware Pre- and Post-processing at the experimental site… Visualization Load Relational Databases MDS+ Experiment simulation Output Database ~1000 -2000 signals {f(t), f(x, t)} Detailed (3 d) time-slice physics simulations: GS 2, ORBIT, M 3 D… Nov. 2, 2009 D. Mc. Cune MDS+ 20 -50 signals {f(t), f(x, t)} Plasma position, Shape, Temperatures, Densities Field, Current, RF and Beam Injected Powers. TRANSP Analysis*: Current diffusion, MHD equilibrium, fast ions, heating, current drive; power, particle and momentum balance. *Fusion. Grid TRANSP on PPPL servers 4

PPPL TRANSP Run Production* Fusion Grid TRANSP *Not shown in graphic: independent run production at JET site (approximately 1000 runs per year). **59 NUBEAM/TRANSP production runs with MPI, up to 16 p. Nov. 2, 2009 D. Mc. Cune 5

PPPL TRANSP SERVICE FY-2005 – FY-2009 (Five Years) 18003 Total Runs Nov. 2, 2009 D. Mc. Cune *~4000 additional JET runs on JET production server over same period. 6

FY-2009 PPPL TRANSP Team Name TRANSP PTRANSP Sci. DAC Other R. Andre Marina G. 70% 100% Kumar Mc. Cune Feibush Ludescher Randerson Total FTE 65% 20% 80% 60% 4. 2 30% 35% 20% 15% 10% 0. 75 0. 70 20% 65% 10% 40% 1. 35 Color code: Physics, Visualization, Engineering/ Operational Support Nov. 2, 2009 D. Mc. Cune 7

TRANSP Staffing Level • Staff levels still insufficient: – Cannot promptly meet short term code development requests of production users; – Long term development projects have lagged. • Main reason: production system support. – Trouble-shooting 100 s of crashed runs. – Detailed answers to 1000 s of user queries. • New hire in FY-2009: Marina Gorelenkova. – Financed by experimental projects: big help. Nov. 2, 2009 D. Mc. Cune 8

A Cautionary Note for FSP • Production support very labor intensive! • Immediate user needs take precedence. – FTEs for long term development = Total FTEs less FTEs for support at best. – It’s worse than this, due to nature of support: • • Repeated interruptions; Numerous short term, urgent, unplanned tasks; User’s research dependent on speedy completion; Highly disruptive to long term planning and effort. • [P]TRANSP code development suffers. Nov. 2, 2009 D. Mc. Cune 9

Separate Production Support and Code Development Roles? • Might work for a sufficiently large team. • But, support personnel need knowledge similar to that of code authors. • Trouble shooting crashed runs = response to user “beta testing” of research code. • Integral part of “Validation & Verification”. • Code developer involvement is important. Nov. 2, 2009 D. Mc. Cune 10

A Day in the Life: Oct. 20, 2009 (zero new lines of code) • Detailed email exchange: • More emails: – JET physicist query, momentum balance terms. – TRANSP extracted Plasma States for GYRO @ GA. – Ideas for NUBEAM simulation of mode induced diffusion of fast ions. – FMCFM module design for PTRANSP, FSP, etc. – NUBEAM development planning (ADAS commit). – SWIM progress report for FY-2009 for DOE. Nov. 2, 2009 D. Mc. Cune – NUBEAM failure in MSTRFP simulations. – Diagnosis of TRANSP output system failure, 160 GB ITER file (output time step too small). – NSTX physicist query, fast ion trapping definition. – Help sought, debug DIII-D “Coppi-Tang” predictions. – Cursory review, two coauthored papers. – Join FSP “frameworks” email list. 11

New for Users This Year • Fusion Grid NUBEAM, MPI up to 16 p: – ADAS atomic physics: state excitation effects on neutral beam deposition. – New current shielding option: Y. R. Lin-Liu, F. L. Hinton, Phys. Plasmas 4 (1997) 4179. – New outputs @ times of interest OUTTIM(…): • • Both GC & FLR computed slowing down ion fb(…); Detailed deposition distribution function data; Distribution data on lost orbits. All available using GET_FBM. • Same NUBEAM in Sci. DACs: CSWIM, FACETs – MPI Scaling to ~1024 p, 1 -4000 ptcls/cpu. Nov. 2, 2009 D. Mc. Cune 12

New for Users This Year • NUBEAM_Fgen – see C. Stark, JP 8. 0060: – Replay NUBEAM calculation over short time window; produce data file for GET_FBM. – El. Vis (Fusion Grid) User Interface; – Start from any existing NBI TRANSP run, any time of interest; replay using Plasma States… – Choose particle statistics (NPTCLS); – Choose number of CPUs up to 16; – SULI summer project. • PTRANSP Free Boundary MHD equilibrium – Details below… Nov. 2, 2009 D. Mc. Cune 13

New for Users This Year • Plasma State based heating and current drive source replay: – TRANSP rerun using sources from prior run. – Sample application (PTRANSP): test transport model variations, holding sources fixed. • “Steerable” namelist quantities. • Two production versions: tshare & pshare: – Latest code in tshare; tested version in pshare. • Prompt input data error checking – At job queuing time; save queue waits. Nov. 2, 2009 D. Mc. Cune 14

Technical Progress (Code Not Ready for Users) • New RF-related Modules: – MPI-TORIC, GENRAY, CQL 3 D. – See GP 8. 00052 (next door). – TORIC Toroidal Mode Spectrum is available. • PTRANSP – Density prediction: tested, some aspects OK, still needs work. – Control Upgrades (being implemented): • Choice of solver algorithm • Choice of summed transport models • Axial, Confinement & Edge plasma regions… Nov. 2, 2009 D. Mc. Cune 15

Details: ADAS in NUBEAM Collaboration: PPPL & IPP/Garching • ADAS atomic physics – correction for atomic state excitation – now available in NUBEAM for H-isotope neutral beams • Significant change from ground state for high target densities (shown). • Improvement over legacy excitation model Janev, C. D. Boley, D. E. Post, Nucl. Fusion Vol. 29, No. 12, pp. 2125 -2140 (1989) Nov. 2, 2009 Asdex-U 17847 n(0) ~ 1. 2 e 20 m**-3 NBI source #1 Figure from “ADAS in Fafner and TRANSP/ NUBEAM” – Michael Kraus et al, October 2009 ADAS Workshop. D. Mc. Cune 16

Details: Isolver Free Boundary Solver • • Isolver is a free boundary equilibrium solver originally written by J. Menard in IDL and rewritten in Fortran for use in TRANSP/PTRANSP. Development of free boundary capabilities proceeding in stages: – EFIT analysis: Completed • Isolver analyzes psi(R, Z) from EFIT time slices used as input to TRANSP to find flux surface representation, separatrix, X points, … – Least squares coil currents (original Isolver mode): Completed • Isolver determines the poloidal field coil currents which will best match a prescribed boundary. – Measured coil currents: 70% Complete • Isolver uses the measured poloidal field coil currents along with a feedback circuit to keep the plasma centered in the vessel. – Passive structures: 30% Complete • The measured poloidal field coil currents are used along with a circuit equation to include the currents induced in the nearby passive structures. – Flux Diffusion: To Be Done • The poloidal field diffusion equation, which is used to evolve the plasma q, is coupled to the equilibrium solver. • The poloidal field coil currents are driven through a circuit equation to best match the measured coil currents. – Plasma current prediction: To Be Done • The plasma current is predicted from the circuit equation. • The drive voltage of the ohmic coils is chosen to match a prescribed plasma current waveform. Nov. 2, 2009 D. Mc. Cune 17

Isolver: EFIT Analysis • EFIT psi(R, Z, t) slice data is input to TRANSP/PTRANSP – LEVGEO=8 mode in TRANSP/PTRANSP for prescribed equilibrium data. – New NLPSIRZ_FLUX=. TRUE. mode uses psi(R, Z, t) input to evaluate flux surfaces. • Isolver determines flux surface based on location of X points and limiter surface. – Flux surface representation returned to TRANSP/PTRANSP is displaced in from separatrix to avoid singularities. – Separate analysis is performed with boundary at separatrix to evaluate enclosed flux and self inductance Li. • • Available for any tokamak with EFIT and limiter data. For application to D 3 D, see Budny et al. APS poster JP 8. 00102, Tuesday afternoon session. Nov. 2, 2009 D. Mc. Cune 18

Isolver: Least Squares Mode • • • Fit coil currents to best match prescribed boundary. Coil currents can be adjusted freely, fixed to measured data or constrained by measured data. Coil geometry currently available for NSTX and D 3 D. – D 3 D example in Budny, et al. JP 8. 00102 poster. • ITER coil geometry awaits TRANSP/PTRANSP output pipeline upgrade to accommodate large amount of psi(R, Z, t) data. NSTX Note: Ohmic coil current KK_OH left unconstrained and became nonzero only during the part of the shot where Isolver had difficulty finding a solution. Nov. 2, 2009 D. Mc. Cune 19

Isolver: Measured Coil Currents Mode • • Measured poloidal coil currents drive Isolver coil currents. Feedback circuits position the plasma. DIII-D – Stabilizes plasma against vertical and horizontal motion. – Feedback circuit currents are added to measured coil currents. – Prescribed boundary is used to identify feedback points for centering the plasma. Vertical feedback point Vertical feedback current source Horizontal feedback point Nov. 2, 2009 Horizontal feedback current source D. Mc. Cune 20

Details: Updatable Namelist Variables Namelist variables designated with tilda (“~”) in TRANSP’s port. spec are “updatable”. For example: beam and RF powers. Hypothetical application: [main namelist] …… nantech=1 ! 1 gyrotron nbeam=2 ! 2 beams powech=0. 0 pinja=2*0. 0 Time NBI ECH Total 80 -130 34 40 74 130 -300 17 30 47 300 -400 17 15 32 400 -500 0 4 4 ~update_time = 80. 0 powech=40. 0 e 6 pinja=17. 0, 17. 0 ~update_time = 130. 0 powech=30. 0 pinja=0. 0, 17. 0 … … [etc. ] Nov. 2, 2009 D. Mc. Cune 21

Details: Running MPI NUBEAM • In Namelist, set: – NBI_PSERVE = 1 • Command line job setup as before – tr_start [tshare] – It will ask you how many CPUs (Np)? – Note: Np > 16 likely to result in long queue wait. • In the near future: – NTORIC_PSERVE = 1 will be available. Nov. 2, 2009 D. Mc. Cune 22

El. Vis/Java Fusion Grid Services • Available Services*: • Recent Changes: – Nubeam_Fgen, get fast ion distrib. fcn. – Get_fbm, view fast ion distrib. fcn. – Cstate, view Plasma State file. – Rplot, view TRANSP run results. – El. Fresco, compute reflectometry *partial list Nov. 2, 2009 D. Mc. Cune – Nubeam_Fgen, Cstate, and Get_fbm are new. – Rplot: • MDS+ run selector panel; • Script editor; • Arrow Key command line recall. – Heavy use of El. Fresco service. 23

El. Vis/Rplot Nov. 2, 2009 D. Mc. Cune 24

I-phone Has an App for that… TRANSP Fusion Grid Users can access these services… Nov. 2, 2009 D. Mc. Cune 25

Future Plans • • • Continue Run Production & Support. Develop RF Modules with MPI. Further Develop MPI NUBEAM. Develop Free Boundary MHD Equilibrium. Develop GUIs– El. Vis/Java Services. PTRANSP / FSP / Sci. DACs: – PTRANSP winding down, much work needed. – Level of effort going forward is unclear. – Note: ITER task agreements need PTRANSP. Nov. 2, 2009 D. Mc. Cune 26

1 d Transport Solver Compatible with PTRANSP, FSP, Sci. DACs… P 1: Plasma State @t 1 • new geometry • new input profiles • missing profiles => equations to solve • new boundary conditions • sources (heating, fueling). Nov. 2, 2009 D. Mc. Cune FMCFM P 0: Plasma State @t 0 • Past time step geometry • Past time step profiles Portable Solver Module Control Object(s): • Choice of Solver • Choice of Transport Models • Options Current Status: PTRANSP generates P 0 and P 1 New profiles Collection of transport models: • NCLASS • GLF 23 • MMM 95 • MMM 08 • TGLF • etc. , etc. Question: Shall we build the rest of this? 27

Opportunity to Take Advantage of TRANSP Legacy Plasma State Time Series from Sci. DAC/FSP Simulation Includes Machine Description TRANSP Namelist Time Dependent Equilibrium & Profiles & Sources TRANSP Analysis TRANSP Output Dataset Question: Shall we build the rest of this? Nov. 2, 2009 Status: Replay Mode demonstrated transfer of sources from Plasma State time series. Access to Physics Modules that start with TRANSP data, including many first principles codes. D. Mc. Cune 28

Summary • Continued heavy research production use of TRANSP and PTRANSP. • Continued Active Development Program – But some developments are “behind schedule”. • PTRANSP funding & work coming to end. – Some goals not met, but, continuation of PTRANSP activity could help FSP & Sci. DACS. – ITER task agreements depend on PTRANSP. • The current TRANSP staff is stretched thin, mainly due to production support. Nov. 2, 2009 D. Mc. Cune 29