УДК 532.5
ОПРЕДЕЛЕНИЕ МОДУЛЯ СДВИГА ЖИДКИХ СРЕД В ЭКСПЕРИМЕНТАХ С КРУТИЛЬНЫМ ВИСКОЗИМЕТРОМ
А.Е. Коренченко, В. П. Бескачко
Конечно-разностными методами получено решение задачи о движении крутильного вискозиметра, заполненного вязкоупругой жидкостью. Исследовано влияние упругих параметров жидкости на характеристики колебаний вискозиметра. Предложен метод идентификации вязкоупругих жидкостей, основанный на наблюдении резонансных эффектов.
Введение. Эксперименты [1] с диском, вращающимся по инерции на поверхности жидкости, указывают на существование упругих свойств у такой традиционно считающейся ньютоновской жидкости, как вода. Эти эксперименты подтверждают известный реологический тезис, что упругими свойствами в той или иной мере обладают все жидкости. Вопрос только в том, как их обнаружить, если они выражены так же слабо, как в воде. Для этого необходимы эксперименты, в которых реализуются весьма малые скорости деформации, например, такие, которые возникают в крутильном вискозиметре по мере затухания его колебаний. Изучение возможностей метода крутильных колебаний для определения упругих свойств «почти ньютоновских» жидкостей и является целью настоящей работы.
Математическая модель. Для описания вязкоупругих эффектов в работе используется представление о жидкости с конвективной упругостью [2]. Такая жидкость определяется как материал, для которого напряжение зависит как от деформации, так и от скорости деформации, причем деформация определяется через различия в конфигурации материала в последовательные моменты времени, а не в сравнении с некоторой предпочтительной формой. Тензор напряжений в такой жидкости может быть представлен в следующем виде: 0=Овяз + 0упр . Упругая часть тензора напряжений изотропной жидкости с конвективной упругостью выражается формулой
^ т ■ О ■ * (г').
В работе это соотношение использовалось в дискретной форме:
•"+1 = о ■ в+я" ■ е /т
упр О • т^упр е ,
где Я"пр - тензор упругих напряжений в момент времени t, О - модуль сдвига, В - тензор малых
деформаций, появившихся за время Дt = - ^ и учтено, что упругие напряжения за каждый интервал Дt не аддитивны, а имеет место максвелловская релаксация [3]. Вязкая часть тензора напряжений удовлетворяет уравнению Ньютона Явяз =-Ц ■ В, В - тензор скоростей деформации. Описанная модель является линейной моделью вязкоупругости и характеризуется помимо вязкости двумя параметрами: модулем сдвига О и временем релаксации упругих напряжений т.
Введем безразмерные переменные, когда все расстояния отнесены к радиусу цилиндра,
скорость - к V /(V - кинематическая вязкость жидкости), давление - к ру2/ 2 ( р - плотность жидкости), время - к ^ц2/у . Тогда в безразмерных координатах ( г, 0, г) движение цилиндра и заполняющей его жидкости, описывается уравнениями:
ди = - и ■ Уи-Ур - Гупр, (1)
V- и = 0, (2)
d соц
=-к-у + Мвяз + Мупр. (3)
-да
Здесь и _ (мг, мо, п2) и р - поля скорости и избыточного по сравнению с pgz давления в жидкости; юб - угловая скорость цилиндра; р - угол поворота цилиндра относительно равновесного К4к
положения; К _ —2— ; к - коэффициент крутильной жесткости нити; 1ц - момент инерции ци-v 1ц
линдра; Мвяз и Мупр - безразмерные моменты вязкой силы трения и упругой силы, приложенные к цилиндру со стороны жидкости. Безразмерная объемная плотность упругой силы Г на каждом временном шаге определяется явным образом по формуле
г г+Аг _ г г 1 упр _ гупр
■ е /т +
Е■Аи _(и-V)Г
г
упр
_Аг/
■ е /т ■ Аг,
здесь Е _
О^Аг
.2
, А - безразмерный интервал, соответствующий временному шагу разностной
схемы
Выражение для моментов М]
5 1
о
вяз и Мупр можно записать в виде
Мвяз _
2пр Д
I
ц
2 дм0_ дz
& + \ г
z _0 0
дмо
дz
н / д,,
& +
z _ н / Д
ц
дыв дг
щ г
dz),
г _1
(4)
2п Д^С -А/
М
t+At
упр
1цV 2
(4
_J г 2
0
дм,
дz
dг + Гг 2 ^
z _0
дz
н / Д,,
dг +
дмв дг
О
dz) + М^ р ■ е
_А//
г_1
Здесь Н - высота цилиндра и предполагается, что сверху жидкость соприкасается с твердой границей - крышкой.
В качестве граничных условий приняты условия прилипания жидкости ко всем твердым границам:
мв (1, z, г) _ Юц; мв (г, 0, г) _ мв (г, Н / Д, г) _ Юц ■ г; мг (г, 0, г) _ мг (г, 0, г) _ мг (г, н / Дц, г) _ м2 (г, н / Дц, г) _ 0, (5)
мг (1, z, г) _ м2 (1, z, г) _ 0.
Начальные условия записываются в виде
р(0) _ 0; Юц (0)_ Ю0, (6)
т.е. при г < 0 вискозиметр покоится вместе с жидкостью в положении равновесия и в момент г _ 0 толчком приобретает некоторую угловую скорость.
Численное решение находилось методом конечных разностей. Использовались равномерные пространственные сетки с максимальной размерностью 50 х 50 в осевом и радиальном направлениях. Дискретизация гидродинамических уравнений производилась по схеме центральных разностей по пространственным переменным и односторонних по времени. Линеаризация полученных в результате нелинейных алгебраических уравнений проводилась по схеме Ньютона, а решение линеаризованных систем уравнений на каждом временном этапе - методом исключения Гаусса [4]. При реализации схемы были приняты меры по предотвращению появления схемной вязкости. Использование для решения системы уравнений метода Гаусса приводит к существенным ограничениям на размеры сеточного разбиения, однако этот метод обладает хорошими консервативными качествами и гарантированной сходимостью. Консервативные свойства численного решения, а именно, равенство нулю интегрального массового потока, проверялись на каждом временном шаге.
Численные расчеты были проведены для цилиндров с различными радиусами, подвешенных на нитях с различной крутильной жесткостью и для модельных жидкостей с различными значениями упругих параметров. Момент инерции цилиндра выбирался так, чтобы в _ I / 1ж _ 10,
здесь 1ж - момент инерции «замороженной» жидкости. Внутренний радиус вискозиметра выбирался в пределах 0,008 ...0,021 м, а отношение высоты цилиндра к внутреннему радиусу у _ н / Дц было равно 1, так как в [5] было показано, что упругие свойства проявляются при про-
0
м
о
г
z _ н / Дц 0
0
Голубев Е.В., Гуревич С.Ю.
Лазерное возбуждение нормальных волн в ферромагнитной пластине
чих равных условиях тем сильнее, чем меньше у. Плотность и вязкость жидкости отвечали воде. Коэффициент крутильной жесткости нити выбирался так, чтобы период колебаний варьировался в пределах 7,5...18 с. Расчеты были проведены для модельных жидкостей, в которых величина
модуля сдвига выбиралась порядка 10-3 Па (слабые вязкоупругие свойства).
Результаты и обсуждение. Результатом расчета для каждого набора параметров был закон движения вискозиметра cp(t). Затем этот закон аппроксимировался функцией вида
<p(t) = A • e~p't sin (f • t + y) (7)
с целью определения параметров колебаний: частоты f и коэффициента затухания р. Подгонка осуществлялась методом наименьших квадратов с минимизацией методом Розенброка [6]. На рис. 1 изображены графики зависимостей коэффициента затухания от периода колебаний для ньютоновской жидкости и для жидкостей с тем же коэффициентом вязкости, но различными значениями упругих параметров. Для получения каждой кривой на рис. 1 варьировалась жесткость нити подвеса, влияющая на частоту колебаний вискозиметра и, таким образом, на период T = 2п / f. Прочие условия численного расчета (масса, внутренний радиус и момент инерции цилиндра, модуль сдвига, вязкость и плотность жидкости) были фиксированы.
Кривые на рис. 1 отличаются тем, что получены для жидких сред с различными временами релаксации. Как видно из рисунка, зависимости p (T)
для жидкостей с временами релаксации т = 30 мин и т = 3 мин различаются мало. Это можно объяснить тем, что характерное время колебаний - их период T не превышает 18 с и поэтому для рассмотренных выше времен релаксации T << т, так что оба эти времени могут считаться «асимптотически большими». Кривые, отвечающие этим временам релаксации, имеют выраженный максимум, располагающийся при T «12 с. Легко показать, что максимум коэффициента затухания наблюдается при совпадении радиуса цилиндра и расстояния, пробегаемого упругой волной за половину периода колебаний, т. е.
Рис. 1. Зависимость коэффициента затухания крутильных колебаний от их периода для модельных вязкоупругих сред с различными значениями времени релаксации; Яц = 0,0103 м,6 = 0,001 Па
T 2
Гупр • - = 1,038• 10 м*
Я,
(8)
где Купр =л]30/р - скорость распространения упругих сдвиговых возмущений в среде.
«Резонансный» характер зависимостей р (Т) является результатом наложения двух процессов : движения среды, обусловленного крутильными колебаниями цилиндра, и волнового движения, связанного с распространением в ней упругой волны. Анализ временной эволюции распределения скоростей в объеме жидкости показывает, что величина азимутальной скорости во внутренних областях вязкоупругой среды превышает соответствующие значения в ньютоновской жидкости и достигает наибольших значений при выполнении условия (8). Интенсификация относительного движения внутренних слоев жидкости увеличивает диссипацию энергии за счет внутреннего трения и приводит к увеличению коэффициента затухания. Вид зависимостей р (Т), полученных при ^ = 17 с и ^ = 3,5 с, т.е. для значений времен релаксации, сравнимых с периодом колебаний и меньших его, существенно изменяется. Во-первых, как видно из рисунка, максимум на кривых смещается в область меньших Т, и может совсем пропасть (как для т = 3,5 с), во-вторых,
коэффициент затухания колебаний изменяется при резонансе на тем меньшую величину, чем меньше время релаксации.
Для идентификации среды как вязкоупругой и определения в рамках данной модели модуля сдвига жидкости следует получить как можно ярче выраженный максимум на кривой зависимости коэффициента затухания от периода колебаний или от другого аргумента, зависящего от периода. Поэтому важно выяснить, как влияют параметры вискозиметра, например, радиус , на характер указанной кривой. На рис. 2 изображены зависимости коэффициента затухания колебаний от безразмерной комбинации £ = Уупр Т / (2) для различных . В соответствии с формулой (8), максимума затухания следует ожидать при значениях £ , близких к 1. Как видно из рисунка, функции р (£), полученные для радиусов = 0,00866 м и
0,01225 м, монотонно убывают и, таким образом, не несут информации о модуле сдвига. Максимум коэффициента затухания выражен тем лучше, чем больше радиус цилиндра вискозиметра. Заметим, что величина, на которую возрастает коэффициент затухания при «резонансе» для = 0,02021 м , достаточна
для уверенной регистрации в экспериментах с крутильным вискозиметром, где относительная погрешность измерения коэффициента затухания может быть сделана поряд-
Рис. 2. Зависимость коэффициента затухания от безразмерного параметра £ = Уупр-Т/2яц ; 0=0,0005 Па, т = 17 с
ка -10"
Таким образом, обнаруженные в рамках рассмотренной реологической модели «резонансные» эффекты позволяют идентифицировать жидкость с конвективной вязкоупругостью и дают возможность определить модуль сдвига по положению максимума коэффициента затухания на его зависимости от периода колебаний вискозиметра. При практической реализации метода период колебаний можно варьировать, изменяя крутильную жесткость нити подвеса, что легко достигается изменением ее длины.
Литература
1. Апакашев Р.А., Павлов В.В. Определение предела прочности и модуля сдвига воды при малых скоростях течения// Механика жидкости и газа. - 1997. - № 1. - С. 3-7.
2. Астарита Дж., Маруччи Дж. Основы гидромеханики неньютоновских жидкостей. - М.: Мир, 1978.-302 с.
3. Бартенев Г.М., Сандитов Д.С. Релаксационные процессы в стеклообразных системах. -Новосибирск: Наука, 1986. - 236 с.
4. Ши Д. Численные методы в задачах теплообмена. - М.: Мир, 1988. - 537 с.
5. Коренченко А.Е., Бескачко В.П. Об определении вязкоупругих свойств жидкости из анализа колебаний крутильного вискозиметра// Труды XI Российской конференции «Строение и свойства металлических и шлаковых расплавов». - Т. 2. Строение и свойства металлических расплавов. - Челябинск: Изд-во ЮУрГУ, 2004. - С. 126-130.
6. Химмельблау Д. Прикладное нелинейное программирование. - М.: Мир, 1975. - 534 с.
Поступила в редакцию 12 августа 2005 г.