Suppressing Instability in a Vlasov-Poisson System by an External Electric Field Through Constrained Optimization

05/29/2023
by   Lukas Einkemmer, et al.
0

Fusion energy offers the potential for the generation of clean, safe, and nearly inexhaustible energy. While notable progress has been made in recent years, significant challenges persist in achieving net energy gain. Improving plasma confinement and stability stands as a crucial task in this regard and requires optimization and control of the plasma system. In this work, we deploy a PDE-constrained optimization formulation that uses a kinetic description for plasma dynamics as the constraint. This is to optimize, over all possible controllable external electric fields, the stability of the plasma dynamics under the condition that the Vlasov–Poisson (VP) equation is satisfied. For computing the functional derivative with respect to the external field in the optimization updates, the adjoint equation is derived. Furthermore, in the discrete setting, where we employ the semi-Lagrangian method as the forward solver, we also explicitly formulate the corresponding adjoint solver and the gradient as the discrete analogy to the adjoint equation and the Frechet derivative. A distinct feature we observed of this constrained optimization is the complex landscape of the objective function and the existence of numerous local minima, largely due to the hyperbolic nature of the VP system. To overcome this issue, we utilize a gradient-accelerated genetic algorithm, leveraging the advantages of the genetic algorithm's exploration feature to cover a broader search of the solution space and the fast local convergence aided by the gradient information. We show that our algorithm obtains good electric fields that are able to maintain a prescribed profile in a beam shaping problem and uses nonlinear effects to suppress plasma instability in a two-stream configuration.

READ FULL TEXT

page 14

page 15

page 16

page 17

page 18

page 19

page 20

page 21

research
08/23/2019

Mathematical and numerical study of a three-dimensional inverse eddy current problem

We study an inverse problem associated with an eddy current model. We fi...
research
05/06/2023

Differential-Equation Constrained Optimization With Stochasticity

Most inverse problems from physical sciences are formulated as PDE-const...
research
04/21/2021

A convergent numerical scheme for a model of liquid crystal dynamics subjected to an electric field

We present a convergent and constraint-preserving numerical discretizati...
research
09/25/2022

Monte Carlo Gradient in Optimization Constrained by Radiative Transport Equation

Can Monte Carlo (MC) solvers be directly used in gradient-based methods ...
research
03/13/2017

A Lagrangian Gauss-Newton-Krylov Solver for Mass- and Intensity-Preserving Diffeomorphic Image Registration

We present an efficient solver for diffeomorphic image registration prob...
research
07/13/2022

Optimal control of dielectric elastomer actuated multibody dynamical systems

In this work, a simulation model for the optimal control of dielectric e...

Please sign up or login with your details

Forgot password? Click here to reset