“你最近,怎么老爱激动呢?”欧叶望向电视机,已切到下一条关于医保改革的新闻。
“北纬32度30分,西经117度19分。”沈奇在谷歌地图中找到这个坐标,位于美国西海岸城市圣迭戈附近的海域:“太玄乎了,真的是匪夷所思,在物理学上,我们称为复杂性,它将看似无关的事情联系起来,分子结构和位错运动、能量、动力、流体、纳维-斯托克斯方程以某种神秘的规律组合在一起,挺复杂的。”
“我俩通过严谨的数学计算得到一个理论上正确,却不存在实际意义的计算结果……这里要修正一下,这个计算结果有可能具备实际意义,但我需要修改数学物理模型,并找到具体的参数值,诸如经纬度、地震强度等等,我也不知道该找些什么参数,我得查询一下资料。”沈奇兴奋的搓搓手,地震学属于物理学的一个分支,他欠缺的是地质学方面的理论知识。
“你……难道要预测地震?”欧叶惊呆了,也有些害怕:“这,真的好吗?”
不怕和人争,就怕和天斗。
历史上大凡对抗老天爷的人物,要么被天收了,要么被其余人类抓起来。
欧叶担心大数学家沈奇深入物理研究后,走上一条疯狂的不归路。
“好了好了,本故事纯属虚构,预测地震什么的理论漏洞百出,洗洗睡了。”沈奇露出一副“我纯属YY你不要太认真”的表情。
“奇,你一定要,忘记这个危险的想法。”欧叶进入卧室,施施然解开了睡袍的带子,半掩着门露出脑袋,红着脸难得主动了一次:“快点冲凉哦,等你……”
他肯定是憋坏了,才会胡思乱想……欧叶褪去睡袍,将卧室灯光调暗,她钻进被窝,身子情不自禁的发烫……
……
睡在同一张床上的男女,做了不同的梦。
女人梦到了她的老爸,老爸特别喜欢她的男票,还将祖传的一块玉送给男票。
男人在半梦半醒的神游状态中,似乎得到了某种指引,他将地理等级升到了5级。
次日清晨沈奇醒来,去到办公室继续他的计算,为此他查询了一些基础性的地质学文献。
“除了数学,其他科目都是坑。”
沈奇为了学好物理,已将化学升到了7级,现在将地理升到5级。
地震学属于物理学旗下的品牌,是固体地球物理学的一个分支。
沈奇将地理升到5级的职业级,多了解一些地质学方面的常识,对于地震学乃至固体地球物理学的深入研究会有帮助。
“地球岩流圈内部的对流取决于地球内核裂变反应产生的能量,这个过程导致板块边缘压力的产生,由于边缘的摩擦,压力能量间歇性释放而非持续性,岩板块相互挤压而非滑动,这种突发的挤压造成了我们所说的地震。”
“我需要一些地震指标参数,以建立一种全新而有效的数学模型。”沈奇有目的性的开始查阅资料。
地震强度定义为=lgs,这个指标非常关键,ark。
在美国南加州地震网的官网上,可以查到1983年至今的地震时序。
对于尺寸为lgs>15的地震,符合一个极其简单明了的幂律衰减规律:N(S>s)∝s-B,s≥1015,B>095
对于菲奖得主沈奇而言,这个定律太过简单,这是著名的关于地震尺寸频率的古登堡-里克特定律,里氏XX级地震中的里氏,就是古登堡-里克特定律中的美国地震学家里克特。
地震强度lgs<2的地震通常被称为微震,人类难以感知,只能通过地震仪记录。
美国南加州地震网上的绝大部分数据采集来自微震,微震也会产生能量,能量是另一个重要指标,沈奇ark能量指标,以及美国南加州地震网所记录的,自1983年以来每一次地震发生的具体时间和坐标经纬度。
孕震过程的流体具有不可忽视的重要性,只不过流体运动在其孕震、发震、余震的过程中,难以采集精确数据,这是困扰地震界及物理界的一大难题。
沈奇忽然间清醒了,以微积分方式描述流体运动的N-S方程不仅适用于天空和海洋,它同样适用于地表之下!
理论上来说是这样的。
沈奇未加证实的理论。
一直到现在,沈奇几乎是依靠直觉在做一些有可能毫无意义的事情。
复杂性的表象需要严密的数学语言给出定量化的内在支撑。
沈奇重新梳理昨天的计算结果,全新而有效的数学模型需要一些核心工具作为灵魂。
这具灵魂应该是一个或一组数学公式、表达式、等等。
N-S方程尚未找到通用的求解公式,以得到通解。
特解倒是有一百多个,基于其中一个特解,沈奇花费数周时间,瞒着欧叶和所有人,一个人推导、计算、验证,最终得到一个复杂的等式:
1/2d/dt∫Ω∣As/2ω∣2dx+v∫Ω∣As+1/2ω∣2dx+α∫Ω∣As/2ω∣2dx+……(B(u,u),Asω)=<χ1,χ2>
根据这个核心等式,建立一个新的数学模型,沈奇编写了一个计算机程序。
系统中没有单独开通计算机科目,数学中的一些分支包含部分计算机知识与技能。
沈奇用基础的Fortran90语言编写程序,分为设定参数个数及其误差范围的主程序和利用模型计算<χ1,χ2>的子程序,来拟合从美国南加州地震网上下载的参数数据。
“激动人心的时刻,到了。”
沈奇深呼吸一口,重重敲击回车键,启动程序。
程序开始执行,沈奇自编自导的构想是否具备实际意义,答案即将揭晓。
程序瞬间生成大量的矩阵、常微分方程和偏微分方程!
沈奇倒吸一口冷气!
靠他自己手动计算这些矩阵和方程,不吃不喝不睡,至少需要连续计算十年以上。
加上欧叶,两人一起干,也得五年。
主机CPU风扇急速转动,CPU快要被玩坏了。
“坚持,再坚持一会儿啊!”
沈奇大吼到,他这台高性能电脑刚换不久,好几万美元的机器不可能这么弱鸡。
好几万美元的机器汗流浃背、泪流满面。
大量的原始数据洪流还在汹涌的倾泻进来,它们是美国南加州地震网自1983年以来采集到的所有真实地震数据,以及沈奇独家发明的基于N-S方程的流体虚拟数据。
呜呜呜!
风扇狂转,它已超越五档,直奔八档而去。
CPU在哭泣、哀嚎,它正处于一生中最难熬的时刻。
艰难时光持续了半个小时,CPU终于解脱了。
机箱猛然抖动了一下,遂恢复平静。
微观、宏观、分子、地壳,流体、能量、间歇性的摩擦,看似无关的事情却又遵循某种神秘规律,沈奇自编自导的这出数理戏剧落下帷幕。
程序启动之前,沈奇甚至怀疑这是荒诞剧。
程序结束之后,沈奇鼓起了掌,这是写实题材,真实感人。
大数据计算最终得出一个极其简单的结果:<55,167>
沈奇早已购置了地球仪,就摆在他的办公桌上。
“纬度55,经度167,目前我只能得到这种程度的计算结果,毕竟我用的是N-S方程的特解。”
很快的,沈奇找到四组坐标:
北纬55度,西经167度,这个坐标位于白令海域,属于美国的领海,在阿留申群岛中的乌尼马克岛西北约200公里处的北太平洋上。
北纬55度,东经167度,这个坐标同样位于白令海域,属于俄罗斯的领海区域,距尼科尔斯科耶岛仅有50公里。
南纬55度,西经167度,这个坐标位于广阔的南太平洋上,离这个坐标最近的大洲是南极洲,最近的国家是新西兰,此处距新西兰约2000公里。
南纬55度,东经167度,这个坐标同样位于南太平洋上,它离新西兰更近了,距新西兰最南端的坎贝尔岛只有100公里。