-
Cooperation in a non-Ergodic World on a Network -- Insurance and Beyond
Authors:
Tobias Wand,
Oliver Kamps,
Benjamin Skjold
Abstract:
Cooperation between individuals is emergent in all parts of society, yet mechanistic reasons for this emergence is ill understood in the literature. A specific example of this is insurance. Recent work has, though, shown that assuming the risk individuals face is proportional to their wealth and optimising the time average growth rate rather than the ensemble average results in a non-zero-sum game…
▽ More
Cooperation between individuals is emergent in all parts of society, yet mechanistic reasons for this emergence is ill understood in the literature. A specific example of this is insurance. Recent work has, though, shown that assuming the risk individuals face is proportional to their wealth and optimising the time average growth rate rather than the ensemble average results in a non-zero-sum game, where both parties benefit from cooperation through insurance contracts. In a recent paper, Peters and Skjold present a simple agent-based model and show how, over time, agents that enter into such cooperatives outperform agents that do not. Here, we extend this work by restricting the possible connections between agents via a lattice network. Under these restrictions, we still find that all agents profit from cooperating through insurance. We, though, further find that clusters of poor and rich agents emerge endogenously on the two-dimensional map and that wealth inequalities persist for a long duration, consistent with the phenomenon known as the poverty trap. By tuning the parameters which control the risk levels, we simulate both highly advantageous and extremely risky gambles and show that despite the qualitative shift in the type of risk, the findings are consistent.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Anticipation of Oligocene's climate heartbeat by simplified eigenvalue estimation
Authors:
Martin Heßler,
Oliver Kamps
Abstract:
The Eocene-Oligocene transition marks a watershed point of earth's climate history. The climate shifts from a greenhouse state to an icehouse state in which Antarctica glaciated for the first time and periodic dynamics arise which are still relevant for our current climate. We analyse a $CaCO_3$ concentration time series which covers the Eocene-Oligocene transition and which is obtained from a Pac…
▽ More
The Eocene-Oligocene transition marks a watershed point of earth's climate history. The climate shifts from a greenhouse state to an icehouse state in which Antarctica glaciated for the first time and periodic dynamics arise which are still relevant for our current climate. We analyse a $CaCO_3$ concentration time series which covers the Eocene-Oligocene transition and which is obtained from a Pacific sediment core at site DSDP1218. Therefore, we introduce a simplified autoregression-based variant of the dominant eigenvalue (DEV) estimation procedure. The DEV works as leading indicator of bifurcation-induced transitions and enables us to identify the bifurcation type. We confirm its reliability in a methodological study and demonstrate the crucial importance of proper detrending to obtain unbiased results. As a remark, we discuss also possible pathways to estimate the stability of limit cycles based on the DEV and the alternative drift slope as a proof of principle. Finally, we present the DEV analysis results of the $CaCO_3$ concentration time series which are reproducible in a wide parameter range. Our findings demonstrate that the onset of Oligocene's periodic dynamics might be announced by a Neimark-Sacker/Hopf bifurcation in course of the Eocene-Oligocene transition 34 mya. (We follow the convention and use mya$\widehat{=}$"million years ago" and Ma$\widehat{=}$"million years" throughout the article.)
△ Less
Submitted 25 September, 2023;
originally announced September 2023.
-
Estimating Stable Fixed Points and Langevin Potentials for Financial Dynamics
Authors:
Tobias Wand,
Timo Wiedemann,
Jan Harren,
Oliver Kamps
Abstract:
The Geometric Brownian Motion (GBM) is a standard model in quantitative finance, but the potential function of its stochastic differential equation (SDE) cannot include stable nonzero prices. This article generalises the GBM to an SDE with polynomial drift of order q and shows via model selection that q=2 is most frequently the optimal model to describe the data. Moreover, Markov chain Monte Carlo…
▽ More
The Geometric Brownian Motion (GBM) is a standard model in quantitative finance, but the potential function of its stochastic differential equation (SDE) cannot include stable nonzero prices. This article generalises the GBM to an SDE with polynomial drift of order q and shows via model selection that q=2 is most frequently the optimal model to describe the data. Moreover, Markov chain Monte Carlo ensembles of the accompanying potential functions show a clear and pronounced potential well, indicating the existence of a stable price.
△ Less
Submitted 28 November, 2023; v1 submitted 21 September, 2023;
originally announced September 2023.
-
Efficient Multi-Change Point Analysis to decode Economic Crisis Information from the S&P500 Mean Market Correlation
Authors:
Martin Heßler,
Tobias Wand,
Oliver Kamps
Abstract:
Identifying macroeconomic events that are responsible for dramatic changes of economy is of particular relevance to understand the overall economic dynamics. We introduce an open-source available efficient Python implementation of a Bayesian multi-trend change point analysis which solves significant memory and computing time limitations to extract crisis information from a correlation metric. Ther…
▽ More
Identifying macroeconomic events that are responsible for dramatic changes of economy is of particular relevance to understand the overall economic dynamics. We introduce an open-source available efficient Python implementation of a Bayesian multi-trend change point analysis which solves significant memory and computing time limitations to extract crisis information from a correlation metric. Therefore, we focus on the recently investigated S&P500 mean market correlation in a period of roughly 20 years that includes the dot-com bubble, the global financial crisis and the Euro crisis. The analysis is performed two-fold: first, in retrospect on the whole dataset and second, in an on-line adaptive manner in pre-crisis segments. The on-line sensitivity horizon is roughly determined to be 80 up to 100 trading days after a crisis onset. A detailed comparison to global economic events supports the interpretation of the mean market correlation as an informative macroeconomic measure by a rather good agreement of change point distributions and major crisis events. Furthermore, the results hint to the importance of the U.S. housing bubble as trigger of the global financial crisis, provide new evidence for the general reasoning of locally (meta)stable economic states and could work as a comparative impact rating of specific economic events.
△ Less
Submitted 25 July, 2023;
originally announced August 2023.
-
Memory Effects, Multiple Time Scales and Local Stability in Langevin Models of the S&P500 Market Correlation
Authors:
Tobias Wand,
Martin Heßler,
Oliver Kamps
Abstract:
The analysis of market correlations is crucial for optimal portfolio selection of correlated assets, but their memory effects have often been neglected. In this work, we analyse the mean market correlation of the S&P500 which corresponds to the main market mode in principle component analysis. We fit a generalised Langevin equation (GLE) to the data whose memory kernel implies that there is a sign…
▽ More
The analysis of market correlations is crucial for optimal portfolio selection of correlated assets, but their memory effects have often been neglected. In this work, we analyse the mean market correlation of the S&P500 which corresponds to the main market mode in principle component analysis. We fit a generalised Langevin equation (GLE) to the data whose memory kernel implies that there is a significant memory effect in the market correlation ranging back at least three trading weeks. The memory kernel improves the forecasting accuracy of the GLE compared to models without memory and hence, such a memory effect has to be taken into account for optimal portfolio selection to minimise risk or for predicting future correlations. Moreover, a Bayesian resilience estimation provides further evidence for non-Markovianity in the data and suggests the existence of a hidden slow time scale that operates on much slower times than the observed daily market data. Assuming that such a slow time scale exists, our work supports previous research on the existence of locally stable market states.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Non-linear, bivariate stochastic modelling of power-grid frequency applied to islands
Authors:
Ulrich Oberhofer,
Leonardo Rydin Gorjão,
G. Cigdem Yalcin,
Oliver Kamps,
Veit Hagenmeyer,
Benjamin Schäfer
Abstract:
Mitigating climate change requires a transition away from fossil fuels towards renewable energy. As a result, power generation becomes more volatile and options for microgrids and islanded power-grid operation are being broadly discussed. Therefore, studying the power grids of physical islands, as a model for islanded microgrids, is of particular interest when it comes to enhancing our understandi…
▽ More
Mitigating climate change requires a transition away from fossil fuels towards renewable energy. As a result, power generation becomes more volatile and options for microgrids and islanded power-grid operation are being broadly discussed. Therefore, studying the power grids of physical islands, as a model for islanded microgrids, is of particular interest when it comes to enhancing our understanding of power-grid stability. In the present paper, we investigate the statistical properties of the power-grid frequency of three island systems: Iceland, Ireland, and the Balearic Islands. We utilise a Fokker-Planck approach to construct stochastic differential equations that describe market activities, control, and noise acting on power-grid dynamics. Using the obtained parameters we create synthetic time series of the frequency dynamics. Our main contribution is to propose two extensions of stochastic power-grid frequency models and showcase the applicability of these new models to non-Gaussian statistics, as encountered in islands.
△ Less
Submitted 22 September, 2023; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Quantifying Tipping Risks in Power Grids and beyond
Authors:
Martin Heßler,
Oliver Kamps
Abstract:
Critical transitions, ubiquitous in nature and technology, necessitate anticipation to avert adverse outcomes. While many studies focus on bifurcation-induced tipping, where a control parameter change leads to destabilization, alternative scenarios are conceivable, e.g. noise-induced tipping by an increasing noise level in a multi-stable system. Although the generating mechanisms can be different,…
▽ More
Critical transitions, ubiquitous in nature and technology, necessitate anticipation to avert adverse outcomes. While many studies focus on bifurcation-induced tipping, where a control parameter change leads to destabilization, alternative scenarios are conceivable, e.g. noise-induced tipping by an increasing noise level in a multi-stable system. Although the generating mechanisms can be different, the observed time series can exhibit similar characteristics. Therefore, we propose a Bayesian Langevin approach, implemented in an open-source tool, which is capable of quantifying both deterministic and intrinsic stochastic dynamics simultaneously. After a detailed proof of concept, we analyse two bus voltage frequency time series of the historic North America Western Interconnection blackout on 10th August 1996. Our results unveil the intricate interplay of changing resilience and noise influence. A comparison with the blackout's timeline supports our frequency dynamics' Langevin model, with the BL-estimation indicating a permanent grid state change already two minutes before the officially defined triggering event. A tree-related high impedance fault or sudden load increases may serve as earlier triggers during this event, as suggested by our findings. This study underscores the importance of distinguishing destabilizing factors for a reliable anticipation of critical transitions, offering a tool for better understanding such events across various disciplines.
△ Less
Submitted 10 November, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
Identifying Dominant Industrial Sectors in Market States of the S&P 500 Financial Data
Authors:
Tobias Wand,
Martin Heßler,
Oliver Kamps
Abstract:
Understanding and forecasting changing market conditions in complex economic systems like the financial market is of great importance to various stakeholders such as financial institutions and regulatory agencies. Based on the finding that the dynamics of sector correlation matrices of the S&P 500 stock market can be described by a sequence of distinct states via a clustering algorithm, we try to…
▽ More
Understanding and forecasting changing market conditions in complex economic systems like the financial market is of great importance to various stakeholders such as financial institutions and regulatory agencies. Based on the finding that the dynamics of sector correlation matrices of the S&P 500 stock market can be described by a sequence of distinct states via a clustering algorithm, we try to identify the industrial sectors dominating the correlation structure of each state. For this purpose, we use a method from Explainable Artificial Intelligence (XAI) on daily S&P 500 stock market data from 1992 to 2012 to assign relevance scores to every feature of each data point. To compare the significance of the features for the entire data set we develop an aggregation procedure and apply a Bayesian change point analysis to identify the most significant sector correlations. We show that the correlation matrix of each state is dominated only by a few sector correlations. Especially the energy and IT sector are identified as key factors in determining the state of the economy. Additionally we show that a reduced surrogate model, using only the eight sector correlations with the highest XAI-relevance, can replicate 90% of the cluster assignments. In general our findings imply an additional dimension reduction of the dynamics of the financial market.
△ Less
Submitted 29 March, 2023; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Efficient Bayesian estimation of a non-Markovian Langevin model driven by correlated noise
Authors:
Clemens Willers,
Oliver Kamps
Abstract:
Data-driven modeling of non-Markovian dynamics is a recent topic of research with applications in many fields such as climate research, molecular dynamics, biophysics, or wind power modeling. In the frequently used standard Langevin equation, memory effects can be implemented through an additional hidden component which functions as correlated noise, thus resulting in a non-Markovian model. It can…
▽ More
Data-driven modeling of non-Markovian dynamics is a recent topic of research with applications in many fields such as climate research, molecular dynamics, biophysics, or wind power modeling. In the frequently used standard Langevin equation, memory effects can be implemented through an additional hidden component which functions as correlated noise, thus resulting in a non-Markovian model. It can be seen as part of the model class of partially observed diffusions which are usually adapted to observed data via Bayesian estimation, whereby the difficulty of the unknown noise values is solved through a Gibbs sampler. However, when regarding large data sets with a length of $10^6$ or $10^7$ data points, sampling the distribution of the same amount of latent variables is unfeasible. For the model discussed in this work, we solve this issue through a direct derivation of the posterior distribution of the Euler-Maruyama approximation of the model via analytical marginalization of the latent variables. Yet, in the case of a nonlinear noise process, the inverse problem of model estimation proves to be ill-posed and still numerically expensive. We handle these complications by restricting the noise to an Ornstein-Uhlenbeck process, which considerably reduces the ambiguity of the estimation. Further, in this case, the estimation can be performed very efficiently if the drift and diffusion functions of the observed component are approximated in a piecewise constant manner. We illustrate the resulting procedure of efficient Bayesian estimation of the considered non-Markovian Langevin model by an example from turbulence.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
Quantifying resilience and the risk of regime shifts under strong correlated noise
Authors:
Martin Heßler,
Oliver Kamps
Abstract:
Early warning indicators often suffer from the shortness and coarse-graining of real-world time series. Furthermore, the typically strong and correlated noise contributions in real applications are severe drawbacks for statistical measures. Even under favourable simulation conditions the measures are of limited capacity due to their qualitative nature and sometimes ambiguous trend-to-noise ratio.…
▽ More
Early warning indicators often suffer from the shortness and coarse-graining of real-world time series. Furthermore, the typically strong and correlated noise contributions in real applications are severe drawbacks for statistical measures. Even under favourable simulation conditions the measures are of limited capacity due to their qualitative nature and sometimes ambiguous trend-to-noise ratio. In order to solve these shortcomings, we analyse the stability of the system via the slope of the deterministic term of a Langevin equation, which is hypothesized to underlie the system dynamics close to the fixed point. The open-source available method is applied to a previously studied seasonal ecological model under noise levels and correlation scenarios commonly observed in real world data. We compare the results to autocorrelation, standard deviation, skewness and kurtosis as leading indicator candidates by a Bayesian model comparison with a linear and a constant model. We show that the slope of the deterministic term is a promising alternative due to its quantitative nature and high robustness against noise levels and types. The commonly computed indicators apart from the autocorrelation with deseasonalization fail to provide reliable insights into the stability of the system in contrast to a previously performed study in which the standard deviation was found to perform best. In addition, we discuss the significant influence of the seasonal nature of the data to the robust computation of the various indicators, before we determine approximately the minimal amount of data per time window that leads to significant trends for the drift slope estimations.
△ Less
Submitted 7 November, 2022; v1 submitted 7 April, 2022;
originally announced April 2022.
-
Bayesian on-line anticipation of critical transitions
Authors:
Martin Heßler,
Oliver Kamps
Abstract:
The design of reliable indicators to anticipate critical transitions in complex systems is an im portant task in order to detect a coming sudden regime shift and to take action in order to either prevent it or mitigate its consequences. We present a data-driven method based on the estimation of a parameterized nonlinear stochastic differential equation that allows for a robust anticipation of crit…
▽ More
The design of reliable indicators to anticipate critical transitions in complex systems is an im portant task in order to detect a coming sudden regime shift and to take action in order to either prevent it or mitigate its consequences. We present a data-driven method based on the estimation of a parameterized nonlinear stochastic differential equation that allows for a robust anticipation of critical transitions even in the presence of strong noise levels like they are present in many real world systems. Since the parameter estimation is done by a Markov Chain Monte Carlo approach we have access to credibility bands allowing for a better interpretation of the reliability of the results. By introducing a Bayesian linear segment fit it is possible to give an estimate for the time horizon in which the transition will probably occur based on the current state of information. This approach is also able to handle nonlinear time dependencies of the parameter controlling the transition. In general the method could be used as a tool for on-line analysis to detect changes in the resilience of the system and to provide information on the probability of the occurrence of a critical transition in future.
△ Less
Submitted 11 August, 2021;
originally announced August 2021.
-
Efficient Bayesian estimation of the generalized Langevin equation from data
Authors:
Clemens Willers,
Oliver Kamps
Abstract:
Modeling non-Markovian time series is a recent topic of research in many fields such as climate modeling, biophysics, molecular dynamics, or finance. The generalized Langevin equation (GLE), given naturally by the Mori-Zwanzig projection formalism, is a frequently used model including memory effects. In applications, a specific form of the GLE is most often obtained on a data-driven basis. Here, B…
▽ More
Modeling non-Markovian time series is a recent topic of research in many fields such as climate modeling, biophysics, molecular dynamics, or finance. The generalized Langevin equation (GLE), given naturally by the Mori-Zwanzig projection formalism, is a frequently used model including memory effects. In applications, a specific form of the GLE is most often obtained on a data-driven basis. Here, Bayesian estimation has the advantage of providing both suitable model parameters and their credibility in a straightforward way. It can be implemented in the approximating case of white noise, which, far from thermodynamic equilibrium, is consistent with the fluctuation-dissipation theorem. However, the exploration of the posterior, which is done via Markov chain Monte Carlo sampling, is numerically expensive, which makes the analysis of large data sets unfeasible. In this work, we discuss an efficient implementation of Bayesian estimation of the GLE based on a piecewise constant approximation of the drift and diffusion functions of the model. In this case, the characteristics of the data are represented by only a few coefficients, so that the numerical cost of the procedure is significantly reduced and independent of the length of the data set. Further, we propose a modification of the memory term of the GLE, leading to an equivalent model with an emphasis on the impact of trends, which ensures that an estimate of the standard Langevin equation provides an effective initial guess for the GLE. We illustrate the capabilities of both the method and the model by an example from turbulence.
△ Less
Submitted 22 July, 2022; v1 submitted 9 July, 2021;
originally announced July 2021.
-
Non-parametric estimation of a Langevin model driven by correlated noise
Authors:
Clemens Willers,
Oliver Kamps
Abstract:
Langevin models are frequently used to model various stochastic processes in different fields of natural and social sciences. They are adapted to measured data by estimation techniques such as maximum likelihood estimation, Markov chain Monte Carlo methods, or the non-parametric direct estimation method introduced by Friedrich et al. The latter has the distinction of being very effective in the co…
▽ More
Langevin models are frequently used to model various stochastic processes in different fields of natural and social sciences. They are adapted to measured data by estimation techniques such as maximum likelihood estimation, Markov chain Monte Carlo methods, or the non-parametric direct estimation method introduced by Friedrich et al. The latter has the distinction of being very effective in the context of large data sets. Due to their $δ$-correlated noise, standard Langevin models are limited to Markovian dynamics. A non-Markovian Langevin model can be formulated by introducing a hidden component that realizes correlated noise. For the estimation of such a partially observed diffusion a different version of the direct estimation method was introduced by Lehle et al. However, this procedure includes the limitation that the correlation length of the noise component is small compared to that of the measured component. In this work we propose another version of the direct estimation method that does not include this restriction. Via this method it is possible to deal with large data sets of a wider range of examples in an effective way. We discuss the abilities of the proposed procedure using several synthetic examples.
△ Less
Submitted 4 March, 2021;
originally announced March 2021.
-
Adaptive stochastic continuation with a modified lifting procedure applied to complex systems
Authors:
Clemens Willers,
Uwe Thiele,
Andrew J. Archer,
David J. B. Lloyd,
Oliver Kamps
Abstract:
Many complex systems occurring in the natural or social sciences or economics are frequently described on a microscopic level, e.g., by lattice- or agent-based models. To analyze the states of such systems and their bifurcation structure on the level of macroscopic observables, one has to rely on equation-free methods like stochastic continuation. Here, we investigate how to improve stochastic con…
▽ More
Many complex systems occurring in the natural or social sciences or economics are frequently described on a microscopic level, e.g., by lattice- or agent-based models. To analyze the states of such systems and their bifurcation structure on the level of macroscopic observables, one has to rely on equation-free methods like stochastic continuation. Here, we investigate how to improve stochastic continuation techniques by adaptively choosing the parameters of the algorithm. This allows one to obtain bifurcation diagrams quite accurately, especially near bifurcation points. We introduce lifting techniques which generate microscopic states with a naturally grown structure, which can be crucial for a reliable evaluation of macroscopic quantities. We show how to calculate fixed points of fluctuating functions by employing suitable linear fits. This procedure offers a simple measure of the statistical error. We demonstrate these improvements by applying the approach in analyses of (i) the Ising model in two dimensions, (ii) an active Ising model, and (iii) a stochastic Swift-Hohenberg model. We conclude by discussing the abilities and remaining problems of the technique.
△ Less
Submitted 6 October, 2020; v1 submitted 5 February, 2020;
originally announced February 2020.
-
Heterogeneities in electricity grids strongly enhance non-Gaussian features of frequency fluctuations under stochastic power input
Authors:
Matthias F. Wolff,
Katrin Schmietendorf,
Pedro G. Lind,
Oliver Kamps,
Joachim Peinke,
Philipp Maass
Abstract:
Stochastic feed-in of fluctuating renewable energies is steadily increasing in modern electricity grids and this becomes an important risk factor for maintaining power grid stability. Here we study the impact of wind power feed-in on the short-term frequency fluctuations in power grids based on an IEEE test grid structure, the swing equation for the dynamics of voltage phase angles, and a series o…
▽ More
Stochastic feed-in of fluctuating renewable energies is steadily increasing in modern electricity grids and this becomes an important risk factor for maintaining power grid stability. Here we study the impact of wind power feed-in on the short-term frequency fluctuations in power grids based on an IEEE test grid structure, the swing equation for the dynamics of voltage phase angles, and a series of measured wind speed data. External control measures are accounted for by adjusting the grid state to the average power feed-in on a time scale of one minute. The wind power is injected at a single node by replacing one of the conventional generator nodes in the test grid by a wind farm. We determine histograms of local frequencies for a large number of one-minute wind speed sequences taken from the measured data and for different injection nodes. These histograms exhibit a common type of shape, which can be described by a Gaussian distribution for small frequencies and a nearly exponentially decaying tail part. Non-Gaussian features become particularly pronounced for wind power injection at locations, which are weakly connected to the main grid structure. This effect is only present when taking into account the heterogeneities in transmission line and node properties of the grid, while it disappears upon homogenizing of these features. The standard deviation of the frequency fluctuations increases linearly with the average injected wind power.
△ Less
Submitted 1 November, 2019; v1 submitted 10 August, 2019;
originally announced August 2019.
-
Bridging between Load-Flow and Kuramoto-like Power Grid Models: A Flexible Approach to Integrating Electrical Storage Units
Authors:
Katrin Schmietendorf,
Oliver Kamps,
Matthias Wolff,
Pedro G. Lind,
Philipp Maass,
Joachim Peinke
Abstract:
In future power systems, electrical storage will be the key technology for balancing feed-in fluctuations. With increasing share of renewables and reduction of system inertia, the focus of research expands towards short-term grid dynamics and collective phenomena. Against this backdrop, Kuramoto-like power grids have been established as a sound mathematical modeling framework bridging between the…
▽ More
In future power systems, electrical storage will be the key technology for balancing feed-in fluctuations. With increasing share of renewables and reduction of system inertia, the focus of research expands towards short-term grid dynamics and collective phenomena. Against this backdrop, Kuramoto-like power grids have been established as a sound mathematical modeling framework bridging between the simplified models from nonlinear dynamics and the more detailed models used in electrical engineering. However, they have a blind spot concerning grid components, which cannot be modeled by oscillator equations, and hence do not allow to investigate storage-related issues from scratch. We remove this shortcoming by bringing together Kuramoto-like and algebraic load-flow equations. This is a substantial extension of the current Kuramoto framework with arbitrary grid components. Based on this concept, we provide a solid starting point for the integration of flexible storage units enabling to address current problems like smart storage control, optimal siting and rough cost estimations. For demonstration purpose, we here consider a wind power application with realistic feed-in conditions. We show how to implement basic control strategies from electrical engineering, give insights into their potential with respect to frequency quality improvement and point out their limitations by maximum capacity and finite-time response.
△ Less
Submitted 5 December, 2018;
originally announced December 2018.
-
The footprint of atmospheric turbulence in power grid frequency measurements
Authors:
Hauke Haehne,
Jannik Schottler,
Matthias Waechter,
Joachim Peinke,
Oliver Kamps
Abstract:
Fluctuating wind energy makes a stable grid operation challenging. Due to the direct contact with atmospheric turbulence, intermittent short-term variations in the wind speed are converted to power fluctuations that cause transient imbalances in the grid. We investigate the impact of wind energy feed-in on short-term fluctuations in the frequency of the public power grid, which we have measured in…
▽ More
Fluctuating wind energy makes a stable grid operation challenging. Due to the direct contact with atmospheric turbulence, intermittent short-term variations in the wind speed are converted to power fluctuations that cause transient imbalances in the grid. We investigate the impact of wind energy feed-in on short-term fluctuations in the frequency of the public power grid, which we have measured in our local distribution grid. By conditioning on wind power production data, provided by the ENTSO-E transparency platform, we demonstrate that wind energy feed-in has a measurable effect on frequency increment statistics for short time scales (< 1 sec) that are below the activation time of frequency control. Our results are in accordance with previous numerical studies of self-organized synchronization in power grids under intermittent perturbation and rise new challenges for a stable operation of future power grids fed by a high share of renewable generation.
△ Less
Submitted 13 March, 2018; v1 submitted 2 February, 2018;
originally announced February 2018.
-
On the stability and quality of power grids subjected to intermittent feed-in
Authors:
Katrin Schmietendorf,
Joachim Peinke,
Oliver Kamps
Abstract:
Feed-in fluctuations induced by renewables are one of the key challenges to the stability and quality of electrical power grids. In particular short-term fluctuations disturb the system on a time scale, on which load balancing does not operate yet and the system is intrinsically governed by self-organized synchronization. Wind and solar power are known to be strongly non-Gaussian with intermittent…
▽ More
Feed-in fluctuations induced by renewables are one of the key challenges to the stability and quality of electrical power grids. In particular short-term fluctuations disturb the system on a time scale, on which load balancing does not operate yet and the system is intrinsically governed by self-organized synchronization. Wind and solar power are known to be strongly non-Gaussian with intermittent increment statistics in these time scales. We investigate the impact of short-term wind fluctuations on the basis of a Kuramoto-like power grid model considering stability in terms of desynchronization and frequency and voltage quality aspects. We compare intermittent feed-in with a realistic power spectrum, correlated Gaussian noise of the same spectrum, and Gaussian white noise. We found out that the likelihood of severe outages is mainly determined by the temporal correlation of the feed-in. The intermittent nature of wind power is transferred into frequency and voltage fluctuations. This establishes a novel type of fluctuations with severe consequences on frequency and voltage quality, which are beyond engineering status of knowledge.
△ Less
Submitted 24 November, 2016;
originally announced November 2016.
-
Self-Organized Synchronization and Voltage Stability in Networks of Synchronous Machines
Authors:
Katrin Schmietendorf,
Joachim Peinke,
Rudolf Friedrich,
Oliver Kamps
Abstract:
The integration of renewable energy sources in the course of the energy transition is accompanied by grid decentralization and fluctuating power feed-in characteristics. This raises new challenges for power system stability and design. We intend to investigate power system stability from the viewpoint of self-organized synchronization aspects. In this approach, the power grid is represented by a n…
▽ More
The integration of renewable energy sources in the course of the energy transition is accompanied by grid decentralization and fluctuating power feed-in characteristics. This raises new challenges for power system stability and design. We intend to investigate power system stability from the viewpoint of self-organized synchronization aspects. In this approach, the power grid is represented by a network of synchronous machines. We supplement the classical Kuramoto-like network model, which assumes constant voltages, with dynamical voltage equations, and thus obtain an extended version, that incorporates the coupled categories voltage stability and rotor angle synchronization. We compare disturbance scenarios in small systems simulated on the basis of both classical and extended model and we discuss resultant implications and possible applications to complex modern power grids.
△ Less
Submitted 10 July, 2013;
originally announced July 2013.
-
The Lundgren-Monin-Novikov Hierarchy: Kinetic Equations for Turbulence
Authors:
R. Friedrich,
A. Daitche,
O. Kamps,
J. Lülff,
M. Voßkuhle,
M. Wilczek
Abstract:
We present an overview of recent works on the statistical description of turbulent flows in terms of probability density functions (PDFs) in the framework of the Lundgren-Monin-Novikov (LMN) hierarchy. Within this framework, evolution equations for the PDFs are derived from the basic equations of fluid motion. The closure problem arises either in terms of a coupling to multi-point PDFs or in terms…
▽ More
We present an overview of recent works on the statistical description of turbulent flows in terms of probability density functions (PDFs) in the framework of the Lundgren-Monin-Novikov (LMN) hierarchy. Within this framework, evolution equations for the PDFs are derived from the basic equations of fluid motion. The closure problem arises either in terms of a coupling to multi-point PDFs or in terms of conditional averages entering the evolution equations as unknown functions. We mainly focus on the latter case and use data from direct numerical simulations (DNS) to specify the unclosed terms. Apart from giving an introduction into the basic analytical techniques, applications to two-dimensional vorticity statistics, to the single-point velocity and vorticity statistics of three-dimensional turbulence, to the temperature statistics of Rayleigh-Bénard convection and to Burgers turbulence are discussed.
△ Less
Submitted 22 October, 2012; v1 submitted 28 September, 2012;
originally announced September 2012.
-
Two-Point Vorticity Statistics in the Inverse Turbulent Cascade
Authors:
R. Friedrich,
M. Voßkuhle,
O. Kamps,
M. Wilczek
Abstract:
We present a statistical analysis of the two-point vorticity probability density of the vorticity field generated in the inverse cascade of stationary two-dimensional turbulence.
We present a statistical analysis of the two-point vorticity probability density of the vorticity field generated in the inverse cascade of stationary two-dimensional turbulence.
△ Less
Submitted 11 December, 2012; v1 submitted 15 December, 2010;
originally announced December 2010.
-
An exact relation between Eulerian and Lagrangian velocity increment statistics
Authors:
O. Kamps,
R. Friedrich,
R. Grauer
Abstract:
We present a formal connection between Lagrangian and Eulerian velocity increment distributions which is applicable to a wide range of turbulent systems ranging from turbulence in incompressible fluids to magnetohydrodynamic turbulence. For the case of the inverse cascade regime of two-dimensional turbulence we numerically estimate the transition probabilities involved in this connection. In thi…
▽ More
We present a formal connection between Lagrangian and Eulerian velocity increment distributions which is applicable to a wide range of turbulent systems ranging from turbulence in incompressible fluids to magnetohydrodynamic turbulence. For the case of the inverse cascade regime of two-dimensional turbulence we numerically estimate the transition probabilities involved in this connection. In this context we are able to directly identify the processes leading to strongly non-Gaussian statistics for the Lagrangian velocity increments.
△ Less
Submitted 1 October, 2008; v1 submitted 25 September, 2008;
originally announced September 2008.
-
Lagrangian statistics in forced two-dimensional turbulence
Authors:
Oliver Kamps,
Rudolf Friedrich
Abstract:
We report on simulations of two-dimensional turbulence in the inverse energy cascade regime. Focusing on the statistics of Lagrangian tracer particles, scaling behavior of the probability density functions of velocity fluctuations is investigated. The results are compared to the three-dimensional case. In particular an analysis in terms of compensated cumulants reveals the transition from a stro…
▽ More
We report on simulations of two-dimensional turbulence in the inverse energy cascade regime. Focusing on the statistics of Lagrangian tracer particles, scaling behavior of the probability density functions of velocity fluctuations is investigated. The results are compared to the three-dimensional case. In particular an analysis in terms of compensated cumulants reveals the transition from a strong non-Gaussian behavior with large tails to Gaussianity. The reported computation of correlation functions for the acceleration components sheds light on the underlying dynamics of the tracer particles.
△ Less
Submitted 9 October, 2007;
originally announced October 2007.
-
Lagrangian Investigation of Two-Dimensional Decaying Turbulence
Authors:
Michael Wilczek,
Oliver Kamps,
Rudolf Friedrich
Abstract:
We present a numerical investigation of two-dimensional decaying turbulence in the Lagrangian framework. Focusing on single particle statistics, we investigate Lagrangian trajectories in a freely evolving turbulent velocity field. The dynamical evolution of the tracer particles is strongly dominated by the emergence and evolution of coherent structures. For a statistical analysis we focus on the…
▽ More
We present a numerical investigation of two-dimensional decaying turbulence in the Lagrangian framework. Focusing on single particle statistics, we investigate Lagrangian trajectories in a freely evolving turbulent velocity field. The dynamical evolution of the tracer particles is strongly dominated by the emergence and evolution of coherent structures. For a statistical analysis we focus on the Lagrangian acceleration as a central quantity. For more geometrical aspects we investigate the curvature along the trajectories. We find strong signatures for self-similar universal behavior.
△ Less
Submitted 5 September, 2007;
originally announced September 2007.
-
A Corrsin type approximation for Lagrangian fluid Turbulence
Authors:
Rudolf Friedrich,
Rainer Grauer,
Holger Homann,
Oliver Kamps
Abstract:
In Lagrangian turbulence one is faced with the puzzle that 2D
Navier-Stokes flows are nearly as intermittent as in three dimensions although no intermittency is present in the inverse cascade in 2D Eulerian turbulence. In addition, an inertial range is very difficult to detect and it is questionable whether it exists at all. Here, we investigate the transition of Eulerian to Lagrangian probabi…
▽ More
In Lagrangian turbulence one is faced with the puzzle that 2D
Navier-Stokes flows are nearly as intermittent as in three dimensions although no intermittency is present in the inverse cascade in 2D Eulerian turbulence. In addition, an inertial range is very difficult to detect and it is questionable whether it exists at all. Here, we investigate the transition of Eulerian to Lagrangian probability density functions (PDFs) which leads to a new type of Lagrangian structure function. This possesses an extended inertial range similar to the case of tracer particles in a frozen turbulent velocity field. This allows a connection to the scaling of Eulerian transversal structure functions.
△ Less
Submitted 22 May, 2007;
originally announced May 2007.