【正文】 未知量恒取上游節(jié)點(diǎn)的函數(shù)值。 CIP 數(shù)值方法是根據(jù)方程（ ） ~（ ）編寫(xiě)的，而上風(fēng)差分格式是在方程（ ）和方程（ ）的基礎(chǔ)上編寫(xiě)的。 因此 通過(guò) CIP 數(shù)值方法算得的三角形的數(shù)值模型的區(qū)域面積與通過(guò)解析解得到的模型的區(qū)域面積近似相等，故 CIP 數(shù)值方法具有較高的守恒性。 CIP 數(shù)值方法是在方程（ ） ~（ ）的基礎(chǔ)上編寫(xiě)的，而上風(fēng)差分格式是根據(jù)方程（ ）和方程（ ）編寫(xiě)的。但 通過(guò) CIP 數(shù)值方法算得的數(shù)值結(jié)果因其計(jì)算數(shù)據(jù)中出現(xiàn)負(fù)值以及比解析解峰值濃度高的數(shù)值，故 CIP 數(shù)值方法存在一定的不穩(wěn)定性。 （ a） 三維 圖 （ b） y=x上的剖面圖 圖 初始條件 大連海洋大學(xué)本科畢業(yè)論文 數(shù)值模型的構(gòu)建及兩種方法的對(duì)比 12 (3)計(jì)算參數(shù)選取 對(duì)流速度 u=， 對(duì)流速度 v=， 空間步長(zhǎng) dx=1m， 空間步長(zhǎng) dy=1m， 時(shí)間步長(zhǎng) dt=。 通過(guò) 上風(fēng)差分格式法算得的圓錐的數(shù)值模型的區(qū)域面積近似為 ，而通過(guò)解析解得到的模型的區(qū)域面積近似為 ，因此 通過(guò)上 風(fēng)差分格式法算得的圓錐的數(shù)值模型的區(qū)域面積與通過(guò)解析解得到的模型的區(qū)域面積相差較大，這是由于數(shù)值在 y=x 剖面的 垂直方向上存在數(shù)值耗散，故上風(fēng)差分格式法在 y=x 剖面上不具有較高的守恒性。但在實(shí)際中，由于受到計(jì)算方法的限制，其結(jié)果發(fā)生了一定的變化，如圖 和 所示 ， 其誤差分析見(jiàn)表 和 。 （ a） CIP 數(shù)值方法 （ b） 上風(fēng)差分格式法 圖 圓錐和圓柱數(shù)值模型數(shù)值解的三維圖（ t=1s） 大連海洋大學(xué)本科畢業(yè)論文 數(shù)值模型的構(gòu)建及兩種方法的對(duì)比 17 （ a） CIP 數(shù)值方法 （ b） 上風(fēng)差分格式法 圖 y=50 剖面上， 描述 圓錐和圓柱 濃度輸運(yùn)的 數(shù)值模型 的數(shù)值解與解析解（ t=1s） 表 圓柱模型在 y=50 剖面上， 數(shù)值格式的區(qū)域面積與解析解之間的相對(duì)誤差 CIP 格式 上風(fēng)格式 解析解 區(qū)域面積 50 相對(duì)誤差 % % 0 表 圓錐模型在 y=50 剖面上， 數(shù)值格式的區(qū)域面積與解析解之間的相對(duì)誤差 CIP 格式 上風(fēng)格式 解析解 區(qū)域面積 14 相對(duì)誤差 % % 0 大連海洋大學(xué)本科畢業(yè)論文 結(jié)論與展望 18 第四章 結(jié)論與 展望 結(jié)論 為了對(duì)解決一些較復(fù)雜的對(duì)流擴(kuò)散問(wèn)題（尤其是對(duì)流問(wèn)題）提供 準(zhǔn)確可靠的 數(shù)值模型 ，本文 研究了一維和二維的 CIP 數(shù)值方法，并利用 C 語(yǔ)言編寫(xiě)程序，建立了基于 CIP 系列模式的物質(zhì)輸移擴(kuò)散模型。 首先我要衷心地感謝我的導(dǎo)師張瑞瑾副教授 ,是她引領(lǐng)我走進(jìn)了這個(gè)既有趣又富有挑戰(zhàn)性的科學(xué)領(lǐng)域。 最后，我要再次感謝大連海洋大學(xué)以及各位老師，感謝我的家人，感謝我的大學(xué)各位同學(xué)，感謝所有的出現(xiàn)在我的生命里的朋友們，你們的 鼓勵(lì)與支持，鞭策與教導(dǎo)都是我人生的寶貴財(cái)富 。授 人以魚(yú)不如授人以漁 ,張瑞瑾導(dǎo)師不僅使我學(xué)習(xí)到了先進(jìn)的科學(xué)知識(shí)，接受了全新的思想觀念 ,學(xué)會(huì)了基本的研究方法以及學(xué)術(shù)素養(yǎng) ,更使我明白了許多為人處世的道理 .張瑞瑾導(dǎo)師對(duì)我成長(zhǎng)的無(wú)微不至的關(guān)懷 ,使我終生難忘。 展望 本文建立 的 是 一維和二維的數(shù)值模型，而在實(shí)際生活中，計(jì)算某種或某些物質(zhì)輸移擴(kuò)散所需要的模型大都是三維 數(shù)值模型，因此今后的研究中應(yīng)著重研究 三維數(shù)值模型。 通過(guò) CIP 數(shù)值方法算得的圓柱的數(shù)值模型的區(qū)域面積（即圖形與 x 軸所圍成的面積，下同）近似為 ，圓錐的數(shù)值模型的區(qū)域面積近似為 ，而通過(guò)解析解得到的圓柱模型的區(qū)域面積近似為 ，圓錐模型的區(qū)域面積近似為 ，因此 通過(guò) CIP 數(shù)值方法算得的圓錐和圓柱的數(shù)值模型區(qū)域面積與通過(guò)解析解得到的模型的區(qū)域面積近似相等，故 CIP 數(shù)值方法在y=50 上具有較高的守恒性。 綜上所述，通過(guò) CIP 數(shù)值方法所得到的圓錐的二維數(shù)值模型與 通過(guò) 上風(fēng)差分格式法所得到的圓錐的二維數(shù)值模型相比，具有較高的準(zhǔn)確性與守恒性。在理論上，描述圓錐濃度輸運(yùn)的數(shù)值模型在純移流條件下形狀和峰值高度都不會(huì)發(fā)生改變。 通過(guò) 上風(fēng)差分格式法算得的矩形的數(shù)值模型的區(qū)域面積近似為，而通過(guò)解析解得到的模型的區(qū)域面積為 ，因此 通過(guò)上 風(fēng)差分格式法算得的矩形的數(shù)值模型的區(qū)域面積與通過(guò) 解析解得到的模型的區(qū)域面積相差較小，故上風(fēng)差分格式法具有較高的守恒性，但不如 CIP 數(shù)值方法守恒性高。 圖 初始條件 (3)計(jì)算參數(shù)選取 對(duì)流速度 u=， 空間步長(zhǎng) dx=， 時(shí)間步長(zhǎng) dt=。 如 圖 （ b）所示，通過(guò) 上風(fēng)差分格式法得到的數(shù)值結(jié)果與解析解的數(shù)值結(jié)果相比，其大連海洋大學(xué)本科畢業(yè)論文 數(shù)值模型的構(gòu)建及兩種方法的對(duì)比 8 模型形狀變化較大，不再呈現(xiàn)三角形形狀，呈現(xiàn)類(lèi)似拋物線形，峰值高度變化也較大，故上風(fēng)差分格式法不具有較高的準(zhǔn)確性。 圖 初始條件 (3) 計(jì)算參數(shù)選取 對(duì)流速度 u=， 空間步長(zhǎng) dx=， 時(shí)間步長(zhǎng) dt=。一階上風(fēng)差分格式在求解線性對(duì)流擴(kuò)散方程時(shí)，具有形式簡(jiǎn)單，計(jì)算方便，數(shù)值穩(wěn)定等優(yōu)點(diǎn)，但由于存在一定的數(shù)值耗散問(wèn)題，故其求解的精度不高，計(jì)算結(jié)果往往不能令人滿意。 二維 CIP 數(shù)值方法 對(duì)于二維的 CIP 數(shù)值方法，其偏微分方程為： 0????????? yfvxfutf （ ） 其中 u 為 x 軸方向的 移動(dòng)速度， v 為 y軸方向的 移動(dòng)速度，且 u和 v 為常數(shù)。 反之，當(dāng) u0 時(shí)，將函數(shù) F（ x）的解析式整理得： niniiini ffbaf ????? ??? 39。? () nifd? () 經(jīng)整理得： niniiini ffbaf ????? ??? 39。 大連海洋大學(xué)本科畢業(yè)論文 CIP 數(shù)值方法與上風(fēng)差分格式 法 2 第二章 CIP 數(shù)值方法與上風(fēng)差分格式法 CIP 數(shù)值方法 一維 CIP 數(shù)值方法 對(duì)于一維的 CIP 數(shù)值方法，其偏微分方程為： 0?????? xfutf （ ） 其中 u 為移動(dòng)速度?；谖锢砩辖?jīng)常遇到的守恒性方程，例如著名的符拉索夫方程，在標(biāo)準(zhǔn) CIP 數(shù)值方法的基礎(chǔ)上，又發(fā)展了守恒型的 CIPCSL2 以及 RCIPCSL2 ]11~9[ 。對(duì) 于 一些簡(jiǎn)單的對(duì)流擴(kuò)散問(wèn)題，可以通過(guò)尋找理論解來(lái)解決。河流排海污染物總量居高不下，陸源入海排污口達(dá)標(biāo)率僅為 52%。 因此， 很有必要研究污染物的在環(huán)境中的運(yùn)動(dòng)狀態(tài)。 11 純移流條件下 描述 圓錐體 濃度輸運(yùn) 模型的構(gòu)建及兩種方法的對(duì)比 9 二維數(shù)值模型的構(gòu)建 4 一維一階上風(fēng)差分格式法 II 第一章 前 言 I Abstract 7 一維數(shù)值模型的構(gòu)建 7 純移流條件下 描述 三角形 濃度輸運(yùn) 模型的構(gòu)建及兩種方法的對(duì)比 18 結(jié)論 19 參考文獻(xiàn) 20 大連海洋大學(xué)本科畢業(yè)論文 摘 要 I 摘 要 隨著現(xiàn)代工業(yè)的發(fā)達(dá)和人類(lèi)人口的增長(zhǎng)，各種污染物都被排放到大氣、河流和海洋中，導(dǎo)致了環(huán)境的破壞。近岸局部海域海水環(huán)境污染依然嚴(yán)重，春季、夏季和秋季劣于第四類(lèi)海水水質(zhì)標(biāo)準(zhǔn) ]2[ 的海域面積分別為 52280、 41140 和 57360 平方公里。在物理運(yùn)動(dòng)中，包括兩種形式的運(yùn)動(dòng)：對(duì)流和擴(kuò)散。 在國(guó)外， CIP 數(shù)值方法發(fā)展較成熟，并得到廣泛應(yīng)用，尤其在物理領(lǐng)域。主要是通過(guò)對(duì)對(duì)流擴(kuò)散方程，尤其是對(duì)流方程進(jìn)行差分離散，利用 C 語(yǔ)言編寫(xiě)程序建立模型來(lái)模擬可 溶性物質(zhì)的輸移擴(kuò)散過(guò)程，并同基于上風(fēng)差分格式給出的數(shù)值結(jié)果進(jìn)行對(duì)比，驗(yàn)證該方法的準(zhǔn)確性和守恒性。21 2)(3 () nifc 39。21 2)(3 () 大連海洋大學(xué)本科畢業(yè)論文 CIP 數(shù)值方法與上風(fēng)差分格式 法 3 其中， tu????