Sparse identification of a predator-prey system from simulation data of a convection model

The use of low-dimensional dynamical systems as reduced models for plasma dynamics is useful as solving an initial value problem requires much less computational resources than fluid simulations. We utilize a data-driven modeling approach to identify a reduced model from simulation data of a convection problem. A convection model with a pressure source centered at the inner boundary models the edge dynamics of a magnetically confined plasma. The convection problem undergoes a sequence of bifurcations as the strength of the pressure source increases. The time evolution of the energies of the pressure profile, the turbulent flow, and the zonal flow capture the fundamental dynamic behavior of the full system. By applying the sparse identification of nonlinear dynamics (SINDy) method, we identify a predator-prey type dynamical system that approximates the underlying dynamics of the three energy state variables. A bifurcation analysis of the system reveals consistency between the bifurcation structures, observ...

[1]  S. Chapman,et al.  Transitions to improved confinement regimes induced by changes in heating in zero-dimensional models for tokamak plasmas , 2014 .

[2]  Metamorphosis of plasma turbulence-shear-flow dynamics through a transcritical bifurcation. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  W. Fundamenski,et al.  Radial interchange motions of plasma filaments , 2006 .

[4]  Wen-Xu Wang,et al.  Predicting catastrophes in nonlinear dynamical systems by compressive sensing. , 2011, Physical review letters.

[5]  Lennart Ljung,et al.  System identification (2nd ed.): theory for the user , 1999 .

[6]  W. Fundamenski,et al.  Interchange turbulence in the TCV scrape-off layer , 2006 .

[7]  B. Wan,et al.  One-dimensional modelling of limit-cycle oscillation and H-mode power scaling , 2015 .

[8]  Steven L. Brunton,et al.  Inferring Biological Networks by Sparse Identification of Nonlinear Dynamics , 2016, IEEE Transactions on Molecular, Biological and Multi-Scale Communications.

[9]  P. Diamond,et al.  Study of the L–I–H transition with a new dual gas puff imaging system in the EAST superconducting tokamak , 2013 .

[10]  S. Chapman,et al.  Robustness of predator-prey models for confinement regime transitions in fusion plasmas , 2013, 1302.0727.

[11]  J. Madsen,et al.  Simulation of transition dynamics to high confinement in fusion plasmas , 2014, 1409.3186.

[12]  Suppression of turbulence at low power input in a model for plasma confinement transitions , 2005 .

[13]  F. Wagner,et al.  A quarter-century of H-mode studies , 2007 .

[14]  P. Diamond,et al.  Direct identification of predator-prey dynamics in gyrokinetic simulations , 2015 .

[15]  Edward A. Spiegel,et al.  Rayleigh‐Bénard Convection: Structures and Dynamics , 1998 .

[16]  S. Benkadda,et al.  Confinement and bursty transport in a flux-driven convection model with sheared flows , 2003 .

[17]  Lennart Ljung,et al.  System Identification: Theory for the User , 1987 .

[18]  Ö. Gürcan,et al.  Characterization of predator–prey dynamics, using the evolution of free energy in plasma turbulence , 2012, 1205.5712.

[19]  W. Horton,et al.  L-h confinement mode dynamics in three-dimensional state space , 1995 .

[20]  P. Diamond,et al.  Weak hysteresis in a simplified model of the L-H transition , 2009 .

[21]  Urbashi Mitra,et al.  Editorial: Inaugural Issue of the IEEE Transactions on Molecular, Biological, and Multi-Scale Communications , 2015, IEEE Trans. Mol. Biol. Multi Scale Commun..

[22]  X. Garbet,et al.  The quasilinear behavior of convective turbulence with sheared flows , 2003 .

[23]  P. Diamond,et al.  Zonal flows and pattern formation , 2015 .

[24]  O. E. Garcia,et al.  Mechanism and scaling for convection of isolated structures in nonuniformly magnetized plasmas , 2005 .

[25]  J. Rasmussen,et al.  Numerical modeling of the transition from low to high confinement in magnetically confined plasma , 2015 .

[26]  P. Diamond,et al.  Mean shear flows, zonal flows, and generalized Kelvin–Helmholtz modes in drift wave turbulence: A minimal model for L→H transition , 2003 .

[27]  Ö. Gürcan,et al.  Spatio-temporal evolution of the L → I → H transition , 2012 .

[28]  O. E. Garcia,et al.  Structures, profile consistency, and transport scaling in electrostatic convection , 2005 .

[29]  S. Brunton,et al.  Discovering governing equations from data by sparse identification of nonlinear dynamical systems , 2015, Proceedings of the National Academy of Sciences.

[30]  M. Brøns,et al.  Bifurcation analysis and dimension reduction of a predator-prey model for the L-H transition , 2013 .

[31]  J. Juul Rasmussen,et al.  Turbulence and intermittent transport at the boundary of magnetized plasmas , 2005 .