УДК 532
DOI 10.18413/2075-4639-2018-50-3-308-316
УРАВНЕНИЯ ТЕЧЕНИЯ МНОГОКОМПОНЕНТНОЙ СМЕСИ ВЯЗКИХ НЕСЖИМАЕМЫХ ЖИДКОСТЕЙ
EQUATIONS OF FLOW MULTICOMPONENT MIXTURE OF VISCOUS
INCOMPRESSIBLE FLUID
А.С. Кравчук, А.И. Кравчук A. S. Kravchuk, A. I. Kravchuk
Белорусский государственный университет, 220030, пр. Независимости, 4, Минск, Республика Беларусь
Belarusian State University, 220030, 4 Nezavisimosti Ave., Minsk, Republic of Belarus
E-mail: as-kravchuk7@yandex.ru
Аннотация
В связи с тем, что при выводе уравнений математической физики используется стандартный прием - вывод уравнений для выделенного элементарного объема среды, то это позволяет, используя методику, разработанную авторами, обобщать известные уравнения на случай композиционных сред с использованием объемных долей компонент, где неявно используется гипотеза о том, что объемные доли компонент являются дискретной случайной величиной, описывающей вероятности присутствия той или иной компоненты неоднородной среды в конкретной точке с заданными координатами как элементарного объема, так и исследуемой геометрической области в целом. Построенная система уравнений позволяет решать задачи гидродинамики, например, для многокомпонентных смесей взаимно-нерастворимых жидкостей.
Abstract
Due to the fact that in the derivation of equations of mathematical physics using the standard method -derivation of the equations for the selected elementary volume of the medium, it allows using the method developed by the authors to generalize known equations to the case of composite media with volume fraction of the component. Implicitly used hypothesis that the volume fractions of the components are discrete random variable that describes the probability of the presence of one component of an inhomo-geneous medium at a particular point with given coordinates at an elementary volume, as well as at geometric area as a whole. The constructed system of equations allows us to solve hydrodynamics problems, for example, for multi-component mixtures of mutually insoluble liquids.
Ключевые слова: несжимаемая жидкость; вязкость, объемная доля компонент; эффективные свойства многокомпонентной смеси; системы уравнений Эйлера для усредненных давлений, векторов скорости и массовых сил; системы уравнений Навье-Стокса для усредненных давлений, вязкостей, векторов скорости и массовых сил.
Keywords: incompressible fluid; the viscosity, the volume fraction of the components; effective properties of a multicomponent mixture; the system of Euler equations for the average pressure, vectors of velocity and the mass forces; Navier-Stokes equations for the averaged pressure, viscosity, vectors of velocity and mass forces.
Введение
Механика жидкостей является одной из основных прикладных дисциплин для многих специальностей физико-математического и инженерного образования. В настоящее время она достаточно хорошо разработана [1, 2].
Однако, как и в других современных разделах математической физики, в этом разделе механики актуальна задача по обобщению известных систем уравнений для однородных жидкостей на многокомпонентные смеси, которые являются аналогами композиционных тел в механике твердого тела [3].
Данная статья демонстрирует возможность прямого переноса методики, разработанной авторами для твердых композиционных тел на многокомпонентные смеси жидкостей [4].
1. Методика исследования
Основным методом вывода систем уравнений динамики является получение уравнений движения выделенного элементарного объема dV = dx • dy • dz (где dx , dy , dz -
малые приращения вдоль осей декартовой системы координат с осями 0X, 0y , 0z ) рассматриваемого объема сплошной среды [1]. Соответственно V - объем, в котором рассматривается движение жидкости.
В настоящее время этот подход получил повсеместное признание и используется при выводе всех уравнений математической физики [5-8].
В связи с гипотезой несжимаемости предполагается, что для плотности жидкости
p(x, y, z) выполняется простейшее характеристическое уравнение р(х, У, z)= const [2].
Особенности вывода уравнений для композиционных сред с использованием
объемных долей компонент
Будем называть представительным объемом композиционной среды наименьший объем, в котором физические характеристики остаются в среднем равными физическим характеристикам всего объема среды.
Отметим, что рассмотрение представительного малого объема среды допускает обобщение вывода данных уравнений практически без изменений на случай соответствующих уравнений для композиционного материала с использованием объемных долей [6, 7].
В отличие от сплошной среды уравнения не будут верны непосредственно при предельном переходе dV ^ 0, т. к. в теории композиционных материалов представительный объем композиционного материала dV ф 0 . Но с другой стороны он бесконечно мал по
сравнению со всем объемом V. Точнее 3dx • dy • dz (характерный размер представительного объема) мал по сравнению с наименьшим характерным объемом композиционного тела. При этом и во всем объеме тела, и в элементарном объеме dV концентрации компонент (объемные доли) равны /к (где k (к = 1 , n) - номер компоненты, при этом выпол-
n
нено условие нормировки ^ / = 1).
к =1
Аналогичная гипотеза о малости представительного объема многокомпонентной смеси жидкостей по сравнению с размерами объема, в котором исследуется течение, будет использоваться в дальнейшем. Очевидно, что для каждой компоненты представительного
объема многокомпонентной смеси жидкостей должны быть известны ее плотность Рк и
вязкость TJk, кроме того действующие на нее вектор массовых сил f к .
При использовании данных гипотез план вывода уравнений математической физики полностью сохраняется с точностью до замены понятия элементарного объема на представительный объем [6, 7].
Отметим также, что в смысле формулировки объемных долей /к их можно рассматривать также как дискретную случайную величину присутствия компоненты с номером к в точке с координатами (X, у, 2) Е V.
Обозначим через (х, у, 2) и Рк (х, у, 2) дискретные значения вектора скорости
и давления на к -ой компоненте среды. Домножая на /к и суммируя по к получаем средние по реализации композиции вектор скорости и давление (матожидания дискретных случайных величин):
n n
v(x,^= • vk&y,z), {p(x,y,z)) = ^yk ■ pk(x,z). (1)
k=1 k=1
Таким образом, стационарность обеспечивается требованием воспроизведения концентраций на всех уровнях объема от представительного dV до общего V. Эргодичность обеспечивается равенством среднего по представительному объему среднему по реализации (1):
(v(x, y, z )) = --JJJ v(x, y, z )dxdydz, (2)
JJJdxdydz (dV)
(dV)
p(x, y, z --JJJ p(x, y, z )dxdydz (3)
JJJ dxdydz {dV)
(dV)
2. Результаты исследования и их обсуждение
Система уравнений Навье-Стокса для безвихревого потока вязкой несжимаемой жидкости имеет вид [1, 2, 9,10]:
— = - — ■ grad (p) + f + — -V2 v, div v = 0, (4)
dt p p
где f = (fx , fy, fz) - вектор массовой силы (размерность ускорения), — - вязкость жидкости.
Из системы (4) при — = 0 можно получить систему уравнений Эйлера движения несжимаемой невязкой жидкости для безвихревого потока [1,11,12].
Усреднение уравнений Эйлера
Будем рассматривать представительный объем dV потока многокомпонентной смеси жидкостей, для которого верна система уравнений (4) при — = 0. По аналогии с гипотезой Фойгта в механике композиционных сред [6, 7] будем предполагать, что распределение градиента давления p в элементарном (представительном) объеме не зависит от
наличия в той или иной точке включений других жидкостей и является непрерывно дифференцируемой функцией. В общем случае будем считать, что вектора массовых сил /к различны для каждой из компонент. Из (4) получаем закон определения вектора скорости жидкостей для к -ой компоненты системы в виде:
к =—— ■ grad (р) + /к , Жу Ук = 0. Ж Рк
(5)
Домножая (5) на У к и суммируя по к, получаем для средних по Фойгту вектора скорости (у), плотности (р)у и вектора массовых сил (/)т/ [6, 7]:
V
V
Р- ■ (р)+(/\г, ^ =(б)
п
V
V
= Х
У к ■ Ук
Р
V
к=1
п
\"1
утк_ V 1=1Рк )
п
V
= X Ук ■ .
(7)
к=1
Отметим в дополнение к (6), что по определению непрерывная функция р = Рк одновременно является и средним значением по Фойгту:
п
п
п
(РV = ХГк • Рк = ХГк • Р = Р ■ХУк = Р . к=1 к=1 к=1
(8)
По аналогии с гипотезой Рейсса в механике композиционных сред [6, 7] будем предполагать, что распределение проекций вектора скорости Ук жидкости элементарного объема не зависит от наличия в той или иной точке включений разных материалов и является непрерывно дифференцируемой по времени функцией. Тогда для к -ой компоненты жидкости, аналогично (5) из (4), получаем:
Рк ' ^ = ~ (Рк ) + Рк • /к ,
div у = 0.
(9)
После домножения на объемные доли У к и суммирования с учетом аналогичного (8)
равенства для вектора скорости, получаем для средних по Рейссу значений систему уравнений:
Ж V
Я
Жг
рг- grad ( р)я )+(7)К,
' Я
Шу У = 0.
(10)
п п
РЯ = ^Ук ■Рк , (Р)я = ^Ук ■ Рк , к=1 к=1
Я
Р
1 п — V"■ХУк ■Рк ■ .
Я к=1
(11)
Далее используется первая гипотеза о том, что эффективный вектор скорости [V элемента гомогенизированной жидкости вычисляется как линейная комбинация ^ = а ■(у + (1 — а) ■ ^ средних по Фойгту ^ и по Рейссу ^^ (где 0 <а< 1 - вещественный коэффициент) при этом предполагается, что для эффективного давления выполнено равенство (р) = (Ру = (Р п . Тогда из (6) и (10) получаем:
({У
(12)
где
Р1а =
а
(1 -а)
4-1
(Ру Р)
Я У
ъ. ам л+(1 а 4)
(13)
Используя обратный набор гипотез ( (р) — а^ Ру +(1 — а)^ Р
'Я
и
V) = ( V) = ^ Уу^ ), получаем еще один вариант оценки эффективных коэффициентов системы уравнений Эйлера течения гомогенизированной несжимаемой жидкости:
л V
Р}а^~ = — grad ((Р )Л^ ^ Р
а
(НУ V = 0,
(14)
где
Р)а = а ■РV + (1 — а>РЯ' </ ■Ра=а{/)у -Р) V + (1 — а){/)Я ■ Р)я
(15)
Интегирируя по а на интервале [0,1] коэффициенты в (13) и (15), получаем вилку Кравчука-Тарасюка оценки эффективной плотности (^РР
несжимаемой многокомпо-
нентной смеси жидкостей и вектора массовых сил ненты [6-8]:
, действующих на ее компо-
Р> =
Р)Я »V _]в(Р)Я
Р) Я —Р) 1Р>
^ IV у
Р) = 1 \Р V +Р Я )
(16)
гиа/ а—'* %' ' "7)
В качестве оценки эффективной плотности жидкости рр и вектора массовых сил
, как и ранее [6, 7], предлагается использовать средние значения вилки Кравчука-Та-расюка:
Р-+ (18)
Очевидно, что если не зависит от номера компоненты к (например, если присутствует только сила тяжести), то соответствующие коэффициенты сильно упрощаются.
Уравнения (13) и (14) определяют средние плотность и вектор массовых сил для системы уравнений Эйлера движения многокомпонентной смеси несжимаемых жидкостей, построенных в смысле средних (матожиданий) по представительному объему вектора скорости и давлений.
Оценка средних коэффициентов вязкости в представительном объеме
Как и ранее будем исходить из анализа элементарного (представительного) объема [1].
Будем использовать уравнения Ньютона для связи касательных напряжений и с компонентами градиента скорости. Опуская очевидные выкладки сразу можно записать, что
при постоянной вязкости ] сила трения Р /г , действующая на все грани выделенного элементарного трехмерного объема, имеет выражение [1]:
р/г — ] • (у2 у) • . (19)
Будем предполагать, что распределение трения на границах выделенного элементарного объема является непрерывным и не зависит от номера к компоненты многокомпонентной жидкости. Тогда из (19) следует, что для к -ой компоненты будет выполняться:
ра —]к ' (у2 Vк ) • dxdyd2. (20)
Разделив (20) на ], домножив на У к и просуммировав по к получаем:
Р/г) — (]р • У2(V) ) • dxdydz, (21)
где
] р -
Г N-1
' П ^
I
Ук
V к —1]к)
(22)
Предполагая непрерывность поля скорости и независимость его от номера компоненты к, не повторяя рассуждения предыдущего параграфа, запишем среднее по Рейссу:
=M r -Г(v
(у 2 ^v^ j • dxdydz,
где
n
M R Mk ■
к=1
(23)
(24)
Дословно используя план (12)-(18), можно записать связь средней силы трения г / (его матожидания в смысле дискретной случайной величины) с матожиданием век-
(25)
тора скорости (V) гомогенизированной жидкости в виде (Рисунок 1):
(Ffr) ~ M ' (У^ V) j • dxdydz :
где
M = \ (M + Mj'
M R• M V ы
M R
'R
M = M V +M) R j(26)
оценки разброса
Рис. 1. Вилка Фойгта-Рейсса , и Кравчука-Тарасюка ]], ]]
вязкости двухкомпонентной смеси несжимаемых жидкостей (] — 0.1 мПа • С, ]2 — 0.1 мПа • с) в зависимости от концентраций компонент У :
1 - (17); 2 - ])К (18); 3 - ] (20); 4 - ]] (20)
Fig. 1. Fig 1. Fork Voigt-Reuss [(7)^, (7)^ J and Kravchuk-Tarasyuk [7/, 7) estimating the viscosity spread of a two-component mixture of incompressible liquids (7 = 0.1 mPa ■ c,
щ = 0.1 mPa • с), depending on the concentrations of the components
Y : 1 - Щ)р (17); 2 - Щ)R (18); 3 - Щ (20); 4 - Щ) (20)
Следуя плану вывода системы уравнений движения представительного объема для однородной жидкости, изложенному в [1], получим эту систему с использованием эффективных коэффициентов (18) и (26):
Р)^ • dV = -grad(ip))• dV + Р) •f • dV + (Ffr). (27)
Подставляя (25) в (27) и сокращая на dV = dxdydz, получаем систему уравнений в терминах матожидания всех физических величин, вычисленных для элементарного объема:
4)
1
grad«p))+( Л + (28)
Кроме того, должно быть выполнено в среднем условие неразрывности ( V) = 0
для многокомпонентной смеси.
Отметим, что коэффициенты (14) и (20) могут использоваться исследователями для решения краевых задач течения несжимаемой жидкости с помощью любых программных средств с учетом необходимости выполнения условия малости минимального линейного размера представительного объема многокомпонентной смеси несжимаемых жидкостей по сравнению с минимальным размером области, в которой исследуется течение.
Выводы
Авторами была обобщена система уравнений Навье-Стокса течения однородной вязкой несжимаемой жидкости на случай использования многокомпонентной смеси. Предполагается, что известны объемные доли компонент.
Также получены эффективные коэффициенты для системы уравнений течения композиционной вязкой несжимаемой жидкости с учетом массовых сил.
Построенное авторами решение позволяет решать задачи гидродинамики, например, для многокомпонентных смесей взаимно-нерастворимых жидкостей.
Кроме того, можно рассмотреть также поведение тонкодисперсных, в том числе магнитно-реологических, суспензий, предполагая, например, для одного твердого вещества, что его доля в составе суспензии равняется yj ф 0 , а его вязкость равна нулю
= 0 ) в коэффициентах (28). Очевидно, этот подход элементарно распространяется на композитные тонкодисперсные суспензии с помощью того же правила: если твердые компоненты имеют номера с j = 1, т, то необходимо, чтобы в коэффициентах (28) У j ф 0,
щj = 0 для любого j (j = 1, m).
Список литературы References
1. Курбатов Ю.Л., Шелудченко В.И., Кравцов В.В. 2003. Механика жидкости и газа. Севастополь: Вебер, 226 с.
Kurbatov Yu.L., Sheludchenko V.I., Kravtsov V.V. 2003. Mechanics of fluid and gas. Seva-stop: Weber, 226 p.
2. Альтшуль А.Д., Киселев П.Г. 1975. Гидравлика и аэродинамика, М.: Стройиздат, 323 с.
Altshul AD, Kiselev P.G. 1975. Hydraulics and aerodynamics. Moscow: Stroiizdat, 323 p.
3. Тарасюк И.А., Кравчук А.С. 2014. Сужение «вилки» Фойгта-Рейсса в теории упругих структурно неоднородных в среднем изотропных композиционных тел без применения вариационных принципов. APRIORI. Серия: Естественные и технические науки, № 3. Режим доступа: http: //apriori -j ournal .ru/seria2/3 -2014/Tarasyuk-Kravchuk.pdf
Tarasyuk IA, Kravchuk A.S. 2014. The narrowing of the Voigt-Reuss "plug" in the theory of elastic structurally inhomogeneous in average isotropic composite bodies without the use of variational principles. APRIORI. Series: Natural and Technical Sciences, № 3. Access mode: http://apriori-journal.ru/seria2/3-2014/Tarasyuk-Kravchuk.pdf
4. Гершуни Г.З., Жуховицкий Е.М. 1972. Конвективная устойчивость несжимаемой жидкости. Москва: Наука, 392 с.
Gershuni GZ, Zhukhovitsky Ye.M. 1972. Convective stability of an incompressible fluid. Moscow: Nauka, 392 p.
5. Араманович И.Г., Левин В.И. 1969. Уравнения математической физики. М.: Наука, 288 с.
Aramanovich I.G., Levin V.I. 1969. Equations of mathematical physics. Moscow: Nauka, 288 p.
6. Кравчук А.С., Кравчук А.И., Тарасюк И.А. 2016. Методика вычисления эффективных коэффициентов в уравнении теплопроводности композиционного тела. Вестник СпбГУ. Серия 4, 2 (60), № 4: 335-341.
Kravchuk A.S., Kravchuk A.I., Tarasyuk I.A. 2016. Method for calculating effective coefficients in the heat equation of a composite body. Bulletin of St. Petersburg State University. Series 4, 2 (60), No. 4: 335-341.
7. Кравчук А.С., Кравчук А.И., Попова Т.С. 2016. Уравнение диффузии композиционной смеси в композиционную среду. Инженерно-физический журнал, 89, № 4: 1041-1046.
Kravchuk A.S., Kravchuk A.I., Popova T.S. 2016. The equation of diffusion of a composite mixture into a composite medium. Engineering and Physics Journal, 89, No. 4: 1041-1046.
8. Eringen A.C. 1983. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J Appl Phys, 54, no. 9: 4703-4710.
9. Kim J.S., Lowengrub J.S. 2006. Interfaces and Multicomponent Fluids. Encyclopedia of Mathematical Physics. New York: Elsevier.
10.Frolov R. 2017. An Efficient Algorithm for the Multicomponent Compressible Navier-Stokes Equations in Low and High Mach Number Regimes. Computational Physics. Available at: https://pdfs.s HYPERLINK "https://pdfs.semanticscholar.org/92a8/b9bb68e98dafa471a76ea4f660871ab7e888.pdf' emanticscholar.org/92a8/b9bb68e98dafa471a76ea4f660871ab7e888.pdf
11.Лойцянский Л.Г. 2003. Механика жидкости и газа. Москва: Дрофа, 840 с.
Loitsyansky L.G. 2003. Mechanics of fluid and gas. Moscow: Drofa, 840 p.
12. Смайлов С.А., Кувшинов К.А. 2012. Механика жидкости и газа. Томск: Издательство Томского политехнического университета, 121 с.
Smailov SA, Kuvshinov K.A. 2012. Mechanics of fluid and gas. Tomsk: Tomsk Polytechnic University Press, 121 p.
13. Ландау Л.Д., Лифшиц Е.М. 1986. Гидродинамика. Москва: Наука, 736 с.
Landau L.D., Lifshitz E.M. 1986. Hydrodynamics. Moscow: Nauka, 736 p.
14.Xu S., Xu Z., Alber M. 2018. Three-phase Model of Visco-elastic Incompressible Fluid Flow and its Computational Implementation. Communications in Computational Physics. Available at: https://www.researchgate.net/publication/322959186_Three-phase_Model_of_Visco-elastic_ Incompress-ible_Fluid_Flow_and_its_Computational_Implementation.
15. Lake L.W. 1989. Enhanced Oil Recovery. Upper Saddle River, New Jersey: Prentice Hall, 550 p.