择要

非线性动力学在诸如双摆和振荡电路等本科物理传授教化中扮演着积极的角色。
雅可比矩阵作为非线性动力学中的主要观点，常被用来剖析不动点的性子。
本文中我们先推导了一个二维系统中雅可比矩阵的详细形式，并将其运用到磁约束核聚变约束模式转换模型，末了延伸求解出 ZCD 模型中各不动点的雅可比矩阵本征值。
我们随后调度了模型中的外界加热功率，并利用雅可比矩阵不雅观察模型中非线性动力学性子的变革。
经研究创造，外界加热功率的改变会使得 ZCD 模型极限环的大小发生变革，即功率越大极限环半径越大。
极限环半径的扩大导致环与鞍点相交为一个同宿轨，诱使系统产生同宿分岔。

AbstractNonlinear dynamics plays an active role in university physics such as double pendulum and resonant circuit. As an important concept in nonlinear dynamics, Jacobian matrix is often adopted to analyze the properties of fixed points. In this paper, we derive the specific form of the Jacobian matrix in a classic two-dimensional system, and apply it to the L-H mode transition model of magnetic confinement fusion, then solve the eigenvalues of the Jacobian matrix at each fixed point in the ZCD model. We adjust the external heating power in the model and analyze the nonlinear dynamic properties of the model by the Jacobian matrix. It is found that the variance of external heating power changes the size of limit cycle in ZCD model; and the higher the power is, the larger the limit cycle will be. The expansion of the limit cycle radius leads itself to intersect with the saddle point into a homologous orbit, which induces the system to occur homoclinic bifurcation.

Key wordsnonlinear dynamics; Jacobian matrix; L-H mode transition; homoclinic bifurcation

非线性动力学（nonlinear dynamics）是研究系统中各种运动状态的定量和定性规律的科学[1-5]。
非线性动力学系统是由一个或一组非线性微分方程描述的随韶光变革的系统。
若想要预测一个别系的运动轨迹，我们首先须要给出它在眇小韶光尺度里的性子并列出动力学方程[6]。

人们对非线性问题的认识始于 17 世纪。
1673 年惠更斯创造了非线性征象[7]。
在 1687 年牛顿创造了运动定律和万有引力之后，数学家和物理学家开始考试测验将其推广到三体问题，但是并没有取获胜利。
直到 19 世纪末，法国科学家庞加莱创立了微分方程的定性理论，为非线性动力学的发展奠定了根本[8]。
他首次创造了系统的混沌行为，混沌是一种确定性系统（deterministic system）展示出敏感依赖于初始条件的非周期行为。
1963 年洛伦兹在研究大气对流模型时，创造了系统蝴蝶形状的运动轨迹，这标志着混沌理论的出身[9]。
混沌理论广泛运用于各个领域，科学家们也因此发展了很多数学方法，在这个过程中 Duffing 方程、Van der Pol 方程、Mathieu 方程等著名的数学模型相继被建立，它们至今仍被人们用来研究非线性系统动力学征象的实质特色[10]。
从 20 世纪 60—70 年代开始，非线性动力学理论成为一门主要的前沿学科，分岔和混沌的研究随之成为新的研究热点。
非线性动力学的思想与方法已被运用到很多的新领域，比如系统生物学、经典力学和社会物理学等。

非线性方程的求解方法相较线性方程更加困难。
定性方法和定量方法是研究非线性动力学的两条路径，二者缺一不可。
定量方法又分为解析法和数值法。
常日情形下，我们很难得到非线性方程的解析解，以是只能求解确定参数时的数值解。
定性方法是利用几何的思维办法，从非线性微分方程入手，利用相图（phase plot）来描述微分方程。
在相图中，我们可以研究不动点的稳定性，并预测系统随韶光蜕变的行为。
在剖析详细问题时，定量地去度量稳定性是很有必要的。
这须要我们在方程的不动点附近，将非线性方程线性化。
在一维向量场中的不动点处，图像的斜率决定了不动点的稳定性。
对付多元系统，则须要求解其雅可比矩阵来剖析不动点的稳定性。
据笔者理解，同时采取定性和定量剖析方法研究详细多元系统非线性动力学性子的中文论文并不多。
以是本文通过剖析磁约束核聚变模式转换中的 ZCD 模型，来做一个传授教化式的推导以期回顾传授教化并延伸到科研。
本文框架如下：

（1） 推导二维非线性动力学系统雅可比矩阵的详细形式。

（2） 大略先容磁约束核聚变装置托卡马克中的L-H约束模式转换模型。

（3） 通过调度ZCD模型中外界加热功率的参数，以实现核聚变约束模式的转换，并利用相图和雅可比矩阵剖析和不雅观察模型中各不动点非线性动力学性子的变革。

1 雅可比矩阵的详细形式

我们考虑一个经典的二维非线性系统

在相平面上同时知足f(m0,n0)=0，g(m0,n0)=0的点（m0,n0）被称为不动点（fixed point）。
我们只要打算出不动点邻域内的轨迹分布及其趋于稳定不动点的衰减率，就可以判断其稳定性。
为此我们设u（t）=m－m0，v（t）=n－n0为不动点附近的微扰，并将其利用泰勒展开可得到微分方程

个中，O（u⊃2;，v⊃2;，uv）为u⊃2;、v⊃2;，uv 的二阶小量。
由于 f（m0，n0）=0，以是

同理可得

因此,对付小的扰动,我们可以得到如下的关系式

个中，矩阵

称为不动点处的雅可比矩阵。

由于 Ο（u2,v2,uv ）为趋近于0的变量,以是 我们得到如下的线性系统

这样我们就可以通过求解雅可比矩阵的本征值和本征向量来判断不动点处的增长趋势，随之判断出不动点的稳定性及系统动力学轨迹的蜕变性子。

为了研究雅可比矩阵的性子，我们设 λ1 和 λ2（λ1≠λ2）为矩阵 J 的本征值，a 和 b 为本征向量，那么在任意初始条件下，我们可以将通解写为

个中，本征值在 rr 的指数位置。
若两个本征值都为负数，则解在两个本征向量方向上指数减少，不动点是稳定节点（stable node）；若两个本征值都为正数，则解在两个本征向量方向上指数增长，不动点是不稳定节点；若本征值一个为正数，一个为负数，则不动点称为鞍点（saddle point），即当 t→∞ 时，轨迹趋近于不稳定流形，而当 t→－∞ 时，轨迹趋近于稳定流形。
若本征值为复数，不动点或为中央（center）或为焦点（focus）。
在本文所研究的 ZCD 模型中，雅可比矩阵有四个本征值。
我们规定当三个本征值为负数，一个为正数时，称为一类鞍点；两个本征值为负数，两个为正数时，称为二类鞍点。
对付个中存在某一个本征值为正数和负数的不动点，我们称其在对应的本征向量方向上表现为源和汇的性子。

2 磁约束核聚变 L-H 转换模型

核聚变是两个轻原子核结合组成一个较重原子核的过程，而个中的 L-H 约束模式转换则是受控核聚变中一个主要的研究分支[11]。
实现核聚变反应须要达到 2 亿摄氏度的高温，在高温下物质表现为等离子体态，等离子体在高温下剧烈运动，以是这对聚变反应堆的设计有着极高的哀求。
科学家们因此构建了多种约束等离子体的装置，个中利用磁场来约束等离子体的托卡马克装置是最为有效的。
托卡马克装置的反应事理为氢元素的同位素氘（21 H 或 D）与氚（31H 或 T）聚合反应天生氦（42He）。
实现核聚变反应须要知足劳森判据（Lawson Criterion）nTeτE＞1021m－3sKeV，个中 n 为等离子体密度，Te 为等离子体温度，τE 为能量约束韶光。
虽然根据爱因斯坦质能方程打算，很少的反应物便可以开释大量的能量，但是目前在实验上等离子体的约束韶光是非常短的，以是人工受控核聚变的实现须要更高的约束水平和更长的约束韶光。

托卡马克装置是由两个紧张磁场合成一个螺旋磁场，使等离子体绕磁力线做拉莫尔回旋运动。
个中，一个磁场是由缠绕在圆环形真空室表面的铜线圈产生的环向磁场，另一个是由等离子体中的带电粒子电流产生的极向磁场，如图1所示。
在非均匀磁场下等离子体会发生漂移[13]，等离子体之间的相互浸染会产生能量输运过程。
由等离子体湍流引起的输运过程称为反常输运。
等离子体边缘湍流较强使得温度梯度较弱，会丢失大量能量，这导致托卡马克处在一个较低的约束水平，即低约束模式（L 模式）。
在 1982 年，ASDEX 托卡马克装置加热期间能量约束问题得到了大幅改进。
Wagner 等人创造在L模的根本上通过提高加热功率，并达到一定的阈值之后会实现 L-H 模式转换，H 模式为高约束模式[14]。
H 模式的边缘区域由于密度和温度溘然增长而涌现边缘输运垒，整体约束水平得到极大的提升。
因此 H 模式是未来反应堆的空想约束模式。

L-H 模式的转换是由E×B剪切流对等离子体湍流的抑制产生的[15]。
剪切流包含有均匀流、带状流以及湍流和压强非线性耦合后产生的测地声模，其对付湍流的抑制浸染是研究的重点。
Kim 和 Diamond 在 2003 年建立了包含等离子体湍流、均匀流、带状流和压力梯度四个变量的零维模型，研究创造通过缓慢增加输入功率改变压力梯度可以很好地抑制湍流，使系统从L模通过振荡行为蜕变到静止的H模。
带状流的紧张浸染是抑制湍流，而均匀流会导致振荡的持续韶光变长。
振荡行为是由于带状流和湍流之间的相互竞争，这符合非线性动力学中经典的猎物捕食者种群竞争模型。
在 2009 年，Malkov 和 Diamond[16]构建了包含等离子体湍流、带状流和温度梯度的模型，个中外界加热功率 q 被包含在温度梯度的表达式中，并作为掌握参数会引发系统的四种状态，分别为 L 模式、T 模式、H 模式、QH 模式。
当 q 超过一定阈值时，T 模式会经历 Hopf 分岔，由焦点变为稳定的极限环。
2013 年 Zhu [5]等人将方程拓展为四变量模型，即 ZCD 模型。
在MD 模型的根本上，引入了第二种非线性构造，即测地声模 GAM，用 U2 表示，见公式（12）。
因此，方程变为两个捕食者，一个猎物的模型。
带状流 U1 和测地声模 U2 都会抑制湍流E的增长，使系统达到高约束模式。

ZCD模型方程为

这是一个四变量耦合的非线性微分方程组，该模型中四个变量分别为等离子体湍流强度 E、温度梯度 N、带状流 U1 和测地声模 U2，外界加热功率 q 为掌握因子。
我们通过前述的方法，可以将系统线性化得到对应的雅可比矩阵 J。

3 外界加热功率引起的同宿分岔

本节我们通过改变外界加热功率 q，用非线性动力学的方法对 ZCD 模型中极限环的性子做进一步研究。
图 2~图 6 为 q 值不同时温度梯度 N、等离子体湍流强度 E、带状流 U1、测地声模 U2 以及加热功率 q 随韶光变革的图像（左侧）和系统在相空间（E, U2, N）中的韶光蜕变（右侧），个中箭头表征蜕变的方向。

图 2 中 q 值为 0.47，其图像包括了四个阶段，分别为 H 模式、L 模式、T 模式和 O 模式。
四个模式分别对应到相空间中我们所求解的四个不动点 A、B、C、D，在图 1（b）中已标出。
我们在图 1（a）中不雅观察到，在 t=6500 附近带状流 U1 逐渐消逝，U2 随之增长，湍流 E 被抑制，从而涌现了一个非常稳定的极限环。
在极限环状态下 U1 为零，E、N、U2 不为零。
表 1 中显示不动点 D（极限环）处本征值包含一对共轭复数，且实部为正数，轨迹螺旋运动至极限环。
实际的实验中系统会通过极限环振荡模式达到 H 模，H 模式为 N 更高的状态。
为研究掌握因子 q 对系统约束水平的影响，我们在图 3~图 6 中改变了外界加热功率 q，其值分别为 q=0.49、q=0.52、q=0.54 和 q=0.549。

我们创造在图 2 至图 5 显示的变量随韶光变革的图像（a）中，随着 q 值的增大，轨迹在第一个不动点 A（H 模）的持续韶光明显增长，q 由 0.47 增至 0.54时，t 由约 400 增至约 2500。
但是轨迹在 B 点（L 模）和 C 点（T 模）处的持续韶光险些不变，即H模对付 q 值的大小更为敏感。
q 值的增大直接导致 N 值增大，这导致带状流 U2 被抑制，在（E,N）平面中 U2 的最大振幅减小，且振荡频率增加。
除此之外，最明显的变革是极限环的周期和半径，它们更加依赖于外界加热功率 q 值的大小。
详细表现为 q 值越大，极限环的周期越长，半径越大。
极限环周期的延长使得系统处在 H 模的韶光更长，即约束韶光更长；极限环半径的增大导致轨迹向（E,N）平面和（U2,N）平面靠近。

极限环周期的变革详细表示为，U2 的振幅变革范围越来越大，最小值越来越趋近于 0，湍流便不能被很好的抑制，在每个周期中湍流存在连忙增长和连忙低落的情形。
在图 5 中，伴随着 U2 和 E 之间的相互浸染，极限环与鞍点 A 和鞍点 B 相重合，系统发生了同宿分岔。

通过数值打算，我们在表 1 中看到当 q=0.54 时，雅可比矩阵的本征值反响不动点 D 是一个鞍点，而并非图 5（b）中的极限环。
在同宿分岔中，随着系统中参数变革，极限环与鞍点可以彼此越来越近直至相交。
环与鞍点在相交后成为一个同宿轨道。
在图 2（b）至图 5（b）相图中我们可以看到在 q 值在达到 0.54 之前，极限环由于周期延长越来越靠近 A、B 这两个鞍点。
当 q 值增加至 0.54 时，极限环膨胀且爆裂为鞍点，同时产生一条同宿轨（见图 5），q=0.54 即为系统产生同宿分岔的一个阈值。
由于鞍点不稳定的性子，当 q 大于 0.54 时，鞍形连接会分裂，导致系统无法回到同宿轨道，即环被摧毁。
同宿分岔是全局分岔的一种类型，其会导致相空间中轨迹的拓扑构造发生变革，而且这种变革不像局部分岔那样局限在一个小的邻域内，拓扑构造的变革可以延伸到任意大的间隔。
如图 6（b）所示，全体系统形成了一条同宿轨道。
在图 6 中，同宿分岔导致 U1 周期振荡，湍流被更好地抑制，约束水平达到最高。

4 结语

在本科物理传授教化中会碰着一系列非线性微分方程的求解和剖析问题，本文回顾并利用了非线性动力学中的雅可比矩阵来剖析磁约束核聚变的约束模式转换模型。
雅可比矩阵的本征值反应了不动点附近变量的衰减率，四维方向上的衰减率共同决定了不动点的性子，这为我们预测模型随韶光的蜕变供应了数值依据。
在ZCD 模型中，我们创造外界加热功率 q 值会影响到极限环的周期，q 值越大极限环的周期越长。
极限环半径的增大导致了极限环与鞍点相重合，诱使系统产生同宿分岔，涌现一条同宿轨。
在产生同宿分岔时，极限环状态的不动点 D 与鞍点 A 非常靠近，导致本征值的打算结果与所绘制相图中不动点的性子不符合，雅可比矩阵可以部分阐明磁约束核聚变中 L-H 转换模式。
在物理传授教化中该当把稳到求解雅可比矩阵的本征值不能准确地描述轨迹的性子，须要结合相图详细剖析。

参考文献

[1]刘锦,张孟.弹簧车摆的动力学仿照研究[J].物理与工程,2021,(5):93-97.

LIU J, ZHANG M. Simulation study on the dynamics of spring car pendulum[J]. Physics and Engineering, 2021, 31(5): 93-97. (in Chinese)

[2]李万祥,何玮,唐恭佩.一类复摆系统的非线性动力学研究[J].华中科技大学学报:自然科学版,2007,(5):27-30.

LI W X, HE W, TANG G P. Research on nonlinear dynamics of a double pendulum system[J]. Journal of Huazhong University of Science and Technology (Nature sci ed), 2007, 35(5): 27-30. (in Chinese)

[3]张连芳,傅敏学,刘滢滢,等.混沌实验传授教化之路[J].物理与工程，2013,(1):21-24.

ZHANG L F, FU M X, LIU Y Y, et al. The way of education for chaos experiment[J]. Physics and Engineering, 2013, 23(1): 21-24.（in Chinese）

[4]李明达,董乔南,杨亚利,等.利用非线性动力学系统研究混沌征象[J].物理与工程，2019,(6):77-84.

LI M, DONG Q N, YANG Y L, et al. Study of chaos with nonlinear dynamical system[J]. Physics and Engineering, 2019, 29(6): 77-84.（in Chinese）

[5]ZHU H, CHAPMAN S C, DENDY R O. Robustness of predator-prey models for confinement regime transitions in fusion plasmas[J]. Physics of Plasmas, 2013, 20(4): 042302.

[6]STROGATZ S H. 非线性动力学与混沌[M]. 孙梅, 汪小帆等, 译. 北京: 机器工业出版社: 2016.

[7]周道其. 科学家揭开惠更斯复摆钟之谜[J]. 当代物理知识, 2002, 14(5):1.

ZHOU D Q. Scientist shave solved the mystery ofHuygens compound pendulum clock[J]. Review of World Invention,2002,14(5): 1. (in Chinese)

[8]陈明晖,邓明立.庞加莱微分方程定性理论研究初探——兼纪念庞加莱诞辰150周年[J].科学.2004,(1):29-31.

CHEN M H, DENG M L. A Preliminary study on the Qualitative Theory of Poincarés differential Equation—Also commemorating the 150th anniversary of Poincarés birth[J]. Science, 2004, 56(1): 29-31. (in Chinese)

[9]LORENZ E N. Deterministic nonperiodic flow[J]. Journal of the Atmospheric Sciences, 1963,20(2):130-141.

[10]叶美盈.一种非线性动力学系统数值仿照的新方法[J].云南师范大学学报:自然科学版,2000,(2):33-36.

YE M Y. A new approach to the numerical simulation of nonlinear dynamics[J]. Journal of Yunan Normal University (Natural sci ed), 2000, 20(2): 33-36. (in Chinese)

[11]MCCRACKEN G H, STOTT P.宇宙能源——聚变[M]. 核工业西南物理研究院翻译组. 北京:原子能出版社，2008.

[12]LI S, JIANG H, REN Z, et al. Optimal tracking for a divergent-type parabolic PDE system in current profile control[C]//ZHANG X G. Abstract and Applied Analysis. New York: Hindawi Publishing Corporation, 2014: 1-8.

[13]刘爱红,佘守宪.轴对称缓变磁场中的等离子体与磁约束事理[J].物理与工程，2002,(5):14-17.

LIU A H, YU SH X. Plasma in slow-varying magnetic field with axial-symmetry and principle of magnetic confinement[J]. Physics and Engineering, 2002, 12(5): 14-17. (in Chinese)

[14]WAGNER F. Regime of Improved Confinement and High Beta in Neutral-Beam-Heated Divertor Discharges of the ASDEX Tokamak[J]. Physical Review Letters, 1982, 49(19): 1408-1412.

[15]BURRELL K H. Effects of E × B velocity shear and magnetic shear on turbulence and transport in magnetic confinement devices[J]. Physics of Plasmas, 1997, 4(5): 1499-1518.

[16]MALKOV M A, DIAMOND P H. Weak hysteresis in a simplified model of the LH transition[J]. Physics of Plasmas, 2009, 16(1): 012504.

基金项目: 北京市自然科学基金帮助项目(1194020，2232045)、北京市教委科技操持一样平常项目(KM201910005001)、国家重点研发操持(2018YFB0703500)。

通讯作者: 朱浩，男，北京工业大学、北京-都柏林国际学院讲师 ，紧张从道理论物理的研究事情，研究方向为等离子体物理与核聚变，hao.zhu@bjut.edu.cn。

引文格式: 艾媛媛，朱浩. 非线性动力学在磁约束核聚变约束模式转换中的运用[J]. 物理与工程，2023，33（5）：95-101,107.

Cite this article: AI Y Y, ZHU H. Application of nonlinear dynamics to L-H mode transition in magnetic confinement fusion[J]. Physics and Engineering, 2023， 33（5）：95-101, 107. (in Chinese)

END

更多精彩文章请点击下面“蓝字”标题查看：

CALL FOR PAPERS｜《物理与工程》征稿量子世纪年中国大学物理教诲MOOC同盟2023年事情会议（扩大）暨范例案例互换会 会议纪要2023年全国高档学校物理根本课程教诲学术研讨会 会议纪要2023 年全国高档学校物理根本课程青年西席讲课比赛在喀什大学举办全国大学物理实验传授教化对口增援（智力援疆）研讨会在新疆师范大学举办王青教授：理解王中林院士“拓展的麦克斯韦方程组”“碰瓷”麦克斯韦：伽利略协变和洛伦兹协变电磁场论趣谈热点：运动介质洛伦兹协变电磁理论2021年《物理与工程》精良论文、精良审稿专家、精良青年学者名单王青教授：源自苏格拉底的问题驱动式教诲——在互动中共同学习和发展读后感：教诲中的现实和远方王青教授：昨晚（6月9日），清华电动力学期末考试朱邦芬院士：“减负”误区及我国科学教诲面临的寻衅《物理与工程》2023年第4期目录乐永康：新冠肺炎疫情防控下美国物理实验传授教化及中美情形比拟顾牡：对付重新制订的《非物理类理工学科大学物理课程传授教化基本哀求》的认识和体会朱邦芬院士：从根本科学班到清华学堂物理班朱邦芬院士：对培养一流拔尖创新人才的思考李学潜教授：物理是一种文化李学潜教授：如何帮助物理系学生迈过从高三到大一这个坎穆良柱：物理课程思政教诲的核心是科学认知能力培养穆良柱：什么是物理及物理文化？穆良柱：什么是ETA物理认知模型穆良柱：什么是ETA物理传授教化法吴国祯教授：我的国外研究生经历印象——应清华大学物理系“基科班20年·学堂班10年纪念活动”而写

陈佳洱，赵凯华，王殖东：面向21世纪，急待重修我国的工科物理教诲王亚愚教授：清华物理系本科人才培养理念与实践葛惟昆教授：关于中外人才培养的几点思考安宇教授：为什么传统的教室讲授模式须要改变安宇教授：实在传授教化便是积累的过程刘玉鑫教授：关于本科生物理根本课程传授教化和教材编著的一些思考沈乾若：重创理科教诲的美加课程改革Henderson C：美国研究基金支持下的物理教诲研究及其对高档物理教诲的影响《物理与工程》期刊是专注于物理教诲传授教化研究的学术期刊，是中国科技核心期刊，1981年创刊，欢迎踊跃投稿，期刊投审稿采编平台：

http://gkwl.cbpt.cnki.net

欢迎关注

《物理与工程》微信"大众年夜众号

连续滑动看下一个轻触阅读原文

物理与工程向上滑动看下一个

原标题：《非线性动力学在磁约束核聚变约束模式转换中的运用》

阅读原文

来源：物理与工程

(本文来自澎湃***，更多原创资讯请***“澎湃***”APP)