Kinetic transport equations are notoriously difficult to simulate because of their complex multiscale behaviors and the need to numerically resolve a high dimensional probability density function. Past literature has focused on building reduced order models (ROM) by analytical methods. In recent years, there is a surge of interest in developing ROM using data-driven or computational tools that offer more applicability and flexibility. This paper is a work towards that direction. Motivated by our previous work of designing ROM for the stationary radiative transfer equation in [30] by leveraging the low-rank structure of the solution manifold induced by the angular variable, we here further advance the methodology to the time-dependent model. Particularly, we take the celebrated reduced basis method (RBM) approach and propose a novel micro-macro decomposed reduced basis method (MMD-RBM). The MMD-RBM is constructed by exploiting, in a greedy fashion, the low-rank structures of both the micro- and macro-solution manifolds with respect to the angular and temporal variables. Our reduced order surrogate consists of: reduced bases for reduced order subspaces and a reduced quadrature rule in the angular space. The proposed MMD-RBM features several structure-preserving components: 1) an equilibrium-respecting strategy to construct reduced order subspaces which better utilize the structure of the decomposed system, and 2) a recipe for preserving positivity of the quadrature weights thus to maintain the stability of the underlying reduced solver. The resulting ROM can be used to achieve a fast online solve for the angular flux in angular directions outside the training set and for arbitrary order moment of the angular flux.
翻译:由于复杂的多尺度行为和需要用数字方式解决高维概率密度功能,因此难以模拟动能传输方程式。过去的文献侧重于通过分析方法建立减序模型(ROM),近年来,对利用数据驱动工具或计算工具开发可提供更多适用性和灵活性的ROM的兴趣激增。本文是朝这个方向发展的一个工作。我们以前为[30] 固定辐射传输方程式设计ROM的工作,利用由角变量引发的在线解决方案的低级别结构,我们在此进一步将方法推向基于时间的模型。特别是,我们采用经庆祝的减序法(RBM)方法,并提议采用新的微缩缩缩式降低基数方法(MMD-RBM)。MD-RBM的构建方式是,以贪婪的方式利用微缩和宏观溶解解方程式结构的低位结构结构,与矩形变量有关,我们缩小的顺序包括:降低排序的基数空间的基数基数基数基础,并降低外部的基数基数的基数,从而将稳定度结构的递减为稳定的螺旋结构。