超音速流動與燃燒的大渦模擬基礎

1、超燃沖壓發(fā)動機燃燒室超燃沖壓發(fā)動機燃燒室流動與燃燒的大渦模擬方法流動與燃燒的大渦模擬方法探討與展望探討與展望 n匯報人：張宏達匯報人：張宏達n院系：熱科學和能源工程系院系：熱科學和能源工程系總結與展望總結與展望第四部分第四部分超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分超燃沖壓發(fā)動機研究概述超燃沖壓發(fā)動機研究概述第一部分第一部分匯報提綱匯報提綱超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分研究背景及意義研究背景及意義超燃沖壓發(fā)動機燃燒室概述超燃沖壓發(fā)動機燃燒室概述湍流數值模擬的發(fā)展現(xiàn)狀湍流數值

2、模擬的發(fā)展現(xiàn)狀在在5050年前，人們意識到高超聲速技術在軍事年前，人們意識到高超聲速技術在軍事和空間技術上的應用價值，便廣泛開展了關于高和空間技術上的應用價值，便廣泛開展了關于高超聲速技術的研究，隨著航天技術的飛速發(fā)展，超聲速技術的研究，隨著航天技術的飛速發(fā)展，世界各國在高超聲速技術上的競爭更加激烈。高世界各國在高超聲速技術上的競爭更加激烈。高超聲速技術主要指研制高超聲速（超聲速技術主要指研制高超聲速（Ma5Ma5）飛行器）飛行器所需要的相關技術。吸氣式高超聲速飛行器無需所需要的相關技術。吸氣式高超聲速飛行器無需自帶氧化劑，有效載荷大，飛行成本低，可控性自帶氧化劑，有效載荷大，飛行成本低，可控

3、性強，安全性高并可以長時間重復使用，將對未來強，安全性高并可以長時間重復使用，將對未來的天地往返運輸系統(tǒng)和軍事攻防對抗體系提供了的天地往返運輸系統(tǒng)和軍事攻防對抗體系提供了非常有價值的新能力。因此目前人們更為關注的非常有價值的新能力。因此目前人們更為關注的是是吸氣式高超聲速飛行器吸氣式高超聲速飛行器相關技術。相關技術。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分在吸氣式推進系統(tǒng)中，最關鍵的技術就在吸氣式推進系統(tǒng)中，最關鍵的技術就是是吸氣式沖壓發(fā)動機吸氣式沖壓發(fā)動機。吸氣式沖壓發(fā)動機。吸氣式沖壓發(fā)動機一般由進氣道、隔離段、燃燒室和尾噴管一般由進氣道、隔離段、燃燒室和尾噴

4、管組成。在沖壓發(fā)動機各部件中，組成。在沖壓發(fā)動機各部件中，燃燒室燃燒室是是能量轉換的主要部件，燃燒室性能對發(fā)動能量轉換的主要部件，燃燒室性能對發(fā)動機研制的成敗起先決作用。機研制的成敗起先決作用。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分由于氣流在發(fā)動機燃燒室內的駐留時間較短，由于氣流在發(fā)動機燃燒室內的駐留時間較短，因此如何因此如何使使燃料的化學能在有限的時間和空間內燃料的化學能在有限的時間和空間內最大限度的轉化為熱能是最大限度的轉化為熱能是超燃沖壓發(fā)動機燃燒室超燃沖壓發(fā)動機燃燒室研究中的一個關鍵技術研究中的一個關鍵技術。通常情況下，采用支板。通常情況下，采用支板構

5、型、后掠斜坡及旋流噴嘴等措施，在燃燒室的構型、后掠斜坡及旋流噴嘴等措施，在燃燒室的局部區(qū)域內產生旋渦，形成回流區(qū)，增加了燃料局部區(qū)域內產生旋渦，形成回流區(qū)，增加了燃料的摻混及混合氣體的停留時間和反應時間，提高的摻混及混合氣體的停留時間和反應時間，提高燃燒效率，并縮短了燃燒區(qū)的長度。燃燒效率，并縮短了燃燒區(qū)的長度。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分研究背景及意義研究背景及意義超燃沖壓發(fā)動機燃燒室概述超燃沖壓發(fā)動機燃燒室概述湍流數值模擬的發(fā)展現(xiàn)狀湍流數值模擬的發(fā)展現(xiàn)狀超燃沖壓發(fā)動機燃燒室內

6、氣體仍為超聲速，超燃沖壓發(fā)動機燃燒室內氣體仍為超聲速，燃燒室內能否實現(xiàn)成功點火與火焰穩(wěn)定是影響整燃燒室內能否實現(xiàn)成功點火與火焰穩(wěn)定是影響整個發(fā)動機工作的關鍵因素。個發(fā)動機工作的關鍵因素。一般來說，火焰一般來說，火焰穩(wěn)定穩(wěn)定可以通過增大火焰?zhèn)鞑ニ俣群蜏p少氣流速度兩種可以通過增大火焰?zhèn)鞑ニ俣群蜏p少氣流速度兩種方法實現(xiàn)。但考慮實際可操作性，一般采用后者。方法實現(xiàn)。但考慮實際可操作性，一般采用后者。通過障礙物在流場中形成回流區(qū)是常用的火焰穩(wěn)通過障礙物在流場中形成回流區(qū)是常用的火焰穩(wěn)定方法。超燃沖壓發(fā)動機中常用的火焰穩(wěn)定器包定方法。超燃沖壓發(fā)動機中常用的火焰穩(wěn)定器包括：壁面燃料橫向噴射、支板、后向臺階和

7、凹腔括：壁面燃料橫向噴射、支板、后向臺階和凹腔火焰穩(wěn)定器等?；鹧娣€(wěn)定器等。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分壁面橫向噴射壁面橫向噴射：橫向射橫向射流會使噴孔附件邊界層分離，流會使噴孔附件邊界層分離，在噴孔上下游均形成回流區(qū)，在噴孔上下游均形成回流區(qū)，部分燃料被卷入回流區(qū)，與部分燃料被卷入回流區(qū)，與回流區(qū)內空氣混合、燃燒，回流區(qū)內空氣混合、燃燒，達到穩(wěn)定火焰作用；另一方達到穩(wěn)定火焰作用；另一方面，由于燃燒室內氣流是超面，由于燃燒室內氣流是超音速，噴流會在主流區(qū)形成音速，噴流會在主流區(qū)形成弓形激波，也具有一定的穩(wěn)弓形激波，也具有一定的穩(wěn)定火焰能力。但是這種穩(wěn)定

8、定火焰能力。但是這種穩(wěn)定方式穩(wěn)焰范圍小，且會在燃方式穩(wěn)焰范圍小，且會在燃燒室內造成過高的總壓損失，燒室內造成過高的總壓損失，一般無法單獨使用。一般無法單獨使用。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分支板支板：具有燃料噴注和火焰穩(wěn)定的作用，還具有燃料噴注和火焰穩(wěn)定的作用，還能起到氣流壓縮作用，主流區(qū)域支板后部形成典能起到氣流壓縮作用，主流區(qū)域支板后部形成典型的低速回流區(qū)，燃料與空氣在回流區(qū)內進行混型的低速回流區(qū)，燃料與空氣在回流區(qū)內進行混合、燃燒，實現(xiàn)穩(wěn)定火焰的目的；另一方面，支合、燃燒，實現(xiàn)穩(wěn)定火焰的目的；另一方面，支板后方形成的高湍流度尾跡有助于提高火焰的傳

9、板后方形成的高湍流度尾跡有助于提高火焰的傳播速度。支板的最大不足是直接置于高速氣流中，播速度。支板的最大不足是直接置于高速氣流中，燃燒室內氣流總溫較高，支板的熱防護問題是影燃燒室內氣流總溫較高，支板的熱防護問題是影響支板應用的主要因素。響支板應用的主要因素。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分后向臺階后向臺階：使用：使用最多的一種火焰穩(wěn)定器，它最多的一種火焰穩(wěn)定器，它形成的回流區(qū)較大，因而火焰穩(wěn)定范圍也比較寬形成的回流區(qū)較大，因而火焰穩(wěn)定范圍也比較寬。凹腔火焰穩(wěn)定器凹腔火焰穩(wěn)定器：新型超聲速火焰穩(wěn)定技術新型超聲速火焰穩(wěn)定技術，一般根據長深比不同，將凹腔火焰穩(wěn)

10、定器分為，一般根據長深比不同，將凹腔火焰穩(wěn)定器分為閉式凹腔（閉式凹腔（L/D10L/D10）和開式凹腔（）和開式凹腔（L/D10L/D10），閉），閉式凹腔剪切層偏向于凹腔內部，極限條件下與凹式凹腔剪切層偏向于凹腔內部，極限條件下與凹腔底壁再附，內陷的剪切層將凹腔內流體分為前腔底壁再附，內陷的剪切層將凹腔內流體分為前后兩部分，并形成兩個回流區(qū)，開式凹腔剛好相后兩部分，并形成兩個回流區(qū)，開式凹腔剛好相反，形成一個大的回流區(qū)。反，形成一個大的回流區(qū)。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分研究背

11、景及意義研究背景及意義超燃沖壓發(fā)動機燃燒室概述超燃沖壓發(fā)動機燃燒室概述湍流數值模擬的發(fā)展現(xiàn)狀湍流數值模擬的發(fā)展現(xiàn)狀數值計算是湍流研究和預測的主要手段，目數值計算是湍流研究和預測的主要手段，目前常用的湍流計算方法可以分為兩類，一類是基前常用的湍流計算方法可以分為兩類，一類是基于于reynoldsreynolds平均的平均的N-SN-S方程的湍流模式理論方程的湍流模式理論（RANSRANS），另一類是湍流的高級數值模擬，包括），另一類是湍流的高級數值模擬，包括大渦模擬（大渦模擬（LESLES）和直接數值模擬（）和直接數值模擬（DNSDNS）。）。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究

12、概述第一部分第一部分RANSRANS：不管哪類具體的雷諾平均方法，平均：不管哪類具體的雷諾平均方法，平均的結果都是將脈動運動時空變化的細節(jié)一概抹平，的結果都是將脈動運動時空變化的細節(jié)一概抹平，丟失了包含在脈動運動中的大量有重要意義的信丟失了包含在脈動運動中的大量有重要意義的信息。息。DNSDNS：直接求解：直接求解N-SN-S方程，理論上這種方法是方程，理論上這種方法是最精確的，但是最大的困難在于：湍流脈動運動最精確的，但是最大的困難在于：湍流脈動運動中包含著大小不同尺度的渦，其最大尺度中包含著大小不同尺度的渦，其最大尺度L L可與平可與平均運動的特征長度相比，而最小尺度取決于粘性均運動的特征

13、長度相比，而最小尺度取決于粘性耗散速度，在湍流統(tǒng)計理論中證明了大小尺度比耗散速度，在湍流統(tǒng)計理論中證明了大小尺度比值與雷諾數的關系：值與雷諾數的關系：超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分3 4LLR為了模擬湍流運動，計算網格的尺度應大到為了模擬湍流運動，計算網格的尺度應大到足以包含最大尺度的渦，應小到足以分辨最小渦足以包含最大尺度的渦，應小到足以分辨最小渦的運動，假設為二維網格，整個計算區(qū)域的網格的運動，假設為二維網格，整個計算區(qū)域的網格點總數應不少于：點總數應不少于：計算要模擬的時間長度應大于大渦的時間尺計算要模擬的時間長度應大于大渦的時間尺度度 ，而計算

14、的時間步長應小于小渦的時間尺度，而計算的時間步長應小于小渦的時間尺度 ，因此所需要的時間步數應不小于：因此所需要的時間步數應不小于：超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分6 4LNR3 4LLR因此，總的計算量正比于因此，總的計算量正比于 ，假定計算，假定計算中每一時間步的每一節(jié)點需要執(zhí)行中每一時間步的每一節(jié)點需要執(zhí)行100100條指令，則條指令，則對于對于 的問題就需要執(zhí)行約的問題就需要執(zhí)行約 條指令條指令，這意味著在運算速度為每秒千萬次量級的計算，這意味著在運算速度為每秒千萬次量級的計算機需要執(zhí)行約機需要執(zhí)行約1212天。天。目前國內外采用目前國內外采用D

15、NSDNS計算的算例基本都是槽計算的算例基本都是槽道流這種低速、低雷諾數、流場較為簡單的情況道流這種低速、低雷諾數、流場較為簡單的情況，對于高速、高雷諾數的模擬只有美國少數幾個，對于高速、高雷諾數的模擬只有美國少數幾個計算中心才有可能。計算中心才有可能。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分9 4LR510LR 13.2510LES：放棄對全部尺度范圍上的渦運動都進行：放棄對全部尺度范圍上的渦運動都進行數值模擬，改為只將比網格尺度大的大渦通過數值模擬，改為只將比網格尺度大的大渦通過N-S方程直接計算，而比網格尺度小的小渦，根據小方程直接計算，而比網格尺度小的小

16、渦，根據小渦對大渦運動的影響建立模型來模擬（即亞格子渦對大渦運動的影響建立模型來模擬（即亞格子模型），這就是大渦模擬的基本思想。模型），這就是大渦模擬的基本思想。超燃沖壓發(fā)動機燃燒室研究概述超燃沖壓發(fā)動機燃燒室研究概述第一部分第一部分超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分超音速流動大渦模擬超音速流動大渦模擬小結小結冷態(tài)流動研究可以將流動過程從燃燒過程中冷態(tài)流動研究可以將流動過程從燃燒過程中解耦出來，能夠對物理過程建立初步而必要的認解耦出來，能夠對物理過程建立初步而必要的認識，是發(fā)動機研究中必不可少的手段之一。眾所識，是發(fā)動機研究中必不可少的手段之一。眾所周知，周知，湍流研

17、究起源于不可壓流體，大渦模擬引湍流研究起源于不可壓流體，大渦模擬引入到超音速中，問題多多。入到超音速中，問題多多。超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分 精彩在下面！精彩在下面！聲明：本匯報所有公式中均忽略體積力項！聲明：本匯報所有公式中均忽略體積力項！超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分連續(xù)性方程連續(xù)性方程0iiutx0iiutxLES中中連續(xù)性方程連續(xù)性方程FavreFavre過濾過濾LESLES的連續(xù)性方程的連續(xù)性方程沒有出現(xiàn)未封閉項！沒有出現(xiàn)未封閉項！超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分動量方程動量方程123ijii

18、jkkijjiju uupSStxxx LES中中動量方程動量方程123ijiijkkijijijjijjuuupSSuuuutxxxx FavreFavre過濾過濾動量動量方程出現(xiàn)不封閉項方程出現(xiàn)不封閉項，亞格子應力需要?；瘉喐褡討π枰；?。sgsijijiju uuu常見的亞格子模型包括常見的亞格子模型包括SmagorinskySmagorinsky渦粘模渦粘模型、尺度相似模型、混合模型和一方程模型等。型、尺度相似模型、混合模型和一方程模型等。目前常用的模型是動態(tài)目前常用的模型是動態(tài)SmagorinskySmagorinsky渦粘模型，渦粘模型，個人認為，在超音速流中，一方程模型更合理。

19、個人認為，在超音速流中，一方程模型更合理。因為在超音速流中，湍流度增強，且具有很強的因為在超音速流中，湍流度增強，且具有很強的各向非均勻性，湍流脈動的局部平衡關系被打破，各向非均勻性，湍流脈動的局部平衡關系被打破，動態(tài)動態(tài)SmagorinskySmagorinsky渦粘模型不再適用渦粘模型不再適用。超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分能量方程能量方程超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分ijiiiiijiiuuEpuETRtxxxxxijiiiiijiiuEupuETRtxxxxx LES中中能量方程能量方程ijiiiiiiijiiiiiijiij

20、iijuEupuETREuEutxxxxxxpupuuuxx？能量方程出現(xiàn)更多的不封閉項，情況更加復雜能量方程出現(xiàn)更多的不封閉項，情況更加復雜。超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分sgsijijiijiuu亞格子粘性項，通常很小，予以忽略。亞格子粘性項，通常很小，予以忽略。 sgsiiiiiHEuEupupu亞格子焓通量，通常被?；瘉喐褡屿释?，通常被?；篜rsgssgsisgsiHHx 超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分超音速流動大渦模擬超音速流動大渦模擬小結小結個人觀點個人觀點：1 1、超音速流中大渦模擬的不封閉項遠多于不、超音速流中大渦

21、模擬的不封閉項遠多于不可壓流體中，因此需要模化，既然要模化，就要可壓流體中，因此需要?；?，既然要?；?，就要引入各種假設，模型的準確性和適用性有待考證引入各種假設，模型的準確性和適用性有待考證，發(fā)展精確且通用的模型難度很大！，發(fā)展精確且通用的模型難度很大！2 2、不封閉項中亞格子應力項使用一方程模型、不封閉項中亞格子應力項使用一方程模型更合理，亞格子焓通量及亞格子粘性項的模型我更合理，亞格子焓通量及亞格子粘性項的模型我也不也不太太清楚，這方面文章較少，值得深入研究。清楚，這方面文章較少，值得深入研究。超音速湍流流動大渦模擬超音速湍流流動大渦模擬第二部分第二部分超音速湍流燃燒大渦模擬超音速湍流燃燒

22、大渦模擬小結小結超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分聲明：本匯報所有公式中均忽略體積力項！聲明：本匯報所有公式中均忽略體積力項！LES中中連續(xù)性方程連續(xù)性方程LES中中動量方程動量方程燃燒場中與冷態(tài)場相同！燃燒場中與冷態(tài)場相同！超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分能量方程能量方程,111ijititiiiiijNNNnT mikkk im in ikmnimmnueu eqputxxxxx DTqh Y VRRTVVxM D 為

23、熱流量，等式右端第三項是熱輻射，為熱流量，等式右端第三項是熱輻射，第四項是第四項是DufourDufour效應，即濃度梯度引起的效應，即濃度梯度引起的熱擴散，這兩項通?？梢院雎?。熱擴散，這兩項通?？梢院雎?。等式右端等式右端第二項的擴散速度有如下計算公式：第二項的擴散速度有如下計算公式：iqiqiq超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分 ,111NNNpkpkT pT kpkpppp kpkkkkpkpkkpX XX XDDPTXV VYXYY ffDPPDYYT等式右端第二項是壓力引起的擴散，第三項是體等式右端第二項是壓力引起的擴散，第三項是體積力引起的擴散，第四項是積力

24、引起的擴散，第四項是SoretSoret效應，即溫度梯效應，即溫度梯度引起的質量擴散，這三項通?？梢院雎?。則得度引起的質量擴散，這三項通?？梢院雎?。則得到簡化的組分擴散公式：到簡化的組分擴散公式：,kkk ikiDYVYx 超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分能量方程能量方程綜上，整理得到能量方程的最終形式：綜上，整理得到能量方程的最終形式：,11Nijiti tikkk ikiiiiijNijiti tkikkkiiiiiijueuepuThYVtxxxxxxueueYpuTh Dtxxxxxxx 或超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分LES中

25、中能量方程能量方程Favre過濾，得到過濾，得到LES中的能量方程：中的能量方程：1111NijiikikkkiiiiiijNijiikikkiikiiiiiijiNkkkkkkkkiiiuEuYpuETh DtxxxxxxxuEuYpuETh DEuEutxxxxxxxxYYh Dh DxxxNiiijiijiijpupuuuxx？相比流動場來說，燃燒場中的能量方程出現(xiàn)非常相比流動場來說，燃燒場中的能量方程出現(xiàn)非常多的不封閉項，情況非常復雜多的不封閉項，情況非常復雜。超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分通常亞格子焓通量被?；ǔ喐褡屿释勘荒；簛喐褡诱承皂椡ǔ：苄?/p>

26、，予以忽略亞格子粘性項通常很小，予以忽略。sgsiiiiiHEuEupupuPrsgssgsisgsiHHx sgsijijiijiuu,11NNsgskkk ikkkkkkiiYYqh Dh Dxx亞格子熱通量通常亞格子熱通量通常也也很小，予以忽略很小，予以忽略。超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分組分方程組分方程kikkkkiiiYuYYDtxxxLES中中組分方程組分方程kikkkkkkikikkkiiiiiiiYu YYYYDu YuYDDtxxxxxxx？超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分組分方程的不封閉項：組分方程的不封閉項：通常通

27、常亞格子亞格子對流組分對流組分通量被模化通量被?；?sgssgskk iikikk sgsiYYu YuYScxkkkkiiYYDDxx通常通常也忽略也忽略kk才是真正的麻煩，才是真正的麻煩，由于過濾后的化學反應由于過濾后的化學反應源項的高度非線性本質而顯得格外復雜，源項的高度非線性本質而顯得格外復雜，如何?；且淮箅y點，值得研究。如何?；且淮箅y點，值得研究。超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬小結小結個人觀點個人觀點：1、超音速燃燒中大渦模擬的不封閉項更多、超音速燃燒中大渦模擬的不封閉項更多，但最難的問題是，但最難的問題是反應源項反應源項的?；哪；?，此問題值得此問題值得深入研究。深入研究。2、我們另辟蹊徑，繞開這個問題，不解能、我們另辟蹊徑，繞開這個問題，不解能量方程，量方程，基于基于層流穩(wěn)態(tài)火焰面假設，然后建立湍層流穩(wěn)態(tài)火焰面假設，然后建立湍流火焰面的庫，通過查庫插值得到過濾的流火焰面的庫，通過查庫插值得到過濾的反應源反應源項、過濾的組分質量分數、過濾的溫度等參數，項、過濾的組分質量分數、過濾的溫度等參數，大大減少了計算量，但計算精度肯定受到影響。大大減少了計算量，但計算精度肯定受到影響。超音速湍流燃燒大渦模擬超音速湍流燃燒大渦模擬第三部分第三部分k1、超音速流動與燃燒相比