譜元方法課程教學(xué)文案

2023年2月7日Tel：82663537e-mail：glqin@11、正交函數(shù)(hánshù)系與譜近似1.1正交函數(shù)系與正交多項(xiàng)式1.2函數(shù)的Fourier展開(zhǎnkāi)1.3Chebyshev譜逼近(離散函數(shù)的Fourier展開(zhǎnkāi))1.4插值函數(shù)的導(dǎo)數(shù)1.5二維函數(shù)的Chebyshev譜近似第一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@21.1正交函數(shù)(hánshù)系與正交多項(xiàng)式設(shè)函數(shù)系{j(x)}，如滿足條件，則稱函數(shù)系{j(x)}，在區(qū)間(qūjiān)[a,b]上關(guān)于權(quán)函數(shù)w(x)正交，當(dāng)j(x)是代數(shù)多項(xiàng)式時(shí)，稱為正交多項(xiàng)式。如Chebyshev多項(xiàng)式，

第二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@3Chebyshev多項(xiàng)式的正交性第三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@41.2函數(shù)(hánshù)的Fourier展開函數(shù)f(x)在區(qū)間[a,b]上滿足(mǎnzú)Dirichlet條件，且加權(quán)平方可積，即對(duì)于權(quán)函數(shù)w(x)，有，則f(x)可在[a,b]上以w(x)為權(quán)函數(shù)，按正交多項(xiàng)式φn(x)展成廣義Fourier級(jí)數(shù)，設(shè){j(x)}為在節(jié)點(diǎn){xk,k=0,1,…,N}上的正交函數(shù)系，權(quán)為{wk>0,k=0,1,…,N}，設(shè)f(x)為在節(jié)點(diǎn)系{xk,k=0,1,…,N}上取值的已知函數(shù)，則函數(shù)f(x)的Fourier展開定義為，第四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@51.3Chebyshev譜逼近(離散函數(shù)(hánshù)的Fourier展開)取節(jié)點(diǎn)系{xk,k=0,1,…,N}為N階Chebyshev多項(xiàng)式的極值(jízhí)點(diǎn)(稱Gauss-Lobatto點(diǎn))，即，權(quán)系數(shù)則上面的展開式變?yōu)?，此式亦可看作f(x)在配置點(diǎn)上的插值,g(x)為插值函數(shù)。第五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@61.4插值函數(shù)(hánshù)的導(dǎo)數(shù)對(duì)微分得：插值函數(shù)在配置點(diǎn)的導(dǎo)數(shù)(dǎoshù)可表示為，第六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@71.5二維函數(shù)(hánshù)的Chebyshev譜近似對(duì)二維函數(shù)u(ξ,η)，在標(biāo)準(zhǔn)(biāozhǔn)正方形單元內(nèi)，定義ξ,η方向節(jié)點(diǎn)系{ξj，j=0,1,…,Nx}和{ηk，k=0,1,…,Ny}，第七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@82譜方法(fāngfǎ)求解微分方程2.1解微分方程的加權(quán)殘量法MWR（MethodofWeightResiduals)2.2譜方法(fāngfǎ)求解微分方程2.3小結(jié)第八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@92.1解微分方程(wēifēnfānɡchénɡ)的加權(quán)殘量法MWR（MethodofWeightResiduals)基本(jīběn)思想MWR的步驟MWR的基本(jīběn)方法第九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@10加權(quán)殘量法基本(jīběn)思想(1)設(shè)微分方程(wēifēnfānɡchénɡ)邊值問題為，加權(quán)余量法，是求微分方程(wēifēnfānɡchénɡ)形式如下的近似解。a=[a0,a1,…,aN]T為待定向量，0,1,…,N為為基函數(shù)，且是一組線性無關(guān)的函數(shù)第十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@11加權(quán)殘量法基本(jīběn)思想(2)“殘量”定義如下，顯然只有當(dāng)試函數(shù)為邊值問題的精確解時(shí)，余量(2.31)才為零。加權(quán)殘量法適當(dāng)選取待定向量a=[a0,a1,…,aN]T，使得(shǐde)“殘量”極小。通常某一權(quán)函數(shù)系{wi(x),i=0,1,…,N}，使得(shǐde)權(quán)函數(shù)與“殘量”正交，來確定待定向量。

第十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@12MWR的步驟(bùzhòu)①選取一組滿足要求的基函數(shù)(hánshù)②構(gòu)造試函數(shù)(hánshù)并滿足所有邊界條件；③選取一組權(quán)函數(shù)(hánshù)④運(yùn)用MWR準(zhǔn)則，得到的代數(shù)方程組；⑤解上述代數(shù)方程組，確定待定參數(shù)。第十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@13MWR的基本(jīběn)方法按權(quán)函數(shù)的不同有以下幾種基本方法Galerkin法(GalerkinMethod)配置(pèizhì)法(CollocationMethod)Tan方法(TanMethod)類似于Galerkin法，Tan方法中，試函數(shù)不需要滿足邊界條件，增加邊界約束系數(shù)來滿足邊界條件

第十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@142.2譜方法(fāngfǎ)求解微分方程以非線性Burger方程為例來說明譜方法的離散過程(guòchéng),并區(qū)別Galerkin法,Tau方法和Collocation方法。Burger方程如下:

第十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@15ChebyshevGalerkin方法(fāngfǎ)(GalerkinMethod方程的近似解表達(dá)(biǎodá)為，由MWR準(zhǔn)則得,第十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@16ChebyshevTau方法(fāngfǎ)(TanMethod方程的近似解表達(dá)為，由MWR準(zhǔn)則(zhǔnzé)得,邊界約束,第十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@17ChebyshevCollocation方法(fāngfǎ)(CollocationMethod)選取節(jié)點(diǎn)(jiédiǎn)系為N階Chebyshev多項(xiàng)式的極值點(diǎn),方程的近似解用下式表示，由Collocation方程得,第十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@182.3小結(jié)(xiǎojié)一種求解偏微分方程的高階精度數(shù)值方法。屬于求解偏微分方程加權(quán)殘量法的特例。譜方法使用高階正交多項(xiàng)式作為展開函數(shù)。譜方法最受人青睞的優(yōu)越性在于它具有“無窮階”的收斂(shōuliǎn)速度，其確切含義為，若原微分方程的解無限可微，則由適當(dāng)?shù)淖V方法所得到的近似解對(duì)原問題的收斂(shōuliǎn)速度比1/N的任何冪都更快，這里N是所取基函數(shù)的個(gè)數(shù)。Kreiss和Oliger,OrszagGottlieb和Orszag第十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@193譜元方法(fāngfǎ)求解微分方程回顧譜方法———使用加權(quán)殘量法譜元方法基本(jīběn)思想常微分方程邊值問題的Galerkin變分原理偏微分方程的Galerkin變分原理Galerkin逼近解常微分方程元素矩陣的形成偏微分方程元素矩陣的形成極坐標(biāo)系下環(huán)形區(qū)域的譜元方法總剛度矩陣的形成第十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@203.1回顧譜方法(fāngfǎ)———使用加權(quán)殘量法

選取一組滿足要求的基函數(shù)

構(gòu)造試函數(shù)，滿足所有的邊界條件

選取一組權(quán)函數(shù)

利用MWR準(zhǔn)則，得到ai的代數(shù)方程

求解上述代數(shù)方程組，確定待定參數(shù)第二十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@213.2譜元方法基本(jīběn)思想使用有限元方法的思想，將求解(qiújiě)域分成若干子域采用譜方法中譜近似技術(shù)代替有限元中的插值函數(shù)采用有限元中等參單元,亦可逼近復(fù)雜求解(qiújiě)邊界基函數(shù)的選用在一個(gè)單元上選取在整個(gè)區(qū)域上選取計(jì)算精度不是太高比較高單元的選擇三角形，四邊形等四邊形提高精度方法網(wǎng)格加密不改變插值次數(shù)增加基函數(shù)多項(xiàng)式的次數(shù)第二十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@223.3常微分方程(wēifēnfānɡchénɡ)邊值問題的Galerkin變分原理考察(kǎochá)二點(diǎn)邊值問題上述問題的Galerkin變分問題是：求，使得，其中第二十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@233.4偏微分方程(wēifēnfānɡchénɡ)的Galerkin變分原理考察二維橢圓邊值問題上述(shàngshù)問題的Galerkin變分問題是：求，使得，第二十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@243.5Galerkin逼近(bījìn)解設(shè)Vh是V的有限維子空間，當(dāng)h->0時(shí)，Vh的維數(shù)無限增加，直到充滿V為止，那么Galerkin逼近解，使得，設(shè)是Vh的基函數(shù)系，那么由于，故有，其中得由于的任意性，得，即得線性代數(shù)(xiànxìnɡdàishù)方程組，可以證明此代數(shù)方程組的解唯一存在，即可解出Galerkin逼近解。第二十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@253.6常微分方程元素矩陣(jǔzhèn)的形成考察(kǎochá)以下微分方程用譜方法求解此微分方程,先將求解區(qū)域分成若干單元,設(shè)單元總數(shù)為N,考察(kǎochá)第i個(gè)單元ei,將ei轉(zhuǎn)化為標(biāo)準(zhǔn)求解單元,設(shè)x為總體坐標(biāo),為標(biāo)準(zhǔn)求解單元內(nèi)的局部坐標(biāo),并設(shè)第二十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@26插值點(diǎn)和插值函數(shù)(hánshù)在第i個(gè)元素內(nèi)選取Nx+1個(gè)點(diǎn)作為(zuòwéi)u(x)和f(x)的插值點(diǎn),本文選取階Chebyshev多項(xiàng)式的極值點(diǎn)作為(zuòwéi)插值點(diǎn),即,u(x)和f(x)可以表達(dá)為，第二十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@27元素(yuánsù)矩陣形成對(duì)一維邊值問題,在ei單元內(nèi)線性方程式變?yōu)槠渲?qízhōng)，上面線性方程組變?yōu)?，第二十七頁，?3頁。2023年2月7日Tel：82663537e-mail：glqin@28的推導(dǎo)(tuīdǎo)過程第二十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@293.7偏微分方程元素矩陣(jǔzhèn)的形成具體考察Helmholtz方程，其中，作以下坐標(biāo)變換在標(biāo)準(zhǔn)(biāozhǔn)正方形單元內(nèi)，使用偽譜Chebyshev譜逼近u(ξ,η)，則元素內(nèi)線性方程組式變?yōu)?，第二十九頁，?3頁。2023年2月7日Tel：82663537e-mail：glqin@30總體(zǒngtǐ)坐標(biāo)系與局部坐標(biāo)系第三十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@31元素(yuánsù)矩陣第三十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@323.8極坐標(biāo)系下環(huán)形區(qū)域(qūyù)的譜元方法在極坐標(biāo)系下，Helmholtz方程表達(dá)為，將極坐標(biāo)系下環(huán)形段單元，變換為直角坐標(biāo)(zhíjiǎozuòbiāo)系下的正方形單元，坐標(biāo)變換關(guān)系如下，元素矩陣的形式為，第三十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@33極坐標(biāo)坐標(biāo)系與局部(júbù)坐標(biāo)系第三十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@34元素(yuánsù)矩陣式中各變量(biànliàng)的表達(dá)式如下，第三十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@353.9總剛度(ɡānɡdù)矩陣的形成節(jié)點(diǎn)影響元素集，節(jié)點(diǎn)所在的元素的并集影響節(jié)點(diǎn)，對(duì)某點(diǎn)，其影響元素內(nèi)的所有節(jié)點(diǎn)相互影響元素集，兩個(gè)節(jié)點(diǎn)和影響元素之交集對(duì)節(jié)點(diǎn)和，總體矩陣元素為，由于在一個(gè)(yīɡè)元素e內(nèi)計(jì)算，故稱它為單元矩陣，上式表明矩陣有貢獻(xiàn)的只是()影響元素集所有單元矩陣相應(yīng)的元素。第三十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@364復(fù)雜區(qū)域(qūyù)的譜元方法采用譜逼近的等參數(shù)單元總體坐標(biāo)系下導(dǎo)數(shù)(dǎoshù)的計(jì)算坐標(biāo)變換矩陣及變換行列式局部坐標(biāo)系下偏導(dǎo)數(shù)(dǎoshù)的計(jì)算總體坐標(biāo)系下偏導(dǎo)數(shù)(dǎoshù)的計(jì)算復(fù)雜求解區(qū)域的譜元方法第三十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@37總體坐標(biāo)(zuòbiāo)與局部坐標(biāo)(zuòbiāo)系第三十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@384.1采用譜逼近(bījìn)的等參數(shù)單元總體坐標(biāo)與局部坐標(biāo)(如圖所示)等參單元，插值函數(shù)表達(dá)式和坐標(biāo)變換(biànhuàn)表達(dá)式一致函數(shù)的插值表達(dá)為，坐標(biāo)變換(biànhuàn)的表達(dá)式為，導(dǎo)數(shù)表達(dá)式第三十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@394.2總體(zǒngtǐ)坐標(biāo)系下導(dǎo)數(shù)的計(jì)算坐標(biāo)變換(biànhuàn)矩陣(Jacobi矩陣)及變換(biànhuàn)行列式(Jacobi行列式)微元面積(miànjī)、微元弧長(為常數(shù))時(shí)，(為常數(shù))時(shí)第三十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@40局部坐標(biāo)系下偏導(dǎo)數(shù)(dǎoshù)的計(jì)算由插值式對(duì)ξ求偏導(dǎo)得，上式在(ξm,ηn)處的值為，比照函數(shù)得，則局部坐標(biāo)系下的偏導(dǎo)數(shù)(dǎoshù)表達(dá)為，第四十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@41總體坐標(biāo)系下偏導(dǎo)數(shù)(dǎoshù)的計(jì)算由插值函數(shù)得，由Jacobi矩陣得總體(zǒngtǐ)坐標(biāo)系下，偏導(dǎo)數(shù)計(jì)算式為，第四十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@424.3復(fù)雜(fùzá)求解區(qū)域的譜元方法考察Helmholtz方程(fāngchéng)則元素矩陣為，第四十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@43元素矩陣(jǔzhèn)各變量表達(dá)式第四十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@445非定常不可壓縮流動(dòng)與自然對(duì)流(duìliú)換熱問題求解不可壓縮流動(dòng)控制方程封閉空腔內(nèi)自然對(duì)流換熱基本(jīběn)控制方程解非定常不可壓縮Navier-Stokes方程的時(shí)間分裂法時(shí)間分裂法求解自然對(duì)流問題時(shí)間分裂法求解Navier-Stokes方程的誤差分析第四十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@455.1非定常不可壓縮流動(dòng)控制(kòngzhì)方程非定常不可壓縮流動(dòng)(liúdòng)無因次Navier-Stokes方程為，以頂蓋移動(dòng)的方腔驅(qū)動(dòng)流為例，定義無因次參量如下，第四十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@465.2封閉空腔(kōnɡqiānɡ)內(nèi)自然對(duì)流換熱基本控制方程封閉空腔(kōnɡqiānɡ)內(nèi)自然對(duì)流換熱基本控制方程第四十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@47無因次量定義(dìngyì)對(duì)正方形封閉空腔(kōnɡqiānɡ)，無因次參數(shù)的定義如下，對(duì)空氣來說，第四十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@485.3解非定常不可壓縮(yāsuō)Navier-Stokes方程的時(shí)間分裂法Navier-Stokes可寫成如下形式(xíngshì)，對(duì)時(shí)間積分得，設(shè)為三步中間值，則分三步求解的過程為：非線性步、壓力步、粘性步第四十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@49非線性步非線性步使用顯式4階Runge-Kutta方法，替代(tìdài)傳統(tǒng)時(shí)間分裂法的Adams-Bashforth方法，使本時(shí)間步僅與上一時(shí)間步有關(guān)，此步中代入初始條件，不使用邊界條件。第四十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@50壓力(yālì)步假設(shè)滿足不可壓縮條件(tiáojiàn)，且在邊界上法向上的投影為零，即，得,此步為解帶有第二類邊界條件(tiáojiàn)的壓力Poisson方程，解出壓力Poisson方程后，再由上式求解。第五十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@51粘性(zhānxìnɡ)步利用Crank-Nicolson格式離散粘性步，得Helmholtz方程，施以(shīyǐ)固體壁面邊界條件即可求解第五十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@525.4時(shí)間分裂法求解自然(zìrán)對(duì)流問題求解(qiújiě)自然對(duì)流問題的過程同Navier-Stokes方程，也是將一個(gè)時(shí)間步分為三步求解(qiújiě)，所不同的是在一個(gè)時(shí)間步內(nèi)除了解Navier-Stokes外，還需求解(qiújiě)能量方程，能量方程中沒有壓力項(xiàng)，因此在第二步壓力Poisson方程求解(qiújiě)中與能量方程無關(guān)。第五十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@53非線性步非線性步仍使用(shǐyòng)顯式4階Runge-Kutta方法，第五十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@54壓力(yālì)步為解帶有第二類邊界條件的壓力Poisson方程(fāngchéng)，解出壓力Poisson方程(fāngchéng)后，再由(5.34)式求解。第五十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@55粘性(zhānxìnɡ)步利用(lìyòng)Crank-Nicolson格式離散粘性步，得Helmholtz方程，施以固體壁面邊界條件即可求解。第五十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@565.5時(shí)間(shíjiān)分裂法求解Navier-Stokes方程的誤差分析利用時(shí)間分裂法求解Navier-Stokes方程基于以下(yǐxià)兩個(gè)假設(shè)，滿足不可壓縮約束，即滿足連續(xù)性方程，在邊界的法線上的投影為零，高精度壓力邊界條件描述，壓力Poisson方程變?yōu)椋谖迨?，?3頁。2023年2月7日Tel：82663537e-mail：glqin@576面向?qū)ο笞V元方法(fāngfǎ)程序設(shè)計(jì)面向?qū)ο蟪绦蛟O(shè)計(jì)思想面向?qū)ο笞V元方法程序設(shè)計(jì)節(jié)點(diǎn)類標(biāo)準(zhǔn)單元類單元類整體區(qū)域類譜元方法程序設(shè)計(jì)中的幾個(gè)問題插值函數(shù)積分的存儲(chǔ)節(jié)點(diǎn)的編碼及總體矩陣的形成(xíngchéng)Navier-Stokes方程和能量方程的求解過程第五十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@586.1面向?qū)ο蟪绦蛟O(shè)計(jì)(chénɡxùshèjì)思想按人們通常的思維方式建立模型(móxíng)，設(shè)計(jì)盡可能自然地表現(xiàn)求解方法的軟件。為此，必須建立直接表現(xiàn)問題中事物與事物相互聯(lián)系的概念，建立適合人們思維方式的描述方法；在面向?qū)ο蟮姆椒ㄖ兄饕幸韵赂拍睿簩?duì)象、消息傳遞、類和繼承、方法等。對(duì)象和消息傳遞是表示事物與事物間的關(guān)系；類和繼承是適應(yīng)人們思維方式對(duì)事物的描述方法；方法是作用在對(duì)象上的各種操作。第五十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@596.2面向?qū)ο笞V元方法(fāngfǎ)程序設(shè)計(jì)面向?qū)ο蟮某绦蛟O(shè)計(jì)包括面向?qū)ο蠓治觥⒚嫦驅(qū)ο笤O(shè)計(jì)和面向?qū)ο髮?shí)現(xiàn)對(duì)譜元方法來說，從譜元方法分析問題的思路出發(fā)，找出描述譜元方法的對(duì)象類（數(shù)據(jù)和方法的結(jié)合體），建立對(duì)象類的關(guān)系，并集成這些對(duì)象。設(shè)計(jì)節(jié)點(diǎn)類、標(biāo)準(zhǔn)單元、單元類、總體區(qū)域類。這些類中用屬性來描述類的特性，用方法來實(shí)現(xiàn)某些操作。有些(yǒuxiē)類中的屬性可能又需要用對(duì)象類來描述，方法的作用對(duì)象也可能是某一類的實(shí)體。第五十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@60節(jié)點(diǎn)(jiédiǎn)類主要的數(shù)據(jù)(shùjù)及方法包含節(jié)點(diǎn)的屬性，如節(jié)點(diǎn)的編號(hào)、節(jié)點(diǎn)的位置坐標(biāo)、節(jié)點(diǎn)的邊界屬性、微分方程右端項(xiàng)及邊界條件值、所要求解的變量及其導(dǎo)數(shù)值。節(jié)點(diǎn)類定義，classCGridPoint{protected://節(jié)點(diǎn)的屬性標(biāo)準(zhǔn)節(jié)點(diǎn)，編碼，坐標(biāo)，邊界屬性，方程右端項(xiàng)；壓力，速度及其導(dǎo)數(shù)；public://節(jié)點(diǎn)類方法屬性設(shè)置，邊界設(shè)置，壓力和速度設(shè)置；屬性取值，邊界屬性取值；}第六十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@61標(biāo)準(zhǔn)單元類包括了在各求解單元中的常量(chángliàng)。類中的數(shù)據(jù)及方法主要有：單元內(nèi)網(wǎng)格的劃分(含兩方向的節(jié)點(diǎn)數(shù)及總節(jié)點(diǎn)數(shù))、標(biāo)準(zhǔn)插值節(jié)點(diǎn)、插值函數(shù)的導(dǎo)數(shù)矩陣及有關(guān)標(biāo)準(zhǔn)插值函數(shù)的積分。標(biāo)準(zhǔn)單元類定義，classCStdElement{protected://標(biāo)準(zhǔn)單元類的屬性網(wǎng)格屬性，標(biāo)準(zhǔn)插值節(jié)點(diǎn)，插值函數(shù)的導(dǎo)數(shù)，插值函數(shù)積分；public://標(biāo)準(zhǔn)單元類方法標(biāo)準(zhǔn)單元類的構(gòu)造，Chebyshev多項(xiàng)式及其導(dǎo)數(shù)；插值函數(shù)的積分，屬性設(shè)置及取值；}第六十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@62單元(dānyuán)類主要數(shù)據(jù)及方法包含元素序號(hào)、單元內(nèi)節(jié)點(diǎn)劃分(和標(biāo)準(zhǔn)單元內(nèi)劃分相對(duì)應(yīng))、標(biāo)準(zhǔn)單元類、節(jié)點(diǎn)類和總體坐標(biāo)與局部坐標(biāo)的變換關(guān)系。單元類定義，classCElement{protected://單元類的屬性(shǔxìng)單元及網(wǎng)格屬性(shǔxìng)，標(biāo)準(zhǔn)單元類，坐標(biāo)變換，單元矩陣；public://單元類方法單元構(gòu)造，標(biāo)準(zhǔn)單元類的構(gòu)造，節(jié)點(diǎn)類構(gòu)造；元素矩陣的形成，屬性(shǔxìng)設(shè)置及取值；}第六十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@63總體(zǒngtǐ)區(qū)域類主要數(shù)據(jù)及方法包含區(qū)域內(nèi)單元數(shù)、總節(jié)點(diǎn)數(shù)、節(jié)點(diǎn)類、標(biāo)準(zhǔn)單元類、總體矩陣等?？傮w區(qū)域類定義，classCDomain{protected://整體區(qū)域類的屬性區(qū)域的屬性，單元及網(wǎng)格屬性，單元類，標(biāo)準(zhǔn)單元類；總體矩陣，對(duì)角元位置，邊界(biānjiè)節(jié)點(diǎn)序號(hào)；public://整體區(qū)域類方法單元構(gòu)造，標(biāo)準(zhǔn)單元類的構(gòu)造，節(jié)點(diǎn)類構(gòu)造；元素矩陣的形成，屬性設(shè)置及取值；}第六十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@646.2譜元方法程序設(shè)計(jì)(chénɡxùshèjì)中的幾個(gè)問題插值函數(shù)積分的存儲(chǔ)節(jié)點(diǎn)的編碼及總體矩陣的形成對(duì)稱帶狀矩陣的一維存儲(chǔ)元素的編碼影響元素和影響節(jié)點(diǎn)對(duì)角(duìjiǎo)元的位置及一維數(shù)組的長度總剛度矩陣第一類邊界條件的處理第六十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@65插值函數(shù)(hánshù)積分的存儲(chǔ)在形成元素矩陣時(shí)，要計(jì)算如下積分，

可以看出，

在(i,j,k)變化時(shí)，不必占用(Nx+1)3個(gè)存儲(chǔ)單元，這樣既節(jié)約了存儲(chǔ)空間又節(jié)省(jiéshěng)了計(jì)算時(shí)間。Bijk所占的空間可由下式計(jì)算，Bijk的檢索按如下方式進(jìn)行，先將i,j,k按由大到小排列，設(shè)i>j>k，則Bijk在數(shù)組中的位置為，第六十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@666.3Navier-Stokes方程和能量方程的求解(qiújiě)過程構(gòu)造標(biāo)準(zhǔn)單元(dānyuán)stdElement，求出對(duì)各單元(dānyuán)不變的量。構(gòu)造求解域pGeometry，對(duì)求解域進(jìn)行網(wǎng)格劃分，同時(shí)構(gòu)造節(jié)點(diǎn)類pDot[NG]，構(gòu)造各求解單元(dānyuán)pElement[Ne]，對(duì)節(jié)點(diǎn)類中U、V、Φ賦初值，賦零值，Φ在熱邊界上賦值1，在冷邊界上賦0，其余兩邊為第二類邊界條件，在形成單元(dānyuán)矩陣時(shí)自動(dòng)滿足。利用Runge-Kutta法求解、、，將求解結(jié)果存于、、解壓力Poisson方程，第二類邊界條件，解出壓力后更新、、求解關(guān)于、、的Helmholtz方程，邊界條件第一或第二類。解出Helmholtz方程后更新、、，即完成了一個(gè)時(shí)間步的求解。第六十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@677計(jì)算(jìsuàn)實(shí)例常微分方程求解(qiújiě)Helmholtz和Poisson的求解(qiújiě)簡單邊界Helmholtz方程求解(qiújiě)復(fù)雜邊界Poisson的求解(qiújiě)移動(dòng)頂蓋方腔驅(qū)動(dòng)非定常流動(dòng)封閉方腔中的自然對(duì)流第六十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@687.1常微分方程(wēifēnfānɡchénɡ)求解【例1】用譜方法求解下面方程該問題(wèntí)的精確解的表達(dá)式為：將求解域[0,-1]等分為10個(gè)求解單元,單元內(nèi)使用5階Chebyshev多項(xiàng)式進(jìn)行譜逼近。數(shù)值計(jì)算結(jié)果與精確解的比較如下,比較表明譜元方法具有相當(dāng)高的精度。第六十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@69例1數(shù)值(shùzí)計(jì)算結(jié)果與精確解的比較如下序號(hào)x數(shù)值解解析解絕對(duì)誤差0123456789100.01-3.83295107661418e-0180.1559246312085910.3082929811551880.4535843017646150.5883485562994740.7092408909868030.8130550556913160.8967554397068940.957507402575330.99270559687679100.1559246312085960.3082929811551960.4535843017646250.5883485562994880.7092408909868240.8130550556913390.8967554397069170.9575074025753540.9927055968768061-3.83295107661418e-018-4.88498130835069e-015-8.04911692853239e-015-9.99200722162641e-015-1.45439216225896e-014-2.04281036531029e-014-2.33146835171283e-014-2.38697950294409e-014-2.34257058195908e-014-1.62092561595273e-014-2.22044604925031e-016第六十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@70【例2】【例2】用譜元方法求解下面方程該問題的精確解的表達(dá)式為：將求解域[0,-1]等分為4個(gè)求解單元,單元內(nèi)使用5階Chebyshev多項(xiàng)式進(jìn)行譜逼近。數(shù)值計(jì)算結(jié)果與精確解的比較如下(rúxià),比較表明譜元方法具有相當(dāng)高的精度第七十頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@71例2數(shù)值(shùzí)計(jì)算結(jié)果與精確解的比較如下0123456789100.01-8.3294087758479e-0190.0147662965508840.0286795455946950.04087816195710.0504833997185180.0565905580149630.0582599254786710.0545073764403450.0442945201363430.0265183090547675.337544225295e-018

00.0147662965612980.0286795455705450.0408781618890680.050483399757920.0565905580149630.0582599255415590.0545073763178480.0442945199587630.0265183091540420

-8.3294087758479e-019-1.041408452529e-0112.4150008381962e-0116.8031476285757e-011-3.9402599238958e-011-3.816391647149e-016-6.2888715068876e-0111.2249663089436e-0101.7757908338245e-010-9.9275830611756e-0115.337544225295e-018

第七十一頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@727.2簡單邊界Helmholtz方程(fāngchéng)求解具體考察Helmholtz方程，精確(jīngquè)解為，采用矩形單元采用曲邊四邊形單元第七十二頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@73兩種單元求解結(jié)果比較(bǐjiào)(簡單邊界)表7-1(采用矩形單元(dānyuán)的求解結(jié)果)表7-2(采用曲邊四邊形單元(dānyuán)的求解結(jié)果)第七十三頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@74復(fù)雜邊界(biānjiè)Poisson的求解考察方程邊界條件由精確解給出采用(cǎiyòng)環(huán)形段單元采用(cǎiyòng)曲邊四邊形單元第七十四頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@75兩種單元求解結(jié)果(jiēguǒ)比較(復(fù)雜邊界)表7-3(采用環(huán)形(huánxínɡ)段單元的求解結(jié)果)表7-4(采用曲邊四邊形單元的求解結(jié)果)第七十五頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@767.3移動(dòng)(yídòng)頂蓋方腔驅(qū)動(dòng)非定常流動(dòng)是計(jì)算流體力學(xué)典型數(shù)值實(shí)驗(yàn)?zāi)Ｐ?，在很多情況下都是用它來檢驗(yàn)數(shù)值方法(fāngfǎ)的有效性；除移動(dòng)頂蓋外，初始速度場(chǎng)為零；移動(dòng)頂蓋上速度u=1,v=0，其余邊界速度為零；取9個(gè)求解單元，使用8階Chebyshev多項(xiàng)式插值；計(jì)算了Re=100的情況，取兩種時(shí)間步長T1=0.004，T2=0.005，計(jì)算2000和1600步，各需12和9.6小時(shí)第七十六頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@77移動(dòng)(yídòng)頂蓋方腔驅(qū)動(dòng)流物理模型第七十七頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@78豎直中心線上的速度(sùdù)U,水平中心線上的速度(sùdù)V圖7-3豎直(shùzhí)中心線上的速度U圖7-4水平(shuǐpíng)中心線上的速度V第七十八頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@79方腔內(nèi)中心點(diǎn)的速度(sùdù)U和速度(sùdù)V隨時(shí)間的變化圖7-5方腔內(nèi)中心點(diǎn)的速度(sùdù)U圖7-6方腔內(nèi)中心點(diǎn)的速度(sùdù)V第七十九頁，共93頁。2023年2月7日Tel：82663537e-mail：glqin@mail.