Coupled effects of mechanics, geometry, and chemistry on bio-membrane behavior

Lipid bilayer membranes are models for cell membranes--the structure that helps regulate cell function. Cell membranes are heterogeneous, and the coupling between composition and shape gives rise to complex behaviors that are important to regulation. This thesis seeks to systematically build and analyze complete models to understand the behavior of multi-component membranes.

We propose a model and use it to derive the equilibrium and stability conditions for a general class of closed multi-component biological membranes. Our analysis shows that the critical modes of these membranes have high frequencies, unlike single-component vesicles, and their stability depends on system size, unlike in systems undergoing spinodal decomposition in flat space. An important implication is that small perturbations may nucleate localized but very large deformations. We compare these results with experimental observations.

We also study open membranes to gain insight into long tubular membranes that arise for example in nerve cells. We derive a complete system of equations for open membranes by using the principle of virtual work. Our linear stability analysis predicts that the tubular membranes tend to have coiling shapes if the tension is small, cylindrical shapes if the tension is moderate, and beading shapes if the tension is large. This is consistent with experimental observations reported in the literature in nerve fibers. Further, we provide numerical solutions to the fully nonlinear equilibrium equations in some problems, and show that the observed mode shapes are consistent with those suggested by linear stability. Our work also proves that beadings of nerve fibers can appear purely as a mechanical response of the membrane.

Response of linear, viscous damped systems to excitations having time-varying frequency

The response of linear, viscous damped systems to excitations having time-varying frequency is the subject of exact and approximate analyses, which are supplemented by an analog computer study of single degree of freedom system response to excitations having frequencies depending linearly and exponentially on time.

The technique of small perturbations and the methods of stationary phase and saddle-point integration, as well as a novel bounding procedure, are utilized to derive approximate expressions characterizing the system response envelope—particularly near resonances—for the general time-varying excitation frequency.

Descriptive measurements of system resonant behavior recorded during the course of the analog study—maximum response, excitation frequency at which maximum response occurs, and the width of the response peak at the half-power level—are investigated to determine dependence upon natural frequency, damping, and the functional form of the excitation frequency.

The laboratory problem of determining the properties of a physical system from records of its response to excitations of this class is considered, and the transient phenomenon known as “ringing” is treated briefly.

It is shown that system resonant behavior, as portrayed by the above measurements and expressions, is relatively insensitive to the specifics of the excitation frequency-time relation and may be described to good order in terms of parameters combining system properties with the time derivative of excitation frequency evaluated at resonance.

One of these parameters is shown useful for predicting whether or not a given excitation having a time-varying frequency will produce strong or subtle changes in the response envelope of a given system relative to the steady-state response envelope. The parameter is shown, additionally, to be useful for predicting whether or not a particular response record will exhibit the “ringing” phenomenon.

Relativistic velocity - potential hydrodynamics and stellar stability

The equations of relativistic, perfect-fluid hydrodynamics are cast in Eulerian form using six scalar "velocity-potential" fields, each of which has an equation of evolution. These equations determine the motion of the fluid through the equation

U_ʋ=µ^-1 (ø,_ʋ + αβ,_ʋ + ƟS,_ʋ).

Einstein's equations and the velocity-potential hydrodynamical equations follow from a variational principle whose action is

I = (R + 16π p) (-g)^1/2 d⁴x,

where R is the scalar curvature of spacetime and p is the pressure of the fluid. These equations are also cast into Hamiltonian form, with Hamiltonian density –T₀⁰ (-g^oo)^-1/2.

The second variation of the action is used as the Lagrangian governing the evolution of small perturbations of differentially rotating stellar models. In Newtonian gravity this leads to linear dynamical stability criteria already known. In general relativity it leads to a new sufficient condition for the stability of such models against arbitrary perturbations.

By introducing three scalar fields defined by

ρ ᵴ = ∇λ + ∇x(xi + ∇xɣi)

(where ᵴ is the vector displacement of the perturbed fluid element, ρ is the mass-density, and i, is an arbitrary vector), the Newtonian stability criteria are greatly simplified for the purpose of practical applications. The relativistic stability criterion is not yet in a form that permits practical calculations, but ways to place it in such a form are discussed.

Application of the fixed point method for solution in time stepping finite element analysis using the inverse vector Jiles-Atherton model

An implementation of the inverse vector Jiles-Atherton model for the solution of non-linear hysteretic finite element problems is presented. The implementation applies the fixed point method with differential reluctivity values obtained from the Jiles-Atherton model. Differential reluctivities are usually computed using numerical differentiation, which is ill-posed and amplifies small perturbations causing large sudden increases or decreases of differential reluctivity values, which may cause numerical problems. A rule based algorithm for conditioning differential reluctivity values is presented. Unwanted perturbations on the computed differential reluctivity values are eliminated or reduced with the aim to guarantee convergence. Details of the algorithm are presented together with an evaluation of the algorithm by a numerical example. The algorithm is shown to guarantee convergence, although the rate of convergence depends on the choice of algorithm parameters. © 2011 IEEE.

A finite element model of magnetization of superconducting bulks using a solid-state flux pump

Superconductors have a bright future; they are able to carry very high current densities, switch rapidly in electronic circuits, detect extremely small perturbations in magnetic fields, and sustain very high magnetic fields. Of most interest to large-scale electrical engineering applications are the ability to carry large currents and to provide large magnetic fields. There are many projects that use the first property, and these have concentrated on power generation, transmission, and utilization; however, there are relatively few, which are currently exploiting the ability to sustain high magnetic fields. The main reason for this is that high field wound magnets can and have been made from both BSCCO and YBCO, but currently, their cost is much higher than the alternative provided by low-Tc materials such as Nb3Sn and NbTi. An alternative form of the material is the bulk form, which can be magnetized to high fields. This paper explains the mechanism, which allows superconductors to be magnetized without the need for high field magnets to perform magnetization. A finite-element model is presented, which is based on the E-J current law. Results from this model show how magnetization of the superconductor builds up cycle upon cycle when a traveling magnetic wave is induced above the superconductor. © 2011 IEEE.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

Three primary kinetic characteristics observed in a pulse-response TAP experiment

The state-by-state transient screening approach based on a pulse-response thin-zone TAP experiment is further developed whereby single-pulse kinetic tests are treated as small perturbations to catalyst compositions and analyzed using integral method of moments. Results on three primary kinetic characteristics, termed basic kinetic coefficients, are presented. These three coefficients were introduced as main observables from experimentally measured TAP-responses in a kinetic-model-free manner. Each was analytically determined from moments of responses with no assumption about the detailed kinetic model. In this paper, the inverse question of how well these coefficients represent the time evolution of the observed responses is addressed. Sets of three basic kinetic coefficients are calculated from model and experimental responses and these calculated values are used to generate 3-coefficient curves in a kinetic-model-free manner. The comparison of these 3-coefficient curves with original responses shows that three basic kinetic coefficients can be sufficient to describe the observed kinetics of exit flow time dependencies with no assumption regarding the detailed kinetic model.

Studies on integrability of some perturbed nonlinear evolution equations

Usually typical dynamical systems are non integrable. But few systems of practical interest are integrable. The soliton concept is a sophisticated mathematical construct based on the integrability of a class ol' nonlinear differential equations. An important feature in the clevelopment. of the theory of solitons and of complete integrability has been the interplay between mathematics and physics. Every integrable system has a lo11g list of special properties that hold for integrable equations and only for them. Actually there is no specific definition for integrability that is suitable for all cases. .There exist several integrable partial clillerential equations( pdes) which can be derived using physically meaningful asymptotic teclmiques from a very large class of pdes. It has been established that many 110nlinear wa.ve equations have solutions of the soliton type and the theory of solitons has found applications in many areas of science. Among these, well-known equations are Korteweg de-Vries(KdV), modified KclV, Nonlinear Schr6dinger(NLS), sine Gordon(SG) etc..These are completely integrable systems. Since a small change in the governing nonlinear prle may cause the destruction of the integrability of the system, it is interesting to study the effect of small perturbations in these equations. This is the motivation of the present work.

Análisis de sensibilidad para un modelo de medición de percepción de marca basado en simulación computacional por agentes

Se presenta el análisis de sensibilidad de un modelo de percepción de marca y ajuste de la inversión en marketing desarrollado en el Laboratorio de Simulación de la Universidad del Rosario. Este trabajo de grado consta de una introducción al tema de análisis de sensibilidad y su complementario el análisis de incertidumbre. Se pasa a mostrar ambos análisis usando un ejemplo simple de aplicación del modelo mediante la aplicación exhaustiva y rigurosa de los pasos descritos en la primera parte. Luego se hace una discusión de la problemática de medición de magnitudes que prueba ser el factor más complejo de la aplicación del modelo en el contexto práctico y finalmente se dan conclusiones sobre los resultados de los análisis.

Estimating climatically relevant singular vectors for decadal predictions of the Atlantic Ocean

A key aspect in designing an ecient decadal prediction system is ensuring that the uncertainty in the ocean initial conditions is sampled optimally. Here, we consider one strategy to address this issue by investigating the growth of optimal perturbations in the HadCM3 global climate model (GCM). More specically, climatically relevant singular vectors (CSVs) - the small perturbations which grow most rapidly for a specic initial condition - are estimated for decadal timescales in the Atlantic Ocean. It is found that reliable CSVs can be estimated by running a large ensemble of integrations of the GCM. Amplication of the optimal perturbations occurs for more than 10 years, and possibly up to 40 years. The identi ed regions for growing perturbations are found to be in the far North Atlantic, and these perturbations cause amplication through an anomalous meridional overturning circulation response. Additionally, this type of analysis potentially informs the design of future ocean observing systems by identifying the sensitive regions where small uncertainties in the ocean state can grow maximally. Although these CSVs are expensive to compute, we identify ways in which the process could be made more ecient in the future.

Response theory for equilibrium and non-equilibrium statistical mechanics: causality and generalized Kramers-Kronig relations

We consider the general response theory recently proposed by Ruelle for describing the impact of small perturbations to the non-equilibrium steady states resulting from Axiom A dynamical systems. We show that the causality of the response functions entails the possibility of writing a set of Kramers-Kronig (K-K) relations for the corresponding susceptibilities at all orders of nonlinearity. Nonetheless, only a special class of directly observable susceptibilities obey K-K relations. Specific results are provided for the case of arbitrary order harmonic response, which allows for a very comprehensive K-K analysis and the establishment of sum rules connecting the asymptotic behavior of the harmonic generation susceptibility to the short-time response of the perturbed system. These results set in a more general theoretical framework previous findings obtained for optical systems and simple mechanical models, and shed light on the very general impact of considering the principle of causality for testing self-consistency: the described dispersion relations constitute unavoidable benchmarks that any experimental and model generated dataset must obey. The theory exposed in the present paper is dual to the time-dependent theory of perturbations to equilibrium states and to non-equilibrium steady states, and has in principle similar range of applicability and limitations. In order to connect the equilibrium and the non equilibrium steady state case, we show how to rewrite the classical response theory by Kubo so that response functions formally identical to those proposed by Ruelle, apart from the measure involved in the phase space integration, are obtained. These results, taking into account the chaotic hypothesis by Gallavotti and Cohen, might be relevant in several fields, including climate research. In particular, whereas the fluctuation-dissipation theorem does not work for non-equilibrium systems, because of the non-equivalence between internal and external fluctuations, K-K relations might be robust tools for the definition of a self-consistent theory of climate change.

DIAGNOSTIC TECHNIQUES OF LOCAL INFLUENCE IN SPATIAL ANALYSIS OF SOYBEAN YIELD

Modeling of spatial dependence structure, concerning geoestatistics approach, is an indispensable tool for fixing parameters that define this structure, applied on interpolation of values in places that are not sampled, by kriging techniques. However, the estimation of parameters can be greatly affected by the presence of atypical observations on sampled data. Thus, this trial aimed at using diagnostics techniques of local influence in spatial linear Gaussians models, applied at geoestatistics in order to evaluate sensitivity of maximum likelihood estimators and restrict maximum likelihood to small perturbations in these data. So, studies with simulated and experimental data were performed. Those results, obtained from the study of real data, allowed us to conclude that the presence of atypical values among the sampled data can have a strong influence on thematic maps, changing, therefore, the spatial dependence. The application of diagnostics techniques of local influence should be part of any geoestatistic analysis, ensuring that the information from thematic maps has better quality and can be used with greater security by farmers.

Semi-hyperbolic mappings in Banach spaces

The definition of semi-hyperbolic dynamical systems generated by Lipschitz continuous and not necessarily invertible mappings in Banach spaces is presented in this thesis. Like hyperbolic mappings, they involve a splitting into stable and unstable spaces, but a slight leakage from the strict invariance of the spaces is possible and the unstable subspaces are assumed to be finite dimensional. Bi-shadowing is a combination of the concepts of shadowing and inverse shadowing and is usually used to compare pseudo-trajectories calculated by a computer with the true trajectories. In this thesis, the concept of bi-shadowing in a Banach space is defined and proved for semi-hyperbolic dynamical systems generated by Lipschitz mappings. As an application to the concept of bishadowing, linear delay differential equations are shown to be bi-shadowing with respect to pseudo-trajectories generated by nonlinear small perturbations of the linear delay equation. This shows robustness of solutions of the linear delay equation with respect to small nonlinear perturbations. Complicated dynamical behaviour is often a consequence of the expansivity of a dynamical system. Semi-hyperbolic dynamical systems generated by Lipschitz mappings on a Banach space are shown to be exponentially expansive, and explicit rates of expansion are determined. The result is applied to a nonsmooth noninvertible system generated by delay differential equation. It is shown that semi-hyperbolic mappings are locally φ-contracting, where -0 is the Hausdorff measure of noncompactness, and that a linear operator is semi-hyperbolic if and only if it is φ-contracting and has no spectral values on the unit circle. The definition of φ-bi-shadowing is given and it is shown that semi-hyperbolic mappings in Banach spaces are φ-bi-shadowing with respect to locally condensing continuous comparison mappings. The result is applied to linear delay differential equations of neutral type with nonsmooth perturbations. Finally, it is shown that a small delay perturbation of an ordinary differential equation with a homoclinic trajectory is ‘chaotic’.

Hierarchy and anarchy in quark mass matrices, or can hierarchy tolerate anarchy?

The consequences of adding random perturbations (anarchy) to a baseline hierarchical model of quark masses and mixings are explored. Even small perturbations of the order of 5% of the smallest non-zero element can already give deviations significantly affecting parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, so any process generating the anarchy should in general be limited to this order of magnitude. The regularities of quark masses and mixings thus appear to be far from a generic feature of randomness in the mass matrices, and more likely indicate an underlying order. (C) 2001 Published by Elsevier B.V. B.V.