Source author record

Oliver Junge

Oliver Junge appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2016arXiv

A fully discrete variational scheme for solving nonlinear Fokker-Planck equations in higher space dimensions

We introduce a novel spatio-temporal discretization for nonlinear Fokker-Planck equations on the multi-dimensional unit cube. This discretization is based on two structural properties of these equations: the first is the representation as a gradient flow of an entropy functional in the $L^2$-Wasserstein metric, the second is the Lagrangian nature, meaning that solutions can be written as the push forward transformation of the initial density under suitable flow maps. The resulting numerical scheme is entropy diminishing and mass conserving. Further, the scheme is weakly stable, which allows us to prove convergence under certain regularity assumptions. Finally, we present results from numerical experiments in space dimension $d=2$.

preprint2016arXiv

Computing Coherent Sets using the Fokker-Planck Equation

We perform a numerical approximation of coherent sets in finite-dimensional smooth dynamical systems by computing singular vectors of the transfer operator for a stochastically perturbed flow. This operator is obtained by solution of a discretized Fokker-Planck equation. For numerical implementation, we employ spectral collocation methods and an exponential time differentiation scheme. We experimentally compare our approach to the more classical method by Ulam that is based on integration of the exact transfer operator.

preprint2015arXiv

On fast computation of finite-time coherent sets using radial basis functions

Finite-time coherent sets inhibit mixing over finite times. The most expensive part of the transfer operator approach to detecting coherent sets is the construction of the operator itself. We present a numerical method based on radial basis function collocation and apply it to a recent transfer operator construction that has been designed specifically for purely advective dynamics. The construction is based on a "dynamic" Laplacian operator and minimises the boundary size of the coherent sets relative to their volume. The main advantage of our new approach is a substantial reduction in the number of Lagrangian trajectories that need to be computed, leading to large speedups in the transfer operator analysis when this computation is costly.

preprint2015arXiv

Pseudo generators for under-resolved molecular dynamics

Many features of a molecule which are of physical interest (e.g. molecular conformations, reaction rates) are described in terms of its dynamics in configuration space. This article deals with the projection of molecular dynamics in phase space onto configuration space. Specifically, we study the situation that the phase space dynamics is governed by a stochastic Langevin equation and study its relation with the configurational Smoluchowski equation in the three different scaling regimes: Firstly, the Smoluchowski equations in non-Cartesian geometries are derived from the overdamped limit of the Langevin equation. Secondly, transfer operator methods are used to describe the metastable behaviour of the system at hand, and an explicit small-time asymptotics is derived on which the Smoluchowski equation turns out to govern the dynamics of the position coordinate (without any assumptions on the damping). By using an adequate reduction technique, these considerations are then extended to one-dimensional reaction coordinates. Thirdly, we sketch three different approaches to approximate the metastable dynamics based on time-local information only.

preprint2015arXiv

Pseudo generators of spatial transfer operators

Metastable behavior in dynamical systems may be a significant challenge for a simulation based analysis. In recent years, transfer operator based approaches to problems exhibiting metastability have matured. In order to make these approaches computationally feasible for larger systems, various reduction techniques have been proposed: For example, Schütte introduced a spatial transfer operator which acts on densities on configuration space, while Weber proposed to avoid trajectory simulation (like Froyland et al.) by considering a discrete generator. In this manuscript, we show that even though the family of spatial transfer operators is not a semigroup, it possesses a well defined generating structure. What is more, the pseudo generators up to order 4 in the Taylor expansion of this family have particularly simple, explicit expressions involving no momentum averaging. This makes collocation methods particularly easy to implement and computationally efficient, which in turn may open the door for further efficiency improvements in, e.g., the computational treatment of conformation dynamics. We experimentally verify the predicted properties of these pseudo generators by means of two academic examples.

preprint2014arXiv

Sparse Control of Alignment Models in High Dimension

For high dimensional particle systems, governed by smooth nonlinearities depending on mutual distances between particles, one can construct low-dimensional representations of the dynamical system, which allow the learning of nearly optimal control strategies in high dimension with overwhelming confidence. In this paper we present an instance of this general statement tailored to the sparse control of models of consensus emergence in high dimension, projected to lower dimensions by means of random linear maps. We show that one can steer, nearly optimally and with high probability, a high-dimensional alignment model to consensus by acting at each switching time on one agent of the system only, with a control rule chosen essentially exclusively according to information gathered from a randomly drawn low-dimensional representation of the control system.

preprint2012arXiv

Lazy global feedbacks for quantized nonlinear event systems

We consider nonlinear event systems with quantized state information and design a globally stabilizing controller from which only the minimal required number of control value changes along the feedback trajectory to a given initial condition is transmitted to the plant. In addition, we present a non-optimal heuristic approach which might reduce the number of control value changes and requires a lower computational effort. The constructions are illustrated by two numerical examples.

preprint2012arXiv

Optimal value functions for weakly coupled systems: a posteriori estimates

We consider weakly coupled LQ optimal control problems and derive estimates on the sensitivity of the optimal value function in dependence of the coupling strength. In order to improve these sensitivity estimates a "coupling adapted" norm is proposed. Our main result is that if a weak coupling suffices to destabilize the closed loop system with the optimal feedback of the uncoupled system then the value function might change drastically with the coupling. As a consequence, it is not reasonable to expect that a weakly coupled system possesses a weakly coupled optimal value function. Also, for a known result on the connection of the separation operator and the stability radius a new and simpler proof is given.

preprint2011arXiv

Computing the Invariant Measure and the Lyapunov Exponent for One-Dimensional Maps using a Measure-Preserving Polynomial Basis

We consider a generalisation of Ulam's method for approximating invariant densities of one-dimensional chaotic maps. Rather than use piecewise constant polynomials to approximate the density, we use polynomials of degree n which are defined by the requirement that they preserve the measure on n+1 neighbouring subintervals. Over the whole interval, this results in a discontinuous piecewise polynomial approximation to the density. We prove error results where this approach is used to approximate smooth densities. We also consider the computation of the Lyapunov exponent using the polynomial density and show that the order of convergence is one order better than for the density itself. Together with using cubic polynomials in the density approximation, this yields a very efficient method for computing highly accurate estimates of the Lyapunov exponent. We illustrate the theoretical findings with some examples.

preprint2011arXiv

Estimating long term behavior of flows without trajectory integration: the infinitesimal generator approach

The long-term distributions of trajectories of a flow are described by invariant densities, i.e. fixed points of an associated transfer operator. In addition, global slowly mixing structures, such as almost-invariant sets, which partition phase space into regions that are almost dynamically disconnected, can also be identified by certain eigenfunctions of this operator. Indeed, these structures are often hard to obtain by brute-force trajectory-based analyses. In a wide variety of applications, transfer operators have proven to be very efficient tools for an analysis of the global behavior of a dynamical system. The computationally most expensive step in the construction of an approximate transfer operator is the numerical integration of many short term trajectories. In this paper, we propose to directly work with the infinitesimal generator instead of the operator, completely avoiding trajectory integration. We propose two different discretization schemes; a cell based discretization and a spectral collocation approach. Convergence can be shown in certain circumstances. We demonstrate numerically that our approach is much more efficient than the operator approach, sometimes by several orders of magnitude.