Математические модели для эффективного управления некоторыми теплофизическими процессами

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

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

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

ДМ5-ТБХ & /СрасТУ-ТЕХ 4 сов. В работах [1] - [3] исследовались вопросы устойчивости стационарных решений смешанных задач для квазилинейных уравнений параболического типа. В [4] был получен результат о стремлении ограниченного решения к стационарному для автономного квазилинейного параболического уравнения второго порядка. В [5] - [10] эти работы были развиты для неавтономных уравнений. В связи с этим автором были получены априорные оценки первой производной при условиях, отличных от условий, при которых получены оценки в [14] - [16].

В [43] - [45] решается задача оптимального управления для уравнения теплопроводности в классической постановке. В [43] доказано существование и единственность решения задачи оптимального управления, в [12] выведены необходимые и достаточные условия оптимальности, в [44] приводится численное решение этой задачи. При этом существенным является то, что решение разлагается в ряд по собственным функциям. Автором в [46] получены необходимые условия оптимальности для почти — линейного уравнения, если /и < к. Если уравнение линейное, то получены необходимые и достаточные условия оптимальности. Позднее задачи управления неадиабатическим трубчатым реактором, используемым в химической технологии, исследовались в работах [47], [49]. Математическая модель реактора задавалась системой интегро — дифференциальных уравнений.

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

Для многих приложений [50] очень важно иметь компактный отпаянный лазер, который не был бы связан с водяным источником охлаждения. Обычно при конструировании радиоэлектронных приборов, требующих для своей нормальной работы развитой внешней поверхности, используется ореб-рение. Автором в работах [51], [53] впервые рассмотрена принципиальная возможность воздушного охлаждения капилляра за счет оребрения внешней поверхности. Вторая работа в России в этом направлении [59] появилась через 5 лет. Автором был рассмотрен стационарный процесс теплопередачи в цилиндрической стенке лазера, который может быть записан в виде системы уравнений теплопроводности. При этом нужно решить следующие вопросы: исследовать принципиальную возможность охлаждения газовой смеси до фиксированной температуры за счет оребренияв тех случаях, когда это возможно, рассчитать площадь поверхности оребрения.

Здесь, кроме теоретических, встречаются чисто практические трудности. Если коэффициент теплопроводности табличная величина, то коэффициент теплоотдачи не является табличной величиной и подлежит либо оценке, либо экспериментальному определению. При расчете мы использовали феноменологический подход. Требующиеся для расчетов параметры либо оценивали на основании физических результатов, либо определяли экспериментально. Расчет системы охлаждения был доведен до построения реального лазера.

В лазерной оптике одной из важнейших задач является усовершенствование методов эффективного преобразования частоты лазерного излучения. а2.

Коэффициент эффективности преобразования вводится как т/ = где а8.

30 амплитуда генерируемого поля.

Процесс сложения частот в условиях двухфотонного резонанса широко используется в ультрафиолетовом (УФ), вакуумно — ультрафиолетовом (ВУФ) и инфракрасном (ИК) диапазонах спектра. При апконверсии частот слабого ИК излучения достигнут коэффициент преобразования (КП) ~ 40 — 65 процентов [80], [81], а при генерации ВУФ излучения [82] - [86] максимально достигнутый коэффициент преобразования составляет ~ 1 процент. Однако вопрос о предельно достижимом коэфициенте преобразования в процессах нелинейного смешения все еще остается открытым. Вопрос о предельных значениях КП в процессах смешения частот типа т3 = 2уо ± IV 2 приведен в работе [86]. Автором анализируется вопрос о предельно достижимом КП при сложении частот самого общего вида ии3 = и) х -{- и)2 и>з при условии, когда сумма частот и>1 + и>2 = юпд, где — частота двухфотонного перехода. Анализ делается на основании математической модели. Заметим, что вследствие того, что сигналы с трудом отличимы от шумов, эксперимент был выполнен после того, как на основании математической модели было указано соотношение на начальные поля и был проведен численный эксперимент.

Изучение закрученных потоков жидкости представляет значительный теоретический и практический интерес. Течение жидкости в вихре имеет весьма сложный характер. При расчетах полей скоростей и давления используется большое количество эмпирических и полуэмпирических зависимостей. Одной из первых работ, посвященных попыткам численного решения уравнений Навье — Стокса для закрученных потоков, является работа [63]. В расчетах вихревого движения вязкой жидкости предложен ряд разностных схем.

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

В данной работе обсуждается возможность описания стационарного осе-симметричного цилиндрического вихря в приосевой зоне на основе системы уравнений Навье — Стокса.

Показана возможность количественного описания приосевой зоны на основе выбранной модели, определяются параметры модели [64], [65], [66], [67].

Результаты экспериментов [69], [70], [72] показывают, что в вихревых камерах цилиндрического вида возможно образование зоны, в которой движение однородной вязкой несжимаемой жидкости является осесимметричес-ким, причем радиальная и тангенциальная компонента вектора скорости зависят лишь от расстояния до оси цилиндра, осевая компонента растет линейно по г вдоль оси цилиндра.

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

В наших работах [64], [65], [66] предложен метод сравнения результатов расчета с экспериментом. В них предлагается вести сравнение результатов расчета для полей скорости в вихрях с сильной закруткой, в которых V2 + ги2 «и2, т. е. для случая, когда величина полной скорости потока велика по сравнению с величиной радиальной компоненты скорости. Таким образом, параметры в задаче выбираются такими, чтобы квадратичное отклонение рассчитанной функции V2 = ао + Ь^(Т')2 от экспериментальной было наименьшим. После того, как параметры таким образом выбраны, можно рассчитать радиальную компоненту скорости, которая в эксперименте не измерялась.

На основании системы уравнений магнитной гидродинамики (МГД — уравнений) автором проведено исследование влияния электромагнитных полей на движение металла и электролита в алюминиевых электролизерах Красноярского алюминиевого завода (КрАЗа). В мировой литературе работы, посвященные математическому моделированию процессов в электролизере, появились сравнительно недавно [90], хотя экспериментальное измерение скорости движения расплава с помощью радиоактивных изотопов было проведено еще в 1973 г. [92], а экспериментальное исследование влияния электромагнитных полей на поверхность расплавленного металла начали проводить в России в 1959 г. [94].

Для различных типов электролизеров математические модели будут разными. Автором исследуются вопросы, связанные с процессом Эру — Холла производства алюминия на необожженных анодах. Следует отметить, что работы отечественных авторов являются наиболее обоснованными и доказательными.

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

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

Процесс Эру — Холла производства алюминия основан на электрохимической реакции превращения глинозема в алюминий.

Электролизер состоит из угольного анада, погруженного в расплавленный криолит. Криолит является электролитом, в котором растворен глинозем. Ванна заполнена расплавленным алюминием, являющимся катодом. Анод взвешен в расплавленном электролите. В процессе работы анодный материал расходуется, поэтому анод должен непрерывно опускаться для сохранения фиксированного расстояния анод — катод.

В электролизере Содерберга анод обжигается на месте, для этого используется тепло от пиролиза угольной пасты. Электрические контакты сделаны посредством вертикальных штырей.

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

На работу электролизера оказывает влияние, например, соседний ряд электролизеров. Таким образом, нужно создать такую модель, которая реагировала бы на изменение конкретных параметров, поддающихся управлению. Па Красноярском алюминиевом заводе (КрАЗе) в настоящее время работают мощные электролизеры, рассчитанные на 165 кА. Как отмечено в [99], влияние на гидродинамические процессы в жидком электролите и металле, а также на физико — химические процессы в них электромагнитных сил, пропорциональных квадрату тока, очень велико. Эффективность работы электролизера оценивается выходом металла по току (77). В [117] анализируется механизм выхода по току в современных ваннах с теоретической и практической точек зрения. Обсуждаются различные факторы, влияющие на 77, такие как циркуляция электролита под действием электромагнитных сил, диффузия металла, температурный баланс и т. д.

В литературе имеется ряд эмпирических и теоретических формул, характеризующих выход по току. Согласно эмпирической формуле Коробова: г] = 1 — 256 700 • 5а0'21 / (/а0'58. Ь • е1294°/т), где.

5а — площадь анода, м2;

За — анодная плотность тока, А / см2;

Ь — междуполюсное расстояние, см;

Т — температура электролита.

Изменение выхода металла по току в зависимости от формы поверхности металла определяется следующим образом: где г]н — номинальный выход по току, — коэффициент, учитывающий влияние перекоса металла, и) — отношение фактического и номинального тока серии. Зависимость выхода по току от скорости дается формулой: где к — коэффициент, учитывающий физико — химические свойства электролита и особенности технологии;

V — скорость движения межфазной жидкости.

По одной из гипотез, снижение выхода по току пропорционально г — ротору в электролите: где /с — ток серии, А;

5 — площадь электролита.

На выход по току металла существенное влияние оказывает форма раздела поверхности металл — электролит и циркуляция металла и электролита [92], [93], а в результате увеличения мощности электролизера эти эффекты стали играть очень важную роль. С одной стороны, необходимо обеспечивать.

Л?7 = (1 ~ Цн,) • / / и) — г] = 1 — к — V.

0,83 3 симметричность и многоконтурность расплава, избегать застойных зон. С другой стороны, с увеличением скорости расплава снижается выход по току, быстрее размываются стенки электролизера.

Необходимость оптимизации магнитных полей от различных токонесущих частей электролизера была осознана уже давно [89], но изучение магнито-гидродинамических (МГД) явлений с целью управления процессом началось сравнительно недавно [90], [91].

В связи со значительной затрудненностью непосредственных измерений на ванне (высокая температура 963° С, химическая агрессивность расплавов, сильные магнитные поля), а также высокой стоимостью различных модификаций конструкции широкое распространение получили методы математического моделирования и вычислительного эксперимента.

Приведем обзор работ, представляющих практический интерес.

Работа [118] содержит подробную математическую модель для описания стационарных магнитогидродинамических состояний электролизера при следующих физических предположениях: жидкости несмешиваемы, несжимаемытечение описывается уравнениями Навье — Стоксаэлектромагнитные поля определяются уравнениями Максвелла, законом Омавыполняется закон сохранения массыдавление считается константой на свободных поверхностяхвязкость рассматривается как параметр.

Выводятся уравнения для нахождения распределения плотности электрического тока, магнитных полей, нахождения скорости движения жидкости, распределение давления. При этом используется предположение о бесконечной длине электролизера. Авторы работы [118] считают, что для реальных электролизеров это предположение ослабляется учетом расхода жидкости. В этих предположениях обсуждаются результаты вычислений для реальных электролизеров. Представлены примеры распределения потенциала, скорости, магнитных полей, определенных в поперечном разрезе. Работы зарубежных авторов, в основном, относятся к электролизерам с обожженными анодами, а не к электролизерам Содерберга.

Гидродинамическое моделирование обслуживания коммерческого электролизера описано в работе [119]. Показаны различные образцы течения и поверхности для разных конструкций электролизера и изменений в управлении. Форма поверхности раздела зависит от дивергенции Лоренцевых сил. Волновое движение может быть дано взаимодействием между током и изменениями магнитных полей. Индуцированный ток частично подавляет этот эффект. Отмечены интересные эффекты: сильное влияние настыли на горизонтальный ток, соответствие ряби на поверхности раздела петлям течения, некоторое соответствие течения электролита (без учета газовыделения) и металла. Рассматриваются расчеты с учетом и без учета настыли. В случае отсутствия настыли скорости течения становятся значительно больше и форма поверхности раздела значительно деформируется по краям. Рассматривается случай попарного отключения анодов. Под отключенными анодами наблюдается возмущение течения (на рисунке незначительное) и «ямки» под отключенными анодами. Рассчитывается развитие событий с интервалом по времени после отключения анодов. Отмечено увеличение вогнутости поверхности сразу после отключения, значительные деформации в углах ванны, стабилизация с течением времени.

Работа [115] подчеркивает значительные экономические успехи в связи с улучшением МГД — ситуации в алюминиевом электролизере и при отливке жидкого металла. В течение последних нескольких лет (с 1973 г.) центр АЪиБХЛЗБЕ — ЬО^А занимался улучшением управления электролизера. За это время потребление энергии было снижено с 15,1 до 14,2 к^Ь/к^А!, выход по току повысился с 88% до 93% и время жизни ванны увеличилось до 2000 дней. Отмечается, что электромагнитные силы могут играть пассивную и активную роль в МГД — процессах. Описываются гипотезы и вычислительная методология, положенные в основу модели: вычисления ведутся в трехмерном пространствеучитывается эффект стальных частейполе скорости вычисляется по уравнениям Навье — Стокса с постоянной эффективной вязкостью в стационарном случаепренебрегается газами и при перерасчете индуктивной составляющей не учитывается форма поверхности раздела.

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

Эффективная вязкость определяется с помощью экспериментальных данных. Модель постоянно согласуется с экспериментами. Рассматривается поле скорости в случае различных токоподводов и также возмущение потока в случае анодной замены. Последний эффект более нагляден, чем в предыдущей работе, наблюдается сильное возмущение под меняемым анодом, это может привести к нарушению технологического режима. Неоднократно подчеркивается необходимость управления течением для получения хороших производственных результатов.

В работе [120] уделяется значительное внимание изменениям течения и скорости в случае замены угловых анадов, а также в случае нарушения связи между различными анодами. Рассчитанные поля скорости сравниваются с замеренными скоростями. Наблюдается некоторое соответствие.

В [93], [90] отмечается, что движение расплавов определяется действием сил, вызванных взаимодействием между электрическим током и магнитным полем, выделением газа на аноде, градиентами температур в электролизере. Одним из первых обратил внимание на циркуляцию электролита и алюминия на электролизерах большой мощности С. М. Мещеряков [121]. Он предположил зависимость циркуляции от сил, вызванных взаимодействием электрического тока и магнитного поля, дал приближенные схемы движения расплавов.

Измерения на работающих электролизерах затрудены из-за высокой температуры 960° С) и химической агрессивности криолитового расплава. Попытка измерить циркуляцию электролита и алюминия на промышленных электролизерах методом радиоактивных изотопов была впервые предпринята М. Ф. Даграмаджи и В. Н. Рудаковым [94]. Ими также были получены приближенные формулы для нахождения колебаний границы раздела металлэлектролит. В [95] также методом радиоактивных изотопов получены схемы движения расплавов, подобные схемам в работах [94], [96]. В автореферате Н. С. Сираева приведены результаты измерения скорости течения электролита и металла, контуры течения. Скорость течения поверхностных слоев электролита в пространстве борт — анод 18 — 30 см/с. Повышенная скорость — правый угол выходного и левый угол входного торцов. Скорость электролита под анодом значительно выше, чем над металлом, из-за газов. Тепловая конвекция играет крайне малую роль — до 1,5 см/с. Циркуляция электролита в междуполюсном зазоре (МПЗ) и металла очень похожи. Толщина газоэлектролитного слоя в МПЗ до 2,2 см. Также имеются работы по измерению скорости металла методом растворения железных стержней [97].

Исследования распределения скоростей расплава алюминиевого электролизера, проводимые на промышленных агрегатах, приведены в работах [95], [97], [122], [123]. Прямые измерения скорости на магнитогидродинамичес-ких моделях [124], [125] дают хорошие результаты, но требуют уникального аппаратурного оформления. В [126], [127] предпринята попытка изучения характера движения на физической модели электролизера с помощью построения распределения электромагнитной силы по данным измерений плотности тока и индукции магнитного поля. Там показана хорошая корреляция между визуально наблюдаемым характером течения и полем электромагнитной силы. Такой подход возможен, если можно наблюдать движение расплава. Если же визуальное наблюдение недоступно, то по полю электромагнитной силы весьма трудно предсказать характер течения. Кроме того, такие исследования малодоступны в условиях промышленного эксперимента из-за невозможности аппаратурного оформления измерения горизонтальных токов в расплаве. В [100] проведено численное исследование течений расплавов в алюминиевом электролизере в поперечном разрезе. В работе [112] рассматривается допустимость того или иного приближения в зависимости от мощности электролизера. В работе [105] была предложена, а в [106], [107] получила свое развитие модель Моро — Эванса, основанная на оценке вклада отдельных членов. В этих работах плотность тока считается линейной функцией координат. Нами в работах [108] - [111] предложена модель для расчета плотности электромагнитных сил, усовершенствована модель [106], [107] для рассматриваемого типа электролизеров, предложены другие краевые условия, компоненты плотности тока не являются линейными функциями координат. Предложена оценка коэффициента турбулентного трения в металле и электролите. Проведены исследования влияния магнитного поля на форму поверхности раздела металл — электролит и на линии течения в металле и электролите.

Диссертация состоит из введения, четырех глав, заключения и библиографии.

Заключение

.

В представленной диссертации рассматривается ряд задач, связанных с математическим моделированием теплофизических процессов, их управлением или конструированием. Задачи оптимизации процесса (физического или технологического) не всегда носят классический характер задач вариационного исчисления или задач оптимального управления.

В первой главе рассматривается задача управления тепловым процессом, описываемым почти — линейным параболическим уравнением, где состояние управляемого объекта характеризуется функцией и (х, которая является решением следующей задачи: где ¡-и < к.

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

Доказано существование такого управления, получено необходимое условие оптимальности для почти — линейного уравнения и получено необходимое и достаточное условие оптимальности в случае, когда /(и) = /(ж, ?) • и.

Лм^-ТЕХ & КрасТУ-ТШ 186.

Щ = ихх + ¡-(и), ?>0, 0 < ж < 1- и (х, 0) = 0, (0<ж<1), 1^(0,*) =0, t>0- иж (1, ?) + аи (1, ?) = ар (Ь), а = const > 0, t>0 1 т о о.

Автор получает результаты, касающиеся оценки их (х, ?)|с (од)> существования решения и устойчивости стационарного решения для квазилинейного параболического уравнения, которые приведены им в § 1 — § 3 первой главы.

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

В § 1 этой главы рассматриваются принципы охлаждения газоразрядных лазеров, анализируются основные характеристики теплосъема с поверхности. Далее проанализирована принципиальная возможность беспринудительного охлаждения за счет оребрения внешней поверхности. В данном случае задача управления ставится так: при заданной вводимой в разряд мощности требуется рассчитать площадь оребрения так, чтобы внутри газоразрядной трубки температура не превышала заданного значения. Предложен метод расчета, рассчитана система воздушного охлаждения для волноводного СС>2 — лазера. На основании этого расчета была создана действующая модель компактного лазера. Приведен расчет системы воздушного охлаждения для N2 — лазера. Полученные формулы позволяют оценить, при какой мощности <5 беспринудительное охлаждение возможно, а при какой — нет. В этой же главе анализируется математическая модель, описывающая че-тырехволновой процесс. В случае фазового согласования, и при наличии так называемого «захвата фаз», систему уравнений можно привести к виду: = —аз кп[ка$а3- а±а2], daz da df da2 as к n [к аз as — ai a2], = a2 к n [к аз as — ai a2], = ai к n [к аз as — ai a2], ai (0) = ai0, а2(0) = а2о, аз (0) = а30, a5(0) = 0.

Требуется найти такие аю, а2о, азо, чтобы коэффициент преобразования г] = принимал максимальное значение. Автором найдены первые интег-азо ралы этой системы, найдена связь между коэффициентом преобразования г), начальными плотностями потоков излучения а2о / аюазо / аю и параметром к. Полученное соотношение позволяет сделать некоторые выводы о коэффициенте преобразования в конкретных ситуациях. В частности при к = 1 при а2о / аю = 1 и а30 / аю = 1, Vmax = 0,2.

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

В случае, когда в вихревой камере цилиндрического вида возможно образование зоны, в которой движение однородной вязкой несжимаемой жидкости является осесимметрическим, причем радиальная и тангенциальная компонента вектора скорости зависит лишь от расстояния до оси цилиндра, а осевая компонента растет линейно по z вдоль оси цилиндра, система уравнений Навье — Стокса, записанная в цилиндрической системе координат, расщепляется, и для Т = — —, где и — радиальная компонента скорости, получаем следующую задачу:

Т" > = i- [(Г')2 — Т’Т + 2) + С], zx.

Т{0) = О, Г (х0) = То, ^(хо) = Ti, (ж =) .

Предложен метод решения краевой задачи. Так как при ж — 0 уравнение имеет особенность, решение ищется в виде ряда оо.

Т = ai х%, ?=i получены рекуррентные формулы для отыскания коэффициентов ряда. Предложена методика сравнения с экспериментом, впервые показана возможность описания эксперимента с помощью системы уравнений Навье — Стокса. Это позволяет рассчитать величины, которые в эксперименте не измерялись.

В четвертой главе исследуется влияние электромагнитных сил на линии течения металла, электролита, скорости течения, форму поверхности раздела металл — электролит. Главным показателем работы алюминиевого электролизера является выход металла по току. Имеется ряд эмпирических формул, связывающих показатель выхода по току с различными характеристиками электролизной ванны, в частности, с формой поверхности раздела металл — электролит, со скоростью, с температурой и т. д.

В основе математической модели лежит полная система уравнений магнитной гидродинамики. В работах [108], [109] предложена математическая модель для расчета (плотности тока) и ~Й (магнитной индукции), основанная на широком использовании натурных измерений, реально учитывающих зависимость и if от скорости. Таким образом, автор рассматривает стационарную систему уравнений Навье — Стокса сIш у*" = 0, где = х ~Й — эффективное распределение плотности электромагнитных сил.

На основании оценки вклада отдельных членов в систему уравнений, показана применимость модели Моро — Эванса к электролизерам рассматриваемого типа.

Получена формула расчета формы поверхности раздела.

Щх’У)=(р2-Рг)9 + ЛХ Р2 > ~ < Р1 > +.

2 + Н г #1 +Р2 д —2— + ^ «У учитьгоающая.

Предложена методика сравнения с экспериментом. Предложена оценка «коэффициента турбулентного трения» в металле и электролите на основании геометрии ванны и распределения электромагнитных сил. Оценка в металле согласуется с экспериментальными данными, получена оценка коэффициента турбулентной вязкости в металле.

На основании проведенных расчетов проанализировано влияние Вх, Ву, Вг (компонент магнитного поля) на линии течения и форму поверхности раздела металл — электролит. Проведен сравнительный анализ опшповок, применяемых на заводе.

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

Получены следушщие’результаты:

1. Рассмотрена задача оптимального управления тепловым процессом, описываемым почти — линейным уравнением параболического типа. Получено необходимое условие оптимальности. В случае, если уравнение линейное, то получено необходимое и достаточное условие.

2. Получена оценка первой производной и теорема об устойчивости стационарного решения первой краевой задачи для квазилинейного параболического уравнения при выполнении условия Л.

Границы области устойчивости позволяют обеспечить устойчивость процесса за счет выбора начального режима.

3. Впервые предложен расчет конструкции компактного СО2 — лазера с воздушным охлаждением, на основании которых создана промышленная модель компактного лазера.

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

5. Анализ математической модели, описывающей четырехволновой процесс, позволяет найти связь между коэффициентом преобразования, начальными плотностями потоков излучения и коэффициентами восприимчивости. Выделена область начальных данных, при которых коэффициент преобразования равен 0,2.

6. На основании исследования приосевой зоны осесимметричного вихря предложен метод описания эксперимента и на его основе — расчет тех величин, которые в эксперименте не измеряются.

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

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