We propose a novel family of asymptotically stable, implicit-explicit, adaptive, time integration method (denoted with the $\theta$-method) for the solution of the fractional advection-diffusion-reaction (FADR) equations. This family of time integration method generalized the computationally explicit $L_1$-method adopted by Brunner (J. Comput. Phys. {\bf 229} 6613-6622 (2010)) as well as the fully implicit method proposed by Jannelli (Comm. Nonlin. Sci. Num. Sim., {\bf 105}, 106073 (2022)). The spectral analysis of the method (involving the group velocity and the phase speed) indicates a region of favorable dispersion for a limited range of Peclet number. The numerical inversion of the coefficient matrix is avoided by exploiting the sparse structure of the matrix in the iterative solver for the Poisson equation. The accuracy and the efficacy of the method is benchmarked using (a) the two-dimensional (2D) fractional diffusion equation, originally proposed by Brunner, and (b) the incompressible, subdiffusive dynamics of a planar viscoelastic channel flow of the Rouse chain melts (FADR equation with fractional time-derivative of order $\alpha=\nicefrac{1}{2}$) and the Zimm chain solution ($\alpha=\nicefrac{2}{3}$). Numerical simulations of the viscoelastic channel flow effectively capture the non-homogeneous regions of high viscosity at low fluid inertia (or the so-called `spatiotemporal macrostructures'), experimentally observed in the flow-instability transition of subdiffusive flows.
翻译:我们提出一个新式的组合,用零星稳定、隐含解释、适应性、适应性、时间整合方法(用美元表示,用美元表示,用美元表示,用美元表示)解决分数对冲反扩散反应(FADR)方程式的解决方案。这个时间整合方法的组合将布伦纳(J.comput.Phys.\bf 229}6613-6622(2010))采用的计算清晰的美元值统一起来,以及Jannelli(Commer.nlin.Sci.Num.Sim.,hbnalalalal-nocial2Sility3,105},106073(2022)提出的完全隐含的方法。该方法的光谱分析(涉及集团速度和阶段速度)表明,一个区域对有限范围的粒子数进行优异差。系数矩阵的翻转,通过利用调式溶液溶液溶液解方程式中的稀薄结构。该方法的准确性和功效是使用(a) 内层(demoal-deal-dealal-dealal-lial-limalial) diaal-deal dival dival romoudal dival) romodal dival romodal romodal romodal romodal romodal romodal ropal routdal rocument rocument rocument rocument rocument se rocument rocument rocumental rocument rocument rocumental rocument rocument se rodre rodre rodal rodre rout rocumental rout rout rodre rocudre rodre rodal se rocudal rocument 区域,由(观察到的两种的双平流,由(一种) rodal rodal rodal rodal rod rodal rodal rodal rodal rodaldal rodal ro) ro) ro) rodal