Mathematics Commons^™

Partial Expansion Of A Lipschitz Domain And Some Applications, Weifeng Qiu, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We show that a Lipschitz domain can be expanded solely near a part of its boundary, assuming that the part is enclosed by a piecewise C¹ curve. The expanded domain as well as the extended part are both Lipschitz. We apply this result to prove a regular decomposition of standard vector Sobolev spaces with vanishing traces only on part of the boundary. Another application in the construction of low-regularity projectors into finite element spaces with partial boundary conditions is also indicated

A Class Of Discontinuous Petrov–Galerkin Methods. Part Iii: Adaptivity, Leszek Demkowicz, Jay Gopalakrishnan, Antti H. Niemi

Mathematics and Statistics Faculty Publications and Presentations

We continue our theoretical and numerical study on the Discontinuous Petrov–Galerkin method with optimal test functions in context of 1D and 2D convection-dominated diffusion problems and hp-adaptivity. With a proper choice of the norm for the test space, we prove robustness (uniform stability with respect to the diffusion parameter) and mesh-independence of the energy norm of the FE error for the 1D problem. With hp-adaptivity and a proper scaling of the norms for the test functions, we establish new limits for solving convection-dominated diffusion problems numerically: for 1D and for 2D problems. The adaptive process is fully automatic and starts …

Benchmark Results For Testing Adaptive Finite Element Eigenvalue Procedures Ii (Cluster Robust Eigenvector And Eigenvalue Estimates), Stefano Giani, Luka Grubisic, Jeffrey S. Ovall

Mathematics and Statistics Faculty Publications and Presentations

As a model benchmark problem for this study we consider a highly singular transmission type eigenvalue problem which we study in detail both analytically as well as numerically. In order to justify our claim of cluster robust and highly accurate approximation of a selected groups of eigenvalues and associated eigenfunctions, we give a new analysis of a class of direct residual eigenspace/vector approximation estimates. Unlike in the first part of the paper, we now use conforming higher order finite elements, since the canonical choice of an appropriate norm to measure eigenvector approximation by discontinuous Galerkin methods is an open problem.

Supporting Implementation Of The Common Core State Standards For Mathematics: Recommendations For Professional Development, Paola Sztajn, Karen A. Marrongelle, Peg Smith, Bonnie L. Melton

Mathematics and Statistics Faculty Publications and Presentations

In 2010, the National Governor’s Association and the Council of Chief State School Officers published the Common Core State Standards for Mathematics (CCSSM) and to date, 44 states, the District of Columbia, and the U.S. Virgin Islands have adopted the document. These content and practice standards, which specify what students are expected to understand and be able to do in K-12 mathematics, represent a significant departure from what mathematics is currently taught in most classrooms and how it is taught. Developing teachers’ capacity to enact these new standards in ways that support the intended student learning outcomes will require considerable …

Short-Term Plasticity At The Schaffer Collateral: A New Model With Implications For Hippocampal Processing, Andrew Hamilton Toland

Dissertations and Theses

A new mathematical model of short-term synaptic plasticity (STP) at the Schaffer collateral is introduced. Like other models of STP, the new model relates short-term synaptic plasticity to an interaction between facilitative and depressive dynamic influences. Unlike previous models, the new model successfully simulates facilitative and depressive dynamics within the framework of the synaptic vesicle cycle. The novelty of the model lies in the description of a competitive interaction between calcium-sensitive proteins for binding sites on the vesicle release machinery. By attributing specific molecular causes to observable presynaptic effects, the new model of STP can predict the effects of specific …

A Review On Convolutions Of Gamma Random Variables, Baha-Eldin Khaledi, Subhash C. Kochar

Mathematics and Statistics Faculty Publications and Presentations

Due to its wide range of applications, the topic of the distribution theory of convolutions of Gamma random variables has attracted the attention of many researchers. In this paper we review some of the latest developments on this problem.

Reliable A-Posteriori Error Estimators For Hp-Adaptive Finite Element Approximations Of Eigenvalue/Leigenvector Problems, Stefano Giani, Luka Grubisic, Jeffrey S. Ovall

Mathematics and Statistics Faculty Publications and Presentations

We present reliable a-posteriori error estimates for hp-adaptive finite element approxima- tions of eigenvalue/eigenvector problems. Starting from our earlier work on h adaptive finite element approximations we show a way to obtain reliable and efficient a-posteriori estimates in the hp-setting. At the core of our analysis is the reduction of the problem on the analysis of the associated boundary value problem. We start from the analysis of Wohlmuth and Melenk and combine this with our a-posteriori estimation framework to obtain eigenvalue/eigenvector approximation bounds.

Biomimetic Broadband Antireflection Gratings On Solar-Grade Multicrystalline Silicon Wafers, Blayne M. Phillips, Peng Jiang, Bin Jiang

Mathematics and Statistics Faculty Publications and Presentations

The authors report a simple and scalable bottom-up technique for fabricating broadband antireflection gratings on solar-grade multicrystalline silicon (mc-Si) wafers. A Langmuir-Blodgett process is developed to assemble close-packed silica microspheres on rough mc-Si substrates. Subwavelength moth-eye pillars can then be patterned on mc-Si by using the silica microspheres as structural template. Hemispherical reflectance measurements show that the resulting mc-Si gratings exhibit near zero reflection for a wide range of wavelengths. Both experimental results and theoretical prediction using a rigorous coupled-wave analysis model show that close-packed moth-eye arrays exhibit better antireflection performance than non-close-packed arrays due to a smoother refractive index …

Determination Of The Electric Field Intensity And Space Charge Density Versus Height Prior To Triggered Lightning, Christopher J. Biagi, Martin A. Uman, Jay Gopalakrishnan, J. D. Hill, Vladimir A. Rakov, T. Ngin, Douglas M. Jordan

Mathematics and Statistics Faculty Publications and Presentations

We infer the vertical profiles of space charge density and electric field intensity above ground by comparing modeling and measurements of the ground-level electric field changes caused by elevating grounded lightning-triggering wires. The ground-level electric fields at distances of 60 m and 350 m were measured during six wire launches that resulted in triggered lightning. The wires were launched when ground-level electric fields ranged from 3.2 to 7.6 kV m⁻¹ and the triggering heights ranged from 123 to 304 m. From wire launch time to lightning initiation time, the ground-level electric field reduction at 60 m ranged from 2.2 …

Integer Optimization And Computational Algebraic Topology, Bala Krishnamoorthy

Systems Science Friday Noon Seminar Series

We present recently discovered connections between integer optimization, or integer programming (IP), and homology. Under reasonable assumptions, these results lead to efficient solutions of several otherwise hard-to-solve problems from computational topology and geometric analysis. The main result equates the total unimodularity of the boundary matrix of a simplicial complex to an algebraic topological condition on the complex (absence of relative torsion), which is often satisfied in real-life applications . When the boundary matrix is totally unimodular, the problem of finding the shortest chain homologous under Z (ring of integers) to a given chain, which is inherently an integer program, can …

Modal Logic And Its Applications, Explained Using Puzzles And Examples, Marek Perkowski

Systems Science Friday Noon Seminar Series

The talk introduces Modal Logic as an extension of classical propositional and First Order Logics. We discuss motivations of Lewis to create modal logic system, axioms and rules of proof. Several examples illustrate deriving theorems from axioms. "Muddy Children" puzzle is used to explain the principles of dealing with uncertainty problems where a temporal lack of response is used as additional information. Other examples include "Narrow Bridge" problem/game which relates to the problem of necessary evil in the world, robot planning and law and robot morality problems, especially related to military robots and use of force by police. Kripke semantics …

Analysis Of The Dpg Method For The Poisson Equation, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We give an error analysis of the recently developed DPG method applied to solve the Poisson equation and a convection-dffusion problem. We prove that the method is quasioptimal. Error estimates in terms of both the mesh size h and the polynomial degree p (for various element shapes) can be derived from our results. Results of extensive numerical experiments are also presented.

A Second Elasticity Element Using The Matrix Bubble, Jay Gopalakrishnan, Johnny Guzmán

Mathematics and Statistics Faculty Publications and Presentations

We presented a family of finite elements that use a polynomial space augmented by certain matrix bubbles in Cockburn et al. (2010) A new elasticity element made for enforcing weak stress symmetry. Math. Comput., 79, 1331–1349 . In this sequel we exhibit a second family of elements that use the same matrix bubble. This second element uses a stress space smaller than the first while maintaining the same space for rotations (which are the Lagrange multipliers corresponding to a weak symmetry constraint). The space of displacements is of one degree less than the first method. The analysis, while similar to …

Commuting Smoothed Projectors In Weighted Norms With An Application To Axisymmetric Maxwell Equations, Jay Gopalakrishnan, Minah Oh

Mathematics and Statistics Faculty Publications and Presentations

We construct finite element projectors that can be applied to functions with low regularity. These projectors are continuous in a weighted norm arising naturally when modeling devices with axial symmetry. They have important commuting diagram properties needed for finite element analysis. As an application, we use the projectors to prove quasioptimal convergence for the edge finite element approximation of the axisymmetric time-harmonic Maxwell equations on nonsmooth domains. Supplementary numerical investigations on convergence deterioration at high wavenumbers and near Maxwell eigenvalues and are also reported.

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 …

Analysis Of Hdg Methods For Stokes Flow, Bernardo Cockburn, Jay Gopalakrishnan, Ngoc Cuong Nguyen, Jaume Peraire, Francisco-Javier Sayas

Mathematics and Statistics Faculty Publications and Presentations

In this paper, we analyze a hybridizable discontinuous Galerkin method for numerically solving the Stokes equations. The method uses polynomials of degree $ k$ for all the components of the approximate solution of the gradient-velocity-pressure formulation. The novelty of the analysis is the use of a new projection tailored to the very structure of the numerical traces of the method. It renders the analysis of the projection of the errors very concise and allows us to see that the projection of the error in the velocity superconverges. As a consequence, we prove that the approximations of the velocity gradient, the …

Balance Systems And The Variational Bicomplex, Serge Preston

Mathematics and Statistics Faculty Publications and Presentations

In this work we show that the systems of balance equations (balance systems) of continuum thermodynamics occupy a natural place in the variational bicomplex formalism. We apply the vertical homotopy decomposition to get a local splitting (in a convenient domain) of a general balance system as the sum of a Lagrangian part and a complemental "pure non-Lagrangian" balance system. In the case when derivatives of the dynamical fields do not enter the constitutive relations of the balance system, the "pure non-Lagrangian" systems coincide with the systems introduced by S. Godunov [Soviet Math. Dokl. 2 (1961), 947–948] and, later, asserted as …

Polynomial Extension Operators. Part Iii, Leszek Demkowicz, Jay Gopalakrishnan, Joachim Schöberl

Mathematics and Statistics Faculty Publications and Presentations

In this concluding part of a series of papers on tetrahedral polynomial extension operators, the existence of a polynomial extension operator in the Sobolev space H(div) is proven constructively. Specifically, on any tetrahedron K, given a function w on the boundary ∂K that is a polynomial on each face, the extension operator applied to w gives a vector function whose components are polynomials of at most the same degree in the tetrahedron. The vector function is an extension in the sense that the trace of its normal component on the boundary ∂K coincides with w. Furthermore, the extension operator is …

A Class Of Discontinuous Petrov–Galerkin Methods. Part Iv: The Optimal Test Norm And Time-Harmonic Wave Propagation In 1d., Jeffrey Zitelli, Leszek Demkowicz, Jay Gopalakrishnan, D. Pardo, V. M. Calo

Mathematics and Statistics Faculty Publications and Presentations

The phase error, or the pollution effect in the finite element solution of wave propagation problems, is a well known phenomenon that must be confronted when solving problems in the high-frequency range. This paper presents a new method with no phase errors for one-dimensional (1D) time-harmonic wave propagation problems using new ideas that hold promise for the multidimensional case. The method is constructed within the framework of the discontinuous Petrov–Galerkin (DPG) method with optimal test functions. We have previously shown that such methods select solutions that are the best possible approximations in an energy norm dual to any selected test …

Supplementary Balance Laws And The Entropy Principle, Serge Preston

Mathematics and Statistics Faculty Publications and Presentations

In this work we study the mathematical aspects of the development in the continuum thermodynamics known as the “Entropy Principle”. It started with the pioneering works of B.Coleman, W.Noll and I. Muller in 60th of XX cent. and got its further development mostly in the works of G. Boillat, I-Shis Liu and T.Ruggeri. “Entropy Principle” combines in itself the structural requirement on the form of balance laws of the thermodynamical system (denote such system (C)) and on the entropy balance law with the convexity condition of the entropy density. First of these requirements has pure mathematical form defining so called …

Reconstructability Analysis Of Elementary Cellular Automata, Martin Zwick, Hui Shi

Systems Science Friday Noon Seminar Series

Reconstructability analysis is a method to determine whether a multivariate relation, defined set- or information-theoretically, is decomposable with or without loss (reduction in constraint) into lower ordinality relations. Set-theoretic reconstructability analysis (SRA) is used to characterize the mappings of elementary cellular automata. The degree of lossless decomposition possible for each mapping is more effective than the λ parameter (Walker & Ashby, Langton) as a predictor of chaotic dynamics.

Complete SRA yields not only the simplest lossless structure but also a vector of losses of all decomposed structures, indexed by parameter, τ. This vector subsumes λ, Wuensche’s Z parameter, and Walker …

Stochastic Orders In Heterogeneous Samples With Applications, Maochao Xu

Dissertations and Theses

The statistics literature has mostly focused on the case when the data available is in the form of a random sample. In many cases, the observations are not identically distributed. Such samples are called heterogeneous samples. The study of heterogeneous samples is of great interest in many areas, such as statistics, econometrics, reliability engineering, operation research and risk analysis.

Stochastic orders between probability distributions is a widely studied concept. There are several kinds of stochastic orders that are used to compare different aspects of probability distributions like location, variability, skewness, dependence, etc.

In this dissertation, most of the work is …

Phase-Linking And The Perceived Motion During Off-Vertical Axis Rotation, Jan Holly, Scott Wood, Gin Mccollum

Mathematics and Statistics Faculty Publications and Presentations

Human off-vertical axis rotation (OVAR) in the dark typically produces perceived motion about a cone, the amplitude of which changes as a function of frequency. This perception is commonly attributed to the fact that both the OVAR and the conical motion have a gravity vector that rotates about the subject. Little-known, however, is that this rotating-gravity explanation for perceived conical motion is inconsistent with basic observations about self-motion perception: (a) that the perceived vertical moves toward alignment with the gravito-inertial acceleration (GIA) and (b) that perceived translation arises from perceived linear acceleration, as derived from the portion of the GIA …

On Residual Lifetimes Of K-Out-Of-N Systems With Nonidentical Components, Subhash C. Kochar, Maochao Xu

Mathematics and Statistics Faculty Publications and Presentations

In this article, mixture representations of survival functions of residual lifetimes of k-out-of-n systems are obtained when the components are independent but not necessarily identically distributed. Then we stochastically compare the residual lifetimes of k-out-of-n systems in one- and two-sample problems. In particular, the results extend some results in Li and Zhao [14], Khaledi and Shaked [13], Sadegh [17], Gurler and Bairamov [7] and Navarro, Balakrishnan, and Samaniego [16]. Applications in the proportional hazard rates model are presented as well.

A Class Of Discontinuous Petrov–Galerkin Methods. Ii. Optimal Test Functions, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We lay out a program for constructing discontinuous Petrov–Galerkin (DPG) schemes having test function spaces that are automatically computable to guarantee stability. Given a trial space, a DPG discretization using its optimal test space counterpart inherits stability from the well posedness of the undiscretized problem. Although the question of stable test space choice had attracted the attention of many previous authors, the novelty in our approach lies in the fact we identify a discontinuous Galerkin (DG) framework wherein test functions, arbitrarily close to the optimal ones, can be locally computed. The idea is presented abstractly and its feasibility illustrated through …

A Projection-Based Error Analysis Of Hdg Methods, Jay Gopalakrishnan, Bernardo Cockburn, Francisco-Javier Sayas

Mathematics and Statistics Faculty Publications and Presentations

We introduce a new technique for the error analysis of hybridizable discontinuous Galerkin (HDG) methods. The technique relies on the use of a new projection whose design is inspired by the form of the numerical traces of the methods. This renders the analysis of the projections of the discretization errors simple and concise. By showing that these projections of the errors are bounded in terms of the distance between the solution and its projection, our studies of influence of the stabilization parameter are reduced to local analyses of approximation by the projection. We illustrate the technique on a specific HDG …

A Class Of Discontinuous Petrov–Galerkin Methods. Part I: The Transport Equation, Leszek Demkowicz, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

Considering a simple model transport problem, we present a new finite element method. While the new method fits in the class of discontinuous Galerkin (DG) methods, it differs from standard DG and streamline diffusion methods, in that it uses a space of discontinuous trial functions tailored for stability. The new method, unlike the older approaches, yields optimal estimates for the primal variable in both the element size h and polynomial degree p, and outperforms the standard upwind DG method.

On The Integrability Of Orthogonal Distributions In Poisson Manifolds, Daniel Fish, Serge Preston

Mathematics and Statistics Faculty Publications and Presentations

We study conditions for the integrability of the distribution defined on a regular Poisson manifold as the orthogonal complement (with respect to some (pseudo)-Riemannian metric) to the tangent spaces of the leaves of a symplectic foliation. Examples of integrability and non-integrability of this distribution are provided.

Hybridization And Postprocessing Techniques For Mixed Eigenfunctions, Bernardo Cockburn, Jay Gopalakrishnan, F. Li, Ngoc Cuong Nguyen, Jaume Peraire

Mathematics and Statistics Faculty Publications and Presentations

We introduce hybridization and postprocessing techniques for the Raviart–Thomas approximation of second-order elliptic eigenvalue problems. Hybridization reduces the Raviart–Thomas approximation to a condensed eigenproblem. The condensed eigenproblem is nonlinear, but smaller than the original mixed approximation. We derive multiple iterative algorithms for solving the condensed eigenproblem and examine their interrelationships and convergence rates. An element-by-element postprocessing technique to improve accuracy of computed eigenfunctions is also presented. We prove that a projection of the error in the eigenspace approximation by the mixed method (of any order) superconverges and that the postprocessed eigenfunction approximations converge faster for smooth eigenfunctions. Numerical experiments using …

Multigrid In A Weighted Space Arising From Axisymmetric Electromagnetics, Dylan M. Copeland, Jay Gopalakrishnan, Minah Oh

Mathematics and Statistics Faculty Publications and Presentations

Consider the space of two-dimensional vector functions whose components and curl are square integrable with respect to the degenerate weight given by the radial variable. This space arises naturally when modeling electromagnetic problems under axial symmetry and performing a dimension reduction via cylindrical coordinates. We prove that if the original three-dimensional domain is convex then the multigrid Vcycle applied to the inner product in this space converges, provided certain modern smoothers are used. For the convergence analysis, we first prove several intermediate results, e.g., the approximation properties of a commuting projector in weighted norms, and a superconvergence estimate for a …