Source author record

Jan Sieber

Jan Sieber 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

27works
10topics
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

27 published item(s)

preprint2022arXiv

Continuation with Non-invasive Control Schemes: Revealing Unstable States in a Pedestrian Evacuation Scenario

This paper presents a framework to perform bifurcation analysis in laboratory experiments or simulations. We employ control-based continuation to study the dynamics of a macroscopic variable of a microscopically defined model, exploring the potential viability of the underlying feedback control techniques in an experiment. In contrast to previous experimental studies that used iterative root-finding methods on the feedback control targets, we propose a feedback control law that is inherently non-invasive. That is, the control discovers the location of equilibria and stabilizes them simultaneously. We call the proposed control zero-in-equilibrium feedback control and we prove that it is able to stabilize branches of equilibria, except at singularities of codimension n+1, where n is the number of state space dimensions the feedback can depend on. We apply the method to a simulated evacuation scenario were pedestrians have to reach an exit after maneuvering left or right around an obstacle. The scenario shows a hysteresis phenomenon with bistability and tipping between two possible steady pedestrian flows in microscopic simulations. We demonstrate for the evacuation scenario that the proposed control law is able to uniformly discover and stabilize steady flows along the entire branch, including points where other non-invasive approaches to feedback control become singular.

preprint2022arXiv

Nonlinear effects of instantaneous and delayed state dependence in a delayed feedback loop

We study a scalar, first-order delay differential equation (DDE) with instantaneous and state-dependent delayed feedback, which itself may be delayed. The state dependence introduces nonlinearity into an otherwise linear system. We investigate the ensuing nonlinear dynamics with the case of instantaneous state dependence as our starting point. We present the bifurcation diagram in the parameter plane of the two feedback strengths showing how periodic orbits bifurcate from a curve of Hopf bifurcations and disappear along a curve where both period and amplitude grow beyond bound as the orbits become saw-tooth shaped. We then `switch on' the delay within the state-dependent feedback term, reflected by a parameter $b>0$. Our main conclusion is that the new parameter $b$ has an immediate effect: as soon as $b>0$ the bifurcation diagram for $b=0$ changes qualitatively and, specifically, the nature of the limiting saw-tooth shaped periodic orbits changes. Moreover, we show $-$ numerically and through center manifold analysis $-$ that a degeneracy at $b=1/3$ of an equilibrium with a double real eigenvalue zero leads to a further qualitative change and acts as an organizing center for the bifurcation diagram. Our results demonstrate that state dependence in delayed feedback terms may give rise to new dynamics and, moreover, that the observed dynamics may change significantly when the state-dependent feedback depends on past states of the system. This is expected to have implications for models arising in different application contexts, such as models of human balancing and conceptual climate models of delayed action oscillator type.

preprint2021arXiv

Controlling coexisting attractors of a class of non-autonomous dynamical systems

This paper studies a control method for switching stable coexisting attractors of a class of non-autonomous dynamical systems. The central idea is to introduce a continuous path for the system's trajectory to transition from its original undesired stable attractor to a desired one by varying one of the system parameters according to the information of the desired attractor. The behaviour of the control is studied numerically for both non-smooth and smooth dynamical systems, using a sof-impact and a Duffing oscillator as examples. Special attention is given to identify the regions where the proposed control strategy is applicable by using the path-following methods implemented via the continuation platform COCO. It is shown that the proposed control concept can be implemented through either using an external control input or varying a system parameter. Finally, extensive numerical results are presented to validate the proposed control methods.

preprint2020arXiv

Calculating the Lyapunov exponents of a piecewise-smooth soft impacting system with a time-delayed feedback controller

Lyapunov exponents are a widely used tool for studying dynamical systems. When calculating Lyapunov exponents for piecewise-smooth systems with time-delayed arguments one faces a lack of continuity in the variational problem. This paper studies how to build a variational equation for the efficient construction of Jacobians along trajectories of the delayed nonsmooth system. Trajectories of the piecewise-smooth system may encounter a so-called grazing event where the trajectory approaches a discontinuity surface in the state space in a non-transversal manner. For this event we develop a grazing point estimation algorithm to ensure the accuracy of trajectories for the nonlinear and the variational equations. We show that the eigenvalues of the Jacobian matrix computed by the algorithm converge with an order consistent with the order of the numerical integration method, therefore guaranteeing the reliability of our proposed numerical method. Finally, the method is demonstrated on a periodically forced impacting oscillator under the time-delayed feedback control.

preprint2016arXiv

DDE-BIFTOOL Manual - Bifurcation analysis of delay differential equations

DDEBIFTOOL is a collection of Matlab routines for numerical bifurcation analysis of systems of delay differential equations with discrete constant and state-dependent delays. The package supports continuation and stability analysis of steady state solutions and periodic solutions. Further one can compute and continue several local and global bifurcations: fold and Hopf bifurcations of steady states; folds, period doublings and torus bifurcations of periodic orbits; and connecting orbits between equilibria. To analyse the stability of steady state solutions, approximations are computed to the rightmost, stability-determining roots of the characteristic equation which can subsequently be used as starting values in a Newton procedure. For periodic solutions, approximations to the Floquet multipliers are computed. The manual describes the structure of the package, its routines, and its data and method parameter structures.

preprint2016arXiv

Early-warning indicators for rate-induced tipping

A dynamical system is said to undergo rate-induced tipping when it fails to track its quasi-equilibrium state due to an above-critical-rate change of system parameters. We study a prototypical model for rate-induced tipping, the saddle-node normal form subject to time-varying equilibrium drift and noise. We find that both most commonly used early-warning indicators, increase in variance and increase in autocorrelation, occur not when the equilibrium drift is fastest but with a delay. We explain this delay by demonstrating that the most likely trajectory for tipping also crosses the tipping threshold with a delay and therefore the tipping itself is delayed. We find solutions of the variational problem determining the most likely tipping path using numerical continuation techniques. The result is a systematic study of the most likely tipping time in the plane of two parameters, distance from tipping threshold and noise intensity.

preprint2016arXiv

Early-warning indicators in the dynamic regime

Early-warning indicators (increase of autocorrelation and variance) are commonly applied to time series data to try and detect tipping points of real-world systems. The theory behind these indicators originates from approximating the fluctuations around an equilibrium observed in time series data by a linear stationary (Ornstein-Uhlenbeck) process. Then for the approach of a bifurcation-type tipping point the formulas for the autocorrelation and variance of an Ornstein-Uhlenbeck process detect the phenomenon `critical slowing down'. The assumption of stationarity and linearity introduces two sources of error in the early-warning indicators. We investigate the difference between the theoretical and observed values for the early-warning indicators for the saddle-node normal form bifurcation with linear drift.

preprint2016arXiv

Generic stabilizability for time-delayed feedback control

Time delayed feedback control is one of the most successful methods to discover dynamically unstable features of a dynamical system in an experiment. This approach feeds back only terms that depend on the difference between the current output and the output from a fixed time T ago. Thus, any periodic orbit of period T in the feedback controlled system is also a periodic orbit of the uncontrolled system, independent of any modelling assumptions. It has been an open problem whether this approach can be successful in general, that is, under genericity conditions similar to those in linear control theory (controllability), or if there are fundamental restrictions to time-delayed feedback control. We show that there are no restrictions in principle. This paper proves the following: for every periodic orbit satisfying a genericity condition slightly stronger than classical linear controllability, one can find control gains that stabilise this orbit with extended time-delayed feedback control. While the paper's techniques are based on linear stability analysis, they exploit the specific properties of linearisations near autonomous periodic orbits in nonlinear systems, and are, thus, mostly relevant for the analysis of nonlinear experiments.

preprint2016arXiv

Relative equilibria and relative periodic solutions in systems with time-delay and $S^{1}$ symmetry

We study properties of basic solutions in systems with dime delays and $S^1$-symmetry. Such basic solutions are relative equilibria (CW solutions) and relative periodic solutions (MW solutions). It follows from the previous theory that the number of CW solutions grows generically linearly with time delay $τ$. Here we show, in particular, that the number of relative periodic solutions grows generically as $τ^2$ when delay increases. Thus, in such systems, the relative periodic solutions are more abundant than relative equilibria. The results are directly applicable to, e.g., Lang-Kobayashi model for the lasers with delayed feedback. We also study stability properties of the solutions for large delays.

preprint2015arXiv

Regular and irregular patterns of self-localized excitation in arrays of coupled phase oscillators

We study a system of phase oscillators with nonlocal coupling in a ring that supports self-organized patterns of coherence and incoherence, called chimera states. Introducing a global feedback loop, connecting the phase lag to the order parameter, we can observe chimera states also for systems with a small number of oscillators. Numerical simulations show a huge variety of regular and irregular patterns composed of localized phase slipping events of single oscillators. Using methods of classical finite dimensional chaos and bifurcation theory, we can identify the emergence of chaotic chimera states as a result of transitions to chaos via period doubling cascades, torus breakup, and intermittency. We can explain the observed phenomena by a mechanism of self-modulated excitability in a discrete excitable medium.

preprint2015arXiv

Shock-Sensitivity in Shell-Like Structures: with Simulations of Spherical Shell Buckling

Under increasing compression, an unbuckled shell is in a metastable state which becomes increasingly precarious as the buckling load is approached. So to induce premature buckling a lateral disturbance will have to overcome a decreasing energy barrier which reaches zero at buckling. Two archetypal problems that exhibit a severe form of this behaviour are the axially-compressed cylindrical shell and the externally pressurized spherical shell. Focussing on the cylinder, a non-destructive technique was recently proposed to estimate the 'shock sensitivity' of a laboratory specimen using a lateral probe to measure the nonlinear load-deflection characteristic. If a symmetry-breaking bifurcation is encountered on the path, computer simulations showed how this can be suppressed by a controlled secondary probe. Here, we extend our understanding by assessing in general terms how a single control can capture remote saddle solutions: in particular how a symmetric probe could locate an asymmetric solution. Then, more specifically, we analyse the spherical shell with point and ring probes, to test the procedure under challenging conditions to assess its range of applicability. Rather than a bifurcation, the spherical shell offers the challenge of a de-stabilizing fold (limit point) under the rigid control of the probe.

preprint2014arXiv

Bifurcation analysis of delay-induced resonances of the El-Nino Southern Oscillation

Models of global climate phenomena of low to intermediate complexity are very useful for providing an understanding at a conceptual level. An important aspect of such models is the presence of a number of feedback loops that feature considerable delay times, usually due to the time it takes to transport energy (for example, in the form of hot/cold air or water) around the globe. In this paper we demonstrate how one can perform a bifurcation analysis of the behaviour of a periodically-forced system with delay in dependence on key parameters. As an example we consider the El-Nino Southern Oscillation (ENSO), which is a sea surface temperature oscillation on a multi-year scale in the basin of the Pacific Ocean. One can think of ENSO as being generated by an interplay between two feedback effects, one positive and one negative, which act only after some delay that is determined by the speed of transport of sea-surface temperature anomalies across the Pacific. We perform here a case study of a simple delayed-feedback oscillator model for ENSO (introduced by Tziperman et al, J. Climate 11 (1998)), which is parametrically forced by annual variation. More specifically, we use numerical bifurcation analysis tools to explore directly regions of delay-induced resonances and other stability boundaries in this delay-differential equation model for ENSO.

preprint2014arXiv

Equation-Free Analysis of Macroscopic Behavior in Traffic and Pedestrian Flow

Equation-free methods make possible an analysis of the evolution of a few coarse-grained or macroscopic quantities for a detailed and realistic model with a large number of fine-grained or microscopic variables, even though no equations are explicitly given on the macroscopic level. This will facilitate a study of how the model behavior depends on parameter values including an understanding of transitions between different types of qualitative behavior. These methods are introduced and explained for traffic jam formation and emergence of oscillatory pedestrian counter flow in a corridor with a narrow door.

preprint2014arXiv

Evolution of the L1 halo family in the radial solar sail CRTBP

We present a detailed investigation of the dramatic changes that occur in the $\mathcal{L}_1$ halo family when radiation pressure is introduced into the Sun-Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. 1991) the families can be fully mapped out as the parameter $β$ is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at $β\approx0.0387$ acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at $β\approx0.289$, resulting in a third new family. The linear stability of the families changes rapidly at low values of $β$, with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter $β$ and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of $β$ there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case.

preprint2014arXiv

Implicit Methods for Equation-Free Analysis: Convergence Results and Analysis of Emergent Waves in Microscopic Traffic Models

We introduce a general formulation for an implicit equation-free method in the setting of slow-fast systems. First, we give a rigorous convergence result for equation-free analysis showing that the implicitly defined coarse-level time stepper converges to the true dynamics on the slow manifold within an error that is exponentially small with respect to the small parameter measuring time scale separation. Second, we apply this result to the idealized traffic modeling problem of phantom jams generated by cars with uniform behavior on a circular road. The traffic jams are waves that travel slowly against the direction of traffic. Equation-free analysis enables us to investigate the behavior of the microscopic traffic model on a macroscopic level. The standard deviation of cars' headways is chosen as the macroscopic measure of the underlying dynamics such that traveling wave solutions correspond to equilibria on the macroscopic level in the equation-free setup. The collapse of the traffic jam to the free flow then corresponds to a saddle-node bifurcation of this macroscopic equilibrium. We continue this bifurcation in two parameters using equation-free analysis.

preprint2014arXiv

Nonlinear dynamic interactions between flow-induced galloping and shell-like buckling

For an elastic system that is non-conservative but autonomous, subjected for example to time-independent loading by a steadily flowing fluid (air or water), a dangerous bifurcation, such as a sub-critical bifurcation, or a cyclic fold, will trigger a dynamic jump to one or more remote stable attractors. When there is more than one candidate attractor, the one onto which the structure settles can then be indeterminate, being sensitive to infinitesimally small variations in starting conditions or parameters. In this paper we develop and study an archetypal model to explore the nonlinear dynamic interactions between galloping at an incipient sub-critical Hopf bifurcation of a structure with shell-like buckling behaviour that is gravity-loaded to approach a sub-critical pitch-fork bifurcation. For the fluid forces, we draw on the aerodynamic coefficients determined experimentally by Novak for the flow around a bluff body of rectangular cross-section. Meanwhile, for the structural component, we consider a variant of the propped-cantilever model that is widely used to illustrate the sub-critical pitchfork: within this model a symmetry-breaking imperfection makes the behaviour generic. The compound bifurcation corresponding to simultaneous galloping and buckling is the so-called Takens-Bodganov Cusp. We make a full unfolding of this codimension-3 bifurcation for our archetypal model to explore the adjacent phase-space topologies and their indeterminacies.

preprint2013arXiv

A method for the reconstruction of unknown non-monotonic growth functions in the chemostat

We propose an adaptive control law that allows one to identify unstable steady states of the open-loop system in the single-species chemostat model without the knowledge of the growth function. We then show how one can use this control law to trace out (reconstruct) the whole graph of the growth function. The process of tracing out the graph can be performed either continuously or step-wise. We present and compare both approaches. Even in the case of two species in competition, which is not directly accessible with our approach due to lack of controllability, feedback control improves identifiability of the non-dominant growth rate.

preprint2013arXiv

Controlling unstable chaos: Stabilizing chimera states by feedback

We present a control scheme that is able to find and stabilize an unstable chaotic regime in a system with a large number of interacting particles. This allows us to track a high dimensional chaotic attractor through a bifurcation where it loses its attractivity. Similar to classical delayed feedback control, the scheme is non-invasive, however, only in an appropriately relaxed sense considering the chaotic regime as a statistical equilibrium displaying random fluctuations as a finite size effect. We demonstrate the control scheme for so-called chimera states, which are coherence-incoherence patterns in coupled oscillator systems. The control makes chimera states observable close to coherence, for small numbers of oscillators, and for random initial conditions.

preprint2013arXiv

Extremum seeking via continuation techniques for optimizing biogas production in the chemostat

We consider the chemostat model with the substrate concentration as the single measurement. We propose a control strategy that drives the system at a steady state maximizing the gas production without the knowledge of the specific growth rate. Our approach separates the extremum seeking problem from the feedback control problem such that each of the two subproblems can be solved with relatively simple algorithms. We are then free to choose any numerical optimization algorithm. We give a demonstration for two choices: one is based on slow-fast dynamics and numerical continuation, the other is a combination of golden-section and Newton iteration. The method copes with non-monotonic growth functions.

preprint2013arXiv

Longtime behavior of coupled wave equations for semiconductor lasers

Coupled wave equations are popular tool for investigating longitudinal dynamical effects in semiconductor lasers, for example, sensitivity to delayed optical feedback. We study a model that consists of a hyperbolic linear system of partial differential equations with one spatial dimension, which is nonlinearly coupled with a slow subsystem of ordinary differential equations. We first prove the basic statements about the existence of solutions of the initial-boundary-value problem and their smooth dependence on initial values and parameters. Hence, the model constitutes a smooth infinite-dimensional dynamical system. Then we exploit the particular slow-fast structure of the system to construct a low-dimensional attracting invariant manifold for certain parameter constellations. The flow on this invariant manifold is described by a system of ordinary differential equations that is accessible to classical bifurcation theory and numerical tools such as AUTO.

preprint2013arXiv

Systematic experimental exploration of bifurcations with non-invasive control

We present a general method for systematically investigating the dynamics and bifurcations of a physical nonlinear experiment. In particular, we show how the odd-number limitation inherent in popular non-invasive control schemes, such as (Pyragas) time-delayed or washout-filtered feedback control, can be overcome for tracking equilibria or forced periodic orbits in experiments. To demonstrate the use of our non-invasive control, we trace out experimentally the resonance surface of a periodically forced mechanical nonlinear oscillator near the onset of instability, around two saddle-node bifurcations (folds) and a cusp bifurcation.

preprint2011arXiv

Characteristic matrices for linear periodic delay differential equations

Szalai et al. (SIAM J. on Sci. Comp. 28(4), 2006) gave a general construction for characteristic matrices for systems of linear delay-differential equations with periodic coefficients. First, we show that matrices constructed in this way can have a discrete set of poles in the complex plane, which may possibly obstruct their use when determining the stability of the linear system. Then we modify and generalize the original construction such that the poles get pushed into a small neighborhood of the origin of the complex plane.

preprint2011arXiv

Nonlinear softening as a predictive precursor to climate tipping

Approaching a dangerous bifurcation, from which a dynamical system such as the Earth's climate will jump (tip) to a different state, the current stable state lies within a shrinking basin of attraction. Persistence of the state becomes increasingly precarious in the presence of noisy disturbances. We consider an underlying potential, as defined theoretically for a saddle-node fold and (via averaging) for a Hopf bifurcation. Close to a stable state, this potential has a parabolic form; but approaching a jump it becomes increasingly dominated by softening nonlinearities. If we have already detected a decrease in the linear decay rate, nonlinear information allows us to estimate the propensity for early tipping due to noise. We argue that one needs to extract information about the nonlinear features (a "softening") of the underlying potential from the time series to judge the probability and timing of tipping. This analysis is the logical next step if one has detected a decrease of the linear decay rate. If there is no discernable trend in the linear analysis, nonlinear softening is even more important in showing the proximity to tipping. After extensive normal form calibration studies, we check two geological time series from paleo-climate tipping events for softening of the underlying well. For the ending of the last ice age, where we find no convincing linear precursor, we identify a statistically significant nonlinear softening towards increasing temperature. The analysis has thus successfully detected a warning of the imminent tipping event.

preprint2011arXiv

On the Stability of a Chain of Phase Oscillators

We study a chain of $N+1$ phase oscillators with asymmetric but uniform coupling. This type of chain possesses $2^{N}$ ways to synchronize in so-called travelling wave states, i.e. states where the phases of the single oscillators are in relative equilibrium. We show that the number of unstable dimensions of a travelling wave equals the number of oscillators with relative phase close to $π$. This implies that only the relative equilibrium corresponding to approximate in-phase synchronization is locally stable. Despite the presence of a Lyapunov-type functional periodic or chaotic phase slipping occurs. For chains of length 3 and 4 we locate the region in parameter space where rotations (corresponding to phase slipping) are present.

preprint2010arXiv

Climate tipping as a noisy bifurcation: a predictive technique

It is often known, from modelling studies, that a certain mode of climate tipping (of the oceanic thermohaline circulation, for example) is governed by an underlying fold bifurcation. For such a case we present a scheme of analysis that determines the best stochastic fit to the existing data. This provides the evolution rate of the effective control parameter, the variation of the stability coefficient, the path itself and its tipping point. By assessing the actual effective level of noise in the available time series, we are then able to make probability estimates of the time of tipping. This new technique is applied, first, to the output of a computer simulation for the end of greenhouse Earth about 34 million years ago when the climate tipped from a tropical state into an icehouse state with ice caps. Second, we use the algorithms to give probabilistic tipping estimates for the end of the most recent glaciation of the Earth using actual archaeological ice-core data.

preprint2010arXiv

Small-scale instabilities in dynamical systems with sliding

We demonstrate with a minimal example that in Filippov systems (dynamical systems governed by discontinuous but piecewise smooth vector fields) stable periodic motion with sliding is not robust with respect to stable singular perturbations. We consider a simple dynamical system that we assume to be a quasi-static approximation of a higher-dimensional system containing a fast stable subsystem. We tune a system parameter such that a stable periodic orbit of the simple system touches the discontinuity surface: this is the so-called grazing-sliding bifurcation. The periodic orbit remains stable, and its local return map becomes piecewise linear. However, when we take into account the fast dynamics the local return map of the periodic orbit changes qualitatively, giving rise to, for example, period-adding cascades or small-scale chaos.