他原本以为叶清河只是在原有一阶模型基础上,做了些高阶补项,微调收敛结构。
作为已经跟着教授快四年的老人,深耕CFD数值这么久,他已默认强激波震荡间断维解,必须人工粘性兜底,这是高超计算不可突破的铁律,是整个行业30年来公认的物理代价,是所有数值格式差分架构无论怎么优化都绕不开的工程死结。
可是,当他逐行扫完这套双空间拓扑介偶体系弱解延拓边界约束、零人工粘性守恒闭环证明,再对照超算输出的全域无震荡、无数值色散、全时空步弱解唯一存在的核验数据时,整个人彻底失神。
瞳孔收紧,呼吸下意识放轻,他伸手快速地翻着报告。
他翻遍报告里的迭代残差曲线、跨子空间通量守恒记录、激波间断处解场连续性判定,却越看越是后背发麻。
目前业内所有的高超CFD地逻辑全是在用数值耗散换稳定,用精度换收敛。
全世界的工程范式都是,激波间断必崩适定性,必生伪解,必靠人工粘性强行抹平,别无二法。
但是叶清河的这份报告,直接从方程存在性层面废掉了这套行业共识。
不需要工程抹平,不需要数值妥协,不需要牺牲激波层精度,直接通过函泛空间重构,让原本无解、不适定、必畸变的强激波间断区域,拥有了严格唯一、全局守恒、完全自洽的弱解体系。
这根本不是优化模型!!
段博文喉咙微微发紧,心底掀起滔天海啸,他抬头看了看一脸平静的叶清河,这么大的事情,他怎么就能够这么平静呢?
这可是推翻了高超声速可压缩流数值求解的底层前提啊!
这要是换任何一个人,别是他,就算是他的导师符嵩,估计也早就兴奋得跳了起来,绝不可能像叶清河这么平静,就好像做了一件可有可无的事情一样。
这几年,他做了太多的流体代码地,调试过上百套差分格式,试过数十种激波捕捉算法,太清楚数值震荡永远是离散格式的固有原罪这个道理了。
这是连续方程转离散迭代天生带的缺陷,是工程人员只能无限修补却永远无法根除的顽疾。
他从来没有怀疑过这个问题。
也从来没想过这个问题是可以被解决的。
就算是他的导师在讲解这个问题的时候,也是一样的态度。
认为这个是绕不开的难题。
可是现在摆在他面前的超算结论,清清楚楚地告诉他,这个问题被叶清河,被这个只是短短接触流体学没多久的叶清河给解决了。
在这一刻,他第一次清晰地认识到,导师在介绍叶清河叶教授的时候自己不如他的那个心情。
也真真正正明白为什么导师在提起叶教授的时候,会是那样的态度。
别人研究了一辈子认为根本无法绕开的问题,他就这么轻松地绕过去了,就这么轻松地解决了。
重构的不是代码,不是格式,也不是参数。
叶教授重构的是NS方程在强间断流场下的解的存在逻辑。
别人都是在离散层面不断地修补误差,他直接在连续理论层面消弭误差根源。
最让他头皮炸裂的是,这套理论完全自洽,超算极限工况核验通过,数学闭环无懈可击。
可以理论层面已经完美无缺,现在交到他手上的任务,仅仅是把绝对真理翻译成计算机能读懂的离散语言。