Applied Mathematics Commons^™

Filtered Subspace Iteration For Selfadjoint Operators, Jay Gopalakrishnan, Luka Grubišić, Jeffrey S. Ovall

Portland Institute for Computational Science Publications

We consider the problem of computing a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. A rational function of the operator is constructed such that the eigenspace of interest is its dominant eigenspace, and a subspace iteration procedure is used to approximate this eigenspace. The computed space is then used to obtain approximations of the eigenvalues of interest. An eigenvalue and eigenspace convergence analysis that considers both iteration error and dis- cretization error is provided. A realization of the proposed approach for a model second-order elliptic operator is based on a …

Computational Algorithms For Improved Representation Of The Model Error Covariance In Weak-Constraint 4d-Var, Jeremy A. Shaw

Dissertations and Theses

Four-dimensional variational data assimilation (4D-Var) provides an estimate to the state of a dynamical system through the minimization of a cost functional that measures the distance to a prior state (background) estimate and observations over a time window. The analysis fit to each information input component is determined by the specification of the error covariance matrices in the data assimilation system (DAS). Weak-constraint 4D-Var (w4D-Var) provides a theoretical framework to account for modeling errors in the analysis scheme. In addition to the specification of the background error covariance matrix, the w4D-Var formulation requires information on the model error statistics and …

Equators Have At Most Countable Many Singularities With Bounded Total Angle, Pilar Herreros, Mario Ponce, J.J.P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

For distinct points p and q in a two-dimensional Riemannian manifold, one defines their mediatrix Lpq as the set of equidistant points to p and q. It is known that mediatrices have a cell decomposition consisting of a finite number of branch points connected by Lipschitz curves. In the case of a topological sphere, mediatrices are called equators and it can benoticed that there are no branching points, thus an equator is a topological circle with possibly many Lipschitz singularities. This paper establishes that mediatrices have the radial …

Computational Methods For Asynchronous Basins, Ian H. Dinwoodie

Mathematics and Statistics Faculty Publications and Presentations

For a Boolean network we consider asynchronous updates and define the exclusive asynchronous basin of attraction for any steady state or cyclic attractor. An algorithm based on commutative algebra is presented to compute the exclusive basin. Finally its use for targeting desirable attractors by selective intervention on network nodes is illustrated with two examples, one cell signalling network and one sensor network measuring human mobility.

Transients In The Synchronization Of Asymmetrically Coupled Oscillator Arrays, Carlos E. Cantos, David K. Hammond, J.J.P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We consider the transient behavior of a large linear array of coupled linear damped harmonic oscillators following perturbation of a single element. Our work is motivated by modeling the behavior of flocks of autonomous vehicles. We first state a number of conjectures that allow us to derive an explicit characterization of the transients, within a certain parameter regime Ω. As corollaries we show that minimizing the transients requires considering non-symmetric coupling, and that within Ω the computed linear growth in N of the transients is independent of (reasonable) boundary conditions.

Accuracy Of Wave Speeds Computed From The Dpg And Hdg Methods For Electromagnetic And Acoustic Waves, Nicole Michelle Olivares

Dissertations and Theses

We study two finite element methods for solving time-harmonic electromagnetic and acoustic problems: the discontinuous Petrov-Galerkin (DPG) method and the hybrid discontinuous Galerkin (HDG) method.

The DPG method for the Helmholtz equation is studied using a test space normed by a modified graph norm. The modification scales one of the terms in the graph norm by an arbitrary positive scaling parameter. We find that, as the parameter approaches zero, better results are obtained, under some circumstances. A dispersion analysis on the multiple interacting stencils that form the DPG method shows that the discrete wavenumbers of the method are complex, explaining …

Introducing Complex Systems Analysis In High School Mathematics Using System Dynamics Modeling: A Potential Game-Changer For Mathematics Instruction, Diana Marie Fisher

Dissertations and Theses

Complex systems abound on this planet, in the composition of the human body, in ecosystems, in social interaction, in political decision-making, and more. Analytical methods allowing us to better understand how these systems operate and, consequently, to have a chance to intervene and change the undesirable behavior of some of the more pernicious systems have developed and continue to be enhanced via quickly changing technology. Some of these analytical methods are accessible by pre-college students, but have not been widely used at that level of education. Jay Forrester, the founder of one of the methodologies, System Dynamics (SD), used to …

Discrete Stability Of Dpg Methods, Ammar Harb

Dissertations and Theses

This dissertation presents a duality theorem of the Aubin-Nitsche type for discontinuous Petrov Galerkin (DPG) methods. This explains the numerically observed higher convergence rates in weaker norms. Considering the specific example of the mild-weak (or primal) DPG method for the Laplace equation, two further results are obtained. First, for triangular meshes, the DPG method continues to be solvable even when the test space degree is reduced, provided it is odd. Second, a non-conforming method of analysis is developed to explain the numerically observed convergence rates for a test space of reduced degree. Finally, for rectangular meshes, the test space is …

The Dc Algorithm & The Constrained Fermat-Torricelli Problem, Nathan Peron Lawrence, George Blikas

Student Research Symposium

The theory of functions expressible as the Difference of Convex (DC) functions has led to the development of a rich field in applied mathematics known as DC Programming.We survey the work of Pham Dinh Tao and Le Thi Hoai An in order to understand the DC Algorithm (DCA) and its use in solving clustering problems. Further, we present several other methods that generalize the DCA for any norm. These powerful tools enable researchers to reformulate objective functions, not necessarily convex, into DC Programs.

The Fermat-Torricelli problem is visited in light of convex analysis and various norms. Pierre de Fermat proposed …

Minimizing Differences Of Convex Functions With Applications To Facility Location And Clustering, Mau Nam Nguyen, R. Blake Rector, Daniel J. Giles

Mathematics and Statistics Faculty Publications and Presentations

In this paper we develop algorithms to solve generalized Fermat-Torricelli problems with both positive and negative weights and multifacility location problems involving distances generated by Minkowski gauges. We also introduce a new model of clustering based on squared distances to convex sets. Using the Nesterov smoothing technique and an algorithm for minimizing differences of convex functions called the DCA introduced by Tao and An, we develop effective algorithms for solving these problems. We demonstrate the algorithms with a variety of numerical examples.

A Multivariate Nonlinear Mixed Effects Model For Longitudinal Image Analysis: Application To Amyloid Imaging, Murat Bilgel, Jerry L. Prince, Dean F. Wong, Susan M. Resnick, Bruno M. Jedynak

Mathematics and Statistics Faculty Publications and Presentations

It is important to characterize the temporal trajectories of disease-related biomarkers in order to monitor progression and identify potential points of intervention. These are especially important for neurodegenerative diseases, as therapeutic intervention is most likely to be effective in the preclinical disease stages prior to significant neuronal damage. Neuroimaging allows for the measurement of structural, functional, and metabolic integrity of the brain at the level of voxels, whose volumes are on the order of mm³. These voxelwise measurements provide a rich collection of disease indicators. Longitudinal neuroimaging studies enable the analysis of changes in these voxelwise measures. However, …

Full State Revivals In Linearly Coupled Chains With Commensurate Eigenspectra, J. J. P. Veerman, Jovan Petrovic

Mathematics and Statistics Faculty Publications and Presentations

Coherent state transfer is an important requirement in the construction of quantum computer hardware. The state transfer can be realized by linear next-neighbour-coupled finite chains. Starting from the commensurability of chain eigenvalues as the general condition of periodic dynamics, we find chains that support full periodic state revivals. For short chains, exact solutions are found analytically by solving the inverse eigenvalue problem to obtain the coupling coefficients between chain elements. We apply the solutions to design optical waveguide arrays and perform numerical simulations of light propagation thorough realistic waveguide structures. Applications of the presented method to the realization of a …

Signal Velocity In Oscillator Arrays, Carlos E. Cantos, David K. Hammond, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We investigate a system of coupled oscillators on the circle, which arises from a simple model for behavior of large numbers of autonomous vehicles. The model considers asymmetric, linear, decentralized dynamics, where the acceleration of each vehicle depends on the relative positions and velocities between itself and a set of local neighbors. We first derive necessary and sufficient conditions for asymptotic stability, then derive expressions for the phase velocity of propagation of disturbances in velocity through this system. We show that the high frequencies exhibit damping, which implies existence of well-defined signal velocities c+>0 and c−f(x−c+t) in the direction …

Geometric Approach To Convex Subdifferential Calculus, Boris S. Mordukhovich, Nguyen Mau Nam

Mathematics and Statistics Faculty Publications and Presentations

In this paper, we develop a geometric approach to convex subdifferential calculus in finite dimensions with employing some ideas of modern variational analysis. This approach allows us to obtain natural and rather easy proofs of basic results of convex subdifferential calculus in full generality and also derive new results of convex analysis concerning optimal value/marginal functions, normals to inverse images of sets under set-valued mappings, calculus rules for coderivatives of single-valued and set-valued mappings, and calculating coderivatives of solution maps to parameterized generalized equations governed by set-valued mappings with convex graphs.

Transients Of Platoons With Asymmetric And Different Laplacians, Ivo Herman, Dan Martinec, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We consider an asymmetric control of platoons of identical vehicles with nearest-neighbor interaction. Recent results show that if the vehicle uses different asymmetries for position and velocity errors, the platoon has a short transient and low overshoots. In this paper we investigate the properties of vehicles with friction. To achieve consensus, an integral part is added to the controller, making the vehicle a third-order system. We show that the parameters can be chosen so that the platoon behaves as a wave equation with different wave velocities. Simulations suggest that our system has a better performance than other nearest-neighbor scenarios. Moreover, …

On The Coupling Of Dpg And Bem, Thomas Führer, Norbert Heuer, Michael Karkulik

Mathematics and Statistics Faculty Publications and Presentations

We develop and analyze strategies to couple the discontinuous Petrov-Galerkin method with optimal test functions to (i) least-squares boundary elements and (ii) various variants of standard Galerkin boundary elements. An essential feature of our method is that, despite the use of boundary integral equations, optimal test functions have to be computed only locally. We apply our findings to a standard transmission problem in full space and present numerical experiments to validate our theory.

On A Convex Set With Nondifferentiable Metric Projection, Shyan S. Akmal, Nguyen Mau Nam, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

A remarkable example of a nonempty closed convex set in the Euclidean plane for which the directional derivative of the metric projection mapping fails to exist was constructed by A. Shapiro. In this paper, we revisit and modify that construction to obtain a convex set with smooth boundary which possesses the same property.

A Tent Pitching Scheme Motivated By Friedrichs Theory, Jay Gopalakrishnan, Peter Monk, Paulina Sepulveda

Mathematics and Statistics Faculty Publications and Presentations

Certain Friedrichs systems can be posed on Hilbert spaces normed with a graph norm. Functions in such spaces arising from advective problems are found to have traces with a weak continuity property at points where the inflow and outflow boundaries meet. Motivated by this continuity property, an explicit space-time finite element scheme of the tent pitching type, with spaces that conform to the continuity property, is designed. Numerical results for a model one-dimensional wave propagation problem are presented.

Dynamics Of Locally Coupled Oscillators With Next-Nearest-Neighbor Interaction, J. Herbrych, A. G. Chazirakis, N. Christakis, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

A theoretical description of decentralized dynamics within linearly coupled, one-dimensional oscillators (agents) with up to next-nearest-neighbor interaction is given. Conditions for stability of such system are presented. Our results indicate that the stable systems have response that grow at least linearly in the system size. We give criteria when this is the case. The dynamics of these systems can be described with traveling waves with strong damping in the high frequencies. Depending on the system parameters, two types of solutions have been found: damped oscillations and reflectionless waves. The latter is a novel result and a feature of systems with …

Monomials And Basin Cylinders For Network Dynamics, Daniel Austin, Ian H. Dinwoodie

Mathematics and Statistics Faculty Publications and Presentations

We describe methods to identify cylinder sets inside a basin of attraction for Boolean dynamics of biological networks. Such sets are used for designing regulatory interventions that make the system evolve towards a chosen attractor, for example initiating apoptosis in a cancer cell. We describe two algebraic methods for identifying cylinders inside a basin of attraction, one based on the Groebner fan that finds monomials that define cylinders and the other on primary decomposition. Both methods are applied to current examples of gene networks.

Stability Of A Circular System With Multiple Asymmetric Laplacians, Ivo Herman, Dan Martinec, J. J. P. Veerman, Michael Sebek

Mathematics and Statistics Faculty Publications and Presentations

We consider an asymptotic stability of a circular system where the coupling Laplacians are different for each state used for synchronization. It is shown that there must be a symmetric coupling in the output state to guarantee the stability for agents with two integrators in the open loop. Systems with agents having three or more integrators cannot be stabilized by any coupling. In addition, recent works in analysis of a scaling in vehicular platoons relate the asymptotic stability of a circular system to a string stability. Therefore, as confirmed by simulations in the paper, our results have an application also …

Periodic State Revivals In Commensurate Waveguide Arrays, Jovan Petrovic, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

Emerging optical and quantum computers require hardware capable of coherent transport of and operations on quantum states. Here, we investigate finite optical waveguide arrays with linear coupling as means of efficient and compact coherent state transfer. Coherent transfer with periodic state revivals is enabled by engineering coupling coefficients between neighbouring waveguides to yield commensurate eigenvalue spectrum. Particular cases of finite arrays have been actively studied to achieve the perfect state transfer by mirroring the input into the output state.

We explore a much wider scope of coherent propagation and revivals of both the state amplitude and phase. We analytically solve …

Nonsmooth Algorithms And Nesterov's Smoothing Technique For Generalized Fermat-Torricelli Problems, Nguyen Mau Nam, Nguyen Thai An, R. Blake Rector, Jie Sun

Mathematics and Statistics Faculty Publications and Presentations

We present algorithms for solving a number of new models of facility location which generalize the classical Fermat--Torricelli problem. Our first approach involves using Nesterov's smoothing technique and the minimization majorization principle to build smooth approximations that are convenient for applying smooth optimization schemes. Another approach uses subgradient-type algorithms to cope directly with the nondifferentiability of the cost functions. Convergence results of the algorithms are proved and numerical tests are presented to show the effectiveness of the proposed algorithms.

Transients In The Synchronization Of Oscillator Arrays, Carlos E. Cantos, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

The purpose of this note is threefold. First we state a few conjectures that allow us to rigorously derive a theory which is asymptotic in N (the number of agents) that describes transients in large arrays of (identical) linear damped harmonic oscillators in R with completely decentralized nearest neighbor interaction. We then use the theory to establish that in a certain range of the parameters transients grow linearly in the number of agents (and faster outside that range). Finally, in the regime where this linear growth occurs we give the constant of proportionality as a function of the signal velocities …

Conditional Tests On Basins Of Attraction With Finite Fields, Ian H. Dinwoodie

Mathematics and Statistics Faculty Publications and Presentations

An iterative method is given for computing the polynomials that vanish on the basin of attraction of a steady state in discrete polynomial dynamics with finite field coefficients. The algorithm is applied to dynamics of a T cell survival network where it is used to compare transition maps conditional on a basin of attraction.

A Posteriori Estimates Using Auxiliary Subspace Techniques, Harri Hakula, Michael Neilan, Jeffrey S. Ovall

Mathematics and Statistics Faculty Publications and Presentations

A posteriori error estimators based on auxiliary subspace techniques for second order elliptic problems in Rd (d ≥ 2) are considered. In this approach, the solution of a global problem is utilized as the error estimator. As the continuity and coercivity of the problem trivially leads to an efficiency bound, the main focus of this paper is to derive an analogous effectivity bound and to determine the computational complexity of the auxiliary approximation problem. With a carefully chosen auxiliary subspace, we prove that the error is bounded above by the error estimate up to oscillation terms. In addition, we show …

Regularity Of Mediatrices In Surfaces, Pilar Herreros, Mario Ponce, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

For distinct points p and q in a two-dimensional Riemannian manifold, one defines their mediatrix Lpq as the set of equidistant points to p and q. It is known that mediatrices have a cell decomposition consisting of a finite number of branch points connected by Lipschitz curves. This paper establishes additional geometric regularity properties of mediatrices. We show that mediatrices have the radial linearizability property, which implies that at each point they have a geometrically defined derivative in the branching directions. Also, we study the particular case of mediatrices on spheres, by showing that they are Lipschitz simple closed curves …

Stability Of Large Flocks: An Example, J. J. P. Veerman, F. M. Tangerman

Mathematics and Statistics Faculty Publications and Presentations

The movement of a flock with a single leader (and a directed path from it to every agent) can be stabilized. Nonetheless for large flocks perturbations in the movement of the leader may grow to a considerable size as they propagate throughout the flock and before they die out over time. As an example we consider a string of N+1 oscillators moving in the line. Each one `observes' the relative velocity and position of only its nearest neighbors. This information is then used to determine its own acceleration. Now we fix all parameters except the number of oscillators. We then …

Asymptotic Reliability Rheory Of K-Out-Of-N Systems, Nuria Torrado, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We formulate a theory that allows us to formulate a simple criterion that ensures that two k-out-of-n systems A and are not ordered. If the systems fail the criterion, it does not follow they are ordered. Thus the theory only serves to avoid some a priori useless comparisons: when neither A nor can be said to be better than the other. The power of the theory lies in its wide potential applicability: the assumptions involve very weak estimates on the asymptotic behavior (as t→0 and as t→∞) of the constituent survival probabilities. We include examples.

Symmetric Nonconforming Mixed Finite Elements For Linear Elasticity, Jay Gopalakrishnan, Johnny Guzmán

Mathematics and Statistics Faculty Publications and Presentations

We present a family of mixed methods for linear elasticity that yield exactly symmetric, but only weakly conforming, stress approximations. The method is presented in both two and three dimensions (on triangular and tetrahedral meshes). The method is efficiently implementable by hybridization. The degrees of freedom of the Lagrange multipliers, which approximate the displacements at the faces, solve a symmetric positive-definite system. The design and analysis of this method is motivated by a new set of unisolvent degrees of freedom for symmetric polynomial matrices. These new degrees of freedom are also used to give a new simple calculation of the …