Skip to content

current_drive_variables

Module containing global variables relating to the current drive system

dx_beam_duct = None module-attribute

width of neutral beam duct where it passes between the TF coils (m) T Inoue et al, Design of neutral beam system for ITER-FEAT, Fusion Engineering and Design, Volumes 56-57, October 2001, Pages 517-521)

big_q_plasma = None module-attribute

Fusion gain; P_fusion / (P_injection + P_ohmic)

f_c_plasma_bootstrap = None module-attribute

bootstrap current fraction (enforced; see i_bootstrap_current)

f_c_plasma_bootstrap_max = None module-attribute

maximum fraction of plasma current from bootstrap; if f_c_plasma_bootstrap_max < 0, bootstrap fraction = abs(f_c_plasma_bootstrap_max)

f_c_plasma_bootstrap_iter89 = None module-attribute

bootstrap current fraction, ITER 1989 model

f_c_plasma_bootstrap_nevins = None module-attribute

bootstrap current fraction, Nevins et al model

f_c_plasma_bootstrap_sauter = None module-attribute

bootstrap current fraction, Sauter et al model

f_c_plasma_bootstrap_wilson = None module-attribute

bootstrap current fraction, Wilson et al model

f_c_plasma_bootstrap_sakai = None module-attribute

Bootstrap current fraction, Sakai et al model

f_c_plasma_bootstrap_aries = None module-attribute

Bootstrap current fraction, ARIES model

f_c_plasma_bootstrap_andrade = None module-attribute

Bootstrap current fraction, Andrade et al model

f_c_plasma_bootstrap_hoang = None module-attribute

Bootstrap current fraction, Hoang et al model

f_c_plasma_bootstrap_wong = None module-attribute

Bootstrap current fraction, Wong et al model

bscf_gi_i = None module-attribute

Bootstrap current fraction, first Gi et al model

bscf_gi_ii = None module-attribute

Bootstrap current fraction, second Gi et al model

f_c_plasma_bootstrap_sugiyama_l = None module-attribute

Bootstrap current fraction, L-mode Sugiyama et al model

f_c_plasma_bootstrap_sugiyama_h = None module-attribute

Bootstrap current fraction, H-mode Sugiyama et al model

cboot = None module-attribute

bootstrap current fraction multiplier

c_beam_total = None module-attribute

neutral beam current (A)

f_c_plasma_diamagnetic_hender = None module-attribute

diamagnetic current fraction, Hender fit

f_c_plasma_diamagnetic_scene = None module-attribute

diamagnetic current fraction, SCENE fit

f_c_plasma_diamagnetic = None module-attribute

diamagnetic current fraction

p_hcd_ecrh_injected_total_mw = None module-attribute

ECH power (MW)

p_ebw_injected_mw = None module-attribute

Electron bernstein power (MW)

p_hcd_ecrh_electric_mw = None module-attribute

ECH wall plug power (MW)

p_hcd_ebw_electric_mw = None module-attribute

Electron bernstein wall plug power (MW)

eta_cd_hcd_primary = None module-attribute

Current drive efficiency of primary HCD system (A/W)

eta_cd_hcd_secondary = None module-attribute

Current drive efficiency of secondary HCD system (A/W)

c_hcd_primary_driven = None module-attribute

Current in plasma driven by primary HCD system (A)

c_hcd_secondary_driven = None module-attribute

Current in plasma driven by secondary HCD system (A)

f_c_plasma_hcd_primary = None module-attribute

Fraction of plasma current driven by primary HCD system

f_c_plasma_hcd_secondary = None module-attribute

Fraction of plasma current driven by secondary HCD system

n_ecrh_harmonic = None module-attribute

cyclotron harmonic frequency number, used in cut-off function

i_ecrh_wave_mode = None module-attribute

Switch for ECRH wave mode : - =0 O-mode - =1 X-mode

e_beam_kev = None module-attribute

neutral beam energy (keV) (iteration variable 19)

eta_hcd_primary_injector_wall_plug = None module-attribute

auxiliary power wall plug to injector efficiency

eta_hcd_secondary_injector_wall_plug = None module-attribute

secondary auxiliary power wall plug to injector efficiency

eta_ecrh_injector_wall_plug = None module-attribute

ECH wall plug to injector efficiency

eta_lowhyb_injector_wall_plug = None module-attribute

lower hybrid wall plug to injector efficiency

eta_icrh_injector_wall_plug = None module-attribute

Ion cyclotron wall plug to injector efficiency

eta_ebw_injector_wall_plug = None module-attribute

Electron bernstein wave wall plug to injector efficiency

eta_beam_injector_wall_plug = None module-attribute

neutral beam wall plug to injector efficiency

f_p_beam_injected_ions = None module-attribute

fraction of beam energy to ions

p_beam_injected_mw = None module-attribute

neutral beam power entering vacuum vessel

f_c_plasma_pfirsch_schluter_scene = None module-attribute

Pfirsch-Schlüter current fraction, SCENE fit

p_beam_shine_through_mw = None module-attribute

neutral beam shine-through power

feffcd = None module-attribute

current drive efficiency fudge factor (iteration variable 47)

f_p_beam_orbit_loss = None module-attribute

fraction of neutral beam power lost after ionisation but before thermalisation (orbit loss fraction)

f_radius_beam_tangency_rmajor = None module-attribute

R_tangential / R_major for neutral beam injection

f_beam_tritium = None module-attribute

fraction of beam that is tritium

eta_cd_norm_hcd_primary = None module-attribute

Normalised current drive efficiency for primary HCD system [(1.0e20 A)/(W m^2)]

eta_cd_dimensionless_hcd_primary = None module-attribute

Dimensionless current drive efficiency for primary HCD system (ζ)

eta_cd_norm_hcd_secondary = None module-attribute

Normalised current drive efficiency for secondary HCD system [(1.0e20 A)/(W m^2)]

eta_cd_dimensionless_hcd_secondary = None module-attribute

Dimensionless current drive efficiency for secondary HCD system (ζ)

eta_cd_norm_ecrh = None module-attribute

User input ECRH gamma (1.0e20 A/(W m^2))

xi_ebw = None module-attribute

User scaling input for EBW plasma heating. Default 0.43

i_hcd_primary = None module-attribute

Switch for current drive efficiency model: - =1 Fenstermacher Lower Hybrid - =2 Ion Cyclotron current drive - =3 Fenstermacher ECH - =4 Ehst Lower Hybrid - =5 ITER Neutral Beam - =6 new Culham Lower Hybrid model - =7 new Culham ECCD model - =8 new Culham Neutral Beam model - =9 RFP option removed in PROCESS (issue #508) - =10 ECRH user input gamma - =11 ECRH "HARE" model (E. Poli, Physics of Plasmas 2019). Removed in #1811. - =12 EBW user scaling input. Scaling (S. Freethy)

i_hcd_secondary = None module-attribute

Switch for 2nd current drive efficiency model: - =0 No fixed current drive - =1 Fenstermacher Lower Hybrid - =2 Ion Cyclotron current drive - =3 Fenstermacher ECH - =4 Ehst Lower Hybrid - =5 ITER Neutral Beam - =6 new Culham Lower Hybrid model - =7 new Culham ECCD model - =8 new Culham Neutral Beam model - =9 RFP option removed in PROCESS (issue #508) - =10 ECRH user input gamma - =11 ECRH "HARE" model (E. Poli, Physics of Plasmas 2019). Removed in #1811. - =12 EBW user scaling input. Scaling (S. Freethy)

i_hcd_calculations = None module-attribute

Switch for current drive calculation: - =0 turned off - =1 turned on

f_p_beam_shine_through = None module-attribute

neutral beam shine-through fraction

dx_beam_shield = None module-attribute

neutral beam duct shielding thickness (m)

p_hcd_primary_extra_heat_mw = None module-attribute

heating power not used for current drive (MW) (iteration variable 11)

p_hcd_secondary_extra_heat_mw = None module-attribute

secondary fixed heating power not used for current drive (MW)

p_hcd_injected_max = None module-attribute

maximum allowable value for injected power (MW) (constraint equation 30)

p_hcd_injected_electrons_mw = None module-attribute

auxiliary injected power to electrons (MW)

p_hcd_injected_ions_mw = None module-attribute

auxiliary injected power to ions (MW)

p_hcd_injected_total_mw = None module-attribute

total auxiliary injected power (MW)

p_hcd_injected_current_total_mw = None module-attribute

total auxiliary injected power (MW)

p_hcd_secondary_injected_mw = None module-attribute

secondary total fixed auxiliary injected power (MW)

p_hcd_primary_injected_mw = None module-attribute

primary auxiliary injected power (MW)

f_c_plasma_internal = None module-attribute

plasma current fraction driven internally (Bootstrap + Diamagnetic + PS)

p_hcd_lowhyb_injected_total_mw = None module-attribute

Total lower hybrid injection power (MW)

p_hcd_icrh_injected_total_mw = None module-attribute

Total ion cyclotron injection power (MW)

p_hcd_ebw_injected_total_mw = None module-attribute

Total electron bernstein wave injection power (MW)

p_beam_plasma_coupled_mw = None module-attribute

Total neutral beam power that is coupled to plasma after losses (MW)

p_hcd_beam_injected_total_mw = None module-attribute

neutral beam injection power (MW)

p_beam_orbit_loss_mw = None module-attribute

neutral beam power lost after ionisation but before thermalisation (orbit loss power) (MW)

f_c_plasma_pfirsch_schluter = None module-attribute

Pfirsch-Schlüter current fraction

p_hcd_lowhyb_electric_mw = None module-attribute

lower hybrid wall plug power (MW)

pwpnb = None module-attribute

neutral beam wall plug power (MW)

radius_beam_tangency = None module-attribute

neutral beam centreline tangency radius (m)

radius_beam_tangency_max = None module-attribute

maximum tangency radius for centreline of beam (m)

n_beam_decay_lengths_core = None module-attribute

neutral beam e-decay lengths to plasma centre

n_beam_decay_lengths_core_required = None module-attribute

permitted neutral beam e-decay lengths to plasma centre

init_current_drive_variables()

Initialise current drive variables

Source code in process/data_structure/current_drive_variables.py 
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
def init_current_drive_variables():
    """Initialise current drive variables"""
    global \
        dx_beam_duct, \
        big_q_plasma, \
        f_c_plasma_bootstrap, \
        f_c_plasma_bootstrap_max, \
        f_c_plasma_bootstrap_iter89, \
        f_c_plasma_bootstrap_nevins, \
        f_c_plasma_bootstrap_sauter, \
        f_c_plasma_bootstrap_wilson, \
        f_c_plasma_bootstrap_sakai, \
        f_c_plasma_bootstrap_aries, \
        f_c_plasma_bootstrap_andrade, \
        f_c_plasma_bootstrap_hoang, \
        f_c_plasma_bootstrap_wong, \
        bscf_gi_i, \
        bscf_gi_ii, \
        f_c_plasma_bootstrap_sugiyama_l, \
        f_c_plasma_bootstrap_sugiyama_h, \
        cboot, \
        c_beam_total, \
        f_c_plasma_diamagnetic_hender, \
        f_c_plasma_diamagnetic_scene, \
        f_c_plasma_diamagnetic, \
        p_hcd_ecrh_injected_total_mw, \
        p_ebw_injected_mw, \
        p_hcd_ecrh_electric_mw, \
        p_hcd_ebw_electric_mw, \
        eta_cd_hcd_primary, \
        eta_cd_hcd_secondary, \
        c_hcd_primary_driven, \
        c_hcd_secondary_driven, \
        f_c_plasma_hcd_primary, \
        f_c_plasma_hcd_secondary, \
        n_ecrh_harmonic, \
        i_ecrh_wave_mode, \
        e_beam_kev, \
        eta_hcd_primary_injector_wall_plug, \
        eta_hcd_secondary_injector_wall_plug, \
        eta_ecrh_injector_wall_plug, \
        eta_lowhyb_injector_wall_plug, \
        eta_icrh_injector_wall_plug, \
        eta_ebw_injector_wall_plug, \
        eta_beam_injector_wall_plug, \
        f_p_beam_injected_ions, \
        p_beam_injected_mw, \
        f_c_plasma_pfirsch_schluter_scene, \
        p_beam_shine_through_mw, \
        feffcd, \
        f_p_beam_orbit_loss, \
        f_radius_beam_tangency_rmajor, \
        f_beam_tritium, \
        eta_cd_norm_hcd_primary, \
        eta_cd_dimensionless_hcd_primary, \
        eta_cd_norm_hcd_secondary, \
        eta_cd_dimensionless_hcd_secondary, \
        eta_cd_norm_ecrh, \
        xi_ebw, \
        i_hcd_primary, \
        i_hcd_secondary, \
        i_hcd_calculations, \
        f_p_beam_shine_through, \
        dx_beam_shield, \
        p_hcd_primary_extra_heat_mw, \
        p_hcd_secondary_extra_heat_mw, \
        p_hcd_injected_max, \
        p_hcd_injected_electrons_mw, \
        p_hcd_injected_ions_mw, \
        p_hcd_injected_total_mw, \
        p_hcd_injected_current_total_mw, \
        p_hcd_secondary_injected_mw, \
        p_hcd_primary_injected_mw, \
        f_c_plasma_internal, \
        p_hcd_lowhyb_injected_total_mw, \
        p_hcd_icrh_injected_total_mw, \
        p_hcd_ebw_injected_total_mw, \
        p_beam_plasma_coupled_mw, \
        p_hcd_beam_injected_total_mw, \
        p_beam_orbit_loss_mw, \
        f_c_plasma_pfirsch_schluter, \
        p_hcd_lowhyb_electric_mw, \
        pwpnb, \
        radius_beam_tangency, \
        radius_beam_tangency_max, \
        n_beam_decay_lengths_core, \
        n_beam_decay_lengths_core_required

    dx_beam_duct = 0.58
    big_q_plasma = 0.0
    f_c_plasma_bootstrap = 0.0
    f_c_plasma_bootstrap_max = 0.9
    f_c_plasma_bootstrap_iter89 = 0.0
    f_c_plasma_bootstrap_nevins = 0.0
    f_c_plasma_bootstrap_sauter = 0.0
    f_c_plasma_bootstrap_wilson = 0.0
    f_c_plasma_bootstrap_sakai = 0.0
    f_c_plasma_bootstrap_aries = 0.0
    f_c_plasma_bootstrap_andrade = 0.0
    f_c_plasma_bootstrap_hoang = 0.0
    f_c_plasma_bootstrap_wong = 0.0
    bscf_gi_i = 0.0
    bscf_gi_ii = 0.0
    f_c_plasma_bootstrap_sugiyama_l = 0.0
    f_c_plasma_bootstrap_sugiyama_h = 0.0
    cboot = 1.0
    c_beam_total = 0.0
    f_c_plasma_diamagnetic_hender = 0.0
    f_c_plasma_diamagnetic_scene = 0.0
    f_c_plasma_diamagnetic = 0.0
    p_hcd_ecrh_injected_total_mw = 0.0
    eta_cd_hcd_primary = 0.0
    n_ecrh_harmonic = 2.0
    i_ecrh_wave_mode = 0
    e_beam_kev = 1.0e3
    eta_hcd_primary_injector_wall_plug = 0.3
    eta_hcd_secondary_injector_wall_plug = 0.3
    eta_ecrh_injector_wall_plug = 0.3
    eta_lowhyb_injector_wall_plug = 0.3
    eta_beam_injector_wall_plug = 0.3
    f_p_beam_injected_ions = 0.5
    p_beam_injected_mw = 0.0
    f_c_plasma_pfirsch_schluter_scene = 0.0
    p_beam_shine_through_mw = 0.0
    feffcd = 1.0
    f_p_beam_orbit_loss = 0.0
    f_radius_beam_tangency_rmajor = 1.05
    f_beam_tritium = 1e-6
    eta_cd_norm_hcd_primary = 0.0
    eta_cd_dimensionless_hcd_primary = 0.0
    eta_cd_norm_ecrh = 0.35
    xi_ebw = 0.8
    i_hcd_primary = 5
    i_hcd_secondary = 0
    i_hcd_calculations = 1
    f_p_beam_shine_through = 0.0
    dx_beam_shield = 0.5
    p_hcd_primary_extra_heat_mw = 0.0
    p_hcd_secondary_extra_heat_mw = 0.0
    p_hcd_injected_max = 150.0
    p_hcd_injected_electrons_mw = 0.0
    p_hcd_injected_ions_mw = 0.0
    p_hcd_injected_total_mw = 0.0
    p_hcd_injected_current_total_mw = 0.0
    p_hcd_secondary_injected_mw = 0.0
    f_c_plasma_internal = 0.0
    p_beam_orbit_loss_mw = 0.0
    f_c_plasma_pfirsch_schluter = 0.0
    pwpnb = 0.0
    radius_beam_tangency = 0.0
    radius_beam_tangency_max = 0.0
    n_beam_decay_lengths_core = 0.0
    n_beam_decay_lengths_core_required = 3.0
    eta_cd_norm_hcd_secondary = 0.0
    eta_cd_dimensionless_hcd_secondary = 0.0
    eta_cd_hcd_secondary = 0.0
    p_ebw_injected_mw = 0.0
    p_hcd_ecrh_electric_mw = 0.0
    p_hcd_ebw_electric_mw = 0.0
    c_hcd_primary_driven = 0.0
    c_hcd_secondary_driven = 0.0
    f_c_plasma_hcd_primary = 0.0
    f_c_plasma_hcd_secondary = 0.0
    eta_icrh_injector_wall_plug = 0.3
    eta_ebw_injector_wall_plug = 0.3
    p_hcd_primary_injected_mw = 0.0
    p_hcd_icrh_injected_total_mw = 0.0
    p_hcd_ebw_injected_total_mw = 0.0
    p_hcd_beam_injected_total_mw = 0.0
    p_hcd_icrh_injected_total_mw = 0.0
    p_hcd_lowhyb_electric_mw = 0.0
    p_hcd_lowhyb_injected_total_mw = 0.0