Skip to content

Fix PF electrical accounting for Ohmic heating #4664

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
grmtrkngtn wants to merge 1 commit into
base: main
Choose a base branch
from
Open

Fix PF electrical accounting for Ohmic heating #4664

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
14 changes: 6 additions & 8 deletions process/models/power.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -594,15 +594,12 @@ def pfpwr(self, output: bool, pulse_timings: PulseTimings):
* abs(self.data.pf_coil.c_pf_coil_turn_peak_input[idx_circuit])
/ self.data.pf_power.pfckts
)

# PF wall plug power dissipated in power supply for ohmic heating (MW)
# This is additional to that required for moving stored energy around
# p_pf_electric_supplies_mw = self.data.physics.p_plasma_ohmic_mw
# / self.data.pf_coil.etapsu
wall_plug_ohmicmw = self.data.physics.p_plasma_ohmic_mw * (
1.0e0 / self.data.pf_coil.etapsu - 1.0e0
# PF wall-plug power required for ohmic heating (MW)
wall_plug_ohmicmw = (
self.data.physics.p_plasma_ohmic_mw / self.data.pf_coil.etapsu

@ajpearcey ajpearcey Oct 8, 2026

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I can be very misinformed, but why does electrical losses in the PF coil depend linearly on the plasma resistivity?

@grmtrkngtn grmtrkngtn Oct 8, 2026

Copy link
Copy Markdown
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This term is intended to represent the power-supply loss associated with transferring the ohmic heating power to the plasma, rather than a resistive loss in the PF coils themselves.

Since the plasma ohmic heating is

P_ohmic = V_loop * I_p = I_p^2 * R_plasma ​

and therefore depends on the plasma resistance, a fixed power-supply efficiency gives a loss

P_ohmic * (1/η_psu − 1)

The PF coil/busbar losses are accounted for separately in pfpowermw.

@chris-ashe chris-ashe Oct 13, 2026

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looking here for some reason we only say there is resistive heating power if their is inductive current drive. What happens in the ST case when there is no inductive flattop?

)
# Total mean wall plug power dissipated in PFC and CS power supplies. Issue #713

# Total mean PF/CS electrical demand
self.data.pf_coil.p_pf_electric_supplies_mw = wall_plug_ohmicmw + pfpowermw

# Output Section
Expand Down Expand Up @@ -1678,6 +1675,7 @@ def plant_electric_production(self):
# by the cooling system (MW)
self.data.heat_transport.p_plant_secondary_heat_mw = (
self.data.power.p_plant_core_systems_elec_mw
- self.data.physics.p_plasma_ohmic_mw
+ self.data.heat_transport.p_hcd_electric_loss_mw
+ self.data.heat_transport.p_coolant_pump_loss_total_mw
+ self.data.heat_transport.p_div_secondary_heat_mw
Expand Down
Loading