Cтраница 4
Его удобно использовать как в простейших моделях взрыва, так и в более сложных случаях, когда приходится решать уравнения в частных производных. В этом случае обычно приходится использовать разностную схему с искусственной вязкостью - которая существенно искажает форму упругой волны. Искусственная вязкость практически не влияет на зависимость положения фронта разрушения от времени и на величину a ( t), которая в основном определяется прочностью среды. [46]
Наиболее распространенным приемом построения однородных разностных схем в газодинамике является метод искусственной вязкости или, иначе, псевдовязкости. В этом методе в систему разностных уравнений газодинамики вводят некоторые дополнительные диссипативные члены, моделирующие действие реальной вязкости в ударных волнах. Введение искусственной вязкости позволяет приближенно описать скачок на фронте ударной волны гладким переходом на протяжении нескольких пространственных шагов сетки. При этом фронт ударной волны проявляется как область относительно быстрого изменения параметров газа. Скорость перемещения этой области совпадает со скоростью ударной волны, а изменения термодинамических величин в ней - со скачками на фронте ударной волны. Но протяженность самой области перехода существенно превышает реальную ширину фронта ударной волны, измеряемую несколькими длинами свободного пробега в газе. [47]
В главе рассматривается построение одномерных дискретных моделей, устанавливаются связи с соответствующими континуальными моделями. Проанализированы численные результаты по распространению одномерных волн в одно -, двух - и трехслойных пакетах. Для сглаживания ударно-волновых профилей использована линейная и квадратичная искусственная вязкость Неймана - Рихтмайера. Рассмотрена модификация схемы распада - разрыва, уменьшающая схемную вязкость. Приведены численные результаты по распространению одномерных волн в слоистых пакетах и моделированию их разрушения. [48]
![]() |
Линии тока и магнитные силовые линии при А 2 в плоскости xz. Угол между BOO и VQO равен 45.| Магнитные силовые линии обволакивают гелиопаузу. [49] |
На рис. 5.44 представлены профили плотности вдоль OCHZ. Разрывы, присущие картине течения разрешены с достаточной точностью, принимая во внимание малое число ячеек в одном из угловых направлений ( R х ф х Х0 99х21х 116), и без паразитных осцилляции. Отметим, что алгоритм не требует введения никакой искусственной вязкости. Ошибка в скачке плотности через большинство разрывов не превышает 3 - 5 % от точного значения. [50]