Data-Driven Dynamic Prediction of Interrelated Heat and Electric Outputs of Microturbines

Microturbines are among the most successfully commercialized distributed energy resources, especially when they are used for combined heat and power generation. However, the interrelated thermal and electrical system dynamic behaviors have not been fully investigated. This is technically challenging due to the complex thermo-fluid-mechanical energy conversion processes which introduce multiple time-scale dynamics and strong nonlinearity into the analysis. To tackle this problem, this paper proposes a simplified model which can predict the coupled thermal and electric output dynamics of microturbines. Considering the time-scale difference of various dynamic processes occuring within microturbines, the electromechanical subsystem is treated as a fast quasi-linear process while the thermo-mechanical subsystem is treated as a slow process with high nonlinearity. A three-stage subspace identification method is utilized to capture the dominant dynamics and predict the electric power output. For the thermo-mechanical process, a radial basis function model trained by the particle swarm optimization method is employed to handle the strong nonlinear characteristics. Experimental tests on a Capstone C30 microturbine show that the proposed modeling method can well capture the system dynamics and produce a good prediction of the coupled thermal and electric outputs in various operating modes.

Reactions of atomic and molecular ions with acetone, 1,1,1- trifluoroacetone, and hexafluoroacetone: An investigation of the effects of molecular structure on the dynamics and kinetics of ion-molecule reactions

A selected ion flow tube study of the reactions of a series of gas-phase atomic cations (S⁺, Xe⁺, O⁺, Kr⁺, N⁺, Ar⁺ and Ne⁺) and molecular ions (SF _n⁺ (n = 1-5), CF_n⁺ (n = 1-3), CF₂Cl⁺, H₃O⁺, NO⁺, N ₂O⁺, CO₂⁺, CO⁺, and N₂⁺) spanning a large range of recombination energies (6.3-21.6 eV), with acetone, 1,1,1-trifluoroacetone, and hexafluoroacetone has been undertaken with the objective of exploring the nature of the reaction ion chemistry as the methyl groups in acetone are substituted for CF₃. The reaction rate coefficients and product ion branching ratios for all 66 reactions, measured at 298 K, are reported. The experimental reaction rate coefficients are compared to theoretically calculated collisional values. Several distinct reaction processes were observed among the large number of reactions studied, including charge transfer (non-dissociative and dissociative), abstraction, ion-molecule associations and, in the case of the reactions involving the reagent ion H₃O⁺, proton transfer. 

Simulation of Distributed Contact in String Instruments: a Modal Expansion Approach

Impactive contact between a vibrating string and a barrier is a strongly nonlinear phenomenon that presents several challenges in the design of numerical models for simulation and sound synthesis of musical string instruments. These are addressed here by applying Hamiltonian methods to incorporate distributed contact forces into a modal framework for discrete-time simulation of the dynamics of a stiff, damped string. The resulting algorithms have spectral accuracy, are unconditionally stable, and require solving a multivariate nonlinear equation that is guaranteed to have a unique solution. Exemplifying results are presented and discussed in terms of accuracy, convergence, and spurious high-frequency oscillations.

Steps towards fully nonlinear simulations of arrays of OWSC

Abandonment of farming systems on upland areas in southwest Britain during the Late Bronze Age – some 3000 years ago – is widely considered a ‘classic’ demonstration of the impact of deteriorating climate on the vulnerability of populations in such marginal environments. Here we test the hypothesis that climate change drove the abandonment of upland areas by developing new chronologies for human activity on upland areas during the Bronze Age across southwest Britain (Dartmoor, Exmoor and Bodmin Moor). We find Bronze Age activity in these areas spanned 3900–2950 calendar years ago with abandonment by 2900 calendar years ago. Holocene Irish bog and lake oak tree populations provide evidence of major shifts in hydroclimate across western Britain and Ireland, coincident with ice rafted debris layers recognized in North Atlantic marine sediments, indicating significant changes in the latitude and intensity of zonal atmospheric circulation across the region. We observe abandonment of upland areas in southwest Britain coinciding with a sustained period of extreme wet conditions that commenced 3100 calendar years ago. Our results are consistent with the view that climate change increased the vulnerability of these early farming communities and led to a less intensive use of such marginal environments across Britain.

The Democratic Dynamics of Government Consultations. Speaking Freely and Listening Properly

There is a growing use of consultation and e-consultation procedures by governments. This chapter seeks to examine the role of consultation as part of a new technology of government. Consultation on policy development can reinvigorate democratic engage- ment but often it can silence views through a sort of participatory disempowerment; it can loosen the democratic anchorage of the public service within the state. The chapter develops a governmentality perspective interrogating what participation, democratic engagement and free speech mean in this context, and how ideas of publicness are constructed, managed and controlled. The focus is on the nature of consultation, its relationship to ideas of free speech and speaking freely, and its potential to empower subaltern counterpublics which can formulate oppositional interpretations and urge alternative conclusions. The aim is to develop an idea of the democratic adequacy of the consultation process and draw out a sense of how democratic engagement here can be structured – for good or ill.

Model reduction for nonlinear aeroelasticity using the Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of computational fluid dynamic (CFD) code for fast nonlinear aerodynamic modeling. First, a nonlinear function is identified on selected interpolation points defined by discrete empirical interpolation method (DEIM). The flow field is then reconstructed by a least square approximation of flow modes extracted by proper orthogonal decomposition (POD). The proposed model is applied in the prediction of limit cycle oscillation for a plunge/pitch airfoil and a delta wing with linear structural model, results are validate against a time accurate CFD-FEM code. The results show the model is able to replicate the aerodynamic forces and flow fields with sufficient accuracy while requiring a fraction of CFD cost.

Ion-acoustic envelope modes in a degenerate relativistic electron-ion plasma

A self-consistent relativistic two-fluid model is proposed for one-dimensional electron-ion plasma dynamics. A multiple scales perturbation technique is employed, leading to an evolution equation for the wave envelope, in the form of a nonlinear Schrödinger type equation (NLSE). The inclusion of relativistic effects is shown to introduce density-dependent factors, not present in the non-relativistic case - in the conditions for modulational instability. The role of relativistic effects on the linear dispersion laws and on envelope soliton solutions of the NLSE is discussed.

Photochemical dynamics of laser-irradiated semiconductor nanostructures

Thesis (Ph.D.)--University of Washington, 2016-08

Big bang bifurcations in von Bertalanffy’s dynamics with strong and weak Allee effects

The main purpose of this work was to study population dynamic discrete models in which the growth of the population is described by generalized von Bertalanffy's functions, with an adjustment or correction factor of polynomial type. The consideration of this correction factor is made with the aim to introduce the Allee effect. To the class of generalized von Bertalanffy's functions is identified and characterized subclasses of strong and weak Allee's functions and functions with no Allee effect. This classification is founded on the concepts of strong and weak Allee's effects to population growth rates associated. A complete description of the dynamic behavior is given, where we provide necessary conditions for the occurrence of unconditional and essential extinction types. The bifurcation structures of the parameter plane are analyzed regarding the evolution of the Allee limit with the aim to understand how the transition from strong Allee effect to no Allee effect, passing through the weak Allee effect, is realized. To generalized von Bertalanffy's functions with strong and weak Allee effects is identified an Allee's effect region, to which is associated the concepts of chaotic semistability curve and Allee's bifurcation point. We verified that under some sufficient conditions, generalized von Bertalanffy's functions have a particular bifurcation structure: the big bang bifurcations of the so-called box-within-a-box type. To this family of maps, the Allee bifurcation points and the big bang bifurcation points are characterized by the symmetric of Allee's limit and by a null intrinsic growth rate. The present paper is also a significant contribution in the framework of the big bang bifurcation analysis for continuous 1D maps and unveil their relationship with the explosion birth and the extinction phenomena.

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

Centre-of-mass motion in multi-particle Schrodinger-Newton dynamics

We investigate the implication of the nonlinear and non-local multi-particle Schrodinger-Newton equation for the motion of the mass centre of an extended multi-particle object, giving self-contained and comprehensible derivations. In particular, we discuss two opposite limiting cases. In the first case, the width of the centre-of-mass wave packet is assumed much larger than the actual extent of the object, in the second case it is assumed much smaller. Both cases result in nonlinear deviations from ordinary free Schrodinger evolution for the centre of mass. On a general conceptual level we include some discussion in order to clarify the physical basis and intention for studying the Schrodinger-Newton equation.

Nonlinear breathing modes at a defect

Recent molecular dynamics (MD) simulations of Cubero et al (1999) of a DNA duplex containing the 'rogue' base difluorotoluene (F) in place of a thymine (T) base show that breathing events can occur on the nanosecond timescale, whereas breathing events in a normal DNA duplex take place on the microsecond timescale. The main aim of this paper is to analyse a nonlinear Klein-Gordon lattice model of the DNA duplex including both nonlinear interactions between opposing bases and a defect in the interaction at one lattice site; each of which can cause localisation of energy. Solutions for a breather mode either side of the defect are derived using multiple-scales asymptotics and are pieced together across the defect to form a solution which includes the effects of the nonlinearity and the defect. We consider defects in the inter-chain interactions and in the along chain interactions. In most cases we find in-phase breather modes and/or out-of-phase breather modes, with one case displaying a shifted mode.