модель трансциллятора. аи. Филиппов, эв. Мухаметзянов, А. Ф. Садриев
 Скачать 1.69 Mb.
|
|
ISSN 1998-4812 Вестник Башкирского университета. 2013. Т. 18. №2 359 * автор, ответственный за переписку УДК 51.73 МАТЕМАТИЧЕСКАЯ МОДЕЛЬ КЛАССИЧЕСКОГО ЛИНЕЙНОГО ТРАНСЦИЛЛЯТОРА © АИ. Филиппов, ЭВ. Мухаметзянов, А. Ф. Садриев*, АИ. Леонтьев, Л. Ф. Садыкова Башкирский государственный университет, Стерлитамакский филиал Россия, Республика Башкортостан, 453103 г. Стерлитамак, пр. Ленина, а. Тел./факс: +7 (3473) 43 22 50. E-mail: sadriev_ainur@mail.ru Создана математическая модель, позволяющая объяснить процессы переноса при колебательных смещениях участков сплошной среды. Установлена зависимость коэффициента трансцилляторного переноса от частоты колебаний. Доказано, что величина коэффициента трансцилляторного переноса пропорциональна квадрату амплитуды колебаний и этот факт является прямым указанием на необратимость процесса, что соответствует природе теплопроводности. Ключевые слова: трансцилляторный перенос теплоты, коэффициент трансциллятор- ного переноса, акустическое поле Введение К настоящему времени накоплено много экспериментальных данных, свидетельствующих о том, что при колебательных смещениях участков сплошной среды возникает дополнительный перенос тепла или массы. Однако теоретически это явление недостаточно изучено. Нами создана математическая модель, позволяющая объяснить процессы переноса в таких системах. Уравнения, описывающие процессы, в общем случае представляют собой нелинейные дифференциальные уравнения с периодическими коэффициентами, при получении аналитических решений этих уравнений возникают большие трудности. Одним из важных достижений в этой области является создание теоретической модели элементарного переносчика – трансциллятора – и исследование переносных свойств. В простейшем случае трансциллятор (transcillator = transfer+oscillator) – пространственно протяженная физическая система, состоящая не менее чем из двух тонких пластин, хотя бы одна из пластин колеблется относительно других. Пластины взаимодействуют друг с другом, обмениваясь скалярным интегральным параметром Q (теплом, массой, зарядом и т.п.). Трансциллятор, пластины которого обмениваются между собой скалярным интегральным параметром по линейному закону Ньютона, назовем линейным. Трансцилляторный перенос, как отмечалось выше, возникает только при колебаниях компонент относительно друг друга. Если компоненты транс- циллятора совершают гармонические колебания, то такой трансциллятор назовем гармоническим. Если скорость колебаний компонент представляет собой сумму гармонических колебаний, то такой транс- циллятор назовем полигармоническим [1–4]. Описание трансцилляторного переноса Поясним существо трансциллятора на примере теплового поля в системе двух соприкасающихся и обменивающихся теплом пластин, одну из которых считаем для простоты неподвижной, а вторая – нижняя – может перемещаться вдоль оси х , расположенной на границе раздела, по колебательной зависимости (рис. 1). Между пластинами происходит теплообмен по закону Ньютона. Существо идеи трансцилляторного переноса заключается в том, что при наличии градиента температуры в системе колебания пластин вдоль оси приводят к дополнительной обычному молекулярному процессу теплопроводности. Для простоты ниже обычной теплопроводностью пренебрегается. Толщина пластин считается намного меньшей, чем амплитуда колебаний. Время релаксации температурного поля между пластинами сопоставимо с периодом колебаний. Рис. 1. Простейшая модель трансциллятора. Нетрудно видеть, что смещение нижней пластины соответствует просто колебательным смещениям участков среды. Поскольку такие смещения широко распространены в природе, то эта модель, называемая трансциллятором, описывает широкий класс явлений переноса тепла, массы, заряда, импульса и т.д. в реальных условиях. Она применима для объяснения влияния акустического поляна явления переноса [5], имеет важное значение для развития фундаментальных идей физики, например, теории турбулентности [6], теории хаоса [7] и т.д. Одним из практических достижений изложенной ниже теории явилось развитие расчетных методов для физических констант явлений переноса, инициированных всплыванием пузырьков в жидко 360 ФИЗИКА сти, особенно при наложении внешнего акустического поля. Всем известно, насколько важно использование технических устройств в технике и химической технологии. Механизм переноса при колебательном движении жидкости в пористой среде проще представить, если предположить, что изменение их взаимного расположения происходит скачкообразно. На рис. 1 изображены четыре состояния колеблющейся среды (нижняя пластина. До начала процесса нагревается часть неподвижной среды (на рис. 1 заштриховано. Из-за разности температур в состоянии происходит переток тепла из неподвижной среды в подвижную в направлении, показанном стрелочкой, соответствующая зона подвижной среды нагревается. В положении 2 после смещения колеблющейся среды наблюдается переток тепла из нагретой части подвижной пластины в более холодную часть неподвижной пластины. Фронт температурного возмущения (пунктирная линия) при этом сдвигается вправо. В положениях 3 и 4 происходит дальнейший сдвиг возмущения. На рис. 1 изображена начальная стадия процесса. Вдаль- нейшем теплообмен, обусловленный колебанием одной из пластин, стабилизируется. В приведенных рассуждениях пренебрегается молекулярной теплопроводностью вдоль пластин. Конвективный перенос тепла существенен только на расстояниях, сравнимых с амплитудой колебаний. Тем не менее, тепловые возмущения распространяются со временем на неограниченное расстояние. Интересно теоретически рассмотреть особенности этого процесса, носящего коллективный характер. При определенных условиях величина транс- цилляторного потока может на несколько порядков превышать молекулярный. Вычисление коэффициента трансцилляторно- го переноса производится следующим образом. Сначала определяется величина конвективного потока тепла с conv. r r ρ = (1) Для определения коэффициента трансцилля- торного переноса осуществляется осреднение выражения по периоду колебаний и по пространственной ячейке. Несмотря на то, что среднее значение скорости при колебательном движении равно нулю 0 v = r , среднее значение конвективного потока, вообще говоря, неравно нулю 0 v ≠ T r , поскольку температурное поле Т зависит от скорости v r , те. они коррелированны. Для определения коэффициента трансцилляторного переноса тр. λ необходимо величину усредненного конвективного потока тепла conv. j r привести к виду, аналогичному закону теплопроводности Фурье T j ∇ λ − = tr. conv. r (2) Одномерный гармонический трансциллятор Для неограниченной в направлении оси х системы уравнения, описывающие перенос локального параметра – температуры, преобразуются к виду ) ( ) 0 1 2 1 2 2 2 Т, (3) ( ) ( ) 0 1 1 2 1 Т) Для простоты сначала будем считать, что в начальный момент температуры пластин одинаковы, а антиградиент температуры в обеих средах совпадает и не зависит от времени и пространственной координаты 0 2 0 1 = = = t t T T , Г 1 , (5) индексы 1 и 2 относятся к верхней и нижней пластинам соответственно. Система уравнений (3) – (4) в этом случае допускает аналитическое решение. Поделив (3) и (4) на ж и ( ) c m − 1 соответственно и вычитая второе из первого, получим следующее уравнение для разности температур ( ) ( ) ( ) Γ = − + − t u T T dt T T d 1 2 1 2 α , (6) где ( ) [ ] ( ) { } 1 2 1 1 1 1 − − + − = α mc c m k . Решение этого уравнения с учетом условий (5) имеет следующий вид ( ) ( ) ∫ τ − τ γ τ Γ = − t d t k u T T 0 1 1 2 exp (7) Такое представление температуры в виде свертки свидетельствует о том, что описываемая система обладает свойствами памяти, следствием которой и является рассматриваемый процесс трансцилляторого переноса. Сложив (3) и (4), получим уравнение ( ) [ ] ( ) 0 1 2 1 1 2 2 = Γ + ∂ − + ∂ t u mc t Т c m Т mc , (8) которое имеет точное решение ( ) ( ) z c d u mc Т c m Т mc t Γ − ∫ τ τ Γ − = − + f 0 2 1 1 2 2 1 (9) Из (7), (9) определяются решения в следующей форме ( ) ( ) [ ] ( ) , exp 0 f 2 0 f 2 1 x d u c mc d t u c mc T t t Γ − Γ + + − Γ − = ∫ ∫ τ τ τ τ α τ (10) ( ) ( ) ( ) [ ] ( ) x d u c mc d t u c c m T t t Γ − Γ + + − Γ − = ∫ ∫ 0 f 2 0 f 1 2 exp 1 τ τ τ τ α τ (11) Здесь ( ) 2 1 f 1 mc c m c + − = – полная объемная теплоемкость трансциллятора. При гармонических колебаниях t A u ω ω = cos выражения для температуры имеют вид , sin f 2 f 2 1 x t A c mc I c mc T Γ − ω Γ + Γ − = (12) ISSN 1998-4812 Вестник Башкирского университета. 2013. Тж где ( ) ( ) [ ] ( ) − − + + = = − = ∫ t t t A d t u I t α ω α ω ω ω α ωα τ τ α τ exp sin cos exp 2 2 0 (13) Вдоль оси x , несмотря на отсутствие теплопроводности, происходит перенос тепла. Он осуществляется за счет конвективного процесса. Величина конвективного потока 2 2 muT c j = (14) в общем случае содержит колебательные компоненты. Несмотря на то, что при колебательном движении среднее значение скорости равно нулю, среднее значение конвективного потока за период колебаний отлично от нуля, поэтому в этом случае возникает дополнительный перенос тепла, называемый трансцилляторным переносом. Вычислим эффективный коэффициент трансцилляторного переноса. Для этого усредним конвективный поток по периоду колебаний в предположении больших времен 2 2 uT m c j = . В результате с учетом (11) имеем ( ) ( ) ( ) [ ] Γ ∫ τ − τ α τ − = t d t u u c cmc m j 0 f 2 exp 1 (15) Сопоставив полученное уравнение с законом Фурье, получим выражение для коэффициента трансцилляторного переноса ( ) ( ) ( ) ( ) [ ] ∫ τ − τ α τ − = λ t d t u t u c cmc m 0 f 2 tr exp 1 (16) Для гармонических колебаний ( ) 0 cos x t A t X + ω − = , ( ) t A t u ω ω = sin подстановка в (16) приводит к следующему результату ( ) ( ) = − − + × × + − = t t t t A c cmc m α ω α ω ω ω ω α α ω λ exp sin cos cos 1 2 2 2 2 f 2 tr ( ) 2 1 2 2 2 2 f 2 ω + α αω − = A c cmc m (17) Выразимтрансцилляторный коэффициент переноса, являющийся аналогом коэффициента теплопроводности, через Φ=α/ω = k1/(γω) ( ) 2 2 1 tr 1 2 Φ + = λ A k (18 ) Из выражения (18) следует, что при коэффициент трансцилляторного переноса достигает наибольшего значения 2 / 2 1 * m A k = λ . Существует характерная частота колебаний α = ω 0 , при которой коэффициент переноса составляет половину максимального. При малых частотах α << ω величина трансцилляторного коэффициента tr λ стремится к нулю. Значение граничной частоты увеличивается с ростом коэффициента ( ) [ ] ( ) { } 1 2 1 1 1 1 − − + − = α mc c m k , те. коэффициента теплообмена между пластинами 1 k и уменьшением эффективной трансцилляторной теплоемкости, определяемой согласно зависимости ( ) [ ] ( ) { } 1 1 2 1 1 1 − − − + − = γ mc c m Трансцилляторный перенос тепла играет определяющую роль в турбулентных потоках жидкости или газа, в пузырьковых жидкостях или в пористых средах, где колебания участков жидкости естественно обусловлены извилистостью пор или могут искусственно вызваны наложением акустического поля. Если для пористой среды – песчаника, насыщенного водой, положить 10 мВт, Дж K м 10 2 3 6 1 ⋅ ⋅ = γ − − , то значения граничной частоты составят кГц 2 / 2 1 1 Причину установленной зависимости транс- цилляторного коэффициента от частоты легко понять, представив разность температур (7) для гармонического колебания с частотой ω и амплитудой А в виде ( ) ( ) − Φ + Φ − − + Φ + Γ = − t t A T T α ϕ ω exp 1 sin 1 2 2 1 2 , (19) где ( ) , / 1 γω = Φ k сдвиг фаз ϕ находится из равенств 2 1 1 cos Φ + = ϕ , 2 Последнее слагаемое в (19) уменьшается со временем, что дает возможность при выполнении неравенств ( ) 1 − Φ ω >> t , 1 1 − γ >> k t перейти к выражению) Эта разность температур определяет тепловой поток между жидкостью и скелетом в установившемся тепловом режиме, когда процесс поглощения тепла средой завершен. Из полученного выражения следует, что между скоростью жидкости и тепловым потоком возникает разность фаз φ. Теплоперенос в пористой среде зависит от сдвига фаз между колебаниями жидкости и теплового потока между насыщающим флюидом и скелетом. Из следует, что на низких частотах указанный сдвиг фаз равен 0, теплоперенос отсутствует. С увеличением частоты тепловой поток вдоль оси х возрастает. Из (18) следует, что величина коэффициента трансцилляторного переноса пропорциональна квадрату амплитуды колебаний. Это означает, что величина tr λ инвариантна по отношению к знаку пространственного параметра – амплитуды колебаний А. Этот факт является прямым указанием на необратимость процесса, что, согласно существующим представлениям, соответствует природе теплопроводности. 362 ФИЗИКА Это обстоятельство оттеняет преимущества развитого нами подхода, поскольку все известные методы в физике приводят к неинвариантным соотношениям по отношению к знаку пространственного параметра. Например, коэффициент теплопроводности газа, согласно молекулярно-кинетической теории, пропорционален первой степени средней длины свободного пробега, те. не инвариантен по отношению к ее знаку. Это находится в серьезном противоречии с необратимостью процессов теплопроводности, согласно которой, в каком бы направлении молекула не двигалась, процесс переноса тепла все равно будет происходить от более нагретых частей к менее нагретым. Видимо, указанное обстоятельство делает необходимым уточнение классической теории переноса. Это легко осуществимо с помощью развитых в данной статье представлений. Действительно, движение молекул всегда можно представить в виде двух взаимопроникающих сред, движущихся навстречу, т.к. центр масс системы неподвижен. Поскольку при этом молекулы разных сред взаимодействуют, то такая система, конечно, может быть описана в рамках трансцилляторной модели. Проблема при этом сводится только к уточнению содержания констант. Это означает, что мы можем стать свидетелями еще одной увлекательной ревизии фундаментальных разделов физической кинетики. ЛИТЕРАТУРА Филиппов АИ, Минлибаев М. Р, Котельников В. А. Некоторые особенности явления вибропереноса тепла в пористых средах Теплофизика высоких температур. Т. 34. № 5. С. 719–723. Филиппов АИ, Котельников В. А, Минлибаев М. Р. Явление вибропереноса в двухкомпонентных осциллирующих взаимодействующих системах. // ИФЖ, 1997. Т. №3. С. 487–492. 3. ФилипповА. И. Особенности теплопереноса в пористой среде при возвратно-поступательном движении жидкости / Башк. гос. унт. Уфа, 1982. С. 10. 4. Буевич АС, Филиппов АИ. К явлениям переноса при колебаниях в двухкомпонентной среде // Инженерно- физический журнал. 1985. Т. XLVIII, № 2. С. 224–230. Филиппов АИ, Филиппов КА. О диффузии подвоз- действием звука // Акустический журнал. 1999. Т. 45. №3. С. 414–417. 6. Зельдович Я. Б.Точное решение задачи диффузии в периодическом поле скорости и турбулентная диффузия // ДАН СССР. 1982. Т. 266. № 4. С. 821–826. 7. Нигматулин Р. И, Филиппов АИ, Ахатов И. Ш, Ниязгу- лов С. А.Уравнения с периодическими коэффициентами и теория хаоса // Статика и динамика упорядоченных сред Межвузовск. научн. сб. Башк. унт. Уфа, 1994. С. 81–93. Поступила в редакцию 31.10.2012 г. |