Add Generomak Equilibrium #335
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
vsnever
reviewed
Nov 29, 2021
| "b_vacuum_radius": 1.48, | ||
| "b_vacuum_magnitude": 2, | ||
| "time": 0.0, | ||
| "f_profile": [ |
There was a problem hiding this comment.
This profile is inconsistent with the Bvac0 value. First of all, it's negative, while Bvac0 is positive. Also, at ψnorm = 1, it should give Bvac0 * R0 = 2.96, but gives -2.93 instead. Can you please renormalize it?
There was a problem hiding this comment.
Thanks for spotting this. I corrected the coordinate system and recalculated f and q profiles.
There was a problem hiding this comment.
Thanks, @Mateasek. Everything seems fine now.
vsnever
approved these changes
Dec 1, 2021
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR adds Generomak Equilibrium
The original equilibrium data were provided by @vsnever.
The equilibrium includes the first wall contoure, which was used to create the first wall meshes. The contour is copying closely the solps simulation provided by @vsnever (at least in the important regions)
The psi grid was clipped to extend only over the vessel.
Positions of the magnetic axis and x-point were found for a denser (interpolated) grid (to reach a better precision)
The position of the x-point determined the new psi at the last closed flux surface (either I misunderstood what the original data meant or the value was completely off)
Last closed flux surface polygon found as a contour given the psi at x-point
The strike-points were found as an intersection of divertor legs contour and the first wall polygon
This PR is connected to issues #211 and #210