• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    蒙特卡羅模擬CFETR中子輸運(yùn)計(jì)算中的全局減方差方法應(yīng)用及對(duì)比

    2016-04-19 05:48:22聶星辰趙平輝祝慶軍中國科學(xué)技術(shù)大學(xué)核科學(xué)技術(shù)學(xué)院合肥3007中國科學(xué)院等離子體物理研究所合肥3003
    核技術(shù) 2016年3期
    關(guān)鍵詞:蒙特卡羅

    聶星辰 李 佳 趙平輝 祝慶軍(中國科學(xué)技術(shù)大學(xué) 核科學(xué)技術(shù)學(xué)院 合肥 3007)(中國科學(xué)院等離子體物理研究所 合肥 3003)

    ?

    蒙特卡羅模擬CFETR中子輸運(yùn)計(jì)算中的全局減方差方法應(yīng)用及對(duì)比

    聶星辰1李 佳1趙平輝1祝慶軍2
    1(中國科學(xué)技術(shù)大學(xué) 核科學(xué)技術(shù)學(xué)院合肥 230027)
    2(中國科學(xué)院等離子體物理研究所合肥 230031)

    摘要在具有全局特性的蒙特卡羅輸運(yùn)精細(xì)計(jì)算的問題中,傳統(tǒng)的MCNP(Monte Carlo N Particle Transport Code)局部減方差方法很難得到理想的計(jì)算結(jié)果,全局減方差方法(Global Variance Reduction,GVR)則是一種有效的解決方法。針對(duì)中國聚變工程試驗(yàn)反應(yīng)堆(Chinese Fusion Engineering Testing Reactor,CFETR)的中子輸運(yùn)過程中減小全局方差的問題,將多種形式的GVR方法應(yīng)用到柱狀CFETR中子學(xué)模型的計(jì)算中。依據(jù)不同的中子分布信息,在算例中應(yīng)用和對(duì)比了6種不同形式的GVR權(quán)窗,并對(duì)不同GVR方法的品質(zhì)因子(FOMG)、標(biāo)準(zhǔn)差(σ)和有效計(jì)數(shù)率(Scoring)進(jìn)行了分析。與AN(MCNP analog method)相比,GVR方法的FOMG有很大的增長,誤差在空間的分布也更加平緩,且具有更高的Scoring。在前人提出的全局減方差的基礎(chǔ)上,在計(jì)算中應(yīng)用一些新的GVR形式(能量、徑跡數(shù)等),計(jì)算結(jié)果表明,基于中子通量的GVR方法的全局計(jì)算效率較AN提高了6.43倍。此外,基于中子能量的全局減方差方法也是一種可行的GVR應(yīng)用形式,其與AN比較,計(jì)算效率提高了5.11倍。綜上,基于中子通量的GVR方法具有最佳的全局減方差效果。

    關(guān)鍵詞全局減方差,權(quán)窗,蒙特卡羅,MCNP,中子學(xué)

    磁約束核聚變國家特殊項(xiàng)目(No.2013GB108004、No.2015GB108002)、國家自然科學(xué)基金(No.11175207)資助

    第一作者:聶星辰,男,1990年出生,2013年畢業(yè)于南華大學(xué),現(xiàn)為碩士研究生,核能科學(xué)與工程專業(yè)

    Supported by National Special Project for Magnetic Confined Nuclear Fusion Energy(No.2013GB108004,No.2015GB108002)and the National Natural Science Foundation of China(No.11175207)

    Frist author:NIE Xingchen,male,born in 1990,graduated from University of South China,master student,major in nuclear science and engineering

    Application and comparison of global variance reduction methods employed
    in Monte Carlo neutron transport for CFETR

    NIE Xingchen1LI Jia1ZHAO Pinghui1ZHU Qingjun2
    1(School of Nuclear Science and Technology,University of Science and Technology of China,Hefei 230027,China)2(Institute of Plasma Physics,Chinese Academy of Sciences,Hefei 230031,China)

    AbstractBackground:For the detailed Monte Carlo radiation transport,which possesses global characteristics,traditional MCNP Localized Variance Reduction(LVR)methods can hardly perform well.Global Variance Reduction(GVR)method is an effective way.Purpose:This study aims to decrease the global variance of Monte Carlo radiation transport calculations for Chinese Fusion Engineering Testing Reactor(CFETR).Methods:According to different information of neutron distribution,six kinds of GVR methods were implied and compared in a cylinder CFETR neutronics model.Global figure of merit(FOMG),standard deviation(σ)and the percentage of scoring were analyzed.Based on predecessors’ research,other kinds of new GVR methods(energy,tracks)were applied in this calculation.Results:In the point of GVR methods,FOMGincreased dramatically,spatial variation inrelative errors was more gradual,and the scoring was higher,compared with AN.The GVR method based on neutron flux achieved the global calculation efficiency of 6.43 times that of AN.Apart from that,the GVR method based on neutron energy can also be a practical GVR kind,achieving the global calculation efficiency of 5.11 times that of AN.Conclusion:The GVR method based on neutron flux can performed best that has a well function of global variance reduction,achieving the global calculation efficiency.

    Key wordsGlobal variance reduction,Weight window,Monte Carlo,MCNP,Neutronics

    針對(duì)中國聚變工程試驗(yàn)反應(yīng)堆(Chinese Fusion Engineering Testing Reactor,CFETR)在整個(gè)空間中的高分辨中子輸運(yùn)計(jì)算和停堆劑量率分析的需求越來越大。在該工程的中子學(xué)分析中,由于CFETR具有龐大的體積和大量的屏蔽材料,而由此引發(fā)的中子深穿透問題會(huì)致使蒙特卡羅模擬(Monte Carlo,MC)很難給出理想結(jié)果。就此類具有全局特性的中子輸運(yùn)問題[1-6],需要給出中子通量的時(shí)間-空間-能量的精細(xì)分布。故而提出了如何獲取可靠全局解的問題。

    MC模擬結(jié)果的精度與樣本數(shù)、方差有關(guān)。通過增大樣本數(shù)或減小方差都可以提高計(jì)算精度。依靠大規(guī)模并行機(jī)的優(yōu)勢,可以增加樣本數(shù),但由于模型的特殊性所帶來的中子深穿透問題,故僅依靠增大樣本數(shù)不能從根本上解決中子通量在全局的精細(xì)分布計(jì)算問題[7]。而提高精度的另一條途徑是減方差。

    通過在MC中使用減方差技巧,可以有效提升中子對(duì)于計(jì)數(shù)器的貢獻(xiàn)率,提高模擬的計(jì)算效率。一種常用和有效的減方差方法是使用權(quán)窗[4-6]。權(quán)窗的作用是依據(jù)重要性調(diào)整中子的權(quán)重,通過設(shè)定權(quán)窗的參數(shù)可以有效控制區(qū)域內(nèi)的中子權(quán)重,進(jìn)而可以避免由輸運(yùn)過程中的權(quán)震蕩而造成的統(tǒng)計(jì)結(jié)果誤差過大的情況。對(duì)于具有全局特性的MC輸運(yùn)計(jì)算問題,傳統(tǒng)的MCNP局部減方差方法(Localized Variance Reduction,LVR)[5-6]中的權(quán)窗技巧由于針對(duì)的是特定的計(jì)數(shù)器或目標(biāo)空間區(qū)域,因而LVR很難得到理想的全局計(jì)算結(jié)果。

    為解決具有全局特性的MC輸運(yùn)計(jì)算問題,本文介紹了全局減方差方法(Global Variance Reduction,GVR)[1-4]。GVR主要功能是將中子輸運(yùn)到模型的全局空間中,使得MC結(jié)果的相對(duì)誤差在全空間中的每處都相對(duì)理想,可以有效減少模型全局的計(jì)算方差。本文依據(jù)MC的中子分布初始信息來獲取GVR權(quán)窗和重要性,并對(duì)比和分析了不同形式的GVR權(quán)窗在CFETR全局特性的MC輸運(yùn)計(jì)算中的效率。

    1 全局減方差原理

    相比于傳統(tǒng)的LVR方法,GVR方法在整個(gè)模型空間中都可以獲得較可信的中子通量精細(xì)分布,具有很好的全局減方差效果。當(dāng)感興趣區(qū)域是全局空間或目標(biāo)計(jì)數(shù)器分布在全空間中時(shí),采用GVR方法會(huì)使得測量結(jié)果的相對(duì)誤差隨空間分布的梯度盡可能小,進(jìn)而提高整個(gè)空間模型中的平均探測效率。為獲取全局空間里中子通量的分布,通過迭代的手法,以提高系統(tǒng)內(nèi)的重要性分布的清晰度,以此指導(dǎo)MC中子輸運(yùn)的徑跡走向,最終達(dá)到全局降方差的目的。

    這里需要指出的是為獲取全空間里可信的中子通量分布,可以針對(duì)全空間的目標(biāo)計(jì)數(shù)器,使用MCNP權(quán)窗生成器來自動(dòng)產(chǎn)生針對(duì)全空間的權(quán)窗。但是對(duì)于一個(gè)占有較大空間的復(fù)合靶柵元,直接使用權(quán)窗生成器具有一定的風(fēng)險(xiǎn)性。根據(jù)權(quán)窗重要性理論,該方法的直接后果是生成的權(quán)窗只會(huì)對(duì)中子較為容易輸運(yùn)到的空間起效果,而對(duì)其他區(qū)域抽樣不均勻。相比之下GVR方法是使用初始的MC模擬的相關(guān)中子信息來設(shè)定權(quán)窗的下閾值,進(jìn)而通過MCNP 內(nèi)置的WWP(Weight window parameter)卡構(gòu)建出GVR權(quán)窗。Wth為權(quán)窗下邊界,Wsplit為權(quán)窗上邊界,定義兩者比值β,如式(1)所示,其在MCNP5中默認(rèn)值是5。

    對(duì)于基于中子通量φ的GVR方法,在低中子通量和高計(jì)算誤差的區(qū)域,權(quán)窗具有較低的下邊界,使得該區(qū)域有更多的中子發(fā)生分裂,增加中子數(shù)的樣本以降低誤差,提高計(jì)算精度。而在高通量和低誤差區(qū)域,權(quán)窗具有相對(duì)高的下邊界,使得在該區(qū)域有更多的中子參與了俄羅斯輪盤賭,降低樣本以減少計(jì)算時(shí)間,而增加計(jì)算效率。所以在GVR權(quán)窗的作用下,中子的密度即會(huì)在低通量和高誤差區(qū)域相對(duì)升高。特別的,對(duì)于固定中子源的計(jì)算模型,其源所在的區(qū)域往往是最高中子通量Max(φ)所在的區(qū)域,如本算例中Max(φ)在等離子體區(qū)內(nèi)。在柵元i處基于中子通量φi的下邊界定義如下[2]:

    對(duì)于基于GVR方法的權(quán)窗或重要性值,其中子分布信息可以從標(biāo)準(zhǔn)MCNP輸出文中所得。使用該分布信息(可以用中子通量、中子數(shù)量、權(quán)重等形式)來產(chǎn)生基于柵元或網(wǎng)格的GVR權(quán)窗,也可以依據(jù)重要性值和權(quán)窗閾值之間的反比關(guān)系,由式(2)反推出重要性分布,如式(3)。將GVR方法應(yīng)用到IMP(Importance)卡中,通過IMP的幾何分裂和俄羅斯輪盤賭原理來實(shí)現(xiàn)全局減方差目的也是種可行的途徑。

    由最初的AN方法來獲取的中子分布信息在中子高衰減區(qū)域可能會(huì)缺失,但它為以后提供了迭代依據(jù)。每次迭代都會(huì)改善中子信息的空間分布情況。在GVR方法中使用多群和能量截?cái)嗟姆椒梢杂行岣叩俾蔥1]?;谥凶油康腉VR方法的實(shí)踐步驟如圖1。通過判斷在目標(biāo)區(qū)域內(nèi)是否有足夠抽樣作為停止迭代的依據(jù)。初始能量截?cái)嗟闹递^高(10-3MeV),以便縮短初始迭代的計(jì)算時(shí)間,而使整體的迭代過程有較快起步。在迭代的進(jìn)程中逐步降低能量截?cái)嗟闹担梢蕴岣呙總€(gè)中子歷史的平均壽命,而獲得空間中更加精細(xì)的中子分布。

    圖1 基于中子通量的GVR方法的驟流程圖Fig.1 Flow chart of GVR method based on neutron flux.

    需注意在低通量區(qū)域,毗連區(qū)域的通量衰減很快,當(dāng)直接采用式(2)算得到的權(quán)窗具有很高的梯度變化,這種非常危險(xiǎn)的減方差行為會(huì)使得中子的權(quán)重變化巨大,出現(xiàn)中子過度分裂的現(xiàn)象,而大大增加了一個(gè)中子歷史的計(jì)算時(shí)間。特別在并行計(jì)算的情況下,該現(xiàn)象會(huì)嚴(yán)重惡化計(jì)算效率。為緩解這種由于過度偏倚現(xiàn)象所帶來的中子輸運(yùn)的長歷史問題,本文通過設(shè)定通量的誤差閾值,而對(duì)通量的結(jié)果進(jìn)行篩選,當(dāng)通量的相對(duì)誤差Re超過限值時(shí)將關(guān)閉權(quán)窗或者統(tǒng)一地設(shè)定權(quán)窗值,而沒有使用全部的通量結(jié)果進(jìn)行計(jì)算。通過數(shù)值實(shí)驗(yàn)表明這種方法可以防止權(quán)窗在空間中有較大的振蕩,有效緩解了長歷史問題。

    2 全局減方差評(píng)估參數(shù)

    MCNP中的 FOM公式作用是衡量給定的計(jì)數(shù)器的計(jì)算效率,而對(duì)于GVR方法的計(jì)算效率則是反應(yīng)在整個(gè)空間分布中的計(jì)數(shù)器的效率。由于和傳統(tǒng)的LVR作用機(jī)制和使用目的不同,這里采用GVR的品質(zhì)因子FOMG[1]。對(duì)于計(jì)算效率較好的GVR方法,其具有較少的計(jì)算時(shí)間和較小的誤差統(tǒng)計(jì),從而具有較高的FOMG值。TCPU代表模擬的CPU計(jì)算時(shí)間,N為柵元或網(wǎng)格的總個(gè)數(shù),Rei代表第i個(gè)柵元的中子通量的相對(duì)誤差。FOMG的定義如下:

    基于相對(duì)誤差的標(biāo)準(zhǔn)差σ[2]可以有效衡量通量的相對(duì)誤差在空間位置上的變化情況。對(duì)于GVR方法的計(jì)算中,在全局空間上的模擬過程應(yīng)盡可能表現(xiàn)出較小的誤差分布梯度,使得隨著中子數(shù)歷史的增加,整個(gè)空間每處的方差以相同的速率收斂。理想GVR方法的σ應(yīng)盡可能趨近0。其定義如下:

    對(duì)于零中子通量的柵元,由于MCNP在該區(qū)域內(nèi)沒有統(tǒng)計(jì)結(jié)果,所以該柵元處的通量誤差是零,該零通量柵元就不適用于FOMG和σ的計(jì)算公式。對(duì)于在計(jì)算過程中出現(xiàn)的大量零通量柵元,其會(huì)影響GVR計(jì)算效率的評(píng)估。對(duì)于零通量柵元,可將相對(duì)誤差由0修改為1,使得相對(duì)誤差在空間中盡可能連續(xù),而后可代入FOMG和σ的計(jì)算公式。

    MCNP的有效計(jì)數(shù)率(Scoring)可以有效反映出沒有計(jì)數(shù)的零通量柵元的份額,故Scoring也是衡量GVR計(jì)算的一個(gè)重要指標(biāo)。綜上所述,當(dāng)FOMG、σ和Scoring的參數(shù)都可以接受時(shí),則該GVR方法的全局減方差效果較好。

    3 GVR方法的應(yīng)用和比較

    本文將多種GVR形式應(yīng)用到CFETR水冷陶瓷增殖堆柱狀中子學(xué)模型的計(jì)算中,對(duì)各種GVR方法的機(jī)理和性能進(jìn)行比較研究,并對(duì)不同GVR方法的FOMG、σ和Scoring進(jìn)行了分析。模型的幾何視圖如圖2所示,其主要參數(shù)見表1。

    圖2 柱狀CFETR中子學(xué)模型Fig.2 Cylinder CFETR neutronics model.

    表2是基于中子通量的GVR方法的迭代過程。第一次計(jì)算的能量截?cái)嘣O(shè)置為ECUT=10-3MeV,而后隨迭代的次數(shù)遞減,遞減到初始的能量截?cái)嘀禐橹?10-9MeV)。隨著GVR方法迭代次數(shù)的上升,平均誤差、σ總體下降,Scoring和FOMG總體上升。數(shù)據(jù)表明每次的迭代都改善了中子的空間分布情況,增大了MCNP有效計(jì)數(shù)率,誤差在空間的分布更加均勻,提高了全局的平均計(jì)算精度,而增加了全局減方差的計(jì)算效率。表2中出現(xiàn)了小部分?jǐn)?shù)據(jù)沒有按趨勢變化的情況,造成數(shù)據(jù)離散的原因可能是因?yàn)橛?jì)算時(shí)間過短(0.5 min),模擬的中子數(shù)太少,而造成數(shù)據(jù)的振蕩。

    表1 中子學(xué)模型主要參數(shù)Table 1 Main parameters for neutronics model.

    表2 基于中子通量的GVR方法的迭代過程(CPU Time=0.5 min)Table 2 Iterations of GVR method based on neutron flux(CPU Time=0.5 min).

    根據(jù)中子的相關(guān)信息分布可以在一定程度上反應(yīng)出系統(tǒng)內(nèi)的重要性分布,參考式(2)的表達(dá)形式,利用不同的中子信息去替換通量可以獲得其他形式的GVR的權(quán)窗。本文應(yīng)用了6種形式的GVR權(quán)窗,并在Davis的基礎(chǔ)上[1],也嘗試使用了其他的GVR形式(能量和中子徑跡)。由于中子的能量在輸運(yùn)的過程是一個(gè)非常重要的參數(shù),且在中子深穿透的過程中其能量也相應(yīng)降低,利用該規(guī)律也可以將中子的能量作為系統(tǒng)內(nèi)的重要性分布的信息來源,MC中子的徑跡也遵循類似的原理。

    表3比較了不同GVR方法的主要參數(shù)。由于計(jì)算模型相對(duì)規(guī)整,表3中的數(shù)據(jù)均是一次迭代的結(jié)果,使用相同的截?cái)鄺l件。對(duì)于沒有采取任何GVR方法的模擬(AN),正如所預(yù)期,AN具有最大的平均誤差、σ和Scoring和最小的FOMG值。而采用多種GVR方法后,相比與AN,F(xiàn)OMG都有很大的增長,相對(duì)誤差在空間的分布也更加平緩,且具有更高的計(jì)數(shù)率。通過對(duì)比多種形式的GVR方法可得,基于中子通量的全局減方差方法(GVR-F)具有最好的全局減方差效果。GVR-F與AN比較,其計(jì)算效率提高了6.43倍,平均誤差降低了1.79倍,σ減少了1.69倍且全局空間都有計(jì)數(shù)率。基于中子權(quán)重的全局減方差方法(GVR-W)也具有很好的全局減方差效率,GVR-W與AN比較,其計(jì)算效率提高了5.98倍,平均誤差降低了1.77倍。除此以外,本文所提出的基于中子能量的全局減方差方法(GVR-E)也是種可選的GVR應(yīng)用形式,其與AN比較,計(jì)算效率提高了5.11倍。GVR-E較Davis[1]所提出的基于中子數(shù)目的全局減方差方法(GVR-P)和 van Wijk[2]的基于中子誤差的全局減方差方法(GVR-R),其全局減方差效率都所提高。

    表3 不同GVR方法的主要參數(shù)(CPU Time=0.5 min)Table 3 Main parameters of different GVR methods(CPU Time=0.5 min).

    圖3展示了AN和GVR-F兩種方法在CFETR外包層和內(nèi)包層的中子通量,計(jì)算誤差和徑跡數(shù)的分布情況。對(duì)于內(nèi)外包層的中子通量分布,AN的中子通量只在模型部分區(qū)域內(nèi)有探測,其MCNP的有效計(jì)數(shù)率較低。在較厚區(qū)域,AN通量的相對(duì)誤差值陡增,且迅速到達(dá)最大誤差值。在使用GVR方法后,通量在整個(gè)模型內(nèi)都有探測結(jié)果,且誤差在空間中的絕大部分基本保持一個(gè)較低的值,變化梯度也較小,減少了全局的方差,進(jìn)而提高了全局的計(jì)算效率。相比于GVR方法,AN的中子徑跡衰減很快,在模型的后半段近乎沒有中子徑跡數(shù),該處的計(jì)數(shù)器的工作效率極低。而對(duì)于采用了GVR方法的中子徑跡分布,其在整個(gè)模型空間中都有較為可觀的徑跡分布,從而使得計(jì)數(shù)器具有較多的信息樣本,其探測結(jié)果也更可信。

    圖3 CFETR包層的中子信息分布Fig.3 Neutron distribution information of CFETR blanket.

    4 結(jié)語

    在解決具有全局特性的MC輸運(yùn)的精細(xì)計(jì)算問題中,本文介紹了多種GVR方法及其應(yīng)用。通過GVR權(quán)窗在CFETR的應(yīng)用并與AN比較可得,GVR方法的FOMG有很大的增長,通量的平均相對(duì)誤差較小,且具有更高的計(jì)數(shù)率。對(duì)于不同的GVR方法比較可得,基于中子通量的全局減方差權(quán)窗(GVR-F)具有最好的全局減方差效果。GVR-F與AN比較,其計(jì)算效率提高了6.43倍。本文所提出的基于中子能量的全局減方差方法(GVR-E)也是種可行的GVR應(yīng)用形式,其與AN比較,計(jì)算效率提高了5.11倍。

    除了本文中已使用的GVR形式以外,還可以考慮使用其他的中子分布信息。此外,將GVR方法所獲得的重要性分布應(yīng)用到IMP卡中來實(shí)現(xiàn)全局減方差的目的也有種可行的途徑。今后也有待將多種GVR方法應(yīng)用到更加復(fù)雜的三維CFETR模型中進(jìn)行對(duì)比和分析。

    參考文獻(xiàn)

    1Davis A,Andrew T.Comparison of global variance reduction techniques for Monte Carlo radiation transport simulations of ITER[J].Fusion Engineering and Design,2011,86:2698-2700.DOI:10.1016/j.fusengdes.2011.01.059

    2van Wijk A J,van den Eynde G,Hoogenboom J E.An easy to implement global variance reduction procedure for MCNP[J].Annals of Nuclear Energy,2011,38:2496-2503.DOI:10.1016/j.anucene.2011.07.037

    3Naish J,Fox F,Ghani Z,et al.Radiation mapping at JET and ITER using advanced computational acceleration techniques and tools[J].Nuclear Technology,2015,192:299-307.DOI:10.13182/nt14-132

    4Cooper M A,Larsen E W.Automated weight windows for global Monte Carlo particle transport calculations[J].Nuclear Science and Engineering,2001,137:1-13.DOI:10.13182/nse00-34

    5Booth,Thomas E.MCNP variance reduction examples[R].LA-UR-12 25907(2012),Los Alamos National Laboratory,2004.DOI:10.2172/1054246

    6李長楷,湯曉斌,岳愛忠.深貫穿問題脈沖高度計(jì)數(shù)減方差技巧[J].核技術(shù),2015,38(3):030501.DOI:10.11889/j.0253-3219.2015.hjs.38.030501 LI Changkai,TANG Xiaobin,YUE Aizhong.Pulse-height tally variance reduction in deep penetration problem[J].Nuclear Techniques,2015,38(3):030501.DOI:10.11889/j.0253-3219.2015.hjs.38.030501

    7許海燕,黃正豐,蔡少輝.蒙特卡羅中子輸運(yùn)問題的全局降方差方法[J].計(jì)算物理,2010,27(5):722-732.DOI:10.3969/j.issn.1001-246X.2010.05.014 XU Haiyan,HUANG Zhengfeng,CAI Shaohui.Global variance reduction method for Monte Carlo particle transport problems[J].Journal of Computational Physics,2010,27(5):722-732.DOI:10.3969/j.issn.1001-246X.2010.05.014

    收稿日期:2015-10-22,修回日期:2015-12-24

    Corresponding author:LI Jia,E-mail:lijia@ustc.edu.cn

    通信作者:李佳,E-mail:lijia@ustc.edu.cn

    DOI:10.11889/j.0253-3219.2016.hjs.39.030501

    中圖分類號(hào)TL62

    猜你喜歡
    蒙特卡羅
    基于復(fù)雜接續(xù)面的蒙特卡羅接續(xù)計(jì)算方法研究
    宮頸癌調(diào)強(qiáng)計(jì)劃在水與介質(zhì)中蒙特卡羅計(jì)算的劑量差異
    利用蒙特卡羅方法求解二重積分
    基于蒙特卡羅法的數(shù)字化電能表計(jì)量誤差評(píng)估
    基于蒙特卡羅模擬的過分相過電壓分析與抑制
    電氣化鐵道(2017年1期)2017-04-16 06:00:20
    脈沖堆弱源點(diǎn)火實(shí)驗(yàn)的蒙特卡羅數(shù)值模擬
    蒙特卡羅有效增殖因子計(jì)算方法研究
    基于蒙特卡羅的戰(zhàn)略投送能力動(dòng)態(tài)評(píng)估方法
    探討蒙特卡羅方法在解微分方程邊值問題中的應(yīng)用
    蒙特卡羅與響應(yīng)面法相結(jié)合的圓柱度公差模型求解
    av国产久精品久网站免费入址| 久久av网站| 国产在视频线精品| 人人妻人人添人人爽欧美一区卜| 一级毛片电影观看| 丰满乱子伦码专区| 免费久久久久久久精品成人欧美视频| 日韩制服骚丝袜av| 观看av在线不卡| 男人操女人黄网站| 欧美日韩成人在线一区二区| 日本-黄色视频高清免费观看| 女性生殖器流出的白浆| 国产伦理片在线播放av一区| 国产精品国产av在线观看| 亚洲精品一区蜜桃| 男人舔女人的私密视频| 亚洲精品中文字幕在线视频| 水蜜桃什么品种好| 五月天丁香电影| 国产片特级美女逼逼视频| 嫩草影院入口| av在线播放精品| 免费观看无遮挡的男女| 久久久久久久久免费视频了| 欧美精品av麻豆av| 2018国产大陆天天弄谢| 欧美精品国产亚洲| av网站在线播放免费| 最新中文字幕久久久久| 五月伊人婷婷丁香| 看非洲黑人一级黄片| 美女主播在线视频| 欧美亚洲 丝袜 人妻 在线| 国产精品偷伦视频观看了| 91aial.com中文字幕在线观看| 午夜日本视频在线| 亚洲国产精品一区二区三区在线| 欧美xxⅹ黑人| 国产熟女欧美一区二区| videosex国产| 在线看a的网站| 亚洲一码二码三码区别大吗| 国产白丝娇喘喷水9色精品| 精品少妇黑人巨大在线播放| 中文字幕人妻丝袜一区二区 | 老司机亚洲免费影院| a 毛片基地| 欧美日韩精品网址| 亚洲三区欧美一区| 777久久人妻少妇嫩草av网站| 五月天丁香电影| 黄色视频在线播放观看不卡| 少妇人妻精品综合一区二区| 少妇人妻精品综合一区二区| 亚洲精品成人av观看孕妇| 日韩一区二区三区影片| 久久久精品区二区三区| 韩国高清视频一区二区三区| 成人二区视频| www.自偷自拍.com| 精品人妻一区二区三区麻豆| av在线app专区| 国产一区有黄有色的免费视频| 日韩欧美一区视频在线观看| 欧美变态另类bdsm刘玥| 久久久精品免费免费高清| 亚洲精品一二三| 宅男免费午夜| 亚洲精品美女久久av网站| 色94色欧美一区二区| 女性生殖器流出的白浆| 人体艺术视频欧美日本| 免费女性裸体啪啪无遮挡网站| 午夜福利网站1000一区二区三区| 丝袜人妻中文字幕| 欧美日韩视频精品一区| 精品国产乱码久久久久久小说| 日日撸夜夜添| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品 欧美亚洲| 久久人人爽av亚洲精品天堂| 毛片一级片免费看久久久久| 春色校园在线视频观看| 日本欧美国产在线视频| 国产片特级美女逼逼视频| 国产男女超爽视频在线观看| 色94色欧美一区二区| 国产av精品麻豆| 欧美精品人与动牲交sv欧美| av在线观看视频网站免费| 久久久久精品久久久久真实原创| 中国三级夫妇交换| 国产亚洲欧美精品永久| 国产亚洲欧美精品永久| 在线观看国产h片| 毛片一级片免费看久久久久| 国产精品熟女久久久久浪| 一区二区三区激情视频| 久热久热在线精品观看| 亚洲少妇的诱惑av| 国产免费又黄又爽又色| 欧美+日韩+精品| 色网站视频免费| 国产福利在线免费观看视频| a级毛片在线看网站| 丝袜人妻中文字幕| 国产熟女午夜一区二区三区| 国产成人免费观看mmmm| 国产精品三级大全| 久久99热这里只频精品6学生| 一区二区三区乱码不卡18| 欧美日韩精品网址| 久久精品久久久久久久性| 国产成人aa在线观看| 999精品在线视频| 国产在线一区二区三区精| 日韩欧美一区视频在线观看| 精品人妻熟女毛片av久久网站| 99精国产麻豆久久婷婷| 纵有疾风起免费观看全集完整版| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区四区五区乱码 | 欧美日韩国产mv在线观看视频| 九色亚洲精品在线播放| 久热久热在线精品观看| 午夜福利视频在线观看免费| 久久久精品免费免费高清| 精品人妻熟女毛片av久久网站| 精品酒店卫生间| 亚洲精品aⅴ在线观看| 日韩不卡一区二区三区视频在线| 日韩大片免费观看网站| 男女边吃奶边做爰视频| 日韩欧美精品免费久久| 亚洲精品成人av观看孕妇| 三上悠亚av全集在线观看| 亚洲欧美色中文字幕在线| 国产成人一区二区在线| 日本免费在线观看一区| 男女边摸边吃奶| 国产精品秋霞免费鲁丝片| 色哟哟·www| 一二三四中文在线观看免费高清| 午夜激情av网站| 精品亚洲乱码少妇综合久久| 极品少妇高潮喷水抽搐| 久久这里有精品视频免费| www.av在线官网国产| 伊人亚洲综合成人网| 久久精品国产自在天天线| 国产男女内射视频| 亚洲 欧美一区二区三区| 久久国内精品自在自线图片| 性色av一级| 丝袜喷水一区| 高清在线视频一区二区三区| 亚洲精品国产av成人精品| 麻豆精品久久久久久蜜桃| 中文乱码字字幕精品一区二区三区| 韩国av在线不卡| 最近手机中文字幕大全| 天堂俺去俺来也www色官网| 午夜影院在线不卡| 久久久久精品人妻al黑| 日韩人妻精品一区2区三区| 色婷婷av一区二区三区视频| 丰满迷人的少妇在线观看| 老司机亚洲免费影院| 熟妇人妻不卡中文字幕| 免费看不卡的av| 久久99一区二区三区| 亚洲国产av影院在线观看| 国产日韩一区二区三区精品不卡| 色播在线永久视频| 青春草国产在线视频| 美女午夜性视频免费| 亚洲av电影在线进入| 日本av免费视频播放| 免费av中文字幕在线| av电影中文网址| 婷婷色av中文字幕| 午夜免费观看性视频| 国产av码专区亚洲av| 国产 一区精品| 搡老乐熟女国产| 中文字幕精品免费在线观看视频| 亚洲国产欧美日韩在线播放| 日韩一区二区视频免费看| 国产有黄有色有爽视频| 久久久久精品久久久久真实原创| 久久影院123| 黄色视频在线播放观看不卡| 久久韩国三级中文字幕| 国产老妇伦熟女老妇高清| 日韩制服丝袜自拍偷拍| 午夜福利,免费看| 精品人妻偷拍中文字幕| 爱豆传媒免费全集在线观看| 国产1区2区3区精品| 丰满饥渴人妻一区二区三| 亚洲精品久久成人aⅴ小说| 午夜福利网站1000一区二区三区| 亚洲av欧美aⅴ国产| 尾随美女入室| 麻豆av在线久日| 亚洲综合色网址| 国产av国产精品国产| 99热全是精品| 久久毛片免费看一区二区三区| 大话2 男鬼变身卡| 久久精品国产自在天天线| 青青草视频在线视频观看| 亚洲欧美日韩另类电影网站| 欧美精品高潮呻吟av久久| 岛国毛片在线播放| 亚洲精品aⅴ在线观看| 精品99又大又爽又粗少妇毛片| 一区二区日韩欧美中文字幕| 亚洲欧美精品自产自拍| 又大又黄又爽视频免费| 男女国产视频网站| 日本欧美国产在线视频| 三上悠亚av全集在线观看| 99re6热这里在线精品视频| 日韩av在线免费看完整版不卡| 99久久中文字幕三级久久日本| 久久青草综合色| 亚洲三区欧美一区| 一区福利在线观看| 最近最新中文字幕大全免费视频 | 不卡av一区二区三区| 免费黄网站久久成人精品| 18禁国产床啪视频网站| 国产免费福利视频在线观看| 久久精品亚洲av国产电影网| 国产男女内射视频| 国产精品熟女久久久久浪| 国产精品久久久av美女十八| 久久久国产精品麻豆| 只有这里有精品99| 成人18禁高潮啪啪吃奶动态图| av女优亚洲男人天堂| 国产xxxxx性猛交| 国产av码专区亚洲av| 精品国产一区二区久久| 卡戴珊不雅视频在线播放| 欧美日韩亚洲高清精品| 免费在线观看完整版高清| 91午夜精品亚洲一区二区三区| 大码成人一级视频| 欧美日本中文国产一区发布| 老鸭窝网址在线观看| 欧美精品av麻豆av| 建设人人有责人人尽责人人享有的| 人妻 亚洲 视频| 亚洲精品aⅴ在线观看| 日韩大片免费观看网站| 男男h啪啪无遮挡| 久久精品国产鲁丝片午夜精品| 亚洲久久久国产精品| 亚洲一区二区三区欧美精品| 精品一区在线观看国产| 麻豆乱淫一区二区| 久久精品aⅴ一区二区三区四区 | 日韩伦理黄色片| 日日摸夜夜添夜夜爱| 婷婷成人精品国产| 亚洲久久久国产精品| 91aial.com中文字幕在线观看| 欧美国产精品一级二级三级| 久久久久网色| 日韩电影二区| 母亲3免费完整高清在线观看 | 免费在线观看视频国产中文字幕亚洲 | 国产精品 欧美亚洲| 亚洲五月色婷婷综合| 亚洲人成网站在线观看播放| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 精品国产超薄肉色丝袜足j| 亚洲精华国产精华液的使用体验| 精品福利永久在线观看| 多毛熟女@视频| 在线免费观看不下载黄p国产| 亚洲,欧美精品.| 午夜日本视频在线| 午夜激情久久久久久久| 不卡av一区二区三区| 乱人伦中国视频| 99国产综合亚洲精品| 女人精品久久久久毛片| 欧美人与性动交α欧美软件| 丝袜喷水一区| 叶爱在线成人免费视频播放| 永久网站在线| 亚洲av福利一区| 亚洲婷婷狠狠爱综合网| 亚洲国产成人一精品久久久| 亚洲av电影在线进入| 久久久久久伊人网av| 在现免费观看毛片| 午夜福利一区二区在线看| av一本久久久久| 在线 av 中文字幕| 午夜免费观看性视频| 日产精品乱码卡一卡2卡三| 国产综合精华液| 国产在线一区二区三区精| 国产精品.久久久| 最新中文字幕久久久久| av卡一久久| 国产免费又黄又爽又色| 电影成人av| 搡老乐熟女国产| 国产精品免费视频内射| 成人午夜精彩视频在线观看| 亚洲精品,欧美精品| 国产av国产精品国产| 精品亚洲成国产av| 亚洲av日韩在线播放| 777米奇影视久久| 97精品久久久久久久久久精品| 黄片小视频在线播放| 日韩中字成人| 午夜91福利影院| 日韩电影二区| 两性夫妻黄色片| 色播在线永久视频| 亚洲国产毛片av蜜桃av| a级毛片在线看网站| 午夜免费鲁丝| 国产乱人偷精品视频| 亚洲欧美精品自产自拍| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| 亚洲在久久综合| 日本wwww免费看| 欧美国产精品va在线观看不卡| 欧美激情 高清一区二区三区| 美女主播在线视频| 美国免费a级毛片| 一区二区日韩欧美中文字幕| 欧美日韩视频精品一区| 大片免费播放器 马上看| 久久久久久久亚洲中文字幕| 纵有疾风起免费观看全集完整版| 久久久久人妻精品一区果冻| 男人舔女人的私密视频| 亚洲一码二码三码区别大吗| 国产高清国产精品国产三级| 97精品久久久久久久久久精品| 亚洲在久久综合| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 免费在线观看完整版高清| 一区二区av电影网| 婷婷色综合www| 国产不卡av网站在线观看| av福利片在线| 色播在线永久视频| 久久婷婷青草| 国产熟女午夜一区二区三区| 咕卡用的链子| 国产色婷婷99| 日韩av不卡免费在线播放| 欧美日韩视频高清一区二区三区二| 熟妇人妻不卡中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 国产男女内射视频| 国产精品久久久久久久久免| 欧美激情极品国产一区二区三区| 免费播放大片免费观看视频在线观看| 日本欧美国产在线视频| 18禁国产床啪视频网站| 日韩不卡一区二区三区视频在线| 五月天丁香电影| 久久av网站| 亚洲伊人久久精品综合| av女优亚洲男人天堂| 日韩,欧美,国产一区二区三区| 日本wwww免费看| 日本vs欧美在线观看视频| 中文乱码字字幕精品一区二区三区| 赤兔流量卡办理| 男人操女人黄网站| 久久精品国产综合久久久| 亚洲精品久久成人aⅴ小说| 欧美bdsm另类| 国产精品久久久久成人av| 亚洲欧美成人精品一区二区| 国产男女超爽视频在线观看| 精品久久蜜臀av无| 日韩中文字幕欧美一区二区 | 成人18禁高潮啪啪吃奶动态图| 精品亚洲成a人片在线观看| 国产成人精品无人区| av在线播放精品| 成年女人在线观看亚洲视频| 精品国产乱码久久久久久小说| 99久久精品国产国产毛片| 十八禁网站网址无遮挡| xxx大片免费视频| 久久人人爽人人片av| 亚洲成色77777| 国产视频首页在线观看| 高清在线视频一区二区三区| 成年美女黄网站色视频大全免费| 91久久精品国产一区二区三区| 亚洲精品国产av蜜桃| 日韩一本色道免费dvd| 大片免费播放器 马上看| 久久久a久久爽久久v久久| 欧美精品一区二区免费开放| 亚洲欧美一区二区三区久久| 两个人免费观看高清视频| 最新中文字幕久久久久| 欧美日韩av久久| 久久鲁丝午夜福利片| 又粗又硬又长又爽又黄的视频| 欧美亚洲日本最大视频资源| 亚洲精品aⅴ在线观看| 久久97久久精品| h视频一区二区三区| 亚洲国产av影院在线观看| 韩国av在线不卡| 亚洲欧美一区二区三区黑人 | 久久这里有精品视频免费| 婷婷色综合大香蕉| a 毛片基地| 亚洲在久久综合| 国产激情久久老熟女| 自线自在国产av| 美女大奶头黄色视频| av卡一久久| 久久久精品免费免费高清| 国产成人91sexporn| 精品少妇一区二区三区视频日本电影 | 国产成人精品婷婷| 中文字幕制服av| 交换朋友夫妻互换小说| 久久 成人 亚洲| videossex国产| 亚洲精华国产精华液的使用体验| 国产精品亚洲av一区麻豆 | 国产精品熟女久久久久浪| www.av在线官网国产| 久久久国产一区二区| 免费黄网站久久成人精品| 免费人妻精品一区二区三区视频| 最近手机中文字幕大全| 可以免费在线观看a视频的电影网站 | 天堂中文最新版在线下载| 亚洲精品av麻豆狂野| 久久毛片免费看一区二区三区| 日韩大片免费观看网站| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 亚洲欧美清纯卡通| 亚洲国产日韩一区二区| 国产精品欧美亚洲77777| 亚洲av在线观看美女高潮| 国产精品偷伦视频观看了| a级毛片黄视频| 午夜福利在线免费观看网站| 黑人欧美特级aaaaaa片| 久久久国产欧美日韩av| 99香蕉大伊视频| 下体分泌物呈黄色| 亚洲成色77777| 久久热在线av| 精品人妻在线不人妻| 啦啦啦在线观看免费高清www| 精品国产超薄肉色丝袜足j| 男的添女的下面高潮视频| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 欧美中文综合在线视频| 97在线人人人人妻| 哪个播放器可以免费观看大片| 侵犯人妻中文字幕一二三四区| 制服丝袜香蕉在线| 亚洲成人手机| 在线精品无人区一区二区三| 午夜激情久久久久久久| 狠狠精品人妻久久久久久综合| 一区二区三区精品91| 人妻少妇偷人精品九色| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利视频精品| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 国产精品国产av在线观看| 高清av免费在线| 女人高潮潮喷娇喘18禁视频| 又大又黄又爽视频免费| 999精品在线视频| 街头女战士在线观看网站| 亚洲国产精品成人久久小说| 国产高清不卡午夜福利| 亚洲色图综合在线观看| 亚洲,欧美,日韩| 久久热在线av| 国产精品不卡视频一区二区| 免费久久久久久久精品成人欧美视频| 天天躁狠狠躁夜夜躁狠狠躁| 菩萨蛮人人尽说江南好唐韦庄| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 妹子高潮喷水视频| 只有这里有精品99| 亚洲精品第二区| 久久免费观看电影| 国产成人aa在线观看| 中文天堂在线官网| 男人爽女人下面视频在线观看| 成年美女黄网站色视频大全免费| av国产久精品久网站免费入址| 日韩熟女老妇一区二区性免费视频| av国产精品久久久久影院| 亚洲精品美女久久久久99蜜臀 | 亚洲男人天堂网一区| 欧美少妇被猛烈插入视频| 亚洲三级黄色毛片| 熟女少妇亚洲综合色aaa.| 国产色婷婷99| 欧美精品国产亚洲| 日本欧美视频一区| 国产一区亚洲一区在线观看| videos熟女内射| 各种免费的搞黄视频| 国产亚洲av片在线观看秒播厂| 桃花免费在线播放| 老女人水多毛片| 一本色道久久久久久精品综合| 91成人精品电影| av免费在线看不卡| 国产野战对白在线观看| 国产精品亚洲av一区麻豆 | 波多野结衣av一区二区av| 女人精品久久久久毛片| 极品少妇高潮喷水抽搐| 成人国产麻豆网| 亚洲伊人色综图| 精品久久久久久电影网| 男人舔女人的私密视频| 看非洲黑人一级黄片| 纯流量卡能插随身wifi吗| 26uuu在线亚洲综合色| 成人二区视频| 观看av在线不卡| 亚洲成人一二三区av| 日本色播在线视频| 女人被躁到高潮嗷嗷叫费观| 免费播放大片免费观看视频在线观看| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久小说| 久久久久久久国产电影| 国产国语露脸激情在线看| 伊人久久大香线蕉亚洲五| 一区二区av电影网| 精品人妻一区二区三区麻豆| 最近最新中文字幕免费大全7| 午夜福利视频在线观看免费| 中文精品一卡2卡3卡4更新| 男女无遮挡免费网站观看| av网站在线播放免费| 久久这里有精品视频免费| 精品一区在线观看国产| 亚洲第一青青草原| 国产成人av激情在线播放| 欧美激情 高清一区二区三区| 欧美日韩视频高清一区二区三区二| 亚洲久久久国产精品| av国产精品久久久久影院| 日韩熟女老妇一区二区性免费视频| 91精品伊人久久大香线蕉| 叶爱在线成人免费视频播放| 日本欧美国产在线视频| 男的添女的下面高潮视频| 三级国产精品片| 亚洲天堂av无毛| 国产深夜福利视频在线观看| 国产精品秋霞免费鲁丝片| 老女人水多毛片| 午夜激情久久久久久久| 免费看不卡的av| 日本爱情动作片www.在线观看| 久久久国产精品麻豆| 国产在视频线精品| 久久精品国产亚洲av高清一级| 久久精品夜色国产| 欧美精品一区二区免费开放| 香蕉精品网在线| 久久久久久伊人网av| 美女午夜性视频免费| 精品一区在线观看国产| 伦精品一区二区三区| 国产片特级美女逼逼视频| 久久久久精品人妻al黑| av卡一久久| 久久99一区二区三区| 韩国av在线不卡| 国产成人a∨麻豆精品| 欧美日韩综合久久久久久| www日本在线高清视频| 久久久亚洲精品成人影院| 国产无遮挡羞羞视频在线观看| 日韩不卡一区二区三区视频在线| 午夜福利视频在线观看免费| 97在线视频观看| 亚洲精品美女久久久久99蜜臀 | 美女高潮到喷水免费观看| 日韩欧美精品免费久久| 人成视频在线观看免费观看|