Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ON THE EVOLUTION OF SOUND WAVES IN WET POROUS MEDIA

Gimaltdinov I.K. 1 Dmitriev V.L. 1 Sitdikova L.F. 1
1 Sterlitamak branch Bashkir state university
1008 KB
Theoretical and experimental study of acoustic wave propagation in porous media is relevant and important for the development of the processes that accompany the use of modern technologies for the use of porous media. In this paper, a theoretical study of wave propagation in wet gas saturated porous media in the interfacial interaction forces, the heat exchange between the skeleton of the porous medium, liquid and gas. In the viscoelastic characteristics recorded skeleton of the porous medium, and it is believed that the thin liquid layer uniformly on the inner surface of its pores. Recorded total system of equations and physical relationships describing the propagation of acoustic waves in a wet porous medium. A dispersion relation, based on which investigated the effect of volume fraction of liquid, gas phase on the propagation of acoustic waves in the medium.
acoustic wave
porous media
the «fast» wave
«slow» wave
1. Zaslavsky Y.M. On the efficiency of excitation of fast and slow waves in the water Bio-and gas-saturated media. Technical Acoustics, 2002, no 2, pp. 1–12.
2. Nigmatulin R.I. Dynamics of multiphase media. Moscow, Nauka, 1987. 464 p.
3. Nigmatulin R.I. Fundamentals of mechanics of heterogeneous media. Moscow, Nauka, 1978. 336 p.
4. Rahmatulin H.A. Fundamentals of gas dynamics of interpenetrating motions of compressible media. J. Appl., 1956, Vol. 20, no 2, pp. 184–195.
5. Shagapov V.S., Husainov I.G., Dmitriev V.L. Propagation of linear waves in gas saturated porous media with the interfacial heat transfer. J. Appl., 2004., Vol. 45, no 4, pp. 114–120.
6. Crandall I.B. Theory of Vibrating System and Sound. N.Y. D. Van Naostrand, 1927. 172 p.

Пористые материалы широко распространены в природе, в современной технике и промышленности. Значительный интерес представляют исследования волновой динамики и акустики дисперсных сред применительно к проблеме подавления звуковых возмущений пористыми средами.

Первая математическая модель пористой среды, состоящей из твердого тела (скелета), пронизанного длинными узкими каналами (порами), заполненными жидкостью, была изложена в [6]. Однако она была весьма идеализированной и неполной, так как позволяла описать только характер затухания акустической волны и не объясняла другие значимые механические эффекты. Позже были разработаны несколько теоретических и полуэмпирических моделей, которые также имели весьма существенные недостатки. С их помощью не удавалось связать частотные характеристики (инерционные и диссипативные), скорость звука и его затухание с измеряемыми структурными элементами сред и физическими свойствами материалов.

Статья Х.А. Рахматулина «Основы газодинамики взаимопроникающих движений сжимаемых сред» [4] положила начало новому этапу развития механики многофазных сред. В ней объединились достижения, идеи и методы газодинамики, механики взаимопроникающих континуумов и реальных процессов межфазного взаимодействия.

Уравнения многоскоростного движения и тепломассообмена в многофазных средах, а также выражения для внутренних и межфазных взаимодействий в плотноупакованных зернистых, порошкообразных и пористых средах были впоследствии получены Р.И. Нигматулиным на основе пространственного осреднения [2, 3].

Теоретический анализ характеристик акустических волн Био, распространяющихся в пористых средах, применительно к случаям насыщения порового пространства жидкостью и в случае газового заполнения пор, выполнен в работе [1]. Численное исследование дисперсии фазовой скорости и коэффициента поглощения продольных волн проведено на основе теории Био для упругих волновых процессов в двухкомпонентной среде.

В сухой насыщенной газом пористой среде влияние теплообменных процессов на распространение волн исследовано в работе [5]. Указаны области частот, когда затухание волн в насыщенной газом пористой среде определяется в основном теплообменными процессами.

pic_46.tif

Рис. 1. Ячейка пористой среды

Основные уравнения

При описании распространения одномерных волн во влажной пористой среде примем следующие допущения: значения длин рассматриваемых волн намного больше размеров пор; скорости жидкой пленки и скелета при прохождении волны равны (υl = υs). В качестве характерных размеров среды примем средний радиус пор a0, среднюю толщину водной пленки h0 и среднюю полутолщину стенок пор b0 (рис. 1).

Система уравнений, описывающих рассматриваемый процесс, имеет следующий вид:

Eqn125.wmf

Eqn126.wmf

Eqn127.wmf (1)

Eqn128.wmf Eqn129.wmf Eqn130.wmf

Здесь pg – давление в газовой фазе; αs, αl и αg – объемные содержания твердой, жидкой и газовой фаз соответственно; σs – напряжение в скелете; Fm – сила присоединенных масс, вызванная инерционным взаимодействием фаз; Fμ – аналог силы вязкого трения Стокса; FB – аналог силы Бассэ, проявляющейся при высоких частотах из-за нестационарности вязкого пограничного слоя около границы с твердой фазой; μg – динамическая вязкость газа; ρj, Eqn131.wmf, υ j, pj, αj – средняя по объему и средняя по фазе плотности, скорость, давление, объемные содержания.

Дополнительным нижним индексом (0) определены параметры, соответствующие невозмущенному состоянию, а параметры без индекса выражают малые возмущения параметров от равновесного значения; верхний индекс (0) соответствует истинному значению параметра.

Для описания неоднородностей температуры произведем схематизацию структуры среды, используя ячеистую схему. При этом пористую среду, насыщенную газом, примем как систему сферических газовых пузырьков, окруженных слоем жидкости и материала скелета. Таким образом, в каждой макроскопической точке, определяемой координатой x, вводим типичную ячейку, состоящую из газового пузырька со слоем жидкости и приходящегося на него скелета. Внутри ячейки имеется распределение микропараметров, а именно температуры Eqn132.wmf и плотности газа Eqn133.wmf где r – микрокоордината, отсчитываемая от центра ячейки. В дальнейшем штрихами наверху снабжены микропараметры.

Для истинных плотностей Eqn131.wmf и объемных содержаний αj справедливы следующие кинематические соотношения:

Eqn134.wmf

Eqn135.wmf

Распределения температур в ячейке пористой среды получим на основе системы уравнений теплопроводности:

Eqn136.wmf

Eqn137.wmf (2)

Eqn138.wmf

где λj и cj – соответственно коэффициенты теплопроводности и удельной теплоемкости при постоянном давлении (j = g, l, s).

Граничные условия для данной системы имеют вид:

Eqn139.wmf Eqn140.wmf

Eqn141.wmf Eqn142.wmf (3)

Eqn143.wmf

Решение системы уравнений будем искать в виде затухающих бегущих волн:

Eqn144.wmf (4)

Eqn145.wmf

где ω – круговая частота; K – комплексное волновое число; δ – коэффициент затухания.

Из условия существования решения вида (4) получим дисперсионное соотношение, на основе которого можно проанализировать характеристики распространения волн:

Eqn146.wmf

где Eqn147.wmf Eqn148.wmf Eqn149.wmf

Eqn150.wmf

Eqn151.wmf

Eqn152.wmf

Eqn153.wmf

Eqn154.wmf Eqn155.wmf Eqn156.wmf Eqn157.wmf Eqn158.wmf Eqn159.wmf

Eqn160.wmf Eqn161.wmf Eqn162.wmf Eqn163.wmf

Eqn164.wmf Eqn165.wmf Eqn167.wmf

Eqn168.wmf Eqn169.wmf

Eqn170.wmf

Результаты расчета

В расчетах параметры фаз взяты при температуре среды 300 К. Для воздуха:

γ = 1,4, Eqn171.wmf, cg = 1006 Дж/(кг∙К), λg = 0,027 Дж/(м∙с∙К), μg = 1,86∙10–5 Па∙с.

Для резины:

μs = 108 Па∙с, λs = 0,15 Дж/(м∙с∙К), cs = 1571 Дж/(кг∙К), Eqn172.wmf, Es = 108 Па.

Для воды:

Eqn173.wmf, cl = 4200 Дж/(кг∙К), λl = 0,6 Дж/(м∙с∙К).

На рис. 2 показаны зависимости коэффициента затухания δ и фазовой скорости Cp «быстрой» и «медленной» волн от частоты ω для системы «резина – вода – воздух»; характерные размеры среды a0 = 103 м. Линия 1 – жидкая фаза отсутствует, αs0 = 0,2, αg0 = 0,8; линия 2 – αl0 = 0,07, αs0 = 0,13, αg0 = 0,8.

pic_47.tif pic_48.tif

Рис. 2. Влияние водной пленки на коэффициент затухания и фазовую скорость «медленной» (сплошные линии) и «быстрой» (пунктирные линии) волн

Видно, что наличие жидкости приводит к тому, что для «медленной» волны в диапазоне частот ω ≤ 80 с–1, а для «быстрой» волны в диапазоне частот ω ≤ 104 с–1 коэффициент затухания больше, чем в случае «сухой» пористой среды. Это связано с тем, что водная пленка является дополнительным «балластом» для скелета и, таким образом, как бы увеличивает его массовую долю. Также из рис. 2 видно, что скорость «быстрой» волны из-за «утяжеления» скелета водой уменьшается на 68 м/с. Для «медленной» волны при частотах ω ≤ 10 с–1 наличие водной пленки приводит к уменьшению фазовой скорости, а для высоких частот к увеличению (около 47 м/с) скорости звука по сравнению со случаем, когда водной пленки нет.

На рис. 3 представлены зависимости коэффициента затухания и фазовой скорости обеих волн от частоты для разных объемных содержаний газовой фазы при постоянном объемном содержании скелета. Линии 1 – a0 = 10–3 м, b0 = 4,5∙10–5 м, h0 = 2,8∙10–5 м αg0 = 0,8, αs0 = 0,13, αl0 = 0,07, линии 2 – a0 = 10–3 м, b0 = 5,2∙10–5 м, h0 = 9,1∙10–5 м αg0 = 0,67 αs0 = 0,13, αl0 = 0,2.

Видно, что изменение объемного содержания жидкой фазы с 0,07 до 0,2 приводит к росту коэффициента затухания «быстрой» и «медленной» волн. При этом скорость «быстрой» волны уменьшается на 60 м/с, что связано с увеличением доли жидкой фазы в пористой среде. Скорость «медленной» волны при меньшем объемном содержании газовой фазы меньше в диапазоне частот ω ≤ 30 с–1, тогда как в области частот ω > 30 с–1 скорость «медленной» волны постепенно растет, достигая на высоких частотах различия до 33 м/с.

pic_49.tif pic_50.tif

Рис. 3. Влияние объемного содержания газовой фазы на коэффициент затухания и фазовую скорость и «медленной» (сплошные линии) и «быстрой» (пунктирные линии) волн

Стоит отметить, что увеличение объемного содержания скелета приводит к росту скорости «быстрой» волны, что связано с упругими свойствами скелета. Коэффициент затухания «медленной» волны при этом выше во всем диапазоне частот. «Быстрая» волна затухает меньше в случае большего объемного содержания скелета в диапазоне частот ω ≤ ω*, тогда как в области частот ω > ω* картина затухания меняется (здесь ω* – некоторая характерная частота, зависящая от объемного содержания скелета среды).

Выводы

В результате проведенных в работе вычислений установлено, что:

  • увеличение объемного содержания водной пленки при постоянном объемном содержании газовой фазы приводит к росту коэффициента затухания и уменьшению скорости «быстрой» волны, что связано с эффектом «утяжеления» скелета водой;
  • увеличение объемного содержания скелета при постоянном объемном содержании жидкости приводит к росту скорости «быстрой» волны, что связано с вязко-упругими свойствами скелета;
  • уменьшение массовой доли газа и увеличение массовой доли воды приводит к снижению фазовой скорости «быстрой» волны, что связано с ролью жидкости как «утяжелителя» скелета. При этом скорость «медленной» волны уменьшается на низких частотах и увеличивается на высоких частотах.

Рецензенты:

Михайлов П.Н., д.ф-м.н., профессор, заведующий кафедрой алгебры, геометрии и методики обучения математике, Стерлитамакский филиал ФГБОУ ВПО «Башкирский государственный университет», г. Стерлитамак;

Биккулова Н.Н., д.ф-м.н., профессор, заведующий кафедрой общей физики, Стерлитамакский филиал ФГБОУ ВПО «Башкирский государственный университет», г. Стерлитамак.

Работа поступила в редакцию 30.10.2013.