亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Benard對流問題五模類Lorenz方程組混沌行為的數(shù)值模擬

        2019-12-03 08:20:10王賀元
        關(guān)鍵詞:軌線板間熱傳導(dǎo)

        王賀元, 張 穎

        (沈陽師范大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 沈陽 110034)

        0 引 言

        Lorenz方程是將滿足一定邊界條件的Navier-Stokes方程和熱傳導(dǎo)方程進(jìn)行傅立葉展開后截斷而產(chǎn)生的[1-2]。其實際背景來源于簡化的大氣對流模型,兩平行板間充滿流體從下板加熱,熱量從下至上傳遞,當(dāng)溫度比較低時,熱量通過熱傳導(dǎo)方式傳遞,流體不發(fā)生運動。當(dāng)溫度增高時,液體就要發(fā)生對流,當(dāng)溫度繼續(xù)升高時,液體會發(fā)生對流。在熱壓力條件下不可壓縮的流體運動可以被描述為如下偏微分方程組:

        (1)

        其中u=u(x,y,z)為流體的速度場,T=T(x,y,z)為流體的溫度場,在這個方程中式中常數(shù)g,ε,v和k分別表示重力加速度,熱膨脹系數(shù),運動粘性系數(shù)和熱傳導(dǎo)系數(shù),ΔT為兩板間的溫度差,p為流體的壓力場,在邊界上u=u(x,y,z)。

        1 豎直截面上的二維對流問題

        僅考慮豎直截面上二維流動問題,引入函數(shù)τ(x,z,t),它的梯度為速度場,同時引入θ(x,z,t)為流體的溫度場,并且在靜態(tài)條件下θ(x,z,t)=T(x,z,t)-T0,這里的T0在兩板間呈線性遞減,假定流體不可壓縮,則方程(1)可表示為如下形式:

        函數(shù)ψ(x,z,t)和θ(x,z,t)在邊界上滿足

        (4)

        對ψ(x,z,t)和θ(x,z,t)傅立葉展開成如下形式:

        (5)

        在具體截斷過程中,將m1m2,n1,n2,賦予具體的值,使得ψ(x,z,t)和θ(x,z,t)成為具體的式子,然后將ψ(x,z,t)和θ(x,z,t)代入式(2)~式(4)中,經(jīng)過復(fù)雜計算通過待定系數(shù)得到截斷方程。

        依照上述截斷的方法,在截取五模類Lorenz方程組時,將ψ(x,z,t)和θ(x,z,t)展式設(shè)為如下形式:

        (6)

        其中x1相當(dāng)于通項中的A1,x2相當(dāng)于通項中的A2,x3相當(dāng)于通項中的A11,x4相當(dāng)于通項中的A12,x5相當(dāng)于通項中的B02,C1,C2為常系數(shù)。

        經(jīng)運算得

        把上述各式代入式(2)和式(3)中,再利用邊界條件(4)經(jīng)整理并根據(jù)物理意義取定其中的參數(shù)得到如下五模類Lorenz方程組:

        (14)

        其中r與溫度有關(guān),稱為Rayleigh數(shù)。

        2 數(shù)值模擬

        上述非線性方程組(14)具有復(fù)雜的動力學(xué)行為,下面我們來數(shù)值模擬當(dāng)參數(shù)r變化時方程組的動力學(xué)行為如下(取坐標(biāo)為x2,x4,x5):

        圖1 r=48.3 繞一點旋轉(zhuǎn)Fig.1 r=48.3 Rotate about a point

        1) 圖1表示當(dāng)r=48.3時方程組的解軌線不斷繞一點(平衡點)旋轉(zhuǎn),且越轉(zhuǎn)越密,此時系統(tǒng)穩(wěn)定,圖2(r=69.2)為軌線繞一點旋轉(zhuǎn)的另外一種形態(tài),軌線不斷增多且向一起靠攏;圖3表示當(dāng)r=71.4時,系統(tǒng)發(fā)生了Hopf分岔[3-6],出現(xiàn)了極限環(huán);圖4所示當(dāng)r=72.3時由原先的繞一點旋轉(zhuǎn)變?yōu)槔@兩點旋轉(zhuǎn)。

        2) 圖5(r=74)、圖6(r=76)為系統(tǒng)不斷分岔出一些新軌線,新軌線繼續(xù)繞兩點旋轉(zhuǎn),即發(fā)生了混沌現(xiàn)象,出現(xiàn)奇怪吸引子[6]。

        圖2 r=69.2軌線聚集Fig.2 r=69.2 Rail line to gather

        圖4 r=72.3繞2點旋轉(zhuǎn)Fig.4 r=72.3 Rotation about two points

        圖5 r=74奇怪吸引子 Fig.5 r=74 Strange attractor

        3) 圖7~圖11分別給出了系統(tǒng)的分岔圖、最大李雅普諾夫指數(shù)、龐加萊截面、返回映射和功率譜,它們展現(xiàn)了系統(tǒng)混沌行為的普適特征。

        圖6 r=76奇怪吸引子Fig.6 r=76 Strange attractor

        圖7 分岔圖Fig.7 Bifurcation diagram

        圖8 最大李雅普諾夫指數(shù)Fig.8 The largest lyapunov index

        圖9 r=84.65龐加萊截面Fig.9 r=84.65 Poincare spectru

        圖10 r=84.65返回映射Fig.10 r=84.65 Return to the map

        圖11 r=84.65功率譜Fig.11 r=84.65 Power spectrum

        3 結(jié) 論

        本文介紹了Lorenz方程截斷的基本方法,對Navier-Stokes方程與熱傳導(dǎo)方程中的速度場和溫度場等變量進(jìn)行二維傅立葉展開,經(jīng)過復(fù)雜運算最后得到新五模方程組,并對其動力學(xué)行為進(jìn)行了數(shù)值模擬。

        猜你喜歡
        軌線板間熱傳導(dǎo)
        時變切換時滯反饋鎮(zhèn)定混沌系統(tǒng)不穩(wěn)定周期軌線1)
        一類三維逆時熱傳導(dǎo)問題的數(shù)值求解
        關(guān)于二維線性自治系統(tǒng)的相圖的教學(xué)
        一類反應(yīng)-擴(kuò)散-對流方程的平衡解分析
        費米氣體光晶格模型的漸近軌線
        電容器的兩種工作狀態(tài)
        關(guān)于電容器中U不變還是Q不變的判斷
        關(guān)于電容器中U不變還是Q不變的判斷
        電容器的兩種工作狀態(tài)
        熱傳導(dǎo)方程解的部分Schauder估計
        亚洲国产av剧一区二区三区| 天天夜碰日日摸日日澡| 人妻丰满熟妇av无码处处不卡| 国产一级在线现免费观看| 加勒比特在线视频播放| 亚洲女人毛茸茸粉红大阴户传播 | 国产精品午夜福利亚洲综合网 | 女人色熟女乱| 国产精品白浆一区二小说| 免费大学生国产在线观看p| 一区二区三区日韩蜜桃| 国产综合色在线精品| 蜜桃成人无码区免费视频网站| 午夜影视啪啪免费体验区入口| 国产午夜激情视频在线看| 曰韩无码av一区二区免费| 久久av高潮av无码av喷吹| 欧美日韩中文亚洲另类春色| 国内精品少妇久久精品| 色偷偷色噜噜狠狠网站30根| 久久伊人色av天堂九九| 日韩精品久久久中文字幕人妻 | 国内精品九九久久精品小草| 最新中文字幕日韩精品| 99视频30精品视频在线观看| 国产爽爽视频在线| 自拍视频在线观看成人| 亚洲视频高清一区二区| 国产98在线 | 日韩| 成人片99久久精品国产桃花岛| 一区二区日本免费观看| 乱色欧美激惰| 粗一硬一长一进一爽一a级| 中文字幕在线一区乱码| 日本熟女中文字幕在线| a级毛片高清免费视频就| 亚洲av成人一区二区三区网址| 日本女优久久精品久久| 人妻少妇中文字幕乱码| 97久久精品人人妻人人| 日本在线观看一区二区三区视频 |