Disruption metadata in the summary IDS

[jorekart]

Dear DD developers,

For integrating the disruption database in SimDB, @olivhoenen suggested to include some metadata now used in the ITER_DISRUPTIONS database, in the summary.disruption. I propose to add the following fields within summary.disruption that may be useful to classify disruptions

• disruption.time (better definition) = min(t_TQ, t_VDE), where t_TQ and t_VDE are the times at which the Thermal Quench (TQ) and the Vertical Displacement Event (VDE) are detected by their respective detectors.

• disruption.time_predisruption = disruption.time - t_pre, where t_pre [s] = 0.0025 * R0[m]^2 and R0 is the major radius where the toroidal field is defined.

• runaway_electrons.current_max = Maximum toroidal current carried by runaway electrons [A]

• runaway_electrons.termination_type = [Benign, non-benign, other]

• halo_current.fraction_poloidal_max = Maximum poloidal halo current normalized by the pre-disruption toroidal plasma current []

• halo_current.fraction_toroidal_max = Maximum toroidal halo current normalized by the pre-disruption toroidal plasma current []

• halo_current.poloidal_tpf = Toroidal peaking factor of the poloidal halo current at the time of maximum poloidal halo current []

• vessel_net_force_max: Maximum of the net force in the vacuum vessel induced by the disruption [m].

  • z-direction
  • x-y plane: Maximum displacement in the x-y plane defined as max(sqrt(Fx(t)^2+Fy(t)^2))

• vessel_net_displacement_max: Maximum of the net displacement of the vacuum vessel due to the electromagnetic forces induced by the disruption [m].

  • z-direction
  • x-y plane: Maximum displacement in the x-y plane defined as max(sqrt(dx(t)^2+dy(t)^2))

• type: [To select among these two options]

  • Major Disruption: Thermal quench precedes the vertical displacement
  • VDE: Vertical displacement precedes the thermal quench

• origin: main cause leading to the disruption [To select among these options]

  • Density limit or MARFE
  • Impurity accumulation in the plasma core
  • Impurity accumulation in the plasma edge
  • Low q_95 or external kink limit
  • Vertical Displacement Event (VDE)
  • Ideal beta limit or internal kink mode
  • Neoclassical Tearing Mode (NTM)
  • Resistive Wall Mode (RWM)
  • Resistive Wall Tearing Mode (RWTM)
  • Error field or low-density limit
  • Mode locking
  • Tearing mode
  • Other

• mitigation_valve: Remove!

• Mitigation

  • mitigation_used: 1 if any disruption mitigation action took place before the end of the current quench, otherwise 0.
  • massive_material_injection [number_of_injection_sequences]:
    • technique: [MGI, SPI, DSP or other…]
    • sequence_type: [pre-TQ staggered, pre-TQ mixed, post-TQ, RE impact mitigation, other]
    • species(n_species):
      • atoms
      • species identifiers (A, Z, etc)
    • pre_GRE_time: Time from the material arrival at the LCFS to the time of the Global Reconnection Event (or GRE)
      • value [s]
      • GRE_time_definition: [time at the top of the Ip-spike, time at the Ip-dip]

• time_radiated_power_max: Move into radiation!

• radiation:

  • fraction:
    • thermal: Radiated energy during the thermal quench normalized by the pre-disruption thermal energy []
    • magnetic: Radiated energy during the current quench normalized by the magnetic energy dissipated within the plasma. []
    • total: Total radiated energy during the disruption normalized by the sum of the pre-disruptive thermal energy and the magnetic energy dissipated within the plasma. []
  • tpf_max: Toroidal peaking factor of the radiation at maximum total radiation power during the disruption
  • time_radiated_power_max

• simulation_phases: describes which disruption phases are covered by the simulations (for simulation data only):

  • disruption_precursors: simulation of evolution plasma parameters and before the thermal quench (1 if included, otherwise 0)
  • thermal_quench: simulation of loss of thermal energy during thermal quench (1 if included, otherwise 0)
  • current_quench: simulation of loss of magnetic energy during current quench (1 if included, otherwise 0)
  • RE_termination: simulation of loss of a runaway electron beam (1 if included, otherwise 0)

In addition, I find useful to filter the database by modelling assumptions. Let's say I am not interested on 2D simulations of the thermal quench, but I want to look at 3D features. For that, we could change the summary.code. to include

• dimensionality: [0,1,2,3,4...] number of spatial coordinates used in the simulation

[jorekart]

@imbeauf please note that @olivhoenen needs the quantities in bold to be implemented in fast track for the metadata. These are "pre_disruptive_values", "type", "runaway_electrons.current_max " and "halo_current.fraction_poloidal_max"

[imbeauf]

I have implemented a first version of the extension, with the quantities in bold.
Are the other quantities ready to be implemented, or do you expect input from others on this discussion ?
Should I create a limited PR right now ?
Please note that there was already a runaways structure in the Summary IDS, independent of the disruption structure, since you have have runaways independently of a disruption. I am thus adding runaways quantities there.

[olivhoenen]

Thank you @imbeauf, I understand that some of the extra definitions still needs to reach consensus in the disruption community, so let's start with a PR on the first version with bold quantities for now.