Maximum-principle preserving space-time isogeometric analysis
This work provides a provably stable and high-order accurate method for convection-dominated problems, which is important for computational fluid dynamics and transport simulations.
The authors propose a nonlinear stabilization technique for convection-diffusion-reaction and pure transport problems using space-time isogeometric analysis, which unconditionally satisfies the discrete maximum principle for arbitrary order and is linearity preserving. The method is validated with numerical experiments in 1D and 2D, showing improved nonlinear convergence and reduced computational cost via a partitioned scheme.
In this work we propose a nonlinear stabilization technique for convection-diffusion-reaction and pure transport problems discretized with space-time isogeometric analysis. The stabilization is based on a graph-theoretic artificial diffusion operator and a novel shock detector for isogeometric analysis. Stabilization in time and space directions are performed similarly, which allow us to use high-order discretizations in time without any CFL-like condition. The method is proven to yield solutions that satisfy the discrete maximum principle (DMP) unconditionally for arbitrary order. In addition, the stabilization is linearity preserving in a space-time sense. Moreover, the scheme is proven to be Lipschitz continuous ensuring that the nonlinear problem is well-posed. Solving large problems using a space-time discretization can become highly costly. Therefore, we also propose a partitioned space-time scheme that allows us to select the length of every time slab, and solve sequentially for every subdomain. As a result, the computational cost is reduced while the stability and convergence properties of the scheme remain unaltered. In addition, we propose a twice differentiable version of the stabilization scheme, which enjoys the same stability properties while the nonlinear convergence is significantly improved. Finally, the proposed schemes are assessed with numerical experiments. In particular, we considered steady and transient pure convection and convection-diffusion problems in one and two dimensions.