170 likes | 187 Views
Collaborative project led by Guido Huysmans and Lynton Appel to develop self-consistent codes for equilibrium reconstruction and MHD stability analysis in plasma physics research.
E N D
IMP#1:Equilibrium and Linear MHD Stability Project leader: Guido Huysmans (CEA) Deputy: Lynton Appel (UKAEA)
Integrated Modelling Project 1 • Objective: • To provide an integrated suite of self-consistent codes (modules) for equilibrium reconstruction and linear MHD stability analysis • Topic 1A : Experimental Equilibrium reconstruction • CEDRES, CLISTE, EFIT, EQUINOX • Topic 1B : Equilibrium codes and linear MHD stability • Equilibrium : CAXE, CHEASE, DIVA, HELENA, VMEC, DINA • Mapping : COTRANS, JMC • MHD Stability : CAS3D, CASTOR, KINX, MISHKA, TERPSICHORE • Also: • Equilibrium toolbox : FLUSH
magnetics MSE equilibrium reconstruction code spec.parameters diagnostic(1)description machine description diagnostic(2)description high resolution equilibrium MHD output description equilibrium description equilibrium description code spec.parameters code spec.parameters MHD stability 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 • 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.
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
Task: ITM-05-IMP1-T1 • 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)
Task: ITM-05-IMP1-T2 • 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 McCarthy • Status : recent work on CLISTE will facilitate move to machine independence • CEDRES : • Persons : Patrick Hertout • Status : work on CEDRES about to be started
Task: ITM-05-IMP1-T3 • 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
Task: ITM-05-IMP1-T4 • Validation of the chain : equilibrium reconstruction -> high resolution equilibrium -> MHD stabilityon 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)
Task: ITM-05-IMP1-T5 • Provide a toolbox for the equilibrium operations (like JET-flush 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
Task: ITM-05-IMP1-T6 • 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
Task: ITM-05-IMP1-T7 • 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 McCarthy • Status: • working prototype (EFIT_ITM, in MATLAB), completely machine independent (Wolfgang Zwingmann) • EFITf90 (Lynton Appel) • CLISTE (Paddy McCarthy) to be combined with EFIT?
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
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/IMP1: 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 now 2006 2007
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.
Agenda IPM#1 session • Introduction : Guido Huysmans • equilibrium toolbox, introduction to flush : Huw Leggate • disussion on toolbox and requirements • EFIT_ITM : Wolfgang Zwingmann • EFITf90 : Lynton Appel • CLISTE: Paddy McCarthy • discussion on EFIT, how to get to one version? • discussion on compatibility framework /IMP1 requirements • progress on the IMP tasks • verification/validation of IMP#1 codes • Physics problem • volunteers
Questions • Storage of flush grid and interpolation results in equilibrium structure • Database requirements for simulation results