Skip to content

Conversation

@chris-ashe
Copy link
Collaborator

@chris-ashe chris-ashe commented Aug 6, 2024

Description

A new model by Sakai et al has been implemented
image

This should be applicable for all aspect ratios.

Unit tests created for large tokamak and ST case

Sakai et al. says at the start of the paper that the dataset used to fit the model includes the toroidal diamagnetic current. Thus we remove the diamagnetic fraction from the calculated value when the function is called

Checklist

I confirm that I have completed the following checks:

  • I have justified any large differences in the regression tests caused by this pull request in the comments.
  • I have added new tests where appropriate for the changes I have made.
  • If I have had to change any existing unit or integration tests, I have justified this change in the pull request comments.
  • If I have made documentation changes, I have checked they render correctly.
  • I have added documentation for my change, if appropriate.

This commit adds a new bootstrap current fraction model based on the Sakai formula. The `bootstrap_fraction_sakai` method is implemented in the `Physics` class in the `physics.py` file. This model calculates the bootstrap fraction using a set of parameters such as plasma beta, bootstrap current coefficient, safety factor, neoclassical viscosity parameters, and epsilon. The formula is derived from the work of Sakai et al. [^1^].

[^1^]: Sakai, R., Fujita, T., Okamoto, A. (2019). Derivation of bootstrap current fraction scaling formula for 0-D system code analysis. Fusion Engineering and Design, 149, 111322. https://doi.org/10.1016/j.fusengdes.2019.111322
@chris-ashe chris-ashe added the Physics Relating to the physics models label Aug 6, 2024
@chris-ashe chris-ashe self-assigned this Aug 6, 2024
@chris-ashe chris-ashe linked an issue Aug 6, 2024 that may be closed by this pull request
Copy link
Collaborator

@stuartmuldrew stuartmuldrew left a comment

Choose a reason for hiding this comment

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

Minor comment about diamagnetic current options

This commit modifies the `Physics` class in the `physics.py` file to update the calculation of the Sakai bootstrap current. The Sakai formula states that the ACCOME dataset used has the toroidal diamagnetic current included in the bootstrap current. Therefore, the commit removes the PROCESS calculated diamagnetic current from the bootstrap current calculation. This change improves the accuracy of the bootstrap current fraction model.
@chris-ashe chris-ashe force-pushed the 954-bootstrap-current-fraction-scaling-formula-for-steady-state-machines branch from ce787b5 to f1414f4 Compare August 7, 2024 09:14
@mkovari
Copy link
Collaborator

mkovari commented Aug 8, 2024

Thanks - a helpful piece of work.
Does it change the results significantly?

Copy link
Collaborator

@stuartmuldrew stuartmuldrew left a comment

Choose a reason for hiding this comment

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

Happy with error message addition.

@chris-ashe chris-ashe merged commit 53f99cd into main Aug 13, 2024
@chris-ashe chris-ashe deleted the 954-bootstrap-current-fraction-scaling-formula-for-steady-state-machines branch August 13, 2024 08:14
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Physics Relating to the physics models

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Bootstrap current fraction scaling formula for steady-state machines

4 participants