Global analysis of dynamical decision-making models through local computation around the hidden saddle

Bistable dynamical switches are frequently encountered in mathematical modeling of biological systems because binary decisions are at the core of many cellular processes. Bistable switches present two stable steady-states, each of them corresponding to a distinct decision. In response to a transient signal, the system can flip back and forth between these two stable steady-states, switching between both decisions. Understanding which parameters and states affect this switch between stable states may shed light on the mechanisms underlying the decision-making process. Yet, answering such a question involves analyzing the global dynamical (i.e., transient) behavior of a nonlinear, possibly high dimensional model. In this paper, we show how a local analysis at a particular equilibrium point of bistable systems is highly relevant to understand the global properties of the switching system. The local analysis is performed at the saddle point, an often disregarded equilibrium point of bistable models but which is shown to be a key ruler of the decision-making process. Results are illustrated on three previously published models of biological switches: two models of apoptosis, the programmed cell death and one model of long-term potentiation, a phenomenon underlying synaptic plasticity. © 2012 Trotta et al.

Global modes with multiple saddle points

Significant progress has been made towards understanding the global stability of slowly-developing shear flows. The WKBJ theory developed by Patrick Huerre and his co-authors has proved absolutely central, with the result that both the linear and the nonlinear stability of a wide range of flows can now be understood in terms of their local absolute/convective instability properties. In many situations, the local absolute frequency possesses a single dominant saddle point in complex X-space (where X is the slow streamwise coordinate of the base flow), which then acts as a single wavemaker driving the entire global linear dynamics. In this paper we consider the more complicated case in which multiple saddles may act as the wavemaker for different values of some control parameter. We derive a frequency selection criterion in the general case, which is then validated against numerical results for the linearized third-order Ginzburg-Landau equation (which possesses two saddle points). We believe that this theory may be relevant to a number of flows, including the boundary layer on a rotating disk and the eccentric Taylor-Couette-Poiseuille flow. © 2014 Elsevier Masson SAS. All rights reserved.

Distributed flow algorithms for scalable similarity visualization

We describe simple yet scalable and distributed algorithms for solving the maximum flow problem and its minimum cost flow variant, motivated by problems of interest in objects similarity visualization. We formulate the fundamental problem as a convex-concave saddle point problem. We then show that this problem can be efficiently solved by a first order method or by exploiting faster quasi-Newton steps. Our proposed approach costs at most O(|ε|) per iteration for a graph with |ε| edges. Further, the number of required iterations can be shown to be independent of number of edges for the first order approximation method. We present experimental results in two applications: mosaic generation and color similarity based image layouting. © 2010 IEEE.

Delayed decision-making in bistable models

Switching between two modes of operation is a common property of biological systems. In continuous-time differential equation models, this is often realised by bistability, i.e. the existence of two asymptotically stable steadystates. Several biological models are shown to exhibit delayed switching, with a pronounced transient phase, in particular for near-threshold perturbations. This study shows that this delay in switching from one mode to the other in response to a transient input is reflected in local properties of an unstable saddle point, which has a one dimensional unstable manifold with a significantly slower eigenvalue than the stable ones. Thus, the trajectories first approximatively converge to the saddle point, then linger along the saddle's unstable manifold before quickly approaching one of the stable equilibria. ©2010 IEEE.

Delayed-decision making in noisy bistable switches

Bistable switches are frequently encountered in biological systems. Typically, a bistable switch models a binary decision where each decision corresponds to one of the two stable equilibria. Recently, we showed that the global decision-making process in bistable switches strongly depends on a particular equilibrium point of these systems, their saddle point. In particular, we showed that a saddle point with a time-scale separation between its attractive and repulsive directions can delay the decision-making process. In this paper, we study the effects of white Gaussian noise on this mechanism of delayed decision-making induced by the saddle point. Results show that the mean decision-time strongly depends on the balance between the initial distance to the separatrix and the noise strength. © IFAC.

Transient Ventilation Dynamics Following a Change in Strength of a Point Source of Heat

We investigate the transient ventilation flow within a confined ventilated space, with high- and low-level openings, when the strength of a low-level point source of heat is changed instantaneously. The steady-flow regime in the space involves a turbulent buoyant plume, which rises from the point source to a well-mixed warm upper layer. The steady-state height of the interface between this layer and the lower layer of exterior fluid is independent of the heat flux, but the upper layer becomes progressively warmer with heat flux. New analogue laboratory experiments of the transient adjustment between steady states identify that if the heat flux is increased, the continuing plume propagates to the top of the room forming a new, warmer layer. This layer gradually deepens, and as the turbulent plume entrains fluid from the original warm layer, the original layer is gradually depleted and disappears, and a new steady state is established. In contrast, if the source buoyancy flux is decreased, the continuing plume is cooler than the original plume, so that on reaching the interface it is of intermediate density between the original warm layer and the external fluid. The plume supplies a new intermediate layer, which gradually deepens with the continuing flow. In turn, the original upper layer becomes depleted, both as a result of being vented through the upper opening of the space, but also due to some penetrative entrainment of this layer by the plume, as the plume overshoots the interface before falling back to supply the new intermediate layer. We develop quantitative models which are in good accord with our experimental data, by combining classical plume theory with models of the penetrative entrainment for the case of a decrease in heating. Typically, we find that the effect of penetrative entrainment on the density of the intruding layer is relatively weak, provided the change in source strength is sufficiently large. However, penetrative entrainment measurably increases the rate at which the depth of the draining layer decreases. We conclude with a discussion of the importance of these results for the control of naturally ventilated spaces.

Computer generated hologram from point cloud using graphics processor

Computer generated holography is an extremely demanding and complex task when it comes to providing realistic reconstructions with full parallax, occlusion, and shadowing. We present an algorithm designed for data-parallel computing on modern graphics processing units to alleviate the computational burden. We apply Gaussian interpolation to create a continuous surface representation from discrete input object points. The algorithm maintains a potential occluder list for each individual hologram plane sample to keep the number of visibility tests to a minimum.We experimented with two approximations that simplify and accelerate occlusion computation. It is observed that letting several neighboring hologramplane samples share visibility information on object points leads to significantly faster computation without causing noticeable artifacts in the reconstructed images. Computing a reduced sample set via nonuniform sampling is also found to be an effective acceleration technique. © 2009 Optical Society of America.

Quasi-static three-point bending of carbon fiber sandwich beams with square honeycomb cores

Sandwich beams comprising identical face sheets and a square honeycomb core were manufactured from carbon fiber composite sheets. Analytical expressions were derived for four competing collapse mechanisms of simply supported and clamped sandwich beams in three-point bending: core shear, face microbuckling, face wrinkling, and indentation. Selected geometries of sandwich beams were tested to illustrate these collapse modes, with good agreement between analytic predictions and measurements of the failure load. Finite element (FE) simulations of the three-point bending responses of these beams were also conducted by constructing a FE model by laying up unidirectional plies in appropriate orientations. The initiation and growth of damage in the laminates were included in the FE calculations. With this embellishment, the FE model was able to predict the measured load versus displacement response and the failure sequence in each of the composite beams. © 2011 American Society of Mechanical Engineers.