Applied Mathematics Commons^™

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.

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 …

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 …

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.

Effect Of Network Structure On The Stability Margin Of Large Vehicle With Distributed Control, He Hao, Prabir Barooah, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We study the problem of distributed control of a large network of double-integrator agents to maintain a rigid formation. A few lead vehicles are given information on the desired trajectory of the formation; while every other vehicle uses linear controller which only depends on relative position and velocity from a few other vehicles, which are called its neighbors. A predetermined information graph defines the neighbor relationships. We limit our attention to information graphs that are D-dimensional lattices, and examine the stability margin of the closed loop, which is measured by the real part of the least stable eigenvalue of …

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 …

Symmetries Of The Central Vestibular System: Forming Movements For Gravity And A Three-Dimensional World, Gin Mccollum, Douglas Hanes

Mathematics and Statistics Faculty Publications and Presentations

Intrinsic dynamics of the central vestibular system (CVS) appear to be at least partly determined by the symmetries of its connections. The CVS contributes to whole-body functions such as upright balance and maintenance of gaze direction. These functions coordinate disparate senses (visual, inertial, somatosensory, auditory) and body movements (leg, trunk, head/neck, eye). They are also unified by geometric conditions. Symmetry groups have been found to structure experimentally-recorded pathways of the central vestibular system. When related to geometric conditions in three-dimensional physical space, these symmetry groups make sense as a logical foundation for sensorimotor coordination.

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.

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 …

Symmetry And Stability Of Homogeneous Flocks, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

The study of the movement of flocks, whether biological or technological, is motivated by the desire to understand the capability of coherent motion of a large number of agents that only receive very limited information. In a biological flock a large group of animals seek their course while moving in a more or less fixed formation. It seems reasonable that the immediate course is determined by leaders at the boundary of the flock. The others follow: what is their algorithm? The most popular technological application consists of cars on a one-lane road. The light turns green and the lead car …

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 …

Finite Sample Properties Of Minimum Kolmogorov-Smirnov Estimator And Maximum Likelihood Estimator For Right-Censored Data, Jerzy Wieczorek

Dissertations and Theses

MKSFitter computes minimum Kolmogorov-Smirnov estimators (MKSEs) for several different continuous univariate distributions, using an evolutionary optimization algorithm, and recommends the distribution and parameter estimates that best minimize the Kolmogorov-Smirnov (K-S) test statistic. We modify this tool by extending it to use the Kaplan-Meier estimate of the cumulative distribution function (CDF) for right-censored data. Using simulated data from the most commonly-used survival distributions, we demonstrate the tool's inability to consistently select the correct distribution type with right-censored data, even for large sample sizes and low censoring rates. We also compare this tool's estimates with the right-censored maximum likelihood estimator (MLE). While …

Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems, Bernardo Cockburn, Jay Gopalakrishnan, Raytcho Lazarov

Mathematics and Statistics Faculty Publications and Presentations

We introduce a unifying framework for hybridization of finite element methods for second order elliptic problems. The methods fitting in the framework are a general class of mixed-dual finite element methods including hybridized mixed, continu- ous Galerkin, non-conforming and a new, wide class of hybridizable discontinuous Galerkin methods. The distinctive feature of the methods in this framework is that the only globally coupled degrees of freedom are those of an approximation of the solution defined only on the boundaries of the elements. Since the associated matrix is sparse, symmetric and positive definite, these methods can be efficiently implemented. Moreover, the …

Stability Of Linear Flocks On A Ring Road, J. J. P. Veerman, Carlos Martins Da Fonseca

Mathematics and Statistics Faculty Publications and Presentations

We discuss some stability problems when each agent of a linear flock on the line interacts with its two nearest neighbors (one on either side).

A Single Particle Impact Model For Motion In Avalanches, J. J. P. Veerman, Dacian Daescu, M. J. Romero-Vallés, P. J. Torres

Mathematics and Statistics Faculty Publications and Presentations

We describe the global behavior of the dynamics of a particle bouncing down an inclined staircase. For small inclinations all orbits eventually stop (independent of the initial condition). For large enough inclinations all orbits end up accelerating indefinitely (also independent of the initial conditions). There is an interval of inclinations of positive length between these two. In that interval the behavior of an orbit depends on its initial condition. In addition to stopping and accelerating orbits, there are also orbits with speeds bounded away from both zero and infinity. A second hallmark of the dynamics is that the orbits going …

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

Mathematics and Statistics Faculty Publications and Presentations

Consider the tangential trace of a vector polynomial on the surface of a tetrahedron. We construct an extension operator that extends such a trace function into a polynomial on the tetrahedron. This operator can be continuously extended to the trace space of H(curl ). Furthermore, it satisfies a commutativity property with an extension operator we constructed in Part I of this series. Such extensions are a fundamental ingredient of high order finite element analysis.

Automated Traffic And The Finite Size Resonance, J. J. P. Veerman, Borko D. Stošić, F. M. Tangerman

Mathematics and Statistics Faculty Publications and Presentations

We investigate in detail what one might call the canonical (automated) traffic problem: A long string of N+1 cars (numbered from 0 to N) moves along a one-lane road “in formation” at a constant velocity and with a unit distance between successive cars. Each car monitors the relative velocity and position of only its neighboring cars. This information is then fed back to its own engine which decelerates (brakes) or accelerates according to the information it receives. The question is: What happens when due to an external influence—a traffic light turning green—the ‘zero’th’ car (the “leader”) accelerates?

As …

A New Elasticity Element Made For Enforcing Weak Stress Symmetry, Bernardo Cockburn, Jay Gopalakrishnan, Johnny Guzmán

Mathematics and Statistics Faculty Publications and Presentations

We introduce a new mixed method for linear elasticity. The novelty is a simplicial element for the approximate stress. For every positive integer k, the row-wise divergence of the element space spans the set of polynomials of total degree k. The degrees of freedom are suited to achieve continuity of the normal stresses. What makes the element distinctive is that its dimension is the smallest required for enforcing a weak symmetry condition on the approximate stress. This is achieved using certain "bubble matrices", which are special divergence-free matrix-valued polynomials. We prove that the approximation error is of order k + …

The Derivation Of Hybridizable Discontinuous Galerkin Methods For Stokes Flow, Bernardo Cockburn, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

In this paper, we introduce a new class of discontinuous Galerkin methods for the Stokes equations. The main feature of these methods is that they can be implemented in an efficient way through a hybridization procedure which reduces the globally coupled unknowns to certain approximations on the element boundaries. We present four ways of hybridizing the methods, which differ by the choice of the globally coupled unknowns. Classical methods for the Stokes equations can be thought of as limiting cases of these new methods.

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

Mathematics and Statistics Faculty Publications and Presentations

In this series of papers, we construct operators that extend certain given functions on the boundary of a tetrahedron into the interior of the tetrahedron, with continuity properties in appropriate Sobolev norms. These extensions are novel in that they have certain polynomial preservation properties important in the analysis of high order finite elements. This part of the series is devoted to introducing our new technique for constructing the extensions, and its application to the case of polynomial extensions from H^½(∂K) into H¹(K), for any tetrahedron K.

Asymptotic And Numerical Techniques For Resonances Of Thin Photonic Structures, Jay Gopalakrishnan, Shari Moskow, Fadil Santosa

Mathematics and Statistics Faculty Publications and Presentations

We consider the problem of calculating resonance frequencies and radiative losses of an optical resonator. The optical resonator is in the form of a thin membrane with variable dielectric properties. This work provides two very different approaches for doing such calculations. The first is an asymptotic method which exploits the small thickness and high index of the membrane. We derive a limiting resonance problem as the thickness goes to zero, and for the case of a simple resonance, find a first order correction. The limiting problem and the correction are in one less space dimension, which can make the approach …

A Mixed Method For Axisymmetric Div-Curl Systems, Dylan M. Copeland, Jay Gopalakrishnan, Joseph E. Pasciak

Mathematics and Statistics Faculty Publications and Presentations

We present a mixed method for a three-dimensional axisymmetric div-curl system reduced to a two-dimensional computational domain via cylindrical coordinates. We show that when the meridian axisymmetric Maxwell problem is approximated by a mixed method using the lowest order Nédélec elements (for the vector variable) and linear elements (for the Lagrange multiplier), one obtains optimal error estimates in certain weighted Sobolev norms. The main ingredient of the analysis is a sequence of projectors in the weighted norms satisfying some commutativity properties.

On The Spectra Of Certain Directed Paths, Carlos Martins Da Fonseca, J. J. P. Veerman

Mathematics and Statistics Faculty Publications and Presentations

We describe the eigenpairs of special kinds of tridiagonal matrices related to problems on traffic on a one-lane road. Some numerical examples are provided.

Multigrid Convergence For Second Order Elliptic Problems With Smooth Complex Coefficients, Jay Gopalakrishnan, Joseph E. Pasciak

Mathematics and Statistics Faculty Publications and Presentations

The finite element method when applied to a second order partial differential equation in divergence form can generate operators that are neither Hermitian nor definite when the coefficient function is complex valued. For such problems, under a uniqueness assumption, we prove the continuous dependence of the exact solution and its finite element approximations on data provided that the coefficients are smooth and uniformly bounded away from zero. Then we show that a multigrid algorithm converges once the coarse mesh size is smaller than some fixed number, providing an efficient solver for computing discrete approximations. Numerical experiments, while confirming the theory, …

Locally Conservative Fluxes For The Continuous Galerkin Method, Bernardo Cockburn, Jay Gopalakrishnan, Haiying Wang

Mathematics and Statistics Faculty Publications and Presentations

The standard continuous Galerkin (CG) finite element method for second order elliptic problems suffers from its inability to provide conservative flux approximations, a much needed quantity in many applications. We show how to overcome this shortcoming by using a two step postprocessing. The first step is the computation of a numerical flux trace defined on element inter- faces and is motivated by the structure of the numerical traces of discontinuous Galerkin methods. This computation is non-local in that it requires the solution of a symmetric positive definite system, but the system is well conditioned independently of mesh size, so it …

The Convergence Of V-Cycle Multigrid Algorithms For Axisymmetric Laplace And Maxwell Equations, Jay Gopalakrishnan, Joseph E. Pasciak

Mathematics and Statistics Faculty Publications and Presentations

We investigate some simple finite element discretizations for the axisymmetric Laplace equation and the azimuthal component of the axisymmetric Maxwell equations as well as multigrid algorithms for these discretizations. Our analysis is targeted at simple model problems and our main result is that the standard V-cycle with point smoothing converges at a rate independent of the number of unknowns. This is contrary to suggestions in the existing literature that line relaxations and semicoarsening are needed in multigrid algorithms to overcome difficulties caused by the singularities in the axisymmetric Maxwell problems. Our multigrid analysis proceeds by applying the well known regularity …

Error Analysis Of Variable Degree Mixed Methods For Elliptic Problems Via Hybridization, Bernardo Cockburn, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

A new approach to error analysis of hybridized mixed methods is proposed and applied to study a new hybridized variable degree Raviart-Thomas method for second order elliptic problems. The approach gives error estimates for the Lagrange multipliers without using error estimates for the other variables. Error estimates for the primal and flux variables then follow from those for the Lagrange multipliers. In contrast, traditional error analyses obtain error estimates for the flux and primal variables first and then use it to get error estimates for the Lagrange multipliers. The new approach not only gives new error estimates for the new …

Incompressible Finite Elements Via Hybridization. Part Ii: The Stokes System In Three Space Dimensions, Bernardo Cockburn, Jay Gopalakrishnan

Mathematics and Statistics Faculty Publications and Presentations

We introduce a method that gives exactly incompressible velocity approximations to Stokes ow in three space dimensions. The method is designed by extending the ideas in Part I (http://archives.pdx.edu/ds/psu/10914) of this series, where the Stokes system in two space dimensions was considered. Thus we hybridize a vorticity-velocity formulation to obtain a new mixed method coupling approximations of tangential velocity and pressure on mesh faces. Once this relatively small tangential velocity-pressure system is solved, it is possible to recover a globally divergence-free numerical approximation of the fluid velocity, an approximation of the vorticity whose tangential component is continuous across …