Source author record

Péter Koltai

Péter Koltai 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

13works
5topics
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

13 published item(s)

preprint2022arXiv

Data-driven modelling of nonlinear dynamics by barycentric coordinates and memory

We present a numerical method to model dynamical systems from data. We use the recently introduced method Scalable Probabilistic Approximation (SPA) to project points from a Euclidean space to convex polytopes and represent these projected states of a system in new, lower-dimensional coordinates denoting their position in the polytope. We then introduce a specific nonlinear transformation to construct a model of the dynamics in the polytope and to transform back into the original state space. To overcome the potential loss of information from the projection to a lower-dimensional polytope, we use memory in the sense of the delay-embedding theorem of Takens. By construction, our method produces stable models. We illustrate the capacity of the method to reproduce even chaotic dynamics and attractors with multiple connected components on various examples.

preprint2022arXiv

Large-scale flow in a cubic Rayleigh-Bénard cell: Long-term turbulence statistics and Markovianity of macrostate transitions

We investigate the large-scale circulation (LSC) in a turbulent Rayleigh-Bénard convection flow in a cubic closed convection cell by means of direct numerical simulations at a Rayleigh number $Ra=10^6$. The numerical studies are conducted for single flow trajectories up to $10^5$ convective free-fall times to obtain a sufficient sampling of the four discrete LSC states, which can be summarized to one macrostate, and the two crossover configurations which are taken by the flow in between for short periods. We find that large-scale dynamics depends strongly on the Prandtl number $Pr$ of the fluid which has values of 0.1, 0.7, and 10. Alternatively, we run an ensemble of 3600 short-term direct numerical simulations to study the transition probabilities between the discrete LSC states. This second approach is also used to probe the Markov property of the dynamics. Our ensemble analysis gave strong indication of Markovianity of the transition process from one LSC state to another, even though the data are still accompanied by considerable noise. It is based on the eigenvalue spectrum of the transition probability matrix, further on the distribution of persistence times and the joint distribution of two successive macrostate persistence times.

preprint2020arXiv

Computation and optimal perturbation of finite-time coherent sets for aperiodic flows without trajectory integration

Understanding the macroscopic behavior of dynamical systems is an important tool to unravel transport mechanisms in complex flows. A decomposition of the state space into coherent sets is a popular way to reveal this essential macroscopic evolution. To compute coherent sets from an aperiodic time-dependent dynamical system we consider the relevant transfer operators and their infinitesimal generators on an augmented space-time manifold. This space-time generator approach avoids trajectory integration, and creates a convenient linearization of the aperiodic evolution. This linearization can be further exploited to create a simple and effective spectral optimization methodology for diminishing or enhancing coherence. We obtain explicit solutions for these optimization problems using Lagrange multipliers and illustrate this technique by increasing and decreasing mixing of spatial regions through small velocity field perturbations.

preprint2020arXiv

Dimensionality Reduction of Complex Metastable Systems via Kernel Embeddings of Transition Manifolds

We present a novel kernel-based machine learning algorithm for identifying the low-dimensional geometry of the effective dynamics of high-dimensional multiscale stochastic systems. Recently, the authors developed a mathematical framework for the computation of optimal reaction coordinates of such systems that is based on learning a parametrization of a low-dimensional transition manifold in a certain function space. In this article, we enhance this approach by embedding and learning this transition manifold in a reproducing kernel Hilbert space, exploiting the favorable properties of kernel embeddings. Under mild assumptions on the kernel, the manifold structure is shown to be preserved under the embedding, and distortion bounds can be derived. This leads to a more robust and more efficient algorithm compared to previous parametrization approaches.

preprint2020arXiv

Memory-based reduced modelling and data-based estimation of opinion spreading

We investigate opinion dynamics based on an agent-based model, and are interested in predicting the evolution of the percentages of the entire agent population that share an opinion. Since these opinion percentages can be seen as an aggregated observation of the full system state, the individual opinions of each agent, we view this in the framework of the Mori-Zwanzig projection formalism. More specifically, we show how to estimate a nonlinear autoregressive model (NAR) with memory from data given by a time series of opinion percentages, and discuss its prediction capacities for various specific topologies of the agent interaction network. We demonstrate that the inclusion of memory terms significantly improves the prediction quality on examples with different network topologies.

preprint2016arXiv

On the numerical approximation of the Perron-Frobenius and Koopman operator

Information about the behavior of dynamical systems can often be obtained by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with a dynamical system. Examples of such operators are the Perron-Frobenius and the Koopman operator. In this paper, we will review different methods that have been developed over the last decades to compute finite-dimensional approximations of these infinite-dimensional operators - e.g. Ulam's method and Extended Dynamic Mode Decomposition (EDMD) - and highlight the similarities and differences between these approaches. The results will be illustrated using simple stochastic differential equations and molecular dynamics examples.

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

Split optimal policy iteration for LQR problems

This technical report is concerned with the convergence properties of what we call the split optimal policy iteration for coupled LQR problems; see section 3.1 in the manuscript. Interestingly, the iteration shows different convergence behavior for continuous and discrete time systems: while global convergence holds for both cases, we have local quadratic convergence for the continuous time case, but only linear convergence for the discrete time case - even though quadratic convergence is retained in the limit as the coupling between the subsystems vanishes.

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.

preprint2012arXiv

Optimizing the stable behavior of parameter-dependent dynamical systems - maximal domains of attraction, minimal absorption times

We propose a method for approximating solutions to optimization problems involving the global stability properties of parameter-dependent continuous-time autonomous dynamical systems. The method relies on an approximation of the infinite-state deterministic system by a finite-state non-deterministic one - a Markov jump process. The key properties of the method are that it does not use any trajectory simulation, and that the parameters and objective function are in a simple (and except for a system of linear equations) explicit relationship.

preprint2011arXiv

Discrete infinitesimal generator of the Frobenius-Perron operator semigroup associated with "outflow systems"

In this technical report the $C_0$ semigroup of Frobenius-Perron operators on $L^1(X)$ is considered, where the underlying dynamical system is such that trajectories may leave the state space $X$ and terminate. We introduce a discrete infinitesimal generator and show, that the operator semigroup generated by this discrete generator converges in $L^1(X)$ pointwise to the Frobenius-Perron operator of the system.

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.