Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения ЛЛу,у 1 = -(1 -Uy)- l -W АА 2 I ЛЛ/,у+1= -(1 - гу)- 1 -и (15.69) У + 1 Уравнения (15.67) - (15.69) получены в результате аппроксимации (15.63) одномерными линейными элементами при разбиении на Аи. На стенке и внешней границе пограничного слоя Таблица 15.4. Параметры, используемые в программе TURVS
имеется лишь один элемент и и/ = О и 1 соответственно. Соответствующие формулы используются в подпрограмме FEl-AR (рис. 15.15). В начале каждого шага вниз по потоку величины и, Kil и uju!l интерполируются по вводимым значениям Ue{l) и dueldl, в подпрограмме ABCD (рис. 15.16) определяются коэффициенты трехдиагональной системы (15.66), а в подпрограммах BANFAC и BANSOL (п. 6.2.3) осуществляется ее решение. В результате определяются значения А9 + SUBROUTINE TURVS(REL,DELTA,IMAX,UE,DUEDX,TAU,TRV,DTRV,АКР, 1APZ,ATR,PC0N) CALCULATE EDDY VISCOSITY, TRV AND DTRV DIMENSION TRV(41),YP{41),TAU(41),DTRV(41),UP(4i) TRVO IS THE EDDY VISCOSITY IN THE OUTER REGION(CLAUSER FORMULATION) RSQ = SQRT{R2L) REQ = SQRT(RSQ) TRVO = ATR*UE*REL*DELTA IMAP = IMAX - 1 AIM = IMAP DU = l./AIM YPd) = 0. UP(1) = 0. UT = SQRT(TAU(1))/REQ PP = -DUEDX/UE/REL/UT/UT/UT EN = 1. + PCON*PP AP = APZ/EN RUT = UT*RSQ UP{2) = DU/UT YP(2) = UP(2) TRV(l) = 0. DO 2 К = 2,IMAP AK = К - 1 и = AK*DU IF(K .EQ. 2)G0 TO 1 UP(K) = U/UT их = l./(l. - и + 2.*DU) UY = l./(l. - и + DU) UZ = l./(l. - U) DEM = DU*(UX/TAU(K-2) + 4.*UY/TAU(K-l) + UZ/TAU(K))/3. YP(K) = YP(K-2) + DEM*RUT 1 DUM = 1. - EXP(-YP(K)/AP) EL = AKP*YP(K)*DUM TRV(K) = EL*EL*TAU(K)M1.-U)/TAU(1) IF(TRV(K) .GT. TRVO)GOTO 3 DTRV(K} = -(l.-U)*EL*EL*TAU(K)*TAU(K)/TAU(l) 2 CONTINUE 3 DO 4 L = K,IMAX DTRV(L) = 0. TRV(L) = TRVO 4 CONTINUE RETURN END Рис. 15.14. Распечатка программы TURVS. Прежде чем получить новое решение 9+ = в AQ +\ рассматривается необходимость изменения шага маршевой переменной Ag в зависимости от значения АбУУй. Если А9 79> > Уу шаг Ag уменьшается вдвое при условии, что полученный шаг будет больше минимального значения Ag. Если А уменьшается наполовину, и т. д. пересчитываются и рассчитывается новое решение А9 + Обычное значение у порядка 0.02. 1 SUBROUTINE FEPAR(IMAX) 3 с COMPUTES VALUES OF СС,АА, AND ЕГ ACROSS BOUNDARY LAYER {15.3ЛЗ) 5 DIMENSION CC(41,5),AA(41,5),EF(4X,5),ADC(5,65),D(65) 6 1,TAU(41),TRV(41),DTRV(41) 7 COMMON CCДА,EF,ABC,D,TAU,TRV,DTRV 8 SI = 1./6. 9 TI = 1./3. 10 TTI = 2.*TI 11 IMAP = IMAX - 1 12 AIM = IMAP 13 DU = l./AIM 14 DUS * DU*DU 15 с 16 с AT WALL 17 С 18 CC(1,1) = 0. 19 CC(1,2) DU/12. 20 CC(1.3) DU/12. 21 AA(1,1) = 0. 22 AA(1,2) = l./DUS + TI 23 AA(1,3) = -AA(1,2) + l./DU 24 EF(1,1) = 0. 25 EF(1,2) = -0.5/DU - TI 26 EF(1,3) = -0.5/DU - SI*(1.-DU) 27 С 28 с AT OUTER EDGE OF B.L. 29 с 30 CC(IMAX,1) = SI*(1. - 0.5*DU) 31 CC(IMAX,2) = TI*(1. - 0.25*DU) 32 CC(IMAX,3) = 0. 33 AA(IMAX,1) = -TI 34 AA(IMAX,2) = TI 35 AA(IMAX,3) 0. 36 EF(IMAX,1) = TI*(1. - 0.5*DU) 37 EF(IMAX,2) = -TI 38 EF(IMAX,3) = 0. 39 с 40 С AT INTERIOR NODES 41 С 42 DO 1 J = 2,IMAP 43 AJ = J-l 44 и = AJ*DU 45 CC(J,1) = SI*(U - 0.5*DU) 46 CC(J,2) = TTI*U 47 CC(J,3) * SIMU + 0.5*DU) 48 DUM = (1. - U)/DU 49 DUS = DUM*DUM 50 AA(J,1) = -DUS - DUM - TI 51 AA(J,2) = 2.*DUS + TTI 52 AA(J,3) = - DUS + DUM - TI 53 DUM = 0.5*(1. - U*U)/DU 54 DUS = 0.5*U - SI 55 EF(J,1) = DUM + DUS - SI*DU 56 EF(J,2) = - TTI 57 EF(J,3) = - DUM + DUS + SI*DU 58 1 CONTINUE 59 RETURN 60 END
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |