Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Original file line number Diff line number Diff line change
Expand Up @@ -114,6 +114,6 @@ The recirculating power comprises the electrical power required to run all of th
electrical systems surrounding the fusion power core, plus the on-site building services, offices,
etc., as shown in Figure 1. Of these, the cryogenic plant power includes the power required to cool
the TF coils from the neutron power absorbed by the coils, the PF coils (as defined by the ratio of
the total PF coil stored energy to the fusion power pulse time `t_pulse_repetition`), and other 'cold' components.
the total PF coil stored energy to the fusion power pulse time `t_plant_pulse_plasma_present`), and other 'cold' components.

[^1]: C. Harrington, *"Development and Implementation of Improved Balance of Plant Models for PROCESS"*, CCFE C5.M15 Milestone Report, August 2014 (copy stored as CCFE internal note T&M/PKNIGHT/PROCESS/027)
2 changes: 1 addition & 1 deletion documentation/io/utilities.md
Original file line number Diff line number Diff line change
Expand Up @@ -109,7 +109,7 @@ Produces a three-page PDF file in the same directory as the input MFILE. The PDF
`Plasma Composition` - Number densities of several ion species relative to the electron density.

`Coil Currents etc` - Peak coil currents of the PF coils in $\text{MA}$, flux swing of the central solenoid
used for startup and total available in $\text{Wb}$. Total burn time `t_burn` in hrs.
used for startup and total available in $\text{Wb}$. Total burn time `t_plant_pulse_burn` in hrs.

`Cost of electricity` - This is the cost of electricity in $ $/ \text{MWh}$. Check the respective cost model
for the reference year of the inflation used.
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -143,7 +143,7 @@ where $\rho_\text{p}$ is the calculated [plasma resistivity](./plasma_resistive_
The total flux required is then simply found by multiplying the loop voltage above by the required duration of the burn phase:

$$
\overbrace{\Phi_{\text{res,burn}}}^{\texttt{vs_burn_required}} = \overbrace{I_{\text{p}} \rho_\text{p} f_{\text{ind}}}^{\texttt{v_burn_resistive}} \times \overbrace{T_{\text{burn}}}^{\texttt{t_burn}}
\overbrace{\Phi_{\text{res,burn}}}^{\texttt{vs_burn_required}} = \overbrace{I_{\text{p}} \rho_\text{p} f_{\text{ind}}}^{\texttt{v_burn_resistive}} \times \overbrace{T_{\text{burn}}}^{\texttt{t_plant_pulse_burn}}
$$

----------------
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -679,9 +679,9 @@ This constraint can be activated by stating `icc = 41` in the input file.

The value of `tohsm` can be set to the required minimum plasma current ramp up time at the start of a pulse. The scaling value `ft_current_ramp_up` can be varied also

The calculated plasma current ramp up time `t_current_ramp_up` is dictated by the [pulsed plant operation configuration](../pulsed-plant.md).
The calculated plasma current ramp up time `t_plant_pulse_plasma_current_ramp_up` is dictated by the [pulsed plant operation configuration](../pulsed-plant.md).

This constraint will ensure that the value of `t_current_ramp_up` is always greater than or equal to `tohsm`
This constraint will ensure that the value of `t_plant_pulse_plasma_current_ramp_up` is always greater than or equal to `tohsm`

--------------------

Expand Down
24 changes: 12 additions & 12 deletions documentation/physics-models/pulsed-plant.md
Original file line number Diff line number Diff line change
Expand Up @@ -29,16 +29,16 @@ change equal to the maximum proposed in [^1], or it can be set by the user. The
constraint is likely to depend on whether the ramp-up is purely inductive or includes current drive,
but this is not taken ito account.

In the steady-state scenario (`i_pulsed_plant` = 0), the plasma current ramp-up time `t_current_ramp_up` is determined as follows.
In the steady-state scenario (`i_pulsed_plant` = 0), the plasma current ramp-up time `t_plant_pulse_plasma_current_ramp_up` is determined as follows.

- If `i_t_current_ramp_up` = 0, the rate of change of plasma current is 0.5 MA/s. The PF coil ramp time `t_precharge`
and shutdown time `t_ramp_down` are (arbitrarily) set equal to `t_current_ramp_up`.
- If `i_t_current_ramp_up` $\neq$ 0, the plasma current ramp-up time `t_current_ramp_up`, and the PF coil ramp
- If `i_t_current_ramp_up` = 0, the rate of change of plasma current is 0.5 MA/s. The PF coil ramp time `t_plant_pulse_coil_precharge`
and shutdown time `t_plant_pulse_plasma_current_ramp_down` are (arbitrarily) set equal to `t_plant_pulse_plasma_current_ramp_up`.
- If `i_t_current_ramp_up` $\neq$ 0, the plasma current ramp-up time `t_plant_pulse_plasma_current_ramp_up`, and the PF coil ramp
and shutdown times are input parameters.

In the pulsed scenario, (`i_pulsed_plant` = 1), the plasma current ramp-up time `t_current_ramp_up` is an input, and it
In the pulsed scenario, (`i_pulsed_plant` = 1), the plasma current ramp-up time `t_plant_pulse_plasma_current_ramp_up` is an input, and it
can be set as an iteration variable (65). The ramp-up and shutdown time in the pulsed case are set
equal to `t_current_ramp_up`. To ensure that the plasma current ramp rate during start-up is prevented from being
equal to `t_plant_pulse_plasma_current_ramp_up`. To ensure that the plasma current ramp rate during start-up is prevented from being
too high, as governed by the requirement to maintain plasma stability by ensuring that the induced
current has time to diffuse into the body of the plasma, constraint equation no. 41 should be
turned on with iteration variable no. 66 `ft_current_ramp_up` and input `t_current_ramp_up_min`, the minimum plasma current
Expand All @@ -55,12 +55,12 @@ constraint equation no. 13 and iteration variable no 21 (`ft_burn_min`).

Over the course of a pulse, the timings are detailed as:

- Precharge (`t_precharge`) - the CS current ramps from zero to maximum value. The other PF coils also ramp from zero to their required values.
- Current ramp-up (`t_current_ramp_up`) - The plasma current ramps up to approx full value. Auxiliary heating is possibly on.
- Fusion ramp (`t_fusion_ramp`) - The plasma temperature and density rise to the full values. The CS and other PF coil currents all change steadily. Auxiliary heating is on.
- Burn time (`t_burn`) - Flat-top duration. The plasma is approximately steady. Fusion power and electricity are produced. The CS and other PF coil currents all change steadily in a pulsed reactor, but are constant for a "steady-state" reactor. Auxiliary heating is on.
- Ramp-down (`t_ramp_down`) - The plasma current, density and temperature all ramp down to zero simultaneously. As a starter for ten we could assume that the CS and PF coil currents also ramp to zero at the same time. Auxiliary heating is possibly on.
- Between pulse (`t_between_pulse`) - CS and PF coil currents are zero - a few minutes may be required to permit vacuum pumping. May be much longer for an experimental device such as DEMO.
- Precharge (`t_plant_pulse_coil_precharge`) - the CS current ramps from zero to maximum value. The other PF coils also ramp from zero to their required values.
- Current ramp-up (`t_plant_pulse_plasma_current_ramp_up`) - The plasma current ramps up to approx full value. Auxiliary heating is possibly on.
- Fusion ramp (`t_plant_pulse_fusion_ramp`) - The plasma temperature and density rise to the full values. The CS and other PF coil currents all change steadily. Auxiliary heating is on.
- Burn time (`t_plant_pulse_burn`) - Flat-top duration. The plasma is approximately steady. Fusion power and electricity are produced. The CS and other PF coil currents all change steadily in a pulsed reactor, but are constant for a "steady-state" reactor. Auxiliary heating is on.
- Ramp-down (`t_plant_pulse_plasma_current_ramp_down`) - The plasma current, density and temperature all ramp down to zero simultaneously. As a starter for ten we could assume that the CS and PF coil currents also ramp to zero at the same time. Auxiliary heating is possibly on.
- Between pulse (`t_plant_pulse_dwell`) - CS and PF coil currents are zero - a few minutes may be required to permit vacuum pumping. May be much longer for an experimental device such as DEMO.
- Pulse repitition (`t_pulse_repitition`) - Sum of all the above times.

A plot showing schematically these timings over a pulse can be found in Figure 2.
Expand Down
20 changes: 10 additions & 10 deletions examples/data/csv_output_large_tokamak_MFILE.DAT
Original file line number Diff line number Diff line change
Expand Up @@ -557,13 +557,13 @@
Required_volt-seconds_during_start-up_(Wb)______________________________ (vssoft)______________________ 2.8401E+02
Available_volt-seconds_during_burn_(Wb)_________________________________ (vsmax)_______________________ 2.8571E+02
# Times #
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_current_ramp_up)________________________ 1.6699E+02
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_plant_pulse_coil_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_plant_pulse_plasma_current_ramp_up)________________________ 1.6699E+02
Heating_time_(s)________________________________________________________ (theat)_______________________ 1.0000E+01
Burn_time_(s)___________________________________________________________ (t_burn)_______________________ 7.1999E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_ramp_down)_______________________ 1.6699E+02
Time_between_pulses_(s)_________________________________________________ (t_between_pulse)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_cycle)______________________ 9.8441E+03 OP
Burn_time_(s)___________________________________________________________ (t_plant_pulse_burn)_______________________ 7.1999E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_plant_pulse_plasma_current_ramp_down)_______________________ 1.6699E+02
Time_between_pulses_(s)_________________________________________________ (t_plant_pulse_dwell)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_plant_pulse_total)______________________ 9.8441E+03 OP
kallenbach_switch_______________________________________________________ (kallenbach_switch)___________ 0
# Radial Build #
TF_coil_radial_placement_switch_________________________________________ (i_tf_inside_cs)____________________ 0
Expand Down Expand Up @@ -1048,8 +1048,8 @@
Chamber_pressure_after_burn_(Pa)________________________________________ (pend)________________________ 1.6766E-01 OP
Chamber_pressure_before_burn_(Pa)_______________________________________ (pstart)______________________ 1.6766E-03
Allowable_pumping_time_switch___________________________________________ (i_vac_pump_dwell)__________________ 0
Dwell_time_between_burns_(s)____________________________________________ (t_between_pulse.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_precharge.)______________________ 5.0000E+02
Dwell_time_between_burns_(s)____________________________________________ (t_plant_pulse_dwell.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_plant_pulse_coil_precharge.)______________________ 5.0000E+02
Allowable_pumping_time_between_burns_(s)________________________________ (tpump)_______________________ 1.8000E+03
Required_D-T_pump_speed_(m3/s)__________________________________________ (s(2))________________________ 5.7790E+00 OP
D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(2))_____________________ 5.6929E+01 OP
Expand Down Expand Up @@ -1631,13 +1631,13 @@ plasma_res_factor = 0.7
i_pulsed_plant = 1

* dwell time [s]
t_between_pulse = 1800.0
t_plant_pulse_dwell = 1800.0

* Switch for pulse timing calculations
pulsetimings = 0

* CS ramp up time [s]
t_precharge = 500.0
t_plant_pulse_coil_precharge = 500.0

* Current drive *
*---------------*
Expand Down
20 changes: 10 additions & 10 deletions examples/data/large_tokamak_1_MFILE.DAT
Original file line number Diff line number Diff line change
Expand Up @@ -555,13 +555,13 @@
Required_volt-seconds_during_start-up_(Wb)______________________________ (vssoft)______________________ 2.8601E+02
Available_volt-seconds_during_burn_(Wb)_________________________________ (vsmax)_______________________ 2.7014E+02
# Times #
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_current_ramp_up)________________________ 1.6631E+02
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_plant_pulse_coil_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_plant_pulse_plasma_current_ramp_up)________________________ 1.6631E+02
Heating_time_(s)________________________________________________________ (theat)_______________________ 1.0000E+01
Burn_time_(s)___________________________________________________________ (t_burn)_______________________ 7.2035E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_ramp_down)_______________________ 1.6631E+02
Time_between_pulses_(s)_________________________________________________ (t_between_pulse)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_cycle)______________________ 9.8463E+03 OP
Burn_time_(s)___________________________________________________________ (t_plant_pulse_burn)_______________________ 7.2035E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_plant_pulse_plasma_current_ramp_down)_______________________ 1.6631E+02
Time_between_pulses_(s)_________________________________________________ (t_plant_pulse_dwell)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_plant_pulse_total)______________________ 9.8463E+03 OP
# Radial Build #
TF_coil_radial_placement_switch_________________________________________ (i_tf_inside_cs)____________________ 0
TF_coil_radial_placement_switch_________________________________________ (i_tf_inside_cs)____________________ 0
Expand Down Expand Up @@ -1043,8 +1043,8 @@
Chamber_pressure_after_burn_(Pa)________________________________________ (pend)________________________ 1.6698E-01 OP
Chamber_pressure_before_burn_(Pa)_______________________________________ (pstart)______________________ 1.6698E-03
Allowable_pumping_time_switch___________________________________________ (i_vac_pump_dwell)__________________ 0
Dwell_time_between_burns_(s)____________________________________________ (t_between_pulse.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_precharge.)______________________ 5.0000E+02
Dwell_time_between_burns_(s)____________________________________________ (t_plant_pulse_dwell.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_plant_pulse_coil_precharge.)______________________ 5.0000E+02
Allowable_pumping_time_between_burns_(s)________________________________ (tpump)_______________________ 1.8000E+03
Required_D-T_pump_speed_(m3/s)__________________________________________ (s(2))________________________ 5.7790E+00 OP
D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(2))_____________________ 5.6700E+01 OP
Expand Down Expand Up @@ -1625,13 +1625,13 @@ plasma_res_factor = 0.7
i_pulsed_plant = 1

* dwell time [s]
t_between_pulse = 1800.0
t_plant_pulse_dwell = 1800.0

* Switch for pulse timing calculations
pulsetimings = 0

* CS ramp up time [s]
t_precharge = 500.0
t_plant_pulse_coil_precharge = 500.0

* Current drive *
*---------------*
Expand Down
20 changes: 10 additions & 10 deletions examples/data/large_tokamak_2_MFILE.DAT
Original file line number Diff line number Diff line change
Expand Up @@ -555,13 +555,13 @@
Required_volt-seconds_during_start-up_(Wb)______________________________ (vssoft)______________________ 2.8601E+02
Available_volt-seconds_during_burn_(Wb)_________________________________ (vsmax)_______________________ 2.7014E+02
# Times #
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_current_ramp_up)________________________ 1.6631E+02
Initial_charge_time_for_CS_from_zero_current_(s)________________________ (t_plant_pulse_coil_precharge)_______________________ 5.0000E+02
Plasma_current_ramp-up_time_(s)_________________________________________ (t_plant_pulse_plasma_current_ramp_up)________________________ 1.6631E+02
Heating_time_(s)________________________________________________________ (theat)_______________________ 1.0000E+01
Burn_time_(s)___________________________________________________________ (t_burn)_______________________ 7.2035E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_ramp_down)_______________________ 1.6631E+02
Time_between_pulses_(s)_________________________________________________ (t_between_pulse)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_cycle)______________________ 9.8463E+03 OP
Burn_time_(s)___________________________________________________________ (t_plant_pulse_burn)_______________________ 7.2035E+03 OP
Reset_time_to_zero_current_for_CS_(s)___________________________________ (t_plant_pulse_plasma_current_ramp_down)_______________________ 1.6631E+02
Time_between_pulses_(s)_________________________________________________ (t_plant_pulse_dwell)______________________ 1.8000E+03
Total_plant_cycle_time_(s)______________________________________________ (t_plant_pulse_total)______________________ 9.8463E+03 OP
# Radial Build #
TF_coil_radial_placement_switch_________________________________________ (i_tf_inside_cs)____________________ 0
TF_coil_radial_placement_switch_________________________________________ (i_tf_inside_cs)____________________ 0
Expand Down Expand Up @@ -1043,8 +1043,8 @@
Chamber_pressure_after_burn_(Pa)________________________________________ (pend)________________________ 1.6698E-01 OP
Chamber_pressure_before_burn_(Pa)_______________________________________ (pstart)______________________ 1.6698E-03
Allowable_pumping_time_switch___________________________________________ (i_vac_pump_dwell)__________________ 0
Dwell_time_between_burns_(s)____________________________________________ (t_between_pulse.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_precharge.)______________________ 5.0000E+02
Dwell_time_between_burns_(s)____________________________________________ (t_plant_pulse_dwell.)_____________________ 1.8000E+03
CS_ramp-up_time_burns_(s)_______________________________________________ (t_plant_pulse_coil_precharge.)______________________ 5.0000E+02
Allowable_pumping_time_between_burns_(s)________________________________ (tpump)_______________________ 1.8000E+03
Required_D-T_pump_speed_(m3/s)__________________________________________ (s(2))________________________ 5.7790E+00 OP
D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(2))_____________________ 5.6700E+01 OP
Expand Down Expand Up @@ -1625,13 +1625,13 @@ plasma_res_factor = 0.7
i_pulsed_plant = 1

* dwell time [s]
t_between_pulse = 1800.0
t_plant_pulse_dwell = 1800.0

* Switch for pulse timing calculations
pulsetimings = 0

* CS ramp up time [s]
t_precharge = 500.0
t_plant_pulse_coil_precharge = 500.0

* Current drive *
*---------------*
Expand Down
Loading
Loading