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

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

07.01.2020 Views

шей, чем скорость роста интегральной энтропии,приведенная в формуле (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

шей, чем скорость роста интегральной энтропии,

приведенная в формуле (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

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

Saved successfully!

Ooh no, something went wrong!