my_system_code/Input.py at main · chencx04/my_system_code · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
# 测试用程序
def assignvalue(input_data,time_data,power_data,hcd_data,prof_data):
    # types import SimpleNamespace

    input_data.i_plasma_geometry=0
    input_data.i_plasma_pedestal=1
    input_data.i_plasma_current = 4
    input_data.aspect=1.7
    input_data.rmajor=3.1
    input_data.kappa=2.5
    input_data.delta=0.6
    input_data.q95=5
    input_data.b_toroidal_rmajor=4
    # input_data.ndensity_e_vol_avg=8.812611330723181e+19
    input_data.ndensity_e_on_axis=1.2e20  # 轴电子密度
    input_data.ndensity_e_ped=3e19
    input_data.ndensity_e_sep=1e19
    # input_data.temp_e_vol_avg=8.386338252709802
    input_data.temp_e_on_axis=12.0  # 轴电子温度
    input_data.temp_e_ped=5
    input_data.temp_e_sep=0.1 # process在输入里的是电子温度，离子温度通过比值计算
    input_data.alpha_temp=2
    input_data.beta_temp=1
    input_data.alpha_ndensity=1
    input_data.rho_ped_ndensity=0.94
    input_data.rho_ped_temp=0.94
    input_data.imp_name=[] # array ["","",""]
    input_data.f_density_alpha_electron=0.1
    input_data.f_density_h_electron=0
    input_data.f_density_beam_electron=0
    input_data.f_density_imp_electron=[] # array
    # 燃料分配比
    input_data.f_density_he3_fuel=0
    input_data.f_density_d_fuel=0.5
    input_data.f_density_t_fuel=0.5
    input_data.f_density_t_beam = 0.5
    input_data.i_plasma_ignited = 1
    input_data.f_temp_electron_ion=1
    input_data.i_plasma_current=4
    input_data.i_q95_as_input_data=1
    input_data.nrho=51
    input_data.i_diamagnetic_current=1
    input_data.i_ps_current=1
    input_data.i_bootsrap_current=1
    input_data.cboot=1 # bootstrap current fraction multiplier
    input_data.q0=1
    input_data.f_current_non_inductive=1
    input_data.i_hcd_secondary=0
    input_data.i_hcd_primary=5  # 不支持 9 和 11 ?
    input_data.i_hcd_calculations=1
    input_data.i_plasma_pedestal=1


    input_data.f_r_beam_tangency_rmajor =1.05
    input_data.e_beam_kev                    =1000
    input_data.feffcd                        =1
    input_data.eta_cd_norm_ecrh              =0.35
    input_data.n_ecrh_harmonic               =2
    input_data.xi_ebw                        =0.8
    input_data.i_ecrh_wave_mode              =0

    input_data.eta_lowhyb_injector_wall_plug =0.95
    input_data.eta_icrh_injector_wall_plug=0.95
    input_data.eta_ecrh_injector_wall_plug=0.95
    input_data.eta_ebw_injector_wall_plug=0.95
    input_data.eta_beam_injector_wall_plug=0.95
    input_data.power_hcd_secondary_injected_mw=0
    input_data.power_hcd_primary_extra_heat_mw=0
    input_data.power_hcd_secondary_extra_heat_mw=0
    input_data.f_power_beam_orbit_loss=0

    input_data.i_alpha_j=0
    input_data.i_ind_plasma_internal_norm=0
    input_data.alpha_j=1
    input_data.ind_plasma_internal_norm=0.9
    input_data.f_current_bootstrap_max=1
    input_data.rho_plasma_core_norm=0.6
    input_data.f_power_plasma_core_rad_reduction=1
    input_data.f_syn_reflect=0.6
    input_data.f_resis_plasma=1
    input_data.beamfus0=1
    input_data.betbm0=1

    input_data.f_alpha_power_deposited=1
    input_data.i_rad_loss=2

    input_data.i_dimension = 0
    input_data.i_plasma_q95 = 0
    input_data.i_pulsed_plant = 0
    input_data.i_t_current_ramp_up = 0
    input_data.i_l_h_threshold = 1
    input_data.f_current_bootstrap = 0.0

    hcd_data.cboot = input_data.cboot
    hcd_data.f_current_bootstrap_max = input_data.f_current_bootstrap_max

    time_data.t_energy_confinement = 0.5
    time_data.pulsetimings = 0.0
    time_data.t_plant_pulse_coil_precharge = 0.0
    time_data.t_plant_pulse_plasma_current_ramp_up = 0.0
    time_data.t_plant_pulse_plasma_current_ramp_down = 0.0
    time_data.t_plant_pulse_fusion_ramp = 0.0
    time_data.t_plant_pulse_burn = 0.5
    time_data.t_plant_pulse_dwell = 0.0