Numerical simulation of a non-isothermal power-law fluid flow in a channel with abrupt contraction

Проводится исследование ламинарного стационарного неизотермического течения степенной жидкости в цилиндрическом канале с внезапным сужением. Формулируется математическая модель течения, которая включает уравнения гидродинамики, записанные в переменных функция токавихрь, и уравнение энергии. Реологические свойства жидкости описываются степенным законом Оствальда де Ваале, в модифицированной форме которого учитывается зависимость эффективной вязкости от температуры. Для решения задачи используется метод установления с последующей реализацией конечноразностного метода на основе схемы переменных направлений. Выполняется оценка влияния вязкой диссипации на структуру потока псевдопластичной, ньютоновской и дилатантной жидкостей. Демонстрируются поля температуры и эффективной вязкости. Приводятся результаты параметрического исследования коэффициента местного гидравлического сопротивления. Viscous fluid flow through a sudden contraction is frequently encountered in a number of industrial equipment dealing with processing and transporting of liquid materials. Generally, the fluid exhibits non-Newtonian behavior and flows under non-isothermal conditions. Such flow is characterized by specific structure, viscous dissipation, and local pressure losses. In this work, a numerical solution to the problem of a non-isothermal power-law fluid flow through a two-to-one axisymmetric abrupt contraction is presented. Mathematical model includes the momentum, continuity and energy equations written in terms of stream function, vorticity and temperature variables. The rheological behavior of the fluid is specified by the Ostwald-de Waele power law. The proposed flow model accounts for viscous dissipation and temperature-dependent rheological properties. To solve the problem, the relaxation method is used, followed by the implementation of the finitedifference method based on the scheme of alternative directions. The equations in a discrete form are solved using the tridiagonal matrix algorithm. It is found that the flow includes both one-dimensional and two-dimensional flow regions. The lengths of these regions are studied versus the Reynolds number, the Peclet numder, and power-law index. Comparing isothermal and non-isothermal cases, it is revealed that an increase in the power-law index leads to a decrease in the downstream two-dimensional zone in the first case, and, in contrast, it provides a significant increase in the second case. The viscosity and temperature distributions are presented to show the effect of the Peclet number for pseudoplastic, Newtonian, and dilatant fluids. The parametric investigation of the local pressure losses is implemented in a wide range of the main parameters.