search

PlasmaControl / DESC

Stellarator Equilibrium and Optimization Suite
MIT License
95 stars 26 forks source link

Rotating ellipse test equilibrium not rotating #1217

Closed jmathewson1869 closed 2 months ago

jmathewson1869 commented 2 months ago

I am able to solve for an equilibrium from an mgrid coil file. However, when plotting the surfaces, they do not seem to vary as phi changes. Any help would be much appreciated. MGRID_COILS.zip

f0uriest commented 2 months ago

It looks like you're just calling eq.solve on the initial guess, which will solve a fixed boundary problem with the initial surface shape you give. After that you'll need to use the {Vacuum}BoundaryError objectives and eq.optimize to actually change the boundary shape to match the external field, see https://desc-docs.readthedocs.io/en/stable/notebooks/tutorials/free_boundary_equilibrium.html

f0uriest commented 2 months ago

Also note that the initial guess for a free boundary run generally needs to have roughly the correct major and minor radius, and fit inside any external coils, so you may need to modify the initial guess depending on the shape of your external field.

jmathewson1869 commented 2 months ago

Ok thank you for the help, I’ll check and see, I think the input has the bounds that I set. But I’ll make sure.

Best, James

On Thu, Aug 22, 2024 at 2:24 PM Rory Conlin @.***> wrote:

Also note that the initial guess for a free boundary run generally needs to have roughly the correct major and minor radius, and fit inside any external coils, so you may need to modify the initial guess depending on the shape of your external field.

— Reply to this email directly, view it on GitHub https://github.com/PlasmaControl/DESC/issues/1217#issuecomment-2305567696, or unsubscribe https://github.com/notifications/unsubscribe-auth/A27HG7RNGSCQX6HFQKFUDVLZSZCJNAVCNFSM6AAAAABM46QEXWVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDGMBVGU3DONRZGY . You are receiving this because you authored the thread.Message ID: @.***>

dpanici commented 2 months ago

Closing as assumed resolved, feel free to reopen if there are still issues @jmathewson1869