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

    對(duì)稱定常微流邊界層

    2013-07-02 12:09:02黃志圣詹華稅
    關(guān)鍵詞:有界邊界條件常數(shù)

    黃志圣,葉 霞,詹華稅

    (集美大學(xué)理學(xué)院,福建 廈門 361021)

    0 引言

    考慮如下對(duì)稱定常微流邊界層方程

    盡管本文的方法是類似于Oleinik的方法,然而,v(y)是y的函數(shù),這使得相應(yīng)的轉(zhuǎn)換和計(jì)算變得更加復(fù)雜,一些限制條件作用在v(y)上使其仍然能夠應(yīng)用Oleinik的線性化方法.此外,如果允許v(y)或其導(dǎo)數(shù)是一個(gè)無界函數(shù),則不能應(yīng)用Oleinik的線性化方法,在這種情況下,得到同樣的結(jié)果看起來非常困難.近年來關(guān)于定常Prandtl系統(tǒng)全局解的存在性方面重要的研究進(jìn)展可參考文獻(xiàn) [5-8],其他相關(guān)工作,可以參考文獻(xiàn) [9].

    本文的目的是研究系統(tǒng) (1)— (2)解的適定性問題.為此引入Crocco變換ξ=x,η=u(x,y)/U(x),w(ξ,η)=uy(t,y)/U(x).經(jīng)過簡單的計(jì)算,可得如下w(ξ,η)的方程和邊界條件:

    與文獻(xiàn) [1]相比,方程 (3)增加了非線性項(xiàng)2vηw2wη+vηηw3,邊界條件 (4)增加了非線性項(xiàng)vηw2.在本文的討論中,要利用Oleinik線性化的方法,討論問題 (1)— (2)解的適定性問題.為此,必須增加以下限制性條件:是有界的,vη≤0,vηη=0.則方程 (3)變?yōu)?

    定義1稱函數(shù)u(x),v(x)為問題(1)的解,如果u在內(nèi)連續(xù),v在D內(nèi)連續(xù)并且在ˉD內(nèi)關(guān)于y連續(xù),u,v的導(dǎo)數(shù)在D內(nèi)連續(xù);方程 (1)對(duì)u,v在D內(nèi)成立,并滿足條件 (2).

    定義2稱函數(shù)w(ξ,η)為問題(4)—(5)的解,如果:w(ξ,η)在內(nèi)連續(xù),并且在Ω內(nèi)有連續(xù)偏導(dǎo)數(shù)wξ,wη,wηη,其中wηη關(guān)于η在η=0處連續(xù);方程 (5)對(duì)w(ξ,η)在Ω內(nèi)成立,并且滿足條件 (4).

    對(duì)任意的函數(shù) f(η,ξ),引進(jìn)如下的記號(hào) fk=fk(η)=f(kh,η),h 〉 0 且為常數(shù),k=0,1,2,…,[X/h].代替方程 (5)和邊界條件 (4),將考察下面的常微分方程組:

    其中,當(dāng)k=0時(shí),μk=0;當(dāng)k≥0時(shí),μk是足夠大的常數(shù),并且U(x)、r(x)是二階連續(xù)導(dǎo)數(shù),v0x(x)有界.在以下的論述中,約定Mi,μ,δ為不依賴于h的正常數(shù).

    1 重要引理

    引理1 設(shè)wk(η)是常微分方程組 (6)的解,當(dāng)0≤η≤1時(shí),wk(η)是連續(xù)的;當(dāng)0≤η〈1時(shí),wk(η)是無窮可微的,vη≤0.wk(η)滿足下面的估計(jì):

    這里,kh≤X,h≤h0,h〉0為常數(shù).

    證明 通過構(gòu)造如下的系統(tǒng) (8)的解,取ε→0的極限來獲得方程組 (6)的解.

    假設(shè)方程組 (8)有一個(gè)解是wk.滿足當(dāng)η=0時(shí),wk為正.令Vk=M3(1-η)e-αkh,α〉0為常數(shù),發(fā)現(xiàn),因?yàn)閁kx〉0,0〈h'〈h,并且所以,當(dāng)η〈1足夠大,M3足夠小時(shí),由假設(shè) vη≤0,有:如果M3足夠小,vη≤0,則令yk=Vk-wk,有:當(dāng)η〈1時(shí),Lε,k(V)-Lε,k(w)〉 0;當(dāng) η =0 時(shí),λε,k(V)- λε,k(w)〉 0 .

    由以上不等式可以推出 yk≤0.事實(shí)上,考慮 Sk=yke-βkh.當(dāng) η〈1,?h″,0〈h″〈h時(shí),由Lε,k(V)-Lε,k(w)〉 0 ,λε,k(V)- λε,k(w)〉 0 ,有

    因?yàn)閁x〉0,當(dāng)η=0時(shí),wk,Vk為正的.由最大值原理,從式 (9)可得Sk≤0.具體來說,Sk在[0,1]的內(nèi)部取不到正的最大,因?yàn)楫?dāng)β〉0足夠大、h足夠小和M3足夠小時(shí),有:這和式 (9)矛盾.在η=0這點(diǎn),因?yàn)镾kη(0)≤0和式 (10)相互矛盾,所以Sk也不可能有正極大值.因此,有Sk≤0,yk=Vk-wk≤0.由此可得,當(dāng) kh≤X時(shí),有wk(η)≥ M3(1- η)e-αkh.

    選擇任意一個(gè)正常數(shù) M4,有下列不等式成立:2/M4.假設(shè)ψ(w)是(-∞,+∞)上的充分光滑函數(shù),使其滿足當(dāng)w≥M4時(shí),ψ(w)=w;當(dāng)w≤M4/4時(shí),ψ(w)=M4/2;當(dāng)M4/4≤w≤M4時(shí),0≤ψ'(w)≤1.

    考慮由方程組 (8)及以下的邊界條件 (11)組成的系統(tǒng)

    當(dāng)ε〉0時(shí),方程 (8)和條件 (11)解的存在性的證明由Leray-Schander定理的證明方法給出.

    Leray-Schauder定理 在Banach空間X中,考慮一族映射y=T(x,k),這里x,y∈X;k是實(shí)參數(shù),a≤k≤b.設(shè):1)T(x,k)定義在x∈X,a≤k≤b;2)對(duì)任意固定的k,T(x,k)在X上連續(xù),即?ε〉0,x0∈X,?δ〉0,如果在X上的有界子集關(guān)于k一致連續(xù),即?ε〉0和任意有界集合X0?X,?δ〉0,當(dāng)時(shí),有是緊集;5)對(duì)x-T(x,k)=0的任意解x,有在X內(nèi)有唯一解.

    現(xiàn)在,做方程 (6)解的存在性證明.考慮下面含參量γ的微分方程系統(tǒng):

    邊界條件

    當(dāng)γ=0時(shí),問題 (12)— (13)為線性的;當(dāng)γ=1時(shí),則該問題轉(zhuǎn)化為方程 (8)和邊界條件(11)的問題.

    對(duì)于這個(gè)系統(tǒng),證明Leray-Schauder定理.考慮T(θ,γ)=w,它把C2([0,1])中的任意矢量函數(shù)θ映射到w≡(w1,…,wm),m=[X/m],其中w是下面的線性微分方程系統(tǒng)的解:

    邊界條件

    當(dāng)γ=0時(shí),問題 (14)— (15)有唯一的解.因?yàn)樵搯栴}是線性的,當(dāng)h足夠小時(shí),在問題(14)中wk的系數(shù)是非正的,并且φ(x)≤0,Ukx〉0.

    非線性問題 (12)— (13)的解wk關(guān)于γ一致有界,并且它的二階導(dǎo)數(shù)也是一致有界的.首先,對(duì)wk做下界的估計(jì).令 V0=M5(1-η)e-αkh.當(dāng) M5足夠小時(shí),可得因?yàn)镸5與γ,h,ε無關(guān).當(dāng)0〈η≤1時(shí),V0代入式 (14),則式 (14)的等式左邊大于零.因此,對(duì)于yk=Vk0-wk,可得,當(dāng)γ,kh≤X,0≤η≤1,0〈ε≤1時(shí),有yk≤0,wk≥Vk0.

    下面對(duì)問題(12)—(13)的解wk做上界的估計(jì).為此把它轉(zhuǎn)化為另一個(gè)未知函數(shù)的形式wk=(M6-eβη)eκkhˉwk,其中κ,β是正常數(shù).選擇適當(dāng)大的M6,β,κ,可以得到以下的方程其中〉0,以及邊界條件

    由式 (12)中關(guān)于wkη的一階方程和邊界條件 (13)中在η=0點(diǎn)wkη的估計(jì)得到wkη關(guān)于γ的一致估計(jì).由相同的方法,可以得出問題 (14)— (15)的解wk(η)和它到三階導(dǎo)數(shù)關(guān)于γ的一致估計(jì).在這些估計(jì)中的常數(shù)依賴于算子 T(θ,η)是 由在 C2([0,1])中函數(shù)θ的有界集映射到w的緊集.因此,在C2([0,1])中關(guān)于方程 (8)— (11)的問題在ε〉0時(shí)解的存在性,即Leray-Schauder定理的結(jié)論.

    前面已經(jīng)證明了問題 (8)的解wk(η)的下界關(guān)于 ε,h的一致估計(jì),即,wk(η)≥ M3(1-η)e-αkh.下面,證明其解的上界的一致估計(jì),即其中,M8,μ是正常數(shù)且與ε,h無關(guān).取μ〈1并且當(dāng)0≤η〈1時(shí),σ〉1.當(dāng)0≤η≤1,M8足夠大時(shí),可以得到Lε,k(V1)〈0.

    另外,對(duì)于μ〈e-1/2,如果M8足夠大,υη≤0,則可以得到λε,k(V1)〈0.由上面的不等式可以得到 Lε,k(V1)-Lε,k(w)〈 0 ,λε,k(V1)- λε,k(w)〈 0 .當(dāng)0 〈 η 〈1,0≤kh≤X 時(shí),設(shè) Sk=Vk1-wk滿足不等式:因?yàn)?,?dāng)η 〈1時(shí),所以當(dāng)0〈η〈1,kh≤X,h足夠小時(shí),由以上的不等式容易推出Sk≥0,Vk1≥wk.

    當(dāng)η=0時(shí),由方程組 (8)結(jié)合wk兩邊的估計(jì),可以得到在0≤η≤1-δ上,wkη,wkηη關(guān)于 ε 的一致估計(jì)成立.對(duì)方程 (8)關(guān)于η微分,可以證明wk的導(dǎo)數(shù)到給定的階數(shù)在0≤η≤1-δ上關(guān)于ε一致有界.因此,問題 (8)的一族解wk依賴于ε,0〈ε≤1,可以抽取一個(gè)序列wk,例如wk和它們?nèi)我饨o定階數(shù)的導(dǎo)數(shù),當(dāng)ε→0時(shí),其在0≤η≤1-δ上一致收斂.因?yàn)閂≤wk≤V1,則這個(gè)有限函數(shù)在0≤η≤1-δ上連續(xù);當(dāng)η〈1時(shí),有限函數(shù)滿足方程組 (6);當(dāng)η=1時(shí),它不存在.顯然,有限函數(shù)wk對(duì)于式 (7)的估計(jì)式成立.存在性證畢.

    對(duì)方程(8)—(11)系統(tǒng)的任意解?wk,下證?wk≥Vk.令?yk=Vk-?wk,選擇適當(dāng)小的M3,有當(dāng)η〈1時(shí),有

    為了得到wk(η)有更好的估計(jì),下面證明引理2.

    引理2 設(shè)wk(η)是問題 (6)的解,有下列的估計(jì)成立:當(dāng)0≤η≤1,0≤kh≤X時(shí),

    當(dāng)0≤η≤1時(shí),

    當(dāng)1-δ≤η≤1時(shí),

    證明 首先,建立估計(jì)式 (18).令Φk=M12(1-η)σe-αkh,則當(dāng)α〉0足夠大,M12足夠小和μ〈e-1/2,有其中另外,如果 M12足夠小,則當(dāng)0≤η〈1,kh≤X時(shí),設(shè) Sk=Φk-wk,因?yàn)?Lk(Φk)-Lk(wk)〉0,λk(Φk)-λk(wk)〉0,則有不等式:

    當(dāng)β〉0為常數(shù)并且其足夠大時(shí),把式(21)、式(22)轉(zhuǎn)化為新函數(shù)由以上類似的證明可以看出,上面的不等式中不可能取正值.因此,Sk=Φk-wk≤0,式(18)估計(jì)成立.用相似的方法證明式 (19).令 φ =K0(1- η)σ ,則 L0(φ)= υK20(1- η)2σ2[-K0/(2σ(1- η))-K0/(4σ3(1- η))]-(1- η)(1+ η)Ux(0)K0(- σ +1/(2σ))+2υηK30(1- η)2σ2(- σ +1/(2σ)).因?yàn)棣訩20=4Ux(0),υη≤0,所以當(dāng)0≤η〈1時(shí),L0(φ)〈0.另外,如果σ(0)足夠大,即σ(0)〉M11,其中M11的取值只與問題 (4)— (5)的數(shù)值有關(guān),可得當(dāng)0≤η〈1時(shí),函數(shù)S=φ-w0滿足不等式因?yàn)棣咋恰?,φηη〈0和S(1)=0,υη≤0,所以S的系數(shù)不可能為正.可以推出當(dāng)0≤η≤1時(shí),S≥0.因此,不等式 (19)成立.

    顯然,如果M10〉0且適當(dāng)大,M10〈σ,而且1-M10/σ≥d〉0,則L0(Q)〉0.取δ〉0,有1-M10/(σ(1-δ))=d.則當(dāng)1-δ≤η〈1時(shí),1-M10/σ≥d,M10〈σ.由此可得,當(dāng)1-δ≤η〈1時(shí),L0(Q)〉0,并且根據(jù)式 (20)和所選擇的d,有并且Q(1)=w0(1)=0.對(duì)于S=Q-w0,當(dāng)1-δ≤η〈1時(shí),得 υ(w0)2Sηη+A0Sη+Qηηυ(w0+Q)S+2υη(w0)2Sη+2υηQη(w0+Q)S 〉 0 ,S(1- δ)≤0 ,S(1)=0.

    因?yàn)镼ηη〈0,當(dāng)M10適當(dāng)大時(shí),Qη〉0.由這些不等式可以推出S≤0,Q≤w0.因此不等式 (20)成立.

    引理3 當(dāng)k=0時(shí),設(shè)w0是問題 (6)的解,則它的導(dǎo)數(shù)w0η滿足估計(jì):當(dāng)0≤η〈1時(shí),-M13σ≤w0η≤M14.

    證明 由引理1和引理2,可以得到:

    令Hk=wk-M2(1-η)σ≤0.當(dāng)η≤1時(shí),Hk≤0,Hk(1)=0;并且當(dāng)η〈1時(shí),Hk(η)是可微的.因此,當(dāng)n→∞,ηkn→1時(shí),Hkη(ηkn)≥0.由此可得

    如果考慮 M9(1-η)σ -wk,當(dāng) n→∞,?ηkn→1時(shí),有

    當(dāng)0≤η〈1時(shí),函數(shù)z0≡w0η滿足方程

    對(duì)w0η做下界的估計(jì).令θ=-M20σ-M21,M20,M21=const〉0.然后,應(yīng)用式 (19),如果M21足夠大,當(dāng)0≤η 〈1時(shí),由假設(shè) υη≤0,有L0(θ)≥-(1-η)2(U2x+2υηK20)M20σ +M21[(1- η)(1+η)U0x-2υηK20(1-η)2]〉0.因此,當(dāng)0≤η〈1時(shí),有L0(θ-z0)〉0.如果M20〉M15,則在η=η0n上θ-z0〈0.因此容易得,當(dāng)0≤η〈1時(shí),有θ-z0〈0.則當(dāng)0≤η〈1時(shí),有w0η=z0≥θ=-M20σ-M21≥-M13σ .

    證明 不等式 (27)— (28)用歸納法證明.當(dāng)k=0時(shí),式 (27)已經(jīng)在上面的引理中證明過.下面要證明如果kh≤X1,當(dāng)k-1對(duì)于所選擇的常數(shù)M22、M23、M24成立時(shí),則k也成立.令wkη=zk,h-1(wk-wk-1)=rk,做出zk,rk,k≥1的方程.對(duì)方程 (6)關(guān)于η求導(dǎo),得

    在η=0對(duì)于方程 (6)的邊界條件,有

    而后,把方程 (6)中wk減去wk-1的差除以h,得:)

    類似地,從方程 (6)的邊界條件可以得到

    考慮函數(shù)φk±rk=qk±,當(dāng)0≤η〈1時(shí),如果kh≤x1,有:

    如果kh≤x2,還可以得到

    因?yàn)樵谑?(34)中,當(dāng)kh≤x1時(shí),qk±的系數(shù)是非正的,由式 (34)— (35)可得,當(dāng)0≤η≤1時(shí),qk±≥0;并且kh≤min(x1,x2)=x3,則有

    在式(28)中令M27等于M28.因?yàn)橛袣w納假設(shè)qk±-1≥0,所以式(34)成立.

    下面對(duì)zk做估計(jì).令F=-M29σ,計(jì)算Pk(F).由wk的估計(jì)式 (18)和式 (7),如果M29適當(dāng)大,當(dāng) kh≤x3時(shí),有:選擇的M29與k,h無關(guān),當(dāng)η=ηkn時(shí),F(xiàn)k-zk≤0;當(dāng)η=0時(shí),F(xiàn)k-zk≤0,其中ηkn對(duì)于不等式 (24)成立.

    當(dāng)0≤η〈1時(shí),對(duì)于差分Fk-zk=Sk,有不等式

    并且有:在式 (37)中Sk的系數(shù)等于Bk+Akη+2υηwkFkη+4υη(Fk+zk).顯然,如果M29足夠大,則其系數(shù)是負(fù)的.用上面的形式取M22等于M29,則Sk-1≤0.因?yàn)樵谑?(37)— (38)中Sk不可能取正值,所取Sk≤0.當(dāng)kh≤x3時(shí),有wkη=zk≥-M22σ,其中M22依賴于問題 (4)— (5)的數(shù)值.可以假設(shè)M22≥M13.

    下面,對(duì)zk做上界的估計(jì).令Ψ=-M30e-αη,如果kh≤x4≤x3,當(dāng)M30足夠小,α足夠大時(shí),不妨設(shè)x4足夠小,所以可以推出 Bkη=-Uk,設(shè) wk-rk〈0,由以上的假設(shè),有 Pk(Ψ)(α 〉0):如果 M30足夠小,kh≤x5,x5=const〉 0,可得當(dāng)0≤η 〈1時(shí),對(duì)于 Sk1=Ψk-zk,有不等式:

    為了證明 Sk1≥0,做新函數(shù) Sk2=Sk1e-βkh,易得如果β〉0足夠大,則當(dāng)M30足夠小時(shí),S2k的系數(shù)在第一個(gè)不等式中是正的.因?yàn)槿绻O(shè)S02≥0,當(dāng)0≤η≤1,kh≤x6時(shí),其中x6=min(x4,x5),則S2k不可能取負(fù)值.因此,當(dāng)kh≤x6時(shí),有S2k≥0,zk≤-M30e-αη.

    在η=1這點(diǎn)附近,可以做進(jìn)一步的估計(jì).令F1=-M31σ并且計(jì)算Pk(F1).如果kh≤x7,0〈1-η 〈δ1,M31和 δ1足夠小時(shí),不妨設(shè) wk-rk〈0,有:Pk1(F)=- υ(wk)2M31[1/(2σ(1- η)2)-1/(4σ3(1- η)2)]+M31(1- η2)Ukx/(2σ(1- η))-BkM31σ +2υwkM231/(2(1- η))-2ηUxM31σ +Bkηwk-Ukrk- υη(wk)2M31/(σ(1- η))+4υηwkM231σ2≤M31[σ(- υM29/2-Bk+ υM2M31-2ηUxM31σ - υη(1-η)M29+4υηwkM31σ)+(υM29/(4σ)+(1+ η)Ukx/(2σ))]+Bkηwk-Ukrk〈 0 .其中為常數(shù).

    考慮差分Sk=Fk1-zk.當(dāng)1-δ1≤η〈1,kh≤x7時(shí),函數(shù)Sk滿足不等式

    如果M31適當(dāng)小,由式 (25)和zk≤M30e-αη,當(dāng)kh≤x6時(shí),有

    因此,選擇的M31只由問題 (4)— (5)的數(shù)值決定.令M23=min{M31,M30e-α(ln(μδ1))-1/2}.設(shè)M23〈M14,則Sk-1≥0.由式(39)、式(40)可得Sk≥0.當(dāng)k≥1時(shí),代入新函數(shù)ˉSk=Ske-βkh,可得:如果 β(h)〉 0 足夠大由此可得,當(dāng)1-δ1≤η 〈1時(shí),有

    不等式 (29)由式 (27)、(28)得來.由方程 (6)可得 υwkwkηη=(ηUk+ μkh)rk/wk-Akzk/wk-Bk-2υηwkzk.因?yàn)锽k≤N3kh,當(dāng)kh足夠小時(shí),有:

    2 解的存在唯一性

    證明 當(dāng)kh〈ξ≤(k+1)時(shí),把wk(η)=w(kh,η)當(dāng)作問題 (6)的解,可以得到一族wh(ξ,η),wh(kh(1- λ)+(k+1)hλ,η)=wk(η)(1- λ)+wk+1(η)λ,0 ≤ λ ≤1,k=0,1,2,… .

    根據(jù)引理1—引理4,wh(ξ,η)關(guān)于ξ滿足Lipschitz條件,并且當(dāng)0〈ξ〈X,0〈η〈1-ε,0〈ε〈1時(shí),wh(ξ,η)在η上關(guān)于h一階一致有界.由Arzela定理知,存在序列hi→0,使得子序列whi(ξ,η)在矩形區(qū)域上一致收斂于某一函數(shù)wh(ξ,η).因?yàn)镸9(1-η)σ≤wh≤M2(1-η)σ,whi可以在Ω內(nèi)一致收斂到w.

    由式 (27)—式 (28)可得w(ξ,η)存在有界的弱導(dǎo)數(shù) wξ,wη,wηη,并且 wwηη在 Ω 內(nèi)有界.不妨設(shè)序列whi(ξ,η)滿足在區(qū)域內(nèi)相應(yīng)的偏導(dǎo)數(shù)在 L2(Ωε)的意義下收斂于 wξ,wη,wηη.

    記wkh=wh(ξ,η)=w(kh,η).由式 (6)的方程可得:

    假設(shè)φ(ξ,η)是光滑函數(shù),且使得suppφ??Ω.令φk(η)=φ(kh,η),在式 (41)兩邊同乘以hφk,將所得的方程關(guān)于η從0到1積分,同時(shí)對(duì)k從1到m(h)求和,有

    為了證明式 (5)中w及其導(dǎo)數(shù)在Ω的子區(qū)域內(nèi)滿足H?lder條件,考慮下面的方程

    下證邊界條件 (4)也是滿足的.由式 (4)的第一個(gè)條件知,wh是一致連續(xù)的,并且當(dāng)η=1時(shí),問題(6)的條件成立.令因?yàn)楫?dāng)η〈1-δ時(shí),在h上一致有界,并且滿足式 (6)的第二個(gè)邊界條件,則:

    當(dāng)hi→0時(shí),在L2(Ω)的意義下 zh?z=υwwη+υηw2-v0w+c.因此,由式 (47)可得即邊界條件 (4)也滿足.

    定理2問題(4)—(5)存在唯一的解ω,有如下的性質(zhì):w在ˉΩ內(nèi)連續(xù),它的導(dǎo)數(shù)wξ,wη,wηη在Ω內(nèi)連續(xù);在Ω內(nèi)wηη≤0,w≥0;當(dāng)η=0時(shí)w〉0;wη關(guān)于η在η=0處連續(xù).

    證明 假設(shè)問題(4)—(5)在Ω內(nèi)的兩個(gè)解w1,w2有如上的性質(zhì).令則可以得到,如果α〉0且α足夠大,當(dāng)η=0時(shí),w1〉0,w2〉0.因此,當(dāng)η≤1時(shí),既沒有正的最大值也沒有負(fù)的最小值,即≡0,w1≡w2,唯一性得證.

    定理1和定理2作為推論,在下面的證明中得到定常對(duì)稱微流邊界層系統(tǒng) (1)— (2)解的存在唯一性.

    定理3問題(1)—(2)在D中的解滿足以下的性質(zhì):u/U,uy/U在中連續(xù)有界;當(dāng)x〉0,y〉0時(shí),u〉0;當(dāng)y→∞時(shí),u→U,u(x,0)=u(0,y)=0;當(dāng)y≥0時(shí),uy/U〉0;當(dāng)y→∞時(shí),uy/U→0,uyy,ux,uy關(guān)于y在中連續(xù)有界;v關(guān)于有界的y是有界的;υ(y)在D內(nèi)有直到三階的有界導(dǎo)數(shù);uyyy在ˉD中有界;uxy關(guān)于有界的y在D內(nèi)有界;uxy,uyyy在內(nèi)連續(xù);uyy/uy關(guān)于y在內(nèi)連續(xù).此外,

    證明 如果w(ξ,η)是問題 (4)— (5)的解,并且有定理1中的性質(zhì),可以利用Crocco變換,把問題 (4)— (5)的解轉(zhuǎn)換成問題 (1)— (2)的解的形式.根據(jù)Crocco變換,可以得到

    由w及其導(dǎo)數(shù)的性質(zhì)知,弱導(dǎo)數(shù)ux,uyy,uyyy在D內(nèi)有界.uxy對(duì)有界的y有界;由w,wξ,wη,wwηη的性質(zhì),知式 (48)成立.且由式 (51)知ux,uyy關(guān)于y連續(xù).

    按如下的方式定義函數(shù)v(x,y)

    關(guān)于y微分式 (52),可得vyuy+vuyy+uyux+uuxy-υu(píng)yyy-2υyuyy-υyyuy=0,即

    函數(shù)w(ξ,η)=uy/U滿足方程 (5),在方程 (5)中用u的導(dǎo)數(shù)代替w的導(dǎo)數(shù),可得

    把式 (54)乘以U后與式 (53)作和,可得uyvy+uyux+ruxuy/r=0,即

    式 (52)和式 (55)就組成了所考慮的系統(tǒng) (1).由式 (4)知,

    定理4 令 u,v是問題 (1)— (2)的解,其導(dǎo)數(shù) ux,uy,vy,uyy,uyyy,uxy在 D中連續(xù),u/U,uy/U在中連續(xù);當(dāng)y≥0,x〉0時(shí),uy〉0;當(dāng)y=0時(shí),uy/U〉0;當(dāng)y→∞時(shí),uy/U→0;uyy/uy,ux關(guān)于y在y=0點(diǎn)連續(xù);(uyyyuy-u2y)/u2>y≤0,則u,v是問題 (1)— (2)的唯一解.證明 此定理的證明與文獻(xiàn)[1]中定理3.1.8的方法一樣,這里不再敘述.

    [1]OLEINIK O A,SAMKHIN V N.Mathematical models in boundary layer theorem [M].Boca Raton:Chapman and Hall/CRC,1999.

    [2]ANDERSON I,TOFTEN H.Numerical solutions of the laminar boundary layer equations for power-law fluids[J].Non-Newton,F(xiàn)luid Mech,1989,32:175-195.

    [3]SCHLICHING H,GERSTEN K,KRAUSE E.Boundary layer theory[M].8th ed.New York:Springer,2004.

    [4]LI L,ZHAN H S.The study of micro-fluid boundary layer theory [J].WSEAS Transaction on Mathematics,2009,8(12):699-711.

    [5]GILBARG D,TRUDINGER N S.Elliptic partial differential equations of second order[M].New York:Springer-Verlag,1977.

    [6]ZHANG J W,ZHAO J N.On the global existence and uniqueness of solutions to non-stationary boundary layer system[J].Science in China,2006,36:870-900.

    [7]E W N.Blow up of solutions of the unsteadly Prandtl's equations[J].Comm Pure Appl Math,1997(1):1287-1293.

    [8]XIN Z,ZHANG L.On the global existence of solutions to the Prandtl's system [J].Adv in Math,2004,191:88-133.

    [9]ZHAN H S.On a quasilinear degenerate parabolic equation [J].Chinese Annals of Math,2006,27A:731-740.

    猜你喜歡
    有界邊界條件常數(shù)
    復(fù)Banach空間的單位球上Bloch-型空間之間的有界的加權(quán)復(fù)合算子
    關(guān)于Landau常數(shù)和Euler-Mascheroni常數(shù)的漸近展開式以及Stirling級(jí)數(shù)的系數(shù)
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
    一類具低階項(xiàng)和退化強(qiáng)制的橢圓方程的有界弱解
    幾個(gè)常數(shù)項(xiàng)級(jí)數(shù)的和
    淺談?wù)?xiàng)有界周期數(shù)列的一些性質(zhì)
    萬有引力常數(shù)的測量
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    基于sub-tile的對(duì)稱有界DNA結(jié)構(gòu)自組裝及應(yīng)用
    亚洲av电影在线进入| 最近最新免费中文字幕在线| 黄色 视频免费看| 成人精品一区二区免费| 亚洲欧美激情综合另类| 视频区欧美日本亚洲| 免费看a级黄色片| 国产精品98久久久久久宅男小说| 欧美另类亚洲清纯唯美| 欧美乱妇无乱码| 我要搜黄色片| 国产主播在线观看一区二区| 日韩精品中文字幕看吧| 精品国内亚洲2022精品成人| 国产视频内射| 久久亚洲真实| 亚洲精品国产精品久久久不卡| 久久精品国产清高在天天线| 亚洲精品一区av在线观看| 一区二区三区激情视频| 看免费av毛片| 一本精品99久久精品77| 最新在线观看一区二区三区| 国产伦人伦偷精品视频| 国产精品日韩av在线免费观看| 少妇人妻一区二区三区视频| 悠悠久久av| 国产成+人综合+亚洲专区| 国产人伦9x9x在线观看| 午夜福利成人在线免费观看| 后天国语完整版免费观看| 香蕉丝袜av| 视频区欧美日本亚洲| 欧美午夜高清在线| 亚洲欧美日韩无卡精品| 成人av一区二区三区在线看| 三级男女做爰猛烈吃奶摸视频| 看片在线看免费视频| 久久久久国产一级毛片高清牌| 嫁个100分男人电影在线观看| 香蕉丝袜av| 琪琪午夜伦伦电影理论片6080| 十八禁网站免费在线| 1024香蕉在线观看| 最新中文字幕久久久久 | 成人欧美大片| 欧美中文日本在线观看视频| 黄色片一级片一级黄色片| 国产高清有码在线观看视频| 国产精品日韩av在线免费观看| 美女高潮喷水抽搐中文字幕| 国产又色又爽无遮挡免费看| 免费看a级黄色片| 男人舔女人的私密视频| 小说图片视频综合网站| 色老头精品视频在线观看| 夜夜躁狠狠躁天天躁| 中文字幕人成人乱码亚洲影| 嫩草影视91久久| 成人三级黄色视频| 日韩 欧美 亚洲 中文字幕| 男女下面进入的视频免费午夜| 精品一区二区三区视频在线观看免费| 国产黄片美女视频| 亚洲av片天天在线观看| 精品一区二区三区av网在线观看| 国产成人一区二区三区免费视频网站| 高清毛片免费观看视频网站| 免费人成视频x8x8入口观看| 91av网站免费观看| 国内毛片毛片毛片毛片毛片| 亚洲精品粉嫩美女一区| 丰满人妻一区二区三区视频av | 亚洲精品乱码久久久v下载方式 | 午夜福利视频1000在线观看| 免费大片18禁| 亚洲熟妇中文字幕五十中出| 香蕉av资源在线| 国内揄拍国产精品人妻在线| 国产野战对白在线观看| 此物有八面人人有两片| 久久久久久人人人人人| 亚洲电影在线观看av| 黑人巨大精品欧美一区二区mp4| 色老头精品视频在线观看| 久久久久免费精品人妻一区二区| 免费观看精品视频网站| 国产又黄又爽又无遮挡在线| 一卡2卡三卡四卡精品乱码亚洲| 国产单亲对白刺激| 午夜福利成人在线免费观看| 久久久久国产一级毛片高清牌| 性欧美人与动物交配| 在线看三级毛片| 无人区码免费观看不卡| 欧美色视频一区免费| 欧美成人一区二区免费高清观看 | 日韩大尺度精品在线看网址| 国产 一区 欧美 日韩| 蜜桃久久精品国产亚洲av| 久久久久国产一级毛片高清牌| 久久久久性生活片| 99在线视频只有这里精品首页| 色精品久久人妻99蜜桃| 久久中文字幕一级| 亚洲狠狠婷婷综合久久图片| 变态另类丝袜制服| 国产高清视频在线观看网站| 国产精品久久久人人做人人爽| 91av网一区二区| 巨乳人妻的诱惑在线观看| 日韩免费av在线播放| 日本免费一区二区三区高清不卡| 国内揄拍国产精品人妻在线| 国产精品女同一区二区软件 | 女警被强在线播放| 18禁黄网站禁片午夜丰满| 精品久久蜜臀av无| 成人18禁在线播放| 制服人妻中文乱码| 99热这里只有精品一区 | 国产成人av激情在线播放| 黄色丝袜av网址大全| 日本一本二区三区精品| 亚洲 欧美一区二区三区| 免费观看精品视频网站| 成人一区二区视频在线观看| 99久久综合精品五月天人人| 狂野欧美白嫩少妇大欣赏| 亚洲五月婷婷丁香| 性色avwww在线观看| 欧美xxxx黑人xx丫x性爽| 99riav亚洲国产免费| 毛片女人毛片| 丝袜人妻中文字幕| 国产真实乱freesex| 国产美女午夜福利| 免费在线观看亚洲国产| 色播亚洲综合网| 日本免费a在线| 亚洲av第一区精品v没综合| 99久久无色码亚洲精品果冻| 一夜夜www| 欧美在线一区亚洲| 麻豆av在线久日| 99国产精品一区二区蜜桃av| 久99久视频精品免费| 黄色丝袜av网址大全| 久久久久久久久中文| 国产亚洲av高清不卡| 最新美女视频免费是黄的| 色av中文字幕| 人妻丰满熟妇av一区二区三区| 国产又黄又爽又无遮挡在线| 亚洲成a人片在线一区二区| 亚洲人与动物交配视频| 精品久久久久久久久久久久久| 国产精品永久免费网站| 国产1区2区3区精品| 国产精品一区二区精品视频观看| 国产一级毛片七仙女欲春2| 不卡一级毛片| 国内少妇人妻偷人精品xxx网站 | 免费高清视频大片| 丰满人妻一区二区三区视频av | 一级黄色大片毛片| 午夜福利在线观看免费完整高清在 | 国产精品av久久久久免费| 成人特级av手机在线观看| 欧美黄色片欧美黄色片| 国产精品免费一区二区三区在线| 精品国产美女av久久久久小说| 久久99热这里只有精品18| 亚洲av美国av| 国产三级黄色录像| 久久久国产精品麻豆| 亚洲精品中文字幕一二三四区| 亚洲aⅴ乱码一区二区在线播放| 国模一区二区三区四区视频 | 激情在线观看视频在线高清| 老司机在亚洲福利影院| 在线观看免费视频日本深夜| 国产99白浆流出| 日韩人妻高清精品专区| 最好的美女福利视频网| 一级a爱片免费观看的视频| 久久这里只有精品中国| 亚洲精品色激情综合| 国产亚洲欧美在线一区二区| 国产精品精品国产色婷婷| 两个人视频免费观看高清| 亚洲精品久久国产高清桃花| 1024手机看黄色片| 精品一区二区三区四区五区乱码| 成年女人毛片免费观看观看9| 国内久久婷婷六月综合欲色啪| 亚洲天堂国产精品一区在线| 欧美最黄视频在线播放免费| 成人精品一区二区免费| 国产熟女xx| 亚洲av五月六月丁香网| 成人特级黄色片久久久久久久| 色哟哟哟哟哟哟| 国产精品影院久久| 国产一区在线观看成人免费| 一进一出好大好爽视频| 国产欧美日韩精品亚洲av| 国产精品爽爽va在线观看网站| 观看美女的网站| 久久亚洲真实| 99精品久久久久人妻精品| www国产在线视频色| 久久久久精品国产欧美久久久| 脱女人内裤的视频| 黑人操中国人逼视频| 欧美不卡视频在线免费观看| 国内精品美女久久久久久| 女人高潮潮喷娇喘18禁视频| 一区福利在线观看| 国产精品久久久久久精品电影| 中文字幕人妻丝袜一区二区| 嫩草影院精品99| 啪啪无遮挡十八禁网站| 午夜精品在线福利| 色视频www国产| 丰满人妻熟妇乱又伦精品不卡| 亚洲 欧美一区二区三区| 特级一级黄色大片| 国产亚洲精品久久久久久毛片| 18禁观看日本| 少妇丰满av| 亚洲精品国产精品久久久不卡| 久久中文字幕人妻熟女| 中亚洲国语对白在线视频| 久久这里只有精品19| 欧洲精品卡2卡3卡4卡5卡区| 国产精品电影一区二区三区| 国产免费男女视频| 色哟哟哟哟哟哟| 搡老岳熟女国产| 日韩欧美三级三区| 制服丝袜大香蕉在线| 免费搜索国产男女视频| 亚洲人成电影免费在线| 亚洲一区高清亚洲精品| 久久热在线av| 亚洲欧美日韩东京热| 人妻夜夜爽99麻豆av| 波多野结衣高清作品| а√天堂www在线а√下载| 天天躁日日操中文字幕| 亚洲成人久久性| 亚洲国产精品999在线| 久久天堂一区二区三区四区| 精品乱码久久久久久99久播| 亚洲性夜色夜夜综合| 一本久久中文字幕| 亚洲精品456在线播放app | 最新在线观看一区二区三区| 成人欧美大片| 色哟哟哟哟哟哟| 国产三级中文精品| 1024香蕉在线观看| 精品国产乱码久久久久久男人| 高清毛片免费观看视频网站| 午夜成年电影在线免费观看| 成人永久免费在线观看视频| 亚洲 国产 在线| 久久精品影院6| 在线观看午夜福利视频| 99热6这里只有精品| 亚洲欧美精品综合久久99| 国产日本99.免费观看| 欧美中文综合在线视频| 香蕉久久夜色| 亚洲午夜精品一区,二区,三区| 三级男女做爰猛烈吃奶摸视频| 成年女人看的毛片在线观看| 身体一侧抽搐| 国产免费av片在线观看野外av| 51午夜福利影视在线观看| 久久午夜亚洲精品久久| 成人亚洲精品av一区二区| 9191精品国产免费久久| 国产高潮美女av| 黄色女人牲交| 国产亚洲欧美在线一区二区| a级毛片在线看网站| 国产成人av教育| 成人亚洲精品av一区二区| 国产成人精品久久二区二区免费| 又紧又爽又黄一区二区| 午夜福利18| 老汉色av国产亚洲站长工具| 俄罗斯特黄特色一大片| 美女cb高潮喷水在线观看 | 精品国产超薄肉色丝袜足j| 亚洲国产精品999在线| 久久香蕉国产精品| 日韩高清综合在线| 婷婷丁香在线五月| 一级毛片精品| 老汉色∧v一级毛片| 国产91精品成人一区二区三区| 最近最新中文字幕大全免费视频| x7x7x7水蜜桃| 日本一本二区三区精品| 成人国产综合亚洲| 午夜免费成人在线视频| 床上黄色一级片| 色播亚洲综合网| 欧美成人免费av一区二区三区| 小蜜桃在线观看免费完整版高清| 一区二区三区激情视频| 免费av毛片视频| 国产欧美日韩精品亚洲av| 亚洲五月天丁香| 欧美日韩亚洲国产一区二区在线观看| 国产毛片a区久久久久| 黄色成人免费大全| 在线视频色国产色| 久久久国产精品麻豆| 欧美日韩一级在线毛片| 啦啦啦韩国在线观看视频| 午夜福利在线观看免费完整高清在 | 中文字幕久久专区| 在线观看免费视频日本深夜| 五月伊人婷婷丁香| 国产一区二区三区在线臀色熟女| 一本一本综合久久| 亚洲在线观看片| 免费av不卡在线播放| 2021天堂中文幕一二区在线观| 亚洲国产欧洲综合997久久,| 99久久无色码亚洲精品果冻| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 亚洲国产高清在线一区二区三| 99热只有精品国产| 这个男人来自地球电影免费观看| 久久伊人香网站| 亚洲色图av天堂| 人人妻,人人澡人人爽秒播| 国产高清视频在线观看网站| 好男人在线观看高清免费视频| 日韩大尺度精品在线看网址| 黄色 视频免费看| 天堂√8在线中文| 国产淫片久久久久久久久 | 国产亚洲精品一区二区www| 亚洲中文字幕一区二区三区有码在线看 | 哪里可以看免费的av片| 国产真实乱freesex| 香蕉久久夜色| 免费av毛片视频| 老鸭窝网址在线观看| 欧美另类亚洲清纯唯美| 全区人妻精品视频| 国产私拍福利视频在线观看| 中国美女看黄片| 人人妻,人人澡人人爽秒播| 美女高潮喷水抽搐中文字幕| 免费观看的影片在线观看| 成人18禁在线播放| 在线观看免费视频日本深夜| 色综合亚洲欧美另类图片| 国产精品 国内视频| 国产久久久一区二区三区| 国产乱人伦免费视频| 久久天堂一区二区三区四区| 免费看十八禁软件| 国产人伦9x9x在线观看| 国产精品一区二区精品视频观看| 国产高清激情床上av| 中亚洲国语对白在线视频| 成人18禁在线播放| 精品不卡国产一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲五月天丁香| 在线播放国产精品三级| 国产视频一区二区在线看| 黄片小视频在线播放| 青草久久国产| 国产在线精品亚洲第一网站| 国产一级毛片七仙女欲春2| 国产精品一及| 毛片女人毛片| 夜夜躁狠狠躁天天躁| 亚洲成人久久性| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区四区五区乱码| 99久久99久久久精品蜜桃| 女人高潮潮喷娇喘18禁视频| 免费看美女性在线毛片视频| 欧美中文日本在线观看视频| 欧美激情久久久久久爽电影| 午夜福利欧美成人| 在线观看午夜福利视频| 国产 一区 欧美 日韩| 精品国产亚洲在线| 日韩三级视频一区二区三区| 国产美女午夜福利| 国产伦一二天堂av在线观看| 18禁黄网站禁片免费观看直播| 成人无遮挡网站| 99精品在免费线老司机午夜| 精华霜和精华液先用哪个| 又大又爽又粗| 亚洲在线自拍视频| 国产av在哪里看| 最近最新中文字幕大全免费视频| 最新美女视频免费是黄的| 成人av一区二区三区在线看| 午夜亚洲福利在线播放| 九色国产91popny在线| 日韩人妻高清精品专区| 欧美日韩乱码在线| 亚洲第一电影网av| av在线蜜桃| 19禁男女啪啪无遮挡网站| 国产日本99.免费观看| 黄色成人免费大全| 国产一区二区激情短视频| 免费在线观看成人毛片| 成人一区二区视频在线观看| 两个人的视频大全免费| 亚洲,欧美精品.| 女警被强在线播放| 国产精品国产高清国产av| 美女黄网站色视频| 午夜激情欧美在线| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| 欧美丝袜亚洲另类 | 欧美极品一区二区三区四区| 草草在线视频免费看| 色综合婷婷激情| 亚洲国产高清在线一区二区三| 噜噜噜噜噜久久久久久91| 亚洲男人的天堂狠狠| 国产乱人视频| 日韩高清综合在线| 在线十欧美十亚洲十日本专区| 精品久久久久久久毛片微露脸| 久久人妻av系列| 草草在线视频免费看| 免费看日本二区| 久久久久久九九精品二区国产| 啪啪无遮挡十八禁网站| 精品国内亚洲2022精品成人| 97人妻精品一区二区三区麻豆| 十八禁人妻一区二区| 中文字幕av在线有码专区| av欧美777| 女警被强在线播放| 国产精品国产高清国产av| 久久精品影院6| 精品久久蜜臀av无| 人妻夜夜爽99麻豆av| 久久中文字幕一级| 最近最新中文字幕大全电影3| 我的老师免费观看完整版| 黄频高清免费视频| 两个人看的免费小视频| 人人妻,人人澡人人爽秒播| 麻豆av在线久日| 久久久久久久久中文| 九九久久精品国产亚洲av麻豆 | 最新在线观看一区二区三区| 亚洲国产高清在线一区二区三| 中文字幕熟女人妻在线| 亚洲国产精品久久男人天堂| 999久久久国产精品视频| 亚洲黑人精品在线| 天堂av国产一区二区熟女人妻| 色噜噜av男人的天堂激情| 亚洲一区高清亚洲精品| 99国产极品粉嫩在线观看| 操出白浆在线播放| 99久久久亚洲精品蜜臀av| 久久草成人影院| 搞女人的毛片| 精品国产美女av久久久久小说| 国产成人av教育| 在线a可以看的网站| 黄色视频,在线免费观看| 国产综合懂色| 精品久久久久久,| 国产淫片久久久久久久久 | 日韩欧美精品v在线| 国产亚洲精品久久久com| 久久久久免费精品人妻一区二区| 美女黄网站色视频| 少妇的丰满在线观看| 午夜福利在线观看免费完整高清在 | 小说图片视频综合网站| 噜噜噜噜噜久久久久久91| 一本精品99久久精品77| 91av网一区二区| 欧美日韩一级在线毛片| 婷婷精品国产亚洲av| 久久久久久大精品| 黄色欧美视频在线观看| 看片在线看免费视频| 热99在线观看视频| 少妇人妻精品综合一区二区| 国产精品1区2区在线观看.| 激情 狠狠 欧美| 欧美一级a爱片免费观看看| 色综合亚洲欧美另类图片| 国产69精品久久久久777片| 午夜福利网站1000一区二区三区| 天天躁日日操中文字幕| 国产伦一二天堂av在线观看| 长腿黑丝高跟| 99久久人妻综合| 国产精品无大码| 热99re8久久精品国产| 永久网站在线| 六月丁香七月| 特级一级黄色大片| 久久精品综合一区二区三区| 国产精品麻豆人妻色哟哟久久 | 亚洲色图av天堂| 精品久久久久久久末码| 日韩,欧美,国产一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 欧美xxxx性猛交bbbb| 寂寞人妻少妇视频99o| 好男人视频免费观看在线| 成人高潮视频无遮挡免费网站| 亚洲最大成人av| 秋霞伦理黄片| 亚洲精品乱久久久久久| 国产久久久一区二区三区| 久久国产乱子免费精品| a级毛色黄片| 欧美又色又爽又黄视频| 能在线免费观看的黄片| 国产成人精品久久久久久| 美女内射精品一级片tv| 伊人久久精品亚洲午夜| 成人漫画全彩无遮挡| 波多野结衣巨乳人妻| 日本猛色少妇xxxxx猛交久久| 乱码一卡2卡4卡精品| 3wmmmm亚洲av在线观看| 美女高潮的动态| 2022亚洲国产成人精品| 精品无人区乱码1区二区| 精品人妻偷拍中文字幕| 国产精品熟女久久久久浪| 小说图片视频综合网站| 国产成人a∨麻豆精品| 久久鲁丝午夜福利片| 免费大片18禁| 亚洲成av人片在线播放无| 91久久精品国产一区二区三区| 国产在线男女| 淫秽高清视频在线观看| 可以在线观看毛片的网站| 久久精品影院6| 精品久久久噜噜| 免费av毛片视频| 级片在线观看| 国产精品一及| 午夜久久久久精精品| www.色视频.com| 久久久久久久国产电影| 看非洲黑人一级黄片| 特大巨黑吊av在线直播| 亚洲成人精品中文字幕电影| 青春草国产在线视频| 亚洲精品日韩在线中文字幕| 51国产日韩欧美| 天天躁夜夜躁狠狠久久av| 亚洲国产色片| 久久久久久九九精品二区国产| 1024手机看黄色片| 午夜福利在线观看免费完整高清在| 男插女下体视频免费在线播放| 精品国产露脸久久av麻豆 | 久久精品国产自在天天线| 成人亚洲精品av一区二区| 国产 一区精品| 如何舔出高潮| 久久精品久久久久久噜噜老黄 | 只有这里有精品99| 久热久热在线精品观看| 成人国产麻豆网| 亚洲av不卡在线观看| 午夜a级毛片| 十八禁国产超污无遮挡网站| 18禁在线播放成人免费| 欧美激情国产日韩精品一区| 嫩草影院新地址| 少妇高潮的动态图| 成人无遮挡网站| 日本-黄色视频高清免费观看| 精品一区二区三区视频在线| 日韩精品青青久久久久久| 久久久久性生活片| 免费看a级黄色片| 免费看日本二区| 午夜a级毛片| 老师上课跳d突然被开到最大视频| 国产成人免费观看mmmm| 99热这里只有是精品50| 床上黄色一级片| 国产在线男女| 深爱激情五月婷婷| 久久久久国产网址| 欧美成人a在线观看| 精品不卡国产一区二区三区|