IMP1 Equilibrium and Linear MHD Stability Project leader

IMP#1: Equilibrium and Linear MHD Stability Project leader: Guido Huysmans (CEA) Deputy: Lynton Appel (UKAEA) G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Integrated Modelling Project 1 • Objective: – To provide an integrated suite of self-consistent codes (modules) for equilibrium reconstruction and linear MHD stability analysis • Topic 1 A : Experimental Equilibrium reconstruction – CEDRES, CLISTE, EFIT, EQUINOX • Topic 1 B : Equilibrium codes and linear MHD stability – Equilibrium: CAXE, CHEASE, DIVA, HELENA, VMEC, DINA – Mapping : COTRANS, JMC – MHD Stability : CAS 3 D, CASTOR, KINX, MISHKA, TERPSICHORE • Also: – Equilibrium toolbox : FLUSH G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

IMP-1: equilibrium and MHD Stability • Standardise contributed codes to become independent of machine /diagnostic data. – Use only external geometry data (from database) – Definition of interfaces between codes and machine and diagnostics magnetics MSE diagnostic(1) description diagnostic(2) description machine description equilibrium reconstruction equilibrium description high resolution equilibrium • Validation and Verification – compare equilibrium and MHD stability codes on benchmark case – Apply codes to a relevant experimental problem/data • MHD Stability limits in plasmas with an internal transport barrier. G. Huysmans ITM workshop: CCP-DCP-IMP-1 code spec. parameters equilibrium description MHD stability code spec. parameters MHD output description 9 -10/11/2005

Standardized Interfaces • Example code interfaces (sketch): Call Equilibrium_Reconstruction(geometry, measurements, equilibrium_out, ERcode_specific_parameters) Call High_Res_Equilibrium(equilibrium_out, equilibrium_high_res, HREcode_specific_parameters) Call MHD_stability(equilibrium_high_res, MHD_stability_structure, MHDcode_specific_parameters) equilibrium_out, equilibrium_high_res : instances of equilibrium structure MHD_stability_structure : to be defined G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 1 • Definition of the data structures to describe a generic machine and its diagnostics. • Definition of the exchange methods of extended equilibrium data and MHD stability results. (collaboration with DCP) – Persons : Jo Lister (DCP), Lynton Appel, Guido Huysmans – Status : • Structures/Objects defined to describe machine, diagnostics and equilibria • Database in MDSplus • close to first frozen version (urgently required) G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 2 • adaptation of equilibrium reconstruction codes to be machine independent, using the data structures of the ITM taskforce. – EQUINOX : • Persons : responsible person left fusion, new person will be appointed • Status : – EFIT : • Persons : Lynton Appel, Wolfgang Zwingmann, Vladimir Drozdov • Status : several versions exist, work going on on multiple fronts – CLISTE : • Persons : Patrick Mc. Carthy • Status : recent work on CLISTE will facilitate move to machine independence – CEDRES : • Persons : Patrick Hertout • Status : work on CEDRES about to be started G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 3 • adaptation of high resolution equilibrium, mapping and linear MHD stability codes to the ITM data structures. – HELENA • First interface with MDSplus database ready using the fortran type definitions using simple euitmopen, euitmget, euitmput routines • Persons : Chistian Konz, Guido Huysmans • Status : in progress – CHEASE • Persons : Olivier Sauter, Hinrich Lutjens • Status : to be done – KINX • Persons : Olivier Sauter, Laurent Villard • Status : to be done – CASTOR, MISHKA • Persons : Christian Konz, Erika Strumberger, Guido Huysmans • Status : to be done G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 4 • Validation of the chain : equilibrium reconstruction -> high resolution equilibrium -> MHD stability on synthetic data and on data from several tokamaks on the MHD stability of internal transport barriers. – 2006 – Looking for persons to collect data (Filomena Nave) G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 5 • Provide a toolbox for the equilibrium operations (like JETflush etc) – – – Provide field magnetic flux and field in space Flux surface averaging Intersections line of sights Finding X-point, plasma boundary … (requirements from IMP’s have been requested) – FLUSH : Huw Leggate – Status : starting up, flush specific requirements to be added to equilibrium structure G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 6 • Create a unique version of the CASTOR and MISHKA MHD stability codes (combining the existing versions into one framework) – Include MISHKA into the CASTOR framework such that one can choose the physics model as an input parameter – Persons : Christian Konz, Erika Strumberger, Guido Huysmans – Status : soon to be started G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Task: ITM-05 -IMP 1 -T 7 • create a unique version of EFIT – Several EFIT versions exist typically adapted to a specific machine • JET, MAST, Tore Supra, … – Persons: Lynton Appel, Wolfgang Zwingmann, Vladimir Drozdov, Emilia Solano, Patrtick Mc. Carthy – Status: • working prototype (EFIT_ITM, in MATLAB), completely machine independent (Wolfgang Zwingmann) • EFITf 90 (Lynton Appel) • CLISTE (Paddy Mc. Carthy) to be combined with EFIT? G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

In addition • Tool(s) for the reconstruction of the plasma boundary – required by EQUINOX – Adapt XLOC to taskforce structures? • Direct free-boundary equilibrium codes (DINA, CEDRES, …) • Extended MHD Stability codes G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

IMP-1 timeline call for participation 17/12/2004 project leader nominated 1/2/2005 Code inventory 11/2/2005 unification of code versions 1/1/2006 DCP/IMP 1: definition of interfaces/data structures 30/9/2005 prototype database/codes 1/1/2006 verification codes 1/4/2006 adaptation codes 1/1/2006 machine/diagnostic data 1/1/2006 validation 1/1/2007 experimental data 1/4/2006 2005 G. Huysmans now ITM workshop: CCP-DCP-IMP-1 2006 2007 9 -10/11/2005

IMP#1 Wishes/Questions • Finalise the schemas definitions – Decide where to put time – How about code specific parameters – How to handle multiple versions (like equilibria made by different codes for the same shot/time) • general solution required (ad hoc for equilibria at the moment) • How to get data into the database? – At least one complete case that can be used as a test case (experimental data and synthetic benchmark) – How to store simulation results? (i. e. results from the whole chain from reconstruction to MHD stability) • IMP#1 request from UAL: – Possibility to read/write structures (or objects) • read and write whole subtrees • Call euitmget(‘equilibrium’, machine, shot, time, an_equilibrium) – Option of read/write structures to file. G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Agenda IPM#1 session • Introduction : Guido Huysmans • equilibrium toolbox, introduction to flush : Huw Leggate • disussion on toolbox and requirements • EFIT_ITM : Wolfgang Zwingmann • EFITf 90 : Lynton Appel • CLISTE: Paddy Mc. Carthy • discussion on EFIT, how to get to one version? • discussion on compatibility framework /IMP 1 requirements • progress on the IMP tasks • verification/validation of IMP#1 codes – Physics problem – volunteers G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005

Questions • Storage of flush grid and interpolation results in equilibrium structure • Database requirements for simulation results G. Huysmans ITM workshop: CCP-DCP-IMP-1 9 -10/11/2005