Распространение коротких акустических импульсов в средах с экспоненциальной и резонансной релаксацией Ларичев Владимир Андреевич

Содержание к диссертации

Введение

Глава 1 Аналитическое описание динамики короткого импульса распространяющегося в среде с произвольным СВР

1.1 Термодинамический подход Мандельштама и Леонтовича к описанию дисперсионно-диссипативных свойств сред в близи состояния термодинамического равновесия.

1.2 Связь динамики короткого импульса с параметрами релаксационной среды 77

1.3 Аппроксимирующее аналитическое выражение, описывающее динамику короткого импульса, распространяющегося в релаксационной среде с произвольным СВР 34

1.4 Примеры численных расчетов динамики импульса 41

Заключение к главе 1 46

Глава 2 Распространение акустического импульса в среде с двумя релаксационными процессами

2.1 Функция Грина задачи о распространении импульса в среде с двумя релаксационными процессами

2.2 Аналитическое выражение для функции Грина среды с двумя релаксационными процессами

2.3 Альтернативное представление функции Грина среды с двумя релаксационными процессами

Заключение к главе 2 66

Глава 3 Распространение импульсов в неоднородных релаксационных средах 68

3.1 Функции Грина линейного и точечного источников в экспоненциально неоднородной релаксационной среде

3.2 Распространение плоского импульса в неоднородной вдоль направления распространения релаксационной среде

Заключение к главе 3 86

Глава 4 Обобщенное локальное уравнение состояния среды с резонансной релаксацией 88

4.1 Обобщение термодинамического подхода на случай среды с резонансной релаксацией

4.2 Дисперсия фазовой скорости и частотная зависимость коэффициента поглощения в среде с резонансной релаксацией

4.3 Механическая интерпретация обобщенной функции отклика сред с резонансной релаксацией

Заключение к главе 4

Глава 5 Распространение коротких импульсов в среде с единственным процессом резонансной релаксации

5.1 Аналитическое пространственно- временное представление функции Грина задачи о распространении плоского импульса в среде с одним обобщенным процессом резонансной релаксации 110

5.2 Асимптотический анализ интегралов 116

5.2.1 Прифронтовые асимптотики 117

5.2.2 Асимптотики больших расстояний. Случай А > 0 119

5.2.3 Асимптотики больших расстояний. Случай А < 0 129

5.3 Основные типы распространения импульса в среде с одним обобщенным процессом резонансной релаксации

Заключение к главе 5 156

Заключение 157

Литература 158

• Аппроксимирующее аналитическое выражение, описывающее динамику короткого импульса, распространяющегося в релаксационной среде с произвольным СВР
• Аналитическое выражение для функции Грина среды с двумя релаксационными процессами
• Распространение плоского импульса в неоднородной вдоль направления распространения релаксационной среде
• Дисперсия фазовой скорости и частотная зависимость коэффициента поглощения в среде с резонансной релаксацией

Введение к работе

Диссертация посвящена проблеме распространения акустических импульсов малой амплитуды в средах, обладающих частотно-зависимым поглощением и дисперсией фазовой скорости. Поскольку поглощение и связанная с ним частотная дисперсия проявляется в той или иной мере практически во всех реальных средах, то это явление играет значительную роль для различных научных и технических приложений, так или иначе связанных с использованием как акустических, так и электромагнитных импульсов. В частности, проявлением дисперсионно-диссипативных свойств является то, что форма и амплитуда импульсов изменяется по мере их распространения, что позволяет использовать это обстоятельство как для диагностики свойств среды, так и для правильной интерпретации передаваемого сигнала.

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

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

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

Другая группа приложений связана с обратной задачей. Поскольку частотная дисперсия в среде определяется релаксационными процессами, происходящими на молекулярном или микроструктурном уровне, решение обратной задачи может быть использовано для определения параметров, характеризующих молекулярные или микроструктурные свойства сред, а также кинетику соответствующих релаксационных процессов. Такими величинами могут быть, например, параметры спектра времен релаксации ([Михайлов, Соловьев, Сырников 1964;Новик, Берри 1975;Кельберт, Чабан 1986;Кельберт, Сазонов 1987;Максимов 1996]). Таким образом, импульсная диагностика ([Кельберт, Сазонов 1987;Нигул 1981;Нигул 1984;Максимов 1996]), основанная на решении обратной задачи может служить инструментом молекулярной акустики и дополнением к традиционной акустической спектроскопии.

Заметим, что в англоязычной литературе использование закономерностей динамики формы импульса для определения каких-либо дисперсионных характеристик среды часто называют спектроскопией во временной области (time domain spectroscopy) (см. например [Roberts, Petropoulos 1996]).

Хотя в обычных условиях импульсная диагностика по-видимому может служить лишь дополнением к традиционным спектроскопическим методам, существуют задачи, в которых у импульсных методов есть значительные преимущества ([Кельберт, Сазонов 1987;Максимов 1996;Roberts, Petropoulos 1996]). Для традиционной спектроскопии обычно требуется сложный комплекс измерений в широком частотном диапазоне. Причем в процессе измерений необходимо поддерживать стационарные условия в течение относительно длительного времени. Это обстоятельство осложняет применение спектроскопии для диагностики быстропротекающих процессов и неустойчивых сред или сред, находящихся под воздействием быстроменяющихся внешних условий. В таких случаях импульсные методы могут иметь существенные преимущества.

Примерами задач, в которых целесообразно использовать импульсные методы диагностики могут служить задачи мезомасштабной диагностики атмосферы ([Иванченко, Николаев 2000]) и океана ([Кельберт, Сазонов 1987]). В этих задачах при проведении спектроскопических измерений возникает весь отмеченный выше комплекс сложностей, в тоже время импульсные измерения могут быть проведены сравнительно просто.

Следует отметить особую роль коротких импульсов в задачах импульсной диагностики. В таких задачах, так же как и в задачах локационного зондирования, важно отделить особенности динамики импульса, определяемые средой, от особенностей определяемых его начальным профилем. В работе [Максимов 1996] показано, что в случае когда длительность импульса много меньше чем характерные времена происходящих в среде релаксационных процессов, его динамика определяется главным образом дисперсионными свойствами среды и практически не зависит от его первоначальной формы. Это обстоятельство делает короткие импульсы удобным диагностическим инструментом.

Заметим также, что корректная постановка и решение обратной задачи - задачи импульсной диагностики - в значительной мере определяется возможностью нахождения решения прямой задачи, явно зависящего от параметров, характеризующих дисперсионно-диссипативные свойства среды. В этом случае обратная задача может быть поставлена как задача отыскания таких значений параметров, при которых решение прямой задачи было наиболее близко к профилю импульса, полученному экспериментально. Общность такой постановки определяется возможностью достаточно общего описания дисперсионно-диссипативных свойств реальных сред.

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

Способы описания дисперсионных свойств сред

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

Остановимся на основных подходах, используемых для описания дисперсионно-диссипативных свойств сред ([Коган 1963]).

Эмпирическое описание

Наиболее простым способом такого описания является прямое использование экспериментальных кривых частотной зависимости коэффициента поглощения и фазовой скорости, полученных в результате прямых спектроскопических измерений в определенной ограниченной области частот. Эти данные в дальнейшем могут быть экстраполированы на всю оставшуюся область частот в рамках той или иной модели, удовлетворяющей принципу причинности, что обычно достигается учетом дисперсионных соотношений Крамерса — Кронинга (см. например [Ландау, Лифшиц Т.VIII1992]). Такой способ не требует информации о физической природе механизмов, ответственных за дисперсионно-диссипативные свойства. Однако свойства среды, выявленные таким способом, являются: чисто эмпирическими, и, не могут быть использованы для другой среды, или той же среды, но находящейся в других условиях.

Примеры теоретического обоснования и использования этого подхода можно найти в работах [Дерягин 1931; Гинзбург 1955; Горшков 1957; Азими и др. 1968; Губкин 1984; Lamb 1962].

Механические модели

Другой подход к описанию дисперсионно-диссипативных свойств основан на моделировании реологии среды, т.е. ее механического поведения под влиянием внешнего воздействия. Несмотря на то что, ни закон Гука, ни закон вязкого течения Ньютона в отдельности в точности не описывают механическое поведение значительной части реальных сред даже в линейном приближении, качественное, а иногда и количественное согласие с экспериментами часто удается получить в предположении, что процессы упругого и вязкого деформирования протекают в среде одновременно. Механическое поведение материала при таком подходе может быть смоделировано совокупностью последовательно и параллельно соединенных упругих элементов (пружинки), подчиняющихся закону Гука и вязких элементов (демпферы, поршни), подчиняющихся закону вязкого течения Ньютона ([Работнов 1977; Новик, Берри 1975; Аскадский 1973]), а также сосредоточенных масс в качестве инерционных элементов ([Achenbach, Chao 1962; Новик, Берри 1975; Николаевский 1985]), если реакция материала имеет также и инерционные свойства.

Наиболее простые модели такого типа это совокупность упругого и вязкого элементов, соединенных последовательно (модель Максвелла) или параллельно (модель Кельвина-Фойгта), а также их обобщение — трехэлементная модель стандартного неупругого тела ([Новик, Берри 1975; Аскадский 1973]), представляющая собой модель Кельвина-Фойгта последовательно соединенную с пружиной. Эти модели как правило лишь качественно описывают поведение неупругих сред. Однако более сложными моделями такого типа, состоящими из большего числа элементов можно эффективно описать поведение целого ряда материалов ([Николаевский и др. 1970; Аскадский 1973; Новик, Берри 1975; Работнов 1977]).

Все возможные механические модели можно рассматривать как частные случаи общего подхода, положенного в основу наследственной теории упругости. Этот подход был предложен Больцманом ([Аскадский 2001]), который исходил из того, что если тело ранее испытывало деформацию, то повторная деформация до той же величины потребует меньшего напряжения, причем это уменьшение напряжения тем больше, чем больше длилась первичная деформация. Далее предполагается, что описанный выше эффект суммируется при многократном воздействии. Таким образом, состояние среды в текущий момент времени оказывается зависимым от суммы возмущений во все предшествующие моменты. Это наиболее общее предположение о локальной реакции линейной среды на возмущение, может быть выражено в виде интегрального уравнения наследственного типа (уравнения Вольтерра) [Работнов 1979; Ландау, Лифшиц Т.VIII 1992] с релаксационными ядрами, описывающими дисперсионные свойства среды.

Вид релаксационных ядер в рамках наследственной теории упругости подбирается, как правило, из эвристических соображений. Такие ядра обычно должны быть положительными, монотонно убывать с ростом времени, удовлетворять критерию «затухающей памяти» [Работнов 1979; Локшин, Суворова 1982], а также аппроксимировать по возможности большее количество экспериментальных зависимостей и быть удобными для интегрирования. В литературе рассматривается ряд модельных релаксационных ядер: экспоненциальное, ядро Абеля, дробно-экспоненциальное ядро Работнова и др. (см. например [Коган 1966; Работнов 1979; Локшин, Суворова 1982, Nigul 1983]). Параметры таких ядер обычно не поддаются явной физической интерпретации. Хотя следует отметить, что определенные соображения в пользу выбора той или иной формы ядра в некоторых случаях могут быть приведены. Механические модели, таким образом, можно рассматривать как метод получения экспоненциальных релаксационных ядер.

Как отмечено в работе [Работнов 1979], способ построения релаксационных ядер, основанный на построении механических реологических моделей, обладает тем преимуществом, что релаксационные ядра построенные на его основе не противоречат законам термодинамики, хотя, по мнению автора [Работнов 1979], «было бы достаточно наивно искать в структуре материала соответствующие упругие и вязкие элементы».

Однако, при всех своих достоинствах подход, основанный на наследственной теории упругости, также является в определенной мере феноменологическим, опосредованным образом учитывающим реальные процессы происходящие в среде при деформировании на молекулярном или микроструктурном уровне. В рамках такого подхода априори трудно определить могут ли свойства конкретной среды быть описаны данным релаксационным ядром. И хотя параметры релаксационных ядер представляют собой некоторые характеристики среды, прямая связь этих параметров с ее микроструктурными характеристиками, как правило затруднена.

Заметим также, что в электродинамике общая линейная связь между напряженностью и электрической индукцией или поляризацией (см. например [Ландау, Лифшиц, Т.VIII 1992; Гинзбург, 1967; Виноградова и др. 1990; Ермаченко 1998] и др.) также описывается уравнением наследственного типа (уравнением Вольтерра). Таким образом, математически задача о распространении электромагнитного импульса в линейной среде, дисперсия в которой определяется уравнением состояния наследственного типа, также оказывается эквивалентной задаче о распространении акустического импульса в среде с уравнением состояния, описываемым в рамках наследственной теории упругости. ([Кельберт, Сазонов 1987; Виноградова и др. 1990] и ДР-)

Микроскопические модели релаксационных механизмов

Качественно иной подход к описанию дисперсионных свойств среды основан на моделировании процессов, происходящих в среде на молекулярном или микроструктурном уровне под влиянием макроскопического внешнего воздействия.

В рамках этого подхода моделируются элементарные процессы, происходящие с микроструктурным элементом среды при макроскопическом воздействии, после чего результат тем или иным образом усредняется по объему среды.

Дисперсионно-диссипативные свойства многих сред были эффективно описаны в рамках этого подхода. В многоатомных газах дисперсия, как показано Кнезером [Kneser 1913] и др., определяются релаксационными процессами передачи энергии между внешними (поступательными) и внутренними (вращательными и колебательными) степенями свободы молекулы. Такие процессы обычно называют кнезеровскими (см. например [Красильников, Крылов 1984]). В газе возможны и другие релаксационные процессы ([Михайлов, Соловьев, Сырников 1964; Красильников, Крылов 1984; Мезон 1969]), в частности, процессы установления химического равновесия и др. Разнообразные релаксационные процессы возможны в жидкостях и в твердых телах ([Красильников, Крылов 1984; Михайлов, Соловьев, Сырников 1964; Кельберт, Сазонов 1991; Новик, Берри 1975]). В частности в жидкостях также возможны кнезеровские релаксационные процессы, различного рода химическая релаксация (например процессы диссоциации молекул электролита), структурная релаксация, связанная с изменением ближнего порядка в расположении молекул, релаксация пузырьков в жидкостях с пузырьками газа, а также другие релаксационные процессы, причем многие из них часто происходят одновременно. Различные релаксационные процессы как на молекулярном, так и на микроструктурном уровне возможны в жидких кристаллах [Капустин, Капустина 1986], стеклах и полимерах [Михайлов, Соловьев, Сырников 1964]. В твердых кристаллических телах ([Новик, Берри 1975; Михайлов, Соловьев, Сырников 1964; Кожевников 1997; Мамин 2001; Ерофеев, Ромашов 2002]) возможны релаксационные процессы, связанные с дислокациями в кристаллах, фазовыми переходами и др.

Однако, если для многоатомных газов как правило удается построить простые и достаточно адекватные микроскопические модели релаксационных процессов, то для жидкостей и твердых тел микроскопическое описание процессов релаксации часто представляет собой очень сложную задачу. ([Михайлов, Соловьев, Сырников 1964; Красильников, Крылов 1984] и др.)

При этом важно также отметить, что во многих жидких и твердых средах (например в полимерах, вязких жидкостях, жидких кристаллах и др.) одновременно может протекать большое количество как различных так и однотипных релаксационных процессов с широким спектром времен релаксации. ([Михайлов, Соловьев, Сырников 1964; Красильников, Крылов 1984; Исакович, Чабан 1988]).

Отметим также, что частотная дисперсия фазовой скорости электромагнитных волн в различных средах может быть интерпретирована в рамках микромоделей релаксационных процессов на молекулярном уровне. Наиболее известными в электродинамике моделями являются модель Дебая полярных диэлектриков [Debye 1929], и модель Лоренца неполярных диэлектриков [Lorentz, 1952].

Модель Лоренца описывает релаксации поляризации молекул неполярного диэлектрика на основе уравнений динамики отдельных связанных электронов [Lorentz, 1952; Гинзбург 1967; Ермаченко 1998]. В этой модели предполагается, что движение электрона может быть описано уравнением гармонического осциллятора с затуханием.

Такое же уравнение описывает релаксацию в жидкости с пузырьками газа или в кристаллах с дислокациями ([Новик, Берри 1975; Накоряков, Покусаев, Шрейбер 1983; Красильников, Крылов 1984; Буланов 2001]).

Дисперсия в полярных диэлектриках была описана Дебаем ([Debye 1929]), как релаксационный процесс теплового разупорядочивания ориентации дипольных молекул, упорядоченных электромагнитным импульсом. В этой модели считается, что тепловая релаксация поляризации происходит по экспоненциальному закону ([Виноградова и др. 1990]), аналогично средам с релаксацией кнезеровского типа. В результате частотные зависимости фазовой скорости и коэффициента поглощения в модели Дебая имеют вид аналогичный тем, что возникают в многоатомных газах вследствие кнезеровских релаксационных процессов.

Иногда модель Дебая используют как феноменологическую модель, аналогично механическим моделям,, рассмотренным выше. Так, например, для аппроксимации экспериментальных данных по диэлектрической проницаемости мускульной ткани использовалась модель Дебая с двумя и пятью эмпирически подобранными временами релаксации. [Hurt 1985].

Адекватные микроскопические модели по-видимому наиболее полно описывают дисперсионные свойства соответствующих сред. В тоже время, как оказалось, разные микроскопические модели приводят в линейном приближении к математически эквивалентным уравнениям, описывающим процессы релаксации (например, отмеченные выше различные по своей природе релаксационные процессы кнезеровского типа и др.). Более того, характеристики релаксационных процессов в линейном приближении как правило не зависят от деталей той или иной модели, а определяются такими комбинациями их микроскопических параметров, которые могут быть интерпретированы, например, как характерные времена соответствующих релаксационных процессов.

Термодинамический подход

Определенной альтернативой микроскопическим моделям является подход, предложенный Мандельштамом и Леонтовичем [Мандельштам, Леонтович 1937]. Он с одной стороны учитывает реальные релаксационные процессы происходящие в среде при возмущении, а с другой позволяет не вникать в микроскопические механизмы этих процессов.

Этот подход основан на описании отклика среды на возмущение в линейном приближении квазиравновесной термодинамики. В этом случае предполагается что состояние термодинамического равновесия среды полностью определяется ее основными термодинамическими переменными (например давлением, плотностью и температурой). Внешнее возмущение (распространяющийся импульс) вьшодит среду из состояния термодинамического равновесия. При этом релаксация внутренних параметров происходит значительно медленнее чем изменение основных термодинамических переменных и можно считать, что среда в каждый момент времени находится в состоянии неполного равновесия, характеризуемого помимо основных термодинамических переменных, еще и текущими значениями внутренних параметров.

Примерами таких внутренних параметров могут быть степень диссоциации молекул, распределение энергии между внутренними и внешними степенями свободы молекул многоатомного газа, концентрация дислокаций и др.

Релаксация среды в целом описывается как совокупность различных релаксационных процессов, каждый из которых, независимо от его физического механизма, характеризуется двумя параметрами: временем релаксации и мощностью, выражаемыми через термодинамические характеристики среды. В этом случае уравнение состояния является уравнением наследственного типа с релаксационным ядром в виде суммы затухающих экспонент, каждая из которых описывает отдельный релаксационный процесс. Это ядро аналитически эквивалентно тому, которое получается из механических моделей, обобщающих стандартное неупругое тело, а также из многих микроскопических моделей. В частности, релаксационное ядро такого типа возникает в модели Кнезера релаксационных процессов в многоатомных газах, поэтому все релаксационные процессы, описываемые экспоненциальными ядрами часто называют кнезеровскими ([Кельберт, Сазонов 1987]).

Термодинамический поход Леонтовича и Мандельштама в настоящее время является общепринятым при описании множества релаксационных процессов в различных средах. Изложению этого подхода и его применению для описания свойств многих сред уделено значительное место в ряде фундаментальных монографий. В частности, в работе [Михайлов, Соловьев, Сырников 1964] проведено описание на основе термодинамического подхода Мандельштама и Леонтовича основных типов молекулярных и микроструктурных релаксационных процессов возможных в жидкостях и газах, а в монографии [Новик, Берри 1975] приведена термодинамическая интерпретация разнообразных релаксаций в кристаллических твердых телах.

Таким образом, в рамках термодинамического подхода Мандельштама и Леонтовича удается описать дисперсионные свойства различных сред в терминах релаксационных процессов, эффективно на макроскопическом уровне учитывающих проявление микроструктуры среды, не углубляясь при этом в микроскопические механизмы этих процессов.

В фундаментальных монографиях [Михайлов, Соловьев, Сырников 1964;Новик, Берри 1975] даже утверждается, основываясь на известных в то время экспериментальных работах, что в рамках подхода Леонтовича Мандельштама могут быть по-видимому описаны практически все релаксационные процессы по крайней мере в газах и жидкостях. Исключение, по мнению этих авторов, составляли только релаксационные механизмы, связанные с дислокациями в кристаллах.

Проблема единого описания релаксационных и резонансных сред в рамках термодинамического подхода

Однако, при всей общности подхода Леонтовича-Манделыптама, существуют среды релаксация которых не описывается суперпозицией экспоненциальных релаксационных процессов. Это, во-первых, микронеоднородные среды, в которых дисперсия является не только частотной, но и пространственной [Исакович 1979; Кельберт, Чабан 1986], а также среды с резонансной релаксацией. В качестве примера можно привести жидкости с пузырьками газа [Накоряков, Покусаев, Шрейбер 1983; Красильников, Крылов 1984; Бескаравайный, Ковалев, Поздеев 1983, Буланов 2001], кристаллы с дислокациями [Новик, Берри 1975, Михайлов, Соловьев, Сырников 1964, Ерофеев, Ромашов 2002], твердые стекла при низкой температуре, системы связанных частиц со спином, помещенных в магнитное поле и др. [Исакович, Чабан 1988], а также, неполярные диэлектрики [Lorentz, 1952] (для электромагнитных волн).

Обычно среды с резонансной релаксацией рассматриваются отдельно от общего термодинамического подхода Мандельштама-Леонтовича. Например, в работе [Кельберт, Чабан 1986] предложено независимо рассматривать три типа релаксации кнезеровскую, резонансную и релаксацию диффузного обмена.

Если ограничиться только средами с локальным откликом на внешнее возмущение (т.е. средами без пространственной дисперсии), то все линейные среды можно разбить на два класса: среды с экспоненциальной релаксацией, свойства которых описывается в рамках подхода Леонтовича-Мандельштама и среды с резонансной релаксацией. К последним относится, в частности, модель Лоренца [Lorentz, 1952], которая по существу является одним из примеров микромоделей, не сводящихся к термодинамическому описанию.

Однако из общих соображений неясно, почему релаксация в резонансных средах не может в линейном приближении быть описана в рамках квазиравновесной термодинамики.

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

Наличие же достаточно общего описания дисперсионно-диссипативных свойств сред, в свою очередь, позволяет в общем виде поставить как прямую задачу - задачу об определении закономерностей распространения импульсов в проивольной линейной среде, так и обратную - задачу об определении параметров такой среды по динамике формы распространяющегося импульса, т.е. задачу импульсной диагностики сред.

Распространение импульсов в диспергирующих средах

Для решения как прямой, так и обратной задачи важно иметь аналитические соотношения, описывающие динамику профиля импульса в пространственно-временной области. Несмотря на то, что в рамках линейной теории импульс во временной области дается интегралом Фурье от его спектральных компонент, выражения для которых сравнительно нетрудно получить аналитически, соотношения между параметрами среды, координатой и временем в пространственно-временном представлении, описывающим динамику импульса, оказываются нетривиальными и представляют основной интерес при решении как прямой, так и обратной задачи.

Для квазимонохроматических импульсов (импульсов огибающей), такие выражения могут быть получены в рамках классической теории дисперсии (см. например [Вайнштейн 1976; Виноградова, и др. 1990] др.). Однако дисперсия квазимонохроматических импульсов определяется главным образом, частотной зависимостью фазовой скорости и коэффициента поглощения лишь вблизи несущей частоты. Короткие же импульсы содержат широкий спектр частот и соответствующая теория для них оказывается неприменимой.

Когда же длительность импульса оказывается сравнимой по величине с характерными временами релаксационных процессов получение явно зависящих от времени аналитических выражений для профиля импульса, вообще говоря, представляет собой весьма сложную задачу. Для решения этой задачи используются как точные, так и различные приближенные методы.

Асимптотические методы

Традиционным способом получения пространственно-временного представления импульса является применение асимптотических методов, в частности, метода перевала. По-видимому, впервые для анализа распространения синусоидальной волны с передним фронтом по среде, описываемой моделью Лоренца, этот метод был применен Зоммерфельдом [Sommerfeld 1914] и Бриллюеном [Brillouin 1914]. Последовательное изложение этого ставшего теперь классическим анализа имеется в монографии Бриллюена [Brillouin 1960], где воспроизведены также и оригинальные статьи. Важные уточнения в методику Зоммерфельда и Бриллюена были внесены Бирвалдом [Bearwald 1930]. Методика теоретического исследования импульсов в диспергирующих средах предложенная в этих работах активно развивается ив настоящее время. [Oughstun, Sherman 1988; Oughstun, Sherman 1989; Shen, Oughstun 1989; Oughstun, Sherman 1990; Oughstun, Laurens 1991; Oughstun, Sherman 1994; Oughstun 1995; Oughstun, Balictsis 1996; Oughstun, Balictsis 1997; Xiao, Oughstun 1998, Barakat, Baumann 1969; Кельберт, Чабан 1986; Кельберт, Сазонов 1988; Кельберт, Сазонов 1991]

Современная техника использования метода перевала для асимптотического анализа динамики импульса в среде Лоренца развита в серии работ Остена, Шермана и их коллег: [Oughstun, Sherman 1988; Oughstun, Sherman 1989; Shen, Oughstun 1989; Oughstun, Sherman 1990; Oughstun, Laurens 1991; Oughstun, Sherman 1994; Oughstun 1995; Oughstun, Balictsis 1996; Oughstun, Balictsis 1997; Xiao, Oughstun 1998]. По сравнению с классическими работами Зоммерфельда и Бриллюена в этих работах основное внимание уделяется более точной аппроксимации зависимости перевальной точки как от пространственной и временной переменных, так и от параметров среды. Кроме того для получения равномерной асимптотики во всей области параметров используется прифронтовое разложение. В целом авторам указанных работ удается получить более точное воспроизведение динамики импульса в области его диспергирования.

В работах [Barakat, Baumann 1969; Кельберт, Чабан 1986; Кельберт, Сазонов 1988; Кельберт, Сазонов 1991] метод перевала использовался для асимптотического анализа распространения импульсов различной начальной формы в среде с одним кнезеровским релаксационным процессом.

Несколько иная идея вычисления обратных преобразований Лапласа или Фурье заключается в разложении подынтегральной функции тем или иным образом в асимптотический ряд с дальнейшим интегрированием нескольких первых членов этого ряда. Этот подход в том или ином виде использовался в работах [Trizna, Weber 1982; Varoquaux, Williams, Avenel 1986; Wyns, Foty, Oughstun 1989; He, Storm 1996; Karlsson, Rikte 1998; Xiao, Oughstun 1998]. Как правило разложение проводилось при частоте стремящейся к бесконечности, и, таким образом, первые проинтегрированные члены ряда описывали прифронтовую часть импульса, определяемую высокочастотными компонентами спектра. Этим обстоятельством ограничивается применимость такого подхода. В работе [Karlsson, Rikte 1998] помимо этого использовалось еще и низкочастотное разложение, что позволило правильно описать относительно низкочастотную часть импульса (предвестник Бриллюена).

Еще один асимптотический подход заключается в замене на больших расстояниях от источника исходного интегро-дифференциального уравнения уравнением более простого диффузионного типа [Roberts, Petropoulos 1996; Roberts, Petropoulos 1999]. Функция Грина для такого уравнения также может быть получена в явном виде и использована на больших расстояниях от источника. Результаты полученные в этом подходе согласуются с результатами, полученными методом перевала. [Roberts, Petropoulos 1999]

Все упомянутые выше работы посвящены распространению импульса в средах с одним или, в некоторых случаях, двумя релаксационными процессами ([Xiao, Oughstun 1998]). Однако, как уже отмечалось, в реальных средах может одновременно протекать множество релаксационных процессов с широким спектром времен релаксации (СВР). Таким образом, для описания распространения импульсов в реальных средах необходимо разработка подходов, учитывающих наличие широких и даже непрерывных СВР.

Ввиду отсутствия до последнего времени общего подхода к проблеме описания распространения импульса в среде с произвольным СВР в ряде работ предпринимались попытки получить описание динамики различных частей импульса в различных асимптотических пределах. Так в работах [Blake 1974; Leander 1991] и [Leander 1993] для среды Максвелла с произвольным СВР получены прифронтовые асимптотики. Асимптотики больших времен и расстояний для импульсов в среде с произвольным СВР рассматривались в работах [Blake 1974; Дунин, Максимов МТТ 1988].

Однако отдельные асимптотики не позволяют описать эволюцию импульса в целом для среды с произвольным СВР, особенно в области его диспергирования. А такая постановка вопроса является весьма актуальной для задач импульсной акустодиагностики релаксационных сред.

В работе [Дунин, Максимов 1988] предложен подход к аналитической аппроксимации профиля короткого импульса в среде с произвольным набором релаксационных процессов, пригодный для всей области распространения. Он основан на построении выражения, переходящего в соответствующих пределах в высокочастотную (прифронтовую) асимптотику и в низкочастотную асимптотику, соответствующую большим расстояниям от источника и от фронта импульса в среде с произвольным набором релаксационных процессов. Коэффициенты этой аналитической аппроксимации выражены через три первых момента СВР < г* >,к = 0,1,2,3.

Однако в работе [Максимов 1996] показано, что информация о параметрах спектра времен релаксации среды (о пяти моментах СВР < г* >,к = -2,-1, 0,1, 2 ) содержится в особенностях профиля импульса как вблизи фронта, исчезающих на больших расстояниях, так и вдали от него. Поэтому важно иметь равномерное описание полного профиля импульса во всей области его диспергирования в терминах указанных параметров СВР. Использование такой аппроксимации позволяет в практической плоскости ставить задачу определения моментов СВР по экспериментально измеренной динамике профиля импульса.

Точные решения

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

Однако до настоящего времени удалось получить только два точных решения для функции Грина короткого импульса для среды с одним экспоненциальным релаксационным процессом (модель стандартного неупругого тела) (см. например [Morrison 1956; Вайнштейн 1976; Дунин 1986], предельными случаями которого являются среды Кельвина-Фойгта [Зверев 1950; Carpenter 1967] и Максвелла [Berry 1958], а также для среды СВР вида 1/г [Дунин, Максимов 1990], выраженную через гипергеометрическую функцию в случае малой дисперсии фазовой скорости. Отметим также, что функции Грина для сред с Е и Ei памятью, полученные в [Нигул 1983], оказываются эквивалентными двум предыдущим случаям. На этом список известных автору точных решений исчерпывается.

Вместе с тем значительный интерес представляла бы функция Грина для среды с двумя релаксационными процессами. Помимо теоретической важности этой задачи для понимания динамики импульса в средах с распределенным СВР, простейшим примером которых является такая среда, эта задача имеет и практическое значение, в частности, для акустики океана, так как в морской воде дисперсия и поглощение, связаны главным, образом с двумя релаксационными процессами при диссоциации солей MgS04 и В(ОН)з. (см. например [Житковский 1995])

Столь же важным является и поиск функции Грина для среды Лоренца с резонансной релаксацией, что позволило бы снять вопрос о равномерном описании динамики короткого импульса в такой среде, возникающий в различных асимптотических подходах [Oughstun, Sherman 1994].

Таким образом, для понимания закономерностей распространения коротких импульсов в релаксационных средах важным является поиск новых аналитических фундаментальных решений в пространственно-временной области.

Распространение импульсов в неоднородных средах

Следует заметить, что на практике создавать одномерные импульсы часто бывает достаточно сложно, и соответствующее описание оказывается приближенным. В тоже время, довольно обычными являются источники, излучающие сферические и цилиндрические волны. Закономерности распространения сферических и цилиндрических импульсов в однородной релаксационной среде рассматривались в работах [Кукуджанов 1963; Blake 1974; Кельберт, Сазонов 1987; Дунин, Максимов МТТ 1988; Кельберт, Сазонов 1991; Рохлин 1995], однако, даже в этом более простом случае найдены лишь асимптотики отдельных частей профиля импульса, а точных решений до последнего времени получено не было.

В силу этого представляет определенный интерес выяснить, как сказывается неодномерность геометрии излучения на динамике изменения профиля импульса, распространяющегося в неоднородной релаксационной среде.

Все упомянутые выше работы, посвященные распространению импульсов, касались только однородных сред. В неоднородных средах помимо частотной дисперсии, связанной с микроструктурой среды, возникает также и пространственная дисперсия, определяемая неоднородностью. Таким образом, как для прямой задачи так и для задачи импульсной акустической диагностики среды на практике важно уметь описьюать совместное влияние на форму распространяющегося импульса диспергирующих свойств среды, связанных с релаксацией и с пространственной неоднородностью. Такая задача может быть актуальна, например, в случае температурно-неоднородной среды, или при моделировании динамики импульсов, распространяющихся в газовой атмосфере, неоднородность плотности в которой определяется полем тяжести.

Идея общего подхода к описанию распространения импульсов в неоднородной среде с дисперсионно-диссипативными свойствами была высказана в работе [Шемякин 1955] и эффективно развита в работах [Дунин, Максимов МТТ 1988; Дунин, Максимов 1990; Максимов 1993; Максимов 1994]. Суть подхода заключается в разделении (факторизации) пространственно неоднородных свойств среды и ее дисперсионно-диссипативных свойств. Если такое разделение оказывается возможным, то решение задачи сводится к свертке решения упругой неоднородной задачи и фундаментального решения плоской задачи для дисперсионно-диссипативной среды.

Этот же подход использован в работе [Максимов 1994] для получения точных решений, описывающих распространение плоского импульса в экспоненциально неоднородной среде, с релаксаций описываемой моделью Максвелла. Представляет интерес обобщить решение этой задачи с учетом реальной неодномерности распространения импульса в атмосфере. Кроме того, интересным представляется оценить влияние температурных неоднородностей в океане по методике описанной в работе [Максимов 1993] на динамику распространения коротких импульсов с целью их дистанционной диагностики импульсными методами.

Экспериментальные результаты

Реальный смысл различные теоретические построения, в том числе и касающиеся распространения импульсов, приобретают тогда, когда они могут быть сопоставлены с экспериментом. В этой связи следует упомянуть несколько экспериментальных работ, которые подтверждают теоретические результаты.

Аналитическое выражение функции Грина во временной области было с успехом использовано для интерпретации результатов экспериментов по распространению коротких импульсов в уксусной кислоте [Андреев, Сапожников, Тимофеев 1994; Зенкова, Зозуля 2001].

Здесь следует заметить, что для интерпретации более ранней экспериментальной работы [Carome, Parks, Mraz 1964] по распространению импульсов конечной длительности в уксусной и пропановой кислоте, терахлориде и дисульфиде углерода использовалось численное суммирование рядов Фурье через которые был выражен источник [Carome, Fleury, Wagner 1964]. А в работе [Moffett, Beyer 1970] использовался подход [Blackstock 1967], основанный на классической теории дисперсии, для интерпретации результатов по распространению конечного импульса с узким спектром в терахлориде углерода. В обоих подходах удалось получить приемлемое согласие с экспериментальными профилями, но оба они, в отличии от [Андреев, Сапожников, Тимофеев 1994; Зенкова, Зозуля 2001] не позволяют получить из экспериментальных профилей важные параметры среды, такие, например, как время релаксации.

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

Цель работы

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

Описание распространения коротких импульсов в однородных средах, обладающих спектром времен релаксации в терминах экспериментально измеряемых параметров СВР.

Поиск точной функции Грина для среды с двумя релаксационными процессами.

Поиск точной функции Грина точечного и линейного источников в изотермической атмосфере с релаксационными свойствами Максвелла.

Единое описание релаксационных и резонансных свойств сред в рамках термодинамического подхода Леонтовича-Манделыптама

Поиск точной функции Грина, описывающей динамику коротких импульсов в средах с резонансной релаксацией и исследование закономерностей распространения короткого импульса в такой среде.

Научная новизна

Научная новизна приводимых в диссертации результатов состоит в том, что впервые получены точные пространственно временные представления функций Грина ряда задач, а именно для среды с двумя релаксационными процессами, среды с единственным процессом резонансной релаксации, изотермической атмосферы с релаксационными свойствами Максвелла для точечного и линейного источников. Для среды с произвольным СВР построена новая аналитическая аппроксимация профиля импульса с использованием экспериментально измеряемых параметров его динамики. Впервые дано обобщение термодинамического подхода Мандельштама-Леонтовича при учете альтернативной формулировки принципа симметрии кинетических коэффициентов Онзагера. На этой основе впервые выведено универсальное уравнение состояния линейных сред с локальным откликом и дана его механическая интерпретация. Для среды с одним процессом резонансной релаксации проведено исследование возможных типов динамики короткого импульса и впервые дана их полная классификация.

Положения выносимые на защиту

На защиту выносятся следующие результаты:

Новая аналитическая аппроксимация пространственно-временной функции Грина одномерной среды с произвольным спектром времен релаксаций (СВР), использующая экспериментально измеряемые моменты СВР.

Новое точное пространственно-временное представление функции Грина одномерной среды с двумя экспоненциальными релаксационными процессами, его анализ и упрощенная аналитическая аппроксимация.

Новые точные пространственно-временные представления функций Грина линейного и точечного источников в изотермической атмосфере с релаксационными свойствами Максвелла и их анализ.

Новое универсальное уравнение состояния сред с линейным локальным откликом, выведенное путем обобщения термодинамического подхода Манделыптама-Леонтовича при учете альтернативной формулировки принципа симметрии кинетических коэффициентов Онзагера, и его механическая интерпретация.

Новое точное пространственно-временное представление функции Грина одномерной среды с одним процессом резонансной релаксации, описываемым в рамках обобщенного уравнения состояния.

Новая полная классификация допустимых форм эволюции короткого импульса, распространяющегося в однородной среде с одним процессом резонансной релаксации, в рамках обобщенного уравнения состояния сред с линейным локальным откликом.

Научная и практическая значимость

Полученные результаты имеют значение как для моделирования динамики профиля импульса в релаксационных средах так и для решения задач импульсной диагностики.

Новая аппроксимация временного представления функции Грина для сред с СВР может быть использована для моделирования распространения акустического импульса в релаксационных средах со сложной реологией. При этом зависимость параметров аппроксимации только от экспериментально определяемых моментов СВР позволяет использовать ее для целей импульсной диагностики сред.

Эта аппроксимация также может быть применена для моделирования динамики коротких импульсов в пространственно неоднородных релаксационных средах, там где получение точных решений представляет принципиальную сложность, а численные методы расчета не достаточно эффективны.

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

Новое точное временное представление функции Грина для сред с двумя экспоненциальными релаксационными процессами также может быть использовано для моделирования динамики импульса и для определения релаксационных параметров сред (времен релаксации и относительных мощностей релаксационных процессов). Это точное временное представление позволяет подробно исследовать эволюцию профиля импульса на расстояниях, сравнимых с дисперсионной длиной. Кроме того, это точное представление позволяет оценить эффективность аппроксимации функции Грина для СВР.

Точные временные представления функций Грина точечного и линейного источника в экспоненциально неоднородной среде с релаксацией, описываемой моделью Максвелла, демонстрируют одновременное влияние пространственной и частотной дисперсии, а также геометрии излучения на динамику формы распространяющегося импульса. В частности показано, что в зависимости от соотношения между релаксационными свойствами среды и пространственной дисперсией, связанной с неоднородностью, динамика импульса может быть либо релаксационной, либо дисперсионной. Более того, при определенных соотношениях параметров, дисперсия, обусловленная неоднородностью среды и дисперсия, связанная с ее релаксационными свойствами могут полностью компенсировать друг друга. Эти временные представления функций Грина могут быть использованы для моделирования распространения импульсов в простой модельной изотермической атмосфере.

Теоретическое значение обобщенного уравнения состояния заключается в том, что им полностью исчерпывается описание линейного локального отклика сред в состояниях близких к термодинамическому равновесию. Это позволяет, помимо прочего, ставить задачу построения общей теории распространения импульсов малой амплитуды в таких средах. В свою очередь на основе такой теории потенциально возникает возможность ставить задачу импульсной диагностики сред с линейным локальным откликом в общем виде.

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

Новое аналитическое представление во временной области функции Грина плоского источника в среде с одним процессом резонансной релаксации, описываемым в рамках обобщенного уравнения состояния, может быть использовано для моделирования распространения импульса в среде с процессом резонансной релаксации наиболее общего вида. В аналитической структуре этого представления удалось как отдельные слагаемые выделить ранее описываемые лишь асимптотически предвестники Зоммерфельда и Бриллюена и описать их формирование на расстояниях сравнимых с дисперсионной длиной.

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

Апробация

Результаты работы докладывались на IV, X, XI сессиях Российского акустического общества (1995,1999,2001); Fourth International Congress on Sound and Vibration (1996); III Международной научно-технической конференции «Современные методы и средства океанологических измерений» (1997); 16 International Congress on Acoustics and 135^th Meeting of Acoustical Society of America (1999); Третьем совещании по магнитной и плазменной аэродинамике в аэрокосмических приложениях (2001); 141 Meeting of Acoustical Society of America (2001); International Photonic Research 2001; семинарах «Акустика неоднородных сред» научной школы проф. Рыбака (2000, 2002); научных сессиях МИФИ 1999,2000,2001,2003.

Публикации

По теме диссертации в научных журналах и трудах конференций опубликовано 26 работ [Maksimov, Larichev 1996; Ларичев, Максимов 1997-1, Ларичев, Максимов 1997-2, Ларичев, Максимов 1997-3, Ларичев, Максимов 1998; Maksimov, Larichev 1998; Larichev, Maksimov 1998; Ларичев, Максимов 1999-1, Ларичев, Максимов 1999-2,

Ларичев, Максимов 1999-3; Larichev, Maksimov 1999-1; Larichev, Maksimov 1999-2; Larichev, Ларичев, Максимов 2000-1; Ларичев, Максимов 2000-2; Ларичев, Максимов 2000-3; Ларичев, Максимов 2001-1; Ларичев, Максимов 2001-2; Ларичев, Максимов 2001-3; Maksimov 2001-1; Larichev, Maksimov 2001-2; Ларичев, Максимов 2002; Ларичев, Максимов 2003-1; Ларичев, Максимов 2003-2] в том числе 7 работ в ведущих научных журналах по данной тематике.

Структура и объем

Структурно диссертация состоит из введения, пяти глав и заключения. Работа содержит 164 страниц, в том числе 115 стр. текста, 54 рисунка, 8 страниц библиографии.

Содержание

В первой главе исследуются общие закономерности распространения коротких импульсов малой амплитуды в средах, совокупность релаксационных процессов в которых при этом может быть описана в рамках линейной неравновесной термодинамики на основе подхода Мандельштама и Леонтовича ([Мандельштам, Леонтович 1937]).

В первом разделе излагается термодинамический подход Мандельштама и Леонтовича для сред с дискретным и непрерывным спектрами времен релаксаций (СВР). В рамках этого подхода воспроизведен вывод известных уравнений состояния сред вблизи термодинамического равновесия (см. например [Новик, Берри 1975; Михайлов, Соловьев, Сырников 1964]).

Во втором разделе из уравнений движения и непрерывности, а также уравнения состояния, вывод которых приведен в первом разделе, получено наследственное уравнение, описывающее распространение волн в средах с СВР кнезеровского типа. Приведено спектральное представление функции Грина плоского источника в данной среде (одномерной задачи Коши для данного уравнения). Приведены высокочастотные и низкочастотные временные асимптотические представления функции Грина. С их помощью показаны некоторые важные закономерности динамики импульса и их связь с моментами СВР ([Максимов 1996]). Также обосновано приближение малой дисперсии для большинства реальных сред (см. например [Новик, Берри 1975; Михайлов, Соловьев, Сырников 1964]).

В третьем разделе, с использованием подхода, предложенного в работе [Дунин, Максимов 1988], получена новая аналитическая аппроксимация временного представления функции Грина задачи Коши сред с произвольным СВР, зависящая только от явно определяемых из профиля импульса моментов СВР. Аппроксимирующее выражение воспроизводит характерные измеримые особенности профиля импульса, характеризуемые этими моментами. Кроме того, для среды с одним релаксационным процессом полученная аппроксимация полностью аналитически воспроизводит точное решение. Однако, в силу своего приближенного характера, данная аппроксимация содержит некоторые ограничения на параметры СВР, связанные с требованием положительности аргументов входящих в нее радикалов. Условия применимости также явно сформулированы.

В четвертом разделе приведены численные расчеты профиля короткого импульса, распространяющегося в средах с различными СВР, выполненные полученной аппроксимации. Также для среды с двумя релаксационными процессами приведено сравнение профилей тела импульса, вычисленных с помощью аппроксимационного выражения и с помощью точного временного представления функции Грина, полученного во второй главе.

Во второй главе получены новые временные представления функции Грина задачи Коши для среды с двумя релаксационными процессами.

В первом разделе второй главы, с использованием теоремы Эфроса об обобщенной свертке, получено временное представление функции Грина, состоящее из трех слагаемых, описывающих соответственно упругий предвестник, высокочастотную и низкочастотную часть тела импульса. Первые два из них выражены через элементарные функции и функции Бесселя, а третье представляет полностью действительный интеграл в конечных пределах, прямой численный расчет которого также не представляет сложностей.

Детальный анализ функции Грина позволил записать ее высокочастотную часть в инвариантном виде, справедливом для среды с произвольным СВР. При этом особенности высокочастотной части тела импульса содержат информацию о первых моментах СВР, определяемых экспериментально.

Во втором разделе, для низкочастотной составляющей тела импульса получено аналитическое выражение, которое позволяет получить ее адекватное описание в широком диапазоне изменения параметров СВР. Справедливость этого выражения подтверждена путем сравнения с результатами прямого численного расчета.

В третьем разделе приведено альтернативное временное представление функции Грина, полученное с помощью обычной теоремы о свертке. Оно состоит из четырех слагаемых, первое из которых описывает упругий предвестник, два других высокочастотные вклады от каждого из релаксационных процессов в отдельности, а четвертое слагаемое представляет собой низкочастотный совокупный вклад от обоих релаксационных процессов.

В третей главе на ряде характерных примеров исследовано влияние пространственной неоднородности релаксационной среды на характер распространения импульсов.

В первом разделе получены точные выражения для функций Грина линейного и точечного источников в неоднородной среде с экспоненциальным распределением плотности вдоль одной из координат для случая, когда дисперсионно диссипативные свойства среды описываются в рамках модели Максвелла. Проанализированы особенности динамики распространения акустического импульса, возбуждаемого в такой среде линейным и точечным источниками.

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

С использованием аппроксимирующего аналитического выражения для функции Грина задачи о распространении импульса в однородной среде, на основании подхода, предложенного в работе [Максимов 1993], проведены оценки влияния неоднородности среды на характер распространения импульса.

В первом разделе четвертой главы дано обобщение термодинамического подхода Мандельштама и Леонтовича, которое приводит к описанию отклика среды в терминах экспоненциальных релаксаций, на случай резонансных релаксаций. В основе этого обобщения лежит альтернативная формулировка принципа симметрии кинетических коэффициентов Онзагера для величин, обладающих различной симметрией по отношению к инверсии времени. С учетом этого обстоятельства выводится обобщенное локальное уравнение состояния среды, находящейся вблизи термодинамического равновесия. Показано, что данное уравнение состояния является исчерпывающим при описании линейного отклика среды в рамках локальных моделей. Показано, что все существующие в настоящее время локальные модели линейного отклика среды, являются частными случаями обобщенного уравнения состояния.

Во втором разделе на примере единственного механизма резонансной релаксации продемонстрировано многообразие дисперсионно-диссипативных свойств, которые описываются в рамках обобщенного уравнения состояния. В этом частном случае получены все возможные типы частотной зависимости фазовой скорости и коэффициента поглощения, и показано, каким образом характер этой зависимости определяется параметрами обобщенного уравнения состояния. Определены множества значений параметров, определяющих различные типы частотной зависимости фазовой скорости и коэффициента поглощения.

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

В пятой главе проведено исследование динамики импульса в среде с единственным процессом резонансной релаксации.

В первом разделе получено новое аналитическое представление функции Грина плоского импульса в среде с одним процессом резонансной релаксации. Это представление является полностью действительным. В нем аналитически выделен упругий предвестник и предвестник Зоммерфельда, а предвестник Бриллюена выражен интегралом от действительной функции в конечных пределах. С помощью этого представления возможен эффективный прямой численный расчет профиля импульса на не слишком больших расстояниях от источника.

Во втором разделе проведен подробный асимптотический анализ пространственно временного представления функции Грина, полученного в первом разделе, на больших расстояниях от источника, что позволило рассчитывать временные профили импульса во всей области распространения.

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

В Заключении сформулированы основные выводы диссертации.

Аппроксимирующее аналитическое выражение, описывающее динамику короткого импульса, распространяющегося в релаксационной среде с произвольным СВР

Полученные результаты имеют значение как для моделирования динамики профиля импульса в релаксационных средах так и для решения задач импульсной диагностики.

Новая аппроксимация временного представления функции Грина для сред с СВР может быть использована для моделирования распространения акустического импульса в релаксационных средах со сложной реологией. При этом зависимость параметров аппроксимации только от экспериментально определяемых моментов СВР позволяет использовать ее для целей импульсной диагностики сред.

Эта аппроксимация также может быть применена для моделирования динамики коротких импульсов в пространственно неоднородных релаксационных средах, там где получение точных решений представляет принципиальную сложность, а численные методы расчета не достаточно эффективны.

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

Новое точное временное представление функции Грина для сред с двумя экспоненциальными релаксационными процессами также может быть использовано для моделирования динамики импульса и для определения релаксационных параметров сред (времен релаксации и относительных мощностей релаксационных процессов). Это точное временное представление позволяет подробно исследовать эволюцию профиля импульса на расстояниях, сравнимых с дисперсионной длиной. Кроме того, это точное представление позволяет оценить эффективность аппроксимации функции Грина для СВР.

Точные временные представления функций Грина точечного и линейного источника в экспоненциально неоднородной среде с релаксацией, описываемой моделью Максвелла, демонстрируют одновременное влияние пространственной и частотной дисперсии, а также геометрии излучения на динамику формы распространяющегося импульса. В частности показано, что в зависимости от соотношения между релаксационными свойствами среды и пространственной дисперсией, связанной с неоднородностью, динамика импульса может быть либо релаксационной, либо дисперсионной. Более того, при определенных соотношениях параметров, дисперсия, обусловленная неоднородностью среды и дисперсия, связанная с ее релаксационными свойствами могут полностью компенсировать друг друга. Эти временные представления функций Грина могут быть использованы для моделирования распространения импульсов в простой модельной изотермической атмосфере.

Теоретическое значение обобщенного уравнения состояния заключается в том, что им полностью исчерпывается описание линейного локального отклика сред в состояниях близких к термодинамическому равновесию. Это позволяет, помимо прочего, ставить задачу построения общей теории распространения импульсов малой амплитуды в таких средах. В свою очередь на основе такой теории потенциально возникает возможность ставить задачу импульсной диагностики сред с линейным локальным откликом в общем виде.

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

Новое аналитическое представление во временной области функции Грина плоского источника в среде с одним процессом резонансной релаксации, описываемым в рамках обобщенного уравнения состояния, может быть использовано для моделирования распространения импульса в среде с процессом резонансной релаксации наиболее общего вида. В аналитической структуре этого представления удалось как отдельные слагаемые выделить ранее описываемые лишь асимптотически предвестники Зоммерфельда и Бриллюена и описать их формирование на расстояниях сравнимых с дисперсионной длиной.

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

Результаты работы докладывались на IV, X, XI сессиях Российского акустического общества (1995,1999,2001); Fourth International Congress on Sound and Vibration (1996); III Международной научно-технической конференции «Современные методы и средства океанологических измерений» (1997); 16 International Congress on Acoustics and 135th Meeting of Acoustical Society of America (1999); Третьем совещании по магнитной и плазменной аэродинамике в аэрокосмических приложениях (2001); 141 Meeting of Acoustical Society of America (2001); International Photonic Research 2001; семинарах «Акустика неоднородных сред» научной школы проф. Рыбака (2000, 2002); научных сессиях МИФИ 1999,2000,2001,2003.

По теме диссертации в научных журналах и трудах конференций опубликовано 26 работ [Maksimov, Larichev 1996; Ларичев, Максимов 1997-1, Ларичев, Максимов 1997-2, Ларичев, Максимов 1997-3, Ларичев, Максимов 1998; Maksimov, Larichev 1998; Larichev, Maksimov 1998; Ларичев, Максимов 1999-1, Ларичев, Максимов 1999-2, Ларичев, Максимов 1999-3; Larichev, Maksimov 1999-1; Larichev, Maksimov 1999-2; Larichev, Ларичев, Максимов 2000-1; Ларичев, Максимов 2000-2; Ларичев, Максимов 2000-3; Ларичев, Максимов 2001-1; Ларичев, Максимов 2001-2; Ларичев, Максимов 2001-3; Maksimov 2001-1; Larichev, Maksimov 2001-2; Ларичев, Максимов 2002; Ларичев, Максимов 2003-1; Ларичев, Максимов 2003-2] в том числе 7 работ в ведущих научных журналах по данной тематике.

Структурно диссертация состоит из введения, пяти глав и заключения. Работа содержит 164 страниц, в том числе 115 стр. текста, 54 рисунка, 8 страниц библиографии.

В первой главе исследуются общие закономерности распространения коротких импульсов малой амплитуды в средах, совокупность релаксационных процессов в которых при этом может быть описана в рамках линейной неравновесной термодинамики на основе подхода Мандельштама и Леонтовича ([Мандельштам, Леонтович 1937]).

В первом разделе излагается термодинамический подход Мандельштама и Леонтовича для сред с дискретным и непрерывным спектрами времен релаксаций (СВР). В рамках этого подхода воспроизведен вывод известных уравнений состояния сред вблизи термодинамического равновесия (см. например [Новик, Берри 1975; Михайлов, Соловьев, Сырников 1964]).

Во втором разделе из уравнений движения и непрерывности, а также уравнения состояния, вывод которых приведен в первом разделе, получено наследственное уравнение, описывающее распространение волн в средах с СВР кнезеровского типа. Приведено спектральное представление функции Грина плоского источника в данной среде (одномерной задачи Коши для данного уравнения). Приведены высокочастотные и низкочастотные временные асимптотические представления функции Грина. С их помощью показаны некоторые важные закономерности динамики импульса и их связь с моментами СВР ([Максимов 1996]). Также обосновано приближение малой дисперсии для большинства реальных сред (см. например [Новик, Берри 1975; Михайлов, Соловьев, Сырников 1964]).

Аналитическое выражение для функции Грина среды с двумя релаксационными процессами

При описании динамики распространения акустических импульсов в реальных средах возникает задача одновременного учета, как релаксационных свойств среды, так и ее пространственной неоднородности.

В первой главе было получена аппроксимация, позволяющая описывать динамику импульса в однородной релаксационной среде. Но с помощью подхода предложенного в работе [Максимов 1993], аппроксимация (1.3.15) может быть использована для описания эволюции профиля одномерного акустического импульса при его распространении в неоднородной среде, когда характерная длинна волны импульса много меньше масштаба пространственной неоднородности среды. В этом случае применимо квазиклассическое описание [Бреховских 1973], которое, однако, не позволяет учитывать диспергирующие свойства неоднородной среды, проявляющиеся при сопоставимости пространственных масштабов импульса и неоднородности среды. Совместное влияние на распространение импульса диспергирующих свойств неоднородной сред и ее релаксационных свойств было исследовано в работе [Максимов 1994] на примере одномерного импульса, распространяющегося в неоднородной максвелловской среде, для которой были получены точные решения. Однако, на практике создавать одномерные импульсы часто бывает достаточно сложно, и соответствующее описание оказывается приближенным. В тоже время, довольно обычными являются источники, излучающие сферические и цилиндрические волны. В силу этого представляет интерес выяснить, как сказывается неодномерность геометрии излучения на динамике изменения профиля импульса распространяющегося в релаксационной среде. Как уже отмечено во Введении, закономерности распространения сферических и цилиндрических импульсов в однородной релаксационной среде рассматривались в работах [Кукуджанов 1963; Blake 1974; Кельберт, Сазонов 1987; Дунин, Максимов МТТ 1988; Кельберт, Сазонов 1991; Рохлин 1995] однако, даже в этом более простом случае однородной среды точных решений получено не было.

В третей главе на ряде характерных примеров исследовано влияние пространственной неоднородности релаксационной среды на характер распространения импульсов. В первом разделе получены точные выражения для функций Грина линейного и точечного источников в неоднородной среде с экспоненциальным распределением плотности вдоль одной из координат для случая, когда дисперсионно диссипативные свойства среды описываются в рамках модели Максвелла. Проанализированы особенности динамики распространения акустического импульса, возбуждаемого в такой среде линейным и точечным источниками. Во втором разделе, на основании асимптотического подхода, предложенного в [Максимов 1993], исследована эволюция профиля одномерного акустического импульса при его распространении в неоднородной среде, когда характерная длинна волны импульса много меньше масштаба пространственной неоднородности среды. С использованием аппроксимирующего аналитического выражения для функции Грина задачи о распространении импульса в однородной среде (1.3.15) на основании подхода предложенного в работе [Максимов 1993] проведены оценки влияния температурной неоднородности среды на характер распространения импульса в морской воде. В этом разделе будут получены точные выражения для функций Грина модельной задачи для линейного и точечного источников в неоднородной среде, дисперсионно-диссипативные свойства которой описываются моделью Максвелла, а плотность экспоненциально распределена вдоль одной из координат В случае единственного релаксационного процесса, предполагая пространственную зависимость равновесной плотности р0 (R) и обратной адиабатической сжимаемости //(R), локальное уравнение состояния (1.2.3) для изотропной релаксационной среды может быть записано в виде: Заметим, что зависимость волнового числа f(R,/?) от переменной р в выражении (3.1.6) проявляется только как функциональная зависимость от К{р). Это обстоятельство позволяет с помощью теоремы Эфроса об обобщенной свертке [Лаврентьев, Ша-бат 1987] представить функцию Грина задачи для уравнения (3.1.5) в бесконечной среде в виде Здесь Gy,(R,p)- функция Грина уравнения Гельмгольца (3.1.5) когда волновое число К. (R,/?) имеет вид: Необходимо отметить, что представление (3.1.6), в котором разделены факторы, ответственные за пространственную неоднородность задачи и релаксационные свойства среды, оказалось эффективным при анализе большой группы задач (впервые такой подход предложен Е.И. Шемякиным в работе [Шемякин 1955]) [Максимов 1993], [Максимов 1994], [Дунин, Максимов МТТ 1988], [Дунин, Максимов ЖПМТФ 1988], особенно в тех случаях, когда функции I(t,) и CL,(R,f) могут быть найдены аналитически. Так, например, для сферических волн в однородной среде, функция Грина которых при отсутствии релаксации имеет вид Gw (R, І) = S(t — R/c)/R , учет релаксационных про цессов в соответствии с (3.1.8) сразу приводит к выражениюG (R,/) = I(t,R) IR, отличающемуся от плоского случая только множителем 1 / R . Одним из немногих случаев (см. литературу в [Максимов 1994]), для которых существует точное решение для фактора I(t,), является среда Максвелла. Она соответствует предельному случаю релаксационной среды с одним временем релаксации г , когда низкочастотная скорость распространения возмущений стремится к нулю CQ — О (Л- 1). В этом случае волновое число К(р) записывается в виде , а выражение для фактора I(t,g) имеет вид [Максимов 1994] Среди ограниченного набора точных решений уравнения Гельмгольца (3.1.5) [Бре-ховских 1973], [Li, Franke, Liu 1993], [Li 1994] достаточно простое для дальнейших вычислений пространственно-временное представление функции Грина G„,(R,f) линейного и точечного источников существует при экспоненциальном распределении плотности вдоль одной из координат и постоянной скорости распространения. Рассмотрим задачу о распространении импульса давления от линейного (двухмерная задача) и точечного источников в среде с экспоненциальным распределением плотности вдоль одной из координат, релаксационные свойства которой описываются моделью Максвелла. Если положить, что плотность распределена по закону p(R) = р(0)ё аг, где а = mg/ kfrT , m-молекулярная масса воздуха,g- ускорение свободного падения, Т - температура, kb - постоянная Больцмана, a ( (R) = const, то эта задача будет описывать распространение импульсов от линейного и точечного источников в изотермической атмосфере с релаксационными свойствами Максвелла (рис.3.1.1), тем не менее следует подчеркнуть, что, конечно, реальную земную атмосферу эта постановка конечно не описывает и выбрана она лишь потому, что в этом случае удается получить точные решения. Такая задача является обобщением результатов, полученных в работе [Максимов 1994] для плоских импульсов.

Распространение плоского импульса в неоднородной вдоль направления распространения релаксационной среде

В предыдущих главах исследовались закономерности распространения импульсов в средах, которые могут быть описаны в рамках термодинамического подхода Леонтовича-Манделыптама. Этот подход позволяет описать большой класс сред, в которых могут протекать релаксационные процессы различной физической природы с экспоненциальной зависимостью от времени. Как показано в главе 1, все многообразие таких сред описывается единым образом в терминах спектра времен экспоненциальных релаксаций. Но в терминах экспоненциальных релаксационных моделей невозможно описывать акустические среды, обладающие резонансными свойствами. Такими свойствами, например, обладает жидкость с пузырьками газа [Красильников, Крылов 1984, Накоряков, Покусаев, Шрейбер 1983]. Существуют и другие примеры проявления резонансных свойств у акустических сред [Новик, Берри 1975, Красильников, Крылов 1984, Кельберт, Чабан 1986].

Однако описания резонансных и релаксационных сред стоят особняком друг от друга. Действительно, в ряде известных монографий ([Новик, Берри 1975], [Михайлов, Соловьев, Сырников 1964], [Ландау, Лифшиц Т.6 1986]) утверждается, что теория внутренних параметров Мандельштама и Леонтовича приводит к совокупности только экспоненциальных релаксаций. Для описания в рамках той же концепции резонансных свойств среды приходится несколько искусственно предполагать, что динамика релаксации внутренних параметров подчиняется кинетическим уравнениям более высокого порядка, чем первый [Кельберт, Чабан 1986]. В этой связи возникает возможность использовать для описания кинетики внутренних параметров уравнения и более высоких степеней с не очень ясной перспективой к каким особенностям это может привести. Тем более что остается открытым вопрос о взаимодействующей системе внутренних параметров, удовлетворяющих системе таких уравнений. Отметим, что в работе [Кельберт, Чабан 1986] сделана попытка на базе теории внутренних параметров рассмотреть три различных типа релаксаций: экспоненциальную (кнезеровскую), резонансную, и релаксацию, связанную с диффузионным обменом. При этом для описания резонансной релаксации приходится писать кинетические уравнения второго порядка, а для релаксации диффузионного обмена - и вовсе использовать уравнение в частных производных, учитывающее нелокальную природу внутреннего параметра, что, вообще говоря, выходит за границы применимости локальной теории Мандельштама и Леонтовича [Мандельштам, Леонтович 1937]. Поэтому подход работы [Кельберт, Чабан 1986], хотя и демонстрирует принципиальную возможность существования релаксаций различного типа, однако не позволяет последовательным образом вывести общую структуру локального уравнения состояния в рамках теории Мандельштама и Леонтовича [Мандельштам, Леонтович 1937].

В первом разделе четвертой главы дано обобщение термодинамического подхода Мандельштама и Леонтовича, которое приводит к описанию отклика среды в терминах экспоненциальных релаксаций, на случай резонансных релаксаций. В основе этого обобщения лежит альтернативная формулировка принципа симметрии кинетических коэффициентов Онзагера для величин, обладающих различной симметрией по отношению к инверсии времени. С учетом этого обстоятельства выводится обобщенное локальное уравнение состояния среды, находящейся вблизи термодинамического равновесия. Показано, что данное уравнение состояния является исчерпывающим при описании линейного отклика среды в рамках локальных моделей. Показано, что все существующие в настоящее время локальные модели линейного отклика среды, являются частными случаями обобщенного уравнения состояния.

Во втором разделе на примере единственного механизма резонансной релаксации продемонстрировано многообразие дисперсионно-диссипативных свойств, которые описываются в рамках обобщенного уравнения состояния. В этом частном случае получены все возможные типы частотной зависимости фазовой скорости и коэффициента поглощения в средах с локальным откликом, и показано, каким образом характер этой зависимости определяется параметрами обобщенного уравнения состояния. Определены множества значений параметров, определяющих различные типы частотной зависимости фазовой скорости и коэффициента поглощения.

В третьем разделе для получения наглядного представления о реологических свойствах среды, описываемой обобщенным уравнением состояния, разработана обобщенная механическая модель релаксационной среды. Для того чтобы обобщить термодинамический подход на случай резонансных сред заметим в начале, что принцип симметрии кинетических коэффициентов Онзагера в форме (1.1.10) из первой главы применим только для тех пар внутренних параметров lt,j\, которые обладают одинаковой четностью при инверсии времени. Если же различные внутренние параметры обладают разной четностью по отношению к инверсии времени, то для таких величин принцип симметрии Онзагера следует писать в виде [Ландау, Лифшиц Т.5 1995]: Поэтому в общем случае матрица у„ уже не будет симметрической и, соответственно,, положительно определенной. Вследствие этого уже нельзя подобрать такое линейное преобразование, как это было сделано в главе 1, при котором одновременно диагонализуются обе матрицы у и и Za» оставляя вдоль диагонали действительные положительные члены. Однако, поскольку в любом случае матрица yt, определяет скорость прироста энтропии [Ландау, Лифшиц Т.5 1995], то все ее собственные числа должны иметь положительную действительную часть. Поскольку произведение действительных матриц у - и Xij также дает матрицу с действительными коэффициентами, то и характеристический многочлен этой матрицы будет иметь действительные коэффициенты. И, как следствие, корни характеристического многочлена будут либо действительными, либо комплексно сопряженными. Если считать среду диссипативной, то действительные части корней должны быть положительны. Однако матрица линейного преобразования {&}— { i = АнА» которая диагонализует матрицу yikXkf в системе (1.1.8), в этом случае уже не будет чисто действительной, и поэтому такое преобразование внутренних параметров, рассматриваемых как действительные физические величины, не корректно. Допустимыми с физической точки зрения могут быть лишь преобразования с действительными коэффициентами. В рамках только таких преобразований матрица уікХи У5106 не может быть полностью диагонализована, однако, она может быть приведена к квазидиагональному (клеточному) виду. В самом деле, в поле действительных чисел простыми делителями характеристического полинома действительной матрицы могут быть лишь полиномы первого и второго порядков. При этом простым делителям первого порядка будут соответствовать диагональные элементы, а делителям второго порядка - клеточные матрицы размера 2x2. Таким образом, в рассматриваемом случае существует такое линейное преобразование переменных {,} - { Я. Л, которое диагонализует матрицу Xij - Zij -\kXuW )у - ijX{u) и приводит к квазидиагональному виду матрицу yij- ij= ikyu X)ij- Матрица Гу=уікХу также будет иметь квазидиагональный вид с тем же количеством клеточных матриц размера 2x2.

Дисперсия фазовой скорости и частотная зависимость коэффициента поглощения в среде с резонансной релаксацией

С учетом альтернативной формулировки принципа симметрии кинетических коэффициентов Онзагера для величин, имеющих различную четность по отношению к инверсии времени, дано обобщение теории внутренних параметров Мандельштама и Ле-онтовича [Мандельштам, Леонтович 1937] на случай резонансной релаксации.

На этой основе получено обобщенное уравнение состояния среды с учетом как резонансных, так и нерезонансных релаксаций. Показано, что в рамках линейного приближения полученное уравнение состояния является общим для сред с локальным откликом на возмущение, находящихся вблизи состояния термодинамического равновесия.

Для среды с единственным процессом резонансной релаксации исследованы частотная дисперсия фазовой скорости и коэффициента поглощения в зависимости от основных параметров уравнения состояния: резонансной частоты, времени релаксации и угла инерционности. Установлены и описаны все возможные особенности поведения коэффициента поглощения и дисперсии фазовой скорости в диссипативных и активных средах и проведена их классификация. На плоскости Q.T — (р определены области параметров, соответствующие выделенным особенностям.

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

Описание линейного отклика среды в терминах совместного спектра собственных частот и времен релаксаций может служить основой для разработки подходов к решению задачи о распространении импульсов малой амплитуды в произвольной однородной среде с линейным локальным откликом на возмущение. А это, в свою очередь, позволяет надеяться решить обратную задачу о диагностике свойств таких сред импульсными методами. Основными параметрами, характеризующими среду и, соответственно, подлежащими определению, являются спектры собственных частот и времен релаксации. Распространение коротких импульсов в среде с единственным процессом резонансной релаксации

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

Для сред, дисперсия в которых может быть описана в рамках экспоненциальных релаксационных процессов, подходы к решению такой задачи описаны в первых трех главах. Они основаны на асимптотическом анализе функций Грина сред с произвольным СВР ([Максимов 1996]) и известных точных пространственно-временных представлениях функций Грина для среды с одним релаксационным процессом [Зверев 1950; Morrison 1956; Berry 1958; Carpenter 1967; Вайнштейн 1976; Дунин 1986 ] и со спектром вида 1/т [Дунин, Максимов 1990].

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

Как уже отмечалось во Введении, теоретические исследования распространения импульса в среде Лоренца с единственным релаксационным процессом насчитывают уже вековую историю, начиная с работ Зоммерфельда [Sommerfeld 1914] и Бриллюена [Brillouin 1914] и заканчивая работами современных исследователей, в частности работами Остена, Шермана и их коллег [Oughstun, Sherman 1988; Oughstun, Sherman 1989; Shen, Oughstun 1989; Oughstun, Sherman 1990; Oughstun, Laurens 1991; Oughstun, Shennan 1994; Oughstun 1995; Oughstun, Balictsis 1996; Oughstun, Balictsis 1997; Xiao, Oughstun 1998], в которых методы, предложенные Зоммерфельдом и Бриллюеном получили наиболее полное развитие. Подробнее обзор исследований в этой области приведен во введении.

Однако упомянутые работы, в том числе и те, в которых делалась попытка получить пространственно-временное представление функции Грина [Не, Storm 1996; Karlsson, Rikte 1998], по существу ограничиваются асимптотическими описаниями различных частей импульса. Правда следует заметить, что в работах Остена и Шермана [Oughstun, Sherman 1989; Oughstun, Sherman 1990; Oughstun, Sherman 1994] делается попытка получения равномерной асимптотики, но и в этом случае речь идет об асимптотическом описании на больших расстояниях от источника. В тоже время существенными для диагностики среды могут оказаться особенности профиля импульса, которые исчезают на больших расстояниях.

Кроме того, как показано в четвертой главе, среда Лоренца является лишь частным случаем среды с резонансной релаксацией. Поэтому возникает вопрос, какие особенности динамики импульса определяются резонансными средами, отличными от среды Лоренца? В данной главе проведено полное исследование динамики импульса в произвольной среде с локальным откликом при наличии единственного процесса резонансной релаксации.

Наиболее простой средой, описываемой в рамках обобщенного термодинамического подхода, и не сводимой к известным, частным случаям среды с экспоненциальной релаксацией и среды Лоренца является, как уже отмечено выше, среда с одним обобщенным резонансным релаксационным процессом, которая описывается уравнением состояния 4.2.1 Далее будет получено пространственно-временное представление функции Грина такой среды. Функция Грина задачи для уравнения (4.2.1) с ядром (4.2.2) в бесконечной среде может быть представлена в виде интеграла Лапласа.

Однако, прямой численный расчет интеграла (5.1.1) по исходному контуру интегрирования представляет значительные сложности, поскольку подынтегральная функция содержит особенности. Асимптотическим, аналитическим и численным оценкам интеграла (5.1.1) в частном случае (р = к12 посвящено большое число работ: [Sommerfeld 1914; Brillouin 1914; Bearwald 1930; Brillouin 1960; Oughstun, Sherman 1988; Oughstun, Sherman 1989; Shen, Oughstun 1989; Oughstun, Sherman 1990; Oughstun, Laurens 1991; Oughstun, Sherman 1994; Oughstun 1995; Oughstun, Balictsis 1996; Oughstun, Balictsis 1997; Xiao, Oughstun 1998; Trizna, Weber 1982; Varoquaux, Williams, Avenel 1986; Wyns, Foty, Oughstun 1989; He, Storm 1996; Karlsson, Rikte 1998].

Похожие диссертации на Распространение коротких акустических импульсов в средах с экспоненциальной и резонансной релаксацией

Динамика распространения предельно коротких импульсов в оптически одноосных средах Соболевский Анатолий Фёдорович

Динамика распространения предельно коротких импульсов в оптически одноосных средахСоболевский Анатолий Федорович

Исследование спектра и релаксации магнитных возбуждений в проводящих средахХалиуллин, Гиниятулла Гатиятович

Исследование влияния самоподобной структуры среды на физические процессы переноса и релаксацииРябов, Ярослав Евгеньевич

Некоторые вопросы теории релаксации разреженных потоков быстрых заряженных частиц в неоднородных средахГаюров, Хаким Шарифович

К теории спин-решеточной релаксации и диффузного затухания стимулированного спинового эха в полимерных системах и неоднородных средахФаткуллин, Наиль Фидаиевич

Особенности распространения сверхкоротких импульсов лазерного излучения в размерноограниченных полупроводниковых структурахКоровай Олеся Васильевна

Особенности распространения ультракоротких лазерных импульсов и их воздействие на простые квантовые системыМихайлов Евгений Михайлович

РАСПРОСТРАНЕНИЕ КОРОТКИХ АКУСТИЧЕСКИХ ИМПУЛЬСОВ В СРЕДАХ С РЕЛАКСАЦИЕЙ И ОБОБЩЕННЫЙ ВАРИАЦИОННЫЙ ПРИНЦИП ДЛЯ ДИССИПАТИВНОЙ МЕХАНИКИ СПЛОШНЫХ СРЕДМаксимов, Герман Адольфович

Распространение термоупругих волн в среде с учетом релаксации теплового потокаБабенков, Михаил Борисович