Расчет полей скорости в соэкструзионных течениях методом конечных элементов
Метод конечных элементов позволяет осуществлять численное моделирование двух - и трехмерных многослойных течений (см. раздел 4.4.5). Отличие от моделирования однослойных течений состоит в том, что в данном случае требуется реализовать еще один алгоритм. Он заключается в том, что сетка конечных элементов сдвигается таким образом, чтобы граница между двумя слоями различных расплавов всегда находилась на границе элемента. Расчеты МКЭ требуют больших объемов дискового пространства, времени для расчетов, поэтому для моделирования необходим, как минимум, микрокомпьютер[28].
Специфические проблемы, связанные с моделированием многослойных течений МКЭ, здесь обсуждаться не будут. Заинтересованные читатели могут найти более подробную информацию по данному вопросу в работе [37]. Однако чтобы продемонстрировать возможности МКЭ, будет произведено сравнение результатов, полученных с помощью численного моделирования и на основе эксперимента.
Изучение области, в которой происходит слияние слоев в многослойном течении, производилось на экспериментальном испытательном стенде, обеспечивающем возможности визуального наблюдения за экспериментом [38]. В рамках проводимых экспериментов изучали формирование границы раздела между слоями и ее положение, а также проводили оценку других явлений, возникающих при совместном течении. В качестве следующего шага была сделана попытка численного моделирования двухслойного плоского течения МКЭ.
На рис. 6.18 представлены фотографии реальных течений и линии тока, полученные расчетным путем для различных соотношений объемных расходов расплавов. Следует отметить, что численное моделирование МКЭ позволяет не только хорошо описать форму и положение границы между слоями, но и дает возможность оценить размеры и положение застойных зон, возникающих в районе слияния слоев в зависимости от соотношения между объемными расходами расплавов. Этот расчет основывается на чисто вязких свойствах материала, описывающихся моделью Карро.
Появление застойной зоны объясняется значительным повышением скорости течения нижнего потока. Более быстрый верхний слой увлекает частицы нижнего слоя в направлении течения, вызывая соответствующее повышение давления в нижней части потока. В результате этого на небольшом расстоянии от места слияния слоев
VJ V2 = 2,5/1 |
Рис. 6.18. Сравнение экспериментальных результатов изучения плоского двухслойного течения с результатами численного моделирования с помощью МКЭ
в нижней части потока формируется локальный максимум давления. Это давление вынуждает расплав течь в обратном направлении, вследствие чего образуется застойная зона. Подобное явление наблюдается при нанесении изоляции на проводники (см. раздел 5.3.3.4).
Таким образом, моделирование с помощью МКЭ позволяет без проблем предсказать появление застойных зон в области слияния потоков различных расплавов в канале экструзионной головки. Это позволяет избежать увеличения времени пребывания расплавов в канале и соответствующие неблагоприятные воздействия на расплавы в этой зоне. Этот пример убедительно показывает, что МКЭ позволяет получить достаточно точное описание даже для весьма сложных процессов течения.