【正文】 前雙側(cè)的主要方程（ 16）中的第二、三、四方程分別表達(dá)了質(zhì)量守恒、動量守恒、 和能量守恒。模擬快速壓縮機(jī)器里氣體運(yùn)動的主要目的是拓展到封閉混合氣體在高溫高壓下可能發(fā)生的自燃現(xiàn)象。 參考文獻(xiàn)： 1. and A. Thomas, An opposed piston rapid pression machine for preflame reaction . Inst. Mech. Engrs. 183 (1969) 365–385. 2. P. Park, Rapid pression machine measurements of ignition delays for primary reference fuels. . thesis, . (1990). 3. . Griffiths and . Barnard, Flame and Combustion. London: Blackie academic and professional (1995) 328pp. 4. L. Brett, J. MacNamara, P. Musch and . Simmie, Simulation of methane autoignition in a rapid pression machine with creviced pistons. Combust. Flame. 124 (2020) 326–329. 5. . Tabaczynski, . Hoult and . Keck, High Reynolds number flow in a moving corner. J. Fluid 42 (1970) 249–255. 6. D. Lee, Autoignition measurements and modeling in a rapid pression machine. . thesis, .(1997). 7. , Combustion Theory. Massachusetts: Perseus Books (1985) 680pp. 8. . Buckmaster(editor), The Mathematics of Combustion (2nd ed.). New York: SIAM (1985) 254pp. 9. I. Barin, Thermochemical Data of Pure Substances, Parts I and II. Weinheim: John Wiley and sons (1993) 1739pp. 10. . Hirschfelder, . Curtiss and . Bird, Molecular Theory of Gases and Liquids. New York: John Wiley (1963) 1249pp. 11. . John, Gas Dynamics (2nd ed.). Boston: Allyn amp。然而，該模型 反映壓縮室 內(nèi)的核心最后的壓力。并且外部條件 ??qx 符合方程（ 15），并在最終保持不 變。嚴(yán)格的數(shù)學(xué)處理這些問題具有重要的價值，但這不是本文分析的目的。當(dāng)波到達(dá)活塞處是，它又反應(yīng)了各自的距離，（ P， ρ， T）又可以計算。鑒于一個事實(shí)，即 θ = 0（ 103）一般情況下，波的速度通常是活塞的 30倍。 數(shù)值方法 方程（ 8）進(jìn)行了數(shù)值綜合的差分算法。然 而，考慮變化的活塞速度使得問題復(fù)雜化了，我們將把活塞的運(yùn)動速度簡單化，因?yàn)榭傮w模型的運(yùn)動一旦完成，活塞的速度基本穩(wěn)定；參見 節(jié)。 V是混合氣體的平均速度。在上面的關(guān)系式中，Ｔｉ和Ｐｉ是初始值，Ｔ和Ｐ是一段時間后的值， γ (s)是絕熱指數(shù)。封閉的壓縮氣體造成 氣體壓力，溫度和密度迅速增加。結(jié)果表明，使用實(shí)驗(yàn)數(shù)據(jù)該模型有有較好的計算和預(yù)測能力 。它來源于文獻(xiàn) (4)，初步壓力和初始溫度分別為 344開爾文。然而， 這里一維空間的研究包括兩方面： ，可以成功抑制活塞運(yùn)動時產(chǎn)生的熱邊界層（文獻(xiàn) 6），從而使結(jié)果更加接近真實(shí)值。只有核心溫度上升到一定的水平，化學(xué) 反應(yīng)可以產(chǎn)生重大的影響，但這段時間通常很短（通常是幾毫秒）。下文討論的計算結(jié)果，我們使用恒定值 λ和 γ ?；钊粋€行程為 80 時間段，40 時間段可以被認(rèn)為是合適的的。應(yīng)該強(qiáng)調(diào)的是，為使 θ接近真實(shí)值 ，波會略增加 (p, ρ , T ) 。算法 我們現(xiàn)在考慮的漸近行為（ 8）在限制 ε → 第一次返回至中心線。 )( 21?O 的方程 ).10)0)0(()(1)0(()010(0,1)111(0,0)10(0,00110xvpxTTxTxTvTxpxvvvvxTTttt???????????????????????????????????????????????????????????????????? 2)11( 中的 1v? 可以寫成： ,0011 0),(sdxttsv ? ???????? ?? 這里我們用了邊界條件 0,01 ?? ?? xv ，把（ 12）代入 4)14( 中，并且有 : )))0)(()0)0(()(1)0((0})0)(({00)(00),(),(sssssdxtPtdtPxtttPtsTxxTTxTxTtsTTT???????????????????????????????????? ,)(0,0,10,10 0????????????xTTxTtTtc 時 )(tTc 隨著條件而變化。 把 (17)和（ 19）帶入（ 16）的第四個方程中，且 ???????? ***0*0*0 ,0,1,0 zzTTv 得出： ),2/0()(0 2*0.*0**0*0.