Разработка и исследование цифровой модели теплового потока при течении вязкой жидкости в канале с внешними нагревающимися элементами - Дипломная работа
Основные понятия теории течения жидкости. Создание математической модели распределения температурного поля в вязкой жидкости. Разработка цифровой модели изменения поля температуры в зависимости от: теплопроводности жидкости и металла, граничных условий.
Аннотация к работе
Такой процесс происходит в движущихся жидкостях и газах и всегда сопровождается передачей тепла посредством теплопроводности от частицы к частице при условии, что вся масса текущей жидкости не обладает одинаковой температурой. До тех пор пока мы рассматриваем области, расположенные внутри потока, и, следовательно, не обращаемся к процессам, происходящим на твердых поверхностях, которые ограничивают поток, или на его свободной поверхности, мы можем обе формы переноса тепла охватить одним понятием - теплопроводность в движущихся средах. В этой области поля скорости и температуры описываются уравнениями, похожими на уравнение, использующееся для задач двумерной теплопроводности. Количество переносимой теплопроводностью энергии, определяемое как плотностью теплового потока QT [Вт/(м2• К)], пропорционально градиенту температуры (закон Фурье): где ? - коэффициент теплопроводности, характеризующий его способность проводить теплоту, Вт/(м•К); знак минус указывает направление переноса теплоты в сторону снижения температуры. Диффузионный член в (2.5) удобно представить в виде: Для двумерной задачи, сформулированной в полярных системах координат, выражение для будет иметь вид: Плотности потоков соотносятся с градиентом ? следующим образом: Уравнение для стационарного медленного течения в канале записывается в виде где ср - удельная теплоемкость при постоянном давлении; k - теплопроводность жидкости.В программе, как уже упоминали раньше, предусмотрен выбор жидкости, с различной вязкостью ? (Па*с); материала кожуха, у которых отличается теплопроводность k (вт/К*м); изменение конфигурации трубы, тем самым придать жесткость канала и увеличить площади теплообмена. В зависимости какие параметры были применены при расчете задачи, получаем различные температурные поля. Решим задачу, используя разные материалы и изменяя жидкость. Так же можем сравнить распределение температуры для кожуха. На рисунке 6.1 и 6.3 приведены результаты нагрева кожуха для случаев, когда материал кожуха чугун и алюминий.Для выполнения дипломной работы был построен теплообменник, который состоит из металлического кожуха и цилиндрической трубы, по которой течет жидкость. Требуется исследовать, как выглядит температурное поле жидкости при изменении настраиваемых параметров. При выполнении поставленной задачи была: - создана математическая модель распределения температуры жидкости с заданной вязкостью в цилиндрическом канале;begin for j:=2 to M2 do for i:=2 to L2 do begin gam[i,j]:=cond; //коэфф-т диффузии Г end; //присваиваем произвольные значения cp:=1; //теплоемкости den:=1; //плотности dpdz:=-1; //градиенту давления t1:=STRTOFLOAT(Form1.Edit1.Text); // температура нагрева t2:=STRTOFLOAT(Form1.Edit2.Text); // температура стенки канала TEMPERGRAF := t2; if (Form1.N41.Checked = true) then trnagrev1 := STRTOFLOAT(Form3.STRINGGRID2.Cells[7,2]) else if (Form1.N61.Checked = true) then begin trnagrev1 := STRTOFLOAT(Form3.STRINGGRID2.Cells[4,2]); MESSAGEBOX(handle, PCHAR(" Температура на границах внутренней трубы должна " #10 "быть не равной и меньше температуры на участке нагрева!"), PCHAR("Ошибка"), MB_OK MB_ICONSTOP) end else begin if STRTOFLOAT(Form1.Edit3.Text) >= STRTOFLOAT(Form1.Edit1.Text) then begin MESSAGEBOX(handle, PCHAR("Температура жидкости должна быть меньше температуры стенки канала!"), PCHAR("Ошибка"), MB_OK MB_ICONSTOP) end else if (STRTOFLOAT(Form1.Edit2.Text) - STRTOFLOAT(Form1.Edit1.
Введение
В учении о теплообмене рассматриваются процессы распространения теплоты в твердых, жидких и газообразных телах. Эти процессы по своей физико-механической природе весьма многообразны, отличаются большой сложностью и обычно развиваются в виде целого комплекса разнородных явлений.
Перенос теплоты может осуществляться тремя способами: теплопроводностью, конвекцией и излучением, или радиацией. Эти формы глубоко различны по своей природе и характеризуются различными законами.
Первая форма теплообмена представляет собой перенос тепла посредством теплопроводности. Она характеризуется тем, что ее возникновение обусловлено наличием вещественной среды, и тем, что теплообмен совершается только между непосредственно соприкасающимися частицами тела. Этот процесс можно себе представить как распространение тепла от частицы к частице.
Второй формой теплообмена является перенос тепла посредством конвекции. Этот перенос осуществляется в тех случаях, когда частицы тела меняют свое местоположение в пространстве и при этом выполняют роль носителей тепла. Такой процесс происходит в движущихся жидкостях и газах и всегда сопровождается передачей тепла посредством теплопроводности от частицы к частице при условии, что вся масса текущей жидкости не обладает одинаковой температурой. До тех пор пока мы рассматриваем области, расположенные внутри потока, и, следовательно, не обращаемся к процессам, происходящим на твердых поверхностях, которые ограничивают поток, или на его свободной поверхности, мы можем обе формы переноса тепла охватить одним понятием - теплопроводность в движущихся средах. Если же рассматриваются и твердые стенки то в общем случае между стенками и движущейся средой наблюдается теплообмен, который обусловлен тем, что частицы среды, соприкасающейся со стенкой, воспринимают от нее тепло и уносят его с собой. Теплообмен между средой и стенкой называют теплоотдачей.
Третья форма теплообмена - перенос тепла посредством излучении. Эта форма характеризуется тем, что некоторая часть внутренней энергии тела преобразуется в энергию излучения и уже в такой форме передается через пространство. Встречая на своем пути другое тело, лучистая энергия полностью или частично вновь превращается в теплоту.
За последние десятилетия интерес к математическому моделированию сложных физических процессов и необходимость в нем в нем заметно возросли. Этому в значительной мере способствует прогресс в развитии компьютерной техники, численных методов решения всех типов задач математической физики и реализуемых на этой основе математических моделей. Любая современная наукоемкая технология, так или иначе, использует результаты вычислительных экспериментов. Сказанное в полной мере относится к исследованиям в области динамики жидкости и тепломассообмена.
Для численного решения практических задач, связанных с теплопереносом, течением жидкости и другими аналогичными явлениями, требуется, как правило, интегрирование системы нелинейных дифференциальных уравнений в частных производных по пространственным координатам и времени. Хотя существуют численные методы для получения такого решения, задача написания и использования общих вычислительных программ для всех практически важных процессов тепломассопереноса достаточно трудна.
Важным классом задач, подобных задачам теплопроводности, являются задачи о полностью развитых течениях и теплопереносе в каналах. Вдали от входного участка обычно существует область, в которой продольная скорость и поле температуры демонстрируют некоторое особое поведение. В этой области поля скорости и температуры описываются уравнениями, похожими на уравнение, использующееся для задач двумерной теплопроводности.
По трубопроводам транспортируют газы и различные жидкости. Встречаются трубопроводы различных сечений, но чаще круглые и полукруглые. Часто к стенкам канала прикрепляют ребра для увеличения площади поверхности, за счет чего улучшается теплообмен. Присутствие ребер приводит также и к увеличения перепада давления по длине канала при том же расходе через его сечение.
Жидкость в трубопроводе подогревают для уменьшения ее вязкости и увеличения скорости перемещения. Важно проанализировать поведение жидкости в трубопроводе, изменение температуры от теплопроводности и количество прикрепленных ребер.
В данном случае исследование актуально именно потому, что определенные аспекты темы изучены не в полной мере и проведенное исследование направлено на преодоление этого пробела.
Целью работы является разработка и исследование цифровой модели теплового потока при течении вязкой жидкости в канале с внешними нагревающимися элементами.
Перед нами стоит ряд задач: - Создание математической модели распределения температурного поля в вязкой жидкости;
- Выбор метода реализации математической модели;
- Разработка алгоритмов реализации модели с помощью выбранных методов;
- Выбор программной среды для его реализации;
- Разработка программы по алгоритму с помощью выбранного программного обеспечения;
- Разработка цифровой модели изменения поля температуры в зависимости от: теплопроводности жидкости и металла, граничных условий;
- Проверка адекватности модели.
1. Постановка задачи
Решим задачу для теплообменника, который состоит из металлического кожуха и цилиндрической трубы, по которой течет жидкость заданной вязкости. Имеется некая труба, в которой сохраняется постоянная температура. В некотором участке эту трубу обхватывает нагревающийся металлический кожух. В нашем случае кожух нагревают с нижней части.
Теплообмен между нагретым кожухом и жидкостью может происходить только через поверхность трубы.
Для исходной области заданы граничные условия на внешней и внутренней границе кожуха. На нижней части кожуха задается постоянная температура T1 (идет нагрев кожуха), на внутренней стенке кожуха задана начальная температура T3. Задается температура жидкости T2, которая втекает в область нагрева (Рисунок 2).
Задается радиус трубы R2 и кожуха R1. Для обеспечения жесткости и увеличения площади теплообмена к внутренней стороне трубы присоединяют ребра жесткости, длиной половины R2.
Задается теплопроводность материала кожуха k и вязкость жидкости ?.
Изза симметрии канала расчетной областью будет половина кольца (Рисунок 3).
Рисунок 3 - Расчетная область канала
Требуется исследовать, как выглядит температурное поле жидкости при изменении: - температуры нагрева кожуха;
- теплопроводность материала кожуха;
- вязкости жидкости;
- конфигурацию трубы (изменение количество ребер);
- температуры стенки кожуха и жидкости.
2. Основные понятия теории течения жидкости
Жидкость - одно из агрегатных состояний вещества, материал, обладающий свойством текучести, т.е. сплошная легкодеформируемая среда, обладающая свойством течь. Основным свойством жидкости, отличающим ее от других агрегатных состояний, является способность неограниченно менять форму под действием касательных механических напряжений, даже сколь угодно малых, практически сохраняя при этом объем [4].
К основным параметрам, характеризующим свойства жидкости, принадлежат: Вязкость ?
Жидкости (как и газы) характеризуются вязкостью. Вязкостью, или внутренним трением жидкости, называется сопротивление, которое жидкость оказывает перемещению одного ее слоя относительно другого под влиянием внешней силы.
Механизм внутреннего трения в жидкостях и газах заключается в том, что хаотически движущиеся молекулы переносят импульс из одного слоя в другой, что приводит к выравниванию скоростей - это описывается введением силы трения.
Различают динамическую вязкость (единицы измерения: пуаз, 0,1Па·с) и кинематическую вязкость (единицы измерения: стокс, м?/с, внесистемная единица - градус Энглера). Кинематическая вязкость может быть получена как отношение динамической вязкости к плотности вещества и своим происхождением обязана классическим методам измерения вязкости, таким как измерение времени вытекания заданного объема через калиброванное отверстие под действием силы тяжести [3].
Плотность ?
Плотностью жидкости или газа называется количество массы, заключающейся в единице объема. Для определения плотности вещества нужно массу тела m разделить на его объем [3].
Размерность плотности - кг/м3 (система единиц СИ).
Теплоемкость Ср
Теплоемкость - физическая величина, определяющая отношение бесконечно малого количества теплоты ?Q, полученного телом, к соответствующему приращению его температуры ?T.
Единицей СИ для удельной теплоемкости является джоуль на килограмм-цельсий. Следовательно, удельную теплоемкость можно рассматривать как теплоемкость единицы массы вещества. На значение удельной теплоемкости влияет температура вещества [3].
Теплопроводность k
Теплопроводность - это молекулярный перенос теплоты между непосредственно соприкасающимися телами или частицами одного тела с различной температурой, при котором происходит обмен энергией движения структурных частиц (молекул, атомов, свободных электронов) [2].
Существует также понятие идеальной жидкости.
Идеальная жидкость - в гидродинамике - воображаемая (идеализированная) жидкость, в которой, в отличие от реальной жидкости, отсутствует вязкость. В идеальной жидкости отсутствует внутреннее трение, то есть нет касательных напряжений между двумя соседними слоями.
Моделью идеальной жидкости пользуются при теоретическом рассмотрении задач, в которых вязкость не является определяющим фактором и ею можно пренебречь. В частности, такая идеализация допустима во многих случаях течения, рассматриваемых гидроаэромеханикой, и дает хорошее описание реальных течений жидкостей и газов на достаточном удалении от омываемых твердых поверхностей и поверхностей раздела с неподвижной средой. Математическое описание течений идеальных жидкостей позволяет найти теоретическое решение ряда задач о движении жидкостей и газов в каналах различной формы, при истечении струй и при обтекании тел [8].
2.2 Общие характеристики течения в каналах
Хоть каналы могут иметь любую геометрическую форму, для простоты рассмотрим прямой канал прямоугольного поперечного сечения, изображенный на рисунке 2.2. Течение в основном направлено вдоль оси z, (назовем ее продольной координатой, в отличие от поперечных координат х и у). Течение вызывается градиентом давления ?р/?z, который обычно отрицателен. Давление практически постоянно в поперечном сечении и изменяется вдоль оси z.
Так как изменения скорости и температуры по оси z малы по сравнению с изменениями их по оси х или у, то очень часто пренебрегают вязким напряжением, вызванным градиентом ?w/?z, и переносом тепла за счет градиента ?Т/?z. При этом обычно член ?р/?z в уравнении для продольной составляющей скорости полагают равным ?р/?z, где р соответствует среднему значению давления в поперечном сечении. Течение в канале может быть стационарным или нестационарным [13].
2.3 Ламинарное движение
Это движение, называют также потенциальным (безвихревым) движением.
При наблюдении за движением жидкости в трубах и каналах можно заметить, что в одном случае жидкость сохраняет определенный строй своих частиц, а в других - перемещаются бессистемно.
Ламинарным называется слоистое течение без перемешивания частиц жидкости и без пульсации скорости и давления. При ламинарном течении жидкости в прямой трубе постоянного сечения все линии тока направлены параллельно оси трубы, при этом отсутствуют поперечные перемещения частиц жидкости. Ламинарное течение - это плавное, упорядоченное, регулярное движение, когда отдельные струйки жидкости, не перемешиваясь, как бы скользят друг по другу. При таком потоке существует лишь молекулярное трение между соседними струйками [13-16].
Ламинарное течение жидкости наблюдается при небольших скоростях ее движения. Внешний слой жидкости, примыкающий к поверхности трубы, в которой она течет, изза сил молекулярного сцепления прилипает к ней и остается неподвижным (рисунок 2.3). Скорости последующих слоев тем больше, чем больше их расстояние до поверхности трубы, и наибольшей скоростью обладает слой, движущийся вдоль оси трубы.
Бугорки и шероховатости обтекаются ламинарным потоком и не оказывают влияние на сопротивление трубы. При увеличении скорости движения толщина ламинарного слоя уменьшается.
При увеличении скорости потока траектории отдельных струек приобретают волнообразный характер, через некоторое время струйка исчезает, перемешиваясь с жидкостью. Характер течения при этом изменился. Траектории отдельных частиц приобретают хаотичный неустановившийся характер. Такое течение называется турбулентным - хаотичным, вихревым [16].
Переход от ламинарного к турбулентному режиму наблюдается при определенной скорости движения жидкости. Эта скорость называется критической ?кр. Для характеристики течения жидкости используют так называемое число Рейнольдса - Re [1]. С ростом Re в некоторый момент происходит потеря устойчивости движения, струйки перемешиваются, в потоке образуются хаотически пульсирующие вихри. Такое неупорядоченное движение вязкой жидкости становится турбулентным. При турбулентном режиме обычно резко растет сопротивление потока жидкости.
Значение этой скорости прямо пропорционально кинематической вязкости жидкости и обратно пропорционально диаметру трубы.
где ? - кинематическая вязкость; k - безразмерный коэффициент; d - внутренний диаметр трубы.
Входящий в эту формулу безразмерный коэффициент k одинаков для всех жидкостей и газов, а также для любых диаметров труб. Этот коэффициент называется критическим числом Рейнольдса Rekp и определяется следующим образом:
Как показывает опыт, ламинарное движение возможно при сравнительно невысоких числах Рейнольдса. Для труб круглого сечения число Рейнольдса Rekp примерно равно 2300. Для нашего случая, когда канал в сечении кольцо и поверхность контакта между жидкостью и стенкой больше, то и больше локальных возмущающих факторов число Рейнольдса будет Rekp =1600 [1].
Критерий подобия числа Рейнольдса позволяет судить о режиме течения жидкости в трубе. При малых значениях числа Рейнольдса ( ) наблюдается ламинарное течение, переход от ламинарного течения к турбулентному происходит в области , а при (для гладких труб) течение - турбулентное. Если число Рейнольдса одинаково для разных жидкостей, то и режим течения таких жидкостей в трубах разных сечений одинаков.
Формула, по которой определяется число Рейнольдса в канале любой геометрии [16],
Гидравлический диметр при конфигурации канала, когда сечение представляет собой кольцо [15]:
где ? - смоченный периметр, D и d - внешний и внутренний диаметр трубы.
2.4 Теплообмен
Теплообмен - это самопроизвольный необратимый перенос теплоты (точнее, энергии в форме теплоты) между телами или участками внутри тела с различной температурой. В соответствии со вторым началом термодинамики теплота переносится в направлении меньшего значения температуры. Теплообмен существен во многих процессах нагревания, охлаждения, конденсации, кипения и оказывает значительное влияние на массообменные процессы [19].
Движущиеся среды, участвующие в теплообмене и интенсифицирующие его, называются теплоносителями (обычно капельные жидкости, газы и пары). Известны два основных способа проведения тепловых процессов: путем теплоотдачи и теплопередачей.
Теплоотдача - теплообмен между поверхностью раздела фаз (чаще твердой поверхностью) и теплоносителем.
Теплопередача - теплообмен между двумя теплоносителями или иными средами через разделяющую их твердую стенку либо.
Различают три разных механизма распространения теплоты: теплопроводность, конвективный и лучистый перенос.
Теплопроводность - перенос энергии от более нагретых участков тела к менее нагретым в результате теплового движения и взаимодействия микрочастиц (атомов, молекул, ионов и др.). В чистом виде теплопроводность может встречаться в твердых телах, не имеющих внутренних пор, и в неподвижных слоях жидкостей, газов или паров. Количество переносимой теплопроводностью энергии, определяемое как плотностью теплового потока QT [Вт/(м2• К)], пропорционально градиенту температуры (закон Фурье):
где ? - коэффициент теплопроводности, характеризующий его способность проводить теплоту, Вт/(м•К); знак минус указывает направление переноса теплоты в сторону снижения температуры.
Численная характеристика теплопроводности материала равна количеству теплоты, проходящей через материал площадью 1 кв.м за единицу времени (секунду) при единичном температурном градиенте.
С повышением температуры теплопроводность жидкостей, за исключением воды, уменьшается, а для всех других тел увеличивается.
Конвективный перенос теплоты - перенос физической теплоты перемещающихся нагретых жидкостей, газов, паров. В наиболее распространенном случае, когда существен лишь перенос внутренней энергии, а переносом механических и потенциальных видов энергии можно пренебречь, плотность теплового потока за счет конвективного переноса составляет:
где w - вектор скорости текучей среды; ?, С, Т - плотность, теплоемкость и температура среды.
Плотностью жидкости или газа называется количество массы, заключающейся в единице объема. Для определения плотности вещества нужно массу тела m разделить на его объем. Размерность плотности - кг/м3 (система единиц СИ).
Теплоемкость - физическая величина, определяющая отношение бесконечно малого количества теплоты ?Q, полученного телом, к соответствующему приращению его температуры ?T.
Единицей СИ для удельной теплоемкости является джоуль на килограмм-цельсий. Следовательно, удельную теплоемкость можно рассматривать как теплоемкость единицы массы вещества. На значение удельной теплоемкости влияет температура вещества.
В большинстве случаев значения w, ?, С и Т потоков теплоносителей таковы, что в направлении движения конвективный перенос преобладает над теплопроводностью. Однако при малых скоростях течения высокотеплопроводных жидкостей (расплавов металлов) может наблюдаться обратное соотношение.
По мере приближения к твердой поверхности, где скорость вязких жидкостей стремится к нулю, qt и qk также становятся сравнимы по величинам. При ламинарном режиме течения в направлении, поперечном движению, конвективный перенос отсутствует.
3. Математическая постановка задачи
3.1 Обобщенная математическая постановка задачи
Применение первого закона термодинамики к бесконечно малому контрольному объему в твердом теле или в неподвижной среде приводит к известному уравнению теплопроводности [1].
Обозначив три пространственные переменные через декартовы координаты х, у, z и время через t, можем записать уравнение теплопроводности в виде где ? - плотность, с - удельная теплоемкость.
В этом уравнении величины ?, с, k могут зависеть от х, у, z, t и быть функциями температуры Т.
Уравнение теплопроводности является прямым следствием закона сохранения, представленного первым законом термодинамики, и пропорциональности плотности потока градиенту температуры. Существует множество других физических процессов, при которых соответствующая плотность потока некоторой величины пропорциональна градиенту этой величины и для которых существует закон сохранения. Отсюда следует, что эти процессы будут описываться дифференциальными уравнениями, аналогичными (3.1).
К подобным процессам можно отнести диффузию химических компонент, движение заряженных частиц в электромагнитном поле, течение в пористых материалах, потенциальные течения, перенос тепла и влаги в почве, а также полностью развитые течение и теплообмен в каналах. Построив вычислительную процедуру для решения уравнения (3.1), мы сможем применить ее и для любого аналогичного процесса, просто придавая новый смысл величинам Т, k и др. Например, можно интерпретировать Т как концентрацию, k как коэффициент диффузии и т.п. Удобнее работать с таким обобщенным дифференциальным уравнением, так как уравнение теплопроводности и другие аналогичные уравнения станут его частными случаями. В дальнейшем будем основываться на подобном обобщенном дифференциальном уравнении [2].
3.2 Дифференциальное уравнение для обобщенной переменной
Обозначим через ? зависимую переменную в обобщенном дифференциальном уравнении. Градиент ? приводит к соответствующему диффузионному потоку. Таким образом, плотность диффузионного потока Jx в направлении х задается в виде:
где Г - обобщенный коэффициент диффузии.
Обобщенное дифференциальное уравнение, представляющее закон сохранения для ?, может быть записано в виде где ? - величина, подобная объемной теплоемкости.
Сравнивая уравнение теплопроводности (2.1) с обобщенным дифференциальным уравнением (2.3), заключаем, что для того чтобы обобщенное уравнение описывало процесс теплопроводности, нужно сделать следующую замену: ? = Т, ? = ?с, Г = k, В общем случае для значения ? существуют соответствующие значения величин ?, Г. В простой задаче эти величины могут быть постоянными, но в общем случае нет ничего необычного в том, что они зависят от пространственных координат и времени, а также от самой переменной ?.
Обобщенное дифференциальное уравнение (2.3) записано в декартовой системе координат. Более компактно оно может быть записано:
где применено соглашение о суммировании по повторяющемуся индексу. Например, для трехмерного случая можно записать здесь х1, х2 и х3 - декартовы координаты (х, у и z) по трем направлениям. Будем называть три члена уравнения (2.5) нестационарным, диффузионным и источниковым членом соответственно.
В более общем виде выражение для плотности диффузионного потока (2.2) может быть записано следующим образом [3]:
где Ji - плотность диффузионного потока в направлении координаты хі.
Диффузионный член в (2.5) удобно представить в виде:
Для двумерной задачи, сформулированной в полярных системах координат, выражение для будет иметь вид:
Плотности потоков соотносятся с градиентом ? следующим образом:
Уравнение для стационарного медленного течения в канале записывается в виде где ср - удельная теплоемкость при постоянном давлении; k - теплопроводность жидкости. Левая часть (2.10) соответствует конвективному переносу в канале, а члены правой части описывают теплопроводность в жидкости. Так как обычно перенос тепла вдоль оси z очень мал по сравнению с переносом в поперечном его сечении, то последним членом в уравнении (2.10) можно пренебречь.
При простом полностью развитом течении в канале поперечные скорости и и v равны нулю, поэтому (2.10) упрощается и принимает вид
Уравнение (2.11) соответствует стационарной форме обобщенного дифференциального уравнения при следующей замене:
4. Методы решения
4.1 Выбор метода решения
Известны [16] несколько основных типов численных методов, применяемых для решения уравнений теплопроводности: - Метод конечных разностей;
- Метод контрольного объема;
- Метод конечных элементов.
Получить аналитическое решение уравнений описывающих явления гидродинамики, как правило, не представляется возможным. В связи с этим используются численные методы. Все методы включают два этапа.
- На первом этапе все дифференциальные уравнения сводятся к алгебраическим уравнениям относительно значений искомых переменных (скорости, давления, температуры и т.п.) в конечном числе точек в пределах области решения. Этот этап называют дискретизацией уравнений
- На втором этапе производится решение алгебраических выражений с помощью соответствующего численного метода.
Различные типы численных методов вычислительной гидродинамики используют свои способы дискретизации основных уравнений.
Каждый из упоминавшихся выше численных методов имеет свои преимущества и недостатки. Преимуществами метода контрольных объемов являются: · Простота применения к расчетным областям произвольной формы
· Наличие отработанных алгоритмов
· Ошибка дискретизации падает с ростом числа вершин
· Основное внимание уделяется балансу потоков через контрольные объемы
· Автоматическое выполнение законов сохранения
Граничные условия
Для численного моделирования требуется задание значений зависимых переменных или их нормальных градиентов на границах расчетной области. Граничные условия для нашей задачи были определены при постановке задачи.
Решение дискретных аналогов уравнений не может быть получено, если граничные условия не заданы. Граничные условия делятся на два типа - с известными значения на границе (задача Дирихле) и известными нормальными градиентами на границе (задача Неймана).
Начальные условия
Для нестационарных задач перед началом расчета необходимо задать значения всех искомых переменных во всей области потока в начальный момент времени.
Источники ошибок
В общем случае существует три основных первопричины ошибок: · Допущения, принимаемые при построении математической модели
· Аппроксимации при дискретизации уравнений
· Отсутствие сходимости в процессе решения
О принимаемых допущениях говорилось ранее - они связаны с неточным описанием реальных физических явлений используемыми уравнениями. К этой же категории можно отнести аппроксимацию геометрии и граничных и начальных условий.
Ошибки дискретизации являются результатом применения аппроксимирующих соотношений для расчетов потоков через грани ячеек и источниковых членов в пределах ячеек.
Под отсутствием сходимости в данном случае подразумевается невозможность получения решения дискретных алгебраических уравнений с желаемой точностью. (Стремление численного решения к точному при последовательном измельчении сетки также часто называется сходимостью). жидкость математический температура цифровой
4.2 Методы решения дифференциальных уравнений
Для решения задач математической физики, не поддающихся аналитическим решениям, используют численные методы, позволяющие аппроксимировать исходные дифференциальные уравнения. Развитие ЭВМ привело к появлению различных численных методов решения уравнений математической физики. Кроме того, за время развития вычислительных методов наметилась их специализация по отраслям приложения. Так, например, при моделировании глобальных атмосферных процессов и при решении акустических задач отдается предпочтение спектральным методам; методы конечных элементов применяются сейчас, в основном, для решения задач механики сплошной среды. Для решения проблем теплообмена и механики жидкости и газа используются как метод конечных разностей, так и метод контрольных объемов [12].
4.2.1 Метод конечных разностей
На рисунке 4.1 изображено разбиение одномерной расчетной области на узловые точки.
Конечно-разностный метод основан на замене производных дифференциального уравнения физического процесса в точке на конечные разности. Производные аппроксимируются разложением в ряды Тейлора. Для аппроксимации производной в точке 2 разложим исследуемую переменную ? в ряд Тейлора в точках 1 и 3
Пренебрегая остаточным членом Rn(x) обоих рядов, вычитая и складывая уравнения, получим:
Заменяя производные в исходном дифференциальном уравнении в соответствии с (2.6), получаем разностный аналог, записываемый в виде системы алгебраических уравнений.
Данный метод не учитывает физической картины процесса, а наличие остаточного члена в ряде Тейлора и предположение о полиномиальном характере изменения зависимой переменной, приводит к нежелательным последствиям, например, для случая экспоненциального изменения переменной. Этим методом решается значительное количество практических задач газовой динамики и диффузии (теплообмена), однако в ряде случаев конечно-разностный метод не подходит для решения уравнений математической физики. Так, для решения дифференциальных уравнений с разрывными коэффициентами в конечно-разностном методе применяют схему с искусственной вязкостью, метод характеристик и дивергентные схемы. Указанные разновидности конечно-разностного метода применяются для задач гидро- и газодинамики [2].
4.2.2 Метод сеток
Дифференциальные уравнения вида решаются преобразованием их в алгебраические уравнения, называемые дискретными аналогами. Эти уравнения содержат в качестве неизвестных значения ? в выбранных дискретных местоположениях, которые образуют сетку и называются расчетными точками. Вокруг каждой расчетной точки строится контрольный объем, и дискретные аналоги получаются интегрированием уравнения по таким контрольным объемам.
Основная идея метода контрольного объема легка и понятна и поддается прямой физической интерпретации. Расчетную область разбивают на некоторое число контрольных объемов, грани которых показаны штриховыми линиями. После в геометрический центр каждого контрольного объема помещается узловая точка.
Видно, что некоторая расчетная точка «сообщается» с четырьмя соседними через четыре грани контрольного объема. Одна из граней приграничного контрольного объема совпадает с границей расчетной области, а граничная точка помещена в центр грани контрольного объема.
Удобно представлять контрольный объем нулевой толщины для граничной точки.
Границы контрольных объемов должны располагаться так чтобы, они совпадали с разрывами в свойствах материала, источниковых членах, граничных условиях и др.
Существует большое число геометрических величин, которые имеют отношение не к расчетным точкам, а к граням контрольных объемов. Грань с номером i лежит между точками i-1 и i. Другими словами, грань имеет тот же номер, что и ближайшая к ней точка в положительном направлении оси координат, т.е. в направлении увеличения i или j. Частным следствием такого построения контрольных объемов и используемой нумерации является то, что грань I = 2 и точка I = 1 совпадают с левой границей расчетной области.
Границы контрольных объемов должны располагаться так чтобы, они совпадали с разрывами в свойствах материала, источниковых членах, граничных условиях и др.
В соответствии с местоположение расчетных точек (рисунок 4.2) в направлении осей параметры ? и у определяются индексами I и J. Индекс I увеличивается вдоль оси ?, а индекс J - вдоль оси у. Значение I=1 задает сеточную линию на левой границе, в то время как I=L1 указывает на правую границу. Аналогично J=1 и M1 соответствуют нижней и верхней границам. На самом деле во многих местах в программе левая и нижняя граница области называются границами I1 и J1 соответственно, а правая и верхняя - границами L1 и M1. Логичность такого обозначения очевидна. Для удобства введены некоторые дополнительные переменные, определяемые так:
Таким образом, сеточные линии I=2 и L2 являются первыми внутренними линиями сетки вблизи соответствующих границ, а линии I=3 и L3 - вторыми внутренними линиями. Подобными свойствами обладают линии J=M2 и МЗ. Таким образом, в диапазонах I=2, L2 и J=2, М2 задаются все внутренние расчетные точки. Зависимая переменная, например температура, в расчетной точке [I, J] будет обозначаться как Т[I, J]. Граничные значения Т обозначены как Т[1, J], Т[L1, J], Т[I, 1] и Т[I, M1].
Программа разработана для полярной системы координат ? - продольная координата, y,r - координаты в радиальном направлении. Разница между ними заключается в том, что координата r должна измеряться от полюса, в то время как координата y может быть выбрана произвольно.
Существует большое число геометрических величин, которые имеют отношение не к расчетным точкам, а к граням контрольных объемов. Грань с номером i лежит между точками i-1 и i. Другими словами, грань имеет тот же номер, что и ближайшая к ней точка в положительном направлении оси координат, т.е. в направлении увеличения i или j. Частным следствием такого построения контрольных объемов и используемой нумерации является то, что грань I = 2 и точка I = 1 совпадают с левой границей расчетной области.
Дифференциальные уравнения законов сохранения массы, энергии, количества движения интегрируют по каждому объему. В результате находят дискретный аналог дифференциального уравнения, в который входят значения в нескольких узловых точках.
Полученный подобным образом дискретный аналог выражает закон сохранения для конечного контрольного объема точно так же, как дифференциальное уравнение выражает закон сохранения для бесконечно малого контрольного объема. Одним из важных свойств метода контрольного объема является то, что в нем заложено точное интегральное сохранение таких величин, как масса, количество движения и энергия на любой группе контрольных объемов и, следовательно, на всей расчетной области. Это свойство проявляется при любом числе узловых точек, а не только в предельном случае очень большого их числа [1].
Именно этот метод и будет использован в дальнейшей работе.
4.3 Полярная система координат
Поставленная задача решена с применением полярной системы координат. Все вычисления произведены в полярных координатах без перевода их в декартовые.
Полярная система координат - двумерная система координат, в которой каждая точка на плоскости определяется двумя числами - полярным углом и полярным радиусом. Полярная система координат особенно полезна в случаях, когда отношения между точками проще изобразить в виде радиусов и углов; в более распространенной, декартовой или прямоугольной системе координат, такие отношения можно установить только путем применения тригонометрических уравнений.
Полярная система координат задается лучом, который называют нулевым или полярной осью. Точка, из которой выходит этот луч, называется началом координат или полюсом. Любая точка на плоскости определяется двумя полярными координатами: радиальной и угловой. Радиальная координата (обычно обозначается r) соответствует расстоянию от точки до начала координат. Угловая координата, также называется полярным углом или азимутом и обозначается ?, равна углу, на который нужно повернуть против часовой стрелки полярную ось для того, чтобы попасть в эту точку.
Определенная таким образом радиальная координата может принимать значения от нуля до бесконечности, а угловая координата изменяется в пределах от 0° до 360°. Однако, для удобства область значений полярной координаты можно расширить за пределы полного угла, а также разрешить ей принимать отрицательные значения, что отвечает повороту полярной оси по часовой стрелке.
Одной из важных особенностей полярной системы координат является то, что одна и та же точка может быть представлена бесконечным количеством способов. Это происходит потому, что для определения азимута точки нужно повернуть полярную ось так, чтобы он указывал на точку. Но направление на точку не изменится, если осуществить произвольное число дополнительных полных оборотов.
Углы в полярных координатах задаются либо в градусах, либо в радианах. Выбор, как правило, зависит от области применения. В навигации традиционно используют градусы, в то время как в некоторых разделах физики, и почти во всех разделах математики используют радианы.
4.4 Обобщенное дискретное уравнение
Значения обобщенной зависимой переменной ? сохраняются в расчетных точках. Дискретное уравнение связывает значение ? в одной расчетной точке со значениями ее в четырех соседних точках. Такое уравнение получается при интегрировании уравнения по контрольному объему, содержащему точку (i, j).
Типичный контрольный объем показан на рисунок 4.5. Проинтегрировав по нему уравнение (2.2), получим
В формуле (4.1) верхний индекс «0» обозначает известное значение ? в начале шага по времени ?t; J - плотность диффузионного потока через грань контрольного объема; - осредненный по контрольному объему источниковый член; A, ?V - площадь грани и контрольный объем.
Диффузионные потоки па гранях контрольного объема e и w могут быть рассчитаны следующим образом:
где De - проводимость между точками P и E, которая вычисляется по значениям Г в этих точках. Зная размеры ?xе- и ?xе , можно рассчитать De, по формуле
Аналогичным образом определяются значения D на других гранях. Проводимость зависит от площади грани контрольного объема .
Так как источниковый член S может зависеть от ?, ж
Вывод
6.1 Исследование влияния различных параметров на физические процессы
В программе, как уже упоминали раньше, предусмотрен выбор жидкости, с различной вязкостью ? (Па*с); материала кожуха, у которых отличается теплопроводность k (вт/К*м); изменение конфигурации трубы, тем самым придать жесткость канала и увеличить площади теплообмена. В зависимости какие параметры были применены при расчете задачи, получаем различные температурные поля.
Приведем несколько экспериментов. Решим задачу, используя разные материалы и изменяя жидкость.
Сравним результаты вычисления для воды и для мазута. На рисунке 6.1 и 6.2 можно увидеть табличные значения нагрева этих жидкостей. В обоих случаях материал кожуха чугун.
Так же можем сравнить распределение температуры для кожуха. На рисунке 6.1 и 6.3 приведены результаты нагрева кожуха для случаев, когда материал кожуха чугун и алюминий.
6.2 Построение и объяснение графиков
Рассмотрим распределение температуры, когда конфигурация трубы разная. Выберем материал кожуха алюминий, в качестве жидкости воду.
Во всех трех случаях температура нагрева кожуха - 100 градусов, а температура жидкости изначально равно 10 градусам. По рисункам можно определить, что площадь нагрева будет больше, конечно же, у канала с шестью ребрами. Температура понижается по направлению к середине канала.
6.3 Программа расчета и моделирования физических полей ELCUT
ELCUT® - это мощный современный комплекс программ для инженерного моделирования электромагнитных, тепловых и механических задач методом конечных элементов. Дружественный пользовательский интерфейс, простота описания даже самых сложных моделей, широкие аналитические возможности комплекса и высокая степень автоматизации всех операций позволяют разработчику полностью сосредоточиться на своей задаче.
Это приложение может быть использовано для проектирования и анализа теплового состояния различных систем. Можно вычислять как установившееся распределение температуры, так и изучать процессы нагрева и охлаждения.
ELCUT широко используется в научных исследованиях, промышленности и образовании. Академические и прикладные научные центры используют моделирование полей для изучения разнообразных физических процессов, а также для разработки приборов и устройств, применяемых в научных исследованиях.
ELCUT становится стандартным инструментом для обучения студентов различных специальностей и все более широко входит в практику преподавания физики, электротехники, прикладных дисциплин во множестве высших учебных заведений России и стран СНГ.
Ведущие электромашиностроительные заводы страны применяют ELCUT для расчета электромагнитных параметров и тепловых режимов электрических машин. Необходимость в анализе поля возникает при расчете для оптимизации существующих машин, либо при расчете машин нетрадиционных конструкций. Тепловые расчеты с помощью ELCUT успешно используют при проектировании систем обогрева с помощью нагреваемых кабелей, специализированных и комбинированных систем контроля температуры, вентилируемых фасадов, ограждающих конструкций в строительстве и архитектурном проектировании.
ELCUT - это программа расчета и моделирования физических полей, и все ее применения так или иначе связаны с наукой. Традиционные направления, где применяется моделирование полей, - физические исследования, электротехника, строительная механика, - в последнее время дополняются многочисленными примерами из области нанотехнологий, биологии, химии.
Моделируемые в ELCUT физические задачи, связанные с электромагнитными полями, электрическими токами, температурными полями, упругими напряжениями и деформациями - сегодня возникают в самых различных направлениях научных исследований. Большая часть применений, описанных на страницах Промышленность и Образование, тоже можно отнести к научным. Популярность программы в ведущих российских и зарубежных научных центрах доказывается ее внедрениями и публикациями, связанными с ELCUT.
6.4 Проверка корректности работы программы
Результаты, полученные в процессе вычислений в дипломной работе, отвечают условиям адекватности. Чтобы убедиться в корректности работы приложения, приведем результаты, вычисленные в программе Elcut.
Температурное поле вычисляется в программе Elcut при таких же граничных условиях и температурных потоках, как и в программе реализованное нами. Если сравнить рисунок 6.7, сделанный в программном комплексе Elcut, и рисунок 6.8, построенный с помощью нашего приложения, то значения температур совпадают.Для выполнения дипломной работы был построен теплообменник, который состоит из металлического кожуха и цилиндрической трубы, по которой течет жидкость.
Требуется исследовать, как выглядит температурное поле жидкости при изменении настраиваемых параметров.
При написании дипломной работы были использованы знания, полученные в процессе учебы.
При выполнении поставленной задачи была: - создана математическая модель распределения температуры жидкости с заданной вязкостью в цилиндрическом канале;
- разработан алгоритм реализации модели с помощью выбранных методов;
- разработана программа по алгоритму с помощью выбранного программного обеспечения;
- разработана цифровая модель изменения поля температуры в зависимости от вязкости жидкости, теплопроводности материала кожуха, конфигурации канала и граничных условий;
- проведена проверка адекватности модели поставленной физической задачи.
В данном проекте также приведены теоретические сведения, необходимые для решения поставленной задачи и рассмотрены результаты выполненной работы.
Все, поставленные ранее задачи, выполнены. Полученный программный продукт готов к использованию.
Список литературы
1. Патанкар С.В., Численное решение задач теплопроводности и конвективного теплообмена при течениях в каналах: пер. с англ. Е.В. Калабина; под ред Г.Г. Янькова. - М.: Издательство МЭИ, 2003. - 312с., ил
2. Чижиумов С.Д., Основы гидродинамики: учеб. пособие - Комсомольск-на-Амуре : ГОУВПО «КНАГТУ», 2007. - 106 с.
3. Патанкар С.В., Численные методы решения задач теплообмена и динамики жидкости: пер. с англ. - М. : Энергоатомиздат, 1984. - 152 с.