Содержание
- 2. Градиент: Градиент функций многих переменных f(x1,…x2) в точке P0 – элемент векторного пространства, компоненты которого суть
- 3. Набла-оператор – дифференциальный оператор, обозначаемый через и символически записываемый в виде «вектора» с компонентами ∂/∂х1,…,∂/∂xn; умножение
- 4. Векторные линии поля gradФ определяются уравнениями или Линии градиента, определяемые этими уравнениями, перпендикулярны к поверхностям уровня
- 5. Ротор Ротор (вращающий, вихрь) векторной функции точки есть векторная функция точки, определяемая формулой: Оператор Лапласа: Оператор
- 6. Внимание! Вывод уравнения движения основывается на законе сохранения и/или изменения импульса. Область , состоящая при всех
- 7. Случай одномерного движения Для одномерного потока газа в соответствии со вторым законом Ньютона уравнение движения может
- 8. Случай многомерного движения Выделим в потоке идеальной жидкости (несжимаемой и без вязкости) элементарную массу в форме
- 9. Случай многомерного движения Таким образом, сумма проекций на боковые грани равна: Проекцию объемной силы определим в
- 10. Случай многомерного движения В результате получаем три уравнения: или где . После деления на dt получаем:
- 11. Случай многомерного движения Замечание: Полная (субстанционная) производная равна сумме локальной и конвективной. При установившемся движении локальная
- 12. Внимание! Вывод уравнения неразрывности основывается на законе сохранения массы. Область , состоящая при всех t из
- 13. Случай одномерного движения Газ занимает весь предоставленный ему объем, что подтверждается и в действительности, поэтому изменение
- 14. Случай многомерного движения Построим внутри потока неподвижный объем, в форме параллелепипеда, и будем рассматривать его как
- 15. Случай многомерного движения Поскольку пренебрегая бесконечно малыми величинами высших порядков, получаем: Если за время dt масса
- 16. Случай многомерного движения В начальный момент времени t масса внутри параллелепипеда равна δmt=ρdxdydz, а в конечный
- 17. Случай многомерного движения В частном случае, при стационарном (установившемся) движении, плотность ρ, как и все остальные
- 18. Замечание: Возможны и другие формы записи этих уравнений, в том числе и для случая динамики переменной
- 19. Для тепловых процессов закон сохранения и превращения энергии выражается в виде первого начала термодинамики, которое для
- 20. Для идеального газа имеем: , что приводит к следующему выражению: Выделяя, объем V, ограниченный поверхностью F,
- 21. Это позволяет записать: Тогда в прямоугольных координатах, принимая гипотезу: с учетом , для количества теплоты Qv,
- 22. Внимание! Уравнение Фурье-Кирхгофа содержит давление P, скорость течения среды w и плотность ρ. Замечание: Для общего
- 23. Отсюда: где dx/dt=wx; dy/dt=wy; dz/dt=wz представляют собой проекции вектора скорости течения среды на соответствующие координаты. С
- 24. Внимание! Для того чтобы подчеркнуть непосредственную связь между производной dφ/dt с движущейся средой (субстанцией), часто её
- 25. Модель невязкого нетеплопроводного газа Поскольку энергия системы меняется за счет работы внешних сил и дополнительного потока
- 26. Внимание! Вывод уравнения распространения теплоты основывается на законе сохранения энергии. Область , состоящая при всех t
- 27. Внимание! Очевидно, что интегральные законы сохранения и балансовые уравнения равносильны, т.к. они выражают одни и те
- 28. Фундаментальное свойство эквивалентности тепловой и механической энергии выражается в виде: , где V=1/ρ – удельный объем.
- 29. Вывод: Даже при решении относительно простой задачи о движении идеальной, гомогенной среды в пространстве необходимо располагать
- 30. Замечание: Оба метода являются кинематическими, они математически связаны друг с другом и возможен переход от уравнений
- 31. Движение одной какой-либо частицы жидкости определяется системой трех уравнений: где t – время, а x, y,
- 32. Тогда, движение всей массы жидкости будет известным, если известна система: где начальные координаты a, b и
- 33. Справочно: Траектория – линия, по которой движется некоторая частица, в частности жидкости. Линия тока – кривая,
- 34. Пример 1. Пусть задана основная система уравнений движения некоторого потока: Определить характер движения и его кинематические
- 35. Пример 2. Составить общие уравнения движения и определить скорости и ускорения для случая вращения жидкости с
- 36. Направляющие косинусы скоростей: Замечание: Для данной частицы они переменные и зависят от угловой скорости ω и
- 37. Основные понятия и определения Динамическая система – совокупность взаимодействующих объектов, состояние которой изменяется во времени, а
- 38. Стационарное состояние (решение) – состояние системы, при котором значения некоторых существенных для него характеристик не изменяются
- 39. Для установившегося движения поле скоростей остается неизменным: Аналитическим условием установившегося движения будет: Замечание: Для того, чтобы
- 40. Внимание! Если основная система уравнений движения известна, то можно определить полную скорость и по направлению и
- 41. Сумма - представляет собой ускорение, связанное с неравномерностью движения при перемещении частицы жидкости в пространстве: Вывод:
- 42. При пространственном движении имеем: или более полно: Внимание! Чтобы получить уравнение линии тока в конечной форме
- 44. Скачать презентацию
Градиент:
Градиент функций многих переменных f(x1,…x2) в точке P0 – элемент векторного
Градиент функций многих переменных f(x1,…x2) в точке P0 – элемент векторного
Градиент (т.е. шагающий) скалярного поля U – векторное поле, компонентами которого служат частные производные ∂U/∂x, ∂U/∂y, ∂U/∂z; обозначается gradU и выражается с помощью набла-оператора как
Градиент скалярной функции точки есть векторная функция точки, определяемая формулой:
где V1 – область, содержащая точку , S1 - замкнутая поверхность, ограничивающая область V1, δ – наибольшее расстояние от точки до точек поверхности S1.
ВВЕДЕНИЕ
Некоторые математические понятия
Набла-оператор – дифференциальный оператор, обозначаемый через и символически записываемый в виде
Набла-оператор – дифференциальный оператор, обозначаемый через и символически записываемый в виде
Пример: В прямоугольных декартовых координатах линейный оператор (набла) определяется формулой:
его применение к скалярным и векторным функциям точки формально соответствует некоммутативной операции умножения на вектор с декартовыми координатами:
Замечание: В каждой точке, где существует вектор , его модуль
равен наибольшей производной по направлению dФ/ds в этой точке, а его направление совпадает с направлением, в котором эта наибольшая производная достигается.
Некоторые математические понятия
Векторные линии поля gradФ определяются уравнениями или
Линии градиента, определяемые этими уравнениями,
Векторные линии поля gradФ определяются уравнениями или
Линии градиента, определяемые этими уравнениями,
Дивергенция:
Дивергенция векторного поля – скалярная характеристика поля V, обозначаемая через div V; в прямоугольных координатах дивергенция имеет вид ∂Vx/∂х + ∂Vy/∂y + ∂Vz/∂z, где Vx , Vy , Vz – компоненты вектора.
Дивергенция (расхождение) векторной функции точки есть скалярная функция точки, определяемая формулой:
Некоторые математические понятия
Ротор
Ротор (вращающий, вихрь) векторной функции точки есть векторная функция точки,
Ротор
Ротор (вращающий, вихрь) векторной функции точки есть векторная функция точки,
Оператор Лапласа:
Оператор Лапласа – дифференциальный оператор, обозначаемый ∆ и задаваемый формулой: ∂2f/∂х21 +…+ ∂2f/∂x2n , где f – функция переменных х1,…,хn.
∆= div grad =
Некоторые математические понятия
Внимание! Вывод уравнения движения основывается на законе сохранения и/или изменения импульса.
Внимание! Вывод уравнения движения основывается на законе сохранения и/или изменения импульса.
Область , состоящая при всех t из одних и тех же частиц, называется движущимся объемом, имеющим поверхность γ(t) с единичным вектором внешней нормали
Если импульс – , то закон сохранения и/или изменения
импульса может быть записан в виде:
Замечание: В подобной постановке силами, действующими на объем ω(t) будут только поверхностные силы давления, направленные по нормали к поверхности рассматриваемого объекта.
Случай одномерного движения
Движение газа происходит в одномерном поле давления p=f(x) с градиентом – dp/dx, поэтому на выделенный элементарный объем dV=Fdx должна действовать сила , сообщающая массе газа ρFdx= ρdV в стационарном потоке
ускорение, равное:
Замечание: Такое ускорение приобретает элементарная масса газа, движущаяся со скоростью w в поле скорости w=f(x) с градиентом dw/dx.
Семинар 1. Уравнения движения невязкой среды (уравнения Эйлера)
Случай одномерного движения
Для одномерного потока газа в соответствии со вторым законом
Случай одномерного движения
Для одномерного потока газа в соответствии со вторым законом
откуда или
и
С учетом соотношения ρ=1/v последнее уравнение получает вид:
где v – удельный объем.
Вывод: Полученное уравнение показывает, что dp и dw имеют всегда разные алгебраические знаки. Это свидетельствует о том, что скорость одномерного потока газа возрастает только в направлении уменьшения давления.
Случай многомерного движения
Выделим в потоке идеальной жидкости (несжимаемой и без вязкости)
Случай многомерного движения
Выделим в потоке идеальной жидкости (несжимаемой и без вязкости)
; ;
Тогда, для первого уравнения имеем:
; ;
Определим силу dFx, полагая, что внешними являются поверхностные силы (силы давления жидкости, окружающей выделенную элементарную массу, на ее боковые грани) и силы объемные (или массовые).
Замечание: Проекции сил давления на боковые грани, перпендикулярные оси ОХ, равны: - по оси (знак плюс),
- против оси (знак минус), где p и p’ - средние гидродинамические давления для граней:
Случай многомерного движения
Таким образом, сумма проекций на боковые грани равна:
Проекцию объемной
Случай многомерного движения
Таким образом, сумма проекций на боковые грани равна:
Проекцию объемной
где X - проекция ускорения объемной силы.
Тогда проекция равнодействующей:
и
и после отнесения к единице массы имеем:
Случай многомерного движения
В результате получаем три уравнения:
или
где . После деления на
Случай многомерного движения
В результате получаем три уравнения:
или
где . После деления на
;
Вывод: Полная производная dwx/dt=ax представляет собой полное ускорение движения данной частицы в направлении оси ОХ, она обычно направлена по скорости (но можно и произвольно) и иногда называется – субстанционной (индивидуальной, материальной и т.п.) производной по времени.
Частная производная ∂wx/∂t, представляет собой ускорение в данной точке, обусловленное неустановившимся движением жидкости ее часто называют локальной производной.
Сумма wx∂wx/∂x+wy∂wx/∂y+wz∂wx/∂z представляет собой ускорение, связанное с неравномерностью движения при перемещении частицы жидкости в пространстве – конвективной производной.
Случай многомерного движения
Замечание: Полная (субстанционная) производная равна сумме локальной и конвективной.
Случай многомерного движения
Замечание: Полная (субстанционная) производная равна сумме локальной и конвективной.
Краткая форма записи уравнения движения
При постоянной массе: При переменной массе:
где Р – тензор напряжений; - плотность распределения реактивных сил, т.е. реактивные силы, отнесенные к единице объема.
Замечание: Возможны два подхода к получению исходных уравнений интегральных законов сохранения. Первый из которых рассматривает изменение во времени массы в переменном объеме ω(t), а второй – в фиксированном (не зависящем от времени) объеме ω, что дает возможность получения законов сохранения в виде балансовых уравнений.
Первый подход: Второй подход:
Вывод: Полученные уравнения относятся к движению как капельной, так и газообразной среды. В систему трех уравнений входят пять неизвестных: составляющие скорости – wx, wy, wz, давление - p и плотность ρ. Поэтому потребуется дополнительно еще два уравнения.
Внимание! Вывод уравнения неразрывности основывается на законе сохранения массы.
Область ,
Внимание! Вывод уравнения неразрывности основывается на законе сохранения массы.
Область ,
Если масса – , то закон сохранения массы может быть записан в виде:
Случай одномерного движения
Скорость прохождения массы газа (массовая скорость) wm через единицу площади поперечного сечения F определяется произведением ρw, где ρ – плотность газа, w – его скорость.
Если учесть, что по мере перемещения газа в канале изменение массовой скорости составляет ∂(ρw)/∂х, то элементарный расход массы dM газа за время dt можно записать в виде соотношения:
Замечание: Изменение F на длине dx при этом не учитывается.
Семинар 2. Уравнение неразрывности
Основные уравнения тепло- и массообмена в биосфере
Случай одномерного движения
Газ занимает весь предоставленный ему объем, что подтверждается и
Случай одномерного движения
Газ занимает весь предоставленный ему объем, что подтверждается и
или
Вывод: Полученное уравнение и является уравнением неразрывности.
Для одномерного потока частную производную ∂(ρw)/∂х можно заменить на полную. Тогда изменение массовой скорости может быть записано в виде: d(ρw)/dх.
Для стационарного потока ρw=const и ρw=М/F.
В рассматриваемых условиях уравнение неразрывности получает вид:
М=ρwF=const
Последнее уравнение можно представить в дифференциальной форме:
или
где v – удельный объем.
Случай многомерного движения
Построим внутри потока неподвижный объем, в форме параллелепипеда, и
Случай многомерного движения
Построим внутри потока неподвижный объем, в форме параллелепипеда, и
Если поток не претерпевает разрыва, то для несжимаемой среды количество массы жидкости, входящей в объем и выходящей из него должны быть равны между собой.
Для сжимаемой среды разность между этими количествами должна равняться приращения массы внутри объема, обусловленному изменением плотности.
В направлении ОХ, через боковые грани, перпендикулярные этой оси, за время внутрь параллелепипеда поступит масса:
а выйдет , где плотность и скорость на выходе
могут быть определены в виде: и
Поскольку втекание и вытекание происходит одновременно, их изменение обуславливается только изменением координаты х.
Следовательно, выходящая масса равна:
Случай многомерного движения
Поскольку
пренебрегая бесконечно малыми величинами высших порядков,
получаем:
Если
Случай многомерного движения
Поскольку
пренебрегая бесконечно малыми величинами высших порядков,
получаем:
Если
Аналогично получаем для других осей:
и полное приращение равно:
Замечание: Это изменение массы при условии неразрывности движения должно равняться изменению массы, обусловленному изменением плотности.
Случай многомерного движения
В начальный момент времени t масса внутри параллелепипеда равна
Случай многомерного движения
В начальный момент времени t масса внутри параллелепипеда равна
Поэтому: ,
и в конечный момент времени масса жидкости в
объеме составит:
Таким образом, приращение массы за время dt равно:
но из условия неразрывности: δm= δM, т.е.
После сокращения окончательно получаем:
Вывод: Это и есть искомое уравнение неразрывности.
Случай многомерного движения
В частном случае, при стационарном (установившемся) движении, плотность ρ,
Случай многомерного движения
В частном случае, при стационарном (установившемся) движении, плотность ρ,
Для несжимаемой жидкости ρ=const и как при установившемся, так и при неустановившемся движении, уравнение неразрывности записывается в виде:
Внимание!!! Поскольку по физическому смыслу частные производные, которые входят в последнее уравнение, представляют собой соответствующие скорости линейных деформаций растяжения или сжатия в направлении осей, это уравнение неразрывности показывает, что в условиях сплошности движения несжимаемой жидкости сумма скоростей относительных линейных деформаций ребер параллелепипеда равна нулю, т.е. коэффициент объемного расширения равен нулю.
Другими словами, это уравнение показывает, что движение жидкости осуществляется таким образом, что данная масса все время занимает один и тот же объем.
Замечание: Возможны и другие формы записи этих уравнений, в том числе
Замечание: Возможны и другие формы записи этих уравнений, в том числе
Для несжимаемой
жидкости:
Краткие формы записи уравнения неразрывности
Для сжимаемой жидкости
(газа):
Вывод: Здесь также возможны два подхода к получению исходных уравнений интегральных законов сохранения. Первый из которых рассматривает изменение во времени массы в переменном объеме ω(t), а второй – в фиксированном (не зависящем от времени) объеме ω, что дает возможность получения законов сохранения в виде балансовых уравнений.
Первый подход Второй подход
Здесь γ(t) – поверхность ω(t) и – единичный вектор внешней нормали к γ(t)
Для постоянной массы:
Для переменной массы:
Для тепловых процессов закон сохранения и превращения энергии выражается в виде
Для тепловых процессов закон сохранения и превращения энергии выражается в виде
где Qv – количества теплоты, поступающей в единицу объема среды в единицу времени, Вт/м3; Lv – работа, совершаемая внешними силами над единицей объема среды за единицу времени, Вт/м3; U – внутренняя энергия одного килограмма среды, Дж/кг; w – скорость движения (течения) среды, м/с.
Замечание: Это уравнение выражает то обстоятельство, что изменение полной энергии тела, которое складывается в данном случае из изменения его внутренней энергии ρU и изменения кинетической энергии ρw2/2, обусловлено количеством теплоты, подводимой к телу, и внешней работой, совершенной над телом.
Поскольку рассматривается движение элемента потока жидкости (газа), то работа в данном случае связана с изменением давления и внутренним трением текущей среды, а изменение внутренней энергии связанно с изменением ее теплосодержания, т.е. энтальпии, уравнением:
Семинар 3. Уравнения распространения теплоты в вещественной среде
Для идеального газа имеем: , что приводит к следующему выражению:
Выделяя, объем
Для идеального газа имеем: , что приводит к следующему выражению:
Выделяя, объем
Первый член этого уравнения представляет собой изменение теплосодержания рассматриваемого объема V; второй – количество теплоты, ушедшее через поверхность F путем теплопроводности; третий – количество теплоты, выделенное внутренними источниками (имеет знак минус, если в теле находятся стоки теплоты).
Замечание: Внутренние источники теплоты могут возникать вследствие излучения, объемных химических реакций, радиоактивного распада вещества, прохождения электрического тока, работы трения и т.д.
Между потоком вектора через замкнутую поверхность F, ограничивающую объем V, и расходимостью (диверегенцией) вектора существует связь, выражаемая формулой Остроградского-Гаусса:
Замечание: Формула справедлива, если в объеме V нет сильных разрывов функции. В тепловой задаче это означает, что область V не должна включать границы раздела фаз.
Уравнения распространения теплоты в вещественной среде
Это позволяет записать:
Тогда в прямоугольных координатах, принимая гипотезу: с учетом ,
Это позволяет записать:
Тогда в прямоугольных координатах, принимая гипотезу: с учетом ,
где λ - коэффициент теплопроводности, qv - интенсивность выделения теплоты внутренним источником.
В результате получаем известное уравнение Фурье-Кирхгофа:
Вывод: Это уравнение параболического типа. Известно, что оно имеет нестационарные решения при допущении о бесконечно большой скорости распространения тепловых возмущений.
Замечание: Реальная скорость распространения тепловых возмущений в вещественной среде имеет порядок не более средней скорости теплового движения структурных частиц, а для излучения равна скорости распространения электромагнитных волн. Однако, во многих практических приложениях, подобное допущение находится в пределах точности инженерных расчетов и использование нестационарных решений уравнения Фурье-Кирхгофа следует признать допустимым.
Уравнения распространения теплоты в вещественной среде
Внимание! Уравнение Фурье-Кирхгофа содержит давление P, скорость течения среды w и
Внимание! Уравнение Фурье-Кирхгофа содержит давление P, скорость течения среды w и
Замечание: Для общего решения задачи о теплообмене в движущейся вещественной среде к этому уравнению необходимо присоединить ещё уравнения, определяющие поле скоростей и связь между термодинамическими параметрами состояния среды.
Вывод: Замыкание системы дифференциальных уравнений теплообмена в движущейся вещественной среде, достигается присоединением к уравнению распространения тепла уравнений движения и неразрывности потока жидкости и уравнения состояния. Тем самым вновь прослеживается связь между гидродинамикой и тепло- и массообменом.
Справочно: Проводя аналогии между этими процессами можно констатировать, что полные производные по времени от давления, теплосодержания (энтальпии), квадрата скорости и плотности являются функциями координат и времени, в связи с чем полный дифференциал любой из них равен:
Уравнения распространения теплоты в вещественной среде
Отсюда:
где dx/dt=wx; dy/dt=wy; dz/dt=wz представляют собой проекции вектора скорости течения среды
Отсюда:
где dx/dt=wx; dy/dt=wy; dz/dt=wz представляют собой проекции вектора скорости течения среды
С физической точки зрения, рассматриваемая полная производная, распадается на две разные части. Первый член правой части характеризует изменение данной величины, связанное с изменением поля этой величины во времени, а сумма остальных трех членов характеризует изменение данной величины, происходящее в связи с перемещением рассматриваемого элемента среды из одной точки пространства в другую.
Иными словами, частная производная ∂φ/∂t представляет собой скорость изменения φ во времени в той точке пространства, в которой элементарный объем среды находится в данный момент времени, а сумма:
представляет собой скорость изменения, которое обусловлено перемещением элементарного объема среды dV из точки со значением φ=φ1 в точку со значением φ=φ2. При этом величина ∂φ/∂t – обычно называется местным (локальным) изменением, а величина - изменением перемещения или конвективным изменением.
Уравнения распространения теплоты в вещественной среде
Внимание! Для того чтобы подчеркнуть непосредственную связь между производной dφ/dt с
Внимание! Для того чтобы подчеркнуть непосредственную связь между производной dφ/dt с
т.к.
Вводя значение субстанционной производной в уравнение Фурье-Кирхгофа, в векторной форме получаем:
В прямоугольной системе координат уравнение Фурье-Кирхгофа имеет вид:
Вывод: Решение системы уравнений, включающей в себя уравнения движения, неразрывности и распространения теплоты не является тривиальным, поэтому практика инженерных расчетов обуславливает необходимость введения обоснованных упрощающих допущений.
Уравнения распространения теплоты в вещественной среде
Модель невязкого нетеплопроводного газа
Поскольку энергия системы меняется за счет работы
Модель невязкого нетеплопроводного газа
Поскольку энергия системы меняется за счет работы
где ω(t) – объем, ε – удельная внутренняя энергия.
В рамках данной модели, газ движется в отсутствии внешних силовых полей и внешних источников энергии. При этом, силами, действующими на объем ω(t) будут только поверхностные силы давления, направленные по нормали к поверхности этого объема. Это позволяет записать:
где γ(t) поверхность ω(t), n – единичный вектор внешней нормали к γ(t).
Кроме того, в этой модели Q=0, а W есть мощность, развиваемая силами давления на γ(t), т.е.
Семинар 4. Частные случаи закона сохранения энергии и
уравнений теплопроводности
Внимание! Вывод уравнения распространения теплоты основывается на законе сохранения энергии.
Область
Внимание! Вывод уравнения распространения теплоты основывается на законе сохранения энергии.
Область
Если энергия – , то закон сохранения энергии может
быть записан в виде:
где полная энергия равна сумме кинетической и внутренней энергии.
Интегральные законы сохранения массы, импульса и энергии в рассматриваемой модели имеют вид:
Соответственно, балансовые уравнения интегральных законов сохранения:
Замечание: Введенное здесь обозначение оператора дифференцирования ∂/∂t (вместо d/dt) подчеркивает, что ввиду независимости ω от t этот оператор может быть внесен под знак интеграла как оператор частной производной по t.
Частные случаи законов сохранения
Внимание! Очевидно, что интегральные законы сохранения и балансовые уравнения равносильны, т.к.
Внимание! Очевидно, что интегральные законы сохранения и балансовые уравнения равносильны, т.к.
Вывод: Представленные системы уравнений являются неопределенными и для их дополнения требуется привлечение термодинамических и характеристических свойств среды.
Термодинамические свойства
При термодинамическом рассмотрении статистически равновесных процессов в атмосфере, наряду с введенными выше параметрами состояния: ρ, р, ε обычно используют еще два основных параметра состояния: температуру T и удельную (отнесенную к единице массы) энтропию S.
Вывод: Фактически при рассмотрении атмосферных явлений, обычно допускают, что атмосферный воздух как термодинамическая система является двухпараметрической средой. Это означает, что его состояние вполне определяется заданием каких-либо двух параметров. Следовательно упомянутые пять параметров должны быть связана тремя соотношениями.
Модель невязкого нетеплопроводного газа
Фундаментальное свойство эквивалентности тепловой и механической энергии выражается в виде: ,
Фундаментальное свойство эквивалентности тепловой и механической энергии выражается в виде: ,
Левая часть данного равенства равна количеству тепла, сообщенному единице массы, а правая – приращению внутренней энергии плюс работа расширения этой порции газа.
Замечание: Для двухпараметрических сред данное соотношение является определением энтропии.
Если рассматривать внутреннюю энергию как функцию удельного объема и энтропии, т.е. положить ε = E(V,S), то данное равенство эквивалентно двум соотношениям: p = Ev(V,S), T = Es(V,S).
Вывод: Задание функции E(V,S) полностью описывает всю термодинамику двухпараметрической среды. Однако на практике такое задание не всегда удобно и термодинамические свойства среды описываются другими соотношениями, в том числе характеристическими, основными из которых являются соотношения:
Идеального газа: pV=RT;
Политропного газа: p=A(S)ρn;
Нормального газа: p=f(ρ,S), ε= e(V,p);
Нормальной жидкости: ρ =φ(p,t);
Несжимаемой изотермической жидкости (газа): ρ=const; T=const.
Первый закон термодинамики
Модель невязкого нетеплопроводного газа
Вывод: Даже при решении относительно простой задачи о движении идеальной, гомогенной
Вывод: Даже при решении относительно простой задачи о движении идеальной, гомогенной
Замечание: В простейшем случае, одной пространной переменной, без введения дополнительных упрощающих допущений, количество таких уравнений может быть сокращено только до трех.
Внимание! Применительно к процессам в гидросфере, вместо уравнений Эйлера зачастую используются уравнения движения вязкой жидкости (уравнения Навье-Стокса):
Вывод: В практике инженерных расчетов для вязкой жидкости в пространстве, необходимо введение в систему либо дополнительных уравнений, либо обоснованных упрощающих допущений.
Справочно: Для реальных процессов тепломассообмена в биосфере, даже без учета дисперсности системы, количество необходимых уравнений может достигать 17 и более (обычно 14).
Модель невязкого нетеплопроводного газа
Замечание: Оба метода являются кинематическими, они математически связаны друг с другом
Замечание: Оба метода являются кинематическими, они математически связаны друг с другом
Основные различия методов:
По методу Лагранжа предусматривается изучение законов движения каждой индивидуальной частицы.
В методе Эйлера задача заключается в изучении поля скоростей, ускорений и других параметров движения. При этом вопрос о том, как движется та или иная индивидуальная частица остается в стороне.
Основное допущение методов - жидкость (газ) рассматривается как легкодеформируемая непрерывная среда. В качестве мельчайшего структурного элемента среды принимается "частица" бесконечно малых размеров.
Семинар 5. Основные уравнения кинематики и динамики
невязкой жидкости
Два метода исследования движения жидкости –
метод Лагранжа и метод Эйлера
Движение одной какой-либо частицы жидкости определяется системой трех уравнений:
где t –
Движение одной какой-либо частицы жидкости определяется системой трех уравнений:
где t –
Внимание! Для описания движения всех частиц потока потребовалось бы бесконечное множество таких уравнений, что невозможно.
Замечание: Это затруднение можно устранить следующим образом: пусть в начальный момент времени t0 (принятый за начало отсчета) каждая частица находится в определенной точке пространства, положение которой определяется координатами x0=a, y0=b, z0=c, где a, b, c - начальные координаты.
Вывод: В результате выбор численных значений a, b и c является выбором конкретной частицы. Её текущие координаты будут другими, нежели текущие координаты частицы с другими начальными координатами. Поэтому для потока в целом текущие координаты зависят не только от времени, но и начальных координат.
Метод Лагранжа
Тогда, движение всей массы жидкости будет известным, если известна система:
где
Тогда, движение всей массы жидкости будет известным, если известна система:
где
Замечание: Если система известна (задана или найдена), то движение определяется в полной мере.
Здесь, составляющие скорости w любой частицы определяются как первые производные по времени от x, y, z и записываются в виде частных производных системы:
ускорения как вторые:
полная скорость и ускорение в виде:
а направление w будет определено по косинусам углов α, β, γ.
Вывод: При исследовании движения по методу Лагранжа геометрические характеристики движения являются траектории и линии отмеченных точек.
Метод Лагранжа
Справочно: Траектория – линия, по которой движется некоторая частица, в частности
Справочно: Траектория – линия, по которой движется некоторая частица, в частности
Линия тока – кривая, проходящая через такие частицы, скорости которых в данный момент времени направлены по касательной к этой линии. Линии тока не пересекаются между собой.
Замечание: При неустановившемся движении линии тока и траектории не совпадают друг с другом. Они также могут не совпадать в гетерогенных системах.
Линия отмеченных точек – линия, на которой в данный момент времени лежат частицы жидкости прошедшие в свое время через одну и ту же начальную точку.
Замечание: При установившемся движении все эти три линии совпадают между собой.
Вывод: В методе Лагранжа траектория выбранной частицы N по выбранным (aN, bN, cN) определяется либо непосредственно исходной системой, путем вычисления x, y, z для ряда моментов времени t1,…t2, либо путем исключения из системы времени t, что приводит к двум уравнениям:
совместное решение которых дает траекторию.
Метод Лагранжа
Пример 1. Пусть задана основная система уравнений движения некоторого потока:
Определить характер
Пример 1. Пусть задана основная система уравнений движения некоторого потока:
Определить характер
Зная положение начальной частицы M(t0), определим для частицы М(t) величину и направление скорости и ускорения.
Проекции полной скорости:
Полную скорость:
Направление движения:
Cоставляющие ускорения:
Вывод: Движение потока является равномерным и прямолинейным.
Замечание: Этот простейший пример иллюстрирует последовательность решения задачи.
Метод Лагранжа
Пример 2. Составить общие уравнения движения и определить скорости и ускорения
Пример 2. Составить общие уравнения движения и определить скорости и ускорения
Выберем произвольную точку М(t0) с координатами x0=a, y0=b, z0=c, таким образом, чтобы ее радиус-вектор равный , был направлен под углом α к оси 0х.
За время t точка М переместилась в положение М(t) и радиус-вектор повернулся на угол α’= ωt. Тогда определим:
Координаты точки:
Замечание: Полученная система трех уравнений - основная система уравнений движения в форме Лагранжа.
Проекции полной скорости:
Полную скорость:
Метод Лагранжа
Направляющие косинусы скоростей:
Замечание: Для данной частицы они переменные и зависят от
Направляющие косинусы скоростей:
Замечание: Для данной частицы они переменные и зависят от
Составляющие ускорения:
Полное ускорение:
Направляющие косинусы ускорений:
Замечание: Для данной частицы они переменные и зависят от угловой скорости ω и времени t.
Вывод: Таким образом, в данном случае полная скорость частицы w=ωr=f(r) является функцией радиуса, а полное ускорение a направлено по радиусу r к центру, т.е. является центростремительным ускорением.
Метод Лагранжа
Основные понятия и определения
Динамическая система – совокупность взаимодействующих объектов, состояние которой
Основные понятия и определения
Динамическая система – совокупность взаимодействующих объектов, состояние которой
Состояние системы - состояние, определяемое совокупностью значений, характерных для данной системы величин, которые обычно называют параметрами состояния.
Переходный процесс – изменение во времени координат динамической системы, возникающее под влиянием воздействий, резко изменяющих её состояние, структуру и параметры.
Нестационарное движение – движение жидкости, при котором в каждой точке потока скорость жидкости, её давление и другие характеристики изменяются во времени.
Стационарное движение – движение жидкости, при котором в каждой точке потока скорость жидкости, её давление и другие характеристики не изменяются во времени.
Замечание: Вследствие не нулевых начальных условий переходный процесс часто представляет собой переход от одного установившегося и/или равновесного состояния системы к другому установившемуся и/или равновесному состоянию.
Семинар 6. Метод Эйлера
Стационарное состояние (решение) – состояние системы, при котором значения некоторых существенных
Стационарное состояние (решение) – состояние системы, при котором значения некоторых существенных
Установившиеся процессы = стационарные процессы – процессы, протекающие в системе, не изменяющей при этом своего состояния.
Установившееся движение = установившейся режим – состояние динамической системы после окончания переходного процесса, когда система находится в равновесии, совершает вынужденные колебания, автоколебания и т.п.
Режим – условие или область существования процесса.
Автомодельный режим – режим, не зависимый от какого-либо параметра.
Метод Эйлера
Рассмотрим движение жидкости в некоторой области. В каждой точки этой области в данный момент времени частицы жидкости имеют скорость w.
Внимание! Вся совокупность векторов, изображающих соответствующие скорости составляет – векторное поле скоростей.
В общем случае (при неустановившемся движении) всё поле изменяется во времени и компоненты скорости являются функциями не только координат, но и времени. Тогда, основная система уравнений имеет вид:
Для установившегося движения поле скоростей остается неизменным:
Аналитическим условием установившегося движения будет:
Для установившегося движения поле скоростей остается неизменным:
Аналитическим условием установившегося движения будет:
Замечание: Для того, чтобы полная скорость оставалась одной и той же как по направлению, так и по значению надо соблюдать аналитическое условие, лишь в этом случае полная скорость остается одной и той же, как по направлению, так и по значению.
Вывод: Равенства нулю частной производной полной скорости по времени ∂w/∂t=0, ещё недостаточно.
Внимание! Строго говоря, для того чтобы движение было установившимся, т.к. частная производная ∂w/∂t определяет лишь изменение скорости во времени и если ∂w/∂t=0, то это обозначает лишь, что w не изменяется по величине.
Однако при этом направление скорости может и изменяться, т.е. движение может быть и неустановившимся, но стационарным.
Метод Эйлера
Внимание! Если основная система уравнений движения известна, то можно определить полную
Внимание! Если основная система уравнений движения известна, то можно определить полную
При этом можно найти полное ускорение:
где проекции ускорения по координатным осям будут равны:
Составим выражение полных производных dwx/dt, dwy/dt, dwz/dt:
После деления на dt получаем:
Полная производная dwx/dt=ax представляет собой полное ускорение движения данной частицы в направлении оси 0x. Она обычно направлена по скорости (но можно и произвольно) и иногда называется – субстанционной (индивидуальной, материальной и т.п.) производной по времени.
Частная производная ∂wх/∂t, представляет собой ускорение в данной точке, обусловленное неустановившимся движением жидкости ее часто называют локальной производной.
Метод Эйлера
Сумма - представляет собой ускорение,
связанное с неравномерностью движения при перемещении
Сумма - представляет собой ускорение,
связанное с неравномерностью движения при перемещении
Вывод: Полная, т.е. субстанционная производная равна сумме локальной и конвективной.
Внимание! При установившемся движении локальная производная ∂wх/∂t=0, но конвективная может и не быть равной нулю – случай неравномерного установившегося движения.
Аналогично для других производных dwy/dt, dwz/dt.
Замечание: При исследовании движения жидкости по методу Эйлера геометрическими характеристиками движения являются линии тока, а не траектории, как в методе Лагранжа.
Пусть в данный момент времени имеем плоскую линию тока ab, её уравнение y=f(x), где x, y - координаты точек данной линии. Тогда wy/wx=tgα, где α - угол касательной к линии тока по отношению к оси 0x. Поэтому:
и или
Вывод: Это и есть дифференциальное уравнение линии тока на плоскости.
Метод Эйлера
При пространственном движении имеем:
или более полно:
Внимание! Чтобы получить уравнение линии
При пространственном движении имеем:
или более полно:
Внимание! Чтобы получить уравнение линии
Принимая время t как параметр и интегрируя, получаем два уравнения:
Вывод: Совместное решение этих уравнений дает уравнение линии тока в конечной форме.
Замечание: Для установившегося движения линии тока совпадают с траекторией, а потому полученные уравнения будут одновременно и уравнениями траекторий.
Внимание! Переход от координат Лагранжа к координатам Эйлера производится дифференцированием, что всегда возможно. Обратный переход должен осуществляться путем интегрирования, что всегда сложнее, а иногда и невозможно.
Метод Эйлера