ГАЗОВОЙ ДИНАМИКИ ЧИСЛЕННЫЕ МЕТОДЫ это:

ГАЗОВОЙ ДИНАМИКИ ЧИСЛЕННЫЕ МЕТОДЫ

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


Здесь :


- одномерный случай, - двумерный случай, - трехмерный случай, - плотность газа, р - давление газа, - скорость газа, Т - температура газа, - полная энергия массы газа, - удельная внутренняя энергия, - коэффициент вязкости сжатия, - коэффициент вязкости сдвига, - коэффициент теплопроводности, db k - символ Кронекера.

Следует различать два основных класса задач газовой динамики:

задача Коши для ограниченной или бесконечной области (нестационарные задачи газовой динамики);

стационарные краевые задачи газовой динамики для конечной или бесконечной области.

В свою очередь эти задачи могут подразделяться на ряд классов в зависимости от физич. свойств течения газа. По этому принципу можно различать:

течения идеального газа ;

течения вязкого несжимаемого газа ( - коэффициент сжимаемости, , - скорость звука), описываемые уравнениями Навье - Стокса;

течения вязкого сжимаемого теплопроводящего газа ( ).

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

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

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

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

4) Аппроксимация уравнений не типа Коши - Ковалевской уравнениями типа Коши - Ковалевской при стремлении соответствующего малого параметра к нулю (уравнения Навье - Стокса, уравнения фильтрации).

5) Построение подвижных разностных сеток как регулярного, так и нерегулярного типов.

6) Разделение в разностной схеме задач аппроксимации во внутренних регулярных и граничных точках.

7) Представление сложной системы линейных или нелинейных алгебраич. уравнений в виде рекуррентных соотношений (метод приближенной или точной факторизации).

8) Метод продолжения краевой задачи за границу и включения ее в краевую задачу с простой областью (метод фиктивных областей).

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

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

Методы 1-го класса основаны на представлении обобщенного решения уравнений газовой динамики как совокупности классич. решений, определенных в некоторых областях, покрывающих фазовое пространство и примыкающих друг к другу через общие границы (линии разрывов) с соблюдением условий примыкания (условия динамич. совместности). В каждой области можно применять разностную схему, пригодную для классич. решений, а условия примыкания должны разрешаться с помощью системы, вообще говоря, нелинейных алгебраич. уравнений. Одним из наиболее распространенных методов дискретного представления классич. решений является метод характеристик. Этот метод используется только для решения задач газовой динамики, описываемых гиперболич. уравнениями, и основан он на свойстве гиперболич. системы уравнений иметь, напр., в случае двух неизвестных функций и двух независимых переменных семейство характеристик, к-рые образуют характеристич. сетку, строящуюся в процессе счета. Метод характеристик появился в газовой динамике сравнительно давно и с успехом применялся для расчета одномерных нестационарных течений с небольшим количеством особенностей, а также расчета двумерных стационарных течений в области гиперболичности уравнений. В расчетах используются также и модификации метода характеристик, в к-рых расчет ведется по слоям, ограниченным фиксированными линиями. В случае двух независимых переменных (одномерные нестационарные задачи или двумерные стационарные задачи, сверхзвуковое обтекание) метод характеристик дает возможность избежать интерполяций и тем самым эффектов сглаживания и ап-проксимационной вязкости. Он позволяет точно определять место возникновения ударных волн внутри поля течения как результат пересечения характеристик одного семейства. При большом количестве неизвестных и независимых переменных начинают появляться недостатки этого метода: возникает аппроксимационная вязкость, при наличии большого числа особенностей алгоритм становится логически сложным. Существенным недостатком метода характеристик является также ограничение на шаг сетки, связанное с критерием устойчивости Куранта, и нестрогое выполнение законов сохранения. Поэтому методом характеристик целесообразно рассчитывать задачи, в к-рых число^разрывов невелико. Для метода характеристик доказана сходимость его решения к решению исходной дифференциальной задачи в случае достаточно гладких течений. С развитием ЭВМ, способных решать сложные логические задачи, метод характеристик будет использоваться более эффективно.

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

Основой для построения разностных схем задач газовой динамики является аппроксимация законов сохранения на заданной подвижной или неподвижной сетке, к-рая приводит к сложной системе нелинейных соотношений явного (явные разностные схемы) или неявного (неявные разностные схемы) типа (см. Гиперболического типа уравнение;численные методы решения).

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

схему в виде:

где . Для замыкания соотношений (*) следует связать величины с величинами Это делается различными способами, к-рые приводят или к явным разностным схемам, если выражены через , или к неявным разностным схемам, если выражены через и нелинейным соотношениям, если применяется один целый шаг, и линейным, если применяется два или несколько дробных шагов (метод предикатор-корректор).

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

Указанная общая методика приводит к разностным схемам сквозного счета, как в случае идеального газа ( ), так и в случае диссипативного процесса ( ). Разностные схемы сквозного счета, вводя аппроксимационную вязкость, сглаживают особенности в переходных областях ширины 0(h).и преобразуют ударные волны в ударные переходы, контактные разрывы в контактные полосы. Аппроксимационная вязкость разностных схем объединяет в себе диссипатнвные свойства уравнений газовой динамики и диссипативные свойства самой разностной схемы. Структура аппрок-симационной вязкости определяется дифференциальным приближением схем, к-рое отличается от исходной системы уравнений членами порядка , где - порядок аппроксимации схемы. Во многих случаях представляется важным насколько разностная схема и ее дифференциальное приближение сохраняют групповые свойства исходной системы дифференциальных уравнений. Сохранение разностной схемой групповых свойств имеет большое значение в практич. счете, особенно в задачах газовой динамики, где, напр., неинвариантность первого дифференциального приближения относительно преобразования Галилея приводит к неприятным счетным эффектам (неустойчивость, немонотонность профилей и т. д.).

Известные разностные схемы сквозного счета имеют на гладких решениях локальную точность, как правило, не выше 3-го порядка и глобальную точность не выше 1-го порядка (учитывая невысокую точность разностной схемы вблизи особенностей). Разностные схемы для уравнений газовой динамики должны удовлетворять, кроме независимых требований аппроксимации и устойчивости, еще ряду практически необходимых требований - дивергентности, экономичности, полной консервативности и т. д. Для многомерных задач строить экономичные разностные схемы позволяет идея расщепления. Дивергентность или консервативность разностной схемы означает выполнение в разностных уравнениях разностных аналогов основных законов сохранения (массы, импульса, полной энергии). Свойство полной консервативности требует выполнения разностных аналогов законов сохранения не только массы, импульса и полной энергии, но и различных видов энергии (кинетической, потенциальной, магнитной).

Рассмотрим несколько конкретных разностных схем вида (*). Полагая в формуле (*)


получим явную разностную схему 1-го порядка аппроксимации:


Здесь


Указанная схема является условно аппроксимирующей при (при схема аппроксимирует систему .уравнений ), дивергентной п условно устойчивой. Условие устойчивости имеет вид: , где с - скорость звука. Полагая в формуле (*)


где


получим абсолютно аппроксимирующую явную разностную схему 2-го порядка аппроксимации:


Это схема - дивергентная и устойчивая при условии Разностная схема


является явной при , неявной при при имеет 2-й порядок аппроксимации, при -1-й порядок аппроксимации; при схема полностью консервативная. Разностная схема


имеет 2-й порядок аппроксимации, неявная, абсолютно аппроксимирующая и недивергентная. Здесь


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

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

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

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

Лит.:[1] Рихтмайер Р. Д., Морток К., Разностные методы решения краевых задач, М., пер. с англ., 1972; [2] Годунов С. К., Рябенький В. С., Разностные схемы, М., 1973; [3] Рождественски и Б. Л., Яненко Н. Н., Системы квазилинейных уравнений ..., М., 1968; [4] Самарский А. А., Введение в теорию разностных схем, М., 1971; [5] Жуков А. И., "Тр. Матем. ин-та АН СССР", 1960, т. 58, [6J Xарлоу Ф. X., в сб.: Вычислительные методы в гидродинамике, М., 1967, 316-42; [7] Дородницын А. А., в кн.: Труды 3 Всесоюзного матем. съезда, т. 3, М., 1958, 447-53; [8] Велоцерковский О. М., в сб.: Численные методы решения задач механики сплошной среды, М., 1969, 101-213; [9] Яненко Н. Н., Анучина Н. Н., Петренко В. Е., Шокин Ю. И., в сб.: Численные методы механики сплошной среды, Новосибирск, 1970, т. 1, 40-62; [10] Яненко Н. Н., Метод дробных шагов решения многомерных задач математической физики, Новосибирск, 1967; [11] Самарский А. А., Попов Ю. П., Разностные схемы газовой динамики, М., 1975; [12] Году нов С. К., 3абродин А. В., Прокопов Г. П., "Ж. вычисл. матем. и матем. физ.", 1961, т. 1, М 6, 1020- 1050. Ю. И. Шокин, Н. Н. Яненко.


Математическая энциклопедия. — М.: Советская энциклопедия. . 1977—1985.

Смотреть что такое "ГАЗОВОЙ ДИНАМИКИ ЧИСЛЕННЫЕ МЕТОДЫ" в других словарях:

  • ГАЗОВОЙ ДИНАМИКИ УРАВНЕНИЯ — математическое выражение основных законов сохранения массы, импульса, энергии газа, описывающих состояние газа. Газ есть совокупность большого числа частиц (молекул, атомов, ионов), находящихся в непрерывном хаотич. движении. Учет взаимодействия… …   Математическая энциклопедия

  • ГИПЕРБОЛИЧЕСКОГО ТИПА УРАВНЕНИЕ — численные методы решения методы решения уравнений гииерболпч. типа на основе вычислительных алгоритмов. Различные математич. модели во многих случаях приводят к дифференциальным уравнениям гиперболич. типа. Такие уравнения имеют точные аиалитич.… …   Математическая энциклопедия

  • ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА — раздел математики, включающий круг вопросов, связанных с использованием ЭВМ. Содержание термина В. м. нельзя считать установившимся, так как эта область математики интенсивно развивается в связи с быстро растущими применениями ЭВМ в новых… …   Математическая энциклопедия

  • КРАЕВАЯ ЗАДАЧА — численные методы решения для уравнений с частными производными приближенные методы решения, в результате к рых решение задачи представляется таблицей чисел. Точно решения (в виде явных формул, рядов и т. п.) К. з. можно построить лишь в редких… …   Математическая энциклопедия

  • Яненко, Николай Николаевич — В Википедии есть статьи о других людях с такой фамилией, см. Яненко. Николай Николаевич Яненко Дата рождения: 22 мая 1921(1921 05 22) Место рождения: город Каинск, Томская губерния (ныне Куйбышев Новосибирской области) Дата смерти …   Википедия

  • СССР. Естественные науки —         Математика          Научные исследования в области математики начали проводиться в России с 18 в., когда членами Петербургской АН стали Л. Эйлер, Д. Бернулли и другие западноевропейские учёные. По замыслу Петра I академики иностранцы… …   Большая советская энциклопедия

  • ГАЗОВАЯ ДИНАМИКА — раздел гидроаэромеханики, в к ром изучается движение сжимаемых сплошных сред (газа, плазмы) и их вз ствие с тв. телами. Как часть физики, Г. д. связана с термодинамикой и акустикой. Св во сжимаемости состоит в способности в ва изменять свой… …   Физическая энциклопедия

  • Крыла теория — математическое описание в рамках определенной схемы течения взаимодействия движущегося крыла летательного аппарата с окружающей средой при заданных внешних условиях, геометрии крыла, законах его движения и деформациях поверхности (упругих или… …   Энциклопедия техники

  • крыла теория — Рис. 1. Основные области в схеме несущей поверхности. крыла теория  математическое описание в рамках определенной схемы течения взаимодействия движущегося крыла летательного аппарата с окружающей средой при заданных внешних условиях, геометрии… …   Энциклопедия «Авиация»

  • крыла теория — Рис. 1. Основные области в схеме несущей поверхности. крыла теория  математическое описание в рамках определенной схемы течения взаимодействия движущегося крыла летательного аппарата с окружающей средой при заданных внешних условиях, геометрии… …   Энциклопедия «Авиация»

Книги

Другие книги по запросу «ГАЗОВОЙ ДИНАМИКИ ЧИСЛЕННЫЕ МЕТОДЫ» >>


Поделиться ссылкой на выделенное

Прямая ссылка:
Нажмите правой клавишей мыши и выберите «Копировать ссылку»