Исследование движения жидкости и механических примесей в ней методом конечных элементов
Движение жидкости будет турбулентным. Турбулентные движения характеризуются тем, что поле истинных скоростей частиц жидкости, которая рассматривается как сплошная среда-континуум, имеет нерегулярный пульсационный характер, является неустановившимся и напоминает хаотическое поле скоростей отдельных молекул, из которых состоят тела. Траектории частиц жидкости при турбулентном движении в высшей степени извилисты. Точные измерения указывают на то, что изменения всех параметров течения при турбулентном движении во времени имеют вид, изображенный на рисунке 2.1, где для примера дано характерное для турбулентных течений изменение плотности в данной точке пространства с течением времени t.
Рис 2.2 Плотность турбулентного течения
На основное изменение плотности (на рисунке- плавная пунктирная линия) накладываются нерегулярные пульсации большой частоты. Во многих вопросах турбулентные течения жидкости целесообразно изучать только в среднем. При исследовании турбулентных течений обычно вводят средние значения компонент скорости Таким образом, в случае турбулентных течений сложное движение жидкости осредняется и при этом возникают проблемы составления полной системы уравнений для определения средних характеристик движения и проблемы изыскания способов экспериментального измерения осредненных характеристик движения. Практика построения моделей для изучения турбулентных движений показывает, что способы введения средних характеристик движения, вообще говоря, несущественны для составления полной системы уравнений теорий турбулентности, но они являются главной основой для разработки методов экспериментальных измерений различных средних величин, проведение которых необходимо для сравнения результатов предложенной теории турбулентности с опытными данными.
Укажем некоторые возможные способы осреднения истинных характеристик движения. Пусть А(х,у,z,t) - некоторая истинная характеристика турбулентного движения. В любой фиксированной точке пространства можно провести осреднение А по времени t. Тогда среднее значение А будет равно
где промежуток времени Т достаточно велик по отношению ко времени отдельных пульсации и мал по отношению ко времени заметного изменения средних характеристик (осредненное движение может быть нестационарным). С другой стороны, в определенный момент времени t осреднение А можно провести по объему; тогда
причем объем V должен удовлетворять условиям, аналогичным условиям, наложенным на промежуток времени Т. Можно провести осреднение по времени и по объему V одновременно. Указанные выше осреднения по времени и по объему возможно проводить с весом, когда среднее значение А определяется, например, следующим образом:
где g(t) - некоторая заданная функция. В различных задачах при выборе V и Т можно руководствоваться различными соображениями, но в имеющихся приложениях результат осреднения рассматривается как не зависящий от V и T. В ряде случаев используются вероятностные способы осреднения, и среднее значение A часто определяется как математическое ожидание A. После введения среднего значения A истинное значение A представляется в виде
A=
где A'-пульсация A; среднее значение пульсации равно нулю,
Потребуем, чтобы операции осреднения во всех случаях обладали следующими свойствами: 1 Среднее значение суммы равняется сумме средних значений:
2 Среднее значение производной от истинной характеристики турбулентного движения равняется производной от среднего значения:
Среднее значение произведения двух сомножителей, из которых только один испытывает турбулентные пульсации, равно произведению средних. В частности,
а равняется сумме произведения средних величин и среднего значения произведений пульсации этих величин:
Заметим, что при определении средних значений с помощью интегрирования по времени или по пространству перечисленные здесь свойства осреднения выполняются лишь приближенно. Пусть проведено осреднение скорости u=
Можно назвать средним значением u2 макроскопическую величину
Для истинных значений u2 в этом случае можно написать
причем Если в механике жидкости для величин ввести средние одинаковым способом, то характеристики движения таким образом осредненного континуума не будут удовлетворять основным законам сохранения и уравнениям состояния, удовлетворяющимся для истинных движений. В различных теориях турбулентности для определенного набора основных величин осреднения вводятся некоторым одинаковым способом, а способы осреднения других величин вводятся по соглашению так, чтобы удовлетворялись основные законы физики, как и при обычном определении этих величин для истинных движений. Рассмотрим турбулентные движения несжимаемой вязкой жидкости. Полная система уравнении движения в этом случае, как известно, состоит из уравнения неразрывности и уравнений импульса, которые в декартовой системе координат имеют вид
где
Для дальнейшего мы не будем фиксировать закон зависимости
При многих теоретических исследованиях турбулентных движений исходят из предпосылки о справедливости уравнений (2.7),(2.8) для истинного неустановившегося пульсирующего движения. Однако ввиду крайней запутанности, извилистости и сложности траекторий частиц жидкости при турбулентном движении получение решений этих уравнений для турбулентных движений представляет собой громоздкую и сложную задачу. Ставится задача о разыскании функциональных соотношений между средними величинами. Уравнения движения для средних, величин получаются путем осреднения уравнений движения (2.7),(2.8) для величин, описывающих мгновенное состояние движения. Осреднив уравнение неразрывности (2.7), на основании свойств операции осреднения легко получим уравнение неразрывности для осредненных величин
оно имеет тот же вид, что и для истинных скоростей (2.7). Осредним уравнения импульса, левую часть которых предварительно запишем в виде (2.10). Так как
(плотность
которые называются уравнениями Рейнольдса. Если связь
Шесть различных величин
вошедших в уравнения Рейнольдса, называются турбулентными напряжениями. Вид зависимости турбулентных напряжений от средних характеристик течения в различных классах задач может быть различным. Таким образом, ввиду нелинейности уравнений истинных движений, после их осреднения получается большее, чем число уравнений, число неизвестных. Следовательно, для математического изучения осредненных турбулентных движений одних уравнений гидромеханики, достаточных для изучения истинных движений, недостаточно. Поэтому полное теоретическое исследование осредненных турбулентных движений возможно только на основании некоторых дополнительных законов или гипотез, справедливость которых может быть в конечном счете установлена только на опыте.
Содержание многих работ по исследованию турбулентных движений [2,10,11,12,13] сводится к изучению справедливости различных простых и естественных гипотез о зависимости турбулентных напряжений от средних скоростей и их градиентов, которые позволяют поставить и решить теоретически основные частные задачи о турбулентном движении. В настоящее время не существует общей математической постановки задачи о произвольных осредненных турбулентных движениях и вообще не выяснена возможность такой формулировки задачи. Иногда по аналогии с законом Навье - Стокса (2.9) полагают, что
где M1=
где с=[0,1:0,2], Dx- шаг сетки (при вычислении методом конечных элементов). Найдем
т.к. по свойствам осредненной функции скорости и учитывая уравнение неразрывности (2.11)
подставим (2.17) и (2.18) в (2.16), получим
подставим (2.19) в (2.13), получим замкнутую систему уравнений
Осталось только добавить начальные и граничные условия
Решить эту систему в общем случае аналитически нельзя, она решается численными методами, например методом конечных элементов. Вычисления проводились в системе ANSYS, картину распределения поля скоростей и давления можно увидеть на рис 2.3 и рис 2.4.
Рис 2.3 Поле распределения скоростей в баке (v0=2,358 м/с)
Рис 2.4 Поле распределения давления в баке Рис 2.5 Поле распределения скоростей в трубопроводе (v0=3,537 м/с)
Начальное значение скорости вычисляется исходя из расхода жидкости
где Q - расход жидкости (м3/с), d - диаметр трубопровода (м). Теперь рассмотрим турбулентное движение жидкости с примесями. Предположим, что движение жидкости не зависит от присутствующих в ней примесей, т.к. дисперсный состав и концентрация примесей очень малы. Считаем частичку примеси шаром радиуса R, тогда, можно рассмотреть силы, действующие со стороны жидкости на эту частицу: FA - выталкивающая сила (Архимедова), P - сила давления, F - сила сопротивления, mg - сила тяжести, U - скорость движения частицы. Считаем, что жидкость линейно-вязкая, т.е. представляет собой сплошную среду, в которой сопротивление сдвигу пропорционально скорости движения.
где V - объем частицы примеси, Движение частицы будет выражаться следующим уравнением
где m - масса частицы (кг). Добавим начальные и граничные условия
Уравнения (2.20), (2.21) и (2.23) позволяют определить характер движения жидкости с примесями. Данные системы решаются с помощью метода конечных элементов. Пример распределения примесей в баке можно видеть на рис 2.7
Рис 2.6 Распределение примесей в баке (v0=2,358 м/с)
Идея метода конечных элементов состоит в том, что любую непрерывную величину, такую, как температура, давление и перемещение, можно аппроксимировать дискретной моделью, которая строится на множестве кусочно-непрерывных функций, определенных на конечном числе подобластей. Кусочно-непрерывные функции определяются с помощью значений непрерывной величины в конечном числе точек рассматриваемой области [9]. В общем случае непрерывная величина заранее неизвестна и нужно определить значения этой величины в некоторых внутренних точках области. Дискретную модель, однако, очень легко построить, если с начало предположить, что числовые значения этой величины в каждой внутренней точке области известны. После этого можно перейти к более общему случаю. Итак, при построении дискретной модели непрерывной величины поступают следующим образом [7]: ·в рассматриваемой области фиксируется конечное число точек. Эти точки называются узловыми точками или просто узлами. ·значения непрерывной величины в каждой узловой точке считается переменной, которая должна быть определена. ·область определения непрерывной величины разбивается на конечное число подобластей называемых элементами. Эти элементы имеют общие узловые точек в совокупности аппроксимируют области. ·непрерывная величина аппроксимируется на каждом элементе полиномом, который определяется с помощью узловых значений этой величины. Для каждого элемента определяется свой полином, но полиномы подбираются таким образом, чтобы сохранялась непрерывность величины вдоль границ элемента (этот полином называется функцией элемента). Рассмотрим трехмерный элемент, представляющий собой параллелограм. Он имеет четыре узла. Обозначим их индексами i, j, k, l, причем обход узлов в том порядке, как они написаны, осуществляется против часовой стрелки. Интерполяционный многочлен имеет вид:
или в матричном виде: Коэффициенты
где вектор Мы получили систему уравнений 4x4. Запишем эту системы в матричной форме:
где Матрица [C] строится для каждого элемента (то есть у каждого элемента своя и зависит от координат узлов конкретного элемента). Запишем матрицу [C]:
Для примера рассмотрим элемент с координатами узлов I=(1,2), J=(2,1), то первые две строки матрицы будут:
Столбец
но так как нам,
Но как нам известно
где
так как
то есть мы получили аппроксимацию Примеры разбиения бака и трубопровода на элементы приведены ниже
Рис 2.7 Пример разбиения трубопровода на элементы
Рис 2.8 Пример разбиения бака на элементы Движение жидкости в трубопроводе аналогично движению в баке-отстойнике, за исключением того, что трубопровод имеет иную геометрию, поэтому будут иные граничные условия. Начальное значение скорости будет вычисляться по формуле (2.22). Как ведут себя примеси в трубопроводе показано на рис 2.8
Рис 2.9 Распределение примесей в трубопроводе
Воспользуйтесь поиском по сайту: ![]() ©2015 - 2025 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...
|