ПОИСК Статьи Чертежи Таблицы Замечания о расчете дозвуковых из "Вычислительная гидродинамика " МОЖНО рассчитать по этой программе, считая число Мо малым, скажем положив Мо = 0.1. Поэтому, казалось бы, программа расчета течений сжимаемой жидкости обладает большей гибкостью. [c.422] В обшем случае, однако, подобный метод будет и весьма неэффективным, и весьма неточным. [c.422] При Мо 0.1 мы имеем а 10 м , и по условию (5.183) At ограничено в основном скоростью звука, что уменьшает максимально допустимые значения А/ в десять раз по сравнению со случаем применения уравнений для несжимаемой жидкости. Более того, паразитные звуковые волны на сетке приведут к возрастанию ошибки, связанной с неразличимостью, и (что, быть может, наиболее важно) ухудшат итерационную сходимость. Черни с соавторами [1950] указал на желательность эффективного отфильтровывания этих паразитных волн и применения поэтому уравнений течения несжимаемой жидкости. [c.422] Ухудшение точности проистекает из неопределенности прп М- 0. Дал е при поверхностном знакомстве с газовой динамикой ясно, что для сверхзвуковых течений используются отношения давлений, например отношение давлений в донной области Рг = Рв/Роо, В то время как для течений несжимаемой жидкости используются разности давлений ), например коэффициент давления в данной области Срв = (Рв — Роо)/дх, где доо — скоростной напор в невозмушенном потоке, т. е. [c.422] При 7 = 1.4 и Моо = 0.1 получаем отношение относительных ошибок г = 0.007. Это означает, что ири вычислении Рг по уравнениям течения сжимаемой жидкости с погрешностью менее 1 °/о получаем значение Срв с погрешностью 100%. [c.423] Для смешанных течений, т. е. для сверхзвуковых течений с областями дозвукового течения, отношения давлений еще имеют смысл и можно пользоваться уравнениями сжимаемой жидкости. Однако задачи с большим интервалом изменения числа Маха по времени очень трудны для расчета. Примерами таких задач являются разгон тела из состояния покоя до сверхзвуковой скорости и расчет взрыва от его начала до поздних стадий. [c.423] Харлоу и Амсден [1968] разработали неявный эйлеров метод расчета движений сплошной среды (метод I E), дающий хорошие результаты от М = О до М 1. Метод I E основан на расчете уравнения неразрывности по неявной схеме, что придает системе уравнений эллиптический характер (см. Фромм [1963] и Руо [1967]). В нем нет ограничения на размер шага по времени, связанного со скоростью звука. Дальнейшие усовершенствования метода I E и его приложения указаны в работах Харлоу и Амсдена [1970], а также Харлоу с соавторами [1971]. [c.423] Вернуться к основной статье