Метод расчета местных характеристик турбулентного теплопереноса крупномасштабными структурами в автомодельной части свободных сдвиговых течений
Диссертация
Представленный в диссертационной работе метод расчета местных характеристик процесса турбулентного переноса импульса и тепла в свободных сдвиговых течениях позволяет учесть и численно исследовать развитие в моногармоническом приближении как крупно-, так и мелкомасштабных пульсаций. Методом расчета можно определить зависимость линейных размеров, частоты, распределения по сечению среднего течения… Читать ещё >
Список литературы
- Петухов B.C. Турбулентность в теории теплообмена. — В кн.: Тепломассообмен- 1. ч. I. Минск, ИТМО АН БССР, 1981, с. 20−51.
- Турбулентность. Современные проблемы прикладной физики/Под ред. Брэдшоу П. и др. М.: Машиностроение, 1980, с. 431.
- Karraan Th.W. fiber laminare und turbulente Riebung. Z. angew. Math. Mech., 1921, Bd. 1, H. 4, S. 233−252.
- Кутателадзе С.С., Леонтьев А. И. Тепломассообмен в турбулентном пограничном слое. М.: Энергия, 1972. — 342 с.
- Tennekes Н., Lurrley J.L. A first cause in turbulence. Gam-bridge: MIT Press, 1972. — 300 p.
- Corrsin S. Advances in geophysics. New-York — San Fransisco- London, Acad. Press, 1974, S8A. p.
- Launder B.B., Spalding D.B. Mathematical models of turbulence.- London: Acad. Press, 1972. 170 p.
- Коловандин Б.А. Моделирование теплопереноса при неоднородной турбулентности. Минск: Наука и техника, 1980. — 183 с.
- Daly B.J., Harlow F.H. Transport equations in turbulence. -Phys. Fluids, 1970, vol. 13, N 11, p. c634−2649.
- Монин А.С., Яглом A.M. Статистическая гидромеханика. 4.1. -М.: Физматгиз, 1965. 639 с.
- Patel V.C., Scheuerer G. Calculation of two-dimensional near and far wakes. AIAA J., 1982, vol. 20, N 7, p. 900−908.
- Ramaprian B.R., Patel V.C. The symmetric turbulent wake of a flat plate. AIAA J., 1982, vol. cO, N 9, p. 1228−1235.
- Ferziger J.H. Turbulent flow simulation: a large eddy simulators viewpoint. Rec. Contr. Fluid Mech., Berlin, 1982, p. 69−72.•14. Таунсенд А. А. Структура турбулентного потока с поперечным сдвигом. М.: ИИЛ, 1959. — 399 с.
- Landahl М.Т. A wave-guide model for turbulent shear flbw.- J. Fluid Mech., 1567, vol. 30, p. 441−459.
- Молло-Кристенсен E. Физика турбулентных течений. Ракетная техника и космонавтика, 1971, т. 9, № 7, с. 3−16.
- Davies P.O.L., Yule A.J. Coherent structures in turbulence.- J. Fluid Mech., 1975, vol. 69, part. 3, p. 513−537.
- Blackwelder R.F. Dynamic measurements in unsteady flows. -In: Proc. Dynamic Flow Gonf. 1978, Marseilles and Baltimore, 1979, P. 173−187.
- Mumford J.C. The structure of the large eddies in fully developed turbulent shear flows. Part 1. The plane jet. J. Fluid Mech., 1982, vol. 118, p. 241−268.
- Hussain A.K.M.I'. Coherent structures-reality and myth.-Phys. Fluids, 1983, vol. 26, N 10, p. 2816−2850.
- Yule A.J. Large-scale struture in the mixing layer of a round jet. J. Fluid Mech., 1978, vol. 89, p. 413−433.
- Hernan M.A., Jimenez J. Computer analysis of a high-speed film of the plane turbulent mixing layer. J. Fluid Mech., 1982, vol. 119, P. 323−345.
- Zaman K.B.M.Q., Hussain A.K.M.F. Vortex pairing in a circular jet under controlled excitation. Part 1. General jet response. J. Fluid Mech., 1980, vol. 101, part 3, p. 449−491.
- Goldschmidt V.W., Young M.F., Ott E.S. Turbulent convective velocities (broadband and wavenumber dependent) in a plane jet. J. Fluid Mech., 1981, vol. 105, p. 327−345.
- Hussain A.K.M.F., Thompson C.A. Controlled symmetric perturbation of the plane jet- an experimental study in the initial region. J. Fluid Mech., 1980, vol. 100, part 2, p. 397 431.
- Hussain A.K.M.F., Clark A.R. On the coherent structure of the axisymmetric mixing layer- a flow-visualization study, J. Fluid Mech., 1981, vol. 104, p. 263−294.
- Hussain A.K.M.F., Zaman K.B.M.Q. The «preferred mode» of the axisymmetric jet. J. Fluid I/lech., 1981, vol. 110, p. 39−71.
- Launder 3.E. Heat and mass transport. In: Heat and mass transport topics in Applied Physics. Berlin, Springer, 1976, vol. 12, p. 231−287.
- Morel R., Key C., Wallace J.M. Evolution of the thermal field in a turbulent wake dounstream from an asymmetrically heated plate. Phys. Fluids, 1983, vol. 26, N 2, p. 416−421.
- Gence J.N., Mathieu J. On the application of successive plane strains to grid-generated turbulence. J. Fluid Mech., 1979, vol. 93, part 3, p. 501−513.
- Gence J.H., Mathieu J. The return to isotropy of an homogeneous turbulence having been submitted to two successive plane strains. J. Pluid Mech., 1980, vol. 101, part. 3, p. 555−566.
- Шарне Ж., Шон Ж.P., Алькарац Е. Тепловые характеристики турбулентного пограничного слоя с изменением знака потока теплана стенке. В кн.: Турбулентные сдвиговые течения I/Под ред. Гиневского А. С. — М.: Машиностроение, 1980. — 431 с.
- Kurbatskii А.P., Yanenko N.N. On modelling of effects of negative production of tempereture fluctuation intensity in the turbulent mixing layer. — J. Pluid Mech., 1983, vol. 13 0, p. 453−462.
- Sreenivasan K.R., Antonia R.A., Britz D. Local izotropy and large structures in a heated turbulent jet. J. Fluid Mech., 1979, vol. 94, part 4, p. 745−775.
- Rajagopalan S., Antonia R. A. Properties <5f the large structure in a slightly heated turbulent mixing layer of a planejet. J. Fluid Mech., 1981, vol. 105, p. 201−281.
- Reynolds A.J. The prediction of turbulent Prandtl and Schmidt numbers. Intern. J. Heat Mass Trans., 1975, vol. 18, N 9, p. 1055−1069.
- Иевлев B.M. Турбулентное движение высокотемпературных сплошных сред. М.: Наука, 1975. — 256 с.
- Петухов Б.С., Медвецкая Н. В. Турбулентное течение и теплообмен. в вертикальных трубах при сильном влиянии подъемных сил. -ТВТ, 1978, т. 16, № 4, с. 778−801.
- LuraLey J.L., Zeman 0., Siess J. The influence of buoyancy on turbulent transport. J. Fluid Mech., Mech., 1978, vol. 84, part 3, p. 531−597.
- Лаундер Б.Е. Модели замыкания для напряжений третье поколение. — В кн.: Турбулентные сдвиговые течения I/Под ред. Гиневского А. С. — М.: Машиностроение, 1980, — 431 с.
- Дайковский А.Г., Полежаев В. И., Федосеев А. И. Численное моделирование переходного и турбулентного режимов конвекции на основе нестационарных уравнений Навье-Стокса. Препринт ИПМ АН СССР, 1978, № 101 — 65 с.
- Полежаев В.И. Численное моделирование турбулентных течений жидкости. В кн.: Численные методы решения задач переноса. Материалы Междунар. школы — семинара. 4.2. Минск, 1979, с. 86−92.
- Orszag S.A., Izraeli М. Numerical simulation of viscous incompressible flows.-Ann.Rev.Fluid Mech., 1974, vol.6, p.281−318.
- Гретцбах Г., Шуман У. Прямое численное моделирование турбулентных полей скорости, давления и температуры течений в каналах. В кн.: Турбулентные сдвиговые течения I/Под ред. Гиневского А. С. — М.: Машиностроение, 1980. — 431 с.
- Kolovandin В.A., Sosinovich V.A., Kravar S.V. To the theory of the mexanism of mixing in turbulent flows with coherent structures. Lett. Heat Mass. Trans., 1978, vol. 5, P" 253 258.
- Хинце И.О. Турбулентность. M.: Физматгиз, 1963. — 680 с.
- Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1974. -711 с.
- Mankbad.i R., Liu. J.Т.С. A study of the interactions between large-scale coherent structures and fine-grained turbulence in a round jet. Philos. trans, roy. soc. London. A. Math. phy3. sci., 1981, vol. 298, N 1443, p. 541−602.
- Юдович В.И. Об устойчивости стационарных течений вязкой несжимаемой жидкости. Докл. АН СССР, 1965, т. 161, № 5,с. I037−1040.
- Гольдштик М.А., Штерн В. Н. Гидродинамическая устойчивость и турбулентность . Новосибирск: Наука, 1977. — 366 с.
- Freymuth P. On Transition in a separated laminar boundary layer. J. Fluid Mech., 1966, vol. 25, part 4, p. 683−704.
- Grow S.G., Champagne F.H. Orderly structure in jet turbulence. J. Fluid Mech., 1971, vol. 48, part 3, p. 547−591.
- Sato H. Stability and transition of a two-dimensional jet. -J. Fluid Mech., i960, vol. 7, part 1, p. 53−80.
- Mack L.M. A numerical study of the temporal eigenvalue spectrum of the Blasius boundary layer. J. Fluid Mech., 1976, vol. 73, part 3, p. 497−520.
- Grosch C.E., Salwen H. The constinuous spectrum of the Orr-Sommerfeld equation. Part 1. The spectrum and the eigenfunc-tions. J. Fluid Mech., 1978, vol. 37, part 1, p. 33−54.
- Legner H.H., Finson M.L. On the stability of fine-scaled turbulent free shear flows. J. Fluid Mech., 1980, vol. 100, part 2, p. 303−319.
- Lumley J.L. Toward a turbulent constitutive relation. J. Fluid Mech., 1972, vol. 41, part 2, p. 413−434.
- Lumley J.L., Khajeh-IIouri B. Computational modeling of turbulent transport. Adv. in Geophysics, ITew-York-San-Francisco-London, Acad. Press, 1974″ S. 13A, p. I69-I92.
- Tavouleris S., Corrsin S. Experiments in nearly homogenous turbulent shear flow with a uniform mean temperature gradient, Part 1. J. Fluid Mech., 1981, vol. 104, p. 3H-347.
- Мачис Г. П., Шланчяускас А. А. Расчет переноса тепла и импульса когерентными структурами в дальнем турбулентном следе за цилиндром. В сб.: Проблемы турбулентного переноса. Минск, ИТМО им. А. В. Лыкова АН БССР, 1979, с. 82−89.
- Tam S.K.W., Chen К.С. A statistical model of turbulence in two-dimensional mixing layer. J. Fluid Mech., 1979, vol. 92, part 2, p. 303−326.
- Liu J.T.C., Alper A. On the large-scale structure in turbulent free shear flows. I Symp. on turb. shear Plows, Pens., 1974, p.11.1−11.11.
- Мачис Г. П. Расчет переноса импульса когерентными структурами в турбулентном следе за цилиндром. В сб.: Исследование по гидродинамике и теплообмену. Новосибирск, ИТФ СО АН СССР, 1980, с. 51−57.
- Rey С., Schon J.P., Mathieu G. Transfert spectral a grand nombre d’ondes der champs turbulents scaleires passifs a nombre de Prandtl tres different de l’unite. Proc. 6th Can. Congre Appl. Mech. Vancouver, 1977, vol. 2, p. 655−656.
- Uberoi M.S., Freymuth P. Spectra of turbulence in wakes behind circular cylinders. Phys. Fluids, 1969, vol. 12, N 7, p.1359−1363.
- Naudescher E. Flow in the wake of self-propelled bodies and related sources of turbulence. J. Fluid Mech., 1965, vol. 22, part 4, p. 625−656.
- Petersen R.A. Influence of wave dispersion on vortex pairing in a jet. J. Iluid Kecli., 1978, vol.89, part 3, p. 469−495.
- Nayfeh A.H. Stability of three-dimensional boundary layers.- AIAA J., 1980, vol. 18, IT 4, p. 406−416.
- Plaschko P. Helical instabilities of slowly divergent jets. -J. Fluid Mech., 1979, vol. 92, part 2, p. 209−215.
- Plaschko P. Stochastic model theory for coherent turbulent structures in circular jets. Phys. Fluids, 1981, vol. 24, U 2, p. 187−193.
- Чэн И. Волнообразные вихри в турбулентной струе. Ракет, техника и космонавтика, 1977, т. 15, № 7, с. II8-I29.
- Chen Y.Y. Spatial waves in turbulent jets. Phys. Fluids, 1974, vol. 17, N 1, p. 46−53.
- Проудин А.П., Фельдман С. Новая модель пульсаций и смешения в турбулентном следе. Ракет. техника и космонавтика, 1965, № 4, с. 33−42.
- Fabris G. Conditional sampling study of the turbulent wake of a cylinder. Part 1. J. Fluid Mech., 1979, vol. 94, part 4, p. 673−709.
- Hussain A.K.M.F., Reynolds W.C. The mechanics of an organised wave in turbulent shear flow. Part 2. Experimental results, — J. Fluid Mech., 1972, vol. 54, part 2, p. 241−263.
- Lessen M., Cheifetz M.G. Stability of plane Couette flow with respect to finite two-dimensional disturbances. Phys. FluidSj1975, vol. 18, N 8, p. 939−950.
- Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М.: Мир, 1980. — с. 279.
- Гиневский А.С., Власов Е. В., Колесников А. В. Аэроакустические взаимодействия. М.: Машиностроение, 1978. — 177 с.
- Wells M.R., Stock D.E. The effects of crossing trajectories on the dispersion of particles in a turbulent flow. J. Fluid Mech. 1983, vol. 136, p. 31−62.
- Либби П.A., Вильяме Ф. А. Основные аспекты проблемы. В кн.: Турбулентные течения реагирующих газов/Либби П., Вильяме Ф., М.: Мир, 1983. — 325 с.
- Broadwell J.П., Breidenthal R.E. A simple model of mixing and chemical rections in a turbulent shear layer. J. F^luid Mech., 1982, vol. 125, p. 397−410.
- Miksad R.W., Jones F.L., Powers E.J. Measurements of nonlinear interactions during natural transition of a symmetric wake. Phys. Fluids, 1983, vol. 26, N 6, p. 1402−1409.
- LaRue J.C., Libby P.A. Statistical properties of the interface in the turbulent wake of a heated cylinder. Phys. Fluids, 1976, vol. 19, N 12, p. 1864−1875.
- Van der Hegge Zijnen B.G. Measurements of the distribution of heat and matter in a plane turbulent jet of air. Appl. Sci. Res., 1958, vol. 7A, p. 277−293.
- Jenkins P.E., Goldschmidt V.W. A study of the intermittent region of a heated two-dimensional plane jet. Rept HL 74−45, Civil. Eng. Dept., Purdue Univ., 1974, p. 1−220.
- Фабри Г. Турбулентная температура и характеристики теплового потока в следе за цилиндром. В кн.: Турбулентные сдвиговые течения. Ч. I. М.: Машиностроение, 1982. — 432 с.
- F’abris G. Higher-order statistics of turbulent fluctuationsin the plane wake. Fhys. Fluids, 1983, vol. 26, N 6, p. 1437−1445.
- Мачис Г. П. Волновой анализ переноса импульса и тепла турбулентными вихрями. В кн.: Структура турбулентных течений, Минск, ИТМО им. А. В. Лыкова АН БССР, 1982, с. 49−56.
- Crighton D.G., Gaster М. Stability of slowly diverging jet flow. J. Fluid Mech., 1976, vol.77, part. 2, p. 397−413.
- Strange P.J.K., Crighton D.G. Spinning modes on axisymmetric jets. Part 1. J. Fluid Mech., 1983, vol. 134, p. 231-^45.
- SUb ROUTINE L12RKV (flOOTrNM, V>l, P4"XtCO, B, C1tC2iC3> ПОДПРОГРАММА HHIfcPftOCA if УПР АВЛЕНИ* Ь С t- -I численным процессом, ПАРАМЕТРЫ ПОДПРОГРАММЫ, ЬЫбиРАе^ЫЕ ПОЛЬ3У 3ATР '!Е М:
- ЙОСГ-ПАРАмЕ! Р ЙЬ. НОРА КОПИЯ Б УР, 2. 1 4ьfacse,-первый корьнь.fKUfc.-rtTOPOaKOpErtb,
- С NM-ЧИ СЛО JНАЧЬНиЙ, ПОМЕЧАЕМЫХ Ь ОПЕРАТИВНОЙ ПАМЯТИ" ВСПОНОГА ТЕЛЬНЫХ
- С КОЭФФИЦИЕНТОВ ПРЯНОЙ ПРиГОИК^.
- С РЦ-ЗИАЧЕЧИК Ui HOC t t Л b Н и й ОШИЬ К И «3 АА АВ, А ЕНО Е ПОЛЬЗОВАТЕЛЕМ. с массивы *, в, to, ci, ьг. сз / переменные г, i ?1, Т4 i, 1 иг, Ti.3, г ч4, гои г
- Г ИСПОЛЬЗУЮТСЯ „НИР* ПОДПРОГРАММЫ АЛЯ О С У лЬС ГЬЛЕи л Я *У6ИЧЕСК0й СПП
- С Ай (ИНТЕРПОЛЯЦИИ ПОМЕЩЁННЫХ i ОПЕРАТИВНОЙ П4"ЯТи РЕЗУЛЬТАТОВ ПРЯ
- С Рииь. П Аг Я Г И КМ (вЫЬиРАЕ*иа ПОЛЬЗОВАГьЛЕМ.
- C0*PLEX*16 v (1 0)“ y к < 1 q)# с 0
* с 2 p с 3 (n м, в) RE* l*b PR, YH. X (hM), g (NMJ, T, T31, TM, TUZ, T*3 RE, TOUT1. GICAL ROOF ^(S = 0 - Г ОПРЕДЕЛЕНИЕ ВХОДНЫА ПАРАМЕТРОВ, ivm, а x = о 1 н t = Р R
- GO ТО (<�», 1 00>, Н&- С РЕЖИМ ПРЯМОЙ НРоГОнКИ. О= YMroutсо. DO GO ТО 1с решим 0е>ра1л0й пригонки.4 Ь U = N S 1 Р1. F ь = * 0 +1 к к л F е +1
- С 0П (- ЬДЕЛЕНИЕ НЛОАНЫА liAPAMtlPOb ИНТЕРПОЛИРОВАНИЯ.
- Х-КАТРИЦА УЗЛОН АЬ ЦИ С Ы. И, СО, С1, С d, СЗ-КО-ЭФФИа»" еН I К| КУБИЧЕСКОГО1. С СПЛАЙНА.1. NM’i =fvM-1 1sX (NO)
- J) / T 4 2−02
/Т439 uo 10 i"nfMM TxX Cl-1)/BCI-1) B (I)"B (i)-T*X (I-1)oO 10 J = 1, й - DO 12 J = *(J, I 4 1 s X (I J
- CI II tJ)-(CO (l+1 (J)-L (Jtl, J > > / I41-T41*+?.0O*C2(i, J >) C3t I, J≥(C2tl+1, J — l 2 (1fJ))/T*1
- C2 (X, J >«3. оо*сг С I, «I)
- C2ll*MrJ)s3.D (j*C?(t m, J)
- C3iN*, J > *C3(KM1, J) 40 = 41.si.00 t 001 = YM
- С ОПРЕДЕЛЕНИЕ НАЧАЛЬНЫХ ЗНАЧЕНИЙ ИНТ ЕГ РИ РУ ЕМЫХ, АБСОЛЮТНОЙ
- OdikhXM И УКАЗАТЕЛИ Р t * И М ^ ИН I ЬГРАТиМД L 1 2 I С V .
- CALL L126CV t ROOT,"5,ГИ, Y) 131=0.001.LAG = 1
- F
GO 10 1: И H l ЕГ РИрОбА Н/t f С11С li, но ЗАКОНЧЕНО, ЕСЛИ ЭТО ЬЫПЖ ОБРАТНАЯ ПРОГОНКА1. С It ОН ELl ЗАДАЧИ. return - PRINT 31, RE, I 31, XFLAb GO Ю 1540 PRINT 41 *M/'X=-NMAX + 11. f (NhAX. L Г. 2) I 0 151. RE I URN 50 131=1.U-9
- PP i NT 31, R E * T31 j I F LAG bo 10 15 fc 0 IF (RE.LE.1,D0> uu TO o2 P R) N I 61 (RE, 1 «T 4 1 RE I URNь2 PRINT 31, RE, i 31,1 FLAG «E=10. OUO*RE I F L A G я 2 60 Ю 10 70 PRINT 71, T, I 41, RE, I 3 1 I FLAfa*2 (it) 10 15 ttO PRINT 81 RE1URN
- С HPv ПРЯМОЙ UPOrOHKt ПЕРЕПОЛНЕНА ОПЕРАТИВНАЯ ПАМЯТЬ, ОТВЕДЕННАЯ ДЛЧ
- С РАЗМЕЩЕНИЯ вСПОМОГА) ЕЛЬНЫХ ФУИKUHй. НОЛЬ3О8АТЕЛЕf ДОЛЖЕН БЫТЬ УВЕ1. С ЛИЧЕН ПА Р AM 11И NM.
- PRIM 91, I, I 41, «М 100 RETURN
- F О RM АТ (10Х"1?Н lOCtRANCtS RtSEf, 2t)12.3,Ib)
- FUKMM I (1 OX, dc HO VfiH f LO T MEMORY VOlliHi10X, 5hA{ Ys, 09,2,2X, 3HUr=, D9t2t1i'HMEM.4A1H.t>I4ENS. =, l4) end
- SUBROUTINE UJbCviKUOl, i S If S Г, Art, С У >
- ИНТЕРПОЛЯЦИЯ ЗНАЧЕНИЯ К и ЭФФ vЦИ EH 1А СОПРОЬ- вЛЕНИЯ ЦИЛИНДРА ПРИ
- CBtj-=E*CB (j) с bbU) eA1*SBlD 00 3 1=10,12
- SB к 1)=E* SB (I> DO 4 1=1,4 U С Y (l)s (O.ODOjO.OOO- СУ (9)в (О.ООО, О. ODOJ i t = с4 + с3
- CO=COSQRT (CE*P+C2) С Y <7 ) = А •> * < 1. О DO+C I i CYt8) s-CQ*CY<7> 1 = 11. J F tROOT) 1 = 21. CC=CDSQRT (CE)1. F=lCC*A*t-CX*A3)/(CC~AL)1. С Y I10) = F*(C1+Ct:/C1>
- CF.= (Cl*A3-AL"A4)/(LL-ALJ1. CY (5) = fc + Ch: t Y ib) s-(*L*tttC*Ct)
- С ПЕЧАТЬ MAPAMfcTPOB в ЗАСОЛОВ K? IАЬЛИОЫ РбЗУЛ^ГАЫо.
- PRIM 26, l, R, bH, P,*F, PT, C, XO, A1, XEO, OTCO. A2,A5,A3,A<., L (., CC GO 10 16 С ОБРАТНАЯ ПРОГОНКА.
- С ВЫЧИСЛЕНИЕ КСРНЬй /РАВНскИй . 1^> .
- If С ROOT) GO ТО 7 00 6 1=5,7,2 J"t 1−3J/2 С Q t С У С I)
- Р1 W) =)» ССУ U-4)) Р2(J > = OKt CQ) с P3tJ) sDI (СО
- F = ^F1*F4-F2*F3)/ <1.00+F*Cl)/ (j О I О 3 6 7 C?*(O,DO,O.»0) 3Ь с i t 1> = c?
- DO 15 1=2,Ь 15 СУ (I)s (V, DO, О,00) f ПЕЧАТЬ ОБОЗНАЧЕНИЙ РЕЗУЛЬТАТОВ НАСЧ с I A.1. HRlNl 25 1 в ftt t URN
- FORMAT (3 в X r ' G Л M V «' 0.3)
- Sub KOU TINE l12tXI (.NSTP, t4Y, OX. X, Y, YO, C1, C2,C3>
- С ПОДПРОГРАММА ьызии* НА ПЕЧАТЬ РЕЗУЛЬТАТОВ OhPAThOU 1 °F О ГО H К И ПО
- С ЗАДАИИО? ПОЛЬЗОЙ, А 1ЬЛЬ Ч ФОРМЕ .nPtA/CMO I Р Е Н A ВОЗМО#НОС1Ь ВЫЧИСЛЕНИЯ
- С АРУ Г tf X ПРОИЗБОДпНХ ВЕЛИЧИН, КАК АМПЛ и Т У Д НЫЕ, Ф, А 30 ЬЫ? ХАРАКТЕРИС1ИКИ
- С РЕШЕНИЙ, ТУРБУЛЕНТНОЕ ЧИ 1, Л0 ПРАНДТЛЯ И ДР. мАС&ИВЫ СХ, F 2 л ПЕРЕМЕННЫ
- С ENERGY, SHEAR ИСПОЛвЗУЮТСЯ ВНУТРИ ПОДПРОГРАММЫ.
- Г N Y-ЧИСЛО ЗНАЧЕН/Й, ПОМЕЩАЕМЫХ 6 ОПЕРАТИВНО* П, А «1Я I И, Р Ё 3 УЛЬ Т, А Г ОВ НРЯ1. С МОР ПРОГОНКИ.
- С * ¦ AL СИЬ ВСПОЧОГ, А 7 tЛЬНЫХ ФУНКЦИЙ F1,., F10.1.f i 1 Р, L>X, X, Y 0, l I СЗ-ПАРАМЕТРЫ, а КОЭФФИЦИЕНТЫ КУЬИЧЕСКОГО СПЛАЛНА.
- С X-ltKyuiht ЗНАЧЕНИЕ ПОПЕРЕЧНО* КООРДИНАТЫ.
- С EV1. ЕУ2-ЭНЕРГИЯ ПРиДОЛЬнОй, ПОПЕРЕЧНОЙ КОМПОНьНТ КРУПНОМАСшIАЬ
- С ПОР (¦УЛьСАЦ** CikOPiiLT/.
- С Еу-ЭНЕРГИЯ КПУПН9ММ iT АЬНЫХ ПУЛЬСАЦИЙ С КОРОС I И,
- С SV-PEBHUflbACOBOfc НАПРЯЖЕНИЕ КРУИНОМаСЩТАЬНЫХ ПУЛЬЦАЦИй СКОРОС1 И.
- С fcd-^НЕРГИЯ КРУПНОМ At 1НА6НЫХ ПУЛЬСАЦИЙ ТЕМПЕРАТУРЫ.1.&н1, ЬН2-ПР0Д0/1ЬнЫй ПОПЕРЕЧНЫЙ По t 0<и Т 6−1Ла,
- С PRI IУРЬУЛЕнIHOt ЧИСЛО ПРАНДТЛЯ .
- С f-→bA30Bbt УГЛЫ ВОЛНОВЫХ АМПЛИТУД. ПЕЧАТАНИЕ ЭМХ ВЕЛИЧИН УПРАВЛЯкТ
- С 1Я ОПЕРА ТО PUf «3 h’xiM 11».
- COF Р LEX *16 V (1 0 > t Ud, С X () f С F1, С > 3 «С Р 5 (CF?, Y0(NY, 8), CI < N Y «8) ,
- C2tNY, 8), C3(NY, 8) R Е * L * 8 X, OX, ENEReY, SHEAR, EVl. EV<>, EVfSV, EH, SH1, SH2,PRT, F (2>, P11(P12 toe MOW / CBtx/ PI 1, P i ←
- EQUIVALENCE (CMrCA (1}),(CF3,CX<2>>t<3>>f<4)> tNERGY (U2)=U2*DC0NjG (U2>
- S>HfcARCU2,CM)=U2*)C0KJG (CF1>+LF | *PC0NJG (U2)oo i i = i, f, a1. J a (X + 1 >/2
- CX
sYOCNS r P, J) + i>A*(C1 INSTP, I)+DX* tCg +DX*C^(NSTP, I>> J U^ r (1 >1. CF1*U2*CF1+CF5 - EV1=tNER6Y (CF1) t tfi: = fcNE"GY (02>5 v = ShE An
- E H = E NERGY (С F3 > 1 = 1
- С F is LiCON J в (CM) EV** PR) = t F 1 + С F 3
- F, 6T, 1, 0−1 0) C.0 10 6 F ®0.006 0 ГО Й
- IF tiJABS (PRI) .(jT.1, J"10) GO lo f > < j)=90.001. GO 10 8i F (1)=5?.29S800"DATAN2(EV, PRT) о 1 = i + 1
- F
- MX. fiT'. 1 *0U~4> nt TO 2 >'Rial. ООО ti 0 TO 31. С bolfcUA HA IE 4 A I b.
- HKIM 11, X,?V1,EV2,EV, SV, EHtSH1, SH2,PRT, F RE i URN
- CA-CA+C3*<1, 000-UMnl'WI)1. Св=р* с A + C2dr2=tp*x*dr1**pt0 R 1 * P U * X *) R 1
- CALL L12T"V<.t:.4li>c., AiU43, I>K
- UCi 1 1 = 13,1 A CF U) = CF (I > /CF15 J = 3 * 1 3 7 С м j > =C F
+CA00 2 1=1,3 IF (I > = CF (X)/CF4 CE = CF (Г) J= I1. CF (J≥CF (J)+CF8*CE J-j - C0J> = CA*CF (- + C8"CFJ + CE*cfAc0u3) = ca*chfc + c"*ch7 + ce*cfa
- C0
- DO A 1 = 2, A, 2 J =ч * I-7
- ClH 1)=CtHI> + CKJ> CO t^> = Cb (4)+CFCS)1. C>0 10 22
- ОБРАТНАЯ ПРОГОНКА t A = С Y С 1 >
- СПлАРН-fH ГЕ^ПОЛЯЦИ* 1АбУ/|ИР0ВАННЫХ ВСПОМОГАТЕЛЬНЫХ ФУНКЦИЯ. if (inpr.eq.o) 1)0 10 2a 60 ГО С t>, 7,6, V, 10>, I NPR oy-oy +h8 GO TO 117 uy=dt+he GO I U 11 Ь О Y = L>Y + 5. 700*F8/1. 3 U Сfau го 11 4 DY = DY+8.00-*tl8/13.DU
- GO ГО 11 1С) l’Y = OY + He 1 1 URIsXl (M S T P >
- С A I L ЦгГКУ1. Г-{иЕ., Х, ОК2"0Ь!3,П4(Сг) DO 12 1=1,13,3 1? Г F l I) = CA*CF< I) + L H * l F (1 + 1) ±C?*CF (I + 2) 00 13 1=10,13,3 J=2*(I+14)/3
- CF4) = CF (I>+CC*CF l>KI=DR2*(?htPGYiCA-+bNEr<6Y(C*>>/4.00 dk3s-or3*shtar (ti, i. h>/4.00
- D R 5 = S H 1 (CF1,CB> + SHKCF4, CA>+SHEAfi (CF4,CF>*$Ml С С F 7, С 6 >1. OR5=-DR5/4.O00Ri*)H2*E4ERGY (CC>/4,0O
- PR4s-UR4*SHEAR (CA (Lt>/4.D00r6"SH1(CF10,CC)*SllEAR (CF13|CG)dr6s-dr6/4,do1. ОС 14 1=2,4 J-С * I-3
- COU> = OCMPIX
- SUB ROUT 1ЧЕ L121fiVlL, X, U0,UI «T1 / С F)
- Г ПОДПРОГРАММА вычисляет КОЭФФИЦИ tH Г hi 6 (К) f H С К.) ПО ЫРАВНЕНИЯМ (2.12)
- С Г (2.В) .ПАРАНЕ 1 Р ГЧЙ ОПРЕДЕЛЯЕТ /НЕТ ПОЛЯ ПУЛЬСАЦИЙ '1ЕЛКОГО МАСщТА
- С ьк / 3 А Д A f» ОПЕРАТОРОМ Си -1UN / С 1 5 / Г К В В ВЫЗЫВАЙТЕ*, 1РUГРAMМе. IRВ =
- С. IKUE.-НЛИЯМИЬ М1=Л*ОМАСЛА6НЬ|Х 11У Л С Л Ц И Й УЧИ Г Ы В, А с ГСЯ, Г R 4 = .FALSE.1. С -чI УЧ/ТЫВАЕ1 С Я .
- С ПОДПРОГРАММА L12T4V, Н С>11НОСТИ, ЯВЛЯ?1СЯ ЧАСТОЙ ПОДПРОГРАММЫ
- С L 1 2 О У V, В КО ГОРОЙ ВЫЧИСЛЯЮТСЯ ВСЕ ОПРЕДЕЛЯЕМЫЕ ПОЛЕМ мелких ПУЛЬС А
- СOF MON /C0b/TH, C0/bO4S/KfPr/C3 245/Cl/C15/TRB/e3oS/U4/C13S/D1fD2i, P t, С H <6 ) , S A v В ( 1 й> J)
- Acl, 1)=SAVBCll*X*SAVb00 4 1=9,10 J = 3» 1−26
- CB (J) = SAVBU)+X*SAvB (I + 2> A (4,1)EX*U1*SAVB (8J+U2*$AVB<14)-*SAV8<4) a (2,1)"X*SAVB1. A<1 , 3) =U1-A(1 ,3)cei3)»Ti*c8(3>
- CBtA)sX*T1*SAVB<12)+T2*CB<3>+SAVat10) С в l 5) = T 1 * С b (5 > Q = L"2*Q**r>1 011SU1*U1 6 = t .365o0*U1
- F (OABS (DKCO)). 6 Т. 1 .00−1 5) 60 ТО 22bi = св
- С0 = «1*<1 .1)0,О.00) 22 С О = W С О00 5 iя 1,3 jo 5 j = 1, 3к = з * с i -1) + jcf (obcq*cc*ftu, j"+c*ft<2.j>+ccxt3)*4<3fj> + cw. o*au"j>> c11sc1*(cq-u0)1F (OabS{DICC11)).Gi.1.0j-15) IjO TO 2 4 Bl=C1 1
- С 1 1 = t>1 * (1 .0 00,0.0 00) с 4 f: s = т 1 / с 11
- C4S0/C11 C12=U1/C11 00 6 1=1,66 cfu+9) = c5*cf (i + 3)"c4*c, ui>00 7 1=10,127 (. Fvl> = CKl>+C12*CF4+3)
- S 1 = 13, 1 a ь cf u +3 > = c4*sa vii w)00 9 Is 16, 1 7 у LFlI> = CF (i)+C12*CFU+2> lFl. L> GO Г0 31 00 11) 1*2,3
- Ca<3M-5) = CQ*tS-,*C(I,1>*S6*CCl,*-)+S?"C (I, 3)+S-J*C (I, 4>> Й1=-0.36 500*01
- C1<�"s (T2 + 2 .0t"0*01*X"T1)/C1 1 111 = U1 M.1 +1. 0 798)0"U2
- F ч OABStDi CC2>). GT. 1. 00−1 5 > GO TO 25 I 2 С 33= I 2* (1. ООО, 0.000)25 !f (pabs (oi (c3>>.gt.1.00−15) go to 2® 1 2 = c22=) 2* (1. 000, 0. 000)
- CI 2*(CI 1-B1 >/CZ Cli=tCl1+rt1)/C3 c2i=b*c2*(2.000*u2+U1*c12) C2^ = -0 2"c2*C3 + c22"(C12 + l.13) C2b=C21
- C2^ = C2^"i-i2 + C2"(C3-(C11ri-A*C13)-4.0l)0*B*ai)
- C31=0.365 00*t2*C2*t02"2.000*01*C12)1. СЗ^ = С32*(CI 3 + 2.*C12)1. C33=C31сз<* = г. уоо*С21
- С2*-2.0Г>0"01 + X-C1 2-C2*CC2 / СО DO 12 1=2,3 00 11 J = 1, 3 * = 3* (I -1)+J
- C-3=C0*(C (It 1) * A 11, j)+C (i, 2)*^C2,J)+C (I, 3)*A (3,J)+C (I,»»)"A (4,J))11 C (i-1,J)"C3+C2*CF<0
- C (I-1rD = CtI-b1)+CH<3*X-5) C3
- C3<*= 1 2M 1. 0 DO , 0. ООО>27 «0 13 1=1,31? t (:, I) = C1*t"CF (I + 6) * (34*-, F (I + 12)+C5"C<2, I> C. S)xC31+C4*SAVB (121 00 1*. 1=1,2 C (1,4)4C1*CF
- С (д, 2) = i { I, 2) С 1 * С t- < J) 1 ** C (i, 3) = C (I, 3)+CF (J + 1) j =600 15 1 = 1 Г 3 JeJ+3*(I-1)
- A (1, 1)=CF (J) C32"C32-C1*CF<13) C33»C33 + CF(lO00 16 1=1,3 LnI):C1«CMDH(1,l) CF(I + 4) = C1*CKI»-3>>C (2,1) CF U"-8)*C1*CF (i"-9)*t<3* t ) J = t* I16 С F (. J ) =A (. I , 1 )1. ОС 18 1 = 13,1418 (,FU)=C1*ef-ti+3>±Cvi-12,<*>
- CfHO^fUt (14)ftF<18)) 00 19 IB 1,619 CFlI)=R*CF(l>ои го is», i3 с f и > = p e * с к i >
- С f (4)s1 .0l>0-CK"t) С F <) 5) = 1. OOd-Pt • Cf 11 9> ои гз i=4,15,111. Cti = CF 11)
- M «ARS (iH
- F i I) St 8 1 • (1. и 0, О. J s)) con г 1nuc RETURN ENDой-is) go 10 23
- SUbftOU f I ll21CV<�«SrNQr'4V(V*.Tf)OUT,ftEL=RR,*B3Ellft.XFL4G,NFEiNSTl>< Or, H, YP, SRE, SAE, KOP, INI I, JFL
- ПОДПРОГРАММА РЕШЕНИЯ ЗАДАЧИ KOUJH ДЛЯ ОБЫ <40 В E НИЬ> X Д ИФ E P с HUH, А ЛЬ h Ы X
- F (IFLAG.EQ.2) GO IF LAGsJ F LAG J F LAG=IF LAG * F A G я 0
- DO ы < = i, a to (f*stp, k) sy (k) GO 10 63 0 Y = 0.00
- CALL L12EXI (NSTP, NT, DY, T, YfC0iC1,C2,C3> CALL L120YV (NS, l, Y"*P"0,'STP, NY"QY"EE, X, C0,C1u 0 to 51. ИНФОРМАЦИИ1. GO TO 20
- А ли. (
, (ABi>E RR Anu. С RElE RR «О TO 30о О TO 4 > SO T (: 4 и A’I у. (ABSfcRRti 1. u. oou) > oj 10 45значён/й10 5 0i pahhl погреjhoci ей - СЕРЬеЗНО ВЫХОДНЫХ, b 0 Tu 5 51. ВЛИЯНИЕ ! u4rk.
- ПА L 1 2 I С V C/i/uXOM ЬОЛЬШОГО ЗА1 О 6и 10 6 5с 2 г С 3) м-1= 1
- X F t. Г. NE. TDU I) lit) i 0 65 GO Ю 3')0 INI = 1 H sOABS (DT) I’lJL^O. 0 DO00 70 ≤ I i Nil
- T0Le"ELE4K"C0AtJS
- J F (N F E. Lfc. Mii) GO i 0 22 J IF L A b = A L A G = A RE I URN I H T = N S T P
- ВНУТРЕННИЙ ОДНО^АГОВЫй ИНТЕГРАТОР. Sa D У1. А=н/A.ООО
- НО 221 К = 1, rjg F5(К)*У (К)+А*УР tK) CALL L12DVV (NS, ri-A, FS, Fl"1,l4T, tiY (S, A, X, C0, C1. C2,C3) E E = A/2•00 A=i, ooo"H/32.0Qo 00 222, U F5tK) sY (K) + A*(YP<((>+3,000*F1
- CALL L120yv<4S, r*3.t)0*Et, FS, F2,2fINT#rJY"S|EEt-<, C0, C1"C2"C3)t = i NT
- A = «i / 21 77.0 00 00 223 K = 1 ,
- F5^K) = YCK> + A*C1932.0 00*YP<7296.0i)0*F2 call L1 2 О У V (4S, Г +1 .2)0*11/1 .J00"t*5#F3"3» I, *Y, YPK, E (-, X, C0,C1, С 2. С 3) As) / A 1 OA. 01)01. DO 22t K*1t4Q
- F54O = Y (O+A*(ttt3Al.000*VP (K)-845.()i)0*F3C<)) + t29AA0.000*F2-<)-- i 2 tv 3 2. О 0 0 * F 1 (K))) CA: L L12OVV A =! / 20 520.0 0 000 225 Kat, HQ
- ЦtiOsY
*A*l (-ftGe0,0D0*YP (<)*<9295.0D)*F3 -5643,QD0"F4(K)))+ ¦ tA←.040"ODO*F1(K)"2e352.000*F2,5,IftlT, r, S, EE,*, C0, C1. C2,C3> Ash/ 7&1 3050 .000 00 230 K=1» - F1lK) = Y (K)+AM (90 2680tOL>0*YP (K)*(3 855 735,)DO*f3(K>-1 371 249.000*3953O6A, 0OO*F2(K)+277O20,000*F5<<))) NFtsNFE+5t EC t I =0 .01.1000 250 K» 1, Nu
- У СПЕШНЫЙ tj, А Г. 1 :) +hgo i О < 263,264), ns iws i p = ns 1 p-1
- CNSTP.GT.O) kJ 10 261 1 F I. A G = 9 RE I URN x ($t p)=-h00 2t>2 n = 1, tt Co (NST P, X) = F1 <. К) GO 10 264 *tS 1 P= I D Г = If P К
- H = t, sigmc0^AX1(ee*0MBS (l-l>, HMIN)/h)1 F (00PUT > GO 10 iOO
- ПРОВЕРКА ПРОДОЛЖЕНИЙ ИНТЕГРИРОВАНИЙ.1. GO 10 10 0t = I OUI
- ИНТЕГРИРОВАНИЕ УСПьШНО ЗАВЕРШЕНО. I F L A G = 2 «е1 urn emufiPkMfcP ВЫЗОВА ПОДПРОГРАММЫ L12HKV. ФИЗИЧЕСКИ.- ПАРАМЕТРАМ tf ОПЕРА 10-Pfc СОМ Мим ПРИСВОЕНО. ЧИСЛОВЫЕ ЗНАЧЕНИЙ.
- К t ^ L * б AL, DtCO, XEO, XD, SH, CI"P, R, PT, PR"W1<300),W3(300) CQMPLEX*1& -/2t300,o>, ui<.c300,8), * 5 t 3 0 0 , в>. v6(330,») COMMON / С 03/A L, О T С 0, ХЕ 0, X D, S H, CI/C04/P/C045/R, PI н I. 5001. SH=0.300 R = 1.03 P=0.700xo®100.00
- XE0=10.00 DTC0*1.0−4 С I = -0. 500 A L = 0. 6 UO
- CALL L1 2hkv t. FALSe f, 300, 6 .000, 1. 00−2, W1, «i2, *3, *4, W5, W6) sfopfc h U1. BLOCK OAlA
- Начальник отделения № 7 Начальник отдела № 76
- Начальник сектора № 761 Экономист отделения № 7У1. Е. П. Марочкин В.Л. Миков
- А.В. Мезенцев Т. Г. Каляскина
- Знв.лабораторией турбулентности, к.т.н.
- Завлабораторией математической теории переноса, д.т.н.1. А. Коловандин1. М. Колесников
- С.н.с.лаборатории турбулентности, к.ф.-м.н.1. В.А.Сосинович