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

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

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

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

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

В случае же течения сжимаемой жидкости уравнение энергии необходимо для решения остальных уравнений, а функция тока в нестационарном случае не определяется. Здесь приходится иметь дело с системой четырех дифференциальных уравнений в частных производных '). ') Если предположить, жо коэффпппепты переноса н удельные тепло. емкостп постоянны, то уравпекия удобно записать через вихрь и эптропкю, кзк зто слолзко в рзоото Цянь Сюэ.сопя [1958].

далее, для более ограни. чеккото класса задач [без учета вязкости ялп эпплопроводкостк и прп отсут. стаяв ударных волн) можно считать знтропкю погтояяяой, что ведет к исключению одной кскомой фуякпкя, Одпзко этот подход яе кспользоввлся тт~яроко ярк численном решении задач газовой лппзмпкк. Согласно другому подходу, развитому в работе Гольдттпз с соавторами [!969), уравнение экергкп, включающее члены с теплопроводкостью. звмеяяют уравнением перекосз энтропии и таким образом жертвуют сохрзнеквем эпергяв для сохраненття энтропии. 4.3. Традиционная форма уравнений З!б 4.2. Традиционная форма уравнений Уравнения неразрывности и количества движения для плоского течения мы приведем в их традиционном размерном виде (см. Шлихтинг [1968]), а уравнение энергии запишем так, как это сделано в книге Овчарека [1964] ') (см.

также любой обычный курс газовой динамики, например Шапиро [1953], Липман и Рошко [1957], Чепмен и Уолкер [197!]). Мы рассматриваем совершенный газ, т. е. считаем, что внутренняя энергия и явчяется только функцией температуры Т, но свойства газа могут быть переменными (черточками сверху, как н ранее, обозначаются размерные величины). Предполагается, что объемные силы отсутствуют, а объемная вязкость и учитывается (до не. которых пор) в коэффициенте л = и — а/з)ь.

Итак, указанные уравнения будут такими: до+а (ру)=б, д[ (4.!) р — "+У.й — 1) (Р У)=б. П[ (4.2) (4.3) (4.4) ~) Вместо малоизвестной у нас книги Овчарека [!964], на которую неоднократно ссылается автор, читатель может каждый раз обращаться к отечественным источникам, например к следующим книгам: Кочин Н. Е., Кибель и. й., Розе Н. В. Тсорстнческая гидромеханнка.

В 2-х частях. — 6-е изд, асправл. и доп. — Мл гэизматгиз, 1963; Седов Л. И. Механика сплошной среды. В 2-х томах. — З.е изд., перераб. — М.: Наука, 1976. — Прим. ред. Уравнения для сверхзвукового течения невязкого газа имеют гиперболический тнп. Конечно-разностные уравнения должны в определенной мере учитывать область зависимости исходных дифференциальных уравнений, что приводит к знаменитому условию устойчивости Куранта — Фридрихса — Леви. Наиболее серьезные вычислительные трудности, присушие расчету сверхзвуковых потоков, связаны с наличием ударных волн. В пределе при больших числах Рейнольдса ударные волны являются разрывалли в решении. Наибольшие усилия специалистов в этой области были направлены на то, чтобы научиться «размазывать» эти разрывы, сохраняя при этом точность расчета на некотором расстоянии от разрывов.

Мы рассмотрим также некоторые методы с выделением скачков, где эти разрывы сохраняются. 4ой Консервативная форма уравнений 317 +.3. Консервативная форма уравнений Уравнения движения невязкой сжимаемой жидкости в консервативной форме были выведены Курантом и Фридрихсоы [1948), однако практически их впервые использовал Лаке [1954) для построения консервативной разностной схемы. Если традиционные дифференциальные уравнения преобразованы таким образом, что основными искомыми переменными становятся консервативные величины р, ри, ро и Е, (величина, которая будет определена ниже), то применение к таким уравнениям консервативных конечно-разностных схем обеспечивает сохранение массы, количества движения и энерп1и, Соотношения Рэнкииа — Гюгонио для прямого скачка') основаны только на этих законах сохранения и не завися~ от деталей внутренней структуры скачка, Отсюда следует, что все устойчивые аппроксимирующие консервативные разностные схемы, примененные ') Некоторые расчеты (Кенпер [19?Об]) были проведены с заменой уравнения (4.4) уравнением для давления др Г др др = = — [ и — + д — + ур??[, д1 дх ду (4,6) где т = с;/с.

— отношение удельных теплоемкостей. При выводе этого урав- нения использовалось уравнение состояния совершенного газа Р= рот. (4.6) Уравнение (4.4) является более общим и более употребительным, так что мы в этой книге будем иметь дело только с ним, ') Значением торможения или полным значением называется то значение какой-либо газодинамической величины, когорое получилось бы при переходе течения обратимым образом в состояние покоя ') См., например, Овчарек 11964] или любой другой куре газовой динамики. Здесь ~) /)/И = д/дг+ йд/дх+ пд/ду — субстанциональная производная, ?Т = Ч Ч вЂ” дивергенция скорости, е, = е + Рз/2— внутренняя энергия торможения') на единицу массы, г) — вектор потока тепла, Т вЂ” тензор полных напряжений.

Для определения этих величин необходимы дополнительные соотношения. В компоненты тензора Т входят как давление, так и вязкие напряжения. Гравитационная постоянная, используемая Шлихтингом [1968], в приведенную систему пе введена явно, но неявно она включается в принятую здесь систему единиц. 4.3. Консервативнал форма уравнений 3! в к уравнениям в консервативной форме, удовлетворяют соотношениям Рэпкина — Гюгонио и, следовательно, дают правильные условия на разрыве '). В работе Лонгли [1960] были опробованы четыре различные разностные схемы, и прн этом оказалось, что пз-за использования уравнений в консервативной форме все они дают правильные значения скорости скачка. Гари [1964] показал, что применение схемы Лакса — Вендроффа к уравнениям в неконсерватпвной форме приводит к значительным погрешностям в величине скорости скачка (хотя волна разрежения рассчитывается несколько точнее).

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

3 (разд. 3.1.3). При такой интерпретации нет необходимости в каких-либо предположениях о непрерывности функций. Поэтому интегральные формы предпочтительнее, н многие полагают, что все физические законы следует записывать в интегральной форме, Конечно-разностные аналоги уравнений Навье — Стокса в интегральной форме выведены в работах Лллена [1968] и Рубина и Прейзера [1968, 19?О]. Для устранения громоздкости опустим черточки над размерными величинами и в уравнениях, связывающих размерные величины, будем ставить над знаком равенства букву р, Уравнение неразрывности (4.!) уже записано в консервативной форме. Консервативную форму других уравнений можно получить путем преобразований, подобных следующим.

Рассмот- '1 Отметим, что использование исходных уравнений в консервативной форме само по себе пе обеспечивает сохранение массы, колгпстпа движении п энергии; необходимо также, чтобы сам конечно-разностиый метод был консервативным. 4.8. Конеервотивнан форма уравнениа рим следующие члены уравнения количества движения в направлении х: д (рн) д (рии) д (рио) д) дх ду др д(ри] д(ро) дн ди ди =и —.+и — +и +р — +ри — +ро — = д( дх ду д) дх ду = и [ д) + Ч (рЧ)1 + р и ' (4'7) Выражение в квадратных скобках равно нулю в силу уравнения неразрывности (4.1); поэтому ~~и о д(ри) + д(рии) + д(рио) " д(ри + Ч ((ри) ц Г)( д( дх ду д) Подставляя (4.8) в (4.2), получаем уравнение количества движения в направлении х в консервативной форме: = — д — Ч.((ри) Ч)+ д [2р д +о01+ д [)х (д +дх))' (4.9а) плп сокращенно О (ри) " др д( = д Ч При)Ч)+ ') +))2' (4.96) Уравнение количества движения в консервативной форме в на- правлении у выводится аналогично: + д [2Р д + И)1, (4,10а) (4.

1Об) — „+Ч ЦЧ)=Ь,, д) (4.1 1) где член 5) состоит из компонент градиентов давления н вязких напряжении, стоящих в правой части уравнения (4.9б). Теперь нли сокращенно — = — д — Ч Про)Ч)+):)в+)~~. д (ро) " дР В уравнениях (4.9) и (4.10) переменные и н о заменены консерватпвнымп переменнымп ри и ро. Теперь в уравнении количества движения в направлении х основной переменной является количество движения ри. Чтобы подчеркнуть это, перепишем уравнение (4.9б), введя обозначение ) = ри: 4 Э Канаараатааная форма урааненаа зао совершенно ясно, что уравнение для переменной / представляет собой уравнение переноса.

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

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

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

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

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