Skip to content
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.

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

🚧 3273 calculated bootstrap current fraction is overridden #3275

Open
wants to merge 8 commits into
base: main
Choose a base branch
from

Conversation

chris-ashe
Copy link
Collaborator

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

Description

A new constraint equation icc=92 and f-value fboot_max have been introduced to take over from:

 if current_drive_variables.bootipf > current_drive_variables.bscfmax:
                current_drive_variables.bootipf = min(
                    current_drive_variables.bootipf, current_drive_variables.bscfmax
                )

The above would override the calculated bootstrap current value from the model if its value ever went greater than bscfmax. Giving possible non-consistent solutions.

This may cause solutions now harder to find due to the initial conditions on the input file which can easily give bootstrap fractions greater than one and causing the injected power constraint icc=30 to crash also. Widening bounds is encouraged and the bootstrap fraction can still be directly set by setting bscfmax to the negative of the required value

The old err242 error catching variable which would output that the bootstrap current fraction limit has been enforced has now been removed from all files.

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.

@chris-ashe chris-ashe added Bug Something isn't working Physics Relating to the physics models labels Aug 9, 2024
@chris-ashe chris-ashe self-assigned this Aug 9, 2024
@chris-ashe chris-ashe linked an issue Aug 9, 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.

In addition to specific comments, you should also check the functionality with cboot.

documentation/proc-pages/physics-models/plasma_current.md Outdated Show resolved Hide resolved
documentation/proc-pages/physics-models/plasma_current.md Outdated Show resolved Hide resolved
@chris-ashe
Copy link
Collaborator Author

cboot is applied to the front of each calculation before assignment. E.g:

current_drive_variables.bscf_sauter = (
            current_drive_variables.cboot * self.bootstrap_fraction_sauter()
        )

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 the changes.

@chris-ashe chris-ashe changed the title 3273 calculated bootstrap current fraction is overridden 🚧 3273 calculated bootstrap current fraction is overridden Aug 21, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Bug Something isn't working Physics Relating to the physics models
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Calculated bootstrap current fraction is overridden
2 participants