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
- AI-Accelerated Gyrokinetic Predictions of Turbulent Transport for Stellarator Design Optimization and Experimental Planning
- PQLS: A Quasilinear Gyrokinetic Transport Solver with a Bayesian Saturation-Rule Closure
- Kinetic modifications to parallel transport in a low-recycling scrape-off layer: A one dimensional gyrokinetic study
- ELMO: An Uncertainty-Aware Simulation-to-Surrogate Workflow for Fast Pedestal Linear-Stability Prediction
- A dynamic pore-network model with discrete microbubble transport
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 -