2025-11-22T04:13:15.983918

A Particle-Field Algorithm with Neural Interpolation for a Parabolic-Hyperbolic Chemotaxis System in 3D

Kim, Xin

Tumor angiogenesis involves a collection of tumor cells moving towards blood vessels for nutrients to grow. Angiogenesis, and in general chemo- taxis, systems have been modeled using partial differential equations (PDEs) and as such require numerical methods to approximate their solutions. Here we study a Parabolic-Hyperbolic Keller-Segel (PHKS) system in three space dimensions. The model arises in the angiogenesis literature. To compute solutions to the PHKS system, we develop a neural stochastic interacting particle-field (NSIPF) method where the density variable is represented as empirical measures of particles and the field variable (concentration of chemoattractant) approximated by a convolutional neural network (CNN). We discuss the performance of NSIPF in computing multi-bump solutions to the system.

academic

A Particle-Field Algorithm with Neural Interpolation for a Parabolic-Hyperbolic Chemotaxis System in 3D

基本信息

论文ID: 2510.13199
标题: A Particle-Field Algorithm with Neural Interpolation for a Parabolic-Hyperbolic Chemotaxis System in 3D
作者: Jongwon David Kim, Jack Xin (University of California, Irvine)
分类: math.NA cs.NA
发表时间: 2025年10月15日 (arXiv预印本)
论文链接: https://arxiv.org/abs/2510.13199v1

摘要

肿瘤血管生成涉及肿瘤细胞向血管移动以获取营养物质的过程。血管生成和一般的趋化性系统已通过偏微分方程(PDEs)建模，因此需要数值方法来近似求解。本文研究三维空间中的抛物-双曲Keller-Segel (PHKS)系统，该模型源于血管生成文献。为求解PHKS系统，作者开发了神经随机相互作用粒子场(NSIPF)方法，其中密度变量表示为粒子的经验测度，场变量(趋化因子浓度)通过卷积神经网络(CNN)近似。文章讨论了NSIPF在计算系统多峰解方面的性能。

研究背景与动机

要解决的问题:
- 肿瘤血管生成中细胞向血管移动的数学建模与数值求解问题
- 三维空间中抛物-双曲Keller-Segel系统的高效数值计算
问题重要性:
- 血管生成是肿瘤转移的关键生物过程，理解其机制对癌症治疗具有重要意义
- 趋化性系统在细菌觅食、组织形成等多个生物过程中起重要作用
现有方法局限性:
- 传统有限差分法(FDM)在三维高分辨率计算中耗时较长
- 经典样条插值方法在三维空间中实现复杂且计算成本高
- 缺乏专门针对抛物-双曲系统的高效粒子方法
研究动机:
- 利用神经网络的高维逼近能力改进插值效率
- 结合粒子方法的无网格优势和CNN的插值性能
- 为三维趋化性系统提供更高效的数值求解方案

核心贡献

提出了神经随机相互作用粒子场(NSIPF)方法，将CNN插值器集成到随机粒子场算法中
建立了PHKS系统的混沌传播理论基础，通过适度相互作用为神经插值提供理论支撑
开发了高效的CNN架构，使用径向对称解训练数据，成功实现2D和3D系统的插值
实现了显著的计算效率提升，相比传统方法在保持解质量的同时大幅降低计算成本

方法详解

任务定义

求解三维空间中的抛物-双曲Keller-Segel系统：

ρₜ = ∇ · (γ∇ρ - χρ∇c)  (1)
cₜ = -cρ                  (2)

其中ρ是细菌密度，c是趋化因子浓度，γ是扩散系数，χ是趋化敏感性。

模型架构

1. 粒子表示

密度函数通过粒子的经验测度近似：

ρₜ ≈ (M₀/P) Σⱼ₌₁ᴾ δ(x - Xⱼₜ)

其中P是粒子数，M₀是守恒质量。

2. SDE表述

粒子运动遵循随机微分方程：

dXᵢₜ = χ∇c(Xᵢₜ,t)dt + √(2γ)dBᵢₜ

3. CNN插值器架构

输入层: 低分辨率浓度场数据
卷积层: 5层3D卷积，通道数为1→16→32→32→32→16→1
激活函数: ReLU
跳跃连接: 增强梯度传播
输出: 高分辨率插值场

4. NSIPF算法流程

Algorithm 1: Neural SIPF
1. 初始化粒子位置基于ρ₀
2. For each time step:
   a. 将粒子分箱得到ρₙ
   b. 更新浓度: cⁿ⁺¹ = cⁿ - δt·cⁿ·ρⁿ
   c. CNN插值得到c(x,tₙ)
   d. 计算∇c并更新粒子位置