электронный журнал открытого доступа Cardiometry - выпуск 14, май 2019
Очередной номер нашего журнала - не совсем обычный. Нами постоянно анализируется не только интерес читателей к журналу, но и то, как новая наука кардиометрия понимается рядовыми врачами и насколько они применяют её на практике. Бесспорно, математические основы гемодинамики очень тяжелы в понимании не только врачу, но даже подготовленному человеку. Очередной номер нашего журнала - не совсем обычный. Нами постоянно анализируется не только интерес читателей к журналу, но и то, как новая наука кардиометрия понимается рядовыми врачами и насколько они применяют её на практике. Бесспорно, математические основы гемодинамики очень тяжелы в понимании не только врачу, но даже подготовленному человеку.
шей, чем скорость роста интегральной энтропии,приведенная в формуле (49.6) в [4]. Действительно, в(49.6) присутствует член, пропорциональный коэффициентувторой вязкости, а в (2.9) такого членауже нет при условии (2.11). Этому уменьшению скоростироста энтропии в (2.9) соответствует при условии(2.11) и относительное уменьшение скоростидиссипации кинетической энергии в (2.5) по сравнениюс выражением (79.1) в [4]. При этом имеет местосходство с принципом минимума производстваэнтропии И. Пригожина (см. в [10]).Итак, для дивергентных течений сжимаемойсреды получено дополнительное уравнение (2.11),которое замыкает эту систему (1.1), (1.2), исходя изтребования положительной определенности скоростироста интегральной энтропии в (2.9) и отрицательнойопределенности скорости диссипацииинтегральной кинетической энергии в (2.5). Поэтомууравнение (2.11) для дивергентных теченийсжимаемой среды должно заменять условие бездивергентности,обычно используемое для замыканиясистемы (1.1), (1.2) в случае приближениянесжимаемой среды.3. Новое дивергентное решениеуравнения НС1. Условие (2.11) определяет точную взаимнуюкомпенсацию нормальных напряжений давленияи нормальных вязких напряжений сжимаемогодивергентного течения. В результате такой компенсацииобращается в нуль второй член в правойчасти уравнения (1.1). Уравнение (1.1) при этомточно совпадает с n- мерным обобщением уравненияБюргерса:∂u∂∂uuii+j= ∆it ∂xjρ(3.1)При этом система (1.2), (3.1) уже является замкнутойи описывает эволюцию плотности и поляскорости движения среды по инерции при наличиизатухания, связанного только с действием сдвиговыхвязких напряжений, соответствующих ненулевойправой части уравнения (3.1).Если же в (3.1) коэффициент вязкости равеннулю, то из (3.1) получается n- мерное уравнениеРХ, для которого в [13] получено точное вихревоерешение, рассмотренное далее и обобщенное дляслучая учета внешнего трения или эффективнойвязкости. Отметим, что в отличие от проведенногоηuздесь и в [13] рассмотрения вихревых решений,ранее исследовались лишь безвихревое решениеуравнения (3.1), соответствующее потенциальномутечению и получаемое при использовании модификациинелинейного преобразования Коула –Хопфа [27, 28].u → u + V()t ,Пусть в (3.1) осуществлена заменагде V i(t) случайное гауссово дельта-коррелированноево времени по ле скорости, для которогоV() tV( τ ) = 2νδ δ( t−τ)i j ijV()t = 0ii i i(3.2)имеют место соотношенияВ (3.2) δ ij– символ Кронекера, δ – дельта-функцияДирака – Хэвисайда, а коэффициент ν характеризуетдействие сил вязкости. Он в общемслучае может зависеть от времени, описывая эффективнуютурбулентную вязкость, но может исовпадать с постоянным коэффициентом кинематическойвязкости, когда рассматриваемое случайноеполе скорости соответствует молекулярнымфлуктуациям. Ограничимся рассмотрениемслучая, когда этот коэффициент в (3.2) являетсяпостоянным, но достаточно большим по величинетак, что имеет место неравенство, позволяющеепренебречь членом в правой части (3.1):ην >>(3.3)min( ρ)где min (ρ) – абсолютный минимум величиныплотности среды в пространстве и во времени.Указанная замена в (3.1) с введением случайногополя скорости, как уже отмечалось во Введении,соответствует применяемому в [15] методуполучения стохастического уравнения НС не засчет использования случайной силы, а путем добавленияслучайной скорости к полю скорости,входящему в обычное детерминированное уравнениеНС. Здесь, в отличие от [15], рассмотренслучай, когда такое случайное поле скорости зависиттолько от времени (в [15] такое поле скоростиназывается дрейфовой частью крупномасштабногонеоднородного случайного поля) и егоучет эквивалентен введению силы объемной вязкости,которая по структуре совпадает с обычнойсилой трения в уравнении НС.Из уравнения (3.1), усредненного с учетом (3.2)при условии (3.3), получаем уравнение28 | Выпуск 14, Май 2019
где угловые скобки соответствуют операции усредненияпо случайному гауссову полю скорости V i(t).При выводе уравнения (3.4) из (3.1), (3.2) кроменеравенства (3.3) используется следующее соотношение(являющееся следствием формулы Фуруцу – Новикова[29–31]):Уравнение (3.4) и без условия (3.3) может соответствоватьуравнению (3.1), если наряду с (3.5) выполненыравенстваu = u ,(см. [32]) и если в (3.1) предварительно провестизаменуη η→ min( ) = ν .ρ ρТакое расцепление корреляций возможно приточном разделении временных масштабов, связанныхс крупномасштабными инерционнымидвижениями и движениями с характерным масштабомвязкой диссипации [32, 23].2. Вместо того чтобы для нахождения среднегополя скорости‹ i› u приближенно решать (см., например,[32]) проблему замыкания при рассмотрениинепосредственно уравнения (3.4) воспользуемсяисходным уравнением, из которого точно следуетименно уравнение (3.4). Это исходное уравнениеимеет вид n- мерного уравнения РХ [10, 13, 22]:Если взять операцию ротора от левой частиуравнения (3.6), то в точности получится уравнениеГельмгольца (1.3), в котором только надо удалитьчлен ν Δ ω ir rи произвести замену u→u+ Vtr().Уравнение (3.6), как показано в [13], имеет следующееточное решение для случая произвольнойразмерности пространства (n = 1, 2, 3 и т. д.):rui(,) xt =r rnr r r rˆ (3.7)= d ξu ( ξδξ ) ( − x + B( t) + tu ( ξ)) det Aгде∫∂u∂uiVk∂x011rk=0k− −∂ui∂ui+ ( uj+ Vj( t))= 0 (3.6)∂t∂xтождеств [13]:jr r r r r∂δξ( − x + B( t) + tu0( ξ))∂xmr r r r−1 ∂δξ− x + B t + tu0≡−Akm∂ξkr r r r r r1tu0 1u0r rδξ (1−ξ)≡det Aˆ0i0t0dV τiτ ˆ∂unA≡ Anm= δnm+ t ,0∂ ξmBti() = ∫ ( )i∂tkji∂uiuj∂x∂ui+ uj= ν∆∂xj=−∆ νiju∂ui= uj∂xiui, (3.4)(3.5)det  – детерминант матрицы Â, а u ( ) 0 ix r – произвольноегладкое начальное поле скорости. Решение(3.7) удовлетворяет уравнению (3.6) только притаких временах, для которых при любых значенияхпространственных координат величина детерминантаматрицы  является положительной, т.е.det  > 0. Поэтому везде далее будем это учитывать и,соответственно, не будет использоваться знак модуляпри написании det Â, если не оговорено обратное.Решение (3.7) только в случае потенциальностиначального поля скорости является потенциальнымбезвихревым, соответствуя нулевому полювихря, для всех последующих моментов времени.Наоборот, оно является вихревым и определяетэволюцию поля вихря при ненулевом начальномполе вихря (см. следующий параграф). Далееограничимся рассмотрением только вихревыхрешений (3.7). Отметим, однако, что в [24] полученоименно потенциальное решение двумерногоуравнения РХ (3.6) (при B r = 0 в (3.6)) в лагранжевомпредставлении, которое точно следует ииз (3.7) при n = 2, как уже отмечено во Введении всвязи с возможностью описания решения уравненияСивашинского (В.2) с помощь потенциальногорешения (3.7).Для одномерного случая с n = 1 имеем в (3.7)det ˆ duA= 1+t dξи решение (3.7) точно совпадает с решениями, приведеннымив [33, 34]. Решение (3.7) получается,если использовать интегральное представлениедля неявного решения уравнения (3.6) в видеr r r ru( xt ,) u ( x Bt () tuxt ( ,))с использованием дельта-функции и следующихr( ( ) ( ξ)) ;δξ ( − ξ+ ( ( ξ) − ( ξ)))≡(3.8)(3.9)где А -1 km – матрица обратная матрице А km .После усреднения по случайному полю B i(t) (сгауссовой плотностью распределения вероятностей)из (3.7) можно получить точное решение≡Выпуск 14, Май 2019 | 29
- Page 1 and 2: Выпуск 14, Май 2019 | 95
- Page 3 and 4: Кардио-окулометрич
- Page 5 and 6: 56 Перспективы прим
- Page 7 and 8: Dr. Hong LeiЧунцинская Ш
- Page 9 and 10: Воронова Ольга Кон
- Page 11 and 12: Дорогой читатель!О
- Page 13 and 14: Рис. 1. Формирование
- Page 15 and 16: QT - длительность ин
- Page 17 and 18: 4. Руденко М.Ю., Зерн
- Page 19 and 20: Закономерности дви
- Page 21 and 22: Рис. 1. Физические м
- Page 23 and 24: ОРИГИНАЛЬНОЕ ИССЛЕ
- Page 25 and 26: Соболева Н 1 (R 3 ) для
- Page 27 and 28: Отметим, что обычно
- Page 29: формулах (49.3) и (49.4),
- Page 33 and 34: времени t. Например,
- Page 35 and 36: использовать явное
- Page 37 and 38: rn 1d ξρ0( ξ ) ⋅(2 )ntπνρ =
- Page 39 and 40: стического функцио
- Page 41 and 42: ОТЧЕТ Подача: 15.1.2019;
- Page 43 and 44: Рис. 2. Изменения ам
- Page 45 and 46: г)Время Лактат КрФ
- Page 47 and 48: Визит Владимира Зе
- Page 49 and 50: Руководству, в том
- Page 51 and 52: методической погре
- Page 53 and 54: ЛЕКЦИИЗаконы и акс
- Page 55 and 56: Использование зако
- Page 57 and 58: ОРИГИНАЛЬНОЕ ИССЛЕ
- Page 59 and 60: ОРИГИНАЛЬНОЕ ИССЛЕ
- Page 61 and 62: эффициент корреляц
- Page 63 and 64: а)б)Рис. 5. Диаграмма
- Page 65 and 66: рисунке 9 представл
- Page 67 and 68: ОРИГИНАЛЬНОЕ ИССЛЕ
- Page 69 and 70: мулов, одна из кото
- Page 71 and 72: Таблица 3. Основные
- Page 73 and 74: ния [14]. Наименьшее
- Page 75 and 76: 14. Rudenko M., Voronova O., Zernov
- Page 77 and 78: ВведениеИспользов
- Page 79 and 80: стимул 7 - «вы пробо
шей, чем скорость роста интегральной энтропии,
приведенная в формуле (49.6) в [4]. Действительно, в
(49.6) присутствует член, пропорциональный коэффициенту
второй вязкости, а в (2.9) такого члена
уже нет при условии (2.11). Этому уменьшению скорости
роста энтропии в (2.9) соответствует при условии
(2.11) и относительное уменьшение скорости
диссипации кинетической энергии в (2.5) по сравнению
с выражением (79.1) в [4]. При этом имеет место
сходство с принципом минимума производства
энтропии И. Пригожина (см. в [10]).
Итак, для дивергентных течений сжимаемой
среды получено дополнительное уравнение (2.11),
которое замыкает эту систему (1.1), (1.2), исходя из
требования положительной определенности скорости
роста интегральной энтропии в (2.9) и отрицательной
определенности скорости диссипации
интегральной кинетической энергии в (2.5). Поэтому
уравнение (2.11) для дивергентных течений
сжимаемой среды должно заменять условие бездивергентности,
обычно используемое для замыкания
системы (1.1), (1.2) в случае приближения
несжимаемой среды.
3. Новое дивергентное решение
уравнения НС
1. Условие (2.11) определяет точную взаимную
компенсацию нормальных напряжений давления
и нормальных вязких напряжений сжимаемого
дивергентного течения. В результате такой компенсации
обращается в нуль второй член в правой
части уравнения (1.1). Уравнение (1.1) при этом
точно совпадает с n- мерным обобщением уравнения
Бюргерса:
∂u
∂
∂u
u
i
i
+
j
= ∆
i
t ∂xj
ρ
(3.1)
При этом система (1.2), (3.1) уже является замкнутой
и описывает эволюцию плотности и поля
скорости движения среды по инерции при наличии
затухания, связанного только с действием сдвиговых
вязких напряжений, соответствующих ненулевой
правой части уравнения (3.1).
Если же в (3.1) коэффициент вязкости равен
нулю, то из (3.1) получается n- мерное уравнение
РХ, для которого в [13] получено точное вихревое
решение, рассмотренное далее и обобщенное для
случая учета внешнего трения или эффективной
вязкости. Отметим, что в отличие от проведенного
η
u
здесь и в [13] рассмотрения вихревых решений,
ранее исследовались лишь безвихревое решение
уравнения (3.1), соответствующее потенциальному
течению и получаемое при использовании модификации
нелинейного преобразования Коула –
Хопфа [27, 28].
u → u + V()
t ,
Пусть в (3.1) осуществлена замена
где V i
(t) случайное гауссово дельта-коррелированное
во времени по ле скорости, для которого
V() tV( τ ) = 2νδ δ( t−τ
)
i j ij
V()
t = 0
i
i i i
(3.2)
имеют место соотношения
В (3.2) δ ij
– символ Кронекера, δ – дельта-функция
Дирака – Хэвисайда, а коэффициент ν характеризует
действие сил вязкости. Он в общем
случае может зависеть от времени, описывая эффективную
турбулентную вязкость, но может и
совпадать с постоянным коэффициентом кинематической
вязкости, когда рассматриваемое случайное
поле скорости соответствует молекулярным
флуктуациям. Ограничимся рассмотрением
случая, когда этот коэффициент в (3.2) является
постоянным, но достаточно большим по величине
так, что имеет место неравенство, позволяющее
пренебречь членом в правой части (3.1):
η
ν >>
(3.3)
min( ρ)
где min (ρ) – абсолютный минимум величины
плотности среды в пространстве и во времени.
Указанная замена в (3.1) с введением случайного
поля скорости, как уже отмечалось во Введении,
соответствует применяемому в [15] методу
получения стохастического уравнения НС не за
счет использования случайной силы, а путем добавления
случайной скорости к полю скорости,
входящему в обычное детерминированное уравнение
НС. Здесь, в отличие от [15], рассмотрен
случай, когда такое случайное поле скорости зависит
только от времени (в [15] такое поле скорости
называется дрейфовой частью крупномасштабного
неоднородного случайного поля) и его
учет эквивалентен введению силы объемной вязкости,
которая по структуре совпадает с обычной
силой трения в уравнении НС.
Из уравнения (3.1), усредненного с учетом (3.2)
при условии (3.3), получаем уравнение
28 | Выпуск 14, Май 2019