Вариационные методы как наиболее эффективный механизм при моделировании взаимосвязанных физических полей в сплошных средах. I. краткий обзор теории Текст научной статьи по специальности «Математика»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2017, том 27, № 1, c. 102-112

МАТЕМАТИЧЕСКИЕ МЕТОДЫ г

И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ

УДК 51-72+51.73+530.1 © Б. П. Шарфарец

ВАРИАЦИОННЫЕ МЕТОДЫ КАК НАИБОЛЕЕ ЭФФЕКТИВНЫЙ МЕХАНИЗМ ПРИ МОДЕЛИРОВАНИИ ВЗАИМОСВЯЗАННЫХ ФИЗИЧЕСКИХ ПОЛЕЙ В СПЛОШНЫХ СРЕДАХ. I. КРАТКИЙ ОБЗОР ТЕОРИИ

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

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

введение

В теории научного приборостроения часто приходится сталкиваться с перекрестным взаимодействием физических полей различной природы (гидродинамических, электромагнитных), различных явлений переноса и т. д., когда действие одного поля на другое порождает обратное воздействие. В случае отсутствия перекрестных взаимодействий полей или их слабого взаимодействия при построении уравнений процесса, как правило, доступно использование разнообразных методов, включая и вариационные. Однако при нетривиальных взаимодействиях полей, которые имеют место, например, для взаимопроникающих сред типа дисперсных смесей или для таких явлений, как акустика пористых сред и т. д., единственным способом построений физически разумных уравнений становится вариационный подход [1, с. 7]. В качестве примеров такого подхода можно привести теоретическое решение Рэлеем задачи Рэ-лея—Плато о капиллярной неустойчивости струй [2] или построение Био взаимосвязанных систем уравнений в теории акустики пористых сред [3, 4] (по этому поводу см. также работу [5], а также обзор [6]). Множество примеров применения вариационного подхода, когда он является единственно возможным, приведено, в частности, в работе [1]. Кроме того, в [1] помещена обширная библиографическая подборка публикаций по рассматриваемому вопросу.

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

[7]). Из длинного ряда работ, затрагивающих применение вариационных подходов к другим областям физики, сошлемся лишь на обзорную работу

[8], классический трактат [9], где вводится понятие диссипативной функции, классические монографии по математической физике [10, 11]. Из узкоспециальных приложений вариационных методов отметим работы [12-14] и еще раз отметим обширный библиографический обзор в [1] по разнообразным приложениям вариационных методов в физике.

постановка задачи

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

предварительные сведения

Говорят, что поведение сплошной среды известно, если известна зависимость от лагранжевых координат и времени некоторой системы величин ик, к = 1,2,..., характеризующих состояние среды и называемых определяющими функциями [1, с. 10]. Построение любой модели сплошной среды начинается с перечисления ее определяющих функций. Так, в частности, в случае движения натянутой гибкой струны такой функцией является ее отклонение от состояния равновесия и = ^( х, t) [11, с. 288], а в случае потенциального движения невязкой жидкости в линейном приближении таковым является потенциал и = р(х,у,z,t) скорости жидкости V(х,у, z,t) , V = Ур [11, с. 294] и т. д.

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

Приведем некоторые общеизвестные определения [1].

Вариацией 8ик (х, t) функции ик (х, t) называют величину

5ик (х, t) =

дик (х,t,s)

де

ds = и dе .

31 = ®

де

dе ,

называемая вариацией функционала 1 .

В вариационных подходах обычно рассматривают интегральные функционалы вида [1, с. 14]

1 (и) = |л(х',ик,ик)с1пх .

(2)

Здесь введено стандартное сокращение для производных функции f ( х1, х2,..., х',...)

/.=д-.

Л' дх'

Согласно (1), 3(ик) = (3ик) . Дифференцирование по е функции 1 (и,е) дает [1, с. 14]

31 = I

дЛ 8ик ип

ди

ди'

3),' )а'

х =

13ЛТ 3ик dnx + Г-^Л Зикп. dn-1 х. \дик -1 дтк '

эу ди,'

(1) Здесь

ЗЛ дЛ д дЛ

С точностью до малых (dе)2 можно написать

ик (х,t,е) = ик (х,t) + Зик (х,t) .

Из определения (1) вытекает, что оператор 3 обладает свойствами оператора дифференцирования:

3( f+g ) = 3( f)+3( g), 3( f • g ) = 3( f )• g+f 3( g).

Если F — заданная функция ик, то

3F = ^3ик. ди

Вариации функционалов. Пусть 1 (и) — некоторый функционал, определенный на множестве М . Берется некоторый элемент и еМ и строится семейство процессов сравнения и(е), и(0) = и , е > 0 . На каждом семействе процессов сравнения и (е) функционал 1 становится функцией е . Вводится величина

5ик дик дх' дик '

п = (п1, п2, п3) — нормаль к границе дУ облас-

ти V . Полагая, что на дУ 5ик =0, имеем

31= 13Л_ 5ик d пх= Г

; зик л

дЛ д дЛ

Л

3и

дик дх' дик

Зи dпх .

Условия стационарности интегрального функционала. Стационарные точки интегрального функционала (2) удовлетворяют [1, с. 15] в области У уравнениям

3Л

37"

= 0.

а на границе области дУ краевым условиям

сл п

—гп = 0.

дик '

(3)

(4)

В случае отсутствия вариаций на границе Зи \ =0 для стационарности функционала (2) достаточно условия (3).

Е = 0

0

вариационные принципы и уравнения

Вариационные принципы

История механики и физики изобилует многочисленными попытками объяснить происходящие явления при помощи небольшого числа универсальных законов и общих принципов. Наиболее удачные попытки связаны с идеей о том, что наблюдаемые явления обладают экстремальными свойствами и искомые общие принципы носят вариационный характер [1, 7, 8 и др.]. Еще Аристотель писал, что природа "идет легчайшим путем, с наименьшими затратами". Первым к физической проблеме этот принцип применил П. Ферма в 1662 г., исследуя закон преломления света на границе двух различных по оптическим свойствам сред. При этом он принял постулат: "Природа действует наиболее легкими и доступными путями. ...Свет выберет такой путь, чтобы время его прохождения оказалось минимальным". Из этого принципа немедленно следует закон Снеллиуса преломления света. В механике первым вариационный принцип сформулировал Мопертюи. Затем усилиями Эйлера, Лагранжа, Гамильтона, Остроградского, Якоби и Пуанкаре этот принцип аналитической механики для консервативных систем с конечным числом степеней свободы был доработан до современного состояния и носит в настоящее время название "принцип наименьшего действия". Интегральный функционал (2) в этом случае имеет вид

Ч

5 = | Ь (I) dt, (2а)

где 5 — действие; Ь ^) = Ь (q ^), q ({), {) —

функция Лагранжа, равная т. н. кинетическому потенциалу Ь = К — П , К и П — соответственно кинетическая и потенциальная энергии системы; q и q — совокупность обобщенных координат и их производных по времени. Уравнение (3) в этом случае имеет вид

дЬ — оЬ = 0 дqk dt дqк = .

Принцип стационарного действия в форме Гамильтона утверждает, что истинная траектория механической системы является стационарной точкой функционала (2а) на множестве траекторий, выходящих в момент времени ^ из точки q1 и попадающих в момент времени t2 в точку q2 [1, с. 23]. Существует также эквивалентный ему принцип стационарного действия в форме Гамильтона—Пуанкаре, когда действие (2а) имеет вид

5 = {( РД' — Н (p, q )) dt, (2б)

tl

где p = ( р1, р2,...) — импульс системы; Н (р, q) — функция Гамильтона системы, определяемая как преобразование Лежандра функции Ь по скоростям [1, с. 95] и равная для консервативных систем сумме кинетической и потенциальной энергий. Стационарные точки функционалов (2а) и (2б) совпадают.

Накопленные в физике результаты показывают [1, с. 24, 33], что вариационные принципы справедливы для всех фундаментальных физических полей, а также для микроявлений, которые можно считать обратимыми. Все элементарные взаимодействия частиц, все поля (электромагнитные, гравитационное, сильные и слабые внутриядерные взаимодействия) подчиняются законам, следующим из условий стационарности некоторого действия. Более того, убеждение в справедливости этого факта стало "общим местом" настолько, что в новых ситуациях, когда закономерности изучаемого явления еще неизвестны, предлагается заранее действие, минимизацией которого находят искомые уравнения. Так был получен ряд уравнений квантовой механики, усовершенствованы уравнения теории относительности. Поэтому в настоящее время утверждается, что все известные фундаментальные физические поля удовлетворяют принципу стационарности действия [1, с. 24].

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

Вариационные уравнения

Исторически первым вариационным уравнением был

принцип возможных перемещений ("золотое правило механики") [1, с. 28].

В современном виде его сформулировал Бер-нулли. Здесь приводится формулировка для системы материальных точек, подчиненных некоторым кинематическим идеальным (т. е. при отсутствии потерь) связям. Пусть х5 — координата 5-й

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

I Г,-Зх, =0.

(5)

При отсутствии кинематических ограничений (Зх, = 0) уравнение (5) тривиально: Г, = 0 .

Принцип Даламбера и уравнение энергии перемещений. В задачах динамики ВУ (5) будет справедливым, если к силам Г, добавить силы инерции — т, а,

1(Г, "т,а,)Зх, = 0.

(6)

Здесь , а, — масса и ускорение ,-й частицы соответственно. ВУ (6) называется принципом Да-ламбера. В случае отсутствия связей из (6) следуют уравнения Ньютона.

ВУ (6) (принцип Даламбера) считают исходным постулатом механики систем с конечным числом степеней свободы.

В работе [15] было показано, что вариационное уравнение механики (6) есть уравнение энергии, записанное для возможных перемещений, которое записывается так:

ЗК = 3Ае + ЗО . (7)

Здесь ЗАе =1 Г, -Зх, — работа внешних сил;

К = I — кинетическая энергия системы;

, 2

функция З О определяется так:

ЗО = -? (т тгЗх, -тv ).

Функционал З О на действительном движении обращается в нуль, поэтому на действительном движении вариационное уравнение (7) переходит в уравнение энергии

ЗК = ЗАе.

ВУ (7) можно рассматривать как запись первого начала термодинамики (уравнения энергии) для возможных перемещений в случае систем с конечным числом степеней свободы [1, с. 30]. Отметим схожесть последнего уравнения и уравнения Лагранжа второго рода (см., например, [16, т. 2, с. 542]).

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

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

П (() , по которым строится функция Лагранжа L ((, () . При переходе к сплошным средам кинетическая К и внутренняя Е энергии становятся функционалами от определяющих функций. Поскольку сплошная среда описывает статистические закономерности для большого числа хаотически двигающихся частиц, возникает такая характеристика, как энтропия S, которая также является функционалом от определяющих функций [1, с. 25]. Большое значение для получения результатов по данной тематике оказали работы Дж. Гиб-бса по термодинамике гетерогенных систем [17], а также работы Л.И. Седова (см. библиографию в [1], под номерами [212-219]).

Предполагается, что энтропия и внутренняя энергия аддитивны по массе:

£ = | р, dv, Е = | ре dv,

У У

а кинетическая энергия по объему:

К = ,

У

где р — плотность массы; dv — элемент объема области У; , и е — плотности энтропии и внутренней энергии на единицу массы соответственно, а к — плотность кинетической энергии в единице объема

1 I |2

к=—рV .

2 1

В процессе при произвольных приращениях З (а не действительных приращениях d) уравнение первого начала термодинамики таково [1, с. 30]:

ЗЕ = ЗАе + ЗФ + ЗО.

(8)

где Е — энергия системы; ЗАе — приращение работы внешних макроскопических сил; З Ф — приток тепла и всех остальных видов энергии; З О — невязка между действительным процессом и случаем произвольного допустимого приращения (в действительном процессе З О 0 ).

Если функционалы Е, ЗАе, ЗФ и ЗО опре-

делены, то уравнение (8) превращается в вариационное уравнение, выражающее первое начало термодинамики для возможных приращений определяющих функций. Функционал 3Е есть вариация

функционала Е. Под 3 Ае понимают значения функционала dАе на значениях 3ик .

Приток тепла 3Ф на любом допустимом процессе определяется так [1, с. 31]:

3Ф = | рТ3я dv — 3Ф',

где 3Ф' — некомпенсированное тепло; Т — абсолютная температура. После подстановки значения 3Ф в (8) получается вариационное уравнение первого начала термодинамики

3Е = 3Ае рТ3яdv — 3Ф'+ 30.

(9)

Отметим, что при выводе (9) использовалось второе начало термодинамики [1, с. 31].

Вариационное уравнение Седова [1, 15]. В ВУ (9) не все величины можно задавать независимо. Часть задается явно, часть определяется. Вводятся следующие выражения:

3Е = { рГ 2у|2 + И dv,

30 = |р^3 — ^3х. |dv .

2

зJJлdvdt + 3Ж * + 3Ж = 0

(10)

где

1 V

л=р| || уГ—е I,

3Ж = —

ч

¡ру13х1 dv +|3а; dt

к ^

3Ж* = Ц |рТ3яdv — 3Ф'+ 3А,

1 V V

dt.

Вариационное уравнение Седова (10) имеет ме-

сто для любого объема V и любого интервала времени t t2 ]. Величина Л называется лагранжианом (плотностью функции Лагранжа, см., например, [11, гл. 3]). Лагранжиан Л и функционал 3Ж * есть задаваемые величины, а функционал 3Ж находят из вариационного уравнения (10) [1, с. 32].

ВУ (10) имеет две отличительные черты. Во-первых, оно записано для любой части сплошной среды (это сближает его по форме с уравнением энергии (9)); во-вторых, оно содержит вклады, связанные с необратимыми процессами. Построение новых моделей в рамках вариационного подхода заключается в фиксировании набора определяющих функций и задании Л и 3Ж*. Вычисление 3Ж соответствует установлению уравнений состояния.

Вариационное уравнение Лагранжа. Это

уравнение используется тогда, когда V — весь объем, занятый сплошной средой, а t1 и t2 — моменты времени, в которые значения определяющих функций заданы. Тогда вариационное уравнение Седова (10) для всего объема перепишется в виде [1, с. 33]

3( I) + 3А = 0,

(11)

где

dt dt

Работа внешних сил Ае представляется в виде суммы работ объемных и поверхностных сил 3Ае = 3А1 + 3А^ . Здесь приняты обозначения V = (v1,v2,Vз) , х = (х1,х2,х3) . После ряда несложных преобразований ВУ (9) переписывается в следующем виде [1, с. 31]:

2 _ _ _

I = JJлdvdt, А = 3Ж*+ 3Ж . (11а)

1 V

Черта сверху говорит о принадлежности величины всему объему V.

Вариационное уравнение (11) обычно называют уравнением Лагранжа. Из уравнения Лагранжа следует замкнутая система уравнений сплошной среды, а также и краевые условия [1, с. 33].

Отметим, что вариационные уравнения первого начала термодинамики Седова и Лагранжа справедливы для обратимых и необратимых процессов.

Условимся, что при упоминании ниже вариационных уравнений вариации функционалов находятся в левой стороне уравнения, приравниваемого к нулю справа, как например в (10), (11).

В некоторых случаях (см. ниже) вариационное уравнение сводится к виду

3(1) = 0

(12)

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

е

ненное описание процессов, в которых отсутствует "необратимая стохастичность", можно дать при помощи непосредственного осреднения микроскопического действия [1, с. 33], [11, с. 289]. Отсюда соответствующие вариационные принципы для различного рода обратимых процессов в теории сплошных сред.

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

и , то система уравнений вариационного типа записывается так (см., например, [1, с. 33], [11, с. 266] и др.):

дЛ д дЛ I д дЛ =0

к = 1, п . (13)

Здесь п — число определяющих функций (полевых переменных); т — пространственная раз-

к

и

ди дх

лономного вариационного уравнения является уравнение (12), где вариация интегрального функционала, будучи равной нулю, превращает вариационное уравнение в вариационный принцип.

Для увеличения прозрачности приведем рассуждения для конечномерного случая с последующим обобщением для бесконечномерного. В конечномерном случае левая часть ВУ представляет собой линейную дифференциальную форму вида Fk (и)Зик , где и = {ик}е Rn. Обозначим ее (по повторяющемуся индексу проводится суммирование) через

ЗО = Fk (и)Зи

(14)

Ясно, что (14) будет вариацией функционала О, если вектор {Fk (и)] потенциален, т. е. существует функция р такая, что

ч (и )=р,

ди,.

к = 1, п .

(15)

к дик

мерность задачи; и г =-

, дt '

Такую форму имеют уравнения электродинамики, гравитации, уравнения механики идеальной жидкости, теории упругости. Во всех этих теориях система (12) полностью определяется заданием только лагранжиана Л .

Важной особенностью явлений, описываемых уравнениями (13), — взаимность физических эффектов. Если рассматривается перекрестное взаимодействие между двумя полями, то действие одного поля на другое автоматически порождает обратное в некотором смысле симметричное воздействие: например, если лагранжиан Л содержит перекрестный член и1^2, то в уравнении для поля и1 появится сила со стороны поля 2, равная и2 и

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

Уравнениям, которые учитывают необратимые процессы, по-видимому, также свойственна специальная структура, однако универсальные обобщения вариационных принципов на необратимые процессы пока неизвестны [1, с. 34].

Условия голономности вариационных уравнений

Вариационное уравнение называется голоном-ным, если его левая часть представляет вариацию некоторого функционала [1, с. 67]. Примером го-

Пусть функции Fk (и) непрерывны и дифференцируемы в области А е Rn. Пусть Зик и З}ик — два бесконечно малых поля в А. Тогда, если обозначить через З'ЗО вариацию формы ЗО вдоль поля З}ик, а через ЗЗ'О — вариацию формы З' О вдоль поля Зик, то можно сформулировать следующую теорему [1, с. 68].

Для того чтобы форма ЗО представляла в области А вариацию некоторой функции, необходимо и достаточно, чтобы во всех точках этой области было выполнено равенство

З'ЗО = ЗЗ' О. (16)

Отсюда следуют хорошо известные условия го-лономности дифференциальных форм (условия

потенциальности вектора (и )}) [1, с. 68 ]

дFk (и )_дFm (и )

ди„

ди,,

к, т = 1, п .

(17)

Условия голономности (16) имеют локальный характер, т. е. их выполнение в окрестности некоторой точки гарантирует существование в этой окрестности функции р(и), такой что ЗО = З (р) . Однозначную функцию р можно указать, если область А односвязна [1, с. 68].

Приведенные рассуждения справедливы и в случае бесконечномерных пространств (в частности, в случае теории сплошных сред), и в качестве критерия голономности также справедливо условие (16) [18, 19].

к

Для иллюстрации приведем пример из [1, с. 69, 70], рассмотренный для функционала

30 = | Е ( х, и, и, и. )3ud"

(18)

где х = {х,}е Rn, и, = и,, и. = и ... Условие (16), (17) для функционала (18) имеет вид

дЕ —дЕ = 0

ди,. и,.,.

(19)

Л(^ у ) = К (у,, ) — и (у,х ) =

=1 р|

ду ( х, t)

дt

— Т

1 п

ду ( х, t)

дх

(20)

Здесь у(х, t) — отклонение от состояния равновесия элемента струны с координатой х во времени /; К и и — линейные плотности соответственно кинетической и потенциальной энергий струны (на единицу ее длины); р — линейная

плотность струны; Т0 — натяжение струны. Найдем 3Л(ул,у х):

3Л(ул у ) = 3У,г 3¥,х =

Легко проверяется, что при Е = 3Л / 3и , где Л = Л( х, и, и ,) — лагранжиан, зависящий от определяющей функции и и ее первых производных, условие (19) выполняется тождественно. Замечание. В Физической энциклопедии [16, т. 1, с. 514] в статье о голономных системах определение голономности дано в менее общем виде, чем приведенное выше, а именно: "В аналитической механике под голономной системой понимается механическая система, в которой все наложенные связи являются геометрическими (голономными). Эти связи налагают ограничения только на возможные положения точек и тел системы, но не на их скорости. Разделение систем на голономные и неголономные весьма существенно, т. к. к голо-номным системам применим принцип наименьшего действия, который несправедлив для неголо-номных систем".

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

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

Натянутая гибкая идеальная струна. При описании этого простейшего примера в [11, с. 288-290] продемонстрирован уже упоминавшийся выше прием усреднения при использовании вариационного подхода в сплошных средах на основе методов аналитической механики дискретных частиц.

Процесс колебаний гибкой струны обратим (консервативен). Поэтому в качестве лагранжиана в [11, с. 289] для случая малых колебаний выбран кинетический потенциал

дЛ д

д¥,х

дЛ д _ 3у +---3у.

ду ( дt дух дх

(21)

д д

Выше учтено, что операторы — и — коммути-

дt дх

руют с оператором 3 . После вариации интеграла

ч ч

I = jjлdvdt получаем 31 = jjзЛdvdt, где 3Л с

^ V t1V

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

(

3Л = —

д дЛ д дЛ +

Л

дх дух дt ду 1

Определим величину 33Л : 3у

3у .

3Л

д дЛ д дЛ

3у дх ду х дt ду { Записываем условие стационарности действия:

3Л

д дЛ д дЛ

3у дх ду х дt ду (

=0

что приводит к уравнению колебания струны [11, с. 290]

С2°у—=0, с 2=Т0

р

дх2 дГ

Исследуем голономность функционала 3Л(уt,у х) . Согласно (14) и (21), имеем

30 = 3Л(уд ,у х )

^ ^ (и)3и1 + Е2 (и)3и2 = 3У,( + 3у,х .

у ду,х

Отсюда следует

х

2

2

2

Fi =

ал

dW,t

F2 =

ал

dW,t

8щ = ôy

Su2 = ôy x .

(22)

Выпишем условие (17) применительно к случаю (22):

а2 л

а2л

= 0.

aV,taV,x av,ßv,t

Отсюда следует, что для вариации лагранжиана (20) условия голономности соблюдаются и мы имеем дело с вариационным принципом, т. е. поиском стационарного значения интегрального функционала.

выводы

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

приложение

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

Диссипативная функция (ДФ)

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

энергии упорядоченного движения в энергию неупорядоченного движения (в конечном счете в тепловую), имеет размерность мощности (и равна потере энергии в единицу времени). ДФ может быть построена для механических систем, у которых скорости макроскопических движений малы настолько, что силы сопротивления движению можно считать линейно зависящими от скоростей [16, т. 1, с. 653], [20, с. 187]. ДФ ввел впервые Рэ-лей [9, гл. 4].

ДФ в случае классической механики является квадратичной формой от значений компонентов скоростей [21, с. 102] или тензоров скорости деформации в случае механики сплошных сред, и, например, для случая движения вязкой теплопроводной жидкости с постоянными коэффициентами вязкости ДФ, отнесенная к единице объема, записывается следующим выражением [13, с. 205] (см. также [16, т. 1, с. 653], [20, с. 187], [22, с. 423] и др.):

Ф = 2/Б: D + Лвг.

Здесь Л, / — параметры Ламе для вязкой изотропной жидкости; Б — тензор скоростей дефор-

маций, D =

s s s

xx xy

s s s

yx yy

s s s

а компоненты тензора

равны

s =

syy =

s„ =

av

_X

ax avy

ay avz az

= s„„ = —

( av ay

avy ï

X + y

1 (avv

S = S = —

yz zy 2

s = s = —

zx xz 2

az av.

-+

ax avz

ay av

+

ax az

в = У-V; V = (vx,vy,vz) — вектор скорости течения жидкости; Б: Б — скалярное произведение двух одинаковых тензоров Б (определение скалярного произведения см. [13, с. 9])

D : D

avx

ax

2 (avyï

+

ay

+(®ü

l az

1

+— 2

^av ay

x + avy ï

2

ax

+— 2

1 ( avy av

-+-az ay

+

1

+— 2

avz

ax

+

avx

az

Окончательно имеем

1

2

2

2

2

Ф = м\2

5x

+

( 5vy V ( 5v_

+

v

+

5y 5x

+

\

5y 5z 5y

+

^ 5z

+

(5v. 5vv Y (5Vy +_5v^2 ^ x21

+ 1 I 5vx

5x 5z

+

+Лв2.

В терминах сдвиговой г/ и объемной д вязко-стей, для которых справедливы соотношения с параметрами Ламе (см., например, [23]) , 2

д = л + — ц, / = ц, последнее выражение для ДФ перепишется в виде

Ф = ц< 2

5v„ Y (5vv Y (5v, л2

5x

+

5v

+ 1 5z

+

+

( 5v 5vv Y (5vv 5v ^ —- +—-

5z 5-

- + -5v 5x

+

(5vz 5vx + 1 — + — I К ^ 5x 5z

+1,-,

(П1)

что полностью совпадает, например, с выражением в работе [24, с. 342].

Представим Ф в переменных £кк :

Ф = 2 Ul X 4+1 [ff - \г,]в

л

(П2)

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

Если для данной среды ДФ известна, то учесть влияние внутреннего трения можно посредством вычисления компонент "диссипативного" тензора стл по формуле ([16, т. 1, с. 653] и [20, с. 188])

' дФ

Тогда результирующий тензор находится как сумма

где в2 = (s2 + s2 + s2 + 2s s + 2s s + 2s s ) ,

^ ^ xx yy zz xx yy xx zz yy zz J '

индексы 1, 2, 3 отвечают соответственно переменным x, y и z. Выражение (П1) для ДФ с точностью до множителя 2 совпадает с соответствующим выражением, представленным в работах [16, т. 1, с. 653] и [20, с. 187]. Аналогичное выражение для ДФ, представленное в работе [25, с. 256], отличается от выражения (П1) только отсутствием в нем коэффициента объемной вязкости g вслед-

2

ствие принятия гипотезы Стокса Л = ——¡и [25,

2

с. 67], что эквивалентно тому, что Л+ — ¡и = 0 или

2

g = Л + — и = 0.

3

Таким образом, Ф из выражения (П2) является квадратичной формой компонент тензора скоростей деформаций с коэффициентами, характеризующими вязкость среды, и представляет собой то

+

где ай — недиссипативная составляющая результирующего тензора.

список литературы

1. Бердичевский В.Л. Вариационные принципы механики сплошной среды. М.: Наука, 1983. 448 с.

2. Rayleigh, Lord. On the Instability of jets // Proceedings of the London mathematical Society. 1878. Vol. 10. P. 4-13.

3. Biot M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range // J. Acoust. Soc. Am. 1956. Vol. 28. P. 168-178. Doi: 10.1121/1.1908239.

4. Biot M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range // J. Acoust. Soc. Am. 1956. Vol. 28. P. 179-191. Doi: 10.1121/1.1908241.

5. Johnson D.L. Recent developments in the acoustic properties of porous media // Frontiers in Physical Acoustics. XCIII. North Holland, 1986. P. 255-290.

6. Князьков Н.Н., Шарфарец Б.П. Акустика пористо-упругих насыщенных жидкостью сред (обзор теории Био) // Научное приборостроение. 2016. Т. 26, № 1. С. 77-84. URL: http://213.170.69.26/mag/2016/full1/ Art11.pdf.

7. Вариационные принципы механики / Сб. статей под ред. Л.С. Полака. М.: Физматгиз, 1959. 952 с.

8. Полак Л.С. Вариационные принципы механики. Их развитие и применение в физике. М.: URSS, 2010. 600 с.

9. Лорд Рэлей. Теория звука. Т. 1. М.: ГИТТЛ, 1955. 504 с.

10. Курант Р., Гильберт Д. Методы математической физики. Т. 1. М.: ГТТИ, 1933. 525 с.

11. Морс Ф.М., Фешбах Г. Методы теоретической физики. Т. 1. М.: Изд-во иностр. лит-ры, 1958. 930 с.

12. Био М. Вариационные принципы в теории теплообмена. М.: Энергия, 1975. 209 с.

13. Серрин Дж. Математические основы классической механики жидкости. М.: Иностранная литература, 1963. 256 с.

2

2

2

2

14. Pierce A.D. Variational formulations in acoustic radiation and scattering // Physical Acoustics. Vol. XXII. San Diego: Acad. Press, 1993. P. 195-394.

15. Седов Л.И. Математические методы построения новых моделей сплошных сред // УМН. 1965. Т. 20, № 5. С. 121-180.

16. Физическая энциклопедия. В 5 томах / Гл. ред. А.М. Прохоров. М.: Советская энциклопедия, 19881998. 704+704+672+704+760 с.

17. Гиббс Дж. Термодинамические работы. М.-Л.: Гос-техиздат, 1950. 492 с.

18. Вайнберг М.М. Вариационные методы исследования нелинейных операторов. М.: Гостехиздат, 1956. 344 с.

19. Вайнберг М.М. Вариационный метод и метод монотонных операторов. М.: Наука, 1972. 416 с.

20. Энциклопедия математической физики / Гл. ред. Л.Д. Фаддеев. М.: Большая российская энциклопедия, 1998. 691 с.

21. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 1. Механика. М.: Наука, 1988. 216 с.

22. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. VI. Гидродинамика. М.: Наука, 1988. 736 с.

23. Шарфарец Б.П., Шарфарец Е.Б., Князьков Н.Н., Па-шовкин Т.Н. Некоторые особенности численного ре-

шения задач термоупругости и гидродинамики тепло-проводящей сжимаемой вязкой жидкости с помощью универсальных пакетов // Научное приборостроение. 2016. Т. 26, № 3. С. 57-63. URL: http://213.170.69.26/ mag/2016/fuU3/Art8.pdf

24. Яворский Б.М., Детлаф А.А., Лебедев А.К. Справочник по физике. Для инженеров и студентов вузов. М.: Оникс: Мир и образование, 2006. 1056 с.

25. Шлихтинг Г. Теория пограничного слоя. М.: Наука, 1974. 712 с.

Институт аналитического приборостроения РАН, г. Санкт-Петербург

Контакты: Шарфарец Борис Пинкусович, sharb@mail.ru

Материал поступил в редакцию 15.12.2016

ISSN 0868-5886

NAUCHNOE PRIBOROSTROENIE, 2017, Vol. 27, No. 1, pp. 102-112

VARIATIONAL METHODS AS THE MOST EFFECTIVE MECHANISM FOR MODELING PHYSICAL INTERRELATED FIELDS IN CONTINUOUS MEDIUM. I. OVERVIEW OF THE THEORY

B. P. Sharfarets

Institute for Analytical Instrumentation of RAS, Saint-Petersburg, Russia

In order to solve the problem of choosing the most effective mechanism for mathematical modeling of Interrelated physical fields in continuous media are considered using the variational approach. It summarizes the current theory on the subject, based on the use of variational principles, or variational equations. Specifies the criteria for their selection, discusses the features of this method. Is provided specific example. In the second part of the work will be done the analysis of the cases, when the variational approach is only possible to produce mathematical models of multiphysical phenomena. The appendix describes the dissipative function (dissipative potential) in the case of a homogeneous viscous heat-conducting fluid, which plays an important role in the modeling of irreversible processes.

Keywords: mechanics of continua, related physical fields, variational principle, the variational equation, terms of holonomic of variational equations, dissipative potential

references

1. Berdichevskij V.L. Variacionnye principy mekhaniki sploshnoj sredy [Variation principles of mechanics of the continuous environment]. Moscow, Nauka Publ., 1983. 448 p. (In Russ.).

2. Rayleigh, Lord. On the Instability of jets. Proceedings of the London mathematical Society, 1878, vol. 10, pp. 4-13.

3. Biot M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range. J. Acoust. Soc. Am., 1956, vol. 28, pp. 168-178. Doi: 10.1121/1.1908239.

4. Biot M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. J. Acoust. Soc. Am., 1956, vol. 28, pp. 179-191. Doi: 10.1121/1.1908241.

5. Johnson D.L. Recent developments in the acoustic properties of porous media. Frontiers in Physical Acoustics, XCIII, North Holland, 1986, pp. 255-290.

6. Knyaz'kov N.N., Sharfarets B.P. [Acoustics of porous-elastic fluid saturated medium (an overview of the Biot theory)]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2016, vol. 26, no. 1, pp. 77-84. URL: Doi: 10.18358/np-26-1-i7784. (In Russ.).

7. Polak L.S., ed. Variacionnye principy mekhaniki [Variation principles of mechanics]. Collection of articles. Moscow, Fizmatgiz Publ., 1959. 952 p. (In Russ.).

8. Polak L.S. Variacionnye principy mekhaniki. Ih razvitie i primenenie v fizike [Variation principles of mechanics. Their development and application in physics]. Moscow, URSS Publ., 2010. 600 p. (In Russ.).

9. Rayleigh, Lord. Teoriya zvuka [The Theory of Sound]. Vol. 1. Moscow, GITTL Publ., 1955. 504 p. (In Russ.).

10. Kurant R., Gilbert D. Metody matematicheskoj fiziki [Methods of mathematical physics]. Vol. 1. Moscow, GTTI Publ., 1933. 525 p. (In Russ.).

11. Mors F.M., Feshbah G. Metody teoreticheskoj fiziki [Methods of theoretical physics]. Vol. 1. Moscow, IIL Publ., 1958. 930 p. (In Russ.).

12. Bio M. Variacionnye principy v teorii teploobmena [The variation principles in the theory of heat exchange]. Moscow, Energy Publ., 1975. 209 p. (In Russ.).

13. Serrin J. Matematicheskie osnovy klassicheskoj mekhaniki zhidkosti [Mathematical principles of classical fluid mechanics]. Moscow, Inostrannaya literatura Publ., 1963. 256 p. (In Russ.).

14. Pierce A.D. Variational formulations in acoustic radiation and scattering. Physical Acoustics, vol. XXII, San Diego, Acad. Press, 1993. pp. 195-394.

15. Sedov L.I. [Mathematical methods of creation of new models of continuous environments]. Uspekhi matemati-cheskikh nauk [Achievements of mathematical sciences], 1965, vol. 20, no. 5, pp. 121-180. (In Russ.).

16. Prohorov A.M., ed. Fizicheskaya ehnciklopediya [Physical encyclopedia]. In 5 vol. Moscow, Sovetskaya ehnciklopediya Publ., 1988-1998. 704+704+672+704+760 p. (In Russ.).

17. Gibbs J. Termodinamicheskie raboty [Thermodynamic works]. Moscow, Leningrad, Gostekhizdat Publ., 1950. 492 p. (In Russ.).

18. Vajnberg M.M. Variacionnye metody issledovaniya neli-nejnyh operatorov [Variation methods of a research of nonlinear operators]. Moscow, Gostekhizdat Publ., 1956. 344 p. (In Russ.).

19. Vajnberg M.M. Variacionnyj metod i metod monotonnyh operatorov [Variation method and method of monotonous operators]. Moscow, Nauka Publ., 1972. 416 p. (In Russ.).

20. Faddeev L.D., ed. Ehnciklopediya matematicheskoj fiziki [Encyclopedia of mathematical physics]. Moscow, Bol'shaya rossijskaya ehnciklopediya Publ., 1998. 691 p. (In Russ.).

21. Landau L.D., Lifshic E.M. Teoreticheskaya fizika. T. 1.

Mekhanika [Theoretical physics. Vol. 1. Mechanics]. Moscow, Nauka Publ., 1988. 216 p. (In Russ.).

22. Landau L.D., Lifshic E.M. Teoreticheskaya fizika. T. VI. Gidrodinamika [Theoretical physics. Vol. VI. Hydrodynamics]. Moscow, Nauka Publ., 1988. 736 p. (In Russ.).

23. Sharfarets B.P., Sharfarets E.B., Knyaz'kov N.N., Pashov-kin T.N. [Some features of the numerical solution of tasks of thermoelasticity and hydrodynamics heat-conducting compressed viscous liquid by means of universal packages]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2016, vol. 26, no. 3, pp. 57-63. Doi: 10.18358/np-26-3-i5763. (In Russ.).

24. Yavorskiy B.M., Detlaf A.A., Lebedev A.K. Spravochnik po fizike. Dlya inzhenerov i studentov vuzov [Reference book on physics. For engineers and students of higher education institutions]. Мoscow, Oniks: Mir i obrazovanie Publ., 2006. 1056 p. (In Russ.).

25. Schlichting H. Teoriya pogranichnogo sloya [Boundary layer theory]. Moscow, Nauka Publ., 1974. 712 p. (In Russ.).

Contacts: Sharfarets Boris Pinkusovich, sharb@mail.ru

Article received in edition: 15.12.2016