An optimization-based positivity-preserving limiter in semi-implicit discontinuous Galerkin schemes solving Fokker-Planck equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Liu, Chen, Hu, Jingwei, Taitano, William T., Zhang, Xiangxiong
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913839563407360
author Liu, Chen
Hu, Jingwei
Taitano, William T.
Zhang, Xiangxiong
author_facet Liu, Chen
Hu, Jingwei
Taitano, William T.
Zhang, Xiangxiong
contents For high-order accurate schemes such as discontinuous Galerkin (DG) methods solving Fokker-Planck equations, it is desired to efficiently enforce positivity without losing conservation and high-order accuracy, especially for implicit time discretizations. We consider an optimization-based positivity-preserving limiter for enforcing positivity of cell averages of DG solutions in a semi-implicit time discretization scheme, so that the point values can be easily enforced to be positive by a simple scaling limiter on the DG polynomial in each cell. The optimization can be efficiently solved by a first-order splitting method with nearly optimal parameters, which has an $\mathcal{O}(N)$ computational complexity and is flexible for parallel computation. Numerical tests are shown on some representative examples to demonstrate the performance of the proposed method.
format Preprint
id arxiv_https___arxiv_org_abs_2410_19143
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle An optimization-based positivity-preserving limiter in semi-implicit discontinuous Galerkin schemes solving Fokker-Planck equations
Liu, Chen
Hu, Jingwei
Taitano, William T.
Zhang, Xiangxiong
Numerical Analysis
For high-order accurate schemes such as discontinuous Galerkin (DG) methods solving Fokker-Planck equations, it is desired to efficiently enforce positivity without losing conservation and high-order accuracy, especially for implicit time discretizations. We consider an optimization-based positivity-preserving limiter for enforcing positivity of cell averages of DG solutions in a semi-implicit time discretization scheme, so that the point values can be easily enforced to be positive by a simple scaling limiter on the DG polynomial in each cell. The optimization can be efficiently solved by a first-order splitting method with nearly optimal parameters, which has an $\mathcal{O}(N)$ computational complexity and is flexible for parallel computation. Numerical tests are shown on some representative examples to demonstrate the performance of the proposed method.
title An optimization-based positivity-preserving limiter in semi-implicit discontinuous Galerkin schemes solving Fokker-Planck equations
topic Numerical Analysis
url https://arxiv.org/abs/2410.19143