Qinchen Song

2papers

2 Papers

73.7NAApr 25
A Filtered MgNet Solver For Radiative Transfer Equations

Qinchen Song, Xinliang Liu, Lei Zhang

Conventional numerical solvers for the radiative transfer equation (RTE) exhibit severe sensitivity to medium parameters. To address this, we propose an operator learning framework that approximates the RTE solution map as a function of material properties. The core architecture, MgNet, preserves the solution operator framework established by recursive skeleton factorization (RSF) but substitutes its coefficient-specific sub-operators (e.g. smoother, prolongation operator and restriction operator) with learnable neural components. This design transcends the the fixed parametric structure of classical schemes, enabling data-driven sub-operator optimization and learning of their medium-parameter dependence. To mitigate spectral bias in operator learning, we introduce an adaptive angular compression technique within the loss function that dynamically suppresses high-frequency modes responsible for training instability. Comprehensive benchmarks demonstrate that, when deployed as a learned preconditioner, MgNet achieves at least 10 times acceleration over conventional preconditioners in the diffusive regime and maintains robust generalization to unseen parameter configurations. By unifying multilevel factorization structure with deep operator learning, this work establishes a physics-constrained operator-learning paradigm for radiative transport simulations.

92.8NAApr 23
Fast Algorithm For Solving Time-dependent Multiscale radiative transport Equation

Qinchen Song, Lei Zhang, Min Tang

When solving the time-dependent radiative transport equation (RTE), implicit time discretization is often employed for its robustness and stability. This results in a sequence of steady-state RTEs with identical cross-sections but varying source terms, whose repeated solution is computationally costly. To address this, we first apply the adaptive tailored finite point scheme (TFPS) for spatial discretization. This scheme exploits prior knowledge of the background media's optical properties to adaptively compress the angular domain, constructing a compressed linear system. A key feature is its ability to reconstruct the layer structure after compression, faithfully capturing the variance at the layer. We then use the Recursive Skeleton Method (RSM) to obtain an explicit multilevel decomposition of the inverse discrete operator, which is reused for all steady-state solutions. Numerical experiments show that our framework achieves high accuracy and significant efficiency across diverse scenarios.