第96章 組會,震撼的各團隊(2/2)
他們————沒有看錯吧?
這些東西,都是這個高中生搞出來的?!
除了李建國和王海峰,剩下三個團隊的人都是滿臉的不可思議。
想想他們最開始還滿是對林葉的不相信,此時只感覺臉有些疼。
至於李建國和王海峰,此時也沉浸在一陣震撼當中。
這個年輕人,可真是————嚇了他們一跳。
這個思路,他們此前竟然從來沒有想到過。
但兩位教授的震撼,持續時間也很短暫。
不多時,作為計算專家的王海峰教授,就直接看出來了這個方案的致命問題。
他坐直了身體,雙手交叉放在了桌上。
「想法不錯,林同學。」王海峰開口了,語氣雖然客氣,但問題直指核心,「但是,你怎麼處理這兩個區域的交界面?顯式算法和隱式算法的通量計算方式完全不同,如果在界面處直接拼接,必然會產生虛假的數值反射和振盪。這對於高超聲速流動來說,任何一點虛假波動都會被激波放大,導致計算發散。」
「王教授說得對。」林葉早有準備,他自信地切換到下一頁PPT,「硬拼接肯定不行。所以,我準備在兩個區域之間,引入一個緩衝過渡層。」
「在這個過渡層里,不做那些非黑即白的切換,而是構造一個基於oid函數的權重因子α,然後在過渡層內求解一個混合通量,大概長這樣——」
【F—hybr id=aF—implicit+(1—a)F—eplicit】
「同時,」林葉指著公式中的一項修正項,「為了保證守恆性,我利用了梅林變換技巧,從數學上推導了一個通量修正項。理論上,這能確保在過渡層內的質量和能量守恆誤差控制在二階精度以內。」
梅林變換!
看著林葉給出的這一步,下面的張濤也露出微微感慨的神情。
他和林葉前幾天用到的方法,終究還是發揮了作用。
而王海峰看著那個精妙的權重函數構造,眼中的懷疑逐漸消退,取而代之的是一絲驚訝。
這個處理手法極其老練,完全不像是個新手。
「理論上————似乎是自洽的。」他沉吟著,腦海中開始推導後面的步驟,思考著具體的可行性。
「等等,我也有個問題。」一直沉默的李建國教授推了推眼鏡,他是做物理模型的,更關心物理真實性,「你這個平衡區的判定,如果太寬泛,會不會漏掉一些雖然緩慢但很關鍵的中間反應?比如一氧化氮的生成反應比較慢,如果你把它當成平衡態處理了,熱流預測可能會偏低。」
「這就是這個判據中漸近項的作用。」林葉解釋道,「李教授,我引入的判據不是固定的閾值,而是與流場的對流時間尺度耦合的。如果流速很快,即使反應慢,判據也會收緊,強制將其劃入剛性區計算。」
「從漸近分析的角度來看,只要判據的截斷誤差小於我們的離散誤差,物理精度就是可以保證的。」
林葉並沒有拿出具體的數據,因為他還沒跑過,但他用嚴密的數學邏輯和物理圖像,構建了一個看起來無懈可擊的理論框架。
會議室里陷入了短暫的沉默。
王海峰教授盯著那頁PPT看了足足半分鐘,然後轉過頭,看向周文淵,語氣中帶著一絲認真:「老周,這路子有點野,但邏輯上是通的。比那個算不動的特徵值分解要有操作性得多。」
李建國教授也微微頷首:「只要那個漸近判據能如他所說的那樣工作,物理上確實沒問題。這實際上是用空間換時間,很聰明的做法。」
周文淵看著兩位被說服的老搭檔,心中一塊大石落地,哈哈大笑了起來,他站起身,環視全場:「還有人有問題嗎?」
三個團隊的人都是你看看我,我看看你。
開玩笑,兩位教授都沒問題了,他們還能有什麼問題呢?
「那,既然各位都覺得沒問題了,那我們就試一試。海峰,這幾天可能要辛苦你們團隊了,把這個算法植入到求解器里去驗證一下。」
「沒問題!」王海峰大手一揮,「只要不再讓我算那什麼讓人頭大的特徵值就行了,寫幾行新代碼算什麼。今晚我就安排人手,咱們爭取三天內跑出個結果!」
>