第512章 心都要碎了(1 / 2)
请关闭浏览器的阅读/畅读/小说模式并且关闭广告屏蔽过滤功能,避免出现内容无法显示或者段落错乱。
“顾老师,第五章这个非交换几何的谱三元组构造,如果把它应用到磁流体系统的非定常边界上,是不是可以替代我们现在用的摄动展开?
摄动展开截断到四阶之后,长时间积分的累积残差一直是我们的痛点,上次肖教授提过可以用李雅普诺夫函数做误差界估计,但我们在行星引力场边界处的重构频率还是太高,算力消耗有点吃不消。
我看肖教授在这一章里用谱三元组处理了非交换边界上的正则化问题,本质上就是把边界的不确定性映射到了一个商空间上,然后通过算子的谱分解来消除冗余自由度。
这个思路如果能嫁接到我们的磁流体模型上,说不定能把摄动展开整个绕过去。”
顾清尘面前的屏幕上,俞巍的脸被分割成两个窗口,一个对着他,另一个共享着一份《数学原理》的电子版,页面停在第五章中间。
俞巍眉头微微拧着,显然已经在这一段上琢磨了好一阵子。
“谱三元组的构造确实可以处理非交换边界,不过,小宿在第五章里用的是静态边界条件。
你们的磁流体系统,边界是动态的,等离子体壁面相互作用本身就是非定常的。
直接把静态构造搬过去,算子的谱分解在时间演化下会失稳。”
他翻到第五章的后半部分,指着一处俞巍刚才没有标注的段落。
“你看这里,他在处理边界条件时引入了一个辅助参数来保证算子的自伴性。这个辅助参数在静态边界下是常数,但是你的边界是时变的,常数假设不成立。
你需要把辅助参数改造成一个动态泛函,让它跟着边界条件一起演化。
动态泛函的形式可以从第六章的和乐群轨道等价类构造里找,那一章专门讨论了非定常边界下的商空间保持问题。”
俞巍快速翻到第六章,目光在顾清尘说的那几页上来回扫了两遍,若有所思地点了点头,随即又问道:
“那重构频率的问题呢?动态泛函能降低重构频率吗?”
“不能直接降低,但是可以让重构变得更平滑。”
“你现在用摄动展开,每次重构相当于在整个商空间上重新做一次截断,截断误差累积就是你们算力消耗的主要来源。
如果用谱三元组加动态泛函,重构就不再是重新截断,而是沿着和乐群轨道做一个连续的等价类滑移。
滑移的计算复杂度比截断低一个数量级,你的算力消耗自然就下来了。”
俞巍的眼睛亮了一下,随即又闪过一丝狡黠的光,追问道:
“那动态泛函的具体形式呢?第六章给出的是一般构造,具体到磁流体系统上,边界条件的时变规律怎么嵌入泛函的参数空间?”
顾清尘看着他,笑了一声。
这个问题问得很有水平,俞巍只是出去历练了几个月,没想到成长了这么多。
“这个问题要看你具体的边界条件了。
你得先把磁流体在等离子体壁面附近的时变规律做一个傅里叶谱分解,找出主频分量,然后再把主频分量作为参数嵌入泛函的核函数里。