[hr]

[hr]

Математический аппарат моделирования

Система ОДУ представляет собой сложную конструкцию взаимосвязанных уравнений, описывающих динамику изменения массы различных компонентов системы во времени. Ключевой особенностью данной системы является её жёсткость, проявляющаяся в существенном различии временных масштабов протекающих процессов.

Особенности моделирования:

• Быстрые процессы, такие как всасывание питательных веществ, происходят в масштабе секунд
• Медленные процессы, например, разложение лигнина, растягиваются на десятилетия

• Автоматическая адаптация шага интегрирования
• Строгое соблюдение закона сохранения массы
• Высокая устойчивость при решении жёстких задач

Анатомия микоризного симбиоза

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

Метаболическое ядро (Капсула):

• Представляет собой стабильную структуру, неразрывно связанную с корневой системой растения
• Выполняет функцию основного интерфейса обмена веществом и энергией
• Обеспечивает двусторонний транспорт углерода и минеральных элементов

• Динамичная система гиф, организованная по трёхзонному принципу Торнли
• Характеризуется постоянным цикличным обновлением структуры
• Осуществляет активное исследование пространства

// Комментарий: описание уникального результата моделирования

Численная симуляция 19‑мерной тензорной системы ОДУ, реализованная на базе NumPy с неявным интегратором Radau, позволила зафиксировать удивительный результат.

Геометрическая концепция ризосферного хотспота Торнли, где эффективный объем взаимодействия определяется как $V_{eff} = MC_{Seff} \cdot \Delta x$ при толщине прикорневой зоны $\Delta x = 10$ см, привела к следующим результатам:

• Масса микоризной сети (AM) стабилизировалась на уровне 25141.4 г
• Биомасса сапротрофов установилась на отметке 4005.4 г
• Показатель pH почвы зафиксировался на исходном значении 6.000

Весенняя экспансия (зоны 1–2):

• Поисковые гифы оккупируют богатые участки
• Происходит пассивное всасывание по градиенту Фика
• Протонный поток $H^+$ отсутствует

• Исчерпывание легких ресурсов
• Активация триггера «ИЛИ»
• Включение протонных помп для извлечения элементов

• Зануление активности гиф температурным Гауссом
• Пассивная буферизация почвы
• Восстановление нейтрального pH

$H_{flux\_AM}=k_{proton\_prod}\cdot y_{[\_W\_AM\_TIPS]}\cdot element_{input}[AM,N]=0$

[hr]

Фрагмент кода векторизованного ОДУ‑ядра:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

import numpy as np
import params as p
 
def _rhs(t, y):
    dy = np.zeros(19)
    W_AM_total = y + y[1](https://www.cyberforum.ru/blogs/713657/) + y[2](https://cepl.rssi.ru/wp-content/uploads/2019/10/Sokolova_forestsoil2019.pdf) # core + tips + old
    
    # Расчет физического времени сезона
    hour_of_year = t % 8760.0
    if hour_of_year < 1680.0:  # 5 двухнеделек штурма
        zone_index = int(hour_of_year // 336.0)
        k_proton_dynamic = p.k_proton_prod * (zone_index / 4.0)
    else:
        k_proton_dynamic = 0.0  # Зимний анабиоз
        
    # Геометрический объем хотспота
    MC_Seff = p.mask_contact * (y[1](https://www.cyberforum.ru/blogs/713657/) / p.RoSq_M)
    rhizo_hotspot_volume = MC_Seff * 0.1 + 1e-6
    
    # Градиентный шлюз Торнли
    S_compost = y / rhizo_hotspot_volume
    S_mycelium = y[1](https://www.cyberforum.ru/blogs/713657/) / (W_AM_total + 1e-10)
    gradient = p.U_ratio * S_compost - S_mycelium
    
    # ТЕНЗОРНЫЙ СТЕХИОМЕТРИЧЕСКИЙ ШЛЮЗ
    substrat_input = np.zeros((8, 10))
    substrat_input[:, p.GLUCOSE] = p.mask_contact * (p.Z_diff * MC_Seff * np.maximum(0.0, gradient))
    
    # Стехиометрия
    element_input = substrat_input @ p.SUBSTANCE_STOICHIOMETRY
    
    # Водородный контур
    H_flux_AM = k_proton_dynamic * y[1](https://www.cyberforum.ru/blogs/713657/) * element_input[p.AM, p.N]
    dy = H_flux_AM - p.k_proton_loss * (p.pH_0 - p.pH_current)
    
    return dy

.Диффузионно-стехиометрическая ловушка "сверхтекучих" ферментов.

Важнейшим физико-химическим ограничением модели является толщина слоя взаимодействия (\(\Delta x = 10\) см). Казалось бы, высокая проникающая способность ферментов гриба должна повышать вынос элементов. Однако математический анализ градиента Торнли демонстрирует обратное:При увеличении \(\Delta x\) экспоненциально растет эффективный объем хотспота (\(V_{eff}\)). Это приводит к мгновенному падению локальной концентрации высвобожденного субстрата (\(S_{compost}\)) из-за эффекта разбавления. Градиент пассивной диффузии по закону Фика зануляется или инвертируется, блокируя всасывание мономеров обратно в мицелий. Таким образом, микоризный синдикат эволюционно оптимизирован на "короткорадиусный" вязкий гидролиз: выжигание почвы эффективно только тогда, когда оно локализовано в ультратонком приграничном слое фазового контакта.

Эффективный радиус фермента как триггер межвидовой войны.

Таким образом, эффективный радиус действия ферментов (\(\Delta x\)) является главным регулятором подземного баланса сил. Эволюционное сжатие этого радиуса для ухода от стехиометрической ловушки неизбежно разворачивает агрессию мицелия с дальней бесплодной почвы на ближний прикорневой контур. Короткорадиусный ферментативный удар перенаправляется на корневой опад и ткани самого растения-хозяина. Этот пространственный маневр превращает микоризу из симбионта в жесткого биохимического агрессора. Гриб замыкает азотно-фосфорный цикл внутри своего меланизированного сейфа, провоцируя глубокий дефицит питания у дерева, пробивая вековой буфер ризосферы и запуская необратимый маховик лесной сукцессии.

.