К вопросу о моделировании термодинамического взаимодействия частиц, взвешенных в двумерной среде

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