2. Первоначальное упрощение исходной модели.

Кинетические уравнения в исходной системе можно сильно упростить, проанализировав значения входящих в них констант и диапазон изменения переменных. В табл. 3 представлены максимальные концентрации активных факторов, достигаемые при моделировании роста тромба по уточненной модели, а также минимальные концентрации их предшественников. Прежде всего видно, что зимогены факторы XI, IX, X и pC действительно практически не исчерпываются, и их концентрации могут быть положены константами, а концентрации их активных форм – факторов XIa, IXa, Xa и pCa – вырастают лишь до чрезвычайно низких абсолютных значений, которыми можно пренебречь в знаменателях выражений для Z и W (см. Приложения А-В).

Табл. 3. Изменения содержания факторов свертывания при пространственном росте тромба (расчет по уточненной модели). Все концентрации выражены в нМ.

Фактор	Исходное содержание	Минимум предшественника	Максимум активного фактора
VII, VIIa	10.1	4.172	5.924
XI,XIa	30	29.99	0.0007066
IX,IXa	90	89.95	0.02225
X, Xa	170	169.1	0.4853
II,IIa	1400	189.9	132.5
VIII, VIIIa	0.7	0	0.3513
V, Va	20	0	19.7
pC, pCa	60	59.51	0.2846

Изменение вида всех уравнений системы на последовательных стадиях упрощения приведены в приложениях Б и В. В Приложении Г приведена система уравнений после перехода к безразмерным концентрациям, изменяющимся в диапазоне от 0 до 1.

3. Вклад диффузионных членов в динамику поведения системы.

Для оценки количественного вклада диффузионных членов была проведена серия расчетов с пренебрежением различными диффузионными членами кинетических уравнений. Расчет проводился в двух режимах. В первом просто пренебрегали одним или несколькими диффузионными членами системы. Во втором режиме расчет проводился так, чтобы пренебречь диффузионными членами в области роста сгустка, сохранив их в приповерхностной области, с тем, чтобы оценить вклад диффузии факторов именно при роста тромба. Для этого диффузионный член домножался на S-образную функцию , равную 1 около активирующей поверхности и резко (на интервале порядка шага по пространству, т.е. 15-30 нм) спадающей до 0 на некотором от нее расстоянии, в районе =₀. Значение ₀=L/20 (L – длина отрезка интегрирования) выбрано исходя из динамики профилей: на этом расстоянии влияние стенки уже довольно мало.

Результаты этих расчетов собраны в табл. 4. Из данных таблицы можно сделать следующие выводы:

1. роль диффузии предшественников наиболее слаба;

2. вклад диффузионного члена Xa фактора заметен только в фазе инициации – до отрыва фронта активации плазмы от активирующей поверхности;

3. в фазе роста важны диффузионые члены IIa, IXa, а при гемофилии А – еще и Va факторов.

Дальнейшую работу проводили при пренебрежении диффузией всех факторов, кроме IIa, IXa и Xa.

Таблица 4. Пренебрежение различными диффузионными членами. Расчет по полной уточненной модели.

Недиффундирующий фактор	Область пространства, где пренебрегается диффузией	Скорость роста на 40-й минуте, нм/мин	Размер сгустка на 40-й минуте, нм	Интегральный тромбиновый потенциал, нМ*мин
Норма
все диффундируют		44.9	1.59	361
Все предшественники	0..L	46.8	1.64	354
	L/20..L	46.9	1.65	379
IIa	0..L	37.5	1.33	336
	L/20..L	37.4	1.3	321
IXa	0..L	30	1.04	170
	L/20..L	30.5	1.15	208
Xa	0..L	0	0	0
	L/20..L	43.6	1.54	346
Va	0..L	43.2	1.53	352
	L/20..L		1.52	349
VIIIa	0..L	43.9	1.55	378
	L/20..L		1.54	372
VIIa	0..L	44.4	1.59	359
	L/20..L	44.3	1.58	356
XIa	0..L	43.3	1.55	348
	L/20..L			347
pCa	0..L	44.9	1.59	368
	L/20..L			370
Диффундируют только IIa, IXa	0..L	0	0	0
	L/20..L	42.4	1.46	364
Диффундируют только IIa, IXa и Xa	0..L	43.5	1.53	365
	L/20..L		1.51	377
1% гемофилия А
все диффундируют		16.0	0.60	3.1
Диффундируют только IIa, IXa	0..L	0	0	0
	L/20..L	5.2	0.3	1.4
Диффундируют только IIa, IXa и Xa	0..L	5.9	0.33	3.3
	L/20..L	5.5	0.31	1.9
Диффундируют только IIa, IXa, Xa, Va	0..L	15.9	0.54	3.5
	L/20..L	16.2	0.53	2.8