下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
有衰減擴散問題的偏微分模型引言物質(zhì)的擴散問題,在石油開采、環(huán)境污染、疾病流行、化學(xué)反應(yīng)、新聞傳播、煤礦瓦斯爆炸、農(nóng)田墑情、水利工程、生態(tài)問題、房屋基建、神經(jīng)傳導(dǎo)、藥物在人體內(nèi)分布以及超導(dǎo)、液晶、燃燒等諸自然科學(xué)與工程技術(shù)領(lǐng)域,十分普遍地存在著。凡與反映擴散有關(guān)的現(xiàn)象,大都能由線性或非線性拋物型偏微分方程作為數(shù)學(xué)模型來定量或定性地加以解決。問題的提出設(shè)有一擴散源,某物質(zhì)從此擴散源向四周擴散,沿x,y,z三個方向的擴散系數(shù)分別為常數(shù),衰減(例如吸收、代謝等)使質(zhì)量的減少與濃度成正比,擴散前周圍空間此物質(zhì)的濃度為零,估計物質(zhì)的分布。理想假設(shè)設(shè)u(x,y,z,t)是t時刻點(x,y,z)處某物質(zhì)的濃度。任取一個閉曲面S,它所圍的區(qū)域是0。模型的建立由于擴散,從t到t+At時刻這段時間內(nèi),通過S流入0的質(zhì)量為TOC\o"1-5"\h\z=ft+Atff(a2色cosa+b2竺cosp+c2色cosy)dSdt
t dx dy dz其中a2其中a2,b2,由高斯公式z方向的擴散系數(shù)。c2分別是沿x,yM=ft+Atfff(a1td2u ,dM=ft+Atfff(a1t2 +b2 +c2 )dxdydzdtdx2dy2 dz20由于衰減,0內(nèi)的質(zhì)量減少為=ft+AtIffk2udxdydzdtt0其中k2為衰減系數(shù)。由物質(zhì)不滅定律,在t到t+At時刻間0內(nèi)由于擴散與衰減的合作用,積存于0內(nèi)的質(zhì)量為M-M。12換一個角度看,在t到t+At時刻間0內(nèi)由于濃度的變化引起的質(zhì)量增加為M=fff[u(x,y,z,t+At)—u(x,y,z)]dxdydz =30=ft"呻dydzdt0顯然,M3=M1-M2,即0=ft+Atfff(a2d2u+b2d2u+c2d2u一k2uWdydzdtt dx2 dy2 dz20由At,t,0的任意性得:du d2u d2u d2u=a2 +b2 +c2 一k2udt dx2 dy2 dz2上述方程是常系數(shù)線性拋物型方程,它就是有衰減的擴散過程的數(shù)學(xué)模型。定解條件的提法設(shè)擴散源在點(x,y,z)處,則此擴散問題滿足Cauchy問000題:du d2u7d2u d2u7——=a2 +b2+c2一k2U<dt dx2dy2 dz2u(x,y,z,0)=m3(x-x)5(y-y)6(z-z)000其中M為擴散源的質(zhì)量。方程的解析解用傅立葉變換可求得Cauchy問題的解析解為( 八 M I(x-x)2(y-y)2(z-z)2u(x,y,z,t)= 嚴exp$— 0-o-”o-k2t卜8Ktabcy/F [ 4a2t 4b2t 4c2t ‘但值得注意的是,在實際應(yīng)用中,參數(shù)a,b,c,k往往是很難獲得的,通常都是利用觀測取樣值進行估計,從而得出u(x,y,z,t)的近似表達式。參數(shù)估計:目的是對上式中出現(xiàn)的參數(shù)a,b,c,k進行估計。已知條件:①點源(擴散源)的質(zhì)量M;點源(擴散源)的位置:(x,y,z);ooot0時刻的觀測取樣值(x,y,z,m),m為t時刻(x,y,z)iiiiio iii處物質(zhì)的濃度,i二1,…。n首先考慮取樣時刻。事實上,取樣時刻是未知的,但若設(shè)取樣時刻為t,o作變量替換t=t「則有t=t/1,從而oodududt du= =tdtdtdtodt即dud2ud2ud2u=ta2 +1b2 +1c2 一tk2udtodx2ody2odz2o上式仍然是常系數(shù)線性拋物型方程,與有衰減的擴散過程的數(shù)學(xué)模型形狀完全一致,故可令觀測取樣值的取樣時刻為t=1。于是,(x,y,z,m)滿足o iiii( [、M J(x-x)2(y-y)2(z-z)272u(x,y,z,1)= —exp]— 0-o- 0-k28兀tabgFt [ 4a2 4b2 4c2 ‘其次考慮參數(shù)估計。對上式兩端取對數(shù),有
a 1 1s=ln(M、-ln(abc)-k2則有關(guān)系式:W=lnu(x,y,z,1)=aX+RY+yZ+s由于我們獲得的觀測取樣值(x,y,z,m)可以轉(zhuǎn)化為相應(yīng)的觀測取樣值iiii+(z-z+(z-z0)24c2+k2(X-x)20-4(y一yo)2(Z-Zo)2(X,Y,Z,W),于是利用多元回歸分析可以求出a、隊丫、s的估計值,從而得iiii到參數(shù)a,b,c,k的
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 冀教版小學(xué)語文五年級下冊知識點
- 公修課先進制造業(yè)基礎(chǔ)知識考題附復(fù)習(xí)資料
- 2024聘用常年法律顧問合同
- 2024網(wǎng)絡(luò)設(shè)備安裝合同
- 2024企業(yè)之間的借款合同范本
- 均值濾波課程設(shè)計
- 土地整治泵房施工方案
- 圍墻滲水施工方案
- 商場店鋪門施工方案
- 咸陽消音片施工方案
- 大學(xué)禮堂、報告廳、會議室使用管理辦法
- 小腸穿孔的護理查房共42張課件
- 天津市新版就業(yè)、勞動合同登記名冊
- 少先隊名師工作室匯報課件
- 樁側(cè)摩阻力計算
- 安全管理人員任命書及安全組織架構(gòu)圖
- 標本留存登記表
- 熱能與動力工程測試技術(shù)- 流量測量
- (中職)導(dǎo)游講解技巧單元一情景一歡迎詞講解技巧教學(xué)課件
- 2023屆高考地理一輪復(fù)習(xí):第4講 地球自轉(zhuǎn)及其地理意義 作業(yè)(人教版)
- 混凝土進場檢驗(澆筑)記錄
評論
0/150
提交評論