07.01.2020 Views

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

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

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

уравнения (3.4) (и уравнения (1.1) при условии

(2.11)) в следующем виде

r

n ˆ 1

ui = ∫ d ξu0i( ξ) det A ⋅

n

(2 πν t )

r r r r (3.10)

2

⎡ ( x −ξ

−tu0( ξ))

⋅exp

⎢−

⎥.

⎣ 4ν

t ⎦

Усредненное решение (3.10), в отличие от (3.7),

является уже сколь угодно гладким на любом неограниченном

отрезке изменения времени, а не только

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

матрицы Â.

r

3. Без учета сил вязкости, когда в (3.7) Bt () = 0,

гладкое решение (3.7), как уже отмечалось, определено

лишь при условии det  > 0 [13]. Оно соответствует

ограниченному интервалу времени 0 ≤ t < t 0

,

где величина предельного минимального времени

существования решения определяется из решения

следующего алгебраического уравнения порядка

n (и последующей минимизации полученного

выражения, зависящего от пространственных

координат, по этим координатам):

ˆ du ( x )

det At ( ) = + t = , n=

dx

где det Û 0

– детерминант трехмерной матрицы

U 0nm

= ∂u 0n

∕ ∂x m

, а

ˆ ∂u01 ∂u02 ∂u01 ∂u02

detU

012

= −

∂x1 ∂x2 ∂x2 ∂x1

детерминант аналогичной матрицы в двумерном

случае для переменных (x 1

, x 2

). При этом det Û 013

, det

Û 023

– детерминанты матриц в двумерном случае

для переменных (x 1

, x 3

) и (x 2

, x 3

), соответственно.

Отметим, что в двумерном случае (3.11) точно

совпадает с условием коллапса, полученным в [24] в

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

исследованной на основании уравнения Сивашинского

(В.2). Для точного совпадения надо при этом

в (3.11) заменить

t→ bt =

01 1

1 0 1

ˆ r

det A( t)

= 1+ tdivu0

+

+ det ˆ = 0,

= 2

2

t U012

n

ˆ r

det A( t)

= 1+ tdivu0

+

2

+ t (detUˆ ˆ ˆ

012

+ detU013 + det U023)

+

+ ˆ = 0 = 3

3

t det U0

, n

Us 0

1

() (exp( γ t ) − )

1

γ

0

.

(3.11)

В одномерном случае при n = 1 из (3.11) имеем минимальное

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

1

t0

= > 0.

du01( x1)

max

dx

В частности, при начальном распределении

2

x1

u01( x1) = aexp( − ), a > 0

2

L

получаем

для значения

x

L

t0

=

a

1 1max

При этом сама реализация сингулярности может

иметь место только при положительных значениях

координаты x 1

> 0, когда уравнение (3.11) имеет

положительное решение для величины времени.

Это означает, что сингулярность (коллапс)

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

случае, когда начальное поле скорости отлично от

нуля только при отрицательных значениях пространственной

координаты x 1

< 0.

Аналогичным образом определяется минимальная

величина времени опрокидывания волны t 0

и

при n > 1. Так для (3.11) в двумерном случае (при

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

для начальной функции тока в виде

0 1 2

x x

= a LL exp( − − ), a>

0

2 2

1 2

1 2 2 2

L1 L2

1

ψ ( x , x ) =

минимальная величина времени существования

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

e LL

1 2

t .

0

=

2a

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

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

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

соответствующих точкам эллипса

x

L

= x =

x

+ = 1

2 2

1 2

2 2

1

L2

Согласно (3.11) необходимым условием реализации

сингулярности является условие существования

действительного положительного решения

квадратного (при n = 2) или кубического

уравнения (при n = 3) относительно переменной

e

2

L

2

.

.

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

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

Saved successfully!

Ooh no, something went wrong!