偏微分方程離散差分格式差分方法等ppt課件

1、三偏微分方程的數(shù)值離散方法三偏微分方程的數(shù)值離散方法 3.1 有限差分法 3.2 有限體積法 (有限元，譜方法，譜元，無(wú)網(wǎng)格，有限解析，邊境元，特征線 3.1 有限差分法 3.1.1 模型方程的差分逼近 3.1.2 差分格式的構(gòu)造 3.1.3 差分方程的修正方程 3.1.4 差分方法的實(shí)際根底 3.1.5 守恒型差分格式 3.1.6 偏微分方程的全離散方法3.1.1 模型方程的差分逼近3.1.2 差分格式的構(gòu)造3.1.3 差分方程的修正方程差分方程所準(zhǔn)確逼近的微分方程稱(chēng)為修正方程 對(duì)于時(shí)間開(kāi)展方程，利用展開(kāi)的方程逐漸消去帶時(shí)間的高階導(dǎo)數(shù)，只留空間導(dǎo)數(shù)。Warming-Hyett方法：差分方程(

2、2)寫(xiě)成算子的方式： (3) 24121616121 )2() 1() 1(! 31! 21) 1(! 31! 21Taylor(2) 22121(1) 04422233233222133322213332221112111xuxutcxuxxuctuttuttuueuuexuxxuxxuxuuuetuttuttutuuuuuuuuuxuctuxxjxxjjttnjnjjjjjjnjnj等價(jià)于：展開(kāi)3.1.3 差分方程的修正方程 (續(xù)) 12220121212112332114444444333344433322223332222(5) 22121183,31,21, 1,1) 1() 1(

3、211) 1( 21216122121) 1( ! 31! 21) 1( (4)u 221)(21) 1(ppppppppkkkkllxxxxxxxxllttlllttllttttttttttxxxxxxxxttxuxuxutueeeebebttbbbbebttetttutetuttutetuttuttutetuttuttuteeeueeue即有最后得到的級(jí)數(shù)表示成可以將則記算子3.1.3 差分方程的修正方程(續(xù)) ? (2) schemefor 1CFL why . min)3(81)(61, 020) 1( : , 0) 1() 1() 1(22224222321221212012212)

4、(122201212121穩(wěn)定性判別條件符合，）：對(duì)于（滿(mǎn)足偶次項(xiàng)系數(shù)件是格式穩(wěn)定的充分必要條基本解為HyettgWarxtctcxtccckkkkeexuxuxutuppppppppppppppikxtippppppppkkkk3.1.4 差分方法的實(shí)際根底 相容性，穩(wěn)定性，收斂性 等價(jià)性定理 Fourier穩(wěn)定性分析3.1.4 差分方法的實(shí)際根底續(xù) Fourier (Von Neumann) 穩(wěn)定性分析(1) 0 , 0)(111cuuxctuunininini1Gfactor ion amplificat )()(: ) 1 ( 111111111nnikxnikxnikxnikxnni

5、nininikxnniikxnniikxnniAAGeAeAeAeAuuuueAueAueAuxtciiiiiii滿(mǎn)足穩(wěn)定性要求的代入誤差的基本解設(shè)3.1.4 差分方法的實(shí)際根底續(xù) Fourier (Von Neumann) 穩(wěn)定性分續(xù) 稱(chēng)為CFL條件 (Courant, Friedrichs, Levy)1 1 1,2sin)1 (41sin)cos1 (1sin)cos1 (1)sin(cos1122222 ifGxkxkxkGxkixkxkixkeGxik3.1.5 守恒型差分格式流膂力學(xué)方程組描畫(huà)物理量的守恒性；守恒律組：定義01diixtfu)(),(:f),(2f 0)(, 212

6、121212111ufuuufuuufflffxtuuxuftunljnljnljnjnjnjnjnjnj滿(mǎn)足相容性條件個(gè)變量的多變量函數(shù)：稱(chēng)為數(shù)值通量，它是其中則為守恒型差分格式。下形式其差分格式如果具有如對(duì)于一維單個(gè)守恒律：3.1.5 守恒型差分格式續(xù)守恒性質(zhì)：非守恒的差分格式普通沒(méi)有對(duì)應(yīng)于原始守恒律的“離散守恒律。0),(),(),(),()0 ,(),(:2/12/1112/12/102102110210210121211dttxufdxtxudttxudttxudxxudxtxutftfxuxuntftfxuxujJJnnJJxxtJtJxxnNkkJNkkJJjJjjJjJjnjn

7、JnJJjJjnjJjJjnj守恒律：律。完全對(duì)應(yīng)于連續(xù)的該積分代表離散的守恒可以看成是積分求和再對(duì)求和守恒型差分格式對(duì)3.1.5 守恒型差分格式續(xù)守恒型差分格式的Lax-Wendroff定理： 假設(shè)守恒型差分格式是和守恒律相容的，且當(dāng)時(shí)間和空間步長(zhǎng)趨于零時(shí)，差分解一致有界，幾乎處處收斂于分片延續(xù)可微的函數(shù)，那么這個(gè)收斂的函數(shù)就是守恒律的一個(gè)弱解。推論：守恒型差分各式的收斂解能自動(dòng)滿(mǎn)足延續(xù)關(guān)系。 用途: (加上熵條件可以得到正確的激波，研討中大量運(yùn)用例如：Lax-Friedrichs 格式，Lax-Wendroff格式，Mac Cormack格式 0)(xuftunjnjnjnjffxtuu2

8、121113.1.6 偏微分方程的全離散方法 對(duì)差分格式的普通要求： 有精度、格式穩(wěn)定、求解效率高 特殊要求 物理定律守恒性、物理特征激波、湍流、旋渦、多介質(zhì)、化學(xué)反響等、有界性(正密度、正溫度、正湍動(dòng)能、正組分濃度等) 主要指非定常方程的時(shí)間離散 3.1.6偏微分方程的全離散方法(續(xù) 兩層格式 Crank-Nicolson格式、P-C格式、Lax-Wendroff格式、MacCormack 格式 Runge-Kutta方法 時(shí)空全守恒：如Godunov格式、central-upwind格式、CESE方法 多層格式 Leap-Frog格式、Adams-Bashforth格式、后三點(diǎn)隱格式 3.

9、1.6.1 兩層格式 Crank-Nicolson格式 Predictor-Corrector格式 Lax-Wendroff 格式 Mac Cormack格式 Runge-Kutta方法stable nalunconditio)(4440)(40)(2011111111111111111nnninininininininininininininnniniBAuBuuuuuuuuuuxctuuxuxuctuuxuctu3.1.6.1 兩層格式(cont. Lax-Wendroff 格式 一步LW格式xtcxtOuuuuuuuxuctuninininininini其中),( ),2(2)(2022

10、11211111) 1(cossin1)2(2)(21:2121GxkxkiAAGeeeeAAFouriernnxikxikxikxiknn穩(wěn)定性3.1.6.1 兩層格式(cont. Lax-Wendroff 格式 兩步LW格式 常系數(shù)Jacobian時(shí)與單步LW等價(jià)。但計(jì)算更簡(jiǎn)單，不涉及矩陣相乘。1 ),( , 0)(10)(12/2/ )(022212121211112121xtOffxtuuffxtuuuxftuninininininininini3.1.6.1 兩層格式(cont. Mac Cormack 格式 (1969) 兩步格式 比LW更簡(jiǎn)單，不需求計(jì)算函數(shù)在半點(diǎn)上的值。 LW兩

11、步格式和MC各式的缺陷：定常解的誤差依賴(lài)于時(shí)間步長(zhǎng)。1 ),(0)(1)(21:0)(1:022*1*11*xtcxtOffxtuuuCffxtuuPxftuiiinininininiiMac Cormack格式的構(gòu)造)(1)(1)(21)()(exp.Taylor )(210:0:0*1*1*1*1*1*1*iinininnavavnininniinininFFxtUFFxtUtUtUtUttUUUUUUxFFtUUCxFFtUUPxFtU3.1.6.2 三層格式 Leap-Frog格式 Adams-Bashforth格式)()()()(12)()()(2)()(2)(3123232221t

12、OtOufufttuftutOufttuftutOtuttutuuuftunnnnnnnnn)()(21)(23211tOufuftuunnnn第二課后閱讀提示 傅德薰，3.1 3.3 水鴻壽3.1 , Ferziger and Peric, Springer Chap. 6作業(yè)2 1.用Fourier法分析 3.1.6.1節(jié)中Crank-Nicolson格式的穩(wěn)定性。 2.分析前面3.1.6節(jié)中Mac Cormack格式是幾階精度。3.2有限體積法 出發(fā)方程為積分型守恒方程直角坐標(biāo)、柱坐標(biāo)、球坐標(biāo) 以控制體為離散量 計(jì)算體積分和面積分需求適當(dāng)?shù)牟逯倒胶头e分公式 (quadrature fo

13、rmula) 適用于任不測(cè)形的網(wǎng)格，復(fù)雜幾何外形 缺陷：難以構(gòu)造大于二階以上的格式3.2.1 定常守恒型方程和控制體dqdsdsSSnnvgrad3.2.2 面積分的逼近 面積分用積分點(diǎn)的值表示(quadrature) 積分點(diǎn)的值用CV的值表示(interpolation) 對(duì)于Simpson公式,對(duì)積分點(diǎn)的插值需求四階精度eSseeneeSeeeSkSSffffdSSfSffdSfdSfdSeek461:order4th :order-2nd3.2.4 體積分的逼近 當(dāng)被積函數(shù)為某種型函數(shù)時(shí)，可以得到準(zhǔn)確的積分，逼近精度取決于型函數(shù)的精度。是精確的積分時(shí)，為常值分布或線性分布幾何中心的值。為

14、PPPPqqCVqqqqd :order-2nd3.2.4 體積分的逼近四階精度：2D直角坐標(biāo)網(wǎng)格最后一式可以四階精度逼近3D的面積分3.2.5 插值和微分 積分點(diǎn)的函數(shù)值和其法向梯度 1st UDS: 取上風(fēng)點(diǎn)的值eeeSudSeSnv插值 2nd order: 向積分點(diǎn)線性插值 等價(jià)于中心差分 (CDS)插值 當(dāng)積分點(diǎn)的函數(shù)是線性插值時(shí) Second order 插值 QUICK (quadratic upwind interpolation for convective kinematics) 插值三階精度，但積分差分往往只需二階精度。插值 高精度： N階精度的quadrture需求N-

15、1階多項(xiàng)式插值公式。 界面上導(dǎo)數(shù)可以用插值公式的微分求出。3.2.5有限體積法的邊境條件 用邊境條件替代面積分 入口：通常給定對(duì)流通量 (mass, momentum, energy, etc.) 壁面和對(duì)稱(chēng)面：通量為零 邊境上函數(shù)值給定：和內(nèi)部CV的值共同構(gòu)建邊境上的導(dǎo)數(shù)FV例子3.2.6 守恒律的有限體積方法 Godunov 格式0),(),(0),(),(: t ,tx,xT0in 0)(21211212111nn1/2i1/2- idttxufdxtxudtdxtxufxdxdttxutuftuiinniinnxxttxxtt或積分在，問(wèn)題的精確解給出。的值是由其中得Riemann),

16、(),(1),(12/121212121112121txudttxuftgdxtxuxuggxtuuittiixxnniiinininnii3.2.6.1 Godunov方法的思想.),( Riemann,),; 0(),(Riemann,Riemann,. ,111 2121121121 2121ninniiinnnRPRLRPiniRinLiinintxxxtxttttttxxxxxxxxtiuuuuuuuuuuuuuu離散分布時(shí)刻的的就構(gòu)成下一輪循環(huán)所需得到的平均值內(nèi)進(jìn)行平均，在每個(gè)小網(wǎng)格區(qū)間精確解在下一時(shí)刻的值都互不干擾，這樣可將，的精確解問(wèn)題相鄰的局部使得每個(gè)小網(wǎng)格內(nèi)來(lái)自取問(wèn)題的精確

17、解并求出該問(wèn)題的考慮初始間斷為內(nèi)的常數(shù)分布小網(wǎng)格區(qū)間看成是時(shí)刻的離散分布量把已知的一階迎風(fēng)格式(CIR格式) 0 xuatu用Godunov思想闡明CIR格式=Godunov格式0 ),(11auutaxuxuninininiRiemann解圖示3.2.6.1 1D Euler方程組的Godunov格式0)(0)(00)()(0)()(0)(22dtpEuEdxdtpuudxudtdxxpEutExputuxut積分方程組微分方程組Godunov格式是基于積分方式的方程組,延續(xù)關(guān)系自動(dòng)滿(mǎn)足，不需求另外思索延續(xù)線上的延續(xù)關(guān)系挪動(dòng)網(wǎng)格上的積分回路11njx11njxnjx1njx1njx1njx1

18、2njxnjx21ntnt挪動(dòng)網(wǎng)格上的Godunov格式 求出。由狀態(tài)方程后，在求出),(,0)()()()(0)()()()(0)()()()(121121121111111121211111211211111112121111121121111121111121epppEupuDuEpuDuEtxxExxEpDuupDuutxxuxxuDuDutxxxxnjnjnjjjjjjjjjjjjjnjnjnjnjnjnjnjnjjjjjjjjjjjnjnjnjnjnjnjnjnjjjjjjjnjnjnjnjnjnj固定網(wǎng)格上的Godunov格式 上為同一個(gè)值）問(wèn)題的解在同一條射線（的射線上的物理量問(wèn)題的解在斜率處局部為（Riemann0Riemann),000, 02111111212112112121111121211211212111211211111wxpuxupuEupuEtEExpuupuutuuxuutxxxxDjjjjjjjjjjjjjjjnjnjnjnjjjjjjjjjjnjnjnjnjjjjjjnjnjnjnjnjn