““
此时此刻。
看着身边出现的这张脸。
毫无防备之下。
叶笃正险些没把手伸进衣兜,掏出自己母亲当年送给自己的本命木牌朝对方甩过去——据说那玩意儿是桃木做的,能驱邪。
当然了。
在掏出木牌之前,叶笃正便先一步反应了过来。
出现在自己面前的这货并不是木乃尹,而是最近一段在基地小有名气的
熟人韩立。
见此情形。
叶笃正在暗自松了口气的同时,也用缓慢的语速掩盖了自己内心最开始的惊恐,开口道:
“韩立同志,我还以为是谁呢,原来是你啊”
说实话。
叶笃正对于徐云的印象其实还是很不错的。
毕竟若是没有他拿出的气象多普勒雷达,气象中心恐怕永远都不会有再次自我证明的机会。
倘若真是如此
那么可以预见,整个气象中心将会在很长的时间里失去斗志。
同时从行业角度来看。
气象多普勒雷达对于整个气象学的帮助也显而易见。
这种设备的出现,很可能为一直看不到未来的气象领域开拓出一条全新的康庄大道——还是华夏占据主导权的那种。
所以无论是从本职工作还是个人情感出发,叶笃正对于徐云的印象都很不错,甚至还带着一丝感激。
因此在被吓了一大跳后。
叶笃正也没表露出丝毫不满,而是笑着对徐云问道:
“韩立同志,你怎么到这儿来了,对了,吃过晚饭了吗?”
此时距离第一批数据出炉已经过去了七八个小时有余,天色早已从白天变成了夜晚,再过几个小时差不多就到十二点了。
就在不久前,基地上还托人送来了晚饭。
“嗯,刚刚喝了些粥。”
徐云朝帐篷外的某个方位歪了歪头,此时依稀可以看到几位副业队员正在忙活着发晚餐。
不过今晚的“标餐”规格并不高,大多都是窝头土豆配上榆树叶的菜叶汤。
窝头硌牙,榆树叶发苦。
徐云能喝到精米粥,主要还是和他病人的身份有关系,恢复期需要调养。
接着徐云又把目光放到了叶笃正的算纸上,认真看了几眼:
“咦?叶主任,这是斯托克斯方程组的变式?”
叶笃正微微一怔,看起来对徐云能够认出方程的内容有些惊讶。
不过他很快便想起了徐云的身份,轻轻点了点头:
“对,正是n-s方程组,在涡度的基础上做了一点改变。”
按照老郭此前的介绍。
徐云是剑桥大学数学系的毕业生,认得出n-s方程组倒也正常,毕竟这个方程可是数学领域的一大难题来着。
或者换个角度说。
以徐云能够拿出气象多普勒雷达理论的能力而言,他认不出n-s方程组才怪呢。
徐云则又转头看了眼噼里啪啦满是算盘声的现场,随口对叶笃正问道:
“叶主任,您现在的进度怎么样了?”
“进度?”
叶笃正抬眼与徐云对视了一秒钟,摇了摇头,嘴角扯出一丝苦笑:
“哪有进度?韩立同志,你现在看到的这些就是全部了——后头该怎么推导我自个儿都不知道呢。”
叶笃正说罢。
手指还捏着圆珠笔前半部分笔头做了个小杠杆,将笔尾在算纸上啪啪的拍了两下,看得出来有些烦闷。
气象数据的计算环节不算什么机密,所以叶笃正倒也没想瞒着什么。
毕竟人都是有倾诉欲的。
接着叶笃正便叹息的摇了摇头,准备老老实实的换个思路——既然他考虑的这种变式没有可行性,那么就只能按照竺老给的方案去计算了。
即便
那个想法大概率存在某些问题。
而就在叶笃正提笔书写之际,他的耳边忽然传来了徐云弱弱的声音:
“叶主任,我有个想法啊在这个变式后面加个伯努利函数,您觉得可行吗?”
叶笃正已经写下了几个字母的笔尖顿时一停。
片刻过后。
叶笃正满是诧异的抬起头,一脸见了木乃尹似的表情看着徐云:
“韩立同志,尼(第四声↓)说嘛?”
情绪激动之下。
叶笃正甚至冒出了老家津门的口音。
而在他对面。
看着眼睛瞪得滚圆的叶笃正,徐云的内心其实同样有些意外——他还以为现在定域分布涡度的概念已经比较完整了呢。
不过很快,他便迅速反应了过来。
也是。
对流-扩散方程的关键人物是苏哈斯·帕坦卡,而此君按照年龄来算,现在才二十岁出头呢。
虽然徐云记不太清他提出sImPLe改进算法的具体时间,但苏哈斯·帕坦卡可不是什么年少成名的天才。
他想要sImPLe改进算法,提出无论如何也要到十多年以后了。
不夸张的说。
这年头整个数学界和物理学界对于纳维-斯托克斯方程的研究,还处在一个非常原始的状态。
就连sImPLe算法也就是求解压力耦合方程的半隐式方法的最初版本,都要在1972年才会被提出。
想到这里。
徐云便决定小小的帮叶笃正一把——虽然他之前确实没有这方面的打算。
但这种能够让兔子赶上甚至反超第一梯队的事儿,他自然还是很乐意为之的。
反正不要钱,多少试一点嘛。
随后徐云顿了顿,飞快的在脑海中组织了一番思路,对叶笃正说道:
“叶主任,我的意思是在这个变式后加个伯努利函数,然后再取个旋度,您觉得可行吗?”
“这是我在剑桥大学那会儿听一位学长说的,当时他们推导的情景恰好也是相同的变式”
唰——
结果徐云话没说完。
叶笃正便低头在纸上写下了一个函数:
c=p/ρ+u2/2。
这个函数来自等式?u2/2=u??u+u×w,也就是伯努利函数。
接着叶笃正又按照徐云的说法取了个旋度,得到了一个新的公式:
?w/?t=?×+v?2w。
别看这个公式瞅起来跟颜文字似的,好像又是 ̄▽ ̄~* ̄▽ ̄/又是w。
而在见到这道算式的时候。
徐云裹在绷带下的表情也随之一松。
呼
他的任务算是完成了
想必聪明的同学也看出来了。
没错!
叶笃正此时写出来的式子,正是涡度拟能方程。
它来自上头对流导数与w的标量积,是对于定域分布的涡度。
其中最右边的散度项通常积分为零,和脑子一样不太需要。
右边剩下的两项分别对应通过涡线拉长产生涡度拟能,以及因为粘滞力损耗的涡动拟能。
从这个式子可以直观看出涡动拟能就像力学能量一样,可以被摩擦力耗散掉。
这个公式在后世讨论湍流的时候会被反复提及,算是一个标识型的公式。
更重要的是
众所周知。
大气扩散属于湍流扩散,目前有三种广泛的应用理论:
梯度输送理论、
湍流统计理论、
相似理论。
而这个式子便是湍流统计理论的重要核心,后世在这个基础上诞生了一种叫做wRF的模型。
没错。
wRF。
这是后世气象数值模拟预报最常见的模型,很多民科在家也用这玩意儿来跑数值。
当然了。
气象领域的民科要远比物理和数学领域的民科高智很多,二者存在很明显的差异。
气象领域的民科与其说是‘民科’。
不如说更像是那些开车载着天文望远镜去看星星的天文爱好者,很少有太多出格的言论。
至少不会动不动就表示自己发明出了永动机,然后一看图纸特么的是太极图
气象领域的民科最喜欢的就是在家里默默跑当地的天气模型,然后巴望着天空看自己的结果准不准确,整体来看还是比较佛系的。
总而言之。
wRF即便是在2023年都属于非常重要的模型,遑论眼下这个时期了。
即便
此时出现的只是一个雏形。
随后叶笃正又把公式引申到了等压面和等密度面领域,进行起了环流的相关计算。
期间乔彩虹这姑娘也兴致冲冲的上前旁观了两分钟,等回到徐云身边的时候表情就变成了这样:
@v@
一个小时后。
叶笃正和陶诗言合力推导出了完整的涡度场,拟合出了一个特殊的数学模型。
从徐云的角度来看。
这个模型和后世的wRF依旧出入较大——毕竟这年头没有后世的算力,但核心逻辑还是类似的。
简单来说就是先采用了圆柱切线空间和水平映射,构建起局部空间并映射其邻居,构建起等轴映射。
https://xianzhe.cc/book/34463/15501527.html