Разделы сайта

Читаемое

Обновления Aug-2018

Промышленность Ижоры -->  Динамика жидкости: уравнения 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 [ 52 ] 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

1 SUBROUTINE POINT(N,FMN)

3 с CALCULATES THE FLOW AT GIVEN POINTS, (XP,YP)

5 DIMENSION X(50)Л(50),XC(50),YC(50),DS(50),FN(50,50)

6 1,FT(50,50),RHS(50),SDE(50),CI(50),SI(50)

7 COMMON X,Y.XC.YC,DS.FN,FT,RHS,PI,CPI CI,SI

8 1,UINF,VINF,SDE

10 FAC = SQRTd. - FMN*FMN)

11 GAM = 1.4

12 CI = 0.5*(GAM-1.)*FMN*FMM

13 C2 = 0.5*GAM*FMN*FMN

14 GMP = GAM/(GAM-1.)

15 1 CONTINUE

16 READ(1,2) XP,YP

17 2 FORMAT(2F8.5)

18 YP = YP*FAC

19 RPS XP*XP + YP*YP

20 IF(RPS .LT. l.OE-04)RETURN

21 С

22 UUUINF

23 VV=VINF

24 DO 3 J=1,N

25* DYJ=SI(J)*DS(J)

26 DXJ=CI(J)*DS(J)

27 SPH=DS(J)*0.5

28 XD4XP-XC(J)

29 YD=YP-YC(J)

30 R=SQRT(XD*XD+YD*YD)

31 BET ATAN2(YD,XD.)

32 ALJ=ATAN2(DYJ,DXJ)

33 GAH=ALJ-BET

34 ZI=R*COS(GAM)

35 ET=-R*SIN(GAM)

36 R1S=((ZI+SPH)**2)+ET*ET

37 R2S ((ZI-SPH)**2)+Et*ET

38 QT=AL0G(R1S/R2S)

39 DEN=ZI*ZI+ET*ET-SPH*SPH

40 GN ET*DS(J)

41 QN 2.0*ATAN2(GN,DEN)

42 UJ=QT*CI(J)-QN*SI(J)

43 Va=QT*SI(J)+QN*CI(J)

44 UU UU-fUJ*SDE(J)

45 VV=VV+VJ*SDE(J)

46 3 CONTINUE

47 С

48 YP = YP/FAC

49 UU = UU/FAC/FAC

50 VV VV/FAC

51 PP=1.-UU*UU-VV*VV

52 IF(FMN .GT. 0.05)PP ((l.+Cl*PP)**GMP-1.)/C2

53 WRITE(6,4)XP,YP,UU,VV,PP

54 4 F0RMAT(/,2X,FL0W AT ХД= ,2Г6.3, U,V=*,2Гб.З. P=

55 1 F6.3)

56 GO TO 1

57 END

Рис. 14.11. Распечатка программы POINT.



Гл. 14. Невязкие течения

Таблица 14.2. Параметры, используемые программой PANEL

Параметр

Описание

Матрица А в (14.6)

Длина малой полуоси эллипса

Число Маха набегающего потока М

FN (К, J)

Нормальная к панели k скорость в точке (Хк,Ук), здаваемая элементом а/

РТ(К, J)

Касательная к панели k скорость в точке (Xk, создаваемая элементом о/

>0, печатается X, У, ХС, ГС, Z)S, С/, SI в BODY =0, печатается FN, FT, RHS в MATELM

Число конечных точек элементов

Число элементов

Вектор R в (14.6)

Вектор плотностей источников а

UINF, VINF

Компоненты скорости набегающего потока f/oe и Voo

X, Y

Координаты конечных точек панелей

ХС, YC

Координаты контрольных точек панелей

CI(J), SI(J)

cos а/, sin а/ на рис. 14.5

DS(J)

Длина панели, ds/

а/, рис. 14.5, MATELM и POINT

kj. рис. 14.5, MATELM; BET в POINT

YAr/, рис. 14.5, MATELM; GAM в POINT

r/, рис. 14.5, MATELM; R в POINT

QN, QT

qlj; q\j в (14.13), MATELM и POINT

QNK, QTK

Компоненты скорости, нормальные и касательные к панели в точке (Xk, ук)

UKJ, VKJ

Ukb ki в (14.12), MATELM; W и У/ в POINT

UU, VV, PP

и, Vy -р в точке (Xk, Ук) в SURVL и в точке {Хр, в POINT

Точное значение касательной составляющей скорости на поверхности эллипса, gt, ex

сравнивается на поверхности эллипса с QEX (рис. 14.12), точным значением тангенциальной компоненты,

.. (1 + 6) у

(14.16)



panel method vith 20 elements, ellipse minor semi-axis .500

onset velocity components = 1.000 .000 freestream mach number = .000

velocity and pressure at the control points

xc,yc= -.976 .077 qn,qt= .000 .4 u,v= .135 .427 p= .800 qex* .453

xc,yc -.880 .224 qn,qt .000 1.067 u,v= .762 .748 p -.139 qex 1.071

xc,yc -.698 .349 qn,qt= .000 1.342 u,v= 1.200 .600 p* -.801 qex 1.342

xc,yc -.448 .440 qn,qt= .000 1.455 u,v= 1.410 .359 p=-1.117 qex= 1.454

xc,yc -.155 .488 qn,qt= .000 1.497 u,v= 1.492 .118 p=-1.240 qex 1.495

xc,yc= .155 .488 qn,qt= .000 1.497 u,v= 1.492 -.118 p=:-1.240 qex= 1.495

xc,yc= .448 .440 qn,qt= .000 1.455 u,v= 1.410 -.359 p=-1.117 qex 1.454

xc,yc= .698 .349 qn,qt=: .000 1.342 u,v 1.200 -.600 p= -.801 qex 1.342

xc,yc= .880 .224 qn,qt= .000 1.067 u,v= .762 -.748 p= -.139 qex* 1.071

xc,yc= .976 .077 qn,qt .000 .448 u,v .135 -.427 p= .800 qex .453

xc.yc= .976 -.077 qn,qt= .000 -.448 u,v= .135 .427 p= .800 qex= -.453

xc,yc .880 -.224 qn,qt= .000-1.067 u,v= .762 .748 p -.139 qex=-1.071

xc,yc .698 -.349 qn,qt= .000-1.342 u,v= 1.200 .600 p= -.801 qex -1.342

xc,yc .448 -.440 qn,qt .000-1.455 u,v= 1.410 .359 p=-1.117 qex -1.454

xc.yc* .155 -.488 qn,qt= .000-1.497 u,v= 1.492 .118 p=-1.240 qex=-1.495

xc,yc -.155 -.488 qn,qt= .000-1.497 u,v= 1.492 -.118 p=-1.240 qex -1.495

xc,yc= -.448 -.440 qn,qt .000-1.455 u,v- 1.410 -.359 p -1.117 qex -1.454

xc,yc= -.698 -.349 qn,qt .000-1.342 u,v= 1.200 -.600 p -.801 qex=-1.342

xc.yc= -.880 -.224 qn,qt .000-1.067 u,v= .762 -.748 p -.139 qex -1.071

xc,yc= -.976 -.077 qn,qt .000 -.448 u,v= .135 -.427 p= .800 qex* -.453

flow at x,y= .000 1.000 u,v= 1.257 .000 p= -.579

Рис. 14.12. Типичная выдача программы PANEL.

Путем введения большого числа панелей, в особенности в окрестности носка (YC 0) и миделя (УС 0.5), можно добиться лучшего соответствия. Кроме того, на рис. 14.12 приведены декартовы компоненты скорости и давления в контрольных точках и в точке (О, 1.0), лежащей вне тела.

14,1.3, Связь с методом граничных элементов

Панельный метод, описанный в п. 14.1.1, особенно эффективен при расчете обтекания однородным потоком изолированных тел. Однако для внутренних течений около препятствия в канале часто оказывается более удобной иная формулировка, основанная на теореме Грина. При таком подходе решение Ф находится непосредственно, без введения промежуточного распределения панельных источников.

Для двумерных задач потенциал в любой точке (х/г, ук) области самой общей формы может быть связан со значениями Ф и дФ/дп на расчетной границе S:

1 Г f /1 дФ , . . с д

(In г,/) {s)ds\ (In r,i) Ф {s) ds

(14.17)



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 [ 52 ] 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182

© 2003 - 2018 Prom Izhora
При копировании текстов приветствуется обратная ссылка
Besucherzahler
счетчик посещений