Квадрат саад тойрч урсах урсгал
Латтис Больцманы
арга нь мөргөлдөөний интегралын ойролцооллоос хамаарч хэд хэдэн төрөлд
хуваагдах ба эдгээрт тайвшралтын хугацаа чухал үүрэг гүйцэтгэдэг. Хамгийн
энгийн мөргөлдөөний операторын ойролцоолол нь ВГК буюу нэг тайвшралтын
хугацаатай мөргөлдөөний илэрхийлэл юм. Үүнээс гадна хоёр тайвшралтын хугацаатай
болон олон тайвшралтын хугацаатай гэж байна. Эдгээр нь ганц тайвшралтын
хугацаатай мөргөлдөөний интегралаас боловсруулахад төвөгтэй боловч тогтвортой
шийдийг өгч чаддаг. ВГК ойролцоололд Больцманы тэгшитгэлийг шийдэх бас олон
янзын зам байна. Олон фазат болон олон бүрэлдэхүүнт урсгал бүрд тусгай тусгай
хандлагууд судлагдсан байна. Мөн Махын тооны хязгаарлалтаас болж турбулент
урсгалыг загварчлах боломж ердийн ВГК-д байхгүй болох ба зарим нэг сайжруулалт
хийгдэх шаардлагатай. Үүнээс гадна тэнцвэрт түгэлтийн функцад нягтыг хаалтанд
байх задаргааны гишүүд дээр шууд нэмж загварчлах боломж байна. Ердийн ВГК болон
нягтын функцыг дотор агуулсан ВГК хоёрын ялгааг квадрат саадыг тойрч урсах
урсгал дээр туршиж Рейнольдсын тооны боломжит утгад хүрье.
Мөргөлдөх
процессыг тооцоолох
Ердийн
|
Агуулсан
|
DO i=0,n
DO j=0,m
t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)
DO k=0,8
t2=u(i,j)*cx(k)+v(i,j)*cy(k)
feq(k,i,j)=w(k)*rho(i,j)*(1.+3.0*t2+4.50*t2*t2-1.50*t1)
f(k,i,j)=omega*feq(k,i,j)+(1.-omega)*f(k,i,j)
END DO
END DO
END DO
|
DO i=0,n
DO j=0,m
t1=u(i,j)*u(i,j)+v(i,j)*v(i,j)
DO k=0,8
t2=u(i,j)*cx(k)+v(i,j)*cy(k)
feq(k,i,j)=w(k)* (rho(i,j)+3.0*t2+4.50*t2*t2-1.50*t1)
f(k,i,j)=omega*feq(k,i,j)+(1.-omega)*f(k,i,j)
END DO
END DO
END DO
|
Харин хурдыг
тооцоолоход
Ердийн
|
Агуулсан
|
DO i=1,n
DO j=1,m-1
usum=0
vsum=0
DO k=0,8
usum=usum+f(k,i,j)*cx(k)
vsum=vsum+f(k,i,j)*cy(k)
END DO
u(i,j)=usum/rho(i,j)
v(i,j)=vsum/rho(i,j)
END DO
END DO
|
DO i=1,n
DO j=1,m-1
usum=0
vsum=0
DO k=0,8
usum=usum+f(k,i,j)*cx(k)
vsum=vsum+f(k,i,j)*cy(k)
END DO
u(i,j)=usum
v(i,j)=vsum
END DO
END DO
|
Тайвшралтын
параметрийг гараас өгч турбулент горимд хурэхийг хичээе. Эхний утга нь анхны
хурд юм. Дараагийн утга нь тайвшралтын параметр юм.
Ердийн ВГК загвар
0.1
1.9
Re=
1139.999 ome 1.900000
tau 0.5263158
Fortran
Pause - Enter command<CR> or <CR> to continue.
Урсгалын шугам |