ПОИСК Статьи Чертежи Таблицы Разностные схемы с искусственной дисперсией из "Разностные методы решения задач газовой динамики Изд.3 " В атом параграфе описапы некоторые приемы работы, которые оказываются полезными прп численных расчетах задач (причем не обязательно газодинамических). [c.234] ВИЯ задачи коэффициенты вязкости, теплопроводности, газовую постоянную, скорость ударной волны, скорость поршня и т. д. При этой в системе уравнений, приведенных к безразмерному виду, появляется ряд безразмерных постоянных коэффициентов, таких, как число Рейнольдса, число Прандтля и т, п. Обычно каждый класс физических задач характеризуется своим набором подобных безразмерных параметров. [c.235] Если предполагается, что программа для ЭВМ, по которой будут осуш ествляться расчеты, должны быть универсальной, тан что с ее помош ью возможно численное исследование задач различных классов, то может оказаться целесообразным другой подход к обезразмериванию задачи. [c.235] Размерности всех остальных величин выражаются через них. [c.235] Для рассматриваемых задач с плоской симметрией л- имеем размерность г см . [c.236] Результаты решения задачи преобразуются к размерному виду также с помощью соотношений (6.2). [c.236] Прп таком способе обезразмеривания переход к расчету задач иного класса связан лишь с измепением величины масштабов и, Хо, Ро, То и никак не отражается па расчетных формулах. [c.236] Опишем ряд типичных тестов, которые используются при отладке программ, предназначенных для расчетов одномерных нестационарных задач газовой динамики в лаграижевых массовых координатах. [c.238] Здесь прежде всего следует назвать задачу о поршне, вдвигаемом в газ с постоянной скоростью и порождающей в нем ударную волну. Постановка этой задачи и ее решение были описаны в 7 гл. I. Остановимся на особенностях ее численного расчета. [c.238] При у = 5/3 имеем Хк = Н11/Ь. [c.238] Шаг сетки по пространству к можпо выбрать, папример, следующим образом Л/ = 3 4, Л = 30-н40, так чго /г 0,1. Заметим, что если при А = 0,1 шаг по времепи взять равным т = 0,01, то это будет соответствовать в условии (6.9) значению 0,22. [c.239] В силу того, что скорость фронта волны 3) равна единице, время расчета можно ограничить величиной = М/Ф = 3 н- 4. [c.239] В более общем случае, когда необходимо пров( рить функционирование всех частей программы прп расчетах задач газовой динамики с теплопроподпостью. в качестве теста можно использовать точное решение для ударно1т полны в среде с теплопроводностью (см. О гл. I). [c.241] Вернуться к основной статье