Call for Participation
2009 Work programme
INTEGRATED TOKAMAK MODELLING TASK FORCE
Integrated Modelling Project 1 (IMP1):
Equilibrium reconstruction and Linear MHD
Deadline December 23, 2008
Project leadership:
Guido Huysmans, [email protected]
Christian Konz, [email protected] (deputy)
Task Force leadership:
Pär Strand, [email protected]
Lars-Göran Eriksson, [email protected] (deputy)
Rui Coelho, [email protected] (deputy)
This Call for Participation aims to establish the IMP1 work programme for 2009 under EFDA Art.5
Task Agreement format.
1
Contents
Contents ...................................................................................................................................... 2
Required Resources ...................................................................................................................... 2
Priority Support............................................................................................................................ 3
Implementation ........................................................................................................................... 3
Intellectual Property Rights Monitoring ........................................................................................ 4
Overview: Equilibrium reconstruction and Linear MHD ................................................................. 4
ITM-09-IMP1-T1: Maintenance and continuing development of the available IMP1 codes .............. 6
ITM-09-IMP1-T2: Free boundary equilibrium and position/shape feedback control ........................ 8
ITM-09-IMP1-T2a: Implementation of a free boundary equilibrium code in the Kepler environment 8
ITM-09-IMP1-T2b: Feedback Control for plasma position and shape control ................................... 9
ITM-09-IMP1-T3: Adaptation of equilibrium reconstruction codes (continuation) ........................ 11
ITM-09-IMP1-T4: Flow extension of equilibrium and MHD stability codes .................................... 12
ITM-09-IMP1-T5: Validation of the equilibrium reconstruction codes available in Kepler .............. 13
ITM-09-IMP1-T6: Definition of 3D data structures for the equilibrium and implementation in 3D
equilibrium codes ....................................................................................................................... 14
Required Resources
The Implementation of the IMP1 work programme for 2009 is estimated to require a minimum of 5
Ppy under Baseline support and 2 Ppy under Priority Support in order to be able to provide a
minimum level of project fulfilment of the deliverables and milestones. The overall resources under
priority support for the ITM Task Force are determined by the Work Programme, but the breakdown
between projects within ITM is left to be optimised by the Task Force Leader in order to satisfy
project requirements while taking into account the resources made available by Associations.
Furthermore the resources under baseline support shown in the Work Programme are meant to be
preliminary. For these reasons, the manpower requirements stated in the project descriptions should
2
be viewed as indicative. Tasks are generally collaborative in character or form parts of a larger
structure where timely delivery and coordination is of importance. Support and training on the ITM
framework and tools will be available and time for the participation shall be allocated in addition to
Task activities.
In the analysis of the answers to the Call, the ITM-TF will do its best to promote collaborative efforts
between associations as well as between individuals in the different tasks. This is to provide internal
peer review, quality of work and to develop a stronger sense of community ownership for the
developed tools. However, and in order to reduce fragmentation and related administrative burden,
if smaller tasks are shared between Associations, care should be taken that no individual’s total
contribution to a project becomes too small. It is suggested that, unless explicitly requested by an
Association and/or technically justified, individual commitments per project should be of at least
3pm. This individual commitment could be shared between different tasks within the same project.
Tasks and individual commitments aimed directly at the maintenance of software already installed in
the ITM-TF platform are excluded from this requirement.
Priority Support
The task:
-
ITM-09-IMP1-T2: Free boundary equilibrium and position/feedback control
is formed under Priority Support.
In addition, there are opportunities for work relating to implementation and development of data
structures, support for development of code interfaces and integration into the KEPLER work flow
platform under the INTEGRATION TEAM. The team will provide and train a set of experts from the
different Modelling project to provide high level support and expertise on application development
and code adaptation for the ITM as a whole.
Implementation
All software development is expected to be implemented on the ITM-TF gateway,
www.efda-itm.eu, under the provisions of the Gateway User Agreement (attached). The latter was
agreed to by the EFDA-SC in June 2008. The Gateway User Agreement details access and sharing
mechanisms for the software developments within ITM-TF. Furthermore, the ITM-TF will provide a
collaborative software development environment, based on Gforge, to support the development of
individual programs and at the same time ensure that Quality Assurance and traceability criteria for
the ITM project are adhered to.
3
Intellectual Property Rights Monitoring
In agreement with the EFDA-SC decision (June 2008), and as described in the
Gateway User Agreement, the ITM-TF will monitor IPR relating to contributed codes. In particular,
the ITM-TF will maintain a record of contributions made to any Software through collaborative Tasks
within the ITM-TF. All such modifications remain fully available to the contributing Associate
provided contributors are acknowledged though the principles stated in the ITM-TF license.
A Rights of access form (attached) is required for all codes being contributed. For any given code, this
document states the current list of contributors, the code's ITM-TF Responsible Officer and technical
reference(s) that should be used in publications involving the code.
Overview: Equilibrium reconstruction and Linear MHD
In 2008, a first full chain of codes from equilibrium reconstruction, high resolution equilibrium and
linear MHD stability have been implemented on the Kepler platform. A new task has been defined for
the maintenance and continued development of these contributed codes.
The validation of the contributed codes, originally foreseen for 2008, will start/continue in 2009
when a stable ITM platform and the data access routines will be available.
The tasks continuing from 2008 are the adaptation of the equilibrium reconstruction codes and the
development of the data structures for a 3D machine description and 3D equilibriums.
New tasks in 2009 concern the free boundary equilibrium codes, the plasma position and shape
controllers and the addition of toroidal/poloidal flow in the equilibrium and linear MHD stability
codes.
2008 Task
ITM-05-IMP1-T3: High resolution
equilibrium codes
Status in 2009
 ITM-09-IMP1-T1: Maintenance and continuing
development of the available IMP1 codes
ITM-05-IMP1-T6: Linear MHD Stability
codes
ITM-05-IMP1-T2: Adaptation of new
equilibrium reconstruction codes
ITM-05-IMP1-T5: ITM equilibrium
toolbox: ITM_FLUSH
ITM-07-IMP1-T1: Full domain predictive
 ITM-09-IMP1-T2: Free boundary equilibrium and
4
equilibrium codes
position/shape feedback control :
ITM-09-IMP1-T2a: Implementation of a free
boundary equilibrium code in the Kepler
environment [Priority Support]
ITM-05-IMP1-T2: Adaptation of new
equilibrium reconstruction codes
 ITM-09-IMP1-T3: Adaptation of equilibrium
reconstruction codes (continuation)
ITM-07-IMP1-T2: Validation and
verification of contributed equilibrium
reconstruction codes
 ITM-09-IMP1-T5: Validation of the equilibrium
reconstruction codes available in Kepler
ITM-07-IMP1-T3: Extension of data
structures to 3D equilibriums.
 ITM-09-IMP1-T6: Definition of 3D data structures
for the equilibrium and implementation in 3D
equilibrium codes
New Tasks in 2009
N
ITM-09-IMP1-T2b: Feedback Control for plasma
position and shape control [Priority Support]
N
ITM-09-IMP1-T4: Flow extension of equilibrium and
MHD stability codes
5
ITM-09-IMP1-T1: Maintenance and continuing development of the
available IMP1 codes
Implementation Scheme: Baseline Support
Description of work:
In 2008, one equilibrium reconstruction code, three high resolution equilibrium codes and two
linear MHD stability codes have been made available as actors within the Kepler environment. A
continuous effort is required to keep the codes up-to-date with respect to the requirements, the
definition of the data structures and the correction of bugs.
Implementation of new features such as additional diagnostics in the equilibrium reconstruction
codes, additional coordinate systems in the high-res equilibrium codes or new solvers in the linear
MHD stability codes also fall under this task.
The task comprises 4 parts corresponding to the applications
1a) ITM equilibrium toolbox (FLUSH)
1b) equilibrium reconstruction codes (EFIT_ITM, EFIT++)
1c) high-res equilibrium codes (CHEASE, HELENA, and CAXE)
1d) linear MHD stability codes (ILSA, KINX)
Deliverables:
-
up-to-date version of the code and documentation available under SVN on the gateway
-
definition of standardised test cases for code verification, publication of verification test on
the IMP1 web pages.
Manpower skills and needs:
Requested manpower/skills:
People with expertise in the contributed IMP1 codes, the IMP1 data structures and the
implementation in the Kepler environment are required. It is estimated that the total manpower
needed for the activities within this task is about 10pm.
Existing Commitments:
This is a new task but continues from several tasks on the adaptation of IMP1 codes to the ITM
requirements. Contributors in 2008 are encouraged to express their ongoing interest.
6
7
ITM-09-IMP1-T2: Free boundary equilibrium and position/shape
feedback control
Implementation Scheme: Priority Support
An important subject for the whole device modelling, one of the main objectives of the ITM
taskforce, are free-boundary equilibrium codes. The plasma profiles are evolved within the separatrix
using 1D transport codes but the plasma shape and the external poloidal flux is evolved quasistatically driven by the currents in the poloidal field coils and possibly the induced currents in the
conducting structures (task 2a). A programmed evolution of the plasma shape requires a feedback
scheme on the coil currents (task 2b).
Given the central importance of the free-boundary equilibrium code(s) with respect to the
community needs, the requirements on the next step developments of the European Transport
solver (IMP3) and phasing with IMP1 general activities, the development of the ITM free-boundary
equilibrium module including the position and shape feedback control will be undertaken with
Priority Support with a ceiling of 24 pm (2 ppy) for 2009.
ITM-09-IMP1-T2a: Implementation of a free boundary equilibrium code
in the Kepler environment
Description of work:
The adaptation of an existing (or newly developed) free boundary equilibrium to the standardised
ITM data structures describing the machine geometry, poloidal field coils and plasma profiles.
The requirements for the direct free boundary code are:
-
given the external PF coil and the plasma profiles for pressure and current (or equivalent),
solve the static Grad-Shafranov equation.
-
secondly, the iron core (and possibly other ferromagnetic structures) and induced currents
in the vessel and metallic structures will have to be included
A free boundary solver which calculates the currents in the poloidal field coils from a prescribed
plasma shape would also be a very valuable contribution.
Within the IMP1 project, support is available for the integration of code modules on the Kepler
platform.
8
Deliverables:
-
A free boundary equilibrium code using the ITM data structures as input and output in one
of the supported programming languages, available as source code on the subversion ITM
server.
o First version of direct free boundary code running under Kepler : 1-8-2009
o Free boundary code with full specification under Kepler :end 2009
-
Documentation of contributed codes available on the subversion server.
Manpower skills and needs:
Requested manpower/skills:
Experience with free boundary equilibrium codes and the ITM data structures.
Existing Commitments:
Task ITM-09-IMP1-T2a is, in part, a continuation of ITM-07-IMP1-T1. Contributors in 2008 are
encouraged to express their ongoing interest.
ITM-09-IMP1-T2b: Feedback Control for plasma position and shape
control
To be able to use the free boundary code(s) contributed under task 2a as part of the transport solver,
a feedback controller is required to control the position and the shape of the plasma.
Description of work:
Adaptation of existing or newly developed feedback controller(s) for the plasma position and shape
into the ITM data structures. As part of the work, the ITM data structures may have to be extended
to take the requirements from the feedback controllers into account.
Deliverables:
(1) Evaluation of the existing ITM data structures with respect to the requirements for feedback
control. If required, extension of the ITM data structures.

Indicative time schedule : 1-4-2009
(2) Modules implementing a feedback controller for the plasma position and shape using the ITM
data structures available on the ITM Kepler platform.

Indicative time schedule : end 2009
9
Manpower skills and needs:
Requested manpower/skills:
Experience with plasma position/shape feedback controllers in modelling of plasma evolution.
Knowledge of the ITM data structures.
Existing Commitments:
This is a new task.
10
ITM-09-IMP1-T3: Adaptation of equilibrium reconstruction codes
(continuation)
Implementation Scheme: Baseline Support
Description of work:
The adaptation of the equilibrium reconstruction codes, newly contributed in 2008, needs to be
continued and completed in 2009.
Deliverables:
-
Adaptation of the equilibrium reconstruction codes to the current ITM data structures.
(Indicative schedule : 1-6-2009)
-
source code and documentation available on the ITM subversion server
-
codes running on the Kepler platform on the ITM gateway. (Indicative schedule : 1-9-2009)
-
verification of standard test cases, results published on the IMP1 website. (Indicative
schedule : end 2009)
Manpower skills and needs:
Requested manpower/skills:
Experience with equilibrium reconstruction codes and the ITM data structures. This activity is
estimated to require about 21 pm to cover the continued the full range of remaining codes from
the 2008 Work Programme.
Existing Commitments:
Task ITM-09-IMP1-T3 is a continuation of ITM-05-IMP1-T2. Contributors in 2008 are encouraged to
express their ongoing interest.
11
ITM-09-IMP1-T4: Flow extension of equilibrium and MHD stability
codes
Implementation Scheme: Baseline Support
Description of work:
The existing ITM data structures describing equilibrium and MHD stability will need to be extended
to include toroidal and possibly poloidal flow.
Adaptation of existing equilibrium with flow and linear MHD stability codes to the extended ITM
data structures.
Deliverables:
-
Extended ITM data structures for the equilibrium and MHD stability including toroidal and
possibly poloidal flow. ( Indicative schedule : 1-4-2009)
-
Adaptation of existing equilibrium codes and/or linear MHD stability codes with toroidal
and/or poloidal flow to the ITM data structures. (Indicative schedule: 1-9-2009)
-
Source code and documentation available on the ITM subversion server (Indicative
schedule: end 2009)
-
Publication of standard test cases on the IMP1 website. (Indicative schedule: end 2009)
Manpower skills and needs:
Requested manpower/skills:
Experience with equilibrium codes and/or linear MHD stability codes with flow and the ITM data
structures. This activity is estimated to require about 15 pm and cover several codes.
Existing Commitments:
This is a new task.
12
ITM-09-IMP1-T5: Validation of the equilibrium reconstruction codes
available in Kepler
Implementation Scheme: Baseline Support
Description of work:
Equilibrium reconstruction codes have been adapted to the ITM requirements and contributed to the
ITM taskforce. The next step is the verification and validation of the equilibrium reconstruction codes
following the V&V procedures as defined within the taskforce. The initial V&V is envisaged to use JET
and AUG data, contributions from other tokamaks are actively encouraged. The V&V exercise
involves the definition of the V&V procedure specific for equilibrium reconstruction, the collection of
suitable data into the ITM database and the use of the equilibrium reconstruction codes within the
ITM framework (code platform).
The end result will be a set of validated and documented codes that can be used on any tokamak for
which a machine description is available
Deliverables:
-
Collection of discharges to be used in the V&V from JET and ASDEX Upgrade available in ITM
database. (Indicative schedule: 1-4-2009)
-
Documentation of the results of the V&V (Indicative schedule: first V&V results , 1-9-2009 &
Completed validation: end 2009)
The deliverables critically depend on the access to JET and AUG equilibrium data. It also depends on
the availability of the machine descriptions from the EDRG and the tools to access the data from JET
and AUG.
Manpower skills and needs:
Requested manpower/skills:
Persons with experience in the reconstruction of equilibriums (but preferably not the original code
authors). This activity is estimated to require about 4 pm in relation to each of the participating
devices.
Existing Commitments:
This task is a continuation of task ITM-07-IMP1-T2. Contributors in 2008 are encouraged to express
their ongoing interest.
13
ITM-09-IMP1-T6: Definition of 3D data structures for the equilibrium
and implementation in 3D equilibrium codes
Implementation Scheme: Baseline Support
Description of work:
To prepare for the inclusion of stellarator codes and for the applications of 3D equilibriums on
tokamaks (like the effect of the ripple on the equilibrium reconstruction) the IMP1 data structures
will need to be extended to 3D. This includes both the machine description and the equilibrium data
structures.
The task on data structures for the 3D machine description is strongly linked with the EDRG that
have a coordinating role with other IMPs working on different aspects of the 3D machine
description (ITM-09-TFL2-EDRG-T2).
Adaptation of existing 3D equilibrium code(s) to the newly defined 3D equilibrium data structures.
Deliverables:
-
Extension of the data structures to 3D machine geometries and equilibriums. (Indicative
schedule: 1-6-2009)
-
3D equilibrium code adapted to the 3D data structures running under Kepler, source code
and documentation available on the ITM subversion server. (Indicative schedule : end 2009)
Manpower skills and needs:
Requested manpower/skills:
Expertise in 3D equilibrium codes is requested. Experience in ITM data structures would be very
useful but support for the ITM data structures is available through the Integration Team. This task is
estimated to require about 8 pm (2 pm for the data structure work).
Existing Commitments:
This is a continuation of task ITM-07-IMP1-T3. People who have expressed their interest in 2008 are
encouraged to confirm their commitment in 2009.
14