Markus Abel

Numerical differentiation of experimental data: local versus global methods

Abstract In the context of the analysis of measured data, one is often faced with the task to differentiate data numerically. Typically, this occurs when measured data are concerned or data are evaluated numerically during the evolution of partial or ordinary differential equations. Usually, one does not take care for accuracy of the resulting estimates of derivatives because modern computers are assumed to be accurate to many digits. But measurements yield intrinsic errors, which are often much less accurate than the limit of the machine used, and there exists the effect of “loss of significance”, well known in numerical mathematics and computational physics. The problem occurs primarily in numerical subtraction, and clearly, the estimation of derivatives involves the approximation of differences. In this article, we discuss several techniques for the estimation of derivatives. As a novel aspect, we divide into local and global methods, and explain the respective shortcomings. We have developed a general scheme for global methods, and illustrate our ideas by spline smoothing and spectral smoothing. The results from these less known techniques are confronted with the ones from local methods. As typical for the latter, we chose Savitzky–Golay-filtering and finite differences. Two basic quantities are used for characterization of results: The variance of the difference of the true derivative and its estimate, and as important new characteristic, the smoothness of the estimate. We apply the different techniques to numerically produced data and demonstrate the application to data from an aeroacoustic experiment. As a result, we find that global methods are generally preferable if a smooth process is considered. For rough estimates local methods work acceptably well.

Cluster-based reduced-order modelling of a mixing layer

We propose a novel cluster-based reduced-order modelling (CROM) strategy of unsteady flows. CROM combines the cluster analysis pioneered in Gunzburgers group (Burkardt et al. 2006) and and transition matrix models introduced in fluid dynamics in Eckhardts group (Schneider et al. 2007). CROM constitutes a potential alternative to POD models and generalises the Ulam-Galerkin method classically used in dynamical systems to determine a finite-rank approximation of the Perron-Frobenius operator. The proposed strategy processes a time-resolved sequence of flow snapshots in two steps. First, the snapshot data are clustered into a small number of representative states, called centroids, in the state space. These centroids partition the state space in complementary non-overlapping regions (centroidal Voronoi cells). Departing from the standard algorithm, the probabilities of the clusters are determined, and the states are sorted by analysis of the transition matrix. Secondly, the transitions between the states are dynamically modelled using a Markov process. Physical mechanisms are then distilled by a refined analysis of the Markov process, e.g. using finite-time Lyapunov exponent and entropic methods. This CROM framework is applied to the Lorenz attractor (as illustrative example), to velocity fields of the spatially evolving incompressible mixing layer and the three-dimensional turbulent wake of a bluff body. For these examples, CROM is shown to identify non-trivial quasi-attractors and transition processes in an unsupervised manner. CROM has numerous potential applications for the systematic identification of physical mechanisms of complex dynamics, for comparison of flow evolution models, for the identification of precursors to desirable and undesirable events, and for flow control applications exploiting nonlinear actuation dynamics.

Soft capacitors for wave energy harvesting

Wave energy harvesting could be a substantial renewable energy source without impact on the global climate and ecology, yet practical attempts have struggled with the problems of wear and catastrophic failure. An innovative technology for ocean wave energy harvesting was recently proposed, based on the use of soft capacitors. This study presents a realistic theoretical and numerical model for the quantitative characterization of this harvesting method. Parameter regions with optimal behavior are found, and novel material descriptors are determined, which dramatically simplify analysis. The characteristics of currently available materials are evaluated, and found to merit a very conservative estimate of 10 years for raw material cost recovery.

Synchronization of Sound Sources

Sound generation and interaction are highly complex, nonlinear, and self-organized. Nearly 150 years ago Rayleigh raised the following problem: two nearby organ pipes of different fundamental frequencies sound together almost inaudibly with identical pitch. This effect is now understood qualitatively by modern synchronization theory M. Abel et al. [J. Acoust. Soc. Am. 119, 2467 (2006)10.1121/1.2170441]. For a detailed investigation, we substituted one pipe by an electric speaker. We observe that even minute driving signals force the pipe to synchronization, thus yielding three decades of synchronization-the largest range ever measured to our knowledge. Furthermore, a mutual silencing of the pipe is found, which can be explained by self-organized oscillations, of use for novel methods of noise abatement. Finally, we develop a reconstruction method which yields a perfect quantitative match of experiment and theory.

Prediction of dynamical systems by symbolic regression

We study the modeling and prediction of dynamical systems based on conventional models derived from measurements. Such algorithms are highly desirable in situations where the underlying dynamics are hard to model from physical principles or simplified models need to be found. We focus on symbolic regression methods as a part of machine learning. These algorithms are capable of learning an analytically tractable model from data, a highly valuable property. Symbolic regression methods can be considered as generalized regression methods. We investigate two particular algorithms, the so-called fast function extraction which is a generalized linear regression algorithm, and genetic programming which is a very general method. Both are able to combine functions in a certain way such that a good model for the prediction of the temporal evolution of a dynamical system can be identified. We illustrate the algorithms by finding a prediction for the evolution of a harmonic oscillator based on measurements, by detecting an arriving front in an excitable system, and as a real-world application, the prediction of solar power production based on energy production observations at a given site together with the weather forecast.

Closed-loop separation control over a sharp edge ramp using genetic programming

Abstract We experimentally perform open and closed-loop control of a separating turbulent boundary layer downstream from a sharp edge ramp. The turbulent boundary layer just above the separation point has a Reynolds number

Maximum entropy analysis of hydraulic pipe networks

Re_{\theta }\approx 3500

Localization in a coupled standard map lattice

Reθ≈3500 based on momentum thickness. The goal of the control is to mitigate separation and early re-attachment. The forcing employs a spanwise array of active vortex generators. The flow state is monitored with skin-friction sensors downstream of the actuators. The feedback control law is obtained using model-free genetic programming control (GPC) (Gautier et al. in J Fluid Mech 770:442–457, 2015). The resulting flow is assessed using the momentum coefficient, pressure distribution and skin friction over the ramp and stereo PIV. The PIV yields vector field statistics, e.g. shear layer growth, the back-flow area and vortex region. GPC is benchmarked against the best periodic forcing. While open-loop control achieves separation reduction by locking-on the shedding mode, GPC gives rise to similar benefits by accelerating the shear layer growth. Moreover, GPC uses less actuation energy.