1,170 research outputs found
A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure
We propose a positivity preserving entropy decreasing finite volume scheme
for nonlinear nonlocal equations with a gradient flow structure. These
properties allow for accurate computations of stationary states and long-time
asymptotics demonstrated by suitably chosen test cases in which these features
of the scheme are essential. The proposed scheme is able to cope with
non-smooth stationary states, different time scales including metastability, as
well as concentrations and self-similar behavior induced by singular nonlocal
kernels. We use the scheme to explore properties of these equations beyond
their present theoretical knowledge
Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization
In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it depends on a shock detector. Further, the resulting method is linearity preserving. The same shock detector is used to gradually lump the mass matrix. The resulting method is LED, positivity preserving, and also satisfies a global DMP. Lipschitz continuity has also been proved. However, the resulting scheme is highly nonlinear, leading to very poor nonlinear convergence rates. We propose a smooth version of the scheme, which leads to twice differentiable nonlinear stabilization schemes. It allows one to straightforwardly use Newton’s method and obtain quadratic convergence. In the numerical experiments, steady and transient linear transport, and transient Burgers’ equation have been considered in 2D. Using the Newton method with a smooth version of the scheme we can reduce 10 to 20 times the number of iterations of Anderson acceleration with the original non-smooth scheme. In any case, these properties are only true for the converged solution, but not for iterates. In this sense, we have also proposed the concept of projected nonlinear solvers, where a projection step is performed at the end of every nonlinear iterations onto a FE space of admissible solutions. The space of admissible solutions is the one that satisfies the desired monotonic properties (maximum principle or positivity).Peer ReviewedPostprint (author's final draft
The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations
The cutoff method, which cuts off the values of a function less than a given
number, is studied for the numerical computation of nonnegative solutions of
parabolic partial differential equations. A convergence analysis is given for a
broad class of finite difference methods combined with cutoff for linear
parabolic equations. Two applications are investigated, linear anisotropic
diffusion problems satisfying the setting of the convergence analysis and
nonlinear lubrication-type equations for which it is unclear if the convergence
analysis applies. The numerical results are shown to be consistent with the
theory and in good agreement with existing results in the literature. The
convergence analysis and applications demonstrate that the cutoff method is an
effective tool for use in the computation of nonnegative solutions. Cutoff can
also be used with other discretization methods such as collocation, finite
volume, finite element, and spectral methods and for the computation of
positive solutions.Comment: 19 pages, 41 figure
- …