2,142 research outputs found
High-order filtered schemes for time-dependent second order HJB equations
In this paper, we present and analyse a class of "filtered" numerical schemes
for second order Hamilton-Jacobi-Bellman equations. Our approach follows the
ideas introduced in B.D. Froese and A.M. Oberman, Convergent filtered schemes
for the Monge-Amp\`ere partial differential equation, SIAM J. Numer. Anal.,
51(1):423--444, 2013, and more recently applied by other authors to stationary
or time-dependent first order Hamilton-Jacobi equations. For high order
approximation schemes (where "high" stands for greater than one), the
inevitable loss of monotonicity prevents the use of the classical theoretical
results for convergence to viscosity solutions. The work introduces a suitable
local modification of these schemes by "filtering" them with a monotone scheme,
such that they can be proven convergent and still show an overall high order
behaviour for smooth enough solutions. We give theoretical proofs of these
claims and illustrate the behaviour with numerical tests from mathematical
finance, focussing also on the use of backward difference formulae (BDF) for
constructing the high order schemes.Comment: 27 pages, 16 figures, 4 table
Adaptive monotone multigrid methods for nonlinear variational problems
A wide range of problems occurring in engineering and industry is characterized by the presence of a free (i.e. a priori unknown) boundary where the underlying physical situation is changing in a discontinuous way. Mathematically, such phenomena can be often reformulated as variational inequalities or related non–smooth minimization problems.
In these research notes, we will describe a new and promising way of constructing fast solvers for the corresponding discretized problems providing globally convergent iterative schemes with (asymptotic) multigrid
convergence speed. The presentation covers physical modelling, existence and uniqueness results, finite element approximation and adaptive mesh–refinement based on a posteriori error estimation. The numerical properties
of the resulting adaptive multilevel algorithm are illustrated by typical applications, such as semiconductor device simulation or continuous casting
- …