Н.А. Слёзкин - Динамика вязкой несжимаемой жидкости (1124064), страница 49
Текст из файла (страница 49)
Эти два обстоятельства и побуждали многих исследователей искать другие приближенные методы решения уравнений для пограничного слоя. Большая группа этих методов, получивших наибольшее применение к решению отдельных задач, основывается на специальном выборе независимых безразмерных переменных, позволяющем аифференциальнью уравнения с частными проивводными (! .13) сводить либо к одному нелинейному обыкновенному дифференциальному уравнению с числовыми коэффициентами, либо к некоторой нос тедовательности обыкновенных дифференциальных уравнений также с числовыми коэффициентами, В этих методах численно решается обыкновенное уравнение или группа уравнений и составляются соответственные таблицы.
Эти таблицы затем могут быть использованы для целой группы соответственных задач (а не Одной какой-либо аадачи). Чтобы показать конкретно сущность этих методов, рассмотрим подробно метод сведения уравнений (1.13) к одному обыкновенному уранению для случая степенного закона изменения скорости внешнего потока, развитый в работах Фокнера и Скэна, з) Хартри з), Л. Г. Лойцянского') и др. Пусть на внешней границе ламинарного пограничного слоя скорость частиц во внешнем потоке распределяется в продольном напра- ') Кочин «!.
Е. я Лойцянский Л. Г„Об одном приближенном Методе расчета ламннариого пограничного слоя, Доклады АН СССР, т. ХХХН1, № 9, !942. В Р з ! и и е г Н. М. а и д 3 и и и 3. '1Н., АКС К АМ № 131 4, 1930. ') Н а г ! г е е ««. й., Ргое. о1 Гае Сзшьгщйе РЫ1, 3ос, 33, !937.
ч) Лей цянс к ий Л. Г., Механика жидкости н газа, Гостехнздат, 1930. 5) пвивлижйнный метод гашения тглвнвний погвлничного слоя 273 еленин по закону и(х)=с», (5. 1) уравнения пограничного слоя (1,13) при использовании (3.12) и (5,1) представятся в виде ди ди „,„г дэи и — + о — = сетха'" ' .+ э —,, дх ду дул' — -'- — = О. дх ' ду (5.3) От размерных координач х и У перейден к безразмерным, полагая х = (х„ г (5.4) У= „— Ут Если в выражении числа Рейнольдса ьт К=в (5. 5) заменить скорость через коэффициент (5.2), то получим: (5.6) Полагая и = (гин и используя (5.6), уравнения (5.3) можно представить в виде ииэ диг „„,, дэиг и, +о,д-.— глх л + —,, х, д"'+ д" =-О. Решения уравнений (5.8), вообще говоря, должны быть функциями бевразмернык координат х ' и у,, Поэтому при переходе к размер- ным координатам размерные решения уравнений (5.3) будут зависеть от масштаба длины Е Так как коэффициенты уравнений (5.3) не (5.8) где с — размерный коэффициент, свяаанный с размерностью длины .и скорости соотношением с= —.
(5. 2) Ги л74 (гл. сш тьогия погяанияиого слоя созсрагат я инион внлс нзсшглбз ллинн ), тз ножно потрсбояать, чгоби бсзрззмсрныс скорости а, и о, лапнсегит не от отъ .и.нюх без. разнерннк коорлинаг .т, и у,. а от такой ич кочбннзцнн. ирн которой исключалась Гнз ззьнсичость рз цгернюз скоростей ог насш таба ллигц» !. Если обрати»»се к псряпчу рзненству (5 4~ и кн и орикс рзвснс гну(5 6). то чо«но ззнстнть ьго ~ анис~и еп»й личГ»ишки~я, ис солертизц(сй нлсшгзбз ллинц нрн нерстоле к ззлчсрнцн лоорлннттзч.
б)лит; « ° 3 ! На»тои оснояаиии ииотин новос нсззн!шнчое безрлзнсрное переменноее т„голл г ли / г.е" оя / ("(л) т,-ь тр - -- — у $г ( 1.0) ! ле з — и»ко!орос постоянное число, 1(з (5.9) бу гьч ниеть д / С'(.ь) дУ 'Р (5.10) дг 1, 1 1 (: глс шгрнс н прзюш иль~и озцюю г лифф»ренцкронзьис но т, 1!з осионзнии успения црилицз,ня и того. ч:о низ,няя грзиипл с,ню ирелстзвляст сосюз линию го»ш 6) теч ичг ть слету»нцнс с(ломия лли неизвестной функции ф(г): при т, = 0 Ф(0) .. О. Ф'(О) — И (5 12) 1'.спользуя спи»кои~с,цгя (5.10) и (5 11).
Г»улсь нче!ь. дц . У (и Ги и) нц С(',,:, <»уз ~ 5. 1 3 г !!з урз»испи» нссьичаснсстгц )ьловня обрзщсшш и нуль скор югн с иа с:енле у -. 0 н пер»иго рззсцстпз ().!О) погбякч » гз / .- ! ди ч геп»рь причем. чго п(толею»игя сосгззляюшзи»слгорл скорости и иограничиои слог чоац г бог ь прс.юглнлс»а и ьи ш и (.т, у) — (,' нт) г! '(т,ц (5. 11) й 5! пяивлижйнный метод гашения твлвняиий погганичного слоя 275 Подставляя значение — из (5,13) и испольауя условия (5.12), будем ди дх иметгп -,—,г 1 и()и . =- — УЪ (7 ~"Ф+ — — Рбф' — Ф)1 2 (Рт (5.14) На основании (5.13) и (5.14) получим: д + д ( )+2 ~У ди ди,,з и ! ()з(ги Следовательно, первое уравнение (5.3) представится в виде Наконец, используя предположение (5.1) и поэагая 2ж т лг+!' (5.15) Ф'(оо) = 1.
(5.1 ?) При каждом отдельном значении постоянного р уравнение (5.!6) можно интегрировать численным методом. В цитированной выше работе Хартри приведена таблица 2 значений функции Ф' при различных аначениях параметра 5 и таблица 3 вспомогательных функций, через которые вычисляются толщина вытеснения и", толщина потери импульса 6*" и напряжение вязкости на стенке. Мы приводим некоторые выдержки из этих таблиц (см.
стр. 276 †2). На основании второго равенства (5.13) получим следуюгпее выражение для напряжения вязкости на стенке: — (д ) — 1(71/ ~~' Фи(6) (5.1 8) Толщина слоя вытеснения 3' была определена выше формулой (2.20). Подставляя в эту формулу значение у из (5.9), получим: о*= ~: :~1 — — ") ду = ~ (1 — Ф'(т))) г(т! ° 1/ — и, = АД) 1/ —, ° (5.!9) о о получим следующее обыкновенное дифференциальное уравнение для неизвестной функции Ф: '„';, +Ф вЂ” ",",., =~~( — ";,)' — ~ (5.16) Для решения уравнения (5,16) к условиям (5.12) необходимо присоединить условие на внешней границе слоя. Если считать слой асимптотическим, то дополнительное условие будет представляться в виде д 5) птивдизкйнный мятод вешания тглвняннй погтхничного слоя 277 условная толщина потери импульса, являющаяся мерой изменения количества движения эа счвт образования пограничного слоя, определяется равенством ~ и'( (7 ) "У а Подставляя значение у из (5.9), получим: 5"= ~ Ф'(О)(! — Ф'(т!))йт, ° ф~ й.', =-В(р) 1/ О!7.
(5,20) Значения функций, входящих з равенства (5.18), (5.19) и (5.20), берут из таблицы 3. Таблица 3 В (5) ! Фж (О) АФ) 1,521 1,Г>87 0,498 Заметим, что случай т = 0 отвечает прямолинейно-параллельному внешнему потоку с постоянной скоростью с, обтекающему продольно- прямолинейную пластинку. )(ля положительных значений показателя в (5,1) мы будем получать так называемые ускоренные потоки, которые имеют место в конфузорных (сходящихся) каналах, а для отрицательных значений т будем иметь замедленные потоки в диффуаорных каналах. Наконец, случай т = 1 мы получаем для пограничного слоя в передней критнческои точне при обтекании внешним потоком выпуклого контура. В этом последнем случае У' = сопз1, и поэтому из (5.19) и д5.20) будет следовать, что обе то.шины не — 0,1988 — О,!9 — 0,18 — 0,16 — 0,14 — 0,10 0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,80 1,00 1',20 2,00 2,359 2,007 1,871 1,708 1,597 0,984 0,911 0,853 0,801 0,764 0,699 0,648 0,607 0,544 0,585 0,677 О,ьбб 0,552 0,739 0,515 0,470 0,435 0,408 0,386 0,367 0,350 0,336 0,312 0,292 0,276 0,250 0,231 0,0000 0,086 0,1285 0,1905 0,2395 0,3191 0,4696 0,5870 0,6869 0,7748 0,8542 0,9277 0,996 1,! 20 1,2326 1,336 ТЕОРИЯ ПОГРАНИЧНОГО СЛОЯ ИГП 1Г.Ч.
чп! будут зависеть от координаты х, т. е. толщина пограничного слои вблизи критической точки будет постоянной. Из этих же формул при учете (5.1) будет следовать, что при т(! толщина пограничного слоя будет расти по течению, как и в случае пограничного слоя на пластинке. Наибольший рост толгцины пограничного слоя по течению будет иметь место в замедленных потоках при лг ( О, тогда как при ускоренном течении при ю,.ь 1 толщина пограничного слоя будет даже убывать по течению. На основании таблицы 1 получается, что при р = — 0,1988, т. е. при Гл = — 0,0904, при убывании скорости внешнего потока по закону и — сх - 0.0000 величина Фь(0) обращается в нуль, и поэтому согласно (5.!8) сила трения будет обращаться в нуль на всей стенке соответствующего канала.
Этот случай можно рассматривать как предельный случай того безотрывного движения в пограничном слое, который может быть изучен этим методом, так как при р ( — 0,1988 пограничный слой либо вообще не может существовать, либо развитый выше метод становится неприменимым. Таким обрааом, задавая различные значения для показателя т, можно получить различные по своему характеру течения в пограничном слое и эти течения буду~ сходны с теми течениями, которые имеют место в отдельных частях действительного пограничного слоя, например на крыле: вблизи критической точки (лг = 1), вблизи точки наименьигего давления (лг = 0) и вблизи точки отрыва (Гл = — 0,0904).