Главная » Просмотр файлов » Роуч П. Вычислительная гидродинамика

Роуч П. Вычислительная гидродинамика (1185913), страница 112

Файл №1185913 Роуч П. Вычислительная гидродинамика (Роуч П. Вычислительная гидродинамика.djvu) 112 страницаРоуч П. Вычислительная гидродинамика (1185913) страница 1122020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 112)

Затем все точки на сторонах, в которых 5 обращается в нуль (будем обозначать их звездочками), графопостроитель соединяет прямымн, как показано на рис. 7.5, б, и просматривается следующая ячейка. Г [/ (Уст+а) с (7+т У+т) Ь (у,а) а (г+т, 7) Рис. 7,5, К постросигпо изолиний иа линейном графопостроителе. В этой простой процедуре предусматривается следующая возможность.

Если изолиния должна проходить через угловую точку ячейки, скажем через точку (7+ 1, У), то эту точку необходимо считать за две точки обращения величины 5 в нуль — одну на стороне и и другую на стороне с(. Поскольку эти две точки слились в одну, то по логике программы они должньг соединяться линией нулевой длины. Приняв такое правило подсчета числа точек изолинпи, легко показать, что на одну ячейку может приходиться либо О, либо 2, либо 4 такие точки. Действительно, при обходе ячейки по замкнутому контуру (а, Ь, с, с() должно быть четное число перемен знака величины 5, а это н означает, что число точек, в которых 5 = О, должно быть четным. Если в ячейке обнаруживаются четыре такие точки, то они соединяются между собой всеми шестью возможными способами, как показано на рис.

7.5, в. На общей картине изолинийтакое построение вполне приемлемо, посколысу невозможно определить форму поверхности внутри ячейки только по значениям в ее угловых точках; эта поверхность может быть ~ребнем, соединяющим две вершины, расположенные по диагонали ячейки; выемкой, соединяющей две впадины, расположенные по диагонали ячейки, и, наконец, поверхностью типа седла. Действительная форма поверхности в ячейке в этом случае подсказывается ее поведением в соседних ячейках, так что неопределенность„имеющая место на рис.

7.5, в, вполне приемлема. 7.3. Обработка информации 498 Нужно особо рассмотреть случай, когда 5 = 0 (с некоторой точностью) в двух соседних узловых точках, скажем в точках (7, 1) и (7 + 1, У). Тогда 5 = 0 вдоль всей стороны а и единственной точки здесь выбрать нельзя. В этом случае на такой линии вообше ие выбирается точек для графопостроителя. В случае, когда область 5 = 0 простирается на много ячеек, такое рассмотрение приведет к обрисовке края, а не внутренности области, что опять вполне приемлемо.

30. Не стройте графиков с более высоким разрешением, чем полученные численные результаты. Вместо линейной интерполяции можно спокойно применять и интерполяции более высокого порядка, обычно используемые при обработке экспериментальных данных, поскольку эти произвольные процедуры не оказывают обратного влияния на основной процесс расчета задачи. Однако они часто вводят взаблуждение, так как сам расчет выполняется не со столь высоким разрешением. Графики с линейной интерполяцией представляют информацию наиболее непосредственно, в то время как сглаженные кривые внешне более привлекательны. Квадратичная интерполяция соответствует схемам второ~о порядка точности.

Особые вопросы возникают при изображении отрывных течений. Положение точек отрыва и вторичного присоединения потока можно определить только при помаши полуаналитических методов (Ундервуд 119б9]). Если такие точки нужно нанести, то надо указать, что их нашли экстраполяцией н т. п.; см. Лаван с соавторами [1969), Роуч и Мюллер 119701, Шаве и Ричарде (1970), (Это же замечание справедливо и для определении положения точек отрыва и вторичного присоединения потока в физическом эксперименте.) В настояшее время в вычислительных центрах для пользователей разрабатываются библиотеки программ построения трехмерных графиков и графиков изолиний. Выбор программ для построения графиков можно поручить профессиональному программисту, если имеется такая возможность.

Во многих вычислительных центрах имеются возможности для изготовления кинофильмов картин течений в виде трехмерных графиков или изолиний, причем каждый кадр дает картину течения в отдельный момент времени. Такие кинофильмы могут быть высокоэффективным средством для передачи динамики иестационарных процессов Пионерами в этой области, так же как и в разработке численных схем, являются исследователи из Лос-Аламосской лаборатории. Опубликованные в отчетах кадры цз таких кинофильмов помогают представить развитие течения Рис. 7.6.

Последовательность кадров кинофильма о нестационарном обтекании препятствия. В изометрической проекции показана поверхность воды вблизи сваи квадратного сечения. Расчет проведен в Лос-Аламосской лаборатории У. Николсом, использовавшим трехмерный алгоритм метода маркеров и ячеек (данные любезно предоставлены С.

Хергом из Лос-Аламосской лаборатопиий зоо 7,о. Обработка информации во времени, как это показано на рис. 7.6. Примеры можно найти в работах Херта [1965], Харлоу и Фромма [1965], Томана и Шевчика [1966], Донована [1968, 1970] и Фромма [1970а].

Лакс [!969] приводит пример, который демонстрирует преимущества вывода информации в виде кинофильма. Длинные волны на воде и некоторые волновые явления в физике плазмы описываются уравнением Кортевега — де Вриза ди ди д'и — + и — + — =-О. дт дх дх' Существование его решения типа уединенной волны открыли Крускал и Забуский, наблюдавшие образование уединенных волн цри изучении кинофильма картины течения. Затем прн более тщательных расчетах это явление было выделено и в конце концов было найдено точное аналитическое решение. Таким образом, при помощи вычислений удалось обнаружить ранее неизвестное свойство нелинейного уравнения в частных производных; это — возможность, о которой мечтал фон Нейман.

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

Из других величин представляют интерес источниковый член в уравнении Пуассона для давления, коэффициент давления торможения, который может служить неким индикатором влияния вязкости (Бургграф [1966], Макано и Хын [1967], Роуч и Мюллер [!970]), дисснпативная функция, энтропия, относительная величина диффузионных и конвективпых членов. Аллен [1968] строил графики величин, показывающих отклонение решения от решения уравнений пограничного слоя. В стационарных задачах с отрывом потока линия тока ф=О, так называемая разделяющая линия тока, играет особую роль, поскольку она отграничивает область возвратного течения, Чтобы подчеркнуть это, ее можно наносить специальными символамн или выделять каким-либо другим способом.

Даусон и Маркус [1970] выделяли линию т[т = О, просто строя две лишш тока к(т = 0 о е, где е много меньше всех расстояний между лнннями тока (они брали е = 0.001 У Е). Поскольку скорости в области возвратного течения малы, обычно рекомендуется брать различные шаги изменения функции тока на графике Ак)т для кр ) 0 и для чт (0 (Роуч и Мюллер [1968], Аллен и Чен [1970]). Направление течения может также обозначаться стрелками, нанесение которых предусмотрено на большинстве линейных графопостроителей.

Можно также менять длину стрелок, поскольку 73.2. Графики и кинофильмы карп н течений Рис. 7.7. Примеры построения линий тока и линий направления течения. На рис. и н б линия тока чр=О отмечена крестиками, выше этой линии ЬО = 0.17, ниже Лю = 0.00!27, а — линии тока при обтекании обратного уступа несжимаемой жидкостью, Ке=!0; б — линии тока и направления течения несжимаемой жидкости около обратного уступа, гте = 1О; в — направление течения в донной области при обтекании обратного уступа сжимаемой жидкостью, Ме = 2.24, Ке ЗОО, у = 1,4. (Иэ работы Роуча н Мюллера )1968).) такая стрелка расположена вдоль линии тока, она показывает направление течения, а ее длина показывает абсолютную величину скорости. Примеры построения линий тока приведены на рис.

7.7, а, б. Картину линий тока, полученную в результате вычислений, можно сравнить с экспериментально наблюдаемой Рис. 7.8а. Линии тока при нестационарном обтекании кругового цилиндра для наблюдателя, движущегося вместе с обтекаемым телом (система координат связана со стенками аэродинамической трубы). Из работы Томана и Шевчика [1966). Ряс. 7.86. Линии тока при нестационарном обтеканви кругового цилиндра для наблюдателя, движущегося вместе с набегающим потоком. (Из работы Томана и Шевчика (1966).) 503 7.3.2. Графики а кинофильмы картин течений картиной течения, которая получается по методике Прандтля с нанесением алюминиевого порошка на поверхности воды. При использовании в расчетах физических переменных функцщо тока чр можно вычислить при помощи криволинейного пни теграла чр, — тра = ~ (и агу — о с(х), который рассчитывается по квадратурпой формуле, соответствующей точности конечно-разностной схемы для определения составляющих скорости.

Так, Рис. 7.9. Положения ударных волн при ааиедлениом движении осесииметрняиого тела: при 1=200 мкс М=2.34 (слева), при )=400 мкс М=1.32 (в неитрс) и при Г = 500 мкс М = 0.80 (справа). (Иа работы Уилкниса 119691.) схемам второго порядка для уравнения переноса соответствует квадратурная формула Симпсона, а не формулы более высокого порядка.

Характеристики

Тип файла
DJVU-файл
Размер
12,43 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее