Karsten Ahnert

Odeint – Solving ordinary differential equations in C++

Many physical, biological or chemical systems are modeled by ordinary differential equations (ODEs) and finding their solution is an every-day-task for many scientists. Here, we introduce a new C++ library dedicated to find numerical solutions of initial value problems of ODEs: odeint (www.odeint.com). odeint is implemented in a highly generic way and provides extensive interoperability at top performance. For example, due to its modular design it can be easily parallized with OpenMP and even runs on CUDA GPUs. Despite that, it provides a convenient interface that allows for a simple and easy usage.

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.

Strong and Weak Chaos in Weakly Nonintegrable Many-Body Hamiltonian Systems

We study properties of chaos in generic one-dimensional nonlinear Hamiltonian lattices comprised of weakly coupled nonlinear oscillators by numerical simulations of continuous-time systems and symplectic maps. For small coupling, the measure of chaos is found to be proportional to the coupling strength and lattice length, with the typical maximal Lyapunov exponent being proportional to the square root of coupling. This strong chaos appears as a result of triplet resonances between nearby modes. In addition to strong chaos we observe a weakly chaotic component having much smaller Lyapunov exponent, the measure of which drops approximately as a square of the coupling strength down to smallest couplings we were able to reach. We argue that this weak chaos is linked to the regime of fast Arnold diffusion discussed by Chirikov and Vecheslavov. In disordered lattices of large size we find a subdiffusive spreading of initially localized wave packets over larger and larger number of modes. The relations between the exponent of this spreading and the exponent in the dependence of the fast Arnold diffusion on coupling strength are analyzed. We also trace parallels between the slow spreading of chaos and deterministic rheology.

Dynamical Thermalization of Disordered Nonlinear Lattices

We study numerically how the energy spreads over a finite disordered nonlinear one-dimensional lattice, where all linear modes are exponentially localized by disorder. We establish emergence of dynamical thermalization characterized as an ergodic chaotic dynamical state with a Gibbs distribution over the modes. Our results show that the fraction of thermalizing modes is finite and grows with the nonlinearity strength.

Scaling of energy spreading in strongly nonlinear disordered lattices

To characterize a destruction of Anderson localization by nonlinearity, we study the spreading behavior of initially localized states in disordered, strongly nonlinear lattices. Due to chaotic nonlinear interaction of localized linear or nonlinear modes, energy spreads nearly subdiffusively. Based on a phenomenological description by virtue of a nonlinear diffusion equation, we establish a one-parameter scaling relation between the velocity of spreading and the density, which is confirmed numerically. From this scaling it follows that for very low densities the spreading slows down compared to the pure power law.

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.

Traveling waves and compactons in phase oscillator lattices

We study waves in a chain of dispersively coupled phase oscillators. Two approaches--a quasicontinuous approximation and an iterative numerical solution of the lattice equation--allow us to characterize different types of traveling waves: compactons, kovatons, solitary waves with exponential tails as well as a novel type of semicompact waves that are compact from one side. Stability of these waves is studied using numerical simulations of the initial value problem.

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.

Metaprogramming Applied to Numerical Problems

From the discovery that the template system of C++ forms a Turing complete language in 1994, a programming technique called Template Metaprogramming has emerged that allows for the creation of faster, more generic and better code. Here, we apply Template Metaprogramming to implement a generic Runge-Kutta scheme that can be used to numerically solve ordinary differential equations. We show that using Template Metaprogramming results in a significantly improved performance compared to a classical implementation.