|
|
|
|
|
ОТДЕЛ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ТЕХНИЧЕСКИХ СИСТЕМ А.А. Белолипецкий (Гипертекстовый выпуск одноимённой статьи из
книги «50 лет
ВЦ РАН: История, люди, достижения». Оглавление |
|
|
|
1 |
|
|
|
2 |
Математическое моделирование технологических процессов производства
и контроля топливных мишеней |
|
|
3 |
|
|
|
4 |
|
|
|
5 |
Моделирование турбулентного потока сжимаемого газа с
ударными волнами |
|
|
6 |
Моделирование газодинамических эффектов при воздействии
СВЧ-энергии на сверхзвуковое обтекание тел |
|
|
7 |
Аппарат обобщённых характеристик для расчёта разрывных
решений уравнений газовой динамики |
|
|
8 |
|
|
|
9 |
Разработка эффективных топологических алгоритмов
построения изоконтуров |
|
|
|
|
|
|
В Математическое моделирование технологических процессов производства и контроля
топливных мишеней С Управляемый
термоядерный синтез с инерциальным удержанием, или ИТС — инерциальный
термоядерный синтез, основан на очень быстром доведении сверхплотного
ядерного топлива до термоядерных температур. Контейнер, содержащий ядерное
топливо (как правило, это дейтерий-тритиевая смесь), называется мишенью. Исследования
показывают, что абляционное сжатие сферически или цилиндрически симметричной
топливной мишени для формирования центра горения является энергетически
оптимальным процессом. Для этого мишень симметрично облучают со всех сторон
лазерным излучением или иными видами излучения (рентгеновское, ионные и
электронные пучки). При этом энергия излучения поглощается во внешней
абляционной оболочке мишени, аблятор испаряется, ионизируется и разлетается
со скоростью порядка 100 км/с, создавая реактивный импульс давления. Этот
импульс сжимает неиспарившееся топливо, которое в центре мишени достигает
необходимой плотности и температуры, что позволяет зажечь D-T термоядерную реакцию.
Образующиеся в результате её протоны с энергией 3,52 Мэв теряют её во внутреннем
слое холодного топлива и нагревают его. Образуется фронт термоядерного
горения, который быстро распространяется в слое холодного топлива. Теоретические
расчёты показывают, что для достижения условия breakeven (при котором энергия выхода в реакции термоядерного
синтеза не меньше вложенной) необходимо существенное сжатие и нагрев вещества
мишени до температур 10 Кэв. Для этого можно использовать одну из трёх схем
облучения топливной мишени. Это прямое облучение, прямое зажигание, или
быстрый поджиг, и непрямое рентгеновское облучение. Хотя структура мишени и
её физические характеристики зависят от схемы облучения, однако неизменной
остаётся одна черта — мишень должна представлять собой сферическую оболочку с
твёрдым слоем D-T топлива внутри неё. Эта нехитрая, на
первый взгляд, конструкция оказалась достаточно сложной при её технической
реализации. Первая проблема состоит в доставке топлива внутрь многослойной
полистироловой оболочки мишени без нарушения целостности последней. Одно из
её решений заключается в том, что оболочка помещается в контейнер с D-T газообразной смесью,
находящейся под давлением. В результате газ проникает через стенку мишени,
постепенно заполняя её, давление его внутри мишени повышается, вместе с этим
следует повышать и внешнее давление, но так, чтобы не разрушить оболочку.
Поскольку в оптимальном по быстродействию режиме разность внешнего и
внутреннего давлений не должна превышать критического значения, при котором
может разрушиться оболочка, необходимо с высокой степенью точности определять
давление газа внутри мишени, измерить которое невозможно. При этом внутреннее
давление может достигать от 300 до 1000 атм. Было необходимо
разработать адекватную математическую модель заполнения газопроницаемых
оболочек до высоких давлений [1-3],
когда состояние газа описывается уравнением Ван-дер Ваальса. Теоретическим
ядром такой модели является постановка и решение нелинейной краевой задачи
для сингулярно возмущённой системы уравнений параболического типа. Решение
было получено в виде асимптотического ряда по степеням малого параметра.
Результаты модельных расчётов были в дальнейшем использованы при
конструировании системы заполнения оболочек в лаборатории термоядерных
мишеней нейтронно-физического отдела ФИАН им. П.Н. Лебедева. После заполнения газом мишень помещается в
криостат, где охлаждается до температур фазовых переходов, в результате
которых газ оседает на внутренней стенке мишени в виде твёрдого слоя.
Поскольку к однородности, шероховатости и разнотолщинности этого слоя
предъявляются высокие технологические требования, процесс охлаждения должен
проходить так, чтобы эти требования выполнялись. Для этого в отделе был
предложен ряд моделей остывания и кристаллизации газа на стенке мишени [4-6].
Расчёты на них позволили определить конструктивные параметры криостатов,
разработанных в вышеназванной лаборатории ФИАН. Контроль качества, или характеризация, мишени
представляет ещё один важный этап в
технологической цепи её производства. Для этого мишень просвечивают
излучением, коэффициенты преломления и поглощения которого зависят или не
зависят от оптической плотности среды, и по теневым изображениям вычисляют
нужные физические и геометрические параметры оболочки и топлива. В отделе разрабатывались два подхода к решению этой
задачи. Первый базировался на методах прямолинейной томографии. Предполагалось наличие источника рентгеновского излучения
и обработка большого массива теневых картинок. Он не получил дальнейшего
развития из-за дороговизны рентгеновского источника и требуемой системы
сканирования, но главным образом из-за неприемлемого времени обработки
изображений. Второй подход, названный пороговым алгоритмом,
ограничивается анализом небольшого числа изображений мишени при облучении её
длинноволновым источником излучения [7, 8]. Он основан на сравнении
распределения освещённости различных участков теневой картинки с неким
эталонным. Если хотя бы для одного участка освещённость выходит за
теоретически вычисленные пороги, мишень бракуется. На предварительном этапе
необходимо вычислить эти пороги, то есть решить большую серию задач
математического программирования. В реальном времени требуется лишь проводить
сравнение полученного изображения с эталонным, что позволяет ускорить процесс
характеризации по сравнению с томографическим методом на три порядка. Ценою
за это ускорение является потеря информации об объёмной структуре топливного
слоя, однако существенный выигрыш во времени компенсирует этот недостаток. В последнее время задачи контроля качества
топливного слоя включают в себя и знание пространственного распределения
концентрации D-T смеси, степень аморфности её и другие вопросы. Один
из подходов к этому кругу проблем видится в исследовании механизмов возникновения
диссипативных структур в физически однородных средах. Асимптотические методы исследования эволюции диссипативных структур в
уравнениях реакции-диффузии и длинных волн С середины 80-х гг. A.M. Тер-Крикоровым и А.А. Белолипецким изучалось поведение
решений полулинейных уравнений параболического типа в окрестности точек
простых бифуркаций. Интерес к этим задачам возник под влиянием исследований
И. Пригожина и его брюссельской школы по неравновесной термодинамике,
синергетике, к ним же примыкали уравнения реакции-диффузии из химической
кинетики, математической биологии и др. Пространственно однородные решения этих уравнений
теряли устойчивость, и возникали неустойчивые пространственно неоднородные решения,
или диссипативные структуры. Находились условия, при выполнении которых
стационарное и однородное по пространству решение теряет устойчивость.
Строились решения, реализующие переход от такого неустойчивого решения к
некоторому устойчивому, стационарному, но неоднородному по пространству решению. Такие решения перехода были названы
фундаментальными. Оказалось, что фундаментальное решение играет роль
аттрактора. Если начальные возмущения расположены в определённой области, то
решение задачи Коши асимптотически ведёт себя как некоторое фундаментальное
решение. Основной аппарат был разработан сначала для
абстрактного нелинейного параболического уравнения [9], а затем
применён к системе уравнений реакции-диффузии [10-12]. Были
найдены условия, при выполнении которых нахождение подобных решений для
уравнений реакции-диффузии сводится к соответствующей проблеме для
абстрактного параболического уравнения [13].
В Исследование внутренних волн в стратифицированных средах Исследование
внутренних волн в стратифицированных средах представляет
значительный научный и практический интерес, поскольку естественными примерами подобных сред являются мировой
океан и земная атмосфера. В работах [15-18]
в линейной постановке исследовались стационарные задачи обтекания в поле силы
тяжести стратифицированным потоком идеальной несжимаемой жидкости различных
препятствий, находящихся как внутри жидкости, так и на горизонтальном дне.
Рассматривался общий закон стратификации, когда жидкость состоит из слоёв,
внутри которых плотность изменяется непрерывно и терпит разрывы первого рода
на границах слоёв. Исследовалось асимптотическое поведение внутренних волн на
больших расстояниях от источников возмущений. Исследования [18-23]
были посвящены нестационарным проблемам внутренних волн в стратифицированных
средах. За счёт специфического выбора эйлерово-лагранжевых переменных удалось
свести общую гидродинамическую задачу к смешанной задаче для системы двух
нелинейных уравнений с граничными условиями на известных границах. Известно,
что в потенциальном поле в стратифицированной среде сохраняется только
проекция вектора вихря скорости на нормаль к поверхности постоянной плотности
(для несжимаемой жидкости) или постоянной энтропии (для идеального газа). Эта
проекция называется потенциальным вихрем. Была решена задача о внутренних
волнах, возникающих в первоначально покоящейся среде от возникшего в
начальный момент потенциального вихря. В чисто линейной постановке
вертикальных колебаний не возникает, но в предположении, что квадрат
потенциального вихря и вертикальные смещения частиц имеют один порядок,
потенциальный вихрь вызывает вертикальные колебания с возрастающей
амплитудой. Было показано, что наличие распределённых источников на оси вихря
усиливает интенсивность потенциального вихря. Этот факт находится в
соответствии с наблюдаемыми в природе явлениями. Зарождающиеся интенсивные
атмосферные вихри сопровождаются интенсивными осадками в виде дождя, града и
снега. Цикл исследований
[24-31]
связан с построением фундаментального решения для уравнения внутренних волн с
коэффициентами, терпящими разрывы на некоторых горизонтальных плоскостях.
Интерес к подобной проблеме связан со структурой реальной атмосферы.
Известно, что есть сравнительно узкие атмосферные слои, в которых частота
Брента-Вяйсяля, являющаяся важной физической константой, резко возрастает,
оставаясь вне этих слоёв практически постоянной. Но именно эти слои играют
роль волноводов, и поэтому исследование колебаний в этих слоях, возникающих
от точечного источника, представляют практический интерес. Изучались
математические модели двухслойной и трёхслойной атмосферы. Решения были
представлены в виде быстро сходящихся рядов. Показано, что только конечное
число членов ряда даёт существенный вклад в асимптотическое поведение решений
на больших временах. Моделирование
турбулентного потока сжимаемого
газа с ударными волнами Моделирование
турбулентного потока сжимаемого газа с ударными волнами тесно связано с
задачами о воздействии ударной волны на газодинамические флуктуации. В
последние годы увеличивается число публикаций, посвящённых исследованию
взаимодействия ударных волн с турбулентностью (см., например, [32-34]).
Однако в этой области остаётся много нерешённых вопросов в связи с тем, что
сама сжимаемая турбулентность является недостаточно изученной. Моделирование
основано на синтезе разностных схем повышенной точности для уравнений Эйлера
и статистических методов [35-39]. Генерация турбулентных
пульсаций производится с помощью прямого
расчёта стохастических полей параметров газа с различными начальными
условиями, характеризующимися случайными амплитудами параметров. Проводится
усреднение по расчётам с различными начальными пакетами параметров со
случайными амплитудами. Разработана пространственная модель сжимаемой
турбулентности [39], которая является обобщением
подхода, принятого в [35-38], на двумерный случай. Для
газодинамических расчётов используется оригинальный двумерный аналог
разностной схемы на минимальном шаблоне второго порядка точности. Проведены
численные исследования взаимодействия однородной изотропной турбулентности с
первоначально плоской ударной волной. Для различных типов турбулентности
получены новые данные о поведении коэффициентов усиления ударной волной
флуктуаций термодинамических параметров газа, завихрённости, а также
кинетической энергии турбулентных пульсаций. Исследованы коэффициенты корреляции. Проведено численное
исследование
зависимости статистических характеристик различных типов турбулентности от коэффициента корреляции флуктуаций
плотности и продольной скорости газа при взаимодействии с ударной волной. Для
флуктуаций плотности получены корреляционные и спектральные функции.
Исследования проводились для чисел Маха ударной волны M =< 2 и параметра возмущения, лежащего в пределах 0,1-0,3. Результаты
численных расчётов согласуются с экспериментальными данными по коэффициенту
усиления пульсаций плотности и виду корреляционной функции [40].
Исследовалось влияние флуктуаций параметров газа на затухание N-волны [38]. Под N-волной понимается квазиустойчивая волновая структура,
состоящая из двух ударных волн и волны разрежения. Получены новые
количественные данные о затухании интенсивности N-волны после прохождения области турбулентных
пульсаций. Проведены
численные расчёты взаимодействия первоначально плоской ударной волны с
вихревыми структурами для режимов взаимодействия, при которых происходит
разрушение отдельных вихрей [41]. Такая структурная перестройка
поля турбулентных пульсаций может являться причиной регистрируемого в
экспериментах существенного изменения линейного масштаба турбулентности при
переходе через фронт ударной волны. На основе разработанной модели было показано, что
турбулентное возмущение в сверхзвуковом потоке может существенно
воздействовать на аэродинамические свойства тел, а также порождать серию
акустических волн вблизи тела. Таким образом, предложенная модель сжимаемой
турбулентности может быть использована для
изучения воздействия турбулентности на
режимы сверхзвукового обтекания тел. Моделирование газодинамических эффектов при воздействии СВЧ-энергии на
сверхзвуковое обтекание тел Эти работы
ведутся совместно с Институтом высоких температур РАН с Систематическое
исследование воздействия на режим сверхзвукового обтекания затупленного тела
тонкого разреженного канала проводится в работах [42, 43].
Выявлено существенное воздействие, оказываемое модифицированным потоком на
лобовое сопротивление затупленных тел. К этой проблематике тесно примыкают
работы по исследованию взаимодействия
ударной волны с бесконечной высокотемпературной нитью (см., например, [45]).
К настоящему времени количество экспериментальных и теоретических работ в
данной области непрерывно растёт. Моделируются
газодинамические последствия наличия в сверхзвуковом потоке, натекающем на
тело, протяжённого по пространству источника энергии в виде канала с
пониженной плотностью и постоянным давлением, равным давлению в невозмущённом
потоке. Такой канал может быть образован, например, в результате нагрева газа
с помощью инжекции в поток СВЧ-энергии (СВЧ-разряда). Результаты расчётов
согласуются с экспериментальными данными исследования взаимодействия
СВЧ-разряда со сверхзвуковыми моделями по динамике отхода головной ударной
волны и эволюции давления в точке торможения [46-50].
Используются оригинальные консервативные разностные схемы второго порядка
точности, включая различные модификации схемы на минимальном шаблоне [51].
Построение методов основано на идее повышения порядка разностных схем [52].
Схемы дополнены методами выделения ударных волн и контактных разрывов, а
также недиссипативными методиками постановки граничных условий.
Рассматривалось обтекание затупленных и заострённых тел. Исследовалась также
ситуация, когда канал постоянно присутствует в потоке, т.е. имеет бесконечную
длину. Получены новые данные об изменении структуры течения, отхода головной
ударной волны, а также сил сопротивления поверхностей тел в ситуациях
отошедшей и присоединённой головной ударной волны. Показана возможность существенного
уменьшения сил сопротивления поверхностей тел за счёт образования вихря в
области между головной ударной волной и телом. При наличии в потоке
бесконечного тонкого разреженного канала рассматривается установление
стационарных режимов обтекания тел и исследуется картина течения в
стационарных режимах. Полученные эффекты могут найти применение для
изменения характеристик обтекаемого тела в целях совершенствования процесса
управления летательными аппаратами. Аппарат обобщённых характеристик для расчёта разрывных решений уравнений газовой
динамики В секторе,
руководимом доктором физ.-матем. наук В.Г. Грудницким, был проведён анализ
численных методов решения гиперболических уравнений схемами высокого порядка
аппроксимации, полученными алгоритмом Ричардсона и аналогичными приёмами [53].
Было показано, что на больших временах эти подходы не повышают точности
решения, то есть фактически не являются схемами высокого порядка. Возможно
эта работа, а также доклад, сделанный на первой конференции Бабенко в
Бендерах, наряду с другими факторами способствовали тому, что бурный интерес
к этим методам быстро угас. Разработаны и применяются в большом количестве
прикладные схемы второго порядка аппроксимации на минимальном шаблоне [54].
Эти схемы, в отличие от большинства схем данного класса, не требуют
применения сложных процедур сглаживания и значительных перестроений шаблона
вблизи границ. Предложен и интенсивно разрабатывается аппарат обобщённых
характеристик [55]. С его помощью впервые были
получены достаточные условия устойчивости для расчёта разрывных решений
законов сохранения газовой динамики (доказана устойчивость расчёта нелинейных
разрывных решений). Установлена прямая связь между устойчивостью и
монотонностью численного решения. Из большого числа практических задач, решение
которых было получено за время работы в ВЦ, можно упомянуть следующие. Была
предложена принципиальная схема, проведены все необходимые расчёты, курировались
создание и успешные испытания устройства, обеспечивающего вывод мощных
электронных пучков из глубокого вакуума в атмосферные условия [56].
Идеология такого устройства использовалась позднее в прикладных институтах,
занимавшихся подобными проблемами. Отдельные элементы устройства отмечались
премиями и медалями на всесоюзных выставках. Проводились обширные численные
исследования прохождения лазерных и электронных импульсов через атмосферу [57].
Особое внимание было уделено режимам многократной самофокусировки, в которых
удалось обосновать качественное изменение режимов прохождения импульсов и
дать разумное объяснение ряду аномальных эффектов. Были проведены численные
исследования взаимодействия ударных волн с телами различной формы. В
частности, установлен факт затухания колебаний при набегании ударных волн на
плоскую и цилиндрическую выемки. Установлены аппроксимирующие зависимости
величины декремента затухания колебаний от определяющих параметров течения,
хорошо совпадающие с экспериментальными результатами. Был выполнен ряд работ
по взаимодействию ударных волн с
затупленными телами, как покоящимися, так и
движущимися со сверхзвуковой скоростью. Исследовались колебания парогазовой полости в воде. Впервые определён декремент затухания колебаний.
Исследовалось прохождение ударных волн через термики различной интенсивности. Методы
решения обратных задач математической физики, возникающих при восстановлении
эмпирической информации о свойствах, структуре и эволюции астрофизических
объектов Методы
и программы предназначаются для: - устранения
априорно неизвестных атмосферных и аберрационных искажений в
(не) когерентных изображениях (не)стационарных астрономических объектов естественного и
искусственного происхождения; - автоматического
адаптивного обнаружения и распознавания как внегалактических объектов, так и
объектов ближнего космоса; - статистически
оптимального (в рамках адаптивного байесовского подхода) оценивания
пространственных, временных, спектральных и поляризационных характеристик
астрономических объектов на основе алгоритмической обработки
восстанавливаемых изображений. Кроме
того, эти методы приложимы в сфере офтальмологии, технической диагностики, оптической и электронной микроскопии. Актуальность
задач подтверждается их принадлежностью к списку проблем, приведённых
в подпрограмме "Распознавание образов и обработка изображений",
включённой Министерством науки и технологий Российской Федерации в
государственную программу "Перспективные информационные
технологии". В
практическом плане методы позволяют: - наделить
локационную систему большей автономностью (без ущерба её эффективности)
вследствие ненужности дополнительного калибровочного измерения мгновенных или
средних передаточных характеристик канала атмосфера-телескоп по сигналу от
яркой эталонной — естественной или лазерной — звезды); - делают
локационную систему более экономичной, поскольку отпадает необходимость
дооснащать телескоп такими уникальными и дорогостоящими вспомогательными
устройствами, как измеритель волнового фронта сигнала от звезды,
деформируемое (адаптивное) зеркало, лазерный передатчик и электронный блок
управления. Разработка эффективных топологических алгоритмов
построения изоконтуров Усложнение
матфизического эксперимента с учетом возрастающих возможностей вычислительной
техники предъявляет новые нестандартные требования к представлению
результатов. Например, для визуализации картины скалярного поля — функции,
заданной в R3, перспективным представляется подход, связанный с
построением графа Риба или графа контуров, позволяющих трассировать
топологическую эволюцию поверхностей уровня. Стандартные алгоритмы визуализации
изображений и машинной графики плохо приспособлены для построения графа
контуров, поскольку не учитывают или даже уничтожают топологическую
информацию, да и весьма трудоемки. Впервые
быстрый алгоритм вычисления графа контуров кусочно-линейной функции, заданной
на односвязной области в R3 был построен в [58]. В настоящее
время его модификации широко используются в Ливерморской национальной
лаборатории для построения графов контуров в R3 и представления
более сложных геометрических объектов. |
|
|
|
1. Белолипецкий А.А. Нелинейная математическая модель заполнения тонкостенных
оболочек газом // Вестник Московского университета. Серия 15. Вычислительная
математика и кибернетика. 2000. № 2.
С. 7-10. 2. Aleksandrova I.V., Belolipetskiy А.А., Koresheva E.R. et al. An efficient method for filling targets with a highly-pressurized gas
fuel // J. of the Moscow Physical Society. № 9. P. 325-335. 3. Aleksandrova I.V., Belolipetskiy A.A., Koresheva E.R. et al. Mathematical models for filling polymer shells with a real gas fuel. // J. Laser and Particle Beams. V. 17,
№ 4. P. 701-712. 4. Aleksandrova I.V., Belolipetskiy А.А., Koresheva E.R. et al. Free-standing
target technologies for ICF // J. Fusion Technology. V. 37, № 1. P. 166-172. 5. Белолипецкий А.А. Математическая модель остывания тонкостенной оболочки
при быстром конвективном перемешивании газа внутри неё // Вестник Московского
университета. Серия 15. Вычислительная математика и кибернетика. 2001.
№ 1. С. 14-17. 6. Белолипецкий А.А. Математическая модель вымерзания газа на внутренней
стенке лазерной мишени // Серия 15. Вычислительная математика и
кибернетика. 2002. С. 23-28. 7. 8. Александрова И.В.,
Белолипецкий А.А., Писарницкая Е.А. Пороговый метод характеризации лазерных мишеней в
проблеме лазерного термоядерного синтеза // М.: ВЦ РАН, 2003. С. 1-75. 9. Белолипецкий А.А.,
Тер-Крикоров A.M. Построение фундаментальных решений абстрактного нелинейного
параболического уравнения в окрестности точки бифуркации // Матем. сб. 1985.
Т. 128. Вып. 3. С. 306-320. 10.
Белолипецкий А.А.,
Тер-Крикоров A.M. Об одном классе решений
абстрактного нелинейного параболического уравнения вблизи точки бифуркации
// Докл. АН СССР. 1984. Т. 279, № 4. С. 777-780. 11.
Белолипецкий А.А.,
Тер-Крикоров A.M. О фундаментальных решениях
нелинейного уравнения теплопроводности // Ж . вычисл. матем. и матем.
физики. 1984. Т. 24, № 6. С. 850-863. 12.
Белолипецкий А.А.,
Стронгина Н.Р., Тер-Крикоров A.M. Некоторые вопросы эволюции диссипативных структур с точки зрения
теории бифуркаций // Математическое моделирование. Методы описания и
исследования сложных систем. / Научн. сб. под редакцией А.А. Петрова и А.А.
Самарского. М.: Наука, 1989. С. 25-54. 13.
Белолипецкий А.А., Тер-Крикоров A.M. Бифуркации в уравнениях реакции-диффузии и связанные с ними
диссипативные структуры // Сборник трудов факультета ВМК МГУ им. Ломоносова.
Численные методы и вычислительный эксперимент. Диалог. М.: Изд-во МГУ, 1998.
С. 15-30. 14.
Белолипецкий А.А.,
Тер-Крикоров A.M. Асимптотические решения типа длинных
волн для одного класса граничных задач математической физики // Нелинейный
анализ и нелинейные дифференциальные уравнения. / Научн. сб. под редакцией В.А. Треногина и А.Ф. Филиппова.
М.: Физматлит, 2003. С. 145-196. 15.
Тер-Крикоров A.M., Бежанов К.А.,
Онуфриев А.Т. Пространственная задача
обтекания неровности дна потоком стратифицированной жидкости конечной
глубины // Известия АН СССР. 1990. № 3. С. 101-111. 16.
Тер-Крикоров A.M., Бежанов К.А. Пространственная задача обтекания источника потоком многослойной
жидкости // Известия РАН. Физика атмосферы и океана. 1993. Т. 29.
№ 6. С. 771-779. 17.
Тер-Крикоров A.M. Исследование задачи обтекания тела плоским стратифицированным
потоком // ПММ. 1993. Т. 57. Вып. 3. С. 41-49. 18.
Тер-Крикоров A.M. О силах, действующих на
препятствие в стратифицированном потоке // ПММ. 1993. Т. 57. выпуск
1. С. 58-64. 19.
Тер-Крикоров A.M. Эволюция вихрей и внутренних волн в стратифицированной
жидкости // ДАН СССР. 1994. Т. 338. № 1. С. 54-56. 20.
Тер-Крикоров A.M. Вихри и внутренние волны в стратифицированной жидкости // ПММ.
1995. Т. 59. Вып. 4. С. 599-606. 21.
Тер-Крикоров A.M., Гордейчик Б.Н. О равномерных аппроксимациях фундаментального решения
уравнения внутренних волн // ПММ. 1996. Т. 60. Вып. 3. С. 443-450. |