Расщепленная схема Кранка-Николсона с прыгающей лягушкой для задачи нестационарного биоконвективного потока с концентрационно-зависимой вязкостью

Основная информация

• ID статьи: 2510.14034
• Название: A decoupled Crank-Nicolson leap-frog scheme for the unsteady bioconvection flows problem with concentration dependent viscosity
• Автор: Chenyang Li (Факультет математических наук, Восточнокитайский педагогический университет)
• Классификация: math.NA cs.NA
• Дата публикации: 15 октября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2510.14034

Аннотация

В данной работе предложена и проанализирована полностью дискретная схема Кранка-Николсона с прыгающей лягушкой (CNLF) для задачи нестационарного биоконвективного потока с концентрационно-зависимой вязкостью. Пространственная дискретизация использует метод конечных элементов Галеркина, временная дискретизация применяет метод CNLF для линейных членов и полуявный метод для нелинейных членов. Доказано, что схема безусловно устойчива, то есть временной шаг не ограничен сверху. Методом энергии получены оптимальные оценки ошибок в норме $L^2$ для скорости и концентрации. Численные эксперименты подтверждают теоретические результаты.

Исследовательский контекст и мотивация

Предметная область

Биоконвекция — это явление конвекции жидкости, вызванное движением микроорганизмов, имеющее важное значение в биологии, экологии и инженерных приложениях. Это явление описывается связанными уравнениями типа Навье-Стокса и уравнением конвекции-диффузии:

1. Уравнения гидродинамики: описывают течение несжимаемой вязкой культуральной жидкости
2. Уравнение переноса микроорганизмов: описывает процесс переноса микроорганизмов

Основные вызовы

1. Концентрационно-зависимая вязкость: в отличие от классических ньютоновских жидкостей, вязкость реальных суспензий зависит от концентрации микроорганизмов
2. Нелинейная связанность: существует сильная связь между полем скорости и полем концентрации
3. Численная устойчивость: необходимо разработать численную схему, которая одновременно устойчива и эффективна

Исследовательская мотивация

Существующие методы имеют следующие ограничения при работе с задачами биоконвекции с концентрационно-зависимой вязкостью:

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

Основные вклады

1. Предложена полностью дискретная схема CNLF: впервые применен метод Кранка-Николсона с прыгающей лягушкой к задаче биоконвекции с концентрационно-зависимой вязкостью
2. Доказана безусловная устойчивость: временной шаг не ограничен сверху, что повышает практическую применимость численной схемы
3. Установлены оптимальные оценки ошибок: получены оптимальные порядки сходимости скорости и концентрации в норме $L^2$
4. Предоставлен расщепленный алгоритм: полуявная обработка позволяет решать только линейные системы на каждом временном шаге, повышая вычислительную эффективность

Описание метода

Постановка задачи

Рассмотрим модель биоконвекции в ограниченной области $\Omega \subset \mathbb{R}^d$ ($d=2$ или $3$):

$\frac{\partial u}{\partial t} - \text{div}(\nu(c)D(u)) + u \cdot \nabla u + \nabla p = -g(1+\gamma c)i_2 + f$

$\nabla \cdot u = 0$

$\frac{\partial c}{\partial t} - \theta\Delta c + u \cdot \nabla c + U\frac{\partial c}{\partial x_2} = 0$

где:

• $u$ — поле скорости, $p$ — давление, $c$ — поле концентрации
• $\nu(c)$ — функция вязкости, зависящая от концентрации
• $D(u) = \frac{1}{2}(\nabla u + \nabla u^T)$ — тензор деформации

Архитектура модели

1. Пространственная дискретизация

Используется смешанный метод конечных элементов:

• Скорость-давление: элементы Mini (P1b-P1)
• Концентрация: кусочно-линейные элементы (P1)

Пространства конечных элементов определяются как: $V_h = \{v_h \in C(\Omega)^2 \cap V | v_h|_K \in (P_1(K) \oplus b(K))^2, \forall K \in T_h\}$$M_h = \{q_h \in C(\Omega) \cap H^1(\Omega) | q_h|_K \in P_1(K), \forall K \in T_h, \int_\Omega q_h dx = 0\}$

2. Временная дискретизация: схема CNLF

Первый шаг (обратный метод Эйлера): $\frac{u_h^1 - u_h^0}{\tau} + \nu(c_h^0 + \alpha)\nabla u_h^1 + B(u_h^0, u_h^1, v_h) - (\nabla \cdot v_h, p_h^1) = \text{ПЧ}$

Последующие шаги (CNLF): $\frac{u_h^{n+1} - u_h^{n-1}}{2\tau} + A(c_h^n, \frac{u_h^{n+1} + u_h^{n-1}}{2}, v_h) + B(u_h^n, \frac{u_h^{n+1} + u_h^{n-1}}{2}, v_h) = \text{ПЧ}$

Технические инновации

1. Стратегия расщепления: полуявная обработка нелинейных членов обеспечивает расщепление уравнений скорости и концентрации
2. Временное интегрирование с прыгающей лягушкой: применяется схема Кранка-Николсона второго порядка точности для линейных членов
3. Обработка переменных коэффициентов: специально разработаны операторы проектирования для работы с концентрационно-зависимой вязкостью

Экспериментальная установка

Набор данных

Используется вычислительная область $\Omega = [0,1] \times [0,1]$ с параметрами:

• $\theta = \gamma = 1$
• Конечное время $T = 1.0$
• Аналитическое решение: $u(x,y,t) = (ye^{-t}(2y-1)(y-1), -xe^{-t}(2x-1)(x-1))^T$$p(x,y,t) = e^{-t}(2x-1)(2y-1)$$c(x,y,t) = e^{-t}\sin(\pi x)\sin(\pi y)$

Метрики оценки

• Ошибка в норме $L^2$: $\|r - r_h\|_{L^2} = \|r(t_N) - r_h^N\|_{L^2}$
• Ошибка в норме $H^1$: $\|r - r_h\|_{H^1}$
• Порядок сходимости: вычисляется путем сгущения сетки

Сравниваемые методы

Рассмотрены три различные модели вязкости:

1. $\nu = 1$ (постоянная вязкость)
2. $\nu = 1 + 0.1c$ (линейная зависимость)
3. $\nu = e^c$ (экспоненциальная зависимость)

Детали реализации

• Временной шаг: $\tau = h$
• Сгущение сетки: $h = 1/4, 1/8, 1/16, 1/32, 1/64, 1/128$
• Инструмент реализации: FreeFEM++

Результаты экспериментов

Основные результаты

Проверка устойчивости

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

Анализ сходимости

Сходимость в норме $L^2$:

• Скорость: достигает второго порядка сходимости для всех моделей вязкости
• Концентрация: достигает второго порядка сходимости
• Давление: достигает первого порядка сходимости

Конкретные численные результаты (для $\nu = 1$):

Абляционные исследования

Сравнение различных моделей вязкости подтверждает:

1. Робастность схемы CNLF для различных функций вязкости
2. Что концентрационно-зависимая вязкость не влияет на сходимость схемы
3. Согласованность теоретически предсказанных порядков сходимости с численными результатами

Экспериментальные выводы

1. Оптимальная сходимость: численные эксперименты полностью подтверждают теоретический анализ с порядком сходимости $O(\tau^2 + h^2)$
2. Робастность: схема демонстрирует хорошую устойчивость и сходимость для различных типов функций вязкости
3. Преимущества в эффективности: расщепленная схема значительно повышает вычислительную эффективность

Теоретический анализ

Анализ устойчивости

Теорема 3.1 (Безусловная устойчивость): $\|u_h^{n+1}\|_{L^2}^2 + \|c_h^{n+1}\|_{L^2}^2 + \kappa\tau\sum_{n=1}^N \|\nabla(u_h^{n+1} + u_h^{n-1})\|_{L^2}^2 \leq C$

Доказательство использует метод энергии, ключевые моменты:

1. Использование свойства кососимметричности для обработки нелинейных членов
2. Применение дискретного неравенства Гронуолла

Оценки ошибок

Теорема 4.1 (Сходимость): При условиях A1 и A2 существует константа $C$ такая, что: $\max_{0 \leq i \leq N}(\|u^i - u_h^i\|_{L^2}^2 + \|c^i - c_h^i\|_{L^2}^2) \leq C(\tau^4 + h^4)$

Доказательство использует математическую индукцию в сочетании с:

1. Оценками ошибок операторов проектирования
2. Анализом ошибок аппроксимации временной дискретизации
3. Тонкой обработкой нелинейных членов

Связанные работы

Численные методы для биоконвекции

1. Случай постоянной вязкости: 23,24 установили существование решений, 27 дали оценки ошибок конечных элементов
2. Случай переменной вязкости: 26 доказали существование и единственность слабого решения, 9 предложили схему BDF2
3. Высокопорядковые методы: 19 разработали линеаризованную схему Кранка-Николсона

Преимущества данной работы по сравнению с существующими

1. Первое применение CNLF: введение метода прыгающей лягушки в задачи биоконвекции
2. Безусловная устойчивость: в отличие от существующих методов, отсутствуют ограничения на временной шаг
3. Расщепленный дизайн: повышает вычислительную эффективность и облегчает параллельную реализацию

Заключение и обсуждение

Основные выводы

1. Успешно построена схема CNLF для задачи биоконвекции с концентрационно-зависимой вязкостью
2. Теоретически доказаны безусловная устойчивость и оптимальная сходимость схемы
3. Численные эксперименты подтверждают корректность теоретических результатов

Ограничения

1. Ограничение по размерности: теоретический анализ в основном ориентирован на двумерный случай
2. Предположения о функции вязкости: требуется липшицева непрерывность и ограниченность
3. Граничные условия: рассмотрены только однородные граничные условия Дирихле

Направления будущих исследований

Авторы предлагают расширить CNLF-фреймворк на:

1. Системы Chemotaxis-Navier-Stokes
2. Системы Patlak-Keller-Segel-Navier-Stokes
3. Системы Chemo-Repulsion-Navier-Stokes

Глубокая оценка

Преимущества

1. Теоретическая строгость: полный анализ устойчивости и сходимости с детальными доказательствами
2. Методологическая инновативность: первое применение метода CNLF к задачам биоконвекции с переменной вязкостью
3. Практическая ценность: безусловная устойчивость обеспечивает большую гибкость при практическом применении
4. Достаточная численная верификация: тестирование на различных моделях вязкости подтверждает робастность метода

Недостатки

1. Сильные теоретические предположения: высокие требования к регулярности решения могут ограничить практическую применимость
2. Отсутствие трехмерного расширения: теоретический анализ в основном ограничен двумерным случаем
3. Недостаточный анализ вычислительной сложности: отсутствует сравнение вычислительной эффективности с другими методами
4. Недостаточное обсуждение чувствительности: недостаточно изучена чувствительность метода к изменениям физических параметров

Влияние

1. Академический вклад: предоставляет новые теоретические инструменты для численных методов биоконвекции
2. Перспективы применения: имеет потенциальное применение в биоинженерии, экологии и других областях
3. Распространение методов: CNLF-фреймворк может быть применим к другим связанным системам

Области применения

1. Моделирование микробных суспензий: применимо к системам, где необходимо учитывать изменение вязкости
2. Экологическая гидродинамика: может использоваться для моделирования течений в естественных водоемах с микроорганизмами
3. Проектирование биореакторов: предоставляет численные инструменты для оптимизации конструкции биореакторов

Дополнительные технические детали

Ключевые математические инструменты

1. Кососимметричная трилинейная форма: $B(u,v,w) = \frac{1}{2}\int_\Omega (u \cdot \nabla v) \cdot w dx - \frac{1}{2}\int_\Omega (u \cdot \nabla w) \cdot v dx$
2. Оператор проектирования с переменными коэффициентами: $\nu(c)(\nabla(u-P_h^{n+1}u), \nabla v_h) + (\nabla \cdot v_h, p-\rho_h^{n+1}p) = 0$
3. Дискретное неравенство Гронуолла: ключевой инструмент для анализа устойчивости

Ключевые моменты численной реализации

1. Обработка начальных условий: первый шаг использует обратный метод Эйлера для обеспечения точности
2. Сохранение массы: выбор функциональных пространств обеспечивает сохранение общей массы микроорганизмов
3. Решение линейных систем: каждый временной шаг требует решения только линейной системы, повышая эффективность

Данная работа вносит важные вклады как в теоретическом, так и в численном аспектах, предоставляя эффективный численный метод для решения сложных задач биологических жидкостей. Несмотря на некоторые теоретические предположения и ограничения по размерности, инновационная схема CNLF и строгий анализ закладывают важную основу для исследований в соответствующих областях.