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

bookmark

第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方程,以及非等溫管狀反應器的溫度與壓力場……

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

“總的來說。”

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

肯定是不成的。

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

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

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

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

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

“我知道了。”

(本章完)

第593章 讓全國產三代機上閱兵!第1249章 標準3,發射!第552章 殺人,還要誅心第439章 安70的命運第857章 大地之眼第1184章 東風,發射!第1276章 責無旁貸第1130章 裝配,出問題了?第1069章 軍用產品線?第1146章 一期工程,4000萬畝!第1275章 讓美國人出錢!(算昨天的)第920章 精確制導時代第614章 華夏戰鬥機的心臟第1112章 我給兩院帶來了一份禮物340.第340章 跑得飛快的西方記者第678章 知識產權還是要保護好第529章 獲獎懸念揭曉!第824章 反航母技術之一(新的一月,求月票390.第390章 三個月兩次專訪?第780章 戰略忽悠的意外效果第459章 我們中出了一個叛徒!第463章 要不要安70?第473章 半場開香檳的麻省理工第409章 飛機出口還是要嚴抓宣傳第1228章 地下堡壘的剋星第117章 來自的需求(爲上月月票金主flanke第9章 學霸之間,亦有差距!第413章 給父母買個房第57章 思路轉變第36章 作繭自縛第202章 回家第873章 打不過就加入第673章 有一句話叫做悶聲大發財第483章 有實力,纔有平等對話的資本!第79章 比幻影2000更強!第843章 我和我的大佬導師第702章 警衛九局第548章 佈局:二桃殺三士第1083章 我很大,但很輕第980章 做一家有夢想的公司271.第271章 獵殺獵人,殲轟電7(日萬)第1075章 走出國門,走向世界323.第323章 願者上鉤第1256章 鷹子,行不行啊?第693章 美國人的報復第1115章 涉核課題?第756章 科工委的進步325.第325章 圖334客機和霹靂12空空導彈第943章 總不能什麼都交給華夏人去幹吧?(4K)第1216章 磁流體發電?第805章 另一種方案的殲10第185章 00批次殲8C,交付部隊第461章 高性能航發的技術拼圖第76章 關注一下發動機(爲盟主flankerr加更第1279章 目標:突破臨近空間389.第389章 拿專利授權換個大獎第1012章 通信中繼平臺第431章 國際軍火市場的攪局者第930章 來自COMSOL的挑戰第1008章 這就開始搞綜合電力系統了?第149章 蘇霍伊的野望第1223章 相愛相殺的歐洲人322.第322章 俄國人也羨慕我們(爲盟主【死跑龍套的】加更55)第565章 真空電子束焊接技術402.第402章 試飛員的“護身符”260.第260章 殲轟7MLU?第591章 連升兩級的戰鬥英雄第133章 百噸級客機的可能性(爲flankerr的第1097章 準備返回第1265章 替NASA宣佈計劃失敗第202章 回家371.第371章 小試牛刀第990章 成年人的選擇是全都要第185章 00批次殲8C,交付部隊第1208章 高超音速競賽第456章 山鷹,兩個半月創造的奇蹟!339.第339章 凱旋358.第358章 國家重點實驗室首席科學家第188章 媽都認不出來的米格21改(二合一)第528章 頂級凡爾賽,江湖上盡是我的傳說第510章 論文又上封面了?第3章 現在的本科生都這麼厲害了?373.第373章 新軟件的第一個用戶!第183章 華夏不同意收購第119章 摸不透的需求373.第373章 新軟件的第一個用戶!第1245章 你有反衛星,我也有反衛星第799章 影響整個數學史的郵件291.第291章 平流層飛艇和國際互聯網第921章 大變局第1285章 波音的困擾第1050章 七艘航母的行蹤第100章 渦噴14的重生(下)(二合一)第100章 渦噴14的重生(下)(二合一)第1185章 有什麼東西飛過去了第57章 思路轉變第665章 現代戰爭是炮戰!第427章 打開歐洲市場的鑰匙第1056章 達索的小心思第1037章 我對錢沒有興趣