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

    基于神經(jīng)網(wǎng)絡(luò)的CFD粗網(wǎng)格模擬優(yōu)化方法研究

    2021-06-16 12:14:42劉曉晶
    核技術(shù) 2021年6期
    關(guān)鍵詞:方根機(jī)器數(shù)值

    靳 爽 劉曉晶 程 旭

    (上海交通大學(xué) 核科學(xué)與工程學(xué)院 上海200240)

    在反應(yīng)堆計(jì)算流體力學(xué)(Computational Fluid Dynamics,CFD)數(shù)值模擬中,計(jì)算精度和計(jì)算時(shí)間是一對突出的矛盾。為了能夠較好地兼顧計(jì)算效率和計(jì)算精度,實(shí)現(xiàn)在較高計(jì)算效率下取得滿意的計(jì)算結(jié)果,研究人員嘗試了多種新技術(shù)方法,比如并行計(jì)算、對網(wǎng)格劃分方式進(jìn)行優(yōu)化[1-2]等。隨著大數(shù)據(jù)和人工智能等領(lǐng)域的發(fā)展,機(jī)器學(xué)習(xí)也成為解決上述問題的一種可行路徑。一方面,大量先前已有的數(shù)值模擬結(jié)果、理論研究成果和實(shí)驗(yàn)結(jié)果積累為數(shù)據(jù)訓(xùn)練提供了廣泛的訓(xùn)練資源;另一方面,隨著反應(yīng)堆熱工水力研究和應(yīng)用的不斷發(fā)展,對“物理過程相似、工況參數(shù)不同的大批量多組工況”的CFD數(shù)值模擬需求日益凸顯,例如研究多種不同入口速度下的棒束流動傳熱CFD結(jié)果、多種搖擺條件(搖擺俯角、搖擺周期的各類組合)下的反應(yīng)堆CFD計(jì)算結(jié)果等,在該類問題中,基于神經(jīng)網(wǎng)絡(luò)的CFD粗網(wǎng)格模擬優(yōu)化方法對于“實(shí)現(xiàn)在較高計(jì)算效率下實(shí)現(xiàn)較高精度的數(shù)值模擬”具有顯著效果。

    根據(jù)美國愛達(dá)荷國家實(shí)驗(yàn)室于2018年發(fā)布的報(bào)告[3],在數(shù)據(jù)驅(qū)動的機(jī)器學(xué)習(xí)架構(gòu)中可以使低分辨率的模型計(jì)算出高分辨率模型下的模擬結(jié)果,并對絕熱腔室中射入空氣湍流模擬結(jié)果進(jìn)行了優(yōu)化。目前,通過機(jī)器學(xué)習(xí)針對反應(yīng)堆熱工水力數(shù)值模擬進(jìn)行優(yōu)化的案例研究還較少,本文擬應(yīng)用機(jī)器學(xué)習(xí)方法,針對粗網(wǎng)格下5×5棒束出口截面冷卻劑子通道溫度的CFD結(jié)果進(jìn)行優(yōu)化,以期為機(jī)器學(xué)習(xí)方法在反應(yīng)堆熱工領(lǐng)域?qū)崿F(xiàn)較高計(jì)算效率和較高計(jì)算精度的數(shù)值模擬應(yīng)用提供參考。

    1 數(shù)值模擬優(yōu)化方法

    1.1 棒束模型

    本文計(jì)算的棒束模型是截短型組件的5×5棒束柵元,其主要參數(shù)參考一體化模塊小型反應(yīng)堆常用的截短型組件參數(shù)而設(shè)定[4-7],如表1所示。

    針對本文計(jì)算的棒束模型,可以劃分為36個(gè)子通道,從左上到右下依次編號為子通道1(sub 1)、子通道2(sub 2)、……、子通道36(sub 36),如圖1所示。

    表1 5×5棒束柵元主要參數(shù)Table 1 Main parameters of 5×5 rod bundle

    圖1 5×5棒束柵元子通道劃分方式Fig.1 Subchannel division method of 5×5 rod bundle

    1.2 問題描述

    本文旨在對不同入口速度工況中出口截面上36個(gè)子通道各自的平均溫度值進(jìn)行修正,粗網(wǎng)格和細(xì)網(wǎng)格均采用STAR-CCM中的結(jié)構(gòu)化定向網(wǎng)格[8],基礎(chǔ)網(wǎng)格為多邊形網(wǎng)格,對于邊界層處的網(wǎng)格,均使用棱柱層網(wǎng)格生成,粗網(wǎng)格和細(xì)網(wǎng)格下均設(shè)置棱柱層層數(shù)為三層,棱柱層增長率取默認(rèn)值1.5,棱柱層總厚度取默認(rèn)值即網(wǎng)格尺寸基數(shù)的33.3%。其網(wǎng)格參數(shù)的具體設(shè)置如表2所示。

    細(xì)網(wǎng)格參數(shù)的選取達(dá)到網(wǎng)格無關(guān)性的要求,基礎(chǔ)尺寸為0.5 mm、0.6 mm、0.7 mm網(wǎng)格下的子通道1、子通道2、子通道16(分別為角通道、邊通道、中心通道)的溫度計(jì)算結(jié)果如表3所示。

    本文根據(jù)入口速度的不同,將20組工況作為訓(xùn)練組用于機(jī)器學(xué)習(xí),另外5組工況作為驗(yàn)證組以檢驗(yàn)修正的效果。其中,5組驗(yàn)證組工況中既包括相對于訓(xùn)練組的內(nèi)含參數(shù)工況,也包括外延參數(shù)工況。訓(xùn)練組工況及驗(yàn)證組工況的入口速度設(shè)定情況如表4所示。鑒于本文計(jì)算的物理過程并不復(fù)雜,在本文涉及的計(jì)算工況中,所選用的湍流模型均為可實(shí)現(xiàn)的k-ε(Realizablek-ε)模型。

    表2 粗網(wǎng)格參數(shù)和細(xì)網(wǎng)格參數(shù)Table 2 Parameters of coarse grid and fine grid

    表3 網(wǎng)格無關(guān)性驗(yàn)證結(jié)果Table 3 Results of grid independence verification

    表4 訓(xùn)練組和驗(yàn)證組工況入口速度設(shè)定情況Table 4 Inlet speed of training group and verification group working conditions

    1.3 基于機(jī)器學(xué)習(xí)的優(yōu)化方法

    數(shù)值模擬的誤差來源主要包括以下三類:第一類是由于物理簡化和(或)數(shù)學(xué)近似導(dǎo)致的模型誤差;第二類是由于網(wǎng)格的存在,導(dǎo)致守恒方程或源項(xiàng)在時(shí)空平均方法下產(chǎn)生信息丟失所引起的網(wǎng)格誤差;第三類是由于迭代收斂、算法選擇等原因引起的其他誤差。

    對于第一類誤差,針對流體和壁面間傳熱、質(zhì)量交換、動量交換等現(xiàn)象存在多種關(guān)系模型進(jìn)行模擬,與網(wǎng)格尺寸相關(guān)的特征長度是這些模型中的一個(gè)重要參數(shù)。對于第二類誤差,同樣由網(wǎng)格尺寸決定。而第三類誤差相較前兩類誤差對計(jì)算結(jié)果的影響要小得多。因此,數(shù)值模擬的誤差與網(wǎng)格尺寸密切相關(guān),且實(shí)際分析時(shí)往往無法區(qū)分模型誤差和網(wǎng)格誤差,應(yīng)將其作為一個(gè)整體來考慮。

    我們可以運(yùn)用模式識別和統(tǒng)計(jì)分析等機(jī)器學(xué)習(xí)方法,從數(shù)據(jù)中直接獲得容易被人忽視的信息,試圖尋找到包括網(wǎng)格信息、模型信息、物理信息在內(nèi)的參數(shù)與誤差之間的定量關(guān)系。在此基礎(chǔ)上,運(yùn)用該定量關(guān)系對粗網(wǎng)格模擬結(jié)果進(jìn)行修正,實(shí)現(xiàn)滿足計(jì)算要求的快速數(shù)值模擬。

    根據(jù)問題的需要,我們可以將模擬對象的幾何參數(shù)、模擬工況的物理參數(shù)及數(shù)值計(jì)算的網(wǎng)格參數(shù)納入機(jī)器學(xué)習(xí)的輸入變量中。針對本文研究的工況,每組工況均采用同一套粗網(wǎng)格和同一套細(xì)網(wǎng)格,各組工況在進(jìn)行CFD計(jì)算時(shí)選用的湍流模型、傳熱模型一致,僅有冷卻劑入口速度設(shè)定不同。因此,本文通過機(jī)器學(xué)習(xí)需要得出的是入口速度、粗網(wǎng)格下出口截面各個(gè)子通道溫度、細(xì)網(wǎng)格下出口截面各個(gè)子通道溫度之間的隱含關(guān)系。

    根據(jù)Chang等[9]的文章,目前數(shù)據(jù)驅(qū)動的熱工水力模擬機(jī)器學(xué)習(xí)架構(gòu)可以依據(jù)數(shù)值模擬是否包含偏微分方程、偏微分方程的形式是否已知、偏微分方程本身是否作為機(jī)器學(xué)習(xí)對象、模型的建立是否需要在不同尺度上分別進(jìn)行等四項(xiàng)標(biāo)準(zhǔn)分成5種類型,本文的研究屬于上述五種架構(gòu)類型中的第二類,其主要步驟如圖2所示。

    圖2 基于機(jī)器學(xué)習(xí)的數(shù)值模擬優(yōu)化框架[9]Fig.2 Optimization framework for numerical simulation based on machine learning[9]

    本文中,定義粗網(wǎng)格下子通道i(subi)的平均溫度為TLi,作為待修正值;細(xì)網(wǎng)格下子通道i的溫度記為THi,作為精確值(修正目標(biāo)值)。定義誤差函數(shù)如下:

    根據(jù)前述的優(yōu)化原理,將誤差函數(shù)E i設(shè)置為如下的形式:

    式中:V為冷卻劑入口速度,子通道編號i=1,2,3,...,36,誤差函數(shù)中的具體函數(shù)形式由機(jī)器學(xué)習(xí)方法得到。將粗網(wǎng)格計(jì)算結(jié)果按照誤差函數(shù)修正后的子通道i的溫度記為TMi,即有:

    在機(jī)器學(xué)習(xí)方法的選擇上,目前主要使用深度學(xué)習(xí)的方法。原則上,任何兩層以上的神經(jīng)網(wǎng)絡(luò)都可以被認(rèn)為是深度學(xué)習(xí)。Hornik等[10]在文章中論證了多層神經(jīng)網(wǎng)絡(luò)是一種通用方法,可以捕捉任何可度量信息的特性。神經(jīng)網(wǎng)絡(luò)包括卷積神經(jīng)網(wǎng)絡(luò)和前饋神經(jīng)網(wǎng)絡(luò)兩種方式。其中,卷積神經(jīng)網(wǎng)絡(luò)通常具有更好的預(yù)測能力,但所需的數(shù)據(jù)訓(xùn)練量也更多。鑒于本文研究的工況相對簡單,本文選用所需數(shù)據(jù)量較少的前饋神經(jīng)網(wǎng)絡(luò)。對于神經(jīng)網(wǎng)絡(luò)層數(shù)和神經(jīng)元數(shù)量的選取,海軍工程大學(xué)曹植珺等[11]的綜述文章中指出目前缺少很好的參數(shù)調(diào)節(jié)的手段,包括訓(xùn)練神經(jīng)網(wǎng)絡(luò)時(shí)隱含層數(shù)、隱含節(jié)點(diǎn)數(shù)、輸入節(jié)點(diǎn)數(shù)等基本依靠經(jīng)驗(yàn)。對于不復(fù)雜的關(guān)系結(jié)構(gòu),神經(jīng)網(wǎng)絡(luò)的層數(shù)通常為2~3層,神經(jīng)元數(shù)量不超過10個(gè),本文選擇的神經(jīng)網(wǎng)絡(luò)層數(shù)為兩層,神經(jīng)元數(shù)量為4個(gè),其結(jié)構(gòu)如圖3所示。

    第一層為輸入變量層,包括工況冷卻劑入口速度V以及粗網(wǎng)格下出口截面子通道溫度TLi(i=1,2,3,...,36)共計(jì)37個(gè)輸入變量;第二層為神經(jīng)元層,共包括4個(gè)神經(jīng)元N1,…,N4;第三層為輸出變量層,包括粗網(wǎng)格與細(xì)網(wǎng)格下出口截面各個(gè)子通道溫度的誤差預(yù)估值E1,E2,…,E36,共計(jì)36個(gè)輸出變量。根據(jù)神經(jīng)網(wǎng)絡(luò)的基本原理,各層元素之間的關(guān)系如下:

    圖3 前饋神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)架構(gòu)Fig.3 Framework of feedforward neural network(FNN)

    式中:f為神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)的函數(shù)形式,為表達(dá)方便,記:

    本文采用均方根誤差來評估優(yōu)化效果,分別定義未經(jīng)優(yōu)化的粗網(wǎng)格下出口截面36個(gè)子通道溫度與精確值(細(xì)網(wǎng)格下相應(yīng)計(jì)算結(jié)果)之間的均方根誤差為RMSEL,經(jīng)過優(yōu)化的粗網(wǎng)格下出口截面36個(gè)子通道溫度與精確值(細(xì)網(wǎng)格計(jì)算結(jié)果)之間的均方根誤差為RMSEM,其形式如下:

    通過比較RMSEM與RMSEL的大小差異來評估優(yōu)化效果,如果0~RMSEM?RMSEL,則可認(rèn)為達(dá)到預(yù)期優(yōu)化效果。

    2 優(yōu)化結(jié)果

    2.1 計(jì)算精度

    使用機(jī)器學(xué)習(xí)得到的誤差函數(shù)對訓(xùn)練組和驗(yàn)證組的粗網(wǎng)格結(jié)果進(jìn)行優(yōu)化前后,出口截面36個(gè)子通道溫度與細(xì)網(wǎng)格下模擬結(jié)果的均方根誤差如圖4所示。

    其中,5組驗(yàn)證組工況的粗網(wǎng)格結(jié)果優(yōu)化前后的均方根誤差如表5所示。

    圖4 粗網(wǎng)格結(jié)果優(yōu)化前后的均方根誤差Fig.4 Root mean square error(RMSE)before and after optimization

    表5 驗(yàn)證組工況粗網(wǎng)格結(jié)果優(yōu)化前后的均方根誤差Table 5 Root mean square error before and after optimization of the verification group results

    通過圖4和表5可以看出,訓(xùn)練組和驗(yàn)證組粗網(wǎng)格下修正優(yōu)化后的結(jié)果相比于優(yōu)化前,與細(xì)網(wǎng)格下結(jié)果的均方根誤差大大降低,5組檢驗(yàn)組工況出口截面36個(gè)子通道溫度均方根誤差平均值從1.16 K降低為9.24×10-3K,滿足0~RMSEM?RMSEL,達(dá)到預(yù)期的優(yōu)化效果。

    以入口速度為2.25 m?s-1的工況為例,其細(xì)網(wǎng)格下、粗網(wǎng)格下結(jié)果修正前后36個(gè)子通道的溫度如圖5所示。通過圖5可以看出,修正優(yōu)化后的各個(gè)子通道溫度值相比于優(yōu)化前,明顯更加接近細(xì)網(wǎng)格下的結(jié)果;且對于邊通道和角通道,優(yōu)化提升的效果更加明顯。

    2.2 計(jì)算效率

    使用本文的方法可以對粗網(wǎng)格下的計(jì)算結(jié)果進(jìn)行修正,取得了良好的優(yōu)化效果,這也使得CFD計(jì)算效率大幅提升。

    對于粗網(wǎng)格、細(xì)網(wǎng)格工況,在運(yùn)算時(shí)均以連續(xù)20個(gè)迭代步內(nèi)出口平均溫度變化不超過0.01℃為判斷計(jì)算達(dá)到收斂的條件。本文涉及的表4中共計(jì)25個(gè)工況在粗網(wǎng)格、細(xì)網(wǎng)格下進(jìn)行CFD運(yùn)算所需要的總CPU系統(tǒng)時(shí)間如表6所示。

    圖5 V=2.25 m?s-1工況下粗網(wǎng)格結(jié)果優(yōu)化前后子通道溫度值與細(xì)網(wǎng)格下結(jié)果的對比Fig.5 Comparison of subchannel temperature before and after optimization under the condition of V=2.25 m?s-1

    從表6可以看出,不同工況下粗網(wǎng)格計(jì)算所需的總CPU系統(tǒng)時(shí)間遠(yuǎn)遠(yuǎn)小于細(xì)網(wǎng)格下所需時(shí)間,上述25個(gè)工況粗網(wǎng)格下總CPU系統(tǒng)時(shí)間的平均值約為細(xì)網(wǎng)格下平均值的2.54%。就本文研究的問題而言,若需要得到入口速度從1.65~2.85 m·s-1間等距100個(gè)工況的計(jì)算結(jié)果,傳統(tǒng)方法需運(yùn)行100個(gè)細(xì)網(wǎng)格工況,而基于神經(jīng)網(wǎng)絡(luò)的粗網(wǎng)格優(yōu)化方法下,僅需計(jì)算25個(gè)細(xì)網(wǎng)格工況和100個(gè)粗網(wǎng)格工況(以表6的CPU系統(tǒng)時(shí)間平均值計(jì)算,約相當(dāng)于三個(gè)細(xì)網(wǎng)格工況的計(jì)算量),總體的計(jì)算效率得到大幅提高。

    3 結(jié)語

    本研究利用前饋神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)方法,對不同入口速度工況下出口截面子通道冷卻劑溫度粗網(wǎng)格下的CFD數(shù)值模擬結(jié)果進(jìn)行了優(yōu)化,得到如下結(jié)論:

    1)采用神經(jīng)網(wǎng)絡(luò)方法對本文工況粗網(wǎng)格下的結(jié)果進(jìn)行修正后,5組檢驗(yàn)組工況出口截面36個(gè)子通道溫度的均方根誤差平均值從1.16 K降低為9.24×10-3K,誤差大大減小,優(yōu)化效果顯著。

    2)在本文研究的問題中,粗網(wǎng)格下總CPU系統(tǒng)時(shí)間的平均值約為細(xì)網(wǎng)格下平均值的2.54%。通過神經(jīng)網(wǎng)絡(luò)對粗網(wǎng)格結(jié)果優(yōu)化的方式,對于批量化的CFD計(jì)算或者已有較為充分的訓(xùn)練數(shù)據(jù)的CFD計(jì)算,其計(jì)算效率將大幅提高。

    3)就本文算例來看,基于前饋神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)框架對于內(nèi)部預(yù)測(即驗(yàn)證工況的速度值在訓(xùn)練集速度值范圍內(nèi))的優(yōu)化效果好于外部預(yù)測。但由于本文算例較為局限,這一結(jié)論還需要更多算例結(jié)果的支持驗(yàn)證。

    本文的研究為利用機(jī)器學(xué)習(xí)框架對粗網(wǎng)格數(shù)值模擬進(jìn)行優(yōu)化,從而兼顧計(jì)算效率與計(jì)算精度提供了可行參照。未來,可以在反應(yīng)堆熱工水力領(lǐng)域的更多問題特別是多耦合或非穩(wěn)態(tài)問題上進(jìn)行更廣泛的驗(yàn)證。

    表6 25個(gè)工況在粗網(wǎng)格、細(xì)網(wǎng)格下運(yùn)行所需的總CPU系統(tǒng)時(shí)間Table 6 The total CPU system time required for 25 operating conditions to run under coarse and fine grids

    猜你喜歡
    方根機(jī)器數(shù)值
    用固定數(shù)值計(jì)算
    方根拓展探究
    機(jī)器狗
    機(jī)器狗
    數(shù)值大小比較“招招鮮”
    未來機(jī)器城
    電影(2018年8期)2018-09-21 08:00:06
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    揭開心算方根之謎
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    無敵機(jī)器蛛
    丰满少妇做爰视频| 免费观看av网站的网址| 免费人成在线观看视频色| 成年人午夜在线观看视频| 亚洲精品色激情综合| 美女脱内裤让男人舔精品视频| 国产欧美日韩一区二区三区在线 | 欧美 日韩 精品 国产| av国产精品久久久久影院| 欧美三级亚洲精品| 亚洲国产精品国产精品| 99热全是精品| 国产精品一区www在线观看| 一区二区三区乱码不卡18| 91久久精品国产一区二区三区| 亚洲精品色激情综合| 日产精品乱码卡一卡2卡三| 大又大粗又爽又黄少妇毛片口| 最近手机中文字幕大全| av在线播放精品| 黑人欧美特级aaaaaa片| 亚洲国产av新网站| 我的老师免费观看完整版| 香蕉精品网在线| 精品国产国语对白av| 久久ye,这里只有精品| 大码成人一级视频| 亚洲精品乱久久久久久| 最新的欧美精品一区二区| 亚洲国产成人一精品久久久| 久久久精品94久久精品| 美女国产视频在线观看| 国产有黄有色有爽视频| 国产日韩欧美在线精品| 免费黄频网站在线观看国产| 亚洲精品一区蜜桃| av在线app专区| 少妇 在线观看| 久久精品久久精品一区二区三区| 久久久久视频综合| 一区二区三区精品91| 国产av精品麻豆| 伦精品一区二区三区| 考比视频在线观看| 欧美97在线视频| 久久人人爽人人片av| 九九久久精品国产亚洲av麻豆| av.在线天堂| 日韩精品免费视频一区二区三区 | 国产在线视频一区二区| 亚洲国产欧美日韩在线播放| 国产精品三级大全| a级毛片在线看网站| 亚洲欧美成人精品一区二区| 亚洲av综合色区一区| av国产精品久久久久影院| 亚洲av男天堂| 午夜福利网站1000一区二区三区| 欧美激情 高清一区二区三区| 老女人水多毛片| 啦啦啦啦在线视频资源| 国产日韩欧美视频二区| 啦啦啦中文免费视频观看日本| 成人手机av| 看十八女毛片水多多多| 91久久精品电影网| 久久这里有精品视频免费| 国产亚洲一区二区精品| 啦啦啦视频在线资源免费观看| 少妇人妻久久综合中文| 亚洲国产精品一区二区三区在线| 我的女老师完整版在线观看| 久久精品熟女亚洲av麻豆精品| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜爱| 男女高潮啪啪啪动态图| 精品人妻一区二区三区麻豆| 亚洲丝袜综合中文字幕| 亚洲五月色婷婷综合| 亚洲高清免费不卡视频| 免费久久久久久久精品成人欧美视频 | 国产 一区精品| 少妇高潮的动态图| 国产白丝娇喘喷水9色精品| 伦理电影大哥的女人| 国产av精品麻豆| 伦精品一区二区三区| 亚洲国产成人一精品久久久| 亚洲精品国产av成人精品| 国产视频首页在线观看| 亚洲欧洲国产日韩| av网站免费在线观看视频| 亚洲经典国产精华液单| 日本黄色日本黄色录像| 亚洲,欧美,日韩| 国产一区二区三区综合在线观看 | 黑人猛操日本美女一级片| 久久99蜜桃精品久久| 久久国产精品大桥未久av| 久久狼人影院| 妹子高潮喷水视频| 精品少妇内射三级| 精品久久久久久久久av| 日本与韩国留学比较| www.av在线官网国产| 精品人妻一区二区三区麻豆| av女优亚洲男人天堂| 亚洲精品久久午夜乱码| 十分钟在线观看高清视频www| 日韩一本色道免费dvd| 99国产精品免费福利视频| 国产69精品久久久久777片| 国产成人a∨麻豆精品| 精品久久国产蜜桃| 观看av在线不卡| 人妻夜夜爽99麻豆av| 不卡视频在线观看欧美| videos熟女内射| 国产精品99久久99久久久不卡 | 两个人的视频大全免费| 狠狠精品人妻久久久久久综合| 成年人午夜在线观看视频| 人妻人人澡人人爽人人| 国产精品久久久久成人av| 欧美日韩av久久| 国产精品国产三级国产专区5o| 亚洲精品成人av观看孕妇| 热re99久久精品国产66热6| 日韩 亚洲 欧美在线| 日韩av免费高清视频| 大香蕉久久成人网| 观看av在线不卡| 成人国产麻豆网| 亚洲美女视频黄频| av免费在线看不卡| 成人毛片60女人毛片免费| 亚洲经典国产精华液单| 你懂的网址亚洲精品在线观看| 亚洲一区二区三区欧美精品| 久久人妻熟女aⅴ| 久久99一区二区三区| 高清黄色对白视频在线免费看| 亚洲精品自拍成人| 视频在线观看一区二区三区| 国产成人一区二区在线| 久久久久久久久久久丰满| 免费看av在线观看网站| 亚洲国产欧美日韩在线播放| 丰满饥渴人妻一区二区三| 美女cb高潮喷水在线观看| 一本久久精品| 国产在视频线精品| a级毛片黄视频| 春色校园在线视频观看| 丰满迷人的少妇在线观看| 人妻一区二区av| 国产亚洲最大av| 啦啦啦在线观看免费高清www| 国产综合精华液| 最近2019中文字幕mv第一页| 午夜福利视频在线观看免费| 国产成人精品福利久久| 国产一区二区三区综合在线观看 | 欧美成人午夜免费资源| 男女无遮挡免费网站观看| 飞空精品影院首页| 免费黄网站久久成人精品| 另类精品久久| 久久99热6这里只有精品| 亚洲国产成人一精品久久久| 亚洲国产毛片av蜜桃av| 亚洲欧美成人精品一区二区| 欧美一级a爱片免费观看看| 伊人久久国产一区二区| 国产视频内射| 国产精品一区www在线观看| 老司机亚洲免费影院| 日韩免费高清中文字幕av| 国产精品熟女久久久久浪| a级毛片黄视频| 飞空精品影院首页| 国产日韩欧美在线精品| 国产色爽女视频免费观看| 狠狠精品人妻久久久久久综合| 性色av一级| 九色亚洲精品在线播放| 国产免费一区二区三区四区乱码| 天美传媒精品一区二区| 国产精品一区二区在线不卡| 久久久国产一区二区| 亚洲精品色激情综合| 精品国产一区二区久久| av福利片在线| 69精品国产乱码久久久| 国产午夜精品久久久久久一区二区三区| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区黑人 | 国产男女内射视频| 亚洲精品视频女| 国产极品天堂在线| 亚洲国产av影院在线观看| 大码成人一级视频| 久久青草综合色| 五月天丁香电影| 另类精品久久| 久久狼人影院| 国产成人精品在线电影| 在线天堂最新版资源| 成人亚洲精品一区在线观看| 欧美激情 高清一区二区三区| 亚洲经典国产精华液单| 大香蕉久久成人网| 一本色道久久久久久精品综合| 综合色丁香网| 午夜福利视频精品| 国产成人a∨麻豆精品| 免费av中文字幕在线| a级毛片免费高清观看在线播放| 免费大片18禁| av女优亚洲男人天堂| 亚洲国产欧美日韩在线播放| 精品亚洲成a人片在线观看| 99久久精品一区二区三区| 亚洲av成人精品一区久久| 青春草国产在线视频| 韩国av在线不卡| 国产69精品久久久久777片| 久久精品久久久久久久性| 久久久久久久久久久丰满| 妹子高潮喷水视频| av在线app专区| 中国美白少妇内射xxxbb| 国产av精品麻豆| 亚洲精品aⅴ在线观看| 男女边摸边吃奶| 少妇丰满av| 亚洲精品国产色婷婷电影| 国产成人av激情在线播放 | 汤姆久久久久久久影院中文字幕| 一边亲一边摸免费视频| 亚洲四区av| 欧美97在线视频| 高清黄色对白视频在线免费看| 最近最新中文字幕免费大全7| 中文字幕av电影在线播放| 丰满饥渴人妻一区二区三| 97超碰精品成人国产| 亚洲欧美一区二区三区黑人 | 亚洲精品456在线播放app| 久久久久久久亚洲中文字幕| 成年人免费黄色播放视频| 男女免费视频国产| 大香蕉久久成人网| 午夜激情久久久久久久| 国产精品久久久久久av不卡| 性色avwww在线观看| 亚洲情色 制服丝袜| 热re99久久国产66热| 亚洲精品av麻豆狂野| 国产亚洲精品久久久com| 久久久久久久大尺度免费视频| 久久毛片免费看一区二区三区| 国产熟女欧美一区二区| 草草在线视频免费看| 一区二区三区四区激情视频| 天堂8中文在线网| 日韩免费高清中文字幕av| 久久午夜福利片| 亚洲精品第二区| 男人爽女人下面视频在线观看| 亚洲美女搞黄在线观看| 如何舔出高潮| 久久精品国产鲁丝片午夜精品| 色5月婷婷丁香| 欧美精品人与动牲交sv欧美| 日韩中字成人| 黑人欧美特级aaaaaa片| 久久久久网色| 日韩视频在线欧美| av.在线天堂| 亚洲怡红院男人天堂| 激情五月婷婷亚洲| 一区二区三区乱码不卡18| 亚洲人成网站在线播| 人妻一区二区av| 这个男人来自地球电影免费观看 | 日本与韩国留学比较| 亚洲国产欧美在线一区| 卡戴珊不雅视频在线播放| 伊人久久国产一区二区| 人体艺术视频欧美日本| 亚洲欧洲国产日韩| 啦啦啦视频在线资源免费观看| 国产69精品久久久久777片| 欧美老熟妇乱子伦牲交| 制服诱惑二区| 免费观看无遮挡的男女| 国产亚洲精品第一综合不卡 | 九九在线视频观看精品| 国产欧美日韩综合在线一区二区| 午夜久久久在线观看| 曰老女人黄片| 国产不卡av网站在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品三级大全| 一区二区三区乱码不卡18| 久久精品国产亚洲av天美| 欧美97在线视频| 国产欧美亚洲国产| 欧美精品人与动牲交sv欧美| 99热国产这里只有精品6| 亚洲欧洲日产国产| 国产 精品1| a级片在线免费高清观看视频| 制服丝袜香蕉在线| 精品久久蜜臀av无| 一本久久精品| 日本91视频免费播放| 欧美日韩一区二区视频在线观看视频在线| 五月玫瑰六月丁香| 老司机影院成人| tube8黄色片| 欧美精品人与动牲交sv欧美| 一边亲一边摸免费视频| 亚洲国产精品成人久久小说| 成人综合一区亚洲| 欧美xxⅹ黑人| 日韩中文字幕视频在线看片| 高清欧美精品videossex| 大又大粗又爽又黄少妇毛片口| 男的添女的下面高潮视频| 亚州av有码| 精品卡一卡二卡四卡免费| 在线观看www视频免费| 国产精品不卡视频一区二区| 久久99精品国语久久久| 搡老乐熟女国产| 国产精品国产三级国产专区5o| 亚洲精华国产精华液的使用体验| 国产白丝娇喘喷水9色精品| 日日撸夜夜添| 99久国产av精品国产电影| 天天影视国产精品| 欧美三级亚洲精品| 国产国拍精品亚洲av在线观看| 一本大道久久a久久精品| av在线老鸭窝| 久久精品人人爽人人爽视色| 纯流量卡能插随身wifi吗| 亚洲在久久综合| 草草在线视频免费看| 能在线免费看毛片的网站| 乱码一卡2卡4卡精品| 高清午夜精品一区二区三区| 我的女老师完整版在线观看| 午夜91福利影院| 国产乱人偷精品视频| 亚洲久久久国产精品| 亚洲av中文av极速乱| 中国国产av一级| 久久精品国产亚洲网站| 女人久久www免费人成看片| freevideosex欧美| 青春草亚洲视频在线观看| 国产无遮挡羞羞视频在线观看| 国产毛片在线视频| 国产欧美日韩综合在线一区二区| 日日爽夜夜爽网站| 99久久人妻综合| 久久国产精品大桥未久av| 看非洲黑人一级黄片| 色视频在线一区二区三区| 热99久久久久精品小说推荐| 午夜激情久久久久久久| 99久久综合免费| 婷婷色av中文字幕| 欧美+日韩+精品| 婷婷色麻豆天堂久久| 各种免费的搞黄视频| 一区二区三区免费毛片| 成年av动漫网址| 精品少妇黑人巨大在线播放| 国产黄色免费在线视频| 国产精品欧美亚洲77777| 大香蕉97超碰在线| 精品亚洲成a人片在线观看| 久久久国产精品麻豆| 少妇丰满av| 日日摸夜夜添夜夜爱| 99久久人妻综合| 亚洲av二区三区四区| 日本色播在线视频| 青青草视频在线视频观看| 少妇人妻精品综合一区二区| 一级毛片aaaaaa免费看小| 久久久亚洲精品成人影院| 国国产精品蜜臀av免费| 黄色视频在线播放观看不卡| 欧美97在线视频| 在线看a的网站| 国产一区亚洲一区在线观看| 亚洲av免费高清在线观看| 丁香六月天网| 国产av码专区亚洲av| 色吧在线观看| 插阴视频在线观看视频| 久久久a久久爽久久v久久| 免费观看在线日韩| 欧美精品国产亚洲| 亚洲国产色片| 亚洲精品久久成人aⅴ小说 | 麻豆成人av视频| 国产黄频视频在线观看| 人人澡人人妻人| 欧美 日韩 精品 国产| 久久久久国产网址| 26uuu在线亚洲综合色| 免费观看无遮挡的男女| 99热这里只有精品一区| 久久国产精品大桥未久av| av免费在线看不卡| 亚洲av免费高清在线观看| 国产日韩欧美亚洲二区| 精品酒店卫生间| 精品久久久久久电影网| 能在线免费看毛片的网站| 国产色爽女视频免费观看| 日本黄大片高清| 亚洲欧洲日产国产| a级毛色黄片| 一个人看视频在线观看www免费| 3wmmmm亚洲av在线观看| 嫩草影院入口| 久久精品国产鲁丝片午夜精品| 精品国产一区二区三区久久久樱花| 成人国语在线视频| 男人爽女人下面视频在线观看| 99久久精品一区二区三区| 免费人妻精品一区二区三区视频| 最近手机中文字幕大全| 看免费成人av毛片| 又粗又硬又长又爽又黄的视频| 婷婷色av中文字幕| 成人亚洲精品一区在线观看| 日韩大片免费观看网站| 在线 av 中文字幕| 少妇被粗大的猛进出69影院 | 精品亚洲乱码少妇综合久久| 你懂的网址亚洲精品在线观看| 777米奇影视久久| 日韩亚洲欧美综合| 十八禁高潮呻吟视频| 精品人妻熟女毛片av久久网站| 少妇被粗大猛烈的视频| 啦啦啦视频在线资源免费观看| 蜜桃久久精品国产亚洲av| 中文字幕人妻丝袜制服| 色哟哟·www| a级毛色黄片| 国产在线免费精品| 精品午夜福利在线看| 人妻系列 视频| 国产成人精品福利久久| 亚洲一级一片aⅴ在线观看| 美女国产高潮福利片在线看| 中国美白少妇内射xxxbb| 有码 亚洲区| 女性被躁到高潮视频| 91精品一卡2卡3卡4卡| 亚洲天堂av无毛| 99久久综合免费| 久久久国产欧美日韩av| 2022亚洲国产成人精品| 高清不卡的av网站| 国产免费又黄又爽又色| freevideosex欧美| 在线观看免费高清a一片| 99九九在线精品视频| 韩国高清视频一区二区三区| 夫妻性生交免费视频一级片| 成人黄色视频免费在线看| 夫妻午夜视频| 亚洲av不卡在线观看| 午夜福利,免费看| 美女主播在线视频| 少妇熟女欧美另类| 成人亚洲精品一区在线观看| 大码成人一级视频| 在线观看www视频免费| 美女大奶头黄色视频| videossex国产| 国产成人精品福利久久| 在线观看国产h片| 考比视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产乱人偷精品视频| 国产午夜精品久久久久久一区二区三区| 免费观看无遮挡的男女| 国产成人一区二区在线| 国产片内射在线| 国产精品国产三级国产专区5o| 熟女电影av网| 青青草视频在线视频观看| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 高清午夜精品一区二区三区| 最近最新中文字幕免费大全7| 免费看av在线观看网站| 久久99一区二区三区| 久久ye,这里只有精品| 久久久久久久国产电影| 国产白丝娇喘喷水9色精品| 亚洲av.av天堂| 精品人妻熟女av久视频| 在线天堂最新版资源| 视频中文字幕在线观看| 另类亚洲欧美激情| 成人国产av品久久久| 亚洲精品乱码久久久久久按摩| 久久国内精品自在自线图片| 少妇高潮的动态图| 在线观看国产h片| 久久国产精品男人的天堂亚洲 | 国产成人精品无人区| 亚洲婷婷狠狠爱综合网| 亚洲在久久综合| 夫妻午夜视频| av在线老鸭窝| 亚洲av.av天堂| 国产精品久久久久久av不卡| 两个人免费观看高清视频| 国产成人91sexporn| 99热6这里只有精品| 国产男女超爽视频在线观看| 国产精品一区二区在线不卡| 久久久久久久亚洲中文字幕| 最黄视频免费看| 99热6这里只有精品| 日本黄大片高清| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 九九久久精品国产亚洲av麻豆| 超碰97精品在线观看| 中文天堂在线官网| 国产精品熟女久久久久浪| 国产男女超爽视频在线观看| 日本爱情动作片www.在线观看| 在线播放无遮挡| 精品久久久久久久久av| 精品久久久噜噜| 亚洲精品乱久久久久久| 亚洲精品乱码久久久v下载方式| 国产一区二区三区av在线| a 毛片基地| 欧美日韩一区二区视频在线观看视频在线| av国产精品久久久久影院| 亚洲欧美成人精品一区二区| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品国产精品| 成人黄色视频免费在线看| 国产视频内射| 亚洲激情五月婷婷啪啪| a级毛片黄视频| 欧美三级亚洲精品| 国产在视频线精品| 人人妻人人澡人人爽人人夜夜| 午夜福利视频精品| 大香蕉久久网| 丰满迷人的少妇在线观看| 赤兔流量卡办理| 国产男人的电影天堂91| 女的被弄到高潮叫床怎么办| 天天躁夜夜躁狠狠久久av| 国产精品一区www在线观看| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 久热这里只有精品99| 亚州av有码| 蜜桃久久精品国产亚洲av| 久久久午夜欧美精品| 欧美日韩av久久| 人体艺术视频欧美日本| 久久精品久久精品一区二区三区| 欧美+日韩+精品| 色哟哟·www| 秋霞伦理黄片| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 亚洲伊人久久精品综合| 久久99热6这里只有精品| 欧美日韩成人在线一区二区| 97精品久久久久久久久久精品| 日韩人妻高清精品专区| a 毛片基地| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 成年人午夜在线观看视频| 国产免费一级a男人的天堂| 国产在视频线精品| 精品熟女少妇av免费看| 国内精品宾馆在线| 亚洲国产精品一区二区三区在线| 亚洲欧美清纯卡通| 亚洲丝袜综合中文字幕| 大片电影免费在线观看免费| 日韩电影二区| 多毛熟女@视频| 赤兔流量卡办理| 亚洲国产精品一区三区| 国产毛片在线视频| 99九九线精品视频在线观看视频| 国产午夜精品一二区理论片|