электронный журнал открытого доступа Cardiometry - выпуск 14, май 2019

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

07.01.2020 Views

бенности (в том числе сингулярность энстрофии)эволюции во времени для завихренности. Полученноев настоящей работе описание эволюциизавихренности в лагранжевом представлениидля двумерного и трехмерного случаев (см. (4.4) и(4.5)) при этом можно рассматривать как конкретизациюполученной в [23] формы решения уравненияЭГ для случая движения жидких частиц поинерции.Кроме того, в настоящей работе получен новыйнеобходимый и достаточный критерий реализациивзрывной сингулярности (коллапса) законечное время (см.(3.11), (3.12)) для невязкогорешения уравнений РХ и ЭГ в одномерном, двумерноми трехмерном случаях. В то же время, в[22] приведен интегральный критерий в виде (3.13)(см. формулу (38) в [22]), который определяет лишьдостаточное условие для реализации коллапсарешения. При этом, например, для случая начальногобездивергентного поля скорости коллапсвозможен только согласно необходимому и достаточномукритерию (3.12), но уже не может бытьустановлен из критерия (3.13). Кроме того, из полученногов [22] рассмотрения взрывного режимадля решения трехмерного уравнения РХ делаетсявывод о невозможности продолжения этого решенияна бесконечное время в пространстве СоболеваН 2 (R 3 ), что отличается от полученного в настоящейработе результата, отмеченного выше.В двумерном же случае имеется точное соответствиекритерия (3.11) с аналогичным критерием,приведенным в [24] (см. формулу (9) в [24])в связи с решением задачи о распространениифронта пламени (генерируемого самоподдерживающейсяэкзотермической химической реакцией)на основе упрощенной версии уравнения Сивашинского[25]:∂f1− U∂t22( ∇f) = γ fs 0В уравнении (В.2) функцияxr= f( x , x ,) t3 1 2(В.2)определяет фронт пламени, представляющийсобой границу раздела между горючим веществом(x 3>0 ) и продуктами сгорания (x 3< 0 ), а U sи γ 0постоянныеположительные величины, характеризующиескорость распространения фронта и интенсивностьгорения, соответственно. При γ 0= 0уравнение (В.2) совпадает с уравнением Гамильтона– Якоби для свободной нерелятивистскойчастицы. Полученное в настоящей работе точноерешение (3.7) уравнения РХ в двумерном случае(точнее в его модификации при учете внешнего тренияс коэффициентом μ и формальном равенствеμ = – γ 0) дает и точное решение уравнения (В.2). Приэтом решение (3.7) описывает потенциальное течениевидаru =−U ∇ rf .5. Важным результатом настоящей работы являетсяполучение замкнутого описания эволюции вовремени энстрофии и любых более высоких моментовполя вихря, а также поля скорости в двумерноми трехмерном случаях. Это достигнуто исходя изсоответствующих аналитических решений уравненийЭГ, РХ и НС для случая нулевой вязкости и дляслучая учета внешнего трения или эффективнойвязкости. В результате точно, а не как обычно приближенно,решена проблема замыкания в теориитурбулентности, которая оставалась нерешенной,несмотря на многочисленные попытки ее приближенногорешения [1]. Здесь же это удается сделатьтолько благодаря относительно простой явнойзависимости от начальных условий для полученныхточных решения уравнений ЭГ и РХ для поляскорости (3.7) и поля вихря ((4.1) и (4.2)), чего нет,например, в известном точном решении уравненияБюргерса, получаемом при использовании нелинейногопреобразования Коу ла – Хопфа.В частности, благодаря этому на основе точныхрешений (4.2) могут быть получены следующиеоценки для интегралов от поля завихренности втрехмерном случае вблизи момента сингулярностирешения:1Ω d r ω ) и∫3 2m3(2m)= x ≅ O( 2m−1( t0− t)3 m 1Ω( m)= ∫ d xr ω ≅ Om( t − t)3(−10) ,когда m = 1, 2, 3... . Отсюда непосредст венно следуетнеравенство,2⎛ 1 ⎞Ω3(2 m) / Ω3( m)≅ O⎜⎟>>1;⎝ t0− t ⎠t → t(В.3)которое свидетельствует о сильной перемежаемостиполя вихря вблизи сингулярности.s024 | Выпуск 14, Май 2019

Отметим, что обычно неравенствоΩдействительно, считается выполненным при сильнойвихревой перемежаемости [15], но ранее егоне было возможности получить из точного решенияпроблемы замыкания в теории турбулентности,как это сделано при получении оценки (В.3).6. В заключительной части настоящей работына основе анализа точного замкнутого решенияуравнения баланса энстрофии (5.6) и скоростиизменения интегральной кинетической энергии в(6.1) – (6.4) обсуждается возможность существованияна неограниченном интервале времени нетолько дивергентных, но и гладких бездивергентныхрешений уравнения НС.1. Уравнение Навье – Стокса (НС) иуравнение Эйлера – Гельмгольца (ЭГ)Уравнение движения сжимаемой среды в общемслучае могут быть записаны в следующем виде [4]:∂u∂tii+ uj=∂xjη 1 ∂ η ∂uk= ∆ui− ( p− ( ζ + )( ));ρ ρ ∂x3 ∂xi>∂u∂Из вида второго члена в правой части (1.1) следует,что для вязкого сжимаемого дивергентноготечения нормальные напряжения определяются нетолько давлением, но и дивергенцией поля скорости.В (1.1), (1.2) u i– скорость среды, по повторяющимсяиндексам здесь и далее подразумеваетсясуммирование от 1 до n (где n – размерность пространстваи далее будут рассмотрены случаи, когдаn = 1, 2, 3), а р, ρ, η, ζ – давление, плотность, постоянныекоэффициенты вязкости и второй вязкостисреды, соответственно [4].Для вязкой несжимаемой среды, имеющей постояннуюплотность ρ = ρ 0из уравнения (1.1) втрехмерном случае (при n = 3) после взятия операцииротора от левой и правой части получаетсяследующее уравнение Эйлера – Гельмгольца (ЭГ)r rВ (1.3) ω = rotu , а ν = η∕ρ0 = const – коэффициентмолекулярной кинетической вязкости.2Ω23(2m)3( m)∆= ∂ xk∂ xk∂ρ∂+ ( ρuk) = 0∂t∂xk,k(1.1)(1.2)∂ωi∂ωi+ uk=∂t∂xk∂u= ω − ωdivu+ ν∆ωik i i∂xk(1.3)Для случая сжимаемой среды уравнение (1.3)также имеет место, но лишь при условиях, чтоη∕ρ = const и равенства нулю ротора от второго членав правой части (1.1). В частности, это имеет местов рассмотренном далее случае равенства нулювторого члена в правой части (1.1), что соответствуетнулевому суммарному балансу нормальных напряжений,обусловленных давлением и вязкостьюдивергентного течения среды.В [13] получено точное вихревое решение трехмерногоуравнения Римана – Хопфа (РХ) (совпадающегос (1.1) при обращении в нуль правой части(1.1)) при произвольных, гладких, обращающихсяв нуль на бесконечности, начальных условиях. Оносовпадает и с точным решением уравнения ЭГ (1.3)для сжимаемой невязкой среды (когда в (1.3) равеннулю последний член в правой части). При этом в[13], в частности, показано, что полученное гладкоерешение может существовать только на ограниченноминтервале времени 0 ≤ t < t 0, (где величина t 0определяется далее из решения уравнения (3.11)).Далее (в третьем параграфе) показано, что длялюбой сколь угодно малой величины эффективнойвязкости (введенной вместо последнего члена вправой части (1.3)) можно получить точное решениеуравнений (1.1) – (1.3), существующее уже нанеограниченном интервале времени.2. Уравнения баланса энергии и энтропии1.Обычно при рассмотрении системы из четырехуравнений (1.1), (1.2) для пяти неизвестных функцийвводят дополнительно условие связи (уравнениесостояния среды) между плотностью и давлениемдля того что бы число уравнений совпало счислом неизвестных функций. Представление обуравнении состояния для неравновесных вихревыхтечений нуждается в уточнении. Вместо этого,для замыкания системы (1.1), (1.2) обычно используютприближение нулевой дивергенции поля скоростидля несжимаемой среды, что, в частности,оправдано при относительно малых (по сравнениюсо скоростью звука) скоростях движения среды.Выведем аналогичное уравнение, замыкающеесистему (1.1), (1.2) для дивергентных теченийсжимаемой среды, которое заменит условиеВыпуск 14, Май 2019 | 25

Отметим, что обычно неравенство

Ω

действительно, считается выполненным при сильной

вихревой перемежаемости [15], но ранее его

не было возможности получить из точного решения

проблемы замыкания в теории турбулентности,

как это сделано при получении оценки (В.3).

6. В заключительной части настоящей работы

на основе анализа точного замкнутого решения

уравнения баланса энстрофии (5.6) и скорости

изменения интегральной кинетической энергии в

(6.1) – (6.4) обсуждается возможность существования

на неограниченном интервале времени не

только дивергентных, но и гладких бездивергентных

решений уравнения НС.

1. Уравнение Навье – Стокса (НС) и

уравнение Эйлера – Гельмгольца (ЭГ)

Уравнение движения сжимаемой среды в общем

случае могут быть записаны в следующем виде [4]:

∂u

∂t

i

i

+ u

j

=

∂xj

η 1 ∂ η ∂uk

= ∆ui

− ( p− ( ζ + )( ));

ρ ρ ∂x

3 ∂x

i

>

∂u

Из вида второго члена в правой части (1.1) следует,

что для вязкого сжимаемого дивергентного

течения нормальные напряжения определяются не

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

В (1.1), (1.2) u i

– скорость среды, по повторяющимся

индексам здесь и далее подразумевается

суммирование от 1 до n (где n – размерность пространства

и далее будут рассмотрены случаи, когда

n = 1, 2, 3), а р, ρ, η, ζ – давление, плотность, постоянные

коэффициенты вязкости и второй вязкости

среды, соответственно [4].

Для вязкой несжимаемой среды, имеющей постоянную

плотность ρ = ρ 0

из уравнения (1.1) в

трехмерном случае (при n = 3) после взятия операции

ротора от левой и правой части получается

следующее уравнение Эйлера – Гельмгольца (ЭГ)

r r

В (1.3) ω = rotu , а ν = η∕ρ0 = const – коэффициент

молекулярной кинетической вязкости.

2

Ω

2

3(2m)

3( m)

∆= ∂ xk

∂ xk

∂ρ

+ ( ρuk

) = 0

∂t

∂xk

,

k

(1.1)

(1.2)

∂ωi

∂ωi

+ uk

=

∂t

∂xk

∂u

= ω − ωdivu

+ ν∆ω

i

k i i

∂xk

(1.3)

Для случая сжимаемой среды уравнение (1.3)

также имеет место, но лишь при условиях, что

η∕ρ = const и равенства нулю ротора от второго члена

в правой части (1.1). В частности, это имеет место

в рассмотренном далее случае равенства нулю

второго члена в правой части (1.1), что соответствует

нулевому суммарному балансу нормальных напряжений,

обусловленных давлением и вязкостью

дивергентного течения среды.

В [13] получено точное вихревое решение трехмерного

уравнения Римана – Хопфа (РХ) (совпадающего

с (1.1) при обращении в нуль правой части

(1.1)) при произвольных, гладких, обращающихся

в нуль на бесконечности, начальных условиях. Оно

совпадает и с точным решением уравнения ЭГ (1.3)

для сжимаемой невязкой среды (когда в (1.3) равен

нулю последний член в правой части). При этом в

[13], в частности, показано, что полученное гладкое

решение может существовать только на ограниченном

интервале времени 0 ≤ t < t 0

, (где величина t 0

определяется далее из решения уравнения (3.11)).

Далее (в третьем параграфе) показано, что для

любой сколь угодно малой величины эффективной

вязкости (введенной вместо последнего члена в

правой части (1.3)) можно получить точное решение

уравнений (1.1) – (1.3), существующее уже на

неограниченном интервале времени.

2. Уравнения баланса энергии и энтропии

1.Обычно при рассмотрении системы из четырех

уравнений (1.1), (1.2) для пяти неизвестных функций

вводят дополнительно условие связи (уравнение

состояния среды) между плотностью и давлением

для того что бы число уравнений совпало с

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

уравнении состояния для неравновесных вихревых

течений нуждается в уточнении. Вместо этого,

для замыкания системы (1.1), (1.2) обычно используют

приближение нулевой дивергенции поля скорости

для несжимаемой среды, что, в частности,

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

со скоростью звука) скоростях движения среды.

Выведем аналогичное уравнение, замыкающее

систему (1.1), (1.2) для дивергентных течений

сжимаемой среды, которое заменит условие

Выпуск 14, Май 2019 | 25

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!