Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения 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.
сравнивается на поверхности эллипса с 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)
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |