Integration of Heating and Fast Particles Models

1. Integration of heating and fast particles models D. Farina, T. Johnson, G. Vlad and IMP5 contributorsthanks to L.-G. Eriksson, D. Kalupin, J. Ferreira, V. Baisuk, P. HuynhTF Leaders : G. Falchetto, R. Coelho, D. Coster EFDA CSU Contact Person: D. Kalupinhttps://www.efda-itm.eu/ITM/html/imp5_public.html 2011-09-13

2. Outline • Overview of the IMP5 activities • Structure of the integrated model • consistency / synergies • Connection between physics models and CPO • connections suggests a design of an integrated workflow • The integrated H&CD actor, IMP5HCD • merging heating from EC/IC/LH/NB/alpha heating • Coupling to the ETS • Strongly coupled model

3. IMP5: Heating, current drive and fast particles • IMP5 organizes the ITM work concerning heating and current drive and fast ions • 38 members (not all active...) , PS ~6.4 ppys, BS ~3.2 ppys • 33 codes (not all at the same status) • Heating schemes: • EC, LH, IC, NBI, and alphas heating • Fast particles: • interactions with MHD waves • run-away electrons, … • In 2010 the activity was focussed on the integration of codes • In 2011: • continue code integration; focus on testing in Kepler and documentation • integrated workflows and coupling to the ETS • 4.08b we had 5 codes running • porting to 4.09a delayed • special support: • synergy between heating schemes • self-consistent wave-particle interaction

4. Structure of integrated model • What are the fundamental modules of an integrated H&CD module? • Two considerations • allow strong coupling, e.g. between wave-field and kinetic solvers • suggests modules based on heating schemes • allow for synergies, e.g. electrons absorb both EC, LH and IC waves and NBI ions are accelerated by IC waves • suggests modules based on physical objects • Thus, we have choosen to modularise the integration into physical objects • we believe this choice will allow easier workflows with both strong coupling and synergies • treatment of strong coupling is discussed later top Physical objects: Wave models Source models Kinetic models Wave models Source models Kinetic models Flexible: EC LH IC NBI alpha ion electron wave heating run-aways

5. IMP5 datastructures (CPOs) • The datastructures follow the same principle; • they are Consistent Physical Objects (CPOs)not organised by the heating schemes • IMP5 CPOs • WAVES: representations for global waves and ray-tracing • DISTRIBUTION: fast/thermal ions, expanded distribution functions, gridded and Monte Carlo representations • DISTSOURCE: sources from NBI, fusion products (alphas) • Machine CPOs: ANTENNAS and NBI

6. Connection: CPOs / models - CPO - Physics model BEAM INJECTION ALPHAS Fusion reaction model Beam deposition model Machine CPOs: INPUT Self-consistency NBI: description of beam injector DISTSOURCE: source for kineticion/electronmodels Machinedescriptions DISTRIBUTION: distribution functions Kineticmodel(e.g. Fokker-Planck) + Experimental data RF HEATING ANTENNAS: boundarycondition at antenna WAVES: Wave field Self-consistency OUTPUT CORESOURCE: input CPO for ETS Wave model Datajoiners and mergers;

7. The composite actor IMP5HCD IMP5 has an actor, IMP5HCD, that combines all heating schemes • Each CPOs is filled by one composite actor svn: http://gforge.efda-itm.eu/svn/keplerworkflows/trunk/4.09a/imp5/imp5hcd/imp5hcd.xml website: https://www.efda-itm.eu/ITM/html/imp5_workflow__imp5hcd.html python graphs ETS output Physics modules CORESOURCE DISTRIBUTION WAVES DISTSOURCE NOTE: this is an actor – to be plugged into any workflow

8. Layered composite actors IMP5HCD has many layers of composite actorse.g. consider the Waves actor… Open composite actor The CPOs needed for GRAY are extracted from the ITM Plasma Bundle, which is a bundle of CPO, and control parameter and the time. see https://www.efda-itm.eu/ITM/html/itm_conventions.html#itm_conventions_20

9. Status of the IMP5HCD • The IMP5HCD is running • in standalone workflow • in the ETS (4.09a for Garching/Lisbon version / tested in 4.08b in Cadarache version) • coupling to other workflows coming… • e.g. integration into MHD-control workflow • In 4.08b we had 5 physics codes running • In 4.09a only one code works so far • Major changes to CPOs has delayed transition to 4.09a for IC and NBI codes • Areas where there are no actors (4.08b or 4.09a): • LH waves • electron Fokker-Planck • sources of fusion products we hope to have these actor during 2011

10. Coupling to the ETS • The IMP5HCD is simple to couple to the ETS • IMP5HCD is included as a single box in> ETS/CONVERGENCE_LOOP/UPDATE_SOURCES/HCD imp5hcd

11. Strongly coupled objects • The IMP5HCD actor is optimised for weakly coupled codes • i.e. slow time evolution / not stiff dynamics • Strongly coupled objects can be handled in local convergence loops (illustrated later)(excluding very strongly coupled models that should be coupled within a single physics module) • Inclusion of fast particle-MHD interaction • at present: fast-particle codes run outside workflow (post processing) • later: ”event models” can be triggered by stability constraints

12. Adopt IMP5HCD to strongly coupled modules • Example problem: coupling between IC wave solver (TORIC) and Fokker-Planck solver (FPsim) • Fast ions are accelerated by wave • the dielectric response depend on fast ions from DISTRIBUTION WAVES Check consistencywave/distribution FPsim TORIC Initial wave field with dielectric respons from distribution at last time step YES Fokker-Planck with wave input NO TORIC Update wave field with dielectric respons from new distribution

13. Fast ion-MHD: MARS-workflow • MARS (initial condition version) is the present field solver and time loop driver for Hybrid MHD-Gyrokinetic code HMAGYC • A standalone workflow has been developed to run MARS • e.g. post processing with data from ETS • To be extended to include the full HMAGYC code mars

14. Summary • IMP5 have developed an actor that merges all heating schemes • Primary application is the ETS • but it could be applied to any problem needing integrated HCD & fast particle modelling • Modularised according to the physics; in line with CPOs • allow treatment of synergies • strongly coupled problems could be included through local convergence loops inside the actor • The actor still lack physics modules for LH and alphas • IC and NBI have been delayed by 4.09a