Главная » Просмотр файлов » Дульнев Г.Н., Парфенов В.Г., Сигалов А.В. Применение ЭВМ для решения задач теплообмена

Дульнев Г.Н., Парфенов В.Г., Сигалов А.В. Применение ЭВМ для решения задач теплообмена (1185899), страница 25

Файл №1185899 Дульнев Г.Н., Парфенов В.Г., Сигалов А.В. Применение ЭВМ для решения задач теплообмена (Дульнев Г.Н., Парфенов В.Г., Сигалов А.В. Применение ЭВМ для решения задач теплообмена.djvu) 25 страницаДульнев Г.Н., Парфенов В.Г., Сигалов А.В. Применение ЭВМ для решения задач теплообмена (1185899) страница 252020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

расщепление по пространственным переменным приводит к погрешности решения, имеющей порядок 0 (Лт), и поэтому допустимо при достаточно малых шагах по времени Лт. Для получения локально-одномерной схемы достаточно провести дискретизацию задачи (3.79), (3.8!) по пространственным переменным с использованием неявных схем: прн!=-! (1а, т=Т»п,т', при!) ! б~:т =ш.'.ю.

(3.8'» м' — иl = — (ы', ь~ — 2ш,', „, 1- ш~,, „,. ~1, (3,901 зт ь» ш< т — бл (3.91) Рассмотрим структуру получившейся системы конечно-разностных уравнений и методику ее решения. Система (3.88), (3.89) для К, „, фактически распадается на не связанные между собой подсистемы, в каждую из которых входят только неизвестные, принадлежащие какой-либо из <горизонтальных прямых». Эти подсистемы решаются путем прогонок «по горизонталям» в направлении оси х, причем на каждом шаге по времени такие прогонки выполняются М раз; и = 1, ..., М.

Аналогично система (3.90), (3.91) распадается на «вертикальные» подсистемы, которые решаются прогонками в направлении у, которые выполняются Л' раз. Таким образом, для определения значений температуры и„',,„ на новом временном слое сначала на основе распределения температуры предыдущего временного слоя и'„ „, прогонками в направлениих находится промежуточное распределение д'„ „„ не имеющее самостоятельного значения, а затем на основе этого промежуточного распределения с помощью вертикальных прогонок вычисляется окончательное распределение нового временного слоя. Процедуре составления системы конечно-разностных уравнений локально-одномерной схемы целесообразно дать следующую физическую интерпретацию. На первом этапе область заменяется набором теплоизолированных между собой горизонтальных стер кней (рис.

3.18, а), для каждого из которых методом баланса записывается соответствующая неявная конечно-разностная схема, учитывающая граничные условия задачи на вертикальных границах х -- 0 и х = 1«как граничные условия для торцов стержня. Подчеркнем, что при составлении уравнений ба.. нса для нижнего и верхнего горизонтальных стержней их боковой теплообмен со средой учитывать не надо, т. е. адиабаты в направлении х проходят и по границам д -О, ч = !». Поэтому система уравнений для первого н последнего го- 1»1 ризонтальных рядов элементарных объемов (л»= ! и т=- М) полно«тью идентична системе для внутреннего ряда. На втором этапе аналогичным путем составляются конечно-разностные уравнения для в«ртикальных стержней (рис.

3.!6, 6). Для простой области прямоугольной формы описанная процедура составления разностной схемы на основе ее физической интерпретации не сильно облегчает работу по сравнению с формальным путем, однако для более сложных областей (см. рис. 3. (2) она оказывается весьма полезной и позволяет избежать ошибок. а) б/ У Миайпл»ы л доадпл»ы Ряс 3.!б Приведенное выше обоснование допустимости применении конечно-разностной схемы расщепления носило полукачественный характер. Однако можно провести и строгое математическое доказательство наличия у схемы (3.88) — (3.9!) свойств аппроксимапии порядка 0 (Лт г)»„'-'- й,',) и безусловной устойчивости (см., например, (24!). Отметим, что прн рассмотрении свойства аппроксимации вводится специальное понятие так называемой суммарной аппроксимации )24), которое заключается в следующем. Каждая из промежуточных систем разностпых уравнений (3.88) или (3.90) в отдельиосзи не обладает свойством аппроксимапии.

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

Заметим, что мощность внутренних источников»( при расщеплении уравнения теплопроводности можно «относить» либо к одному из направлений, как это было сделано выше, либо распределять с некоторыми весовыми коэффициентами между от- дельными направлениями. Например, для трехмерной задачи при составлении каждого из одномерных уравнений можно записать в него член д '3 )4роме локально-одномерной существуют и другие экономичные схемы. В частности, для двумерных задач получила распространение схема переменных направлений (4, 241. 4 зхь ЛРОГРАммнАЯ РеАлизАция численнОГО Решения МНОГОМЕРНЫХ ЗАДАЧ С ПОМОЩЬЮ ЛОКАЛЬНО-ОДНОМЕРНОЙ СХЕМЫ В качестве примера рассмотрим программу для решения по локально-одномерной схеме нестационарного трехмерного уравнения теплопроводности дтя параллелепипеда (рис.

3.17): дТ ( д»Т доТ си'Т дт 1 д»о ' дсо ' д»о (3.92) с граничными условиями дТ (3.93) х„.-=. (л — 1) йх, У„=(т — 1) йс, со — (Я вЂ” 1) )Г». Локально-одномерная схема при просзранственной сетке равно- мерной по каждой из координат х, у, г имеет вид: уравнения для первой промежуточной сеточной функции 1';, л ! аохлх сра»".

и чо»Л» У..,о — (1+ + . 2хдт/ ' >. А» ~'чо, ср à — — — и,, ) =-О. 2Ь(, 3 Ат ио» /4» ср -1- — ( — + — и,,о1=0, и=2 ..., У вЂ” 1, х (, 3 Дт о'"'/ о а1» Их сра 1 Г Г ( Х 2ХДт / ~""' о+ , 41 ах ао 1'4. сп 2» 1 3 дт 123 н начальным условием Т(х, у, г, т) 1»=.о Т,.

(3.94) ПРи численном Решении вводитсЯ сещ1чньи фУнкциа и'„, „, ы соответствующая температуре Т (х„у„, го, т,): 1 2 3 4 5 8 7 8 3 19 11 12 13 14 15 !6 !7 !8 13 29 21 22 23 24 25 28 27 28 28 39 31 32 33 34 35 36 37 38 39 49 41 42 43 44 48 47 48 48 59 51 52 53 ПРОГРАМ)й РАСЧЕТА НЕСТАЦИОНАРНОГО ТРЕХМЕРНОГО твакРАтутного поля ПАРллеелкпипЕ(Ы по ЛОКАЛЬНО-ОДНОМЕРНОЙ СХЕМЕ РХМЮВХОИ АЬР(8),(Щ(6),ХУ(19),)П(У(!9), «4[29),ВХ(29),с(29),Р(29),с(29), «с(29,29,29),вх(29),82(ге),и(29) !. ВВОД ИСХОДНЫХ ДАННйХ РХ,РХ,Р2 - РАЗМЕРЫ ПАРАЛЛЕЛЕПИПЕДА Аь,ся - теплспгсжщнссть, ав'ещия ТАплощя(ость ОЧ вЂ” РАСПРЕДЕЛЕНИЕ МОИНССТИ ВНУТРЕННИХ ИСТОЧНИКОВ АРР(6) - новиициенхы тхплсотдхчи нА ГРАнях (Х 9, Х РХ, Х 9, Х РХ, 2 9, « Р2) СЗ(6) - ПЛСТНССТИ ТЕПЛОВОГО ПОТОНА ПОВЕРХНОСТНЫХ ИСТОЧНИКОВ НА ГРАнЯХ ТЭ вЂ” НАЧАЛЬНАЯ ТЕМПЕРАТУРА НХ,КХ,Ж - ЧИСЛО ТОЧЕК ПО ОСЯМ ТАП вЂ” ИАГ ПО ВРЕМЕНИ ТМАХ вЂ” МАКСИМАЛЬНОЕ ВРЕМЯ зу - число щ(водов нА печхть в течение вины тч(зч) - мОменты ВРемени ВИВОдА нА печАть КУ вЂ” ЧИСЛО СЕЧЕНИИ, ПЕРПЕНДИКУЛЯРНЫХ ОСИ 2, в котоРых выводятся нА печхть тюя)ХРАТХРЫ КУЧ(КУ) - НОМЕРА СЕЧЕНИИ ВЫВОДА НА ПЕЧАТЬ НЕАЭ 1, РХ,РХ,Р2,АЬ,СН,ОУ 1 РОЯМАТ(6Р!9.3) РЯ1ИТ З,ЭХ,ЭХ,РХ,АЬ,СН,ОУ 2 РОНМАТ(' РАЗМЕРЫ:',ЗС19.3/ «' А(~',С!9.3,' СН ',С19.3,' ОЧ ',С19.3) НЕАЭ 1, (АЬР(Н),Н 1,6),(Ф(Н),н 1,8) РЯ!ЬТ 3,(А)Р(н),Н 1,6),(ЧЗ(И),И 1,8) 3 РОНМАТ(' КЮЭФИ~ИЕЗПЫ ТЕПЛООТДАЧИ:'/6319.3/ «' ПЛОТНОСТИ ТЕПЛОВОГО ПОТОКА:'/6019.3) яехэ 1, тэ,тАО,ТИАЛ РЯХИТ 4,ТЭ,ТАС,ТМАХ 4 РОЯМАТ(' Те ',РВ.З,' ТАС ',С!9.3,' ТМ(С( ',С19.3) ВЕАР 5, НХ,ИХ,ИЗ,ЪЧ,КУ 5 РОЯМАТ( 1915) РН1ИТ 6,ИХ,НХ,Н2 6 РОНМАТ(' ЧИСЛО УЗЛОВ:',315) НЕАЭ 5, (КУУ(И),Н 1,КЧ) РН1ЬТ 7,(КЧЧ(н),И 1,НУ) 7 РОНМАТ(' СЕЧЕННЯ ВЫВОДА:'/(1915)) ЯЕАР 1, (ТУ(Н),н (,ЗУ) РЯ1ЬТ 8,(ТУ(И),И 1,ЗУ) 8 РОЯМАТ(' МОМЕНТЫ ВРЕМЕНИ ВЫВОДА.''/(6С19.3)) 2.

ЗАПАНИЕ НАЧАЛЬНОГО РАСПРЕЛЕЛЕННЯ ТЕМПЕРАТУРЫ РО 1Э К 1,Н2 РО !9 М 1,ИХ ЭО 19 И !,ИХ !9 С(И,М,К) ТЭ 3. ВЫЧИСЛЕНИЕ ПОСТОЯННЫХ ВЕЩРИИ( Рис. 3.17 уравнения для второй промежуточной сеточной функции (Р'„ (ниже для сокрашения числа формул приводим лишь уравнения для внутренних узлов, а уравнения для граничных узлов т 1 н л( == М запишутся аналогично уравнениям для К„, »„» при л 1, л й() срд„' ( х! /д«ср Л)з дт уравнения для искомой сеточной функции и'„, „,, ы соответствующей температуре на /-ч временном слое, и«.»!.»«-! '(2 / ° ° /пп, «. » цл.~п,» ! + лдт / (!» / Чс ср — ' + )Р',.) =0. ~(з лт На каждом шаге по времени для расчета разностного решении и,', „, » сначала «прогонками по направлениям, пар,!,юельным осн «», определяется первая промежуточная сеточная функция 1»'„„, !атем «прогонками по направлениям, параллельныч оси 1(», вычисляется вторая промежуточная сеточная функция (т'„»,, » и, наконец, «прогонками по направлениям, параллельным оси г», находит.

ся искомое решение и„' „,, „. При этом сеточная функция 1'„', „ «вписывается на местофункции и,', „, », функция («",, „, » — на ме сто '»",»,„„а и,', „, » на место («';, „, ». Поэтому для хранения ссгочной функции температур используется лишь один трехчерный массив () размерами !» М '' К. Структура программы расчета трехмерной задачи во чногоч похожа на рассмотренную в э 3.5. структуру програччы д.!я одномерной задачи.

Основное отличие состоит в организации расчетов внутри цикла по времени. Каждый шаг по времени состои! нз трех частей, в которых определяются сеточные функции К,', » (операторы 92 — 105), »Г,', „,, (операторы 106 — 121), и,'...,, (операторы 122 — 137).

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

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

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