虽然徐云记不太清他提出SImPLE改进算法的具体时间，但苏哈斯·帕坦卡可不是什么年少成名的天才。
他想要SImPLE改进算法，提出无论如何也要到十多年以后了。
不夸张的说。
这年头整个数学界和物理学界对于纳维-斯托克斯方程的研究，还处在一个非常原始的状态。
就连SImPLE算法.....也就是求解压力耦合方程的半隐式方法的最初版本，都要在1972年才会被提出。
想到这里。
徐云便决定小小的帮叶笃正一把——虽然他之前确实没有这方面的打算。
但这种能够让兔子赶上甚至反超第一梯队的事儿，他自然还是很乐意为之的。
反正不要钱，多少试一点嘛。
随后徐云顿了顿，飞快的在脑海中组织了一番思路，对叶笃正说道：
叶主任，我的意思是在这个变式后加个伯努利函数，然后再取个旋度，您觉得可行吗？
这是我在剑桥大学那会儿听一位学长说的，当时他们推导的情景恰好也是相同的变式.......
唰——
结果徐云话没说完。
叶笃正便低头在纸上写下了一个函数：
C=p/+u2/2。
这个函数来自等式?(u2/2)=(u??)u+u，也就是伯努利函数。
接着叶笃正又按照徐云的说法取了个旋度，得到了一个新的公式：
?/?t=?[u]+v?2。
别看这个公式瞅起来跟颜文字似的，好像又是(￣▽￣)~*(￣▽￣)／又是~(￣▽￣)~*。
对于叶笃正而言。
在见到它的一瞬间，他的心脏便狠狠漏跳了一大拍！
这是......
的演化方程！
同时由于?(u)=(??)u?(u??)的缘故，所以这个演化方程还可以改写为对流导数的形式：
d。
写到这里。
叶笃正再次一停顿，扭头又
看向了徐云，迫不及待的问道：
韩立同志，后面呢？后面的思路是什么？
此时此刻。
叶笃正仿佛回到了自己在芝加哥读书的日子。
当时他在追一本连载于芝加哥日报的推理，每每看完一章时便迫不及待的想要疯狂进行催更。
如果不是怕失去留学海外的宝贵资格。
叶笃正甚至考虑过要不要把作者绑到小黑屋去更新——一天必须要更新个五万字，要不然当天不能吃饭！
而在他对面。
徐云则示意乔彩虹将自己的轮椅再朝叶笃正靠近了一些。
随后他从叶笃正手中接过纸和笔，一边写一边解释道：
叶主任，这个方程想要继续推导下去，首先就要明白这个变式的物理意义。
我们在这里再导入一个角动量方程做个对比...你看，物理意义应该就很明显了吧？
叶笃正认真看了小半分钟，很快哦了一声：
哦，我懂了。
右边描述的是因为流体元的拉长，体元惯量矩的改变，还有就是粘性力矩作用在体元上，没错吧？
徐云点了点头。
这个变式的物理意义，差不多可以算是后世涡度的入门级概念。
也就是流体块的涡度可能因为它的拉长而改变，引起惯量矩的改变，或者因为粘性应力加速或者减速。
紧接着。
徐云又写了个佩克来数。
也就是Pe=ud/α，又在上头换了个圈，带入回了原式。
看到这里。
叶笃正的鼻翼中忽然传出了一声带着意外的鼻音，眉头骤然一扬。
他发现了一个此前从未意识到的问题：
根据变式来看。
二维流中涡度是对流，并且像热量一样可以扩散，那么关于佩克来特数的类比就是.....
Re=u?/v。
这意味着涡度像热量一样，在二维流内部不能凭空产生或毁灭。
并且它可以通过对流从一个地方移动到另一个地方。
但另一方面。
∫dV对于所有定域的涡度团是守恒的。
也就是说......
漩涡通过速度场对流，通过扩散传播，但是每个漩涡内总的涡度保持不变。
换而言之.....
边界正是涡度的来源！
这是一个叶笃正从未想过的概念，这代表着他之前的很多思路都是错误的，他确实低估了边界的深度。
但这也同样代表着.....
一个新模型的可能！
准确来说应该是......