Ламинарные и турбулентные течения полидисперсной газовзвеси с взаимодействием частиц и межфазным массообменом тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Наумов, Владимир Аркадьевич АВТОР
доктора технических наук УЧЕНАЯ СТЕПЕНЬ
Киев МЕСТО ЗАЩИТЫ
1993 ГОД ЗАЩИТЫ
   
01.02.05 КОД ВАК РФ
Автореферат по механике на тему «Ламинарные и турбулентные течения полидисперсной газовзвеси с взаимодействием частиц и межфазным массообменом»
 
Автореферат диссертации на тему "Ламинарные и турбулентные течения полидисперсной газовзвеси с взаимодействием частиц и межфазным массообменом"

Академия наук Увраины Кнститу? . гидромеханики

На права;? рукописи

Каумоз Владимир Аркадьевич

УД1-: 5Е2.525 .

ДАШШВДВ И ТЛСТШНЫЕ ТЕЧЕНИЯ ПОШЩКЖРСНОЯ ГАЗОВЗВЕСИ а БЗАИМОДЗЙСГВШЗЛ ЧАСТЩ И МЕЕЙАЗНЫМ МАСШШНОМ

Специальность. 01.02.05 - механика кидкоотея» ..

газа и плазш

Автореферат

диссертации на соискание ученой степени доктора' технических наук

РГ6 од

1 НЮН 1933

Киев 199Э-

Работа выполнена в Квдстккградском' техничгском институте рыбное промышленности и хозяйства и Институте проблем энергосбережения АН Украины.

Научный-консультант: догаор техничесшвс наук А.АЛРАгЗзЕР Официальные оппоненты: ;

■ . доктор технических наук | С.И.КРИДЬ

,у доктор технических наук, профессор ВЛ.ТИМСНШКО

/ доктор физико-математическ:х наук, профессор Ю.И.ШМАКОВ

Ведущая организация:

- Институт технической теплофизик« АН Утраты

■ Защита состоится " , Ю^5 года заседании спе-циапизировешого ученого совета Д О!.. ОЙ. 01. ПРИ Институте гидромеханики АН Украины з^_час. в зале заседаний ученого совета

по адресу: ■■ \ г \

252057 Киев 57, ул. ГеляЗсва, 8/4. " \

. С диссертацией ког-но ознакомиться е', библиотека жститута.

Автореферат разослан Х^"_" " \ 1593г.

!Ученый секретарь - ^ - \

. . специализированного совета . \

¡¡.ср.-**. ч.^/Ц, V0. Хм^от^'^у^оЬ

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Газовые потоки, несущие дисперсные твердые частицы или капли, широко распространены в природе (дзи-кение водяных капель в облаках и туманах, песчаные и пылевые бури, рассеяние- вредных примесе.: в атмосфере и др.), и используется практически во всех отраслях техники .(системы газоочистки и пневмотранспорта, проточные химические реакторы со взвешенными частицами, аппараты для термической и механической обработки.сыпучих материалов, распылительные сушилки и теплообменники, различные энергетические установки с двухфазным рабочим телом или с газодисперсными теплоносителями и т.д.). Для разработки новых :: совершенствования существующих технологических устройств, связанных с течением газовзвеси, необходимо' детальное математическое моделирование физико-химических процессов, протекавших з тагсгг течениях. Теоретическому и экспериментальному изучению различны:-: аспектов гидродинамики и теплообмена'газовзвесей посвящены работы многих исследователей. Вместе с тем ввиду исключительной сложности проблемы общая теория движения газовзвесей далека от своего завершения. Многие актуальные вопросы гидродинамики двухфазных течений изучены недостаточно. Так, турбулентные течения газовзвеси исследованы в основном в простейшем случае: при люнод::-сперсных твердых частицах без их взаимодействия и других физико-химических явлений. На практике часто турбулентное движение газовзвеси осложнено■коагуляцией и дроблением частиц (капель), фазовыми переходами и химическими реакциями. В исследовании таких потоков сделаны только первые шаги. Важно подчеркнуть' неаддитивность указанных .явлений, их нельзя изучать раздельно из-за существенного-взаимного влияния. В связи с этим проблема построения моделей турбулентных течений полидисперсной газовзвеси с сильным взаимодействием частиц, фазовыми переходами и химическими реакциями является весьма актуальной.

Цель работы: ' V

- создание основ теории потокой полидисперсной газовзвесз с взаимодействием капель и межфазным массообменом;

- исследование с ее помощью газодисперсных течений в струях, пограничных слоях' и каналах при взаимном влиянии турбулентности, коагуляции и дробления частиц, фазовых переходов и химических реакций; .

- разработка методов и-программ расчета на ЭВМ турбулентных

- С ~

газовзвесп пркг-знительно к устройствам современной тех-

цуг-г-я ксвизка. В «гассертацгонной работе '. Газрг'стаки основы теории движения полидисперсной газовзвеск е 'сххънам взаимодействием частиц, фазовые; переходами к хши-чезкпгл: реакгд^яиг:

- голу^ена система уравнений турбулентных течений полндисйер-сной газовзвеси с коагуляцией и дроблением частгц, меьфаз-нык кассообмезо!.:; . •

- предложена модель турбулентного переноса в сдвиговых тече-:-:газовзвеси с сильным взаимодействием частиц и глекфаз-ншл кассообменом, Еключашая'уравнения для вторых моментов пульсашюншх.скоростей обеих фаз.

Впервые в рамках.гидродинамического направления решен ряд задач о течении газовзвесн в различных условиях:

- в газодкспероном турбулентном пограничном слое (ГГПС), в т. ч. в соплах Лаваля; Г

- в газоди^персных турбулентных струях с переменным Фракционным составом, взаимодействием' частиц, их плавлением, испарением к горением;

- в канале с учетом дисперсий скоростей отдельных фракций;

- в газодиспзрсном ламинарном пограничном слое (ГЛПС) с учетом подъемных сил, действующих на частицы, их испарения, коагуляции и дробления, образования гадкой плевки.

2. ¡¡случены новые результаты в исследовании физических явлений, проксходящих с отдельной частицей в потоке (удар частила о Еероховат; ю стенку, движение испарякащейс.: капли, действие сели Магнуса на вращающуюся частицу и др.).

■ Пгактгч'зская ценность работы заключается в определении ос-коекш. закономерностей ламинарного, и турбулентного движения, тепле- и кассопероноса полидкеперекой газо?~веси при наличии коа-гулпдп: и дробления частиц, фазовых переходов и хешчеекпх реакций. Предложенные М£т51.:ат*£ческие модели к методы расчета могут быть использованы при проектировании энергетических установок п технологически аштзрзтев, з которых рабочим телом является по-лидкеперзная газовзвесь.

На основе диссертационных исследований были разработаны к Еяедргнн елгдущкз. пакеты прикладных программ:

- о -

1. Турбулентной высокотегятературной струи с плавязцэлисл частицами при газопламеннсм нанесении покрытий. - з Белорусском Н~: порошковой металлургии (г. Минск);

2. ГТПС з каналах пэременного сечения - з Московском институте теплотехники:

3.' Полядисперсного течения -с коагуляцией и дроблением капель з нолем массообменном скруббере - з ЕГО "Государственной институт прикладной химии" (г. Санкт-Петербург);

4. Высокоскоростного течения газовзвеси з. фурме длл подачи инертных газов и пероашов з расплав на предприятиям концерна "Азовмаш", Мариупольском металлургическом комбинате т.'.. Ильич;.

Екедреннне программные комплексы используются при проектировании изделий новой техники, а такае для ссверпенствсзанпл технологических процессов.

Результаты работы могут быть использованы в вузах при поучении механики многофазных сред.

Автор защищает: .

1. Основы теории течений . полидисперснсй газовзвсси с сильный взаимодействием частиц а меафазным массосбменем;

2. Комплекс математических моделей и методов расчета

чных типов турбулентных газодисперсных течений с коагулт^ю-й п дроблением капель, фазовыми переходами и зддспэскпмз реакциями;

3. Результаты расчетов течений газовзвеси в туреулентнпд струях, пограничных слоях и каналах в условиях 1аз8пгзсго влняннл протекающих в них физических процессов.

Апробация работы. Основные результаты диссертационной рз,': -ты были доложены и обсуздены на:

- I конференции по механика Академий нар: сса. згрзл. - "загг, ! 537;.

- II Минском международном форуме по тзшюмасзсесггэну,

- Ш Всесоюзном съезде по теоретической и прикладной Москза, 1991;

- I и II Всесоюзных конференциях ГТеплсфизика и гндрогазедппз:'.:-ка процессов кипения и конденсации". - .Рига, 1952, 1983;

- ХГЗГ, ХГ Всесоюзных конференциях "Актуальные вопросы аз -родисперсных систем". - Одесса, 1936, 1989;.

- Всесоюзной конференции по кинетической теории разрегзнных :: плотеых газовых смесей, и механике неоднородных сред. - .Тзнин-

- 4 - ■■

град, 193?; . •

- Г!И л IX Всесоюзных симпозиумах по горению и взрыву. - Таш-:-:ei::s 1936; Суздаль,-198а;

- Г.", Y £ YI Всесоюзк-х научных совещаниях по теоретическим и прикладным аспектам турбулентных течений.-Таллинн, 1982,1985,

• OQC

Публикации. Материалы диссертации опубликованы в 57 печат-:пл: работах (в т.т. монографии, переведенной на английский язык).

Структура и объем работы. Диссертация изложена на 275 страница;:., содержит 14Г рисунков, список литературы из 282 наименовании, ваего 348 страниц. Работе состоит из введения,'восьми раздало в и заключения.

СОДЕРЖАНИЕ РАБОТЫ

Во введении обосновывается актуальность проблемы, формулируется цель работы, кратко описывается ее структура.

Разделы 1 - 3 объединены в часть I - Основы теории. Первый раздел посвящен выводу системы уравнений .турбулентного .течения полидисперсной газовзвеси с коагуляцией и дроблением частиц. В пгдрззделе I.I приведен обзор методов построения математических моделей газодисперсных потоков; весомый вклад в развитие теории таких течений внесли Бусройд, А.Н.Крайко.С.И.Криль, Кроу, Карк-ер, Ю.П.Лунькин, Р.И.Нигматулин, Х.А.Рахматулин, Coy, В.В:Стру-miftckíií, Трусделл, H.A.Фукс, Ф.И.Франк ль и др. Анализируются работы Г.Л.Баоухи и А.А.Шрайбера, И.М.Васенина с соавторами, 'Джиле спи, Марбла, Л.Е.Стернина, А.П.Тишина с соавторами и др., посвященные моделированию течений с леременным фракционным составом дисперсной фазы. Обосновывается целесообразность гибридного подхода, сочетающего элементы гидродинамического и статистического описания эволюции состояния дисперсной фазы и позволяющего, в отливе от традициошшх методов, корректно учитывать различие параметров (скоростей, температур и т.д.) частиц данной фракции в зависимости от их предыстории.

В подразделах 1.2 и 1.3 сформулировано кинетическое уравнение для полидисперсного ансамбля капель с учетом мевфазного мас-сообмена из-за химичеыих реакций или фазовых переходов, и с его помощью получены гидродинамические уравнения для актуальных зна-

чэгзЗ тр5".зтроз састоянпя а-5 фр^ащз. В нщгезлегэ 1.4 к уравнениям дзуя££2нсго псдядзсдзрсного потока приманена процедура . Рзйгольдса я по^/чена система уразнэянй турбулентного течения ГЕБОвзвеся с коагуляте! я дроблением частиц, мзгфазнжа кассооб-

мзксм '

а<с > а

+--«р ><¥ > + (р"/'» = ^ > « о* >, (1)

ел е г * * " я дв

5<р > а н '

—,9 + «о ><•/ > + <р' У >) = •-£<.!>, (2) а г а г '3 3 3 » ..1 3' '

а а •

--«о XV > + <р'У>) +--«р ><У XV > + <Р'Ч"><У > +

ф £ ® Л "ев д Г я. в ' г» в и

. V <р хУ*У> + <У хр'У> + <р'УУ > + <0 ХО > +-<Р'0'>) =

г» ее в ~в с Ге вв .- * в в в

= <р ><Р > <р,р,> + <1 > + <*Г ><У > + <«1* V > , (3)

г в * гв в а ^в д« дю 3 '

а а

- «о XV > + <р'У >) -г —— «р XV Х'У > •+ <о'У х¥ > +

д '9 9 ГЭ а ' а г' в .■ 3 ' Э в . 3

• а<р> а<т > + <о хУУ>••+ ЧУ хр'У> + <р'УУ>) ■

зд В гд д г9 в 9

а г а г

- Т!(<о ><Р > + <р'Р'> + ^ ><У > + <.Г V >) , ' - . (4) -

• в я ~ в в са да д* с« 5

* £ ■ ^ » _ , * ^ _ ^

- а ' . а : "" , V ■ ' • . . •:

- «р ><а > + <р'1г'>) 4---(<р хй ><у > + <р'1г'хУ > + .

а г г * а г , "•■-•"'..•

+ <о х1г'У> + <Ь ><р'У'> '+ "<р,П,У'>.+ <р ><!• > + <р'*Г>) = ,-

'в а « в гв а гв в а Гв в ж '

= <о х(3 > + <о'СГ> + <П > + ' с-Г х1г > • + '<«Г'Ь,"У + -

• а а » & 9 а - да да . да 9а

+ 1/2<^ >С«У >- <У >)г+ <(У' - У')г>] + 1У2<«Г' (У - У)г> +

■. . д». * ... -я»—»' . э» 9» • ; ■

чсЛГ'-СУ* -У )>• (<•/ >-<У >>;'. -V : - (5)

дв э» . в' ч > д* » . 1 • , ...

а в :

N

+ <p xh'V> + <h xp*V'> + <p'h'V'>) = 2 [<p xF > • (<V >-<V >)+

•g g g .9 s rg g g . s ® з *

+ <p'F'>-«V >-<V >) + <pxF'-(V'-V')> + <F >-<p'(V'-V')> +

"в » g ■ в g з s 1 я g s

■ 5<p> ~ d<q >

+ <p'F' • (V- V' )> - <p ><Q' > - <p'Q'>) + —- - —— - +

° з u . 4» » . g t а г

>+ <x>'V<V:> > <%' ■ v V) -g Я 9 9

+ Í<V >--l<p> + <ÍV'--1 D'

I » зг> 19 a r-1 •

- £ C<J xii > т <J' h* > + t/2<J >C«V >-<V >)2 + <(V' -

93 с;* ge g» vja • g ga

3 — 1

- V')2>] т. 1 /2<J' (V - У'Г> + <J' ' (V - V )>-«V >-<V >)}., (6) g g» g* g g* g» g a< g '

a a

— (<p ><ь,> + <о'ь;>) + — л 9 1 . 91 a r

<— (<p ><ь,> + <о'ь;>) + — (<p ><ь, ><v > + <p*b*><v > т

a»- rg i 'gi a_ rg i g . . g i 9

a<q, > N

+ <p ХЬ,'У'> + <b. ><p'V'> + <p'b'V'>) --1 - 2 <J, > , (7)

ig irss rgig я „ _ r. i» '

- «p ><L > + <p'L'>) +—- (<p ><L XV > + .'<p'L'xV > -i-

з t V ' • a a r * * 1 *

<p ><L'V><L xp'V'> + <p'L'V> + <p ><C > + <p'C'>) =

' в ел 8 * 0 s ' Q & 8 r в Л ' a

= <p XK > + <p'M'> + <B > + <J ><L > + <J' l' > (S)

' я a r6 * «i 0** ал . gt. з» ? л '

a a

~ «Pex5a> + <p:a;>) + — «p.xo.xv + <p;e;xv.> * <p.><«:v:> + <e.><p:v:> + <p:e;v:» = <ss> <I/S<J"!> +

4 4/3<J9s> + <Jo>) + + 4/3Jce+ JjS, (9)

- <p> - pgc *n<%> e <v/Mt, x¡ <jle> = >; (Ю)

- где ; V - скорость; h. - энталышя; L - кдееавдесаспа моменг; t -

время; г - радиус-вектор; Т - температура; р - давление; р, р° -распределенная и истинная плотности; е - диаметр частиц; индекс g относится к параметрам газа, з - з-й фракции. дисперсной фазы; ~ - сила межфазного взаимодействия, отнесенная к массе частицы; С - интенсивность мекфазного теплообмена; И - момент сопротивления несущей среды;■ , , Пв, В=- члены, обусловленные коагуля-' лкей и дроблением частиц при их взаимодействии; .1 - межфазны..: мзссообменом; - часть .Г.,' связанная с изменением диаметра

капель 2 при их взаимодействии с мелкими каплями; J) - интенси-

з 1С

вность перехода массы из 1-й компоненты газа с. массовой долей Ь1 в 8-и-фракцию частиц; Удв ,1гдв, скорость, знтальпия, кинетический момент массы, совершающей фазовый переход £ —» з; й -универсальная газозая постоянная; •М1- молекулярная масса 1-й компоненты газовой фазы (1 =1,2,...,п); N - количество фракций дисперсной фазы; , - векторы теплового потока в- газе и потока массы 1-й компонента; угловые - скобки о означают осреднение по времени, ' - пульсационная составляющая. Т.к. объемная доля дисперсной примеси ар мала, тензор касательных напряжений в несущей среде определяется как в чистом газе.

Осредненнке уравнения движения газовзвеси, кроме корреляций, отгасываюЕззс турбулентный перенос массы, импульса, тепла з потоке

и др., содержат корреляции, обусловленные отклонением скоростей, энтальпий, кинетических моментов дисперсных частиц от средних для з-й фракции (эти отклонения обозначаются "), 0 = УУ", % = 1г"У", С = ЬТ", где верхняя черта обозна-

' . в в « ' в а 6 в в в * - г

чает осреднение по частицам данной фракции. Корреляция 0и обусловлена разными скоростями частиц данной фракции в' некоторой точке пространства из-за разной предыстории. Корреляции вида <У!|У^> обусловлены отклонением в некоторой точке' потока мгновенного значения, в данном случае скорости дисперсной газы, от осрзднен-нсго по времени. Если предположить, что все частицы з-й фракции имеют в данной точке одинаковые скорости, то Оя = 0, & <Ф О, т.к. причина поязлснкя последней корреляции - в турбулентном характере потока. В ламинарном течении <У^> = 0, а 0.

Второй раздел посвящен основным физико-химическим процессам, происходящпг.! о частицей (каплей) в потоке газа. В подразделе 2.1 приведены выражения для сил, действующих на частицу: аэродинамического сопротивления, термофорзза, подъемных сил Саффмэка и Ма-

•гнуса.. Представлены рвоульйаш оазхоа по опрэдзлэЕЯ» те С в.фор&улв Рубшов'а-Комзра да скш Ыашуса

V-1/6 см р;68 >< V, V V* V = % -л*(11).

где О , 0,- векторы угловой скорости частяда и окрргагсего ее элемента газа, соответственно (при йеыэ|Пи|5г/у—• 0, Ие з|Уа

о Сы= %). Ка лабораторной установке измерялось отклонение падающего перика при различных начальных угловых скоростям. В результате обработки опытных данных савиэстао с денная другзк авторов получека формула

1; при Не^< 1 СО ,

о, яп> . " * ^

6,05/Es ; при 100 ^ Но < 3,6-10

В том кз подразделе вычислена величина реактивной еяза, действующей на сферическую каплю, испаряхщуюся з неоднородном ламинарном потоке собственного пара с Линейным профзлеы температуры (на внешней границе слоя Т ~ Тш).и достпгсуз темпзра-гура ш-

пения т ... •'

з

гв • % б

Зр'е

~ г»

»a-v*

JrpJ Cív -,IS1.S) (1 - Y/a).,. (13)

где У - расстояние до стенки, а - .шбйицаеет тетглоотдзчя поверхности капли, fw- температура стенки,.а - толщина пристеночного слоя4 Iй- скрытая теплота парообразования, рпг - плотность пара у поверхности капли.

В подразделе 2.2 приведены выражения для интексинзоотн тепло- и массообмена. одиночной частицы в потоке газа, б т.ч. при испарении и гетерогенной .горении.

Взаимодействие частиц (капель) с твердой поверхностью ц мэ-зду собой рассмотрена в подразделе 2.3. На основе сравнения результатов расчетов с опытныгш дёенхгя' прп поыоци наглядного' графического метода' Рауса показано, что при ударе сферической частицы о шероховатую стенку с небольшой скоростью в точешз всего 'времени взаимодействии вшшняэтся закон трения Кулона. При ssa-^чктельнну^скоростях, Еесферичноетя частиц, Ее^ашх .размерах шэ-, ржоватостк' говераноск: следует пользоваться', гипотезой о поста-

Р

ялстзе (для данных материалов часта* к стешс?) отношения танген-оадышх скоростей точки касания Б пссле U* и до удара и , 0^/0r= const. Разработана блок-схема для-игрецеления занекаязг»-с:'лх характеристик несферических чзсткц после /дара о шероховатую стенку. 5 тс:-: т подразяелз приведен:: формул! лля опредеде-::::?; кинематических параметров Бзаккодрйст5у:гя;:х частяд (капель).

еодер:::нт рэоекяя дгу\ задач в прибдпгопя'. единичной пст^ци: 1) честна .во зрздаюцейся полости; С) частица в слое Кузтта, унесенная с поверхности. Ка рис. 1 представлены за-знек'ости безразмерного времени достхкснйя частицей Бкоггней сте-itkil полости t,Q (ft - угловая скорость грхлшш 1юлостй) от расстояния до оси г = й /R (R - рэдаус полости) к момент t = О для различных Stic, = бгП/ (i 3vxp), рс/рр°. Величина • tx быстро падает при увеличение г . Аналитически -и чпелзкзо исследованы еозкоккооти пренебрежения отдельными силаig:, действукгцхмк ка частицу во Бргцакдейся полости. Для слоя Куэтта установлены условия, при которых унесенная с твердой псьерхностп частица не вернется на нее.

Третий раздел гссгяден разработке модели турбулентности для сдвиговых течений гззсдзвеси с кассообмеком, коагуля-

цией к дробленом частиц. В подразделе 3.1 приведен сбзор мето-дез модедцроваяпя турбулентного переноса з газсдисперсных потокам. Ог:.-зчается вклад а их развитие F.К.Абрамовича и Т.А.Гпршо-впч, Ю.А.Зуегича, Л.Б.Гавина, З.Р.Горбиса, Л.Я.Зайчика, М.К.Лаа-тса. Zraue. Эльхобапп и др. В работах указанных авторов примесь - ösrar?5CKS .'¿окодисперсна, взаннодейехгзе частщ отсутствует'. В не;.'.;Юг.»Е работа:«:, посвященных созданию ¡¿оделдй турбулентных течений газовззест: о коагуляцией и дроблением капель- (Ю.В.Зуев и H.A. ле!7сн~:лск::Л, С.В.Ь&дзедев), используется упрощенная схема образования осколков, параметры бзажодействия капель определяются по продольным осреднением скоростям фаз, корреляция, связанные с коагуляцией, не учитываются, отсутствует кэгфазннй массообмен.

• 3 работах .1.11.Зайчика с соавторами и Л.В.Кондратьева доказано, что длл корректного описания турбулентных потоков массы, газульса и тепла дисперсной фазы в существенно неоднородных по- • токах (например, у стенки) необходимо -использовать уравнения пз- -рзноса зторых моментов не только.газа, но л частиц. В подразделе 3.2 получен;? дифференциальные уравнения для вторых моментов пуль-ездкошае скоростей обеих фаз. Уравнения для турбулентной энер-

гпп газа II н квадрата пу.чьсаиконноа поперечной скорости дисперсной фазы имзкт вид-

д к ,д }: д г <Э К р' .

- + 2 -Е — " 4- — )>| -

<3 V д - д I. ; 6 Ъ. ■■ д .к . * о -

> а V а V

- 2<7'.7\>-£гО,-£-г-£,, 5 = т " < -^ —> (14)

<3 X " р эх о X '

1 4 .Л

<:;;-> £ = £ £о'7' г'.> + <о ><?'.У > - ср'У х? .>]. (15)

3 Р . 'г с, 1 в I ' я 3 1. ( л % х, з г -

<о > в, = 7 [<7'. (V . - 7'. )> + ><?'. (V .- 7' )> +

' 5 -3 51 дв дя д1 дв доь д>

+ <7'. «Г >(<7 -.>-<?>)]. (16).

в' г»' <?.<- ' .

о<7'г> , а 7'г д<7 > а<7'гх, <0 >| -2- н- 7! |<7' . -2. > + 2<У?7'. >--5, 4- <7 .> -У +

в I с? г Т 1 51 5 1 ■ э х. " а х. -1-'

1 ! J

а ¥,Г с;<7 > г 3 7'-* ,

- О'--2- > + 2<р'7'> 4- у1 !<о'7' . -2. > + 2<7 .> 5<0'7' > •

•• 3 I " * д X Т 31 ч I • »

а<7 > о ' а<7

----± + __1 >» + 2<р'7'7' .>----= г + + а- . (17)

О X. в а х 31^ г •

; 1 '

=2(<р XV' > + <р' 7'ХР > + <р' 7' Г' >) . По)

е 1 С ? «?у • Й £ £5 у г8 « ' * ,

я- 2 (<!' V' > - а'у-хУ > - <3 ><7*г> - С.Г '7 ,г>) (19)

£ я з в« а * • 5 ' '

= 2(<1' 7!> - <«Т' 7*XV > - <3 ■ XV'2> - \7,г>). (2С*

дву е аз э я дз ' с ;*е в

Урпвнс~5я (14), (17'( стдеч&йтся 01 аналога'Пзкх урагнекл."!, получениях други?® авторами, уч-зтсм ме?;фазного масссо?ме:-:г (сл&-геекые £т), коагуляггли и дро-Злепия капель (член ).. Слагав-¡ста е., ? обусловлены г.:е:лр2зной силой, б - генералкой турзулен-пгосхи в следе за частицей. Получено также уравнение для е.:' • Представление остальных коррелят в уравнениях (1 )-(Э) рас-

'.смотрено в подраздел» 3.3. Пренебрегая в дифференциальных уравнениях для корреляций <tTV'>, <h'VV> конвектквнгаля и -диффузлон-кыми членачн, получены выражения ■

<Ц* Ч'> = - v ■ а<и >/ ó Y + Ф <"J '/V . (21

3 s st 5 ' "> - 9 s *

j 2

<7 > 7 7

ГЗ . < V 1 •

V -----5 =

1-7

Г5 Ч ' SV

'ах

<n' V'> = - Л. / (<0.>С ) díh.V <3 Y + ч>; V >

9

- .. „. с /с

<Р.Ч 7.т+ Тау : 7.т+ 7яу 1 7.,+ % 7sy + %

1X2

7 =-у + 7 , . 7 = 7 , 7 = р (1 г а Р.г +• a,Re" ), а = 0,1 '79,

'вх "а 'ео 'а у 'а 'с s í s 2 -з 1

а = 0,013; 7 = <з (0.5 a Re* + a Re*), Re*= |<V„ > |S /v,

2 ' 'eo s * i s 2 b s 1 rs 1 s ?

где vsi, - коэффициенты турбулентной вязкости и теплопроводности дисперсной фазы; Са, Ср - теплое?,кости газа при постоянном давлении и материала частиц; ц>ку, Ф„ - показатели затухания корреляций .параметров газа вдоль траектории двякеная частицы; 7=-т -величина обратная времени тепловой релаксации частицы s-й фракции; 6 = 18-уЛ /бг.

' з р а

Формулы (21), (23) обобщают гипотезы, ранее попользованные другим азтората. Первое слагаемое б (21) зависит не только от профиля осредкенной скорости дисперсной фазы б данной точке, но и от предыстории, т.к. v (а также Хз1 к коэффициент турбулентной диффузии частиц вычисляется по квадрату пульсацданной скорости частиц, a <Vs2> находится из дифференциального уравнения (17). В связи с этим указанные формулы могут быть использованы и при моделировании неоднородных турублентшх потоков, в т. . ч. пристеночных течений газовзвеси. Второе слагаемое обусловлено локальными явлениями. Для неинерционшх частиц 7 =¡ 7з = ► со, ... фя—* 1» vol— 0, тогда <U'V;> —> <и<У>. Для инерционных частиц 7ВХ> 7„у 0. Фв —- 0; коэффициент vai определяется отношением двух малых величин- <V¿2> и (7ах+ 7„у )> поэтом^- он кокет быть ко-

г

нечным.

Показано, что коэффициент .турбулентной вязкости газа в выражениях

- <1ГГ> = V д<и >/ЭУ, - <Ь'Чг> = у./Рг. а<11 >/ЭУ (25)

д д * д д д <■ I д

целесообразно вычислять по гипотезе Н.И.Акатнова, согласно которой для (к-е) и к моделей турбулентности можно, соответственно,

записать • . • '

уг = с^ к*/(е + ер)} ' (26)

1/2 а/х

■ V = си к 1/(1 + е 1 /к ). (27)

I (Л Ь р I

(27) следует из (26), если положить е = Спк3/'г/11, Со = 1, макромасштаб турбулентности, С^, Св- эмпирические константы. Формулы (26),(27) отличаются от обычно используемых т°г= С^/е, •у"^ = С^к1'2 тем,- что учитывается не только диссипация турбулентной энергии под действием вязких напрякений в, но-и диссипация за счет силы меафазного взаимодействия ер. На^рнс. 2'пред-ставлены рассчитанные в турбулентной струе профили ■ ^/(П^й.) (кривые: 2- 3- vt по формуле (27)) и опытные точкн 1 (Олек-хауз и др., 1387) при 2/(1=20, Пд_- скорость газа на оСн струн на срезе трубы диаметром й. Видно, что результат расчета по формуле (27) лучше соответствует экспериментальным данным. ■ '.

В подразделе 3.4 для различных условий вычислены осреднен-ные по времени интенсивности мекфазного массообмена и исследовано его влияние на пульсацнонные скорости обеих, фаз. На рис. 3 в качестве примера показана зависимость 6^.= ер+- eJ 0,5ег.гЕ./к, тЕ= 1/фху) в уравнении (14) от относительной инерционности монодисперсных капель при их испарении в изобарическом турбулентном однородном потоке газа; Рг = 2/3, зе = <рр>/<рд> = 1, В = Сп(Т -- Т,) / Iй, равновесная температура испарения при данном давлении, Сп- теплоемкость пара. При ртЕ —» о сильно инерционные капли не откликаются на пульсации скорости газа,, поэтому е^ —► 0. При —* оо пульсационная энергия дисперсной фазы такая же, как несущей, вновь е^—► О. В этих предельных случаях влияние дисперсной фазы; в т.ч. и массообмена, на Турбулентную энергию несущей среда мозшо не учитывать. При- промеяуточных значениях р*^

*

больше значения В приводят к•увеличении ■ £г за счет возрастать интенсивности перехода слабо турбулентной массы дисперсной фазы в несущую. При еток уровень к понижается.

В П2ЖЩШЬ5Э.-1кь> найдены осреднэнные по времени интенсивности тзшюмассооомзна и обмена импульсом при коагуляции к дробле-тша капель. Налраизр, составляющая связанная с изменением

деашзтоа капель б при их взаимодействии с меньшими каплями, он* в А

ределяегся .

- з<р > / (2р«е=) "ё* г фг <р>.(1 + на ) (28)

в- ' £ ' р в г в г в 1 г ге ;

г ш±

- = (в-- + О Г Э 7 , V ^ {<"/ — V >|. (29)

гв - ^ г а' г« гв га » г »•" х • '

«рв><р;у;> -<рг><р;у;».<уге>/(у^<рг><ре>), <30>

где Згя - коэффициент осаждения мелких частиц на крупных; Фгв -отношение изменения массы.мишени размером бе к общей массе попавших, в нее снарядов размером бг(Са > 6J. Проведенной анализ показал, 'что в течениях типа пограничного с..оя Иагв<< 1; величины Js, Пь, 1в.мошо вычислять по осреднепннм значениям параметров (подставляя в формулу для 1в скорости <Ve>*= <VE>+<p^V^>/<pi>).

.. Получено выражение для зедзникы , характеризующей измзпе-: кие <7^2> при взашодейст! л капель (см. уравнение (17)); в частном случае, при полной коагуляции оно икеет следующий в ид

" %/<Рв> = з/.(р°е ) кг.<рг> (<v;v«> - <v;2>). <31 >

.Ев рис.. 4 показано изменение отношения ф = от числа Рей-

нольдов' Rere = У^б/v при <pr>/<pg> = 0,125 и различных 6гв-' = Qr/5s. Видео," что при реальных для турбулентных потоков значениях параметров величины Ф. и Фг одного пгрядка'. Например, Hers= = 20 в воздухе при 5 = 150 мкм и V =2 м/с. При увеличении Re

tv G Г» А rs

и ¡уменьшении &гт величина $6 возрастает из-за усиления интенсив--носр: ЕзакмодейстЕия капель и коагуляции. При Heri_— 0, Sri>—» 1 взаимодействие"* отсутствует, —0.

.Разделы л 4-^8 составляют часть II - Прикладные' задачи.. Чет Ее-' ртый раздел, посвяаан высоко скоростник турбулентным течениям г&-.;.зовйвэск с твердым пайвдасторсяаа материалом в канале. Б родраа-; деле 4.1 'в''одномерном чисто /пщромв&тческои прибликекки полу-'

- 15 - • • / .

чэна система уравнений указанного течения без шафазнзго ■. кассо-обмена, коагуляции и дробления частиц. Анализ'граничных условий на стенке к уравнению сохранения вшульса дисперсной фазы для турбулентного течения з трубе позволил шчислить величину коэффициента трения (коэффициента.гидравлических потерь твердой при-кео:;). Pia рис. 5 показано, что результаты расчета течения газо-езезси u канале по такой одномерной модели хорошо'согласуются с опытная: точками (Вебер, 1988).

В подоазделе 4.2 предложен метод расчета высокоскоростного течения полидисперсной газовзвеси в канале при нестационарном тзплоподзоде через стенки. Для учета тешгаподзода к по.току газовзвеси вместе с системой уравнений,, ошсыващих ее течение, решается уравнение теплопроводности цилиндрического канала. IIa основе метода стрельбы разработан алгоритм, позволяющей проводить расчеты и при дозвуковом истечении штока из трубы, и при звуковом. В подразделе 4.3 проведены расчеты течений', газовзвеси применительно к продувке на металлургических производствах инертных газов и порошков через фурму. Такой процесс - неотъемлемый элемент новых методов получения высококачественных сталей. Фурма представляет собой трубу з защитной футеровке, опущеннуз в рас-плазленкый металл. . : .

На рис. б\7 представлены результаты расчетоз при длине фурма Н = 4 м, внутреннем диаметре & = 20 мм, длине части фурмы в расплаве 3,5 м, противодавлении 3,57 ат. На рис.6 показано изменение давления вдоль канала при различных объемных расходах газа G^ и порошка (ар- отношение массовых расходов дисперсной и несу-ущей фаз). Видно, что при неизменном ßg и возрастании требуется увеличить рй (кривые 3,4). Еа рас. 7 при 0,165 "м3/с, жр= = 2 представлено изменение скоростей газа (нулевая крпзая) й разных Фракций частил вдоль канала. Штриховая линия - результат расчета скорости баз учета соударения частиц.

Расчеты показали, что на металлургических комбинатах из-за неправильного подхода к конструированию фурм для продувки расплава в ковшах от 50 до S5 % энергии скатого газа теряется в се-гулируншх устройствах. При уменьшении -дааиьтра фургн и увеличении тешература ее .стэнки потэрн энергии схатого газа ейзгказте-ся, она расходуется в основном на увеличение импульса струи, что валю для перемешивания расплава, рассредоточения газа я порззез в объеме металла.

ct к

а /

1 ор / X

« ./л i s* \ /1 ^ V1

Г" —---- х \\ 1

6-

i «/с

4 ко

IZO

ÍDO

¿00

KM

Рис. 5 U, н fe

к • n К ï .»?

{ 0,QS5 J »

s a, H 3 !

s Q,1 es 2 §

h O/Cà' «

\ \

1

ÍOO

YS

SO

/ 4 /71 / / /!

<! 0,1

2 0,1

3 0,4

4 0,&

S vÁ

О - f V 2 5 o i 2

. РйС. б -' РиС. 7

Х,М

- 17 -

Дятай раздел посддщзн исследованию го.шдпоперсных газока-тт ESicifl з какала. В подразделе 5.1 проведено мпсле:.тое исслэдог-эйне полидисперсного ггзокапгльиого течения о коггудяпа-er¡ и дроблением с по;лощью слете»::? урденений, полученной з чисто ГйдродпЕатзпческом одномерном псдблипенит:. При эток учнтнзготся дробление капель не только при их дзан^одолзтдии, ко п газсг-'нг-кпческое.

В подраздел? 5.2 дли списания дгпзонпя дисперсной фаза написана систз;.та уравнений. которая следует к:: (t )-(íQ) п одпснл-р-изотор:«д1ч?скоа прибзп'кзнш: и учетпгеае? различие скоростей капель, пр;-каддена:;;их одной <тракцпи,

-1- (р, и ) = J (р ü S ) = 6 (1/3 У*' + J ) (32)

dX ■ 5 8 й 7.

— Со (и "+ о')] - о íü - U )/т + I т р (33)

q ' s 15 з ' ' <s ' д s s у ' £ /

rf - г 2 . ЙС,

а (р ü а ) = -- г с р U/% + .—3 + Э

■ * » « ' * I-S (, d ^ )

а

(3<;)

где г? ip* - дгясп-зрсия скорости з-й фракции, тв - гремя динамической ролакспцки частиц s-й фракция, g>. -- проекция ускорения свободного падения на ось канала д. Л - интенсив-кость изменения вздкшзн р о| з 'рзсультате коагуляции w дробления капель.

Получено выражение для D.; например, составляющая из-за'пополнения s-й фракции 'осколкамп имеет следующей вид

С 3 > 3 ь ¡"r г - 4 г г

D = —— 7 — V - К + (1 - Ф ) а | с +

s о , э ^ q qr .¡г чг» i qrs.

- r = "1 Q 4 ■

+ (ü - и ) 1 . . - (35)

' • j >

где ü , ог - сралняя скорость л дисперсия скорости осколков размером 5 • прп "ссг'лдадействип q-снарчдов и r-кшеней. Из (35} видно, нто D¿3> ¡завися? не только от дисперсий скоростей, но и от квадрата разности средних скоростей. Проведен анализ оздиокнос-та от различпгх факторов штзнсгзности Бзатдсодейстзил <-зстзд одной фргхспкя, rajErac разные скорости.

На ряс. 8-1Q представлена результата расчетов нисходящего •водоеоздуаного течения в гэрткк&шюм к&калс при следующих зке-чо'&ях параметров г. верхней сечении: U.0~ 1 м/с, о.0= о, TJc0- 5 расход вода - 7 ms '(¡«'час) . Исходный фракционный состав капель: 1 е. = 50 ЫКМ. 2 - ICO, 3 - 200, 4 - 400, 5 •- 800, 6 -1500, Кодера кришх на рисунках соотвотствуат номера?/, фрак^Ш..

4;i'V¡eh:;jo исследование показало, что г.эла-'кна о* определяется в ochoí"-:íом --гшерацкеЗ за счет оскс^ксг., приходящих в з-ю фрахщis. ¿; ¿шссилатяваым членом, обусловлю пм:-.: ск^оЯ метанного .(2 - ? р oVr^) в урашоши (34). Интенсивность сбр&зозания крупных оскодк-лз (е данном расчете O.S ш) мала, поэтому в крупных фракциях о* близка к нулю. Вместе с тем количество осколкоз, поступашпх в мелкие фракции, велико, поэтому о; ьамзтно возрастает по X. -к для' сашх мелких частиц суце-ственную роль в уравнении (34) играет дцссипативннй член Z_. Из-за этого на рис. 8 кривая 1 проходят ниже кривой 2. В соответствии с рис.З величина о. мелких фракций не только имеет такой se порядок, что к средняя скорость, ко монет превышать ее в 2 pasa.

На рис. 9,10 штриховые кривые - результаты расчета з чисто гидродинамическом приближении (o s 0). Как видно из рис. 9, ско-. рость 6-й фракции под действием силы тяжести увеличивается вниз по колонке.-Сила аэродина:/,ического сопротивления, действующая со стороны встречного потока' газа па мелкие капли, значительно ■ больше, поэтому на начальном участке кх скороот:-; ггад-ят (линии 2-4), концентрации - возрастают (рис. 10). По горе увеличения разности скоростей фракций начинает сказываться взаимодействие капель. Вначале преобладает процесс осазденпя мелких капель на крупных, в соответствии с эти?,! концентрации мелких капель (линии 1,3) интенсивно убывают от X и 0,5 м до X <и 2 м. При болышсс 1 этот процесс компенсируется образованием осколков при столкновениях каггель, и pt- ра меняются незначительно. На средних, скорое- -тях и концентрациях крупных фракций (кривые 5,6) учет сf практически но сказывается, т.к. для них oj; очень мала. Расчет средних скоростей и ра мелких капель в гидродинагличеасом приближена!, как видно из рисунков, приводит к существенной погрешности для течений в вертикальной колонне.

Подразделы 5.3. 5.4 посвящены разработке модели газокапель-еого течения с коагуляцией и дроблением в противотоке и применению ее для расчета очистки газов от вредной компоненты в полых -

Ц. . м/с

rte,

?.. ; кг/n3

i П 2 ! J í 4 6

8,HW 50 seo tsoo

/ / rh '/V

< o u- Э

o s У

y s

S

Рис. 10

' Рис. 11

кзоеооЗменнах скрубберах. При етом капли подаются сверху, а газ снизу, U < 0. В связи с тем, что з рассматриваемых течениях мелки капли могут выноситься вверх встречным потоком газа, каждая s-я фракция разбивается на две группк:

1) капли, шзодке скорость, нанравленнуэ вниз, и,_> 0;

2) капли, имзщае скорость, налразлзннуз вверх, U < 0. Уравнения, в отличие от -(£2) -(34), осредкяются по каплям,- при-надлззазгш указанным группам, и содержат члены, обусловленные переходом капель из группы "ннсходкщих" в группу "восходящих". На рис. 11 сравниваются опытные данные (Э.Я.Терат, 1Э76) и результаты расчета Обычного коэффициента кассоггзрздачи, равного Kv= = BJU, j/H, Бк = - Ia(bK/be), где br, "оо- концентрация зредкой пршгэси в зергнзм сечении колонны (при выходе rasa) и з никнем (при входе), соответственно. Высота рабочего участка скруббера -3,3 м; плотность орокения 9 м3/(мгчас), скорость капель в исходном верхнем сечении Uj_= 8,6 м/с. Штриховая кривая - расчет без учета коагуляции и дробления капель. При росте скорости встречного потока газа примерно до 9 м/с сказываются в оснозном два фактора, улучшающие массообмзн. Во-первых, из-за увеличения разности скоростей возрастает число Шервуда, вп-вторих, усщщшд?» ся дробление при взаимодействии квдедь g образованием моадцс соколков, т.е. увеличивается-поверхность массооемена, За счет первого из . зтнх факторов увеличивается Kv и при расчете без учета взаимодействия капель (штриховые кривые). При увеличения скорости газа более 9 м/с наряду с указанными факторами все более заме тнув роль играет уменьшение времени пребывания капель в копон-

Из рис. 11 видно, что расчет без учета коагуляции и дробления капель существенно занижает величину Kv Однако и сплошная кривая ланит' несколько нике опытных точек. Видимо, здесь сказывается наличие, в экспериментальной установке зон завихрениям рециркуляции, rasa, вызывающих усиление иасссобмена. , . ' Раздел 6 посвящен исследованию ITC. В подразделе 6 Л из (1) -(10) в приближении гурублентгой'струи получена система уравпв-. ний оСесимметричной ГТС с переменным фракционный! составом и мек-* фазным массосбыеном.

, /-г'. В, подразделе 6.2 описан метод расчета ГТС и проведено спав; не^ве результатов • расчетов с опытными данными для изотермических . г монодисперсных струй, деленные исследования показали, что одно; параметрическая k-иодель с использованием формулы" (27) для вычи-

слэкия V дает при. сравнении с экспериментальными данными результаты не хуке к-е модели; значения <7^2>, рассчитанные с помощью ' уравнения (17) и вычисленные в локально-однородном прибли-кении, различаются в основной части струи не более чем на 15-20 %. Влияние метода расчета <У'2> на осреднешше параметры - менее 3%. В дальнейшем в данном разделе уравнения для б и <7^г> не используются.

Подраздел 6.3 посвящен экспериментальному и численному исследованию неизотермической ГТС с плавящимися частицами в установках для газопламенного нанесения покрытий из порошков на поверхность поврежденных деталей и изделий'. Процесс напыления покрытий газопламенным методом включает в себя нагревание и ускорение дисперсных частиц в факеле сгоревшего газа и последующее нанесение этих частиц, нагретых до "расплавленного .состояшя, на предварительно подготовленную поверхность.

Экспериментальные исследования проводились на установке газопламенного напыления Белорусского республиканского.НПО порошковой металлургии. Температура в струе измерялась-с помощью термопары, скорость - измерителем скоростей светящихся объектов. Профиль потока массы дисперсной фазы оценивался по шлифу нанесенного покрытия на образцах с помощью микроскопа. В качестзе . дисперсной фазы в опытах использовались узкие фракции порошков с формой частиц, близкой к сферической.

В связи с тем, что длина рабочего участка струи при газопламенном напылении много больше ее радиуса, область влиянйя обрабатываемой поверхности на газовую фазу сравнительно невелика. На параметры частиц поверхность оказывает очень'слабое влияние из-за юс существенной инерционности. Задача решалась для области течения, в которой завершилось горение и сфор\шровалась изобарическая струя. Сделанные допущения позволили использовать систему уравнений неизотермическбй ГТС, полученную в подразделе 6.1. Доля массы частиц з видком состояний в данной точке течения при численном исследовании определялась с помощью величины

°> ТР < ти

^ = -ЬР3)> Vе V (3б)

11.

?де 1г*в-, ^ - энтальпия дисперсной фазы при- теипемпературе ила-

влекия в твердом и кидком состоянии, соответственно; 1т*-п*а-= Iм,1й - удельная теплота плавления.

На рис..12 сплошные'кривые - рассчитанные осевые скорости газа, страховые, - дисперсной фазы, точки - скорости частиц, -,оп-ределекнке экспериментально при следующих размерах: 1 - б = 30 • мкм, 2-80. Здесь расход горючей смеси - 0,75 г/с; медного порошка - 1 г/с; координата X отсчитквается от сопловой головки. Наблюдается хорошее согласие опытных точек к расчетных кривых. На рис. 13 показано изменение величины вдоль оси струи: 1 -б = 30 мкм; 2 - 60; 3,4 - 80. Кривая 4 получена при увеличенном в два раза расходе порошка (т.е. 2 г/с). При чрезмерном увеличении расхода дисперсной примеси (кривая 4) снижается доля массы частиц в расплавленном состоянии,, что отрицательно влияет на качество напыления. К аналогичному результату приводит уменьшение расхода горючего газа при неизменном расходе порошка. Созданная на основе'данной модели программа расчета ГТС с плавящимися частицами внедрена и используется для совершенствования технологии нанесения покрытий газопламенным методом.

В подразделе 6.4 проведен расчет полидисперсной турбулентной -струи с гетерогенно горящими частицами. При зтом процесс горения рассматривался на основе* предположения сб эффективной кинетике (реакция между горючим и- окислителем полагалась односту-•пенчатсй). Продукты реакции имели' температуру и скорость частиц, теплота реакции шла на нагрев частиц, а затем за. счет мегфазкого теплообмена на нагрев газа. В качестве примера приведены результаты расчета ГТС с горением углеродных частиц. Показана принципиальная возможность расчета горения пылевидного топлива с помощью модели неизотермической ГТС с ыегфазкым массообкеногл.

В подразделе 6.5 представлено сравнен:» результатов численного исследования полидасперсных турбулентных струй при наличии коагуляции, дробления, испарения капель с "опытными данными. На рис.14 показано изменение диаметра Заутера вдоль оси струи. Точки - опытные данные (Дкул, Ераут, Унгут, 1533), кривые - результаты расчета (штриховая - при полной коагуляции взаимодействующих капзль). Видно,'что расчет без учета дробления капель сущес-.твеЕно завышает их средний размер. И в опытных данных, и в расчете медленный рост б32 на начальном.участке струи затем сменяется, значительно более интенсивным, это объясняется данными рис. . 15, на котором представлено изменение на оси струи продольных

<Щ>,<ир>,н/с

- or¡t

55

ÍO

ä Jr / г- □ i О 9 \ \

/ / î^l /

0,5

h 1 I tZ "Ч 5 ч \ i .

SO {SO X,Mñ -Óо /50 250 К}нц

Рис. 12

50

26

/ / л О

i Г; / , /

/ i- . / / 1 О ° i

Рис, 13

9 20

/0

■' / </ /о % & ^ s о i Ä 2

/ / а/ / J il

■/г? д

а

¿0 ... 400 Х,мм О Рис. 14

SO <00 Х/-М

JO

- 24 - '

■скоростей двух фракций <1 - с начальным диаметром б = 15 мкм, 285), сплошная кривая -.скорость газа, точки - опытные данные. На срезе сопла.капли всех фракций имели низкие скорости, близкие по величине; интенсивность их взаимодействия была' невелика. Из рис. 15 видно, что по X скорости мелких фракций возрастают значительно быстрее, чем крупных; разности |<ие>-<11г>| растут по X, из-за увеличения величина ,б32(Х) быстро растет. При Х>110 мм продольные скорости разных фракций сближаются, поэтому темпы роста 5эг(Х) замедляются. .Численные исследования ГТС с существенным испарением капель показали, что расчеты без учета взаимодействия капель приводят к значительному занижению их размеров.

'Раздел 7 посвящен исследованию ГЛПС. В подразделе 7.1 методом сращиваемых асимптотических разложений, впервые примененным для ГЛПС А.Н.Осипцовым, получена система уравнений движения газовзвеси в ЛПС на пластине с учетом подъемных сил, действующих на частицы. Проекция уравнения сохранения импульса э-й фракции дисперсной примеси на поперечную ось У для неравновесной по скоростям фаз области ГЛПС имеет следующий вид (здесь Зи= 0, 18= 0)

д 7в д Ув д и и -" + V —" = (V - V )/т + ЪЛ -9 1 (и - и ). (37)

- в а х * а у 9 а " 11 а у 9.

•где Ь4= 3,075 \р т/'2/ 8. В уравнении (37)-последнее слагаемое представляет собой подъемную силу Саффмэна При продольном

обтекании пластины равновесным потоком газовзвеси в ЛПС > поэтому ^ 0, т.е. сила Саф^ыэна направлена к пластине. В результате в области, неравновесной по скоростям фаз, поперечная скорость частиц 7я< 0, что вызывает осаждение примеси на пластине. На рис.16 показаны рассчитанные по монодисперсной модели пб-перечные скорости частиц у пластины (г;°= V ух. / и , х =Х/1 ,

Р* Р Р ОЭ р

1р = П^ гр- длина динамической релаксации частицы) при различных Не6 = и^б/у. Модуль г/Д, достигает максимума и затем стремится к .нулю в -дальней (по х> равновесной то скоростям области ГЛПС. На рис. 17 показаны профили распределенной плотности дисперсной фазы в различных сечениях ГЛПС (С = У[П00/(Хг')]1^г); профиль 1 соответствует области осаждения частиц (неравновесной по скоростям фаз); 2,3 - равновесной,, у стенки происходит накопление примеси; Установлено, что осаждающаяся примесь увеличивает коэффициент . трения . ■ '''•'.:'• ■ '

X |

y

¿ 3,s\

3

ÍW» ~ ~ i o.. A 1 o . J

IVc. 16

1 2 o/o r~c. V

Mr

i}2\

S32 /<?■

e>»

'32

/o'^e.

i 0,1

z o/l

3 0,5

%

1

1 í

j X j

jírK^i/^i-f

! > i/ r" T

1 / A , /Va

// /7' A-

O 3 6 X^-

- / (0,5 ; (33)

Причем upa уззлнчззна Es¿ взла^пна максимума- 0f (х) уменьшается.

2 подразделе 7.2 аналитаческя показано, что езда О&фгмзна аз влияет на структуру ГЛПС з"дальней по х области пластины. 3 частности коэффициент трзния Cf= С°(1 + гг^)1'2. где 0° - коэф|к-цпенэг трзния з однофазном ЛПО. Однако разновесная область расположена по s дальше, чем без учета Р., ,. г Елияет на данами-

о. i s cl i

ку позэдзнгм 'скоростей фаз в области их слабой неравновесноста. 2 атом se подразделе -получен интегральный закон сохранения им-ншудьса в ГШ10

0,5 С. = —5 х -~р -г (п;л„ -1) р Л / (р „U ), (39) f б Z d X ' р* а ' ' pv pv рСО со'} '

гдз 3*. О* - тогдана-потери акпульса несущей и дисперсной фаз. Аналогичное ссотноазниэ, подученное в работах других авторов (И. И.Друж5НЕНская. ;С.Соу), не содержит последнего слагаемого, т.к. в зтях .работах нэ учтены подьвггннв силы, действующе на частицы в ГЛПС. '

Лог-раздел 7.3 носвлцзн- расчету ГЖС с чаотнцегг:, покидагца-поверхность, з двух случаях: 1 ) ЛЕС с пршесьа, которая уносится с поверхности потоком; 2) ЖС с отскоком частиц. В первом ■ случае задача решается с немощью уравнений, подучешшх в подразделе 7.1; в качестве .граничного условия на поверхности задастся .-по кширзчэскоа .формуле поток -уносимой примеси. Во втором случае лаош'летрн частзц рассчитываются по лагрангсевым уравнениям (траз-К'горяый мзтед), а газа - по уравнениям гидродинамшет.

В подразделе 7.4 исследовал полидисперсный ЖС на пластине с испарением, коагуляцией и дроблением капель, образованием сс-кой пленка. Т.к. толщкни ôf и скорость пленки ïïf, образующейся при выпадении.капель на поверхность, невелики, профили Uf, Tf . считались линейными; .0f .определялась из интегрального закона сохранения массы жидкости в плеше. На поверхности пленки задава-■ ззхсь условия.сохранения массы, импульса к тепла. Были проведены : расчеты .npsí'условиях опытов (Хшшда, Ыаеда, Екай, 1980) с малым . влагосодерканкем в набегающем потоке зе^ = 0,01. На рис. 18 показана, шхфшг диаме тра Саутера. По направлению к стенке ô3z умонь-.кается ьз-за испарения капель,- а у стенки значительно возрастает

- £í -

вследствие ш. коггуляшгл. Псе~^,щ:ез обусловлена тем, .450'медкге кеют уке прошли область оегзщенкя (so i) я происходи г-гаргста-.чпе их концентрации у стенка (см, рис. 17).

На рис. 19 представлено сравненаз опжгЗая данных (точки) с рззультатамн расчетов локального коэффициента теплоотдачи пдас-тшш в ГЛЮ а Q/(TV- Т ), -•де Q - лскалыэ! тепловоз поток; а0 - хоэ^геиент теплоотдачи в однофазном ЯСС. Штриховая кривая - ь расчете кгпдп. считались конодгсперс-н^:сг с диамзтром, рагннгг дга-аетру Згутв-ра (в.заякзде£стзпе кшеяь не учитывалось). Результаты расчете?, хорошо согласуются с оютш&тз деннн:,<п, еащ нз ез jm-тытывать коагуляцию в пристенной sens, то вычислэнэое значение с сказывается существенно завысснлсл.

3 разделе-.8 исследован ГШС. Изучению пристенных турбуяэет-пнх потеков газовзвеси посвящено большое количество работ, однако в подавляющем большинстве рассматриваются течения в длинных, узких каналах-(трубах). Анализ немногочисленных опубликованных работ, в которых для изучения ГШС используются уравнения Рей-' нольдса обеих фаз (Возе, Ы.М.Исмаалов), показываю?, что на пути создания модели ГТПС сделаны только" первне шага. Нз иесдадовезв ГШС в соплах, влияние на структуру ГТПС различных ослоезязеех факторов (ускорения внешего течения, полндяспэрскости частиц).

3 подразделе 8.1 метод сраппЕвеках асимптотических разложений, являвшийся обобщением подхода В.В.Сычева на случай запыленного потока, применяется для анализа уравнений: Рейвольдса в ГШС на пластине. При этом считается, что параметр! обеих фаз на внешней границе ТПС определены из ранения систеггн уравнений "невязкого" течения. Показано, что в случае инерционного' осаздения частиц (7sga~Us00) з первом -Приближении их параметры моено считать неизменныш поперек ТПС, равнкш V , U , р ; тогда систзка уравнений изотермического ГТПС имеет еид (<р> =.р00(х))-

9<и > а<и > i a<d> а , а<п >

—а + <у ; - п

9 Я Т 3

<tj > -г + <v > —-2- ---- + - iv -- xu'v^i -

a X Э Y ' 3 Y 5 9 i

ах 3 а т <рд> ах а у 4 а у - 5 Рв«>«ид> - ^«К^р^Ь - а<и;> / + з<?д> / зу = о. (40)

Здесь движение газовзвзеи описывается з рамках модели "замороженного" течения. Частицы, попадая в ТПС с большой поперечной скоростью, "протыкает" его, почти не успевая изменить свои пара-

?.г;трн пг-д влияние:.! газа.- ТЬзсуцал ке среда изняет параметры своз: -. я од действием дисперсной фазы.

иг. роЕЗНИЯ "невязкой" задачи без учета ТПС определено, иг. ¿'-г ;1и"--р;шонного выпадения частиц на поверхность, необходимо учитывать пз»знвЕЕе параметров обеих фаз в ТЕС.

5 подразделе 8.2 проведен расчет ТПС с частицами, осаэдаю-га г;.-:оской поверхности, с помощью уравнений (405,(14) и ур;.г.г-::г.:г; г, г: е. .Последние два уравнения шдифндаровгна в сост-с моделью Лаундера - Шарш для . пристенных течений. Ка с;:с.20 показана профили безразмерной пульсацкоккой энергии несущей сгйд;: г.гл двух значениях йе6 (Ь:+= к/иД, у =- ид/г>, и - ди-насг-юскся скорость), нулевая кривая - однофазный ТПС. Генерация т.. р'улектной энергии в следе за инерцгонншл частпцаш С;к приводит к возрастанию к в пристенной зоне по сравнению с однофазным случае;.;. При уменьшении Не6 возрастает дассапагоя пульсыщонаоЗ энергия на частицах е , поэтому уровень к сшскаетаг.

В подразделе 8.3 получена система уравнений криволинейного ГТЛН с коагуляцией и дроблением частиц,- кекфазнтл,: массообкенок. Результата численного решения система уравнений с соответствующими граннчныш условиями для ГПи сопла Лазали (Б.Н.МасловД.Е. Стернян, А.А.Шрайбер, 1982) .представлены в подразделе 8.4. Ка рис. 21 показано изменение•толшзны ГТПС, рассчитанной по 95-про-ьг-нтному уровню и : 0- агр = 0, 1-4.- гер= 3; 1 - пржлесь полкдк-сперснв; 2 - - келодкепзрена (штриховые крпв'-с); 2 - 105 5/^= ■= 3,0, й#- радиус сопла в критическом сечении при z = г/йх= 0, 1 -„координата, отсчитываемая по оси сопла: 3 - 1,72; 4- - 0,92 (т. е. вторая кривая есть результат расчета по шнолксперской коделк с разгром а<в- орэдвим для полидиспербного ансамбля). На рвевх-рящемся участке сопла (а >0) частицы разных размеров начинает выпадать на .стенку"при разных а; мелкие - благе к кркткч5скг.:гу • сеченюо, чем более краткие. При расчетах с конодисперсной примесью, начиная о точки выпадения частиц соответствующего размера .на стенку, относительная толдкна пограни«шого слоя заметно укень-•:раетсяпо-сравнению с однофазным ТПС (особенно резко при небольших. б - из-за.'большой силы мезкфазного взаимодействия). Полндис-.'•пэрсность • пргмеси приводит к более. плавному уменьшен:» 5.Л1.

(кривая- ' '

>На- тис.22 представлены результаты расчета интегрального :ю-.~эф$щиен?а потерь на трение газа с; стенки с-тда

х

■Сс= -^¡К К С? <ие> / 3 1)У,6Х / ф а® и*). . (41)

о .

Обозначения как на рис. 21. Градиент (5<и.>/бх)у. у стенки увенчивается в двухфазном ТПС до сравнению с'однофазным, поэтому и потери на трениз е ГТПС возрастают. При расчетах с монодисперсными частицами уменьшение б приводит к увеличению £ .

В подразделе 8.5 проведен расчет ГТПС в сопле Лаваля с образованием кидкой пленки из осаждающихся частиц. В расчете считалось, что при 2 = 0 происходит срыв пленки со-стенок сопла. В расширяющейся части сопла пленка появляется после качала осааде-ния частиц наименьшей фракции. Рассчитаны тепловые потоки в стенку: при однофазном течении

= (л. 8 <Тд>/5 , (42)

б ГТПС без шинки

, с<т >-ч

= Iх--9 I —'Е Н<р XV > + <р'7'»С <Т >]„.(43)

* д У ■'V в р - 6

в ГТПС с пленкой'

. о, = а Т, / в -_')у. (44)

На рис. 23 показаны отношения тепловых потоков в двухфазном и однофазном ТПС: '1 - без пленки 0^ = 0р/0°; 2,3 - с пленкой 01о = 2 - по голпдисперсной модели, 3 - монодксперекой

1СР 3. В двухфазном штоке в отсутствие пленки тепловой

поток в.стенку на небольшом участке возрастает более чем в 2 раза по сравнению с однофазным потоком за счет осаздения на стенке частиц'с высокой температурой (второе слагаемое в правой часта (43)).' При налички пленки вблизи входа сопла <3, нарастает по той ке причина, но по мере увеличения бг тепловой поток в стенку . уменьшается, величина ...0^ становится меньше единицы. Пленка за за счет, своего термического сопротивления уменьшает тепловой поток в стенку.

— Г'Ь заключении перечислены основные результаты диссертации

1. Разработаны основы теории .ламинарных и турбулентных течений полидиснерсной газовззеси с сильным взаимодействием частиц ! 2; кеафазным массообменом:

• - 1.1. В рамках ! непрерывного подхода к описанию взакмодейст-вия-частиц (метод Лагранка) получена система уразненкй эволюции Ьсостояния^двухфазной среды. .Предложен гибридный подход, сочетай-

ций элемента гидродинамического к. статистического ошсания движения газовзвеси я позволяющий учесть различие параметров отдельных частиц одной фракции путем введения соответствующих корреля-ционннх моментов.

1.2.'Предложена модель турбулентного переноса в сдвиговых "течениях полидисперсной газовзвеси, включающая

- уравнения для пульсационной энергии газа К и дисперсной примеси отличающиеся учетом меафазкого масообмена, коагуляции и дробления частиц; , -

- представление корреляционных моментов, связанных с турбулентным переносом массы, импульса и тепла обеих фаз;

- учет взаимного влияния турбулентное эффектов и других явлений, имеющих место в потоке.' - .

2. На базе те ори п'.1 получены системы уравнений, проведены аналитические и численные' исследования различных тшов турбулен-.тных и ламинарных течений газовзвеси. В результате' установлено следующее:

2.1. Для турбулентных газодиспзрсшх струй

- применение гипотезы Н.И.Акатнова для (27) позволяет с помощью однопараметрической к-модэли получить согласие с опытными данными не хуке, чем по (к-е)-модели;

- вычисление корреляций, .связанных с дисперсной фазой, в рамках локально-однородного подхода не вносит заметной погрешности в определение осреднениях параметров по сравнению с использованием дифференциальных уравнений для указанных корреляций;-

- расчет испарения капель при отношении массовых расходов фаз аер ~ 1 без учета их коагуляции и дробления приводит к существен. кому закиаэнню среднего размера капель;

2.2. Для газокапельных течений в каналах'. ' .

- рост дисперсии скоростей отдельных фракций обусловлен в основном образованием осколков со скоростями, отличными от скоростей ""старых" капель того не размера;

- в вертикальных колоннах, в отличие от ускоряющихся течений (на- ' пример, в соплах), ведагиин ав кеДкзх фракций сравнимы с вели. чинами средних скоростей ж даш могут превышать их;'' ■ .--2.3. Для газодисперсного*ламинарного-пограничного слоя-- • ^ '.-..',

- при продольном обтекании пластины равновесным-,потоком- .газовз-" -веси подменная сила Сефф^эна, направленная к поверхности,_ при-, :

• вода? к ссэг<двяж>: частиц у. пврадзвй" 'кромки и---к:;расзнревщ) зй- ';;:

гтя та-

равновесной дю скоростям -фаз области. вниз по- сэгоку; -

- обнаруженное-з опнтах существенное уведзчзнгв ксзф&щга нлсютдача нагретой пяасгкна.при ее обтекания потоком с мангн содзрзанЕЭМ каш ль - объясняется жпарзшз-.л л коагудяцизй кэлпз Зракцкй дасшрсной фазы, накаплгсващаися возле пластины.

2.'4. Для турбулзнкюго пограничного слоя

- с вошцьа катода сргщваекшс всзмтотгческпх раздогэвзй показано, 'что з первой прз&шэавя форма безразмерного уравнения двеезеея газа з вязком подслое нз зависят нз только от градиента давления, ко и от сзпн'кеафазного взаимодействия;

- при шардаогЕоа .осгздвния частиц ГШС можно сшспзать з ргкках , лрайзшзэзгя "аакоракешого" течения;

. - £оздэ£ствиэ дисперсной ирзнзси уменьшает инщгцу ТПС и увелн-' чнвеет ковфрщевт'трения. .

3. Получены новые результаты при .изучении физических лро-Шсбов/ярат&ДЕШ с одиночной частицей (кашей);

- получена формула, для реактивной саны, действующей на кспаряз-щуюся дашш в неоднородней'температурном пола;.

- путем сравнения с. 0IHTHÄI3 EäHEKJSi пс газано, что при ударе • ; сфзричаской; -частицы с небольшой скоростью ^ еерагозатуи стенку ,

. * ифавэдужз.-закон 'трения Кулана;

щшлэрщенгаяьно установлено, чго при больпих числах ■ Райноль-~ сагл.Шщуса 0и -укэкызаешея,. прэдгшана

.ешнр^ская-зазискмость ' См(Кеы);

-шц^чааа-уоювяя,-щр ъощщ-чёхящвгтюсж,сдува о даверхнос-.

' ЩЛОКЕНЭТ .ПрЕСХЗНОЧНЫЙ СЛОЙ. ' ' .'. "1 _____ .

.. .V 4/..".;?азр^ташг:Ч2Слэнвыз .штода,-.алгоритмы 'г пакеты пршела-. дннх жротремк'лм расчета: . •"-'/- -/--ГК5 .сгсдатужщэй згдробленизм часищ, их-плавлением,' кшаре-

; -; ГЯГО ,с испарением, коагуляцией и дроблением капель;

- Щ1С с учетов обр^ваная -ездкой плзнки из осаздазддхся капель; ; г-^течаний вгетжлах^с-твэр^щггн'частзцаш.ила кашшш.

: проведены шрокне;численные

; ^оцвсса:в-'ЕЭЕЬто^\уотрса1я®ах .соврекен-

■; А-'з^-техида^зЬя^^ выводов и рзюмсецдацйа.пр -----

'• диавмннх'-даргэдкфс»:'

■ - диааатра .цроюшяетзше фура да-''

'! /котораяг'зэрйзтея ц регулировочном

- 33 - v

клапане, использовать для перемешивания расплава, рассредоточения инертного газа и порошка в объеме ковша;

5.2. В установках для газопламенного нанесения покрытий при уменьшении размеров частиц порошка необходимо сникать расход горючей смеси; в противном случае может произойти перегрев частиц значительно выше температуры плавления, что приводит к разбрызгивание капель, взаимодействующих с поверхностью;

5.". В рамках разработанных моделей находит обьяснгтше.экспериментальный' факт увеличения объемного коэффициента массопере-дачи при очистке газа в скруббере с ростом скорости газа;

5.4. В ТПС сопел Лаваля при наличии жидкой плеши из осаждавшихся капель существенно снижается тепловой поток в стенку.

Основные публикации по теме диссертации:

1. Шрайбер A.A., Гавин Л.В., Наумов В.А., Яценко В.П. Турбулентные течения газовзвеси.- Киев: Наукова думка, 1989.- 240с.

2. Shralber A.A., Ga~In L.B., NaumoV V.A., Yatsenko 7.-P. Turbul jii X flows 1л gas suspensions. - Hemisphere Publishing Corporation. - New York, 1S90. - 248p.

3. Гавин Л.Б., Медведев С.В., Наумов В.А. Модель турбулентной пароквдкостЕой струн и ее численное исследование// Материалы Есес. конференции "Теплофизика и гидрогазодинамика процессов кипения и конденсации". Рига, 1932. - Т.2,Л.2. - С.41-47.

4. Гавин Л.В., Наумов В.А., Шор В.В. Модель (к - а) двухфазной турбулентной струи и ее численное исследование'// й-хзико-химнческне процессы в энергетических установках. - Минск: КТМО анбссг, ;:зз. - c.n-i5. ' ..

5. Гавин JE.Б., Наумов В.А. Турбулентная двухфазная струя и ее численное исследование//"НЕЕ. т 1983.- Т.44, Н б.- С.927-932.

6. Гавин Л.Б., Наумов В.А. корреляция пульсационннх скоростей дисперсной фа_ы в струйных течениях// Кзв. АН СССР. -1933. - N 4. - С.61-65.

Гавин Л.Б., Наумов В.А., Шор В.В. Численное исследование газовой струи с тяжелыми частицам, на основе двухпараметрптеской модели турбулентности// ПМТФ. 1984. - N 1. - С.62-57.

8. Наумов В.А. Турбулентная структура двухфазной струи ■ в условиях динамической неравновесности фаз// Турбулентные двухфазные течения и техника эксперимента.- Таллинн, 1985. - С.83-87.

9. Гавин Л.Б., Наумов В.А. Численное исследование распреде-

лешя лиспесной примеси в турбулентной струе// ИФЖ. - 1985'. - Т. •49, N 4. - C.54S-555. .

10. Газдн Л.Б., Наумов В.А. Влияние дисперсной примеси на турбулентную структуру струи// Доклады АН СССР'.- 1985. - Т. 283, II 2. - С.336-339.

11. Гавин Л.Б..Медведев В.А., Наумов В.А. Модель двухфазной турбулентной струи с горящими в диффузионном режиме дисперсными частицами// Материалы 8 Всес.симпозиума по горения и взрыву. Горение гетерогенных и газовых систем.- Ташкент, 1986. - С.79-82.

12. Гавин Л.Б., Наумов В.А. Турбулентная структура двухфазной струи// Проблемы турбулентных течений. - М.: ' Наука, 1986. -С.135-141.

13. Шрайбер A.A., Гавин Л.Б., Наумов В.А.,'Яценко В.П. Тур-(' буленткне течения газовзвеси// Материалы первой конференции но механике Академий наук социалистических стран. - Прага, 1987. -7.3. - С.249-252.

14. Гавин Л.Б..Медведев В.А., Наумов В.А. Модель двухфазной турбулентной струи с учетом гетерогенного горения частиц// Физика горения и взрыва. - 1988. - Т. 24, N 3. - С.12-17.

15. Наумов В.А. Расчет ламинарного пограничного слоя на пластине с учетом подъемных сил, действующих на дисперсную при-у.зсь// Изв. АН СССР. МЯГ. - 1988. - N6. - С.171-173.

.15. Наумов В.А. Ламинарный пограничный слой с осаждающейся

■ дисперсной примесью// ИФЖ. - 1983. - Т. 55, К 4. - С.555-565.

17. Наумов В.А. Расчет теплообмена и осаждения дисперской приглеси в ламинарном пограничном слое на пластине// 'Моделирование з механике. - Новосибирск, 1989. - Т. 3(20), N 5. - С.71-78.

.13. Наумов В.А. 0 методе аналитического исследования газо-днспэрспого ламинарного двухфазного слоя// Изв. АН СССР. Энергетика и транспорт. - 1989. - IT 4. - С. 170-171.

• 19. Наумов В.А. Два подхода к описанию удара сферического тела о шероховатую стенку// Прикладная механика. - 1989. - Т.25, N 5. - С.116-119.

20. Наумов В.А., Кутузова Т.А. Расчет ламинарного пограничного слоя с инерционной примесью в сужающемся канале// Гидротра-

■ нспорт и судовые системы: Сб.науч. трудов КТИРШС.'- Калининград, 1939.'-- С.25-31.

21. Наумов В.А. Сопротивление продольно обтекаемой пластины в потоке двухфазной жидкости// Гидротранспорт и судовые системы:

- 35 -

Сб. науч. трудов КТИРПХ. - Калининград, 1989. -С.31-36.

22. Наумов З.А. Об уравнениях турбулентного газодисперсного пограничного слоя// Турбулентные течения и техника эксперимента.

- Таллинн, 1939. - С.186-183.

• 23. Наумов В.А. Расчет некзотермического ламинарного"погра-, яичного слоя на пластине с учетом подъемных сил, действувщгх на дисперсную примесь// Теплофизика высоких температур. - 19S0.- Т. 28, К 4. - С.314-815.

24. Каукоз. В.А. Влияние мекфазного массообмена на турбулентную энергию течения газовзвеси// ШТФ. - 1990. - К 1. - С.66-69.

25. Наумов В.А. Об одной ошибке расчета ' движения частицы при сушке// ЙФЖ. - 1990. - Г. 5S, Кб. - С.1С30.

26. Куземко Р.Д., Наумоз В.А. Математическая модель- течения з фурмах для глубинной продувки расплавов порошками- и инертными газами// Вопросы теории и практики сталеплавильного производства. - М: Металлургия, 1991. - С.103-125.

27. Наумов В.А., Кондратьев Л.В., Шор В.В. Расчет газодисперсного турбулентного пограничного слоя в comiáx Лаваля// Моделирование в механике. - Новосибирск, 1991. - Т. 5(22), N 1, С. 153-161.

28. Наумоз В.А. Расчет неизотермического газодисперсного ламинарного пограничного слое на пластине с учетом отскока частиц// Моделирование з механике. - Новосибирск,. 1991. - Т. 5(22), N 4. - С.96-101.

29. Наумов В.А., Кондратьев Л.З. Математическое моделирование тепломассопереноса при полидисперсном газокапельном течении ■ в полом скруббере// Теплометрия и теплосберекение: Сб.науч. трудов Института проблем энергосбережения АН Украины. - Киев, 1991.

— С.69—7о. ■

30. Броунштейн Б.Кг,Кондратьев Л.3.,Наумов З.А.«Поляков О.Л. Экспериментальное и численное исследование полидисперсных газокапельных течений в полнх скрубберах// Аннотации докладов-VII Всесоюзного съезда по теоретической и-прикладной механике.- Москва, 1991. - С.63.

31. Кондратьев Л.В., Наумов В.А., Тихонович Я.З.,'Шор З.В. Экспериментальное и численное исследование .двухфазной турбулентной струи при газопламенном напылении// Теплофизика высоких температур. - 1932. - Т. 30, N 1. - С.139-144. .

32. Броунштейн Б.И., Кондратьев Л.В., Наумов В.А., Поля-

ков О.Л. Тешгомассоперенос при полидисперсном газокапвльном течении в полом скруббере'/ Тепломассообмен-ММФ-92. Тепломассообмен в химико-технологических устройствах. - Т. 11-.- Минск, 1992. - С.7-10.

33.'Кондратьев Л.В., Наумов В.А., Тихонович Я.3., Шор В.В.' Экспериментальное и численное исследование тепломассопереноса, плавлегая и распределения дисперсных частиц е турбулентной струе //Тзплсмассообмен-ММФ-92. Тепломассобмен в дисперсных системах.-Т. 5. - Минск, 1932.. - С.141-144.

34. Наумов" В.А. Динамика тяжелой частицы во вращающейся полости// Теоретические основы химической технологии. - 1992. - Т. £6, II 2. - С.300-304.

35. Наумов В.А. Удар выпуклого твердого тела о шероховатую стопку// Прикладная механика. - 1992. - Т.28; К 4. - С.38-42.

. 36.-Наумов В.А. Расчет■газодисперсного ламинарного, пограни-'кого слоя на пластине с учетом образующейся жидкой пленки//Йзв." Российской АН. МЖГ. - 1992. - N 2. - С. 179-181.

■ 37. Наумов В.А. Ламинарный пограничный слой с . дисперсной примесью, уносимой с поверхности пластины// ИФЖ. - 1992. - Т.63, N 2, - С.189-193. , . "

38. Наумов В.А. Тешгомассоперенос в ламинарном пограничном слое на пластине с коагуляцией,- дроблением, испарением купель и • образованием жидкой пленки// Докл. АН Украины. Математика, естествознание,' технические науки. - 100,? - " 4 - ^ »у_сн

Подписано в печать 15._04.SS... Формат 60x84/16. Букага офсетная 3.' Печать офсетная. Усл. печ. л. 2,0,- /сл. вр.-отт» 2,0. /ч.-изд. л. 2,0__ Тир&к /00 зет. Заказ от го of. Бесплатно.

-^».m. i,. ...............................'..........................ми.

Подготовлено и отпечатан» .п Институт® проблем зкергсс5ерс;сния АК /храины. 254070 Киеп 70, ул. Покровская, 11.