動(dòng)力學(xué)蒙特卡羅方法在多尺度模擬中的研究及應(yīng)用.pdf

1、先進(jìn)的能源系統(tǒng)具有高度復(fù)雜性和多尺度現(xiàn)象的特點(diǎn)，并要求具有極高效率和近零排放。一般情況下，應(yīng)用于宏觀現(xiàn)象的連續(xù)性數(shù)學(xué)描述在微觀尺度上具有一定的局限性，而在全部宏觀域中采用微觀尺度上的進(jìn)化規(guī)則如動(dòng)力學(xué)蒙特卡羅（kMC）、格子-Boltzmann（LB）或分子動(dòng)力學(xué)（MD）等方法則計(jì)算代價(jià)太過(guò)高昂。合理劃分系統(tǒng)尺度并各自采用有效的模擬算法能有效解決復(fù)雜系統(tǒng)存在的計(jì)算代價(jià)和組織復(fù)雜性這兩大瓶頸問(wèn)題。先進(jìn)能源利用系統(tǒng)在實(shí)現(xiàn)能量轉(zhuǎn)換和環(huán)境控制中大

2、多要涉及多種類型的化學(xué)反應(yīng)過(guò)程，對(duì)于在微觀尺度上模擬反應(yīng)-擴(kuò)散過(guò)程，kMC是一種計(jì)算量適中的有效算法。kMC方法從微觀尺度出發(fā)，采用隨機(jī)的方式模擬系統(tǒng)中微粒的運(yùn)動(dòng)，通過(guò)微觀物質(zhì)量的變化體現(xiàn)系統(tǒng)的宏觀性質(zhì)，為宏觀模擬提供重要的數(shù)據(jù)信息。
　　 本論文研究了基于微觀模擬技術(shù)的動(dòng)力學(xué)蒙特卡羅方法（kMC）的機(jī)理，并重點(diǎn)研究了kMC的一種具體實(shí)現(xiàn)形式——Gillespie隨機(jī)仿真算法。以固體氧化物燃料電池（SOFC）為研究對(duì)象，首先采用

3、Gillespie算法模擬了電化學(xué)反應(yīng)過(guò)程，然后采用kMC模擬了氧化釔穩(wěn)定氧化鋯（YSZ）電解質(zhì)中的氧空位的擴(kuò)散過(guò)程。kMC對(duì)SOFC的應(yīng)用是研究的重點(diǎn)，主要包括以下幾方面的內(nèi)容：
　　 1、研究了kMC微觀模擬方法中的Gillespie隨機(jī)仿真算法，采用此方法跟蹤SOFC電化學(xué)反應(yīng)中各物質(zhì)分子數(shù)目的變化，從而體現(xiàn)電化學(xué)反應(yīng)的動(dòng)力學(xué)特征。
　　 2、沿YSZ電解質(zhì)厚度方向建立了一維kMC格子模型，由隨機(jī)數(shù)確定氧空位在格子

4、上躍遷的步進(jìn)時(shí)間和發(fā)生的事件，跟蹤每層格子氧空位數(shù)目的變化，并不斷更新系統(tǒng)狀態(tài)。著重分析了溫度、電解質(zhì)初始結(jié)構(gòu)分布和電極兩端電壓頻率對(duì)電流密度、（由正電荷和電子在電極兩端聚集產(chǎn)生的第N層格子處的電勢(shì)）、（第N層格子處的總電壓）和氧空位分布的影響。結(jié)果具有概率統(tǒng)計(jì)的特點(diǎn)，準(zhǔn)確描述了系統(tǒng)的動(dòng)力學(xué)特征。
　　 本文初步嘗試了應(yīng)用kMC模擬能源利用系統(tǒng)中的微觀過(guò)程，以SOFC為例研究了電化學(xué)反應(yīng)和氧空位在電解質(zhì)中的擴(kuò)散過(guò)程，仿真結(jié)果表明