Search

We investigate flow-driven amoeboid motility as exhibited by microplasmodia of Physarum polycephalum. A poroelastic two-phase model with rigid boundaries is extended to the case of free boundaries and substrate friction. The cytoskeleton is modeled as an active viscoelastic solid permeated by a fluid phase describing the cytosol. A feedback loop between a chemical regulator, active mechanical deformations, and induced flows gives rise to oscillatory and irregular motion accompanied by spatio-temporal contraction patterns. We cover extended parameter regimes of active tension and substrate friction by numerical simulations in one spatial dimension and reproduce experimentally observed oscillation periods and amplitudes. In line with experiments, the model predicts alternating forward and backward ectoplasmatic flow at the boundaries with reversed flow in the center. However, for all cases of periodic and irregular motion, we observe practically no net motion. A simple theoretical argument shows that directed motion is not possible with a spatially independent substrate friction.

Traveling localized spots represent an important class of self-organized two-dimensional patterns in reaction-diffusion systems. We study open-loop control intended to guide a stable spot along a desired trajectory with desired velocity. Simultaneously, the spot's concentration profile does not change under control. For a given protocol of motion, we first express the control signal analytically in terms of the Goldstone modes and the propagation velocity of the uncontrolled spot. Thus, detailed information about the underlying nonlinear reaction kinetics is unnecessary. Then, we confirm the optimality of this solution by demonstrating numerically its equivalence to the solution of a regularized, optimal control problem. To solve the latter, the analytical expressions for the control are excellent initial guesses speeding-up substantially the otherwise time-consuming calculations., Comment: 21 pages, 8 figures , 6 videos

We investigate optimal control of dynamical systems which are affine, i.e., linear in control, but nonlinear in state. The control task is to enforce the system state to follow a prescribed desired trajectory as closely as possible, a task also known as optimal trajectory tracking. To obtain well-behaved solutions to optimal control, a regularization term with coefficient $\varepsilon$ must be included in the cost functional. Assuming $\varepsilon$ to be small, we reinterpret affine optimal control problems as singularly perturbed differential equations. Performing a singular perturbation expansion, approximations for the optimal tracking of arbitrary desired trajectories are derived. For $\varepsilon=0$, the state trajectory may become discontinuous, and the control may diverge. On the other hand, the analytical treatment becomes exact. We identify the conditions leading to linear evolution equations. These result in exact analytical solutions for an entire class of nonlinear trajectory tracking problems. The class comprises, among others, mechanical control systems in one spatial dimension and the FitzHugh-Nagumo model with a control acting on the activator., Comment: 5 pages, 1 figure, 11 pages appendix

Trajectory tracking of nonlinear dynamical systems with affine open-loop controls is investigated. The control task is to enforce the system state to follow a prescribed desired trajectory as closely as possible. We introduce exactly realizable desired trajectories as these trajectories which can be tracked exactly by an appropriate control. Exactly realizable trajectories are characterized mathematically by means of Moore-Penrose projectors constructed from the input matrix. The approach leads to differential-algebraic systems of equations and is considerably simpler than the related concept of system inversion. Furthermore, we identify a particularly simple class of nonlinear affine control systems. Systems in this class satisfy the so-called linearizing assumption and share many properties with linear control systems. For example, conditions for controllability can be formulated in terms of a rank condition for a controllability matrix analogously to the Kalman rank condition for linear time-invariant systems., Comment: 10 pages, no figures

This thesis investigates optimal trajectory tracking of nonlinear dynamical systems with affine controls. The control task is to enforce the system state to follow a prescribed desired trajectory as closely as possible. The concept of so-called exactly realizable trajectories is proposed. For exactly realizable desired trajectories exists a control signal which enforces the state to exactly follow the desired trajectory. For a given affine control system, these trajectories are characterized by the so-called constraint equation. This approach does not only yield an explicit expression for the control signal in terms of the desired trajectory, but also identifies a particularly simple class of nonlinear control systems. Based on that insight, the regularization parameter is used as the small parameter for a perturbation expansion. This results in a reinterpretation of affine optimal control problems with small regularization term as singularly perturbed differential equations. The small parameter originates from the formulation of the control problem and does not involve simplifying assumptions about the system dynamics. Combining this approach with the linearizing assumption, approximate and partly linear equations for the optimal trajectory tracking of arbitrary desired trajectories are derived. For vanishing regularization parameter, the state trajectory becomes discontinuous and the control signal diverges. On the other hand, the analytical treatment becomes exact and the solutions are exclusively governed by linear differential equations. Thus, the possibility of linear structures underlying nonlinear optimal control is revealed. This fact enables the derivation of exact analytical solutions to an entire class of nonlinear trajectory tracking problems with affine controls. This class comprises mechanical control systems in one spatial dimension and the FitzHugh-Nagumo model., Comment: 240 pages, 36 figures, PhD thesis

This work deals with the position control of selected patterns in reaction-diffusion systems. Exemplarily, the Schl\"{o}gl and FitzHugh-Nagumo model are discussed using three different approaches. First, an analytical solution is proposed. Second, the standard optimal control procedure is applied. The third approach extends standard optimal control to so-called sparse optimal control that results in very localized control signals and allows the analysis of second order optimality conditions., Comment: 22 pages, 3 figures, 2 tables

In two-dimensional reaction-diffusion systems, local curvature perturbations in the shape of traveling waves are typically damped out and disappear in the course of time. If, however, the inhibitor diffuses much faster than the activator, transversal instabilities can arise, leading from flat to folded, spatio-temporally modulated wave shapes and to spreading spiral turbulence. For experimentally relevant parameter values, the photosensitive Belousov-Zhabotinsky reaction (PBZR) does not exhibit transversal wave instabilities. Here, we propose a mechanism to artificially induce these instabilities via a wave shape dependent spatio-temporal feedback loop, and study the emerging wave patterns. In numerical simulations with the modified Oregonator model for the PBZR using experimentally realistic parameter values we demonstrate the feasibility of this control scheme. Conversely, in a piecewise-linear version of the FitzHugh-Nagumo model transversal instabilities and spiral turbulence in the uncontrolled system are shown to be suppressed in the presence of control, thereby stabilising flat wave propagation., Comment: 22 pages, 12 figures

We present a method to control the two-dimensional shape of traveling wave solutions to reaction-diffusion systems, as e.g. interfaces and excitation pulses. Control signals that realize a pre-given wave shape are determined analytically from nonlinear evolution equation for isoconcentration lines as the perturbed nonlinear phase diffusion equation or the perturbed linear eikonal equation. While the control enforces a desired wave shape perpendicular to the local propagation direction, the wave profile along the propagation direction itself remains almost unaffected. Provided that the one-dimensional wave profile and its propagation velocity can be measured experimentally, and the diffusion coefficients of the reacting species are given, the new approach can be applied even if the underlying nonlinear reaction kinetics are unknown., Comment: 12 pages, 4 figures

Using the Schl\"{o}gl model as a paradigmatic example of a bistable reaction-diffusion system, we discuss some physically feasible options of open and closed loop spatio-temporal control of chemical wave propagation., Comment: 22 pages, 7 figures

We consider the stability of position control of traveling waves in reaction-diffusion system as proposed in {[}J. L\"ober, H. Engel, arXiv:1304.2327{]}. Instead of analyzing the controlled reaction-diffusion system, stability is studied on the reduced level of the equation of motion for the position over time of perturbed traveling waves. We find an interval of perturbations of initial conditions for which position control is stable. This interval can be interpreted as a localized region where traveling waves are susceptible to perturbations. For stationary solutions of reaction-diffusion systems with reflection symmetry, this region does not exist. Analytical results are in qualitative agreement with numerical simulations of the controlled Schl\"ogl model.

We present a method to control the position as a function of time of one-dimensional traveling wave solutions to reaction-diffusion systems according to a pre-specified protocol of motion. Given this protocol, the control function is found as the solution of a perturbatively derived integral equation. Two cases are considered. First, we derive an analytical expression for the space ($x$) and time ($t$) dependent control function $f\left(x,t\right)$ that is valid for arbitrary protocols and many reaction-diffusion systems. These results are close to numerically computed optimal controls. Second, for stationary control of traveling waves in one-component systems, the integral equation reduces to a Fredholm integral equation of the first kind. In both cases, the control can be expressed in terms of the uncontrolled wave profile and its propagation velocity, rendering detailed knowledge of the reaction kinetics unnecessary.

We propose a non-perturbative attempt to solve the kinematic equations for spiral waves in excitable media. From the eikonal equation for the wave front we derive an implicit analytical relation between rotation frequency $\Omega$ and core radius $R_{0}$. For free, rigidly rotating spiral waves our analytical prediction is in good agreement with numerical solutions of the linear eikonal equation not only for very large but also for intermediate and small values of the core radius. An equivalent $\Omega\left(R_{+}\right)$ dependence improves the result by Keener and Tyson for spiral waves pinned to a circular defect with radius $R_{+}$ with Neumann boundaries at the periphery. Simultaneously, analytical approximations for the shape of free and pinned spirals are given. We discuss the reasons why the ansatz fails to correctly describe the result for the dependence of the rotation frequency on the excitability of the medium., Comment: 12 pages, 6 figures

Wave propagation in one-dimensional heterogeneous bistable media is studied using the Schl\"ogl model as a representative example. Starting from the analytically known traveling wave solution for the homogeneous medium, infinitely extended, spatially periodic variations in kinetic parameters as the excitation threshold, for example, are taken into account perturbatively. Two different multiple scale perturbation methods are applied to derive a differential equation for the position of the front under perturbations. This equation allows the computation of a time independent average velocity, depending on the spatial period length and the amplitude of the heterogeneities. The projection method reveals to be applicable in the range of intermediate and large period lengths but fails when the spatial period becomes smaller than the front width. Then, a second order averaging method must be applied. These analytical results are capable to predict propagation failure, velocity overshoot, and the asymptotic value for the front velocity in the limit of large period lengths in qualitative, often quantitative agreement with the results of numerical simulations of the underlying reaction-diffusion equation. Very good agreement between numerical and analytical results has been obtained for waves propagating through a medium with periodically varied excitation threshold., Comment: 13 pages, 9 figures