Mathematics Commons^™

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 …

Inverse Limits With Upper Semi-Continuous Set Valued Bonding Functions: An Example, Christopher David Jacobsen

Masters Theses

"While there is a wealth of information pertaining to inverse limits with single valued bonding maps, comparatively little is known about inverse limits with upper semi-continuous set valued bonding functions. In order to add somewhat to the communal knowledge on the subject, this paper provides an example of an inverse limit with a single upper semi-continuous set valued bonding function. It is then shown that the space is a continuum, and its structure is examined via its arc components and through various of its properties, such as dimension and decomposability"--Abstract, page iii.

Analysis Of Transient Growth In Iterative Learning Control Using Pseudospectra, Douglas A. Bristow, John R. Singler

Mechanical and Aerospace Engineering Faculty Research & Creative Works

In this paper we examine the problem of transient growth in Iterative Learning Co ntrol (ILC). Transient growth is generally avoided in design by using robust monotonic convergence (RMC) criteria. However, RMC leads to fundamental performance limitations. We consider the possibility of allowing safe transient growth in ILC algorithms as a means to circumvent these limitations. Here the pseudospectra is used for the first time to study transient growth in ILC. Basic properties of the pseudospectra that are relevant to the ILC problem are presented. Two ILC design problems are considered and examined using pseduospectra. The pseudospectra provides new results …

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 …

Risk Matrix Input Data Biases, Eric D. Smith, William T. Siefert, David Drain

Engineering Management and Systems Engineering Faculty Research & Creative Works

Risk matrices used in industry characterize particular risks in terms of the likelihood of occurrence, and the consequence of the actualized risk. Human cognitive bias research led by Daniel Kahneman and Amos Tversky exposed systematic translations of objective probability and value as judged by human subjects. Applying these translations to the risk matrix allows the formation of statistical hypotheses of risk point placement biases. Industry-generated risk matrix data reveals evidence of biases in the judgment of likelihood and consequence-principally, likelihood centering, a systematic increase in consequence, and a diagonal bias. Statistical analyses are conducted with linear regression, normal distribution fitting, …

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 …

Holomorphic Extensions In Toric Varieties, Malgorzata Aneta Marciniak

Doctoral Dissertations

"The dissertation describes the Hartogs and the Hartogs-Bochner extension phenomena in smooth toric varieties and their connection with the first cohomology group with compact support and sheaf coefficients. The affirmative and negative results are proved for toric surfaces and for line bundles over toric varieties using topological, analytic, and algebraic methods"--Abstract, page iii.

The Kalman Filter On Time Scales, Nicholas J. Wintz

Doctoral Dissertations

"In this work, we study concepts in optimal control for dynamic equations on time scales, which unfies the discrete and continuous cases. After a brief introduction of dynamic equations on time scales, we will examine controllability and observability for linear systems. Then we construct and solve the linear quadratic regulator for arbitrary time scales. Here, we seek to find an optimal control that minimizes a given cost function associated with a linear system. We will find such an input under two different settings; when the final state is fixed and when it is free. Later, we extend these results to …

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.