Skip to content
DraconiticResearch
Search Sources API
186,625 studies31 sourcesUpdated 58 min ago
PreprintarXiv licencearXiv

Dynamic Optimization Workflow for Gyrokinetic Stability and Transport

Abstract

A dynamic optimization workflow for linear gyrokinetic stability that couples the gyrokinetic code CGYRO with the libEnsemble framework for ensemble management has been developed. The workflow is designed to efficiently explore multi-dimensional plasma parameter spaces and to address the practical challenges of simulation-driven optimization in fusion plasma physics, including expensive objective evaluations, unavailable gradient information, non-smoothness arising from transitions between dominant instability branches, and simultaneous multi-core simulations on high-performance computing resources. Rather than optimizing the growth rate at a single wavenumber, the objective maximizes the growth rate-to-wavenumber ratio over the ion-scale spectrum. A derivative-free, model-based optimization algorithm employs manifold sampling to handle the nonsmooth composite structure of the objective function. The approach is applied to the optimization of local magnetic flux-surface shaping parameters, with an equilibrium baseline drawn from DIII-D negative triangularity experimental parameters. The workflow is general and extensible and provides a foundation for future optimization studies involving nonlinear transport objectives and higher-dimensional parameter spaces.

Related

  1. AI-Accelerated Gyrokinetic Predictions of Turbulent Transport for Stellarator Design Optimization and Experimental Planning
    R. Michael Churchill, Matt Landreman, Jong Youl Choi +10
    arXivPreprintarXiv licenceFull textPlasma physicsarXiv
  2. PQLS: A Quasilinear Gyrokinetic Transport Solver with a Bayesian Saturation-Rule Closure
    F. Wilms, A. Agrawal, J. J. Freigang +3
    arXivPreprintCC BY 4.0Full textPlasma physicsarXiv
  3. Kinetic modifications to parallel transport in a low-recycling scrape-off layer: A one dimensional gyrokinetic study
    A. Maan, E. Perez, M. Francisquez +5
    arXivPreprintCC BY 4.0Full textPlasma physicsarXiv
  4. ELMO: An Uncertainty-Aware Simulation-to-Surrogate Workflow for Fast Pedestal Linear-Stability Prediction
    Nami Li, X. Q. Xu, T. Osborne +5
    arXivPreprintarXiv licenceFull textPlasma physicsarXiv
  5. A dynamic pore-network model with discrete microbubble transport
    Kristoffer Skjelanger, Anna Mareike Kostelecky, Maziar Veyskarami +1
    arXivPreprintCC BY-NC-ND 4.0Full textFluid dynamicsarXiv
BibTeX
@misc{belli2026dynamic,
  title = {{Dynamic Optimization Workflow for Gyrokinetic Stability and Transport}},
  author = {E. A. Belli and J. Larson and J. Candy and I. Sfiligoi},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.17496},
  eprint = {2610.17496},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.17496},
}
RIS
TY  - UNPB
TI  - Dynamic Optimization Workflow for Gyrokinetic Stability and Transport
AU  - E. A. Belli
AU  - J. Larson
AU  - J. Candy
AU  - I. Sfiligoi
PY  - 2026
DA  - 2026/10/12
PB  - arXiv
DO  - 10.48550/arxiv.2610.17496
UR  - https://arxiv.org/abs/2610.17496
AB  - A dynamic optimization workflow for linear gyrokinetic stability that couples the gyrokinetic code CGYRO with the libEnsemble framework for ensemble management has been developed. The workflow is designed to efficiently explore multi-dimensional plasma parameter spaces and to address the practical challenges of simulation-driven optimization in fusion plasma physics, including expensive objective evaluations, unavailable gradient information, non-smoothness arising from transitions between dominant instability branches, and simultaneous multi-core simulations on high-performance computing resources. Rather than optimizing the growth rate at a single wavenumber, the objective maximizes the growth rate-to-wavenumber ratio over the ion-scale spectrum. A derivative-free, model-based optimization algorithm employs manifold sampling to handle the nonsmooth composite structure of the objective function. The approach is applied to the optimization of local magnetic flux-surface shaping parameters, with an equilibrium baseline drawn from DIII-D negative triangularity experimental parameters. The workflow is general and extensible and provides a foundation for future optimization studies involving nonlinear transport objectives and higher-dimensional parameter spaces.
ER  -