Optimal control using flux potentials: A way to construct bound-preserving finite element schemes for conservation laws
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel pith:7VEFAORHrecord.jsonopen to challenge →
read the original abstract
To ensure preservation of local or global bounds for numerical solutions of conservation laws, we constrain a baseline finite element discretization using optimization-based (OB) flux correction. The main novelty of the proposed methodology lies in the use of flux potentials as control variables and targets of inequality-constrained optimization problems for numerical fluxes. In contrast to optimal control via general source terms, the discrete conservation property of flux-corrected finite element approximations is guaranteed without the need to impose additional equality constraints. Since the number of flux potentials is less than the number of fluxes in the multidimensional case, the potential-based version of optimal flux control involves fewer unknowns than direct calculation of optimal fluxes. We show that the feasible set of a potential-state potential-target (PP) optimization problem is nonempty and choose a primal-dual Newton method for calculating the optimal flux potentials. The results of numerical studies for linear advection and anisotropic diffusion problems in 2D demonstrate the superiority of the new OB-PP algorithms to closed-form flux limiting under worst-case assumptions.
This paper has not been read by Pith yet.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.