。 不夸张的说。 这年头整个数学界和物理学界对于纳维-斯托克斯方程的研究,还处在一个非常原始的状态。 就连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对于所有定域的涡度团是守恒的。 也就是说...... 漩涡通过速度场对流,通过扩散传播,但是每个漩涡内总的涡度保持不变。 换而言之..... 边界正是涡度的来源! 这是一个叶笃正从未想过的概念,这代表着他之前的很多思路都是错误的,他确实低估了边界的深度。 但这也同样代表着..... 一个新模型的可能! 准确来说应该是...... 气象学中第一个真正可行的新模型! 要知道。 虽然挪威学派在数值天气预测这方面贡献很大很大,但即便是到现在,整个气象行业也依旧没有一个真正的模型。 事实上。 按照正常历史发展。 气象学要到1971年才会由拉苏尔建立出第一个气候模型。 并且拉苏尔建立的模型预测的还不是局部天气,而是与全球变暖有关的气候模型。 而眼下...... 叶笃正的面前出现了一条新路。