306.第306章 非線性偏微分方程組降維問題

雖然常浩南的獲獎讓學校多了不少宣傳口的工作要做,但唐林天自然也不會忘了前幾天答應的事情。

沒花幾天功夫,學校裡接的外網就已經可以訪問那幾個比較主流的學術數據庫了。

常浩南從大概半年前,剛開始認真考慮編寫一個全新的仿真建模軟件時就非常清楚地意識到,多物理場,尤其是強耦合多物理場問題的研究,本質上是對非線性偏微分方程組的求解。

但這種事情,落實到工程領域的操作上,往往就不是一句話那麼簡單了。

尤其是考慮到現如今的超級計算機運算速度並不樂觀的情況下。

絕大多數偏微分方程都沒有辦法求得解析解,至少短時間內,只能從數值解的方向下功夫。

很多在數學上十分具有美感的解法未必實用。

傳統上對於非線性偏微分方程動態系統的降維主要採用基於變量離散的方法,典型的比如有限元法,有限體積法和有限差分法,堪稱這一領域的御三家。

但也不是沒有其它的思路。

就比如常浩南某天晚上在機房休息時,無意中看到的這篇論文。

儘管是應用數學領域的文章,但卻發在了一份看上去毫無關係的化學工程領域期刊上。

Chemical Engineering Journal

一份在十幾年後算是聲名赫赫,但這功夫只是剛剛創刊,並不起眼的雜誌。

之所以會吸引他點進去,並用每秒幾Kb的速度下載下來,主要是因爲摘要寫得太有吸引力了。

“目前通用的有限差分法和有限元方法對非線性偏微分方程動態系統進行降維只能得到維數很高的常微分方程系統,在近四十年的時間裡,基於變量分離的系統降維方法得到了飛速發展,在滿足一定的條件下能避免基於空間離散方法帶來的一些本質問題,將一類非線性偏微分方程動態系統降至較低的維數,便於快速分析計算、優化及主動控制器的實現,可以應用於對化學工程領域內常見的力熱耦合問題進行數值分析……”

儘管涉及到的具體問題和飛行器設計風馬牛不相及,但裡面提到的力熱耦合本來也是常浩南目前面對的最基礎,也是最緊迫的問題。

這段摘要簡直說到了他的心坎裡。

他相信幾個月之前,當盧育英在蓉城第一次看到自己那篇論文的時候,內心的通透感也就不外乎如此。

幾分鐘的下載時間從未如同現在這般漫長。

常浩南緊盯着屏幕上面一格一格的進度條,幾乎在下載完成的一瞬間就點開了那份文檔。

“衆所周知,任何一個連續函數能被傅里葉級數序列的展開式近似表示,基於上述原理,非線性偏微分方程中的時空親和變量,能夠展開成一個無限維空間基函數集合和其對應的時間係數的級數和的形式:

X(z,t)=(i=1,∞)∑φi(z)xi(t)

其中xi(t)表示每個基函數φi(z)對應的時間係數……”

確實很基礎。

時空變量分離技術並不是什麼新鮮玩意,任何一本數學物理方法或者類似的教材上都能找到,只是一般認爲適合使用分離變量法的偏微分方程應該具有一定的形式和特徵,如線性、齊次、可分離、係數只依賴於一個變量等等,這極大地限制了此類方法的應用。

因此常浩南迅速略過了這部分內容,直接看向了第三節,往往也是正文的第一節:

爲了詳細和清楚地闡述非線性偏微分方程動態系統降維的方法,本小節釆用拋物型非線性偏微分方程系統作爲對象進行闡述……

“來了!”

看到感興趣內容的他精神一振,就連剛剛的些許睏意都瞬間煙消雲散。

邊界條件和初始條件分別爲:

其中x(z,t)表示時空狀態變量,且爲定義在空間區域[a,b]上的無窮維希爾伯特空間上的連續函數。表示空間座標,z∈[a,b]表示空間座標,爲過程定義的實數域上的子空間,t∈[0,∞)表示時間變量……

……

最終,可以得到希爾伯特空間H([a,b])中上述非線性偏微分方程系統的表達形式:

x(z,t)/t=Ax(z,t)+Bu(z,t)+(x,z,t)

x(z,0)=x0(z)

下面給出兩個仿真實例,分別是一維空間的無量綱Kuramoto-Sivashinsky方程,以及非等溫管狀反應器的溫度與壓力場……

“嗯……有點東西……”常浩南看到後面,內心瞭然地點了點頭。

“總的來說。”

他從旁邊的打印機裡面抽出一張紙,開始自言自語地總結起來,

“首先,選擇合適的空間正交基函數且採用時空分離技術對非線性偏微分方程動態系統進行時空變量分離,即將系統的時空親合變量在選定或求得的正交空間基函數上展開,將展開式代入原系統後結合非線性伽遼金方法……”

一個小時的時間很快在他的寫寫畫畫中過去了。

雖然文章中用於闡述理論的對象只是個非常簡單的拋物型系統,但後面舉出來的兩個應用算例確實還算可以,配得上作者在摘要裡面吹出來的牛逼。

這篇文章甚至值得投一個更高影響力的期刊,之所以出現在這裡,大概率是因爲作者和主編出自同一個學校,收到了約稿的邀請。

實際上,常浩南總結到最後,還發現了作者本人都沒有寫出來的部分。

文章裡面的這套方法不僅可以應用於傳熱和流場計算,只要稍加修改,甚至可以用於處理傳質問題和化學反應過程本身。

換句話說,化工生產過程涉及到的全部特徵“三傳一反”,都可以被囊括進去。

當然,沒寫出來未必是作者沒發現,很可能是留着東西準備再發一篇文章出來……

“不過麼……問題也還是有的。”

常浩南看着面前已經被寫滿的三張草稿紙。

雖然可以應用的領域非常廣,但並不意味着這篇文章裡提到的方法就是什麼萬能鑰匙,可以直接搬到常浩南所需要的場景下面來。

“釆用特徵函數爲空間基函數結合權重殘差方法對非線性偏微分方程動態系統進行降維可以得到有限維常微分方程動態系統來近似原系統的無窮維動態,但本質上還是採用線性手段近似,對於真正的強非線性問題而言仍然不夠,可是如果在降維過程中採用別的空間基函數,如傅里葉序列函數和正交基函數等,又可能與非線性偏微分方程動態系統本身的特徵毫無關聯……”

想到這裡,他側過頭看了一眼旁邊擺着超算的另外一個房間。

在理論上當然沒什麼問題,不過真要是開始計算的話……

由於這個自己負責的計算中心還只是剛剛啓用,因此目前用到它的項目不多,儘管如此,也已經讓這臺超算的負荷拉到一個不低的程度了。

如果擱在十年之後,按照文章裡面的思路硬算未必不可行,但以如今國內的超算水平,恐怕算個相控陣雷達陣面的力熱電耦合高低得花上幾年時間……

有這時間都造個測試版本出來了。

肯定是不成的。

“如果採用平衡截斷方法或者最優化方法呢……”

常浩南手中的筆尖重新開始在紙上滑動起來。

很快,第四頁和第五頁草稿紙也被寫滿。

機房中傳來吱嘎吱嘎的設備運行聲。

窗外的月亮從地平線爬到半空,又逐漸落下,最終迎來一輪朝陽。

“我知道了。”

第47章 主動顫振控制,成功!第467章 日本軸承產業感到有些不安第172章 大推力渦扇發動機的暢想第650章 之前換裝殲8C的部隊都後悔了第710章 巴爾幹,終局第480章 國產化渦扇9!第741章 數學不會就是不會第665章 現代戰爭是炮戰!第131章 展示一下殲83的爬升率364.第364章 超越時代的仿真建模工具!第127章 我怎麼突然多了篇論文出來?(二合第718章 新裝備太多,也是一種煩惱第88章 進度飛快的新舟60第202章 回家第54章 組織已經決定了227.第227章 大飛機產業園!協議達成!第182章 殲轟7的未來第189章 DSI進氣道的設想第417章 正式開工!第656章 年終總結,與未來規劃287.第287章 該省的省,該花的花,轟6再改改268.第268章 空軍的財神爺第40章 好好聽講!第664章 組織已經決定了,準備主動出擊!第575章 給小日子挖個坑第427章 打開歐洲市場的鑰匙275.第275章 航空工業今年賺麻了第160章 阿古斯塔要深度合作!第35章 保研名額第755章 分頭行動381.第381章 設計一個DSI進氣道練練手229.第229章 給年輕同志壓擔子第526章 佈局工業控制領域第503章 常浩南同志成了香餑餑第31章 全新的工程模型第774章 數學界的震動348.第348章 江城重型機牀廠第130章 航展準備(爲月票金主flankerr加更第654章 組織上已經給你準備好了326.第326章 客機,C808和C909第109章 抉擇第635章 機彈彈彈組網協同無源定位第408章 軍宣拯救計劃第721章 南亞風雲第428章 逐漸成型的航空產業第99章 渦噴14的重生(上)(二合一)第49章 殲7F306.第306章 非線性偏微分方程組降維問題第176章 喝水帶來的靈感(第600章 軸承,壞了第680章 準備搞事,來自東方的支援第534章 半導體產業鏈,得從頭抓起333.第333章 人民海軍,向前進第589章 航校的老熟人第203章 在做了在做了(新建文件夾)第555章 法國艦隊來訪第604章 躬匠精神就完事了?第697章 雷達抓鬼,高炮驅邪252.第252章 搞預警機也要有備無患第550章 回擊挑釁!第5章 殲7的改造方案第78章 殲83,第二次首飛!第461章 高性能航發的技術拼圖第601章 NSK集團,東窗事發第204章 挑戰米格25!378.第378章 通往高性能航發之路!385.第385章 “息壤”戰洪魔第779章 先把華夏航發的名聲給打出去第86章 遠程偵察第525章 Nh90直升機,看向東方381.第381章 設計一個DSI進氣道練練手第416章 抉擇,渦扇的10技術路線233.第233章 俘獲電子戰吊艙了?第605章 《一個軸承引發的悲劇》第658章 把航母項目整個搬回來?第481章 無所適從的英國人第199章 高級軟臥第551章 裝備工業司的開門紅第147章 俄國人有點懵了第20章 中華兒女多奇志第187章 殲7,超進化!(二合一)第576章 世界,正悄無聲息地發生變化第647章 側衛進化論第176章 喝水帶來的靈感(第732章 法國人被卡脖子了?277.第277章 海空對抗演習第40章 好好聽講!第478章 怎麼,你又要改國軍標?第436章 有人要買我們的客機?288.第288章 我一個本科生,怎麼就成了超算第134章 衆生相(二合一)第674章 夢想集團,不太對勁第786章 你這不是撞槍口上了麼?第764章 去京城,看閱兵!第445章 竟然真有用?第591章 連升兩級的戰鬥英雄第547章 合資,但不能當買辦第644章 起飛吧,太行!第671章 專職秘書第765章 十月一日,晴
第47章 主動顫振控制,成功!第467章 日本軸承產業感到有些不安第172章 大推力渦扇發動機的暢想第650章 之前換裝殲8C的部隊都後悔了第710章 巴爾幹,終局第480章 國產化渦扇9!第741章 數學不會就是不會第665章 現代戰爭是炮戰!第131章 展示一下殲83的爬升率364.第364章 超越時代的仿真建模工具!第127章 我怎麼突然多了篇論文出來?(二合第718章 新裝備太多,也是一種煩惱第88章 進度飛快的新舟60第202章 回家第54章 組織已經決定了227.第227章 大飛機產業園!協議達成!第182章 殲轟7的未來第189章 DSI進氣道的設想第417章 正式開工!第656章 年終總結,與未來規劃287.第287章 該省的省,該花的花,轟6再改改268.第268章 空軍的財神爺第40章 好好聽講!第664章 組織已經決定了,準備主動出擊!第575章 給小日子挖個坑第427章 打開歐洲市場的鑰匙275.第275章 航空工業今年賺麻了第160章 阿古斯塔要深度合作!第35章 保研名額第755章 分頭行動381.第381章 設計一個DSI進氣道練練手229.第229章 給年輕同志壓擔子第526章 佈局工業控制領域第503章 常浩南同志成了香餑餑第31章 全新的工程模型第774章 數學界的震動348.第348章 江城重型機牀廠第130章 航展準備(爲月票金主flankerr加更第654章 組織上已經給你準備好了326.第326章 客機,C808和C909第109章 抉擇第635章 機彈彈彈組網協同無源定位第408章 軍宣拯救計劃第721章 南亞風雲第428章 逐漸成型的航空產業第99章 渦噴14的重生(上)(二合一)第49章 殲7F306.第306章 非線性偏微分方程組降維問題第176章 喝水帶來的靈感(第600章 軸承,壞了第680章 準備搞事,來自東方的支援第534章 半導體產業鏈,得從頭抓起333.第333章 人民海軍,向前進第589章 航校的老熟人第203章 在做了在做了(新建文件夾)第555章 法國艦隊來訪第604章 躬匠精神就完事了?第697章 雷達抓鬼,高炮驅邪252.第252章 搞預警機也要有備無患第550章 回擊挑釁!第5章 殲7的改造方案第78章 殲83,第二次首飛!第461章 高性能航發的技術拼圖第601章 NSK集團,東窗事發第204章 挑戰米格25!378.第378章 通往高性能航發之路!385.第385章 “息壤”戰洪魔第779章 先把華夏航發的名聲給打出去第86章 遠程偵察第525章 Nh90直升機,看向東方381.第381章 設計一個DSI進氣道練練手第416章 抉擇,渦扇的10技術路線233.第233章 俘獲電子戰吊艙了?第605章 《一個軸承引發的悲劇》第658章 把航母項目整個搬回來?第481章 無所適從的英國人第199章 高級軟臥第551章 裝備工業司的開門紅第147章 俄國人有點懵了第20章 中華兒女多奇志第187章 殲7,超進化!(二合一)第576章 世界,正悄無聲息地發生變化第647章 側衛進化論第176章 喝水帶來的靈感(第732章 法國人被卡脖子了?277.第277章 海空對抗演習第40章 好好聽講!第478章 怎麼,你又要改國軍標?第436章 有人要買我們的客機?288.第288章 我一個本科生,怎麼就成了超算第134章 衆生相(二合一)第674章 夢想集團,不太對勁第786章 你這不是撞槍口上了麼?第764章 去京城,看閱兵!第445章 竟然真有用?第591章 連升兩級的戰鬥英雄第547章 合資,但不能當買辦第644章 起飛吧,太行!第671章 專職秘書第765章 十月一日,晴