【VIP專享】第四章--進排氣管內非定常流的數(shù)值分析2課件

1、第四章 進、排氣管內非定常流的數(shù)值分析4.1 非定常一維流動的基本方程4.2 特征線法的基本概念4.3 特征線方程和相容性方程4.4 特征線方程的數(shù)值解法4.5 邊界條件概述4.6 排氣閥的邊界條件4.7 渦輪邊界條件4.8 分支管接頭的邊界條件概述內燃機的進、排氣過程決定了氣缸充氣量的多少、多缸機各缸換氣的均勻性，是決定內燃機性能優(yōu)劣的關鍵因素之一，所以內燃機中對進排氣非定常流動特性的研究很有必要。流體力學中早已闡明了非定常流場中壓力波傳播原理，并推導了計算流場狀態(tài)參數(shù)的偏微分基本方程組。1885年黎曼（Riemann)提出求解雙曲型偏微分方程的特征線原理與解法之后，非定常流的工業(yè)應用研究才

2、得以開始。20世紀50年代海勒（Haller）和詹尼（Jenny）將特征線圖解法應用于內燃機中，可較精確模擬管道中實際情況，并據(jù)以改進和優(yōu)化設計，但仍因計算工作量過大而難以推廣應用。直到20世紀60年代電子計算機廣泛推廣應用。20世紀60年代以來，英國曼切斯特大學的本森（R. S. Benson）教授及其工作組長期致力于內燃機進排氣及增壓系統(tǒng)的非定常流研究和特征線模擬計算工作，70年代開發(fā)了UMIST程序，可用于24缸以下機型。1982年，他的同事以本森的名義出版內燃機的熱力學與空氣動力學，系統(tǒng)總結了以特征線法為基礎的這一領域的研究和計算成果前后30年1960-1990年代初，特征線法一直在內

3、燃機一維非定常流研究和工業(yè)應用中處于主導支配地位，期間德國賽福爾特教授1978年完成了Promo模擬計算程序，國內清華大學80年代初引進UMIST程序移植到微機上MK-14程序。以特征線法為基礎的內燃機一維非定常氣流的研究，主要集中于兩個方面：一是進行循環(huán)模擬計算，用以指導進氣、排氣、增壓系統(tǒng)及整機的實驗、設計，改善和提高內燃機的綜合性能；二是利用壓力波傳播特性，開發(fā)和創(chuàng)建新型的裝置和系統(tǒng)，比如，利用壓力波傳輸能量的脈沖增壓系統(tǒng)、調諧增壓系統(tǒng)以及將排氣壓力波能量直接送到新鮮充氣中的氣波增壓裝置，這些都是在充分掌握非定常流中壓力波傳播原理之后的創(chuàng)新之作。近一二十年來，這一領域的研究工作發(fā)展趨勢：

4、特征線法逐漸被求解方便、精度更高、速度更快的直接求解偏微分方程組的有限差分法和有限體積法所取代，一維向考慮介質真實的粘性、傳熱和湍流等無形特點，特別是多管接頭還使用了多維模型；過渡過程的動態(tài)模擬及非定常流的反饋控制等新領域。盡管模擬計算已由特征線法轉為直接針對偏微分方程的更精確的有限差分法，但特征線法原理及解法在一維非定常流研究中仍具有重要而不可取代的作用。一方面，在管網(wǎng)系統(tǒng)的各個準定常流的邊界處，仍需利用進入邊界的壓力單波等單波線方程來聯(lián)立求解；另一方面，特征線法的物理概念及圖解分析法，是直觀而形象地描述和理解管內壓力過程的有效途徑，也是學習和一維非定常流原理，并進行創(chuàng)新性工作不可缺少的工具

5、。正是基于這一原因，包括本書在內的有關內燃機數(shù)值模擬書籍都對特征線原理和解法進行了詳細的介紹。4.1 非定常一維流動的基本方程進排氣管內流動一般作一維流動處理計算時假定（1）管軸向幾何尺寸比徑向的大得多，軸向流動效應比徑向的大得多，可以認為每一截面上流動參數(shù)是均勻的，僅隨軸向坐標和時間變化，準一維流動（2）管壁剛性，截面積變化緩慢（3）用一維模型模擬摩擦、傳熱（4）管內流體為完全氣體，不計重力4.1 非定常一維流動的基本方程取控制容積ABCD，根據(jù)上述假定，就可從 質量守恒定律 動量定理和 能量守恒定律出發(fā)，來推導管內不定常非等熵流動的基本方程 一、非定常一維不等熵流動1連續(xù)方程2動量方程2動

6、量方程（續(xù)）G的物理意義：單位流體質量分配到沿x方向的比摩擦力。3能量方程4熵方程（129頁）（4.1.31）考察流體質點通過控制表面時的熵值變化沿軌跡線代入能量方程熱力學第二定律有二、非定常一維等熵流動對于等截面非定常一維等熵流動，S=常數(shù)，沒有摩擦與傳熱（q=0，G=0）連續(xù)方程動量方程能量方程等熵條件下，能量方程為音速方程，自動滿足p，可由熱力學公式由a表示，只有兩個獨立參數(shù)a，u。這組方程的特點是：1.一階偏微分方程組2.非線性 (系數(shù)中含未知函數(shù))3.系數(shù)中不含偏導數(shù)， 一階擬線性偏微分方程組4.雙曲線型偏微分方程組，可用特征線求解二、非定常一維等熵流動（續(xù)）等熵條件狀態(tài)方程利用上述

7、公式可將p，用a來表示對上式先取對數(shù)再求偏導又有將以上四式代入連續(xù)方程和動量方程，得連續(xù)方程動量方程4.2 特征線法的基本概念特征線定義：從物理觀點來看 特征線是物理擾動的傳播路徑，馬赫線是超聲流場中特征線從數(shù)學觀點來看， 沿特征線能將偏微分方程換成全微分方程，穿過特征線時，物理參數(shù)的偏導可以不連續(xù)（不能唯一確定），而物理參數(shù)本身卻保持連續(xù)。 利用這一特性可以用來尋求特征線下面用一柯西問題的例子，進一步說明特征線的概念4.2 特征線法的基本概念（續(xù)）4.2 特征線法的基本概念（續(xù)）4.2 特征線法的基本概念特征線定義：從物理觀點來看 特征線是物理擾動的傳播路徑從數(shù)學觀點來看， 沿特征線能將偏微

8、分方程換成全微分方程，穿過特征線時，物理參數(shù)的偏導可以不連續(xù)（不能唯一確定），而物理參數(shù)本身卻保持連續(xù)。 利用這一特性可以用來尋求特征線用一柯西問題的例子，進一步說明特征線的概念4.2 特征線法的基本概念（續(xù)）特征線方程相容性方程4.2 特征線法的基本概念（續(xù)）考慮非定常一維等熵、等截面流動聯(lián)立r1方向u的全微分r2方向a的全微分連續(xù)方程動量方程4.2 特征線法的基本概念（續(xù)）考慮非定常一維等熵流動（續(xù)）黎曼變量非定常一維等熵流動相容性方程解的積分常數(shù) 4.2 特征線法的基本概念（續(xù)）考慮非定常一維等熵流動（續(xù)）123物理意義4.3 特征線方程和相容性方程一、特征線方程和相容性方程的推導4.3

9、 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）一、特征線方程和相容性方程的推導（教材將方程組作線性組合） 00根據(jù)特征線定義，若沿特征線有全微分特征線：相容性方程： 4.3 特征線方程和相容性方程（續(xù)）二、特征線方程實用化轉換二、特征線方程實用化轉換（續(xù)）消去消去4.3 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）4.3 特征線方程和相容性方程（續(xù)）U ，A沿特征線將偏微分方程求解問題轉換成常微分方程求解問題最后由，AA求得p，u 待求點P點的參數(shù)

10、由其依賴區(qū)決定已知量由d d dAA方程確定4.4 特征線方程的數(shù)值解法特征線解法：（1）圖解法：在X-Z，U-A平面上畫圖求解（2）數(shù)值解法 用有限差分方法對特征方程和相容性方程進行離散化，然后對有限差分方程數(shù)值求解 數(shù)值求解又包括（1）直接步進法；（2）逆步進法 不同方法采用不同差分網(wǎng)格例1 如圖，有一根一端封閉，另一端全部開口的等截面管道。設在開口端有一個正弦波形的壓力波進人管道，現(xiàn)要求求出壓力波在封閉管端經(jīng)一次反射后，在靠近開口管端處的壓力波形。解 按以下步驟在XZ平面上構筑特征線：首先，對X和Z軸選取作圖比例尺，即X和Z坐標軸的單位長度分別為s和r (mm)，于是對某一網(wǎng)格點，在圖上

11、量出的X的實際長度msX，同樣Z在圖上的實際長度為n=rZ，圖中s/r =2圖解法（等熵等截面）（1）壓力波未進入管道前，管道內的起始值為U=0，A1, 可知=1（2）正弦波進入管內發(fā)出特征線，由于此時無左行波，因此保持起始值 =1，特征線值的求法如下： 取Z=0.2，即2A點，在壓力波圖上量出p/pR=1.14， 故A=a/aR=(p/pR)-1/(2)=(1.14)0.143=1.0189 =2A-=21.0189-1=1.038 tan +=2(3 -2 )=2(31.038-21)=2.228 +=6583（與Z軸夾角） 由此可以作出2A至2,1點的特征線。（3）右管端為一閉端，U=0

12、，因此=，由此可求出邊界點1,1、2,2、3,3、4,4、5,5點上的值。（4）特征線相交點的計算：例如3,2，其值與3,1點相同，值與2,2點相同，即3,2= 3,1=1.052, 3,2= 2,2=1.038。（5）左端管口處B1至B5點處壓力波的計算：由于此處不存在右行波，故4.4 特征線方程的數(shù)值解法特征線解法：（1）圖解法：在X-Z，U-A平面上畫圖求解（2）數(shù)值解法 用有限差分方法對特征方程和相容性方程進行離散化，然后對有限差分方程數(shù)值求解 數(shù)值求解又包括（1）直接步進法；（2）逆步進法 不同方法采用不同差分網(wǎng)格一、直接步進法（特征線網(wǎng)格） 二、逆步進法（正交網(wǎng)格法）二、逆步進法（

13、正交網(wǎng)格法）（續(xù)）三、計算步長的選取和穩(wěn)定性準則由左圖，s點的斜率為三、計算步長的選取和穩(wěn)定性準則（續(xù)）由于沿x軸可以分為許多網(wǎng)格，對每一網(wǎng)格點均有組A和U值，因此對整個X抽滿足穩(wěn)定性準則的Z值為用圖示的方法來說明 Z最小值的選取方法，圖上過P4 點的特征線形成最小的Z，以此作為計算的時間步長，就可保證在求m+1行上各點的和值時，各交點均可落在m行相應的間隔(如圖上虛線和點劃線所示)之內，以保證計算的穩(wěn)定。三、計算步長的選取和穩(wěn)定性準則（續(xù)）不同流動狀態(tài)下特征線的相互位置（145圖4.4.4）四、特征線的斜率四、特征線的斜率U0U0超聲速流： U0dAA第四章 進、排氣管內非定常流的數(shù)值分析4

14、.1 非定常一維流動的基本方程4.2 特征線法的基本概念4.3 特征線方程和相容性方程4.4 特征線方程的數(shù)值解法4.5 邊界條件概述4.6 排氣閥的邊界條件4.7 渦輪邊界條件4.8 分支管接頭的邊界條件沿X軸將計算區(qū)間分為間隔為X網(wǎng)格，由每一網(wǎng)格點的A和U值及穩(wěn)定性準則，確定下一時間步長Z值4.4 特征線方程的數(shù)值解法AAAA(dAA)T(d)R(d)S4.5 邊界條件概述 前面介紹已知初值條件求解管道內網(wǎng)格點流動參數(shù)的數(shù)值方法，對于管端邊界上的點，只有附加給出邊界約束條件，才能求得確定解 4.5 邊界條件概述（續(xù)）4.6 排氣閥的邊界條件一、排氣正流一、排氣正流（續(xù)）思路：由缸內滯止狀態(tài)求得排氣管狀態(tài)參數(shù)（p，u）方程 消去喉口中間變量(at,ut)（氣缸至喉口為等熵 流動）氣缸到喉口處的流動為等熵流兩邊乘p/pc/u得到未知， 二、排氣閥倒流（亞聲速）二、排氣閥倒流（亞聲速）二、排氣閥倒流（亞聲速）三、排氣閥無流