Statistics and Probability Commons^™

Balanced Pod Algorithm For Robust Control Design For Linear Distributed Parameter Systems, John R. Singler, Belinda A. Batten

Mathematics and Statistics Faculty Research & Creative Works

A mathematical model of a physical system is never perfect; therefore, robust control laws are necessary for guaranteed stabilization of the nominal model and also "nearby" systems, including hopefully the actual physical system. We consider the computation of a robust control law for large-scale finite dimensional linear systems and a class of linear distributed parameter systems. The controller is robust with respect to left coprime factor perturbations of the nominal system. We present an algorithm based on balanced proper orthogonal decomposition to compute the nonstandard features of this robust control law. Numerical results are presented for a convection diffusion partial …

Computational Issues In Sensitivity Analysis For 1d Interface Problems, L. G. Davis, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

This paper is concerned with the construction of accurate and e cient computational algorithms for the numerical approximation of sensitivities with respect to a parameter dependent interface location. Motivated by sensitivity analysis with respect to piezoelectric actuator placement on an

Boundary Data Maps For Schrödinger Operators On A Compact Interval, Stephen L. Clark, Fritz Gesztesy, M. Mitrea

Mathematics and Statistics Faculty Research & Creative Works

We provide a systematic study of boundary data maps, that is, 2 x 2 matrix-valued Dirichlet-to-Neumann and more generally, Robin-to-Robin maps, associated with one-dimensional Schrödinger operators on a compact interval [0, R] with separated boundary conditions at 0 and R. Most of our results are formulated in the non-self-adjoint context. Our principal results include explicit representations of these boundary data maps in terms of the resolvent of the underlying Schrödinger operator and the associated boundary trace maps, Krein-type resolvent formulas relating Schrödinger operators corresponding to different (separated) boundary conditions, and a derivation of the Herglotz property of boundary data maps …

High Accuracy Combination Method For Solving The Systems Of Nonlinear Volterra Integral And Integro-Differential Equations With Weakly Singular Kernels Of The Second Kind, Xiaoming He, Lu Pan, Tao Lü

Mathematics and Statistics Faculty Research & Creative Works

This paper presents a high accuracy combination algorithm for solving the systems of nonlinear Volterra integral and integro-differential equations with weakly singular kernels of the second kind. Two quadrature algorithms for solving the systems are discussed, which possess high accuracy order and the asymptotic expansion of the errors. By means of combination algorithm, we may obtain a numerical solution with higher accuracy order than the original two quadrature algorithms. Moreover an a posteriori error estimation for the algorithm is derived. Both of the theory and the numerical examples show that the algorithm is effective and saves storage capacity and computational …

The Hodrick-Prescott Filter: A Special Case Of Penalized Spline Smoothing, Robert Paige L., A. A. Trindade

Mathematics and Statistics Faculty Research & Creative Works

We prove that the Hodrick-Prescott Filter (HPF), a commonly used method for smoothing econometric time series, is a special case of a linear penalized spline model with knots placed at all observed time points (except the first and last) and uncorrelated residuals. This equivalence then furnishes a rich variety of existing data-driven parameter estimation methods, particularly restricted maximum likelihood (REML) and generalized cross-validation (GCV). This has profound implications for users of HPF who have hitherto typically relied on subjective choice, rather than estimation, for the smoothing parameter. By viewing estimates as roots of an appropriate quadratic estimating equation, we also …

Incorporating Genome Annotation In The Statistical Analysis Of Genomic And Epigenomic Tiling Array Data, Gayla R. Olbricht

Mathematics and Statistics Faculty Research & Creative Works

"A wealth of information and technologies are currently available for the genomewide investigation of many types of biological phenomena. Genomic annotation databases provide information about the DNA sequence of a particular organism and give locations of different types of genomic elements, such as the exons and introns of genes. Microarrays are a powerful type of technology that make use of DNA sequence information to investigate different types of biological phenomena on a genome-wide level. Tiling arrays are a unique type of microarray that provide unbiased, highdensity coverage of a genomic region, making them well suited for many applications, such as …

Gronwall-Ouiang-Type Integral Inequalities On Time Scales, Ailian Liu, Martin Bohner

Mathematics and Statistics Faculty Research & Creative Works

We present several Gronwall-OuIang-type integral inequalities on time scales. Firstly, an OuIang inequality on time scales is discussed. Then we extend the Gronwall-type inequalities to multiple integrals. Some special cases of our results contain continuous Gronwall-type inequalities and their discrete analogues. Several examples are included to illustrate our results at the end.

Weyl-Titchmarsh Theory For Hamiltonian Dynamic Systems, Shurong Sun, Shaozhu Chen, Martin Bohner

Mathematics and Statistics Faculty Research & Creative Works

We establish the Weyl-Titchmarsh theory for singular linear Hamiltonian dynamic systems on a time scale T , which allows one to treat both continuous and discrete linear Hamiltonian systems as special cases for T= ℝ and T= ℤ within one theory and to explain the discrepancies between these two theories. This paper extends the Weyl-Titchmarsh theory and provides a foundation for studying spectral theory of Hamiltonian dynamic systems. These investigations are part of a larger program which includes the following: (i) M(λ) theory for singular Hamiltonian systems, (ii) on the spectrum of Hamiltonian systems, (iii) on boundary value problems for …

Approximating Stationary Statistical Properties, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

It is well-known that physical laws for large chaotic dynamical systems are revealed statistically. Many times these statistical properties of the system must be approximated numerically. the main contribution of this manuscript is to provide simple and natural criterions on numerical methods (temporal and spatial discretization) that are able to capture the stationary statistical properties of the underlying dissipative chaotic dynamical systems asymptotically. the result on temporal approximation is a recent finding of the author, and the result on spatial approximation is a new one. Applications to the infinite Prandtl number model for convection and the barotropic quasi-geostrophic model are …

A Comparison Of Balanced Truncation Methods For Closed Loop Systems, John R. Singler, Belinda A. Batten

Mathematics and Statistics Faculty Research & Creative Works

Real-time control of a physical system necessitates controllers that are low order. In this paper, we compare two balanced truncation methods as a means of designing low order compensators for partial differential equation (PDE) systems. The first method is the application of balanced truncation to the compensator dynamics, rather than the state dynamics, as was done in cite{Skelton:1984}. The second method, LQG balanced truncation, applies the balancing technique to the Riccati operators obtained from a specific LQG design. We discuss snapshot-based algorithms for constructing the reduced order compensators and present numerical results for a two dimensional convection diffusion PDE system.

Bilinear Immersed Finite Elements For Interface Problems, Xiaoming He

Mathematics and Statistics Faculty Research & Creative Works

In this dissertation we discuss bilinear immersed finite elements (IFE) for solving interface problems. The related research works can be categorized into three aspects: (1) the construction of the bilinear immersed finite element spaces; (2) numerical methods based on these IFE spaces for solving interface problems; and (3) the corresponding error analysis. All of these together form a solid foundation for the bilinear IFEs.

The research on immersed finite elements is motivated by many real world applications, in which a simulation domain is often formed by several materials separated from each other by curves or surfaces while a mesh independent …

A Snapshot Algorithm For Linear Feedback Flow Control Design, Benjamin T. Dickinson, Belinda A. Batten, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

The control of fluid flows has many applications. For micro air vehicles, integrated flow control designs could enhance flight stability by mitigating the effect of destabilizing air flows in their low Reynolds number regimes. However, computing model based feedback control designs can be challenging due to high dimensional discretized flow models. In this work, we investigate the use of a snapshot algorithm proposed in Ref. 1 to approximate the feedback gain operator for a linear incompressible unsteady flow problem on a bounded domain. The main component of the algorithm is obtaining solution snapshots of certain linear flow problems. Numerical results …

A Proper Orthogonal Decomposition Approach To Approximate Balanced Truncation Of Infinite Dimensional Linear Systems, John R. Singler, Belinda A. Batten

Mathematics and Statistics Faculty Research & Creative Works

We extend a method for approximate balanced reduced order model derivation for finite dimensional linear systems developed by Rowley (Int. J. Bifur. Chaos Appl. Sci. Eng. 15(3) (2005), pp. 997-1013) to infinite dimensional systems. The algorithm is related to standard balanced truncation, but includes aspects of the proper orthogonal decomposition in its computational approach. The method can be also applied to nonlinear systems. Numerical results are presented for a convection diffusion system.

Effect Of Dlk1 And Rtl1 But Not Meg3 Or Meg8 On Muscle Gene Expression In Callipyge Lambs, Jolena N. Fleming-Waddell, Gayla R. Olbricht, Tasia M. Taxis, Jason D. White, Tony Vuocolo, Bruce A. Craig, Ross L. Tellam, Mike K. Neary, Noelle E. Cockett, Christopher A. Bidwell

Mathematics and Statistics Faculty Research & Creative Works

Callipyge sheep exhibit extreme postnatal muscle hypertrophy in the loin and hindquarters as a result of a single nucleotide polymorphism (SNP) in the imprinted DLK1-DIO3 domain on ovine chromosome 18. The callipyge SNP up-regulates the expression of surrounding transcripts when inherited in cis without altering their allele-specific imprinting status. The callipyge phenotype exhibits polar overdominant inheritance since only paternal heterozygous animals have muscle hypertrophy. Two studies were conducted profiling gene expression in lamb muscles to determine the down-stream effects of over-expression of paternal allele-specific DLK1 and RTL1 as well as maternal allele-specific MEG3, RTL1AS and MEG8, using Affymetrix bovine expression …

A Finite Element Splitting Extrapolation For Second Order Hyperbolic Equations, Xiaoming He, Tao Lü

Mathematics and Statistics Faculty Research & Creative Works

Splitting extrapolation is an efficient technique for solving large scale scientific and engineering problems in parallel. This article discusses a finite element splitting extrapolation for second order hyperbolic equations with time-dependent coefficients. This method possesses a higher degree of parallelism, less computational complexity, and more flexibility than Richardson extrapolation while achieving the same accuracy. By means of domain decomposition and isoparametric mapping, some grid parameters are chosen according to the problem. The multiparameter asymptotic expansion of the d-quadratic finite element error is also established. The splitting extrapolation formulas are developed from this expansion. An approximation with higher accuracy on a …

A Uniformly Dissipative Scheme For Stationary Statistical Properties Of The Infinite Prandtl Number Model, Wenfang (Wendy) Cheng, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

The purpose of this short communication is to announce that a class of numerical schemes, uniformly dissipative approximations, which uniformly preserve the dissipativity of the continuous infinite dimensional dissipative complex (chaotic) systems possess desirable properties in terms of approximating stationary statistics properties. in particular, the stationary statistical properties of these uniformly dissipative schemes converge to those of the continuous system at vanishing mesh size. the idea is illustrated on the infinite Prandtl number model for convection and semi-discretization in time, although the general strategy works for a broad class of dissipative complex systems and fully discretized approximations. as far as …

A Semi-Implicit Scheme For Stationary Statistical Properties Of The Infinite Prandtl Number Model, Wenfang Cheng, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

We propose a semisecret in time semi-implicit numerical scheme for the infinite Prandtl model for convection. Besides the usual finite time convergence, this scheme enjoys the additional highly desirable feature that the stationary statistical properties of the scheme converge to those of the infinite Prandtl number model at vanishing time stop. One of the key characteristics of the scheme is that it preserves the dissipativity of the infinite Prandtl number model uniformly in terms of the time stop. So far as wo know, this is the first rigorous result on convergence of stationary statistical properties of numerical schemes for infinite …

Confluent Mappings And Arc Kelley Continua, W. J. Charatonik, Janusz R. Prajs, J. J. Charatonik

Mathematics and Statistics Faculty Research & Creative Works

A Kelley continuum X, also called a continuum with the property of Kelley, such that, for each p X, each subcontinuum K containing p is approximated by arc-wise connected continua containing p, is called an arc Kelley continuum. A continuum homeomorphic to the inverse limit of locally connected continua with confluent bonding maps is said to be confluently LC-representable. The main subject of the paper is a study of deep connections between the arc Kelley continua and confluent mappings. It is shown that if a continuum X admits, for each ε > 0, a confluent ε-mapping onto a(n) (arc) Kelley continuum, …

Generically There Is But One Self Homeomorphism Of The Cantor Set, Ethan Akin, Eli Glasner, Benjamin Weiss

Mathematics and Statistics Faculty Research & Creative Works

We describe a self-homeomorphism R of the Cantor set X and then show that its conjugacy class in the Polish group H(X) of all homeomorphisms of X forms a dense Gδ subset of H(X). We also provide an example of a locally compact, second countable topological group which has a dense conjugacy class. © 2008 American Mathematical Society.

The Detection Of Unsteady Flow Separation With Bioinspired Hair-Cell Sensors, Benjamin T. Dickinson, John R. Singler, Belinda A. Batten

Mathematics and Statistics Faculty Research & Creative Works

Biologists hypothesize that thousands of micro-scale hairs found on bat wings function as a network of air-flow sensors as part of a biological feedback flow control loop. In this work, we investigate hair-cell sensors as a means of detecting flow features in an unsteady separating flow over a cylinder. Individual hair-cell sensors were modeled using an Euler-Bernoulli beam equation forced by the fluid flow. When multiple sensor simulations are combined into an array of hair-cells, the response is shown to detect the onset and span of flow reversal, the upstream movement of the point of zero wall shear-stress, and the …

Approximate Low Rank Solutions Of Lyapunov Equations Via Proper Orthogonal Decomposition, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

We present an algorithm to approximate the solution Z of a stable Lyapunov equation AZ + ZA* + BB* = 0 using proper orthogonal decomposition (POD). This algorithm is applicable to large-scale problems and certain infinite dimensional problems as long as the rank of B is relatively small. In the infinite dimensional case, the algorithm does not require matrix approximations of the operators A and B. POD is used in a systematic way to provide convergence theory and simple a priori error bounds.

Stationary Statistical Properties Of Rayleigh-Bénard Convection At Large Prandtl Number, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

This is the third in a series of our study of Rayleigh-Bénard convection at large Prandtl number. Here we investigate whether stationary statistical properties of the Boussinesq system for Rayleigh-Bénard convection at large Prandtl number are related to those of the infinite Prandtl number model for convection that is formally derived from the Boussinesq system via setting the Prandtl number to infinity. We study asymptotic behavior of stationary statistical solutions, or in-variant measures, to the Boussinesq system for Rayleigh-Bénard convection at large Prandtl number. in particular, we show that the invariant measures of the Boussinesq system for Rayleigh-Bénard convection converge …

Evaluating Statistical Methods Using Plasmode Data Sets In The Age Of Massive Public Databases: An Illustration Using False Discovery Rates, Gary L. Gadbury, Qinfang Xiang, Lin Yang, Stephen Barnes, Grier P. Page, David B. Allison

Mathematics and Statistics Faculty Research & Creative Works

Plasmode is a term coined several years ago to describe data sets that are derived from real data but for which some truth is known. Omic techniques, most especially microarray and genome wide association studies, have catalyzed a new zeitgeist of data sharing that is making data and data sets publicly available on an unprecedented scale. Coupling such data resources with a science of plasmode use would allow statistical methodologists to vet proposed techniques empirically (as opposed to only theoretically) and with data that are by definition realistic and representative. We illustrate the technique of empirical statistics by consideration of …

Property Of Kelley For The Cartesian Products And Hyperspaces, W. J. Charatonik, J. J. Charatonik

Mathematics and Statistics Faculty Research & Creative Works

A continuum X having the property of Kelley is constructed such that neither X × [0, 1], nor the hyperspace C(X), nor small Whitney levels in C(X) have the property of Kelley. This answers several questions asked in the literature.

Iterated Oscillation Criteria For Delay Dynamic Equations Of First Order, B. Karpuz, O. Öcalan, Martin Bohner

Mathematics and Statistics Faculty Research & Creative Works

We obtain new sufficient conditions for the oscillation of all solutions of first-order delay dynamic equations on arbitrary time scales, hence combining and extending results for corresponding differential and difference equations. Examples, some of which coincide with well-known results on particular time scales, are provided to illustrate the applicability of our results.

On The Asymptotic Integration Of Nonlinear Dynamic Equations, Elvan Akin, Smail Djebali, Toufik Moussaoui, Martin Bohner

Mathematics and Statistics Faculty Research & Creative Works

The purpose of this paper is to study the existence and asymptotic behavior of solutions to a class of second-order nonlinear dynamic equations on unbounded time scales. Four different results are obtained by using the Banach fixed point theorem, the Boyd and Wong fixed point theorem, the Leray-Schauder nonlinear alternative, and the Schauder fixed point theorem. For each theorem, an illustrative example is presented. The results provide unification and some extensions in the time scale setup of the theory of asymptotic integration of nonlinear equations both in the continuous and discrete cases

Transition To Turbulence, Small Disturbances, And Sensitivity Analysis Ii: The Navier-Stokes Equations, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

Recent research has shown that small disturbances in the linearized Navier-Stokes equations cause large energy growth in solutions. Although many researchers believe that this interaction triggers transition to turbulence in flow systems, the role of the nonlinearity in this process has not been thoroughly investigated. This paper is the second of a two part work in which sensitivity analysis is used to study the effects of small disturbances on the transition process. In the first part, sensitivity analysis was used to predict the effects of a small disturbance on solutions of a motivating problem, a highly sensitive one dimensional Burgers' …

Transition To Turbulence, Small Disturbances, And Sensitivity Analysis I: A Motivating Problem, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

For over 100 years, researchers have attempted to predict transition to turbulence in fluid flows by analyzing the spectrum of the linearized Navier-Stokes equations. However, for many simple flows this approach fails to match experimental results. Recently, new scenarios for transition have been proposed that are based on the interaction of the linearized equations of motion with small disturbances to the flow system. These new "mostly linear" theories have increased our understanding of the transition process, but the role of nonlinearity has not been explored in detail. This paper is the first of a two part work in which sensitivity …

Asymptotic Behavior Of The Global Attractors To The Boussinesq System For Rayleigh-Bénard Convection At Large Prandtl Number, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

We study asymptotic behavior of the global attractors to the Boussinesq system for Rayleigh-Bénard convection at large Prandtl number. in particular, we show that the global attractors to the Boussinesq system for Rayleigh-Bénard convection converge to that of the infinite-Prandtl- number model for convection as the Prandtl number approaches infinity. This offers partial justification of the infinite-Prandtl-number model for convection as a valid simplified model for convection at large Prandtl number even in the long-time regime. © 2006 Wiley Periodicals, Inc.

Double Integral Calculus Of Variations On Time Scales, Gusein Sh. Guseinov, Martin Bohner

Mathematics and Statistics Faculty Research & Creative Works

We consider a version of the double integral calculus of variations on time scales, which includes as special cases the classical two-variable calculus of variations and the discrete two-variable calculus of variations. Necessary and sufficient conditions for a local extremum are established, among them an analogue of the Euler-Lagrange equation.