Файл: Щербинин Э.В. Струйные течения вязкой жидкости в магнитном поле.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 180

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Т а к им

образом, при сравнительно общем предположении о

х а р а к т е р е

приближенного з а д а н и я профиля скорости

основным

результатом решения является полный р а з м ы в струи

на конеч­

ном расстоянии от источника.

 

1.2.СТРУЯ У ПЛОСКОЙ ПОВЕРХНОСТИ

Пр и расчете плоской струи, распространяющейся вдоль плос­ кой поверхности, коэффициенты р я д а (4.3) найдем, подчинив

функцию

-^j = / ( г | )

следующим условиям:

 

 

 

 

и = 0

при у=0;

 

и = 0

при у=8;

ди

 

при у = 8;

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.10')

д2и

= 0

при

 

у = б ;

 

 

д2и

= 0

при у =

0.

 

 

 

ду2

 

 

ду2

 

 

 

 

 

r

J

 

 

 

 

 

 

 

 

 

 

Последнее

 

условие

получено

непосредственно

из

уравнения

д в и ж е н и я

(3.1) с использованием

условия

непроницаемости

по­

верхности

V = 0 . К а к видно из (4.10'), электромагнитный член

не

входит

в граничные

условия вплоть до вторых производных ско­

рости по поперечной

координате.

 

 

 

 

 

 

Влияние электромагнитных сил на профиль скорости учиты­

вается

при использовании

производных более высокого порядка .

Е р м о л а е в а и Соковишин

 

[10] ввели та к называемое

условие сов­

местности:

 

 

 

 

 

 

 

 

 

 

 

 

 

д*и

 

аВ2 ди

 

 

п

 

 

 

 

 

4.10"

- г — =

 

 

при 0 = 0,

 

 

 

 

ду3

 

pv

ду

 

 

 

 

 

 

 

 

 

 

 

получающееся

 

путем дифференцирования

(3.1) по у

и использо­

вания

(4.10'), связывающее третью и первую производные на

стенке.

Условие совместности

приводит к

появлению

параметра

 

аВ2

 

 

 

 

 

 

 

 

 

 

 

 

 

pv

б 2 .

 

 

 

 

 

 

 

 

 

 

(4.11)

П о м и м о (4.10)

 

можно

з а д а т ь

еще

 

 

 

 

 

дЧі

 

при у=6,

 

 

 

 

 

 

 

 

 

 

— = 0

 

 

 

 

 

 

 

 

(4.10'")

10 — 2274


а т а к ж е значения

производных

порядка выше третьего

на внеш­

ней и внутренней

границах слоя.

 

 

Последовательно

привлекая

все новые и новые

условия и з

серии (4.10), получим д л я F(r| ) и С7(г|) из /(т)) =F(r\)

+XG(r\)

в ы р а ж е н и я , приведенные

в табл . 4.2. П р и построении

f(n,) один

из коэффициентов

ряда

(4.3) Л,-, вообще говоря,

я в л я ю щ и й с я

функцией х, остается произвольным. Ег о можно ввести в функ­ цию U (х).

Т а б л и ц а 4.2

Fin)

про­

филя

 

1'

11(1- - ч

2

• ї ї ) 2

3

т| f i ­• ї ї ) 3

4ll f i ­ 11 ) 2 (1 +2г| )

5ll ( 1 - i l ) 3 ( H - 3 i i )

6

4 ( 1 - - i l ) 2 X

 

X (l+211 + Зті 2 )

7

• 1 ( 1 - т і ) 4 ( 1 + 4 п )

8

11(1- - i l ) 3 X

 

X ( 1 + З Ї ] + 6.І2 )

движения (3.1) в точке

 

 

К а к

видно

из

приведен­

 

 

ной таблицы, влияние элек ­

 

 

тромагнитных

сил учитыва ­

 

 

ется

лишь

при использова­

 

с (и)

нии

полиномов

не ниже ш е ­

 

 

стой

степени.

Ка к и в

пре­

0

 

дыдущей задаче,

эт о приво ­

 

дит

к

увеличению

объема

0

 

 

вычислительной работы.

 

0

 

 

0

 

Д л я

определения

двух

0

 

характеристик

течения

U(x)

П 3 ( 1 - і і ) 2

и Ь(х)

привлечем

уравнение

 

 

0

 

импульсов

(3.7)

и

уравне ­

1

1 ^ ( 1 - г , ) 3

ние

(3.9). Д ж а у г а ш т и н

[11]

 

 

вместо

уравнения

 

(3.7) ис­

 

 

пользовал

ви д

уравнения

максимального значения

скорости

у—ут-

Однако

в х о д я щ а я в эт о уравнение

величина

 

пр и при -

ближенном з а д а н и и профиля м о ж е т привести

°У

v=ym

неточ­

к большим

ностям

[10]. Так , при применении полинома

третьей степени чис­

ленный

коэффициент

в в ы р а ж е н и и

н а п р я ж е н и я

трения

(3.57)

становится равным 0,556 (вместо точного значения 0,221).

 

Подставив

(4.3)

в (3.7) и

(3.9),

получим

систему нелиней­

ных уравнений:

 

 

 

 

 

 

 

 

d

U4a{X)~-

U

аВ2

 

 

 

 

 

 

 

dx

-V—8

Ubb(X)

;

 

 

 

 

 

d

 

 

оВ2

 

 

 

 

 

 

 

(4.12)

U42c(l)

U282ab2(X)

,

 

 

 

 

 

 

dx

2 Р

 

 

 

 

 

 

в

которой

 

 

 

 

 

 

 

 

 

 

1

 

 

I

 

1

 

 

11

 

 

а=

J

P(r\)di\,

Ь = J

f(n)dn,

с=

j

( f2(y\)

J

fd4)

dr\,

 

 

0

 

0

 

 

0

 

0

 

 

 

a X з а д а н о соотношением (4.11).


П р и

использовании

полиномов,

д л я которых- G(т))

= 0,

ре­

шение

системы

(4.12)

можно найти

в

квадратурах, в

частности

 

/

Чс/—

1

\

12(За-1)

 

 

 

 

 

/

 

+ - VР ~ A ' )

dX=Cx,

 

 

 

 

b " ,

A ( 1

 

 

(4.13)

г д е а = ^ , р = - ^ .

 

 

 

 

 

 

Сходимость этого интеграла при Л->-оо обеспечивается

при

условии

^

> | .

 

 

 

 

 

 

 

§ 2. А С И М П Т О Т И Ч Е С К И Й М Е Т О Д

 

 

 

 

 

Представление

о конечной толщине

пограничного слоя не яв­

л я е т с я

принципиальной особенностью

приближенных

методов.

К качественно тем

ж е

результата м приводит

и представление

об

асимптотическом

х а р а к т е р е слоя. П р и

выборе

профиля

скорости,

асимптотически убывающег о к нулю на бесконечности и имею­

щего

один

максимум

в сечении

#=const ,

неизбежен

произвол .

Н а л и ч и е в определенных случаях

в н е м а г н и т н о й

гидродина­

мике

точных

решений

уравнений

П р а н д т л я

струйного

погранич­

ного слоя наводит на мысль выбрать из семейства профилей именно тот, который следует из точного решения, и затем ис­

пользовать его д л я

приближенного

описания струи

в

м а г н и т ­

н о й

гидродинамике .

 

 

 

 

 

 

Конечно, такой профиль автоматически удовлетворяет соот­

ветствующим граничным условиям

(за исключением тех из них,

которые

вытекают

из

уравнений

движения, например

(4.10")),

а з а д а ч а

сводится

к

выяснению,

в

какое

сечение x = c o n s t сле­

д у е т

переместить профиль скорости

и как

изменить

м а к с и м а л ь ­

ное её значение при наложении магнитного поля, чтобы безраз ­ мерный профиль остался тем ж е , что и в отсутствие поля.

2.1. ПЛОСКА Я З А Т О П Л Е Н Н А Я СТРУЯ

 

З а д а д и м с я точным профилем скорости, известным

из реше­

ния з а д а ч и о непроводящей струе:

 

U = Umf(y\)=Um^^—.

(4.14)

сп Зг|

 

ю-

 


П о д с т а в л я я

(4.14)

в уравнения

(3.1) и

(3.7), придем к

системе

 

 

 

 

 

 

 

 

 

 

d

3

 

 

 

(4.8), в которой следует положить

-=~7г>

Я = —18. Вычислив

интеграл в

(4.9), имеем решение

 

 

 

 

 

 

 

72\'СУ/з

 

 

 

 

 

 

 

 

 

 

 

б2 = . ,

r M

v

 

 

 

 

 

 

 

 

 

(4.15)

 

1 - C N r ' 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ст52

 

в котором

С — постоянная

интегрирования,

a N =

. Постоян­

ную

интегрирования

определим

из

условия, что в предельном

случае при N - й З

 

 

 

 

 

 

 

 

 

б2 +

(10,91)2

( - ^ - 2

V V / =

 

 

 

 

 

 

 

 

 

 

 

 

Jo

 

 

 

 

 

 

 

 

 

Это

дает

 

 

 

 

 

 

 

 

 

 

 

 

С =

1,655

 

 

 

 

 

 

 

 

 

 

 

 

где J0 — импульс струи, заданный в начале координат

х=0 и

полученный

при предельном

переходе при N-»-0 из (3.7). Соответ­

ственно д л я ит

имеем

 

 

 

 

 

 

 

 

Hh l =0,454 ( и у > - ™ * \

 

 

 

 

 

 

( 4 , 6 >

Теперь

можно

определить остальные

характеристики

струи.

Р а с х о д

в ы р а ж а е т с я формулой

 

 

 

 

 

 

Q = p /

udy = 3,302 ( — ) ' / з

( l - C N x V 3 ) ' / 2

)

 

 

(4.17)

 

—со

 

 

 

 

 

^

 

 

 

 

 

 

 

а импульс

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

/ = р J

u*dy = J0(l-CW>yb.

 

 

 

 

 

 

(4.18)

Естественно, результаты

(4.15) — (4.18) полностью совпадают

с точным

решением,

приведенным

в § 3 главы

I I I . Н е

будем

поэтому

более

подробно останавливаться

на этом

решении.


2.2. СТРУЯ У ПЛОСКОЙ ПОВЕРХНОСТИ 1121

 

 

 

 

 

 

Д л я

решения задачи, сформулированной

в § 1 главы

I I I , вве­

дем функцию тока из уравнения

неразрывности

(3.2)

и

опреде­

лим ее следующим

образом:

 

 

 

 

 

 

 

и=—^',

 

и = - — І ,

 

ср = г|)т Дг,),

 

 

 

 

(4.19)

 

оу

 

 

ох

 

 

 

 

 

 

 

 

 

 

где, ка к

и

прежде, r

l=

"^"> а

б (л;)

имеет смысл поперечной ши­

рины струи. С введением функции тока по (4.19)

уравнения (3.7)

и

(3.9)

перепишутся

к а к

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

оо

 

 

 

 

 

/ ^ • ^ ( ^ )

= - v ^ r ( 0 ) - N ^ / ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.20)

со

 

 

 

 

 

 

со

 

 

 

 

v

'

 

Д л я

определения

интегралов,

входящих

в уравнение

(4.20) f

зададимся

функцией

/(т)),

определяемой

решением

А к а т н о в а

[13] д л я непроводящей

струи:

 

 

 

 

 

 

 

/ ' = - ! - ( w = - f 2 ) ,

 

 

 

 

 

 

 

 

 

и, соответственно,

ar\=

^ 3/а^|/а

^2 .

 

 

 

 

 

 

Если при этом воспользоваться еще

значениями

функции

тока на

бесконечности

/\>о = 2,515

и второй производной в нуле

/"(0) =0,221, то после вычисления интегралов система

(4.20)

пе­

рейдет в

 

 

 

 

 

 

 

 

 

 

 

 

( - ^ У = - 0 ^ - ^ - 2 , 8 6 г ^ т

;

 

 

 

 

 

V

6

'

 

 

6

 

 

 

 

 

 

 

(4.21>

{ ^ y

 

= - 3 , 1 8 N ^ .

 

 

 

 

 

 

 

 

В

целях

упрощения

выражений

(4.21) перейдем

к функции

Г =

= i p m 6 . Тогда из

 

 

 

 

 

 

 

 

 

 

7-6'=2,22N 63 + 0,756v6;

 

 

 

 

 

 

 

^

r = l ,9N6 2 + l ,008v