2020/08/16 by Chenghua Duan, Duan, Chenghua, Wenbin Chen +7 · 2 citations
Computer Science · Engineering · Mathematics · #Applied mathematics #Convergence (economics) #Degenerate energy levels #Dissipation #FOS: Mathematics #Fluid Dynamics and Turbulent Flows #Fokker–Planck equation #Fractional Differential Equations Solutions #Mathematical Biology Tumor Growth #Mathematical analysis #Mathematics #Nonlinear system #Numerical Analysis (math.NA) #Numerical analysis #Numerical stability #Order of accuracy #Partial differential equation #Physics #Singularity #Smoothness #Uniqueness #cs.NA #math.NA
paper · pdf · doi:10.48550/arxiv.2008.06903
published in arXiv (Cornell University) (Cornell University)
arxiv created 2020/08/16 · openalex publication_date 2020/08/16 · arxiv updated 2020/08/18 · openalex created_date 2020/08/21 · openalex updated_date 2026/08/06
In this work, we develop novel structure-preserving numerical schemes for a class of nonlinear Fokker--Planck equations with nonlocal interactions. Such equations can cover many cases of importance, such as porous medium equations with external potentials, optimal transport problems, and aggregation-diffusion models. Based on the Energetic Variational Approach, a trajectory equation is first derived by using the balance between the maximal dissipation principle and least action principle. By a convex-splitting technique, we propose energy dissipating numerical schemes for the trajectory equation. Rigorous numerical analysis reveals that the nonlinear numerical schemes are uniquely solvable, naturally respect mass conservation and positivity at fully discrete level, and preserve steady states. Under certain smoothness assumptions, the numerical schemes are shown to be second order accurate in space and first order accurate in time. Extensive numerical simulations are performed to demonstrate several valuable features of the proposed schemes. In addition to the preservation of physical structures, such as positivity, mass conservation, discrete energy dissipation, blue and steady states, numerical simulations further reveal that our numerical schemes are capable of solving degenerate cases of the Fokker--Planck equations effectively and robustly. It is shown that the developed numerical schemes have convergence order even in degenerate cases with the presence of solutions having compact support, can accurately and robustly compute the waiting time of free boundaries without any oscillation, and can approximate blow-up singularity up to machine precision.