Signalverarbeitung im Nervensystem hängt sowohl von der Netzwerkstruktur, als auch den zellulären Eigenschaften der Nervenzellen ab. In dieser Abhandlung werden zwei zelluläre Eigenschaften im Hinblick auf ihre funktionellen Anpassungsmöglichkeiten untersucht: Es wird gezeigt, dass neuronale Morphologie die Signalweiterleitung unter Berücksichtigung energetischer Beschränkungen verstärken kann, und dass selbst kleine Änderungen in biophysikalischen Parametern die Aktivierungsbifurkation in Nervenzellen, und damit deren Informationskodierung, wechseln können. Im ersten Teil dieser Abhandlung wird, unter Verwendung von mathematischen Modellen und Daten, die Hypothese aufgestellt, dass Energie-effiziente Signalweiterleitung als starker Evolutionsdruck für unterschiedliche Zellkörperlagen bei Nervenzellen wirkt. Um Energie zu sparen, kann die Signalweiterleitung vom Dendrit zum Axon verstärkt werden, indem relativ kleine Zellkörper zwischen Dendrit und Axon eingebaut werden, während relativ große Zellkörper besser ausgelagert werden. Im zweiten Teil wird gezeigt, dass biophysikalische Parameter, wie Temperatur, Membranwiderstand oder Kapazität, den Feuermechanismus des Neurons ändern, und damit gleichfalls Aktionspotential-basierte Informationsverarbeitung. Diese Arbeit identifiziert die sogenannte "saddle-node-loop" (Sattel-Knoten-Schlaufe) Bifurkation als den Übergang, der besonders drastische funktionale Auswirkungen hat. Neben der Änderung neuronaler Filtereigenschaften sowie der Ankopplung an Stimuli, führt die "saddle-node-loop" Bifurkation zu einer Erhöhung der Netzwerk-Synchronisation, was möglicherweise für das Auslösen von Anfällen durch Temperatur, wie bei Fieberkrämpfen, interessant sein könnte. / Signal processing in nervous systems is shaped by the connectome as well as the cellular properties of nerve cells. In this thesis, two cellular properties are investigated with respect to the functional adaptations they provide: It is shown that neuronal morphology can improve signal transmission under energetic constraints, and that even small changes in biophysical parameters can switch spike generation, and thus information encoding. In the first project of the thesis, mathematical modeling and data are deployed to suggest energy-efficient signaling as a major evolutionary pressure behind morphological adaptations of cell body location: In order to save energy, the electrical signal transmission from dendrite to axon can be enhanced if a relatively small cell body is located between dendrite and axon, while a relatively large cell body should be externalized. In the second project, it is shown that biophysical parameters, such as temperature, membrane leak or capacitance, can transform neuronal excitability (i.e., the spike onset bifurcation) and, with that, spike-based information processing. This thesis identifies the so-called saddle-node-loop bifurcation as the transition with particularly drastic functional implications. Besides altering neuronal filters and stimulus locking, the saddle-node-loop bifurcation leads to an increase in network synchronization, which may potentially be relevant for the initiation of seizures in response to increased temperature, such as during fever cramps.

It is now widely assumed that many artefacts found in the prehistoric archaeological record were not casually discarded as unwanted material but were deposited in features and contexts with structure and meaning. This appears to include saddle and rotary querns for they are often found whole and apparently still usable or, conversely, deliberately broken. Analysis of the structured deposition of querns in the south-west of England shows that they were deposited in features on both domestic and non-domestic sites. Furthermore, the location and state of the querns, together with the artefacts found in association with them, indicates that they were deposited with different levels and layers of meaning, even within the same type of feature. The deposition of querns appears to have pervaded all aspects of prehistoric life and death suggesting that they played a role above, but nevertheless related to, their prime task of milling. An exploration of the object biography of querns demonstrates the importance of what are often considered to be mundane tools to subsistence communities. Each quern has its own unique life history, its meaning and value determined by the reasons that gave cause for its manufacture, the material from which it was made, the use(s) to which it was put and who used it. However, all querns share points of commonality, related to their function as milling tools, their role as transformers of raw material(s) into usable products (s), their association with women and the production of food, and the movement of the upper stone. Through these, symbolical links can be made between querns and agricultural, human and building life cycles, gender relations and the turning of the heavens. The reason for a quern’s deposition in the archaeological record may have drawn upon one or more unique or common values.

Cette thèse vise à concevoir des solveurs efficaces pour résoudre des systèmes linéaires, résultant des problèmes d'optimisation sous contraintes dans certaines applications de dynamique des structures et vibration (la corrélation calcul-essai, la localisation d'erreur, le modèle hybride, l'évaluation des dommages, etc.). Ces applications reposent sur la résolution de problèmes inverses, exprimés sous la forme de la minimisation d'une fonctionnelle en énergie. Cette fonctionnelle implique à la fois, des données issues d'un modèle numérique éléments finis, et des essais expérimentaux. Ceci conduit à des modèles de haute qualité, mais les systèmes linéaires point-selle associés, sont coûteux à résoudre. Nous proposons deux classes différentes de méthodes pour traiter le système. La première classe repose sur une méthode de factorisation directe profitant de la topologie et des propriétés spéciales de la matrice point-selle. Après une première renumérotation pour regrouper les pivots en blocs d'ordre 2. L'élimination de Gauss est conduite à partir de ces pivots et en utilisant un ordre spécial d'élimination réduisant le remplissage. Les résultats numériques confirment des gains significatifs en terme de remplissage, jusqu'à deux fois meilleurs que la littérature pour la topologie étudiée. La seconde classe de solveurs propose une approche à double projection du système étudié sur le noyau des contraintes, en faisant une distinction entre les contraintes cinématiques et celles reliées aux capteurs sur la structure. La première projection est explicite en utilisant une base creuse du noyau. La deuxième est implicite. Elle est basée sur l'emploi d'un préconditionneur contraint avec des méthodes itératives de type Krylov. Différentes approximations des blocs du préconditionneur sont proposées. L'approche est implémentée dans un environnement distribué parallèle utilisant la bibliothèque PETSc. Des gains significatifs en terme de coût de calcul et de mémoire sont illustrés sur plusieurs applications industrielles. / This thesis aims at designing efficient numerical solution methods to solve linear systems, arising in constrained optimization problems in some structural dynamics and vibration applications (test-analysis correlation, model error localization,hybrid model, damage assessment, etc.). These applications rely on solving inverse problems, by means of minimization of an energy-based functional. This latter involves both data from a numerical finite element model and from experimental tests, which leads to high quality models, but the associated linear systems, that have a saddle-point coefficient matrices, are long and costly to solve. We propose two different classes of methods to deal with these problems. First, a direct factorization method that takes advantage of the special structures and properties of these saddle point matrices. The Gaussian elimination factorization is implemented in order to factorize the saddle point matrices block-wise with small blocks of orders 2 and using a fill-in reducing topological ordering. We obtain significant gains in memory cost (up to 50%) due to enhanced factors sparsity in comparison to literature. The second class is based on a double projection of the generated saddle point system onto the nullspace of the constraints. The first projection onto the kinematic constraints is proposed as an explicit process through the computation of a sparse null basis. Then, we detail the application of a constraint preconditioner within a Krylov subspace solver, as an implicit second projection of the system onto the nullspace of the sensors constraints. We further present and compare different approximations of the constraint preconditioner. The approach is implemented in a parallel distributed environment using the PETSc library. Significant gains in computational cost and memory are illustrated on several industrial applications.

This thesis presents numerical techniques for solving problems of incompressible flow coupled to scalar transport equations using finite element discretizations in space. The two applications considered in this thesis are multi-phase flow, modeled by level set or phase field methods, and planetary mantle convection based on the Boussinesq approximation. A systematic numerical study of approximation errors in evaluating the surface tension in finite element models for two-phase flow is presented. Forces constructed from a gradient in the same discrete function space as used for the pressure are shown to give the best performance. Moreover, two approaches for introducing contact line dynamics into level set methods are proposed. Firstly, a multiscale approach extracts a slip velocity from a micro simulation based on the phase field method and imposes it as a boundary condition in the macro model. This multiscale method is shown to provide an efficient model for the simulation of contact-line driven flow. The second approach combines a level set method based on a smoothed color function with a the phase field method in different parts of the domain. Away from contact lines, the additional information in phase field models is not necessary and it is disabled from the equations by a switch function. An in-depth convergence study is performed in order to quantify the benefits from this combination. Also, the resulting hybrid method is shown to satisfy an a priori energy estimate. For the simulation of mantle convection, an implementation framework based on modern finite element and solver packages is presented. The framework is capable of running on today's large computing clusters with thousands of processors. All parts in the solution chain, from mesh adaptation over assembly to the solution of linear systems, are done in a fully distributed way. These tools are used for a parallel solver that combines higher order time and space discretizations. For treating the convection-dominated temperature equation, an advanced stabilization technique based on an artificial viscosity is used. For more efficient evaluation of finite element operators in iterative methods, a matrix-free implementation built on cell-based quadrature is proposed. We obtain remarkable speedups over sparse matrix-vector products for many finite elements which are of practical interest. Our approach is particularly efficient for systems of differential equations.

This thesis deals with the design of a robust and safe control algorithm to aim at an artificial pancreas. More precisely we will be interested in controlling the stabilizing part of a classical cure. To meet this objective, the design of a robust nonlinear model predictive controller based on the solution of a saddle point optimization problem is considered. Also, to test the controller performances in a realistic case, numerical simulations on a FDA validated testing platform are envisaged.In a first part, we present an extension of the usual nonlinear model predictive controller designed to robustly control, in a sampled-data framework, systems described by nonlinear ordinary differential equations. This controller, which computes the best control input by considering the solution of a constrained saddle point optimization problem, is called saddle point model predictive controller (SPMPC). Using this controller, it is proved that the closed-loop is Ultimately Bounded and, with some assumptions on the problem structure, Input-to State practically Stable. Then, we are interested in numerically solving the corresponding control problem. To do so, we propose an algorithm inspired from the augmented Lagrangian technique and which makes use of adjoint model.In a second part, we consider the application of this controller to the problem of artificial blood glucose control. After a modeling phase, two models are retained. A simple one will be used to design the controller and a complex one will be used to simulate realistic virtual patients. This latter is needed to validate our control approach. In order to compute a good control input, the SPMPC controller needs the full state value. However, the sensors can only provide the value of blood glucose. That is why the design of an adequate observer is envisaged. Then, numerical simulations are performed. The results show the interest of the approach. For all virtual patients, no hypoglycemia event occurs and the time spent in hyperglycemia is too short to induce damageable consequences. Finally, the interest of extending the SPMPC approach to consider the control of time delay systems in a sampled-data framework is numerically explored.

One of the most important design issues for phase change materials is to engineer the phase transition process. The challenge of accurately predicting a phase transition is estimating the true value of transition rate, which is determined by the saddle point with the minimum energy barrier between stable states on the potential energy surface (PES). In this thesis, a new algorithm for searching the minimum energy path (MEP) is presented. The new algorithm is able to locate both the saddle point and local minima simultaneously. Therefore no prior knowledge of the precise positions for the reactant and product on the PES is needed. Unlike existing pathway search methods, the algorithm is able to search multiple transition paths on the PES simultaneously, which gives us a more comprehensive view of the energy landscape than searching individual ones. In this method, a Bézier curve is used to represent each transition path. During the searching process, the reactant and product states are located by minimizing the two end control points of the curve, while the shape of the transition pathway is refined by moving the intermediate control points of the curve in the conjugate directions. A curve subdivision scheme is developed so that multiple transitions paths can be located. The algorithm is demonstrated by examples of LEPS potential, LEPS plus harmonic oscillator potential, and PESs defined by Rastrigin function and Schwefel function.

Plasma enhanced chemical vapor deposition systems are massively deployed to grow numerous thin film coatings including hydrogenated amorphous silicon. A new deposition chamber was designed, procured, and constructed to investigate the plasma properties of a 100% silane (SiH4) glow discharge with varying chamber pressure and inter-electrode spacing. A Hiden EQP1000 ion mass spectrometer sampled the plasma from the substrates point of view. Ion energy distributions were obtained using four different excitation sources +DC, –DC, radio frequency (at 13.56 MHz), and the DC Saddle Field (DCSF) in the tetrode configuration. The shape of the ion energy distributions was constant for the capacitively coupled +DC, –DC, and rf (at higher pressures of 75 and 160 mTorr) glow discharges. The shape of the ion energy distributions for the DCSF plasma exhibited a double peak or saddle structure analogous to radio frequency plasmas. The width between the peaks (peak separation) was controlled by the pressure and the semi-transparent cathode to semi-transparent anode distance. Ion energy distributions from the DCSF plasma concurred with rf and +DC ion energy distributions at specific pressures and inter-electrode distances. This result demonstrates the versatility of the DCSF glow discharge system. Moreover, control of the peak separation is modeled to be iii equivalent to controlling the critical ratio (ion transit time in the sheath to the electron oscillating period), and/or the inferred electron oscillating sheath potential. The DCSF possesses a fusion of rf and +DC methods. The long high energy tail or constant background are indicative of a +DC high voltage sheath in which there is an increasing fraction of collisionless ions as the anode-cathode distance increases. These collisionless ions are provided by the oscillating electrons (or rf nature) of the DCSF method. Higher order silane (silicon containing) ions increase in relative intensity with increasing inter-electrode spacing for the +DC, –DC, and rf plasmas. These higher order silane ions are also detected in the DCSF plasma, and can be reduced at either lower pressure or lower cathode to anode or cathode to substrate distances.

Plasma enhanced chemical vapor deposition systems are massively deployed to grow numerous thin film coatings including hydrogenated amorphous silicon. A new deposition chamber was designed, procured, and constructed to investigate the plasma properties of a 100% silane (SiH4) glow discharge with varying chamber pressure and inter-electrode spacing. A Hiden EQP1000 ion mass spectrometer sampled the plasma from the substrates point of view. Ion energy distributions were obtained using four different excitation sources +DC, –DC, radio frequency (at 13.56 MHz), and the DC Saddle Field (DCSF) in the tetrode configuration. The shape of the ion energy distributions was constant for the capacitively coupled +DC, –DC, and rf (at higher pressures of 75 and 160 mTorr) glow discharges. The shape of the ion energy distributions for the DCSF plasma exhibited a double peak or saddle structure analogous to radio frequency plasmas. The width between the peaks (peak separation) was controlled by the pressure and the semi-transparent cathode to semi-transparent anode distance. Ion energy distributions from the DCSF plasma concurred with rf and +DC ion energy distributions at specific pressures and inter-electrode distances. This result demonstrates the versatility of the DCSF glow discharge system. Moreover, control of the peak separation is modeled to be iii equivalent to controlling the critical ratio (ion transit time in the sheath to the electron oscillating period), and/or the inferred electron oscillating sheath potential. The DCSF possesses a fusion of rf and +DC methods. The long high energy tail or constant background are indicative of a +DC high voltage sheath in which there is an increasing fraction of collisionless ions as the anode-cathode distance increases. These collisionless ions are provided by the oscillating electrons (or rf nature) of the DCSF method. Higher order silane (silicon containing) ions increase in relative intensity with increasing inter-electrode spacing for the +DC, –DC, and rf plasmas. These higher order silane ions are also detected in the DCSF plasma, and can be reduced at either lower pressure or lower cathode to anode or cathode to substrate distances.

Diss. (sammanfattning) Uppsala : Sveriges lantbruksuniv., 2002. / Härtill 6 uppsatser.

Orientador: Vanderlei Minori Horita / Banca: Ali Tahzibi / Banca: Paulo Ricardo Silva / Resumo: Estudamos intersecções homoclínicas de variedades estável e instável de pontos peródicos. Toda intersecção homoclínica produz um comportamento curioso na dinâmiôa. Nosso modelo de tal fenômeno é a famosa ferradura de Smale, a qual é um conjunto hiperbólico para um difeomorfismo. Além disso, estudamos dinâmica não hiperbólica cuja perda de hiperbolicidade é divido à tangências homoclínicas. Elas tem um papel central na teoria de sistemas dinâmicos. O desdobramento de uma tangência homoclínica produz dinâmicas muito interessantes. Neste trabalho estudamos a criação de cascatas de bifurcações de duplicação de período e um esquema de renormalização para uma tangência homoclínica. / Abstract: We study homoclinic intersection of stable and unstable manifolds of periodic points. Every homoclinic intersection produce a intricate behavior of the dynamics. Our model of such phenomena is the so called Smalesþs horseshoe, which is a hyperbolic set for a di eomorphism. We also study non hyperbolic dynamics whose lack of hyperbolicity is due to homoclinic tangencies. They play a central role in the theory of dynamical systems. The unfolding of a homoclinic tangency produce many interesting dynamics. In this work we study creation of cascade of period doubling bifurcations and a renormalization scheme for a homoclinic tangency. / Mestre