Adaptive Mixture Modelling Metropolis Methods for Bayesian Analysis of Non-linear State-Space Models.

We describe a strategy for Markov chain Monte Carlo analysis of non-linear, non-Gaussian state-space models involving batch analysis for inference on dynamic, latent state variables and fixed model parameters. The key innovation is a Metropolis-Hastings method for the time series of state variables based on sequential approximation of filtering and smoothing densities using normal mixtures. These mixtures are propagated through the non-linearities using an accurate, local mixture approximation method, and we use a regenerating procedure to deal with potential degeneracy of mixture components. This provides accurate, direct approximations to sequential filtering and retrospective smoothing distributions, and hence a useful construction of global Metropolis proposal distributions for simulation of posteriors for the set of states. This analysis is embedded within a Gibbs sampler to include uncertain fixed parameters. We give an example motivated by an application in systems biology. Supplemental materials provide an example based on a stochastic volatility model as well as MATLAB code.

Nonholonomically Constrained Dynamics and Optimization of Rolling Isolation Systems

Rolling Isolation Systems provide a simple and effective means for protecting components from horizontal floor vibrations. In these systems a platform rolls on four steel balls which, in turn, rest within shallow bowls. The trajectories of the balls is uniquely determined by the horizontal and rotational velocity components of the rolling platform, and thus provides nonholonomic constraints. In general, the bowls are not parabolic, so the potential energy function of this system is not quadratic. This thesis presents the application of Gauss's Principle of Least Constraint to the modeling of rolling isolation platforms. The equations of motion are described in terms of a redundant set of constrained coordinates. Coordinate accelerations are uniquely determined at any point in time via Gauss's Principle by solving a linearly constrained quadratic minimization. In the absence of any modeled damping, the equations of motion conserve energy. This mathematical model is then used to find the bowl profile that minimizes response acceleration subject to displacement constraint.

The synchronization of superparamagnetic beads driven by a micro-magnetic ratchet.

We present theoretical, numerical, and experimental analyses on the non-linear dynamic behavior of superparamagnetic beads exposed to a periodic array of micro-magnets and an external rotating field. The agreement between theoretical and experimental results revealed that non-linear magnetic forcing dynamics are responsible for transitions between phase-locked orbits, sub-harmonic orbits, and closed orbits, representing different mobility regimes of colloidal beads. These results suggest that the non-linear behavior can be exploited to construct a novel colloidal separation device that can achieve effectively infinite separation resolution for different types of beads, by exploiting minor differences in their bead's properties. We also identify a unique set of initial conditions, which we denote the "devil's gate" which can be used to expeditiously identify the full range of mobility for a given bead type.

Discharge competence and pattern formation in peatlands: a meta-ecosystem model of the Everglades ridge-slough landscape.

Regular landscape patterning arises from spatially-dependent feedbacks, and can undergo catastrophic loss in response to changing landscape drivers. The central Everglades (Florida, USA) historically exhibited regular, linear, flow-parallel orientation of high-elevation sawgrass ridges and low-elevation sloughs that has degraded due to hydrologic modification. In this study, we use a meta-ecosystem approach to model a mechanism for the establishment, persistence, and loss of this landscape. The discharge competence (or self-organizing canal) hypothesis assumes non-linear relationships between peat accretion and water depth, and describes flow-dependent feedbacks of microtopography on water depth. Closed-form model solutions demonstrate that 1) this mechanism can produce spontaneous divergence of local elevation; 2) divergent and homogenous states can exhibit global bi-stability; and 3) feedbacks that produce divergence act anisotropically. Thus, discharge competence and non-linear peat accretion dynamics may explain the establishment, persistence, and loss of landscape pattern, even in the absence of other spatial feedbacks. Our model provides specific, testable predictions that may allow discrimination between the self-organizing canal hypotheses and competing explanations. The potential for global bi-stability suggested by our model suggests that hydrologic restoration may not re-initiate spontaneous pattern establishment, particularly where distinct soil elevation modes have been lost. As a result, we recommend that management efforts should prioritize maintenance of historic hydroperiods in areas of conserved pattern over restoration of hydrologic regimes in degraded regions. This study illustrates the value of simple meta-ecosystem models for investigation of spatial processes.

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

Dynamical Heterogeneity and Nonlinear Susceptibility in Short-Ranged Attractive Supercooled Liquids

Recent work has demonstrated the strong qualitative differences between the dynamics near a glass transition driven by short-ranged repulsion and one governed by short-ranged attraction. Here, we study in detail the behavior of non-linear, higher-order correlation functions that measure the growth of length scales associated with dynamical heterogeneity in both types of systems. We find that this measure is qualitatively different in the repulsive and attractive cases with regards to the wave vector dependence as well as the time dependence of the standard non-linear four-point dynamical susceptibility. We discuss the implications of these results for the general understanding of dynamical heterogeneity in glass-forming liquids.