Mathematics Commons^™

Roger Temam On The Occasion Of His 70th Birthday, Claude Michel Brauner, Danielle Hilhorst, Alain Miranville, Shouhong Wang, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

No abstract provided.

Approximation Of Stationary Statistical Properties Of Dissipative Dynamical Systems: Time Discretization, Xiaoming Wang

Mathematics and Statistics Faculty Research & Creative Works

We consider temporal approximation of stationary statistical properties of dissipative infinite-dimensional dynamical systems. We demonstrate that stationary statistical properties of the time discrete approximations, i.e., numerical scheme, converge to those of the underlying continuous dissipative infinite-dimensional dynamical system under three very natural assumptions as the time step approaches zero. the three conditions that are sufficient for the convergence of the stationary statistical properties are: (1) uniform dissipativity of the scheme in the sense that the union of the global attractors for the numerical approximations is pre-compact in the phase space; (2) convergence of the solutions of the numerical scheme to …

Finite Element Approximations For Stokes-Darcy Flow With Beavers-Joseph Interface Conditions, Yanzhao Cao, Max Gunzburger, Xiaolong Hu, Fei Hua, Xiaoming Wang, Weidong Zhao

Mathematics and Statistics Faculty Research & Creative Works

Numerical solutions using finite element methods are considered for transient flow in a porous medium coupled to free flow in embedded conduits. Such situations arise, for example, for groundwater flows in karst aquifers. the coupled flow is modeled by the Darcy equation in a porous medium and the Stokes equations in the conduit domain. on the interface between the matrix and conduit, Beavers-Joseph interface conditions, instead of the simplified Beavers-Joseph-Saffman conditions, are imposed. Convergence and error estimates for finite element approximations are obtained. Numerical experiments illustrate the validity of the theoretical results. © 2010 Society for Industrial and Applied Mathematics.

Dynamical Laws Of The Coupled Gross-Pitaevskii Equations For Spin-1 Bose-Einstein Condensates, Weizhu Bao, Yanzhi Zhang

Mathematics and Statistics Faculty Research & Creative Works

In this paper, we derive analytically the dynamical laws of the coupled Gross- Pitaevskii equations (CGPEs) without/with an angular momentum rotation term and an external magnetic field for modelling nonrotating/rotating spin-1 Bose-Eintein condensates. We prove the conservation of the angular momentum expectation when the external trapping potential is radially symmetric in two dimensions and cylindrically symmetric in three dimensions; obtain a system of first order ordinary differential equations (ODEs) governing the dynamics of the density of each component and solve the ODEs analytically in a few cases; derive a second order ODE for the dynamics of the condensate width and …

Nonlinear Model Reduction Using Group Proper Orthogonal Decomposition, Benjamin T. Dickinson, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

We propose a new method to reduce the cost of computing nonlinear terms in projec- tion based reduced order models with global basis functions. We develop this method by extending ideas from the group nite element (GFE) method to proper orthogonal decomposition (POD) and call it the group POD method. Here, a scalar two-dimensional Burgers' equation is used as a model problem for the group POD method. Numerical results show that group POD models of Burgers' equation are as accurate and are computationally more e cient than standard POD models of Burgers' equation.

Optimality Of Balanced Proper Orthogonal Decomposition For Data Reconstruction, John R. Singler

Mathematics and Statistics Faculty Research & Creative Works

Proper orthogonal decomposition (POD) finds an orthonormal basis yielding an optimal reconstruction of a given dataset. We consider an optimal data reconstruction problem for two general datasets related to balanced POD, which is an algorithm for balanced truncation model reduction for linear systems. We consider balanced POD outside of the linear systems framework, and prove that it solves the optimal data reconstruction problem. the theoretical result is illustrated with an example.

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.

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 …

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 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.