@book{feec-cbms, author = {Arnold, Douglas N.}, title = {Finite {E}lement {E}xterior {C}alculus}, series = {CBMS-NSF Regional Conference Series in Applied Mathematics}, volume = {93}, publisher = {Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA}, year = {2018}, pages = {xii+120}, isbn = {978-1-611975-53-6}, keywords = {finite element method, mixed method, PDE, finite element exterior calculus, de Rham complex, Hodge theory, thermodynamics, fluid flow, solid deformation, elasticity, electricity and magnetism, numerical methods}, abstract = {Computational methods to approximate the solution of differential equations play a crucial role in science, engineering, mathematics, and technology. The key processes that govern the physical world--wave propagation, thermodynamics, fluid flow, solid deformation, electricity and magnetism, quantum mechanics, general relativity, and many more--are described by differential equations. We depend on numerical methods for the ability to simulate, explore, predict, and control systems involving these processes. The finite element exterior calculus, or FEEC, is a powerful new theoretical approach to the design and understanding of numerical methods to solve partial differential equations (PDEs). The methods derived with FEEC preserve crucial geometric and topological structures underlying the equations and are among the most successful examples of structure-preserving methods in numerical PDEs. This volume aims to help numerical analysts master the fundamentals of FEEC, including the geometrical and functional analysis preliminaries, quickly and in one place. It is also accessible to mathematicians and students of mathematics from areas other than numerical analysis who are interested in understanding how techniques from geometry and topology play a role in numerical PDEs.} }

@proceedings{compatible, title = {Compatible Spatial Discretizations}, booktitle = {Compatible Spatial Discretizations}, editor = {Arnold, D. and Bochev, P. and Lehoucq, R. and Nicolaides, R. and Shashkov, M.}, series = {IMA Vol. Math. Appl.}, volume = {142}, publisher = {Springer, Berlin}, year = {2006}, abstract = {This volume constains papers presented at a highly successful IMA Hot Topics Workshop: Compatible Spatial Discretizations for Partial Differential Equations, held on May 11-15, 2004.}, }

@unpublished{curved_hhj, author = {Arnold, Douglas N. and Walker, Shawn W.}, title = {The {Hellan-Herrmann-Johnson} method with curved elements}, year = {2019}, note = {Submitted. arXiv preprint 1909.09687}, abstract = {We study the finite element approximation of the Kirchhoff plate equation on domains with curved boundaries using the Hellan-Herrmann-Johnson (HHJ) method. We prove optimal convergence on domains with piecewise C^{k+1} boundary for k \geq 1 when using a parametric (curved) HHJ space. Computational results are given that demonstrate optimal convergence and how convergence degrades when curved triangles of insufficient polynomial degree are used. Moreover, we show that the lowest order HHJ method on a polygonal approximation of the disk does not succumb to the classic Babuska paradox, highlighting the geometrically non-conforming aspect of the HHJ method. }, keywords = {Kirchhoff plate, simply-supported, parametric finite elements, mesh-dependent norms, geometric consistency error, Babuska paradox}, subjclass = {65N30, 35J40, 35Q72}, pdf = {http://umn.edu/~arnold/papers/curved_hhj.pdf} }

@article{agmonweyl, author = {Arnold, Douglas N. and David, Guy and Filoche, Marcel and Jerison, David and Mayboroda, Svitlana}, title = {Localization of eigenfunctions via an effective potential}, journal = {Communications in Partial Differential Equations}, volume = {44}, number = {11}, pages = {1186-1216}, year = {2019}, skipnote = {To appear. arXiv preprint 1712.02419}, abstract = {We consider the Neumann boundary value problem for an elliptic operator L = -div A grad + V on a bounded bi-Lipschitz domain in R^n. More generally, we will treat closed manifolds and manifolds with boundary. The eigenfunctions of L are often localized, as a result of disorder of the potential V, the matrix of coefficients A, irregularities of the boundary, or all of the above. In earlier work, two of us introduced the function u solving Lu=1, and showed numerically that it strongly reflects this localization. In this paper, we deepen the connection between the eigenfunctions and this landscape function u by proving that its reciprocal 1/u acts as an effective potential. The effective potential governs the exponential decay of the eigenfunctions of the system and delivers information on the distribution of eigenvalues near the bottom of the spectrum. }, keywords = {Agmon distance, Schr\"odinger equation, spectrum, the landscape of localization}, subjclass = {35P,35Q40}, pdf = {http://umn.edu/~arnold/papers/agmonweyl.pdf}, doi = {10.1080/03605302.2019.1626420} }

@article{specpred, author = {Arnold, Douglas N. and David, Guy and Filoche, Marcel and Jerison, David and Mayboroda, Svitlana}, title = {Computing spectra without solving eigenvalue problems}, journal = {SIAM Journal on Scientific Computing}, volume = {41}, number = {1}, year = {2019}, pages = {B69-B92}, skipnote = {arXiv preprint 1711.04888}, abstract = { The approximation of the eigenvalues and eigenfunctions of an elliptic operator is a key computational task in many areas of applied mathematics and computational physics. An important case, especially in quantum physics, is the computation of the spectrum of a Schr\"odinger operator with a disordered potential. Unlike plane waves or Bloch waves that arise as Schr\"odinger eigenfunctions for periodic and other ordered potentials, for many forms of disordered potentials the eigenfunctions remain essentially localized in a very small subset of the initial domain. A celebrated example is Anderson localization, for which, in a continuous version, the potential is a piecewise constant function on a uniform grid whose values are sampled independently from a uniform random distribution. We present here a new method for approximating the eigenvalues and the subregions which support such localized eigenfunctions. This approach is based on the recent theoretical tools of the localization landscape and effective potential. The approach is deterministic in the sense that the approximations are calculated based on the examination of a particular realization of a random potential, and predict quantities that depend sensitively on the particular realization, rather than furnishing statistical or probabilistic results about the spectrum associated to a family of potentials with a certain distribution. These methods, which have only been partially justified theoretically, enable the calculation of the locations and shapes of the approximate supports of the eigenfunctions, the approximate values of many of the eigenvalues, and of the eigenvalue counting function and density of states, all at the cost of solving a single source problem for the same elliptic operator. We study the effectiveness and limitations of the approach through extensive computations in one and two dimensions, using a variety of piecewise constant potentials with values sampled from various different correlated or uncorrelated random distributions.}, keywords = {localization, spectrum, eigenvalue, Schr\"odinger operator}, subjclass = {65N25, 81-08, 82B44}, pdf = {http://umn.edu/~arnold/papers/specpred.pdf}, doi = {10.1137/17M1156721} }

@article{feecparabolic, author = {Arnold, Douglas N. and Chen, Hongtao}, title = {Finite element exterior calculus for parabolic problems}, journal = {ESAIM Math. Model. Numer. Anal.}, fjournal = {ESAIM. Mathematical Modelling and Numerical Analysis}, volume = {51}, number = {1}, year = {2017}, pages = {17-34}, skipnote = {arXiv preprint 1209.1142}, abstract = { In this paper, we consider the extension of the finite element exterior calculus from elliptic problems, in which the Hodge Laplacian is an appropriate model problem, to parabolic problems, for which we take the Hodge heat equation as our model problem. The numerical method we study is a Galerkin method based on a mixed variational formulation and using as subspaces the same spaces of finite element differential forms which are used for elliptic problems. We analyze both the semidiscrete and a fully-discrete numerical scheme. }, keywords = {finite element exterior calculus, mixed finite element method, parabolic equation, Hodge heat equation}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/feecparabolic.pdf}, doi = {10.1051/m2an/2016013} }

@article{feeclot, author = {Arnold, Douglas N. and Li, Lizao}, title = {Finite element exterior calculus with lower-order terms}, journal = {Math. Comput.}, fjournal = {Mathematics of Computation}, volume = {86}, year = {2017}, pages = {2193-2212}, skipnote = {arXiv preprint 1509.06463}, abstract = { The scalar and vector Laplacians are basic operators in physics and engineering. In applications, they show up frequently perturbed by lower-order terms. The effect of such perturbations on mixed finite element methods in the scalar case is well-understood, but that in the vector case is not. In this paper, we first show that surprisingly for certain elements there is degradation of the convergence rates with certain lower-order terms even when both the solution and the data are smooth. We then give a systematic analysis of lower-order terms in mixed methods by extending the Finite Element Exterior Calculus (FEEC) framework, which contains the scalar, vector Laplacian, and many other elliptic operators as special cases. We prove that stable mixed discretization remains stable with lower-order terms for sufficiently fine discretization. Moreover, we derive sharp improved error estimates for each individual variable. In particular, this yields new results for the vector Laplacian problem which are useful in applications such as electromagnetism and acoustics modeling. Further our results imply many previous results for the scalar problem and thus unifies them all under the FEEC framework. }, keywords = {finite element exterior calculus, lower order terms}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/feeclot.pdf}, doi = {10.1090/mcom/3158} }

@article{localization, author = {Arnold, Douglas N. and David, Guy and Jerison, David and Mayboroda, Svitlana and Filoche, Marcel}, title = {Effective confining potential of quantum states in disordered media}, journal = {Physical Review Letters}, fjournal = {Physical Review Letters}, volume = {116}, number = {5}, year = {2016}, abstract = { The amplitude of localized quantum states in random or disordered media may exhibit long range exponential decay. We present here a theory that unveils the existence of an effective potential which finely governs the confinement of these states. In this picture, the boundaries of the localization subregions for low energy eigenfunctions correspond to the barriers of this effective potential, and the long range exponential decay characteristic of Anderson localization is explained as the consequence of multiple tunneling in the dense network of barriers created by this effective potential. Finally, we show that the Weyl’s formula based on this potential turns out to be a remarkable approximation of the density of states for a large variety of one-dimensional systems, periodic or random. }, keywords = {Anderson localization, effective potential, Weyl's law, disordered media}, subjclass = {PACS: 71.23.An, 72.15.Rn, 03.65.Ge}, pdf = {http://umn.edu/~arnold/papers/localization.pdf}, doi = {10.1103/PhysRevLett.116.056602} }

@article{quadfeec, author = {Arnold, Douglas N. and Boffi, Daniele and Bonizzoni, Francesca}, title = {Finite element differential forms on curvilinear cubic meshes and their approximation properties}, year = {2015}, volume = {129}, issue = {1}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, pages = {1--20}, abstract = { We study the approximation properties of a wide class of finite element differential forms on curvilinear cubic meshes in n dimensions. Specifically, we consider meshes in which each element is the image of a cubical reference element under a diffeomorphism, and finite element spaces in which the shape functions and degrees of freedom are obtained from the reference element by pullback of differential forms. In the case where the diffeomorphisms from the reference element are all affine, i.e., mesh consists of parallelotopes, it is standard that the rate of convergence in L2 exceeds by one the degree of the largest full polynomial space contained in the reference space of shape functions. When the diffeomorphism is multilinear, the rate of convergence for the same space of reference shape function may degrade severely, the more so when the form degree is larger. The main result of the paper gives a sufficient condition on the reference shape functions to obtain a given rate of convergence.}, keywords = {65N30}, subjclass = {mixed finite elements, finite element differential forms, finite element exterior calculus, cubical meshes, cubes}, pdf = {http://umn.edu/~arnold/papers/quadfeec.pdf}, skipnote = {arXiv preprint 1212.6559 }, doi = {10.1007/s00211-014-0631-3} }

@article{quadelas, author = {Arnold, Douglas N. and Awanou, Gerard and Qiu, Weifeng}, title = {Mixed finite elements for elasticity on quadrilateral meshes}, journal = {Adv. Comput. Math.}, fjournal = {Advances in Computational Mathematics}, year = {2015}, volume = {41}, issue = {3}, pages = {553--572}, abstract = { We present stable mixed finite elements for planar linear elasticity on general quadrilateral meshes. The symmetry of the stress tensor is imposed weakly and so there are three primary variables, the stress tensor, the displacement vector field, and the scalar rotation. We develop and analyze a stable family of methods, indexed by an integer r > 1 and with rate of convergence in the L2 norm of order r for all the variables. The methods use Raviart-Thomas elements for the stress, piecewise tensor product polynomials for the displacement, and piecewise polynomials for the rotation. We also present a simple first order element, not belonging to this family. It uses the lowest order BDM elements for the stress, and piecewise constants for the displacement and rotation, and achieves first order convergence for all three variables. }, keywords = {mixed finite element method, linear elasticty, quadrilateral elements}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/quadelas.pdf}, skipnote = {Published online September 9, 2014. arXiv preprint 1306.6821}, doi = {10.1007/s10444-014-9376-x} }

@incollection{stability, author = {Arnold, Douglas N.}, title = {Stability, consistency, and convergence of numerical discretizations}, booktitle = {Encyclopedia of Applied and Computational Mathematics}, editor = {Engquist, Bj\"orn}, publisher = {Springer}, year = {2015}, pages = {1358-1364}, abstract = {This expository article discusses the meaning of consistency, stability, and convergence of numerical discretizations of differential equations. It provides a general framework for quantifying them in which the fundamental theorem that consistency and stability imply convergence can be rigorously stated and proved. These concepts are illustrated with examples coming from both finite difference methods and finite element methods.}, keywords = {consistency, stability, convergence, discretization}, pdf = {http://umn.edu/~arnold/papers/stability.pdf}, doi = {10.1007/978-3-540-70529-1} }

@incollection{golf-flight, author = {Arnold, Douglas N.}, title = {The flight of a golf ball}, booktitle = {The Princeton Companion to Applied Mathematics}, editor = {Nicholas J. Higham and Mark R. Dennis and Paul Glendinning and Paul A. Martin and Fadil Santosa and Jared Tanner}, year = {2015}, publisher = {Princeton University Press}, address = {Princeton, NJ, USA}, pages = {746-749}, abstract = {This expository article describes the modeling of the flight of a golf ball, with particular attention to the resolution of the drag crisis, the role of the dimples, and the optimization of the ball.}, keywords = {golf, drag crisis}, pdf = {http://umn.edu/~arnold/papers/golf-flight.pdf} }

@article{periodic-table, author = {Arnold, Douglas N. and Logg, Anders}, title = {Periodic table of the finite elements}, journal = {SIAM News}, fjournal = {SIAM News}, volume = {47}, year = {2014}, number = {9}, abstract = {This article presents and explains the periodic table of finite elements.}, keywords = {finite element, exterior calculus}, pdf = {http://umn.edu/~arnold/papers/periodic-table.pdf} }

@article{mixed-elastodynamics, author = {Arnold, Douglas N. and Lee, Jeonghun J.}, title = {Mixed methods for elastodynamics with weak symmetry}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {52}, year = {2014}, number = {6}, pages = {2743--2769}, abstract = { We analyze the application to elastodynamic problems of mixed finite element methods for elasticity with weak symmetry. Our approach leads to a semidiscrete method which consists of a system of ordinary differential equations without algebraic constraints. Our error analysis, which is based on a new elliptic projection operator, applies to several mixed finite element spaces developed for elastostatics. The error estimates we obtain are robust for nearly incompressible materials.}, keywords = {mixed finite element, elastodynamics, weak symmetry}, subjclass = {65N30, 74H15, 74S05}, pdf = {http://umn.edu/~arnold/papers/mixed-elastodynamics.pdf}, skipnote = {arXiv preprint 1312.6319}, doi = {10.1137/13095032X} }

@article{cubicderham, author = {Arnold, Douglas N. and Awanou, Gerard}, title = {Finite element differential forms on cubical meshes}, journal = {Math. Comput.}, fjournal = {Mathematics of Computation}, year = {2014}, volume = {83}, issue = {288}, pages = {1551--1570}, abstract = { We develop a family of finite element spaces of differential forms defined on cubical meshes in any number of dimensions. The family contains elements of all polynomial degrees and all form degrees. In two dimensions, these include the serendipity finite elements and the rectangular BDM elements. In three dimensions they include a recent generalization of the serendipity spaces, and new H(curl) and H(div) finite element spaces. Spaces in the family can be combined to give finite element subcomplexes of the de Rham complex which satisfy the basic hypotheses of the finite element exterior calculus, and hence can be used for stable discretization of a variety of problems. The construction and properties of the spaces are established in a uniform manner using finite element exterior calculus. }, keywords = {65N30}, subjclass = {mixed finite elements, finite element differential forms, finite element exterior calculus, cubical meshes, cubes}, pdf = {http://umn.edu/~arnold/papers/cubicderham.pdf}, skipnote = {arXiv preprint 1204.2595}, doi = {10.1090/S0025-5718-2013-02783-4} }

@article{dstar, author = {Arnold, Douglas N. and Falk, Richard S. and Guzm{\'a}n, Johnny and Tsogtgerel, Gantumur}, title = {On the consistency of the combinatorial codifferential}, year = {2014}, journal = {Trans. Amer. Math. Soc.}, fjournal = {Transactions of the American Mathematical Society}, volume = {366}, issue = {10}, pages = {5487-5502}, abstract = { In 1976, Dodziuk and Patodi employed Whitney forms to define a combinatorial codifferential operator on cochains, and they raised the question whether it is consistent in the sense that for a smooth enough differential form the combinatorial codifferential of the associated cochain converges to the exterior codifferential of the form as the triangulation is refined. In 1991, Smits proved this to be the case for the combinatorial codifferential applied to 1-forms in two dimensions under the additional assumption that the initial triangulation is refined in a completely regular fashion, by dividing each triangle into four similar triangles. In this paper we extend Smits's result to arbitrary dimensions, showing that the combinatorial codifferential on 1-forms is consistent if the triangulations are uniform or piecewise uniform in a certain precise sense. We also show that this restriction on the triangulations is needed, giving a counterexample in which a different regular refinement procedure, namely Whitney's standard subdivision, is used. Further, we show by numerical example that for 2-forms in three dimensions, the combinatorial codifferential is not consistent even for the most regular subdivision process. }, keywords = {consistency, combinatorial codifferential, Whitney form, finite element}, subjclass = {58A10, 65N30; 39A12, 57Q55}, pdf = {http://umn.edu/~arnold/papers/dstar.pdf}, skipnote = {arXiv preprint 1212.4472}, doi = {10.1090/S0002-9947-2014-06134-5} }

@article{nonconfelas3d, author = {Arnold, Douglas N. and Awanou, Gerard and Winther, Ragnar}, title = {Nonconforming tetrahedral mixed finite elements for elasticity}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, year = {2014}, volume = {24}, issue = {4}, pages = {783-796}, abstract = { This paper presents a nonconforming finite element approximation of the space of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for stress approximation together with the full space of piecewise linear vector fields for displacement, this gives a stable mixed finite element method which is shown to be linearly convergent for both the stress and displacement, and which is significantly simpler than any stable conforming mixed finite element method. The method may be viewed as the three-dimensional analogue of a previously developed element in two dimensions. As in that case, a variant of the method is proposed as well, in which the displacement approximation is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the linear convergence is retained. }, keywords = {mixed method; finite element; linear elasticity; nonconforming}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/nonconfelas3d.pdf}, skipnote = {arXiv preprint 1210.6256}, doi = {10.1142/S021820251350067X} }

@incollection{fe-families, author = {Arnold, Douglas N.}, title = {Spaces of finite element differential forms}, booktitle = {Analysis and Numerics of Partial Differential Equations}, editor = {Gianazza, U. and Brezzi, F. and Colli Franzone, P. and Gilardi, G.}, publisher = {Springer}, year = {2013}, pages = {117--140}, abstract = {We discuss the construction of ﬁnite element spaces of differential forms which satisfy the crucial assumptions of the ﬁnite element exterior calculus, namely that they can be assembled into subcomplexes of the de Rham complex which admit commuting projections. We present two families of spaces in the case of simplicial meshes, and two other families in the case of cubical meshes. We make use of the exterior calculus and the Koszul complex to deﬁne and understand the spaces. These tools allow us to treat a wide variety of situations, which are often treated separately, in a uniﬁed fashion.}, keywords = {finite element differential form, finite element exterior calculus}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/fe-families.pdf}, skipnote = {arXiv preprint 1208.2041}, doi = {10.1007/978-88-470-2592-9_9} }

@article{boycott, author = {Arnold, Douglas N. and Cohn, Henry}, title = {Mathematicians take a stand}, journal = {Notices Amer. Math. Soc.}, fjournal = {Notices of the American Mathematical Society}, volume = {59}, number = {6}, year = {2012}, pages = {828--833}, abstract = { We survey the reasons for the ongoing boycott of the publisher Elsevier. We examine Elsevier's pricing and bundling policies, restrictions on dissemination by authors, and lapses in ethics and peer review, and we conclude with thoughts about the future of mathematical publishing. }, pdf = {http://umn.edu/~arnold/papers/boycott.pdf}, skipnote = {arXiv preprint 1204.1351}, doi = {10.1090/noti857} }

@article{vectorlaplacian, author = {Arnold, Douglas N. and Falk, Richard S. and Gopalakrishnan, Jay}, title = {Mixed finite element approximation of the vector {L}aplacian with {D}irichlet boundary conditions}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, volume = {22}, number = {9}, year = {2012}, pages = {26 pages}, abstract = { We consider the ﬁnite element solution of the vector Laplace equation on a domain in two dimensions. For various choices of boundary conditions, it is known that a mixed ﬁnite element method, in which the rotation of the solution is introduced as a second unknown, is advantageous, and appropriate choices of mixed ﬁnite element spaces lead to a stable, optimally convergent discretization. However, the theory that leads to these conclusions does not apply to the case of Dirichlet boundary conditions, in which both components of the solution vanish on the boundary. We show, by computational example, that indeed such mixed ﬁnite elements do not perform optimally in this case, and we analyze the suboptimal convergence that does occur. As we indicate, these results have implications for the solution of the biharmonic equation and of the Stokes equations using a mixed formulation involving the vorticity. }, keywords = {vector Laplacian, Hodge Laplacian, mixed finite elements}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/vectorlaplacian.pdf}, skipnote = {arXiv preprint 1109.3668}, doi = {10.1142/S0218202512500248} }

@incollection{integrityconf, author = {Arnold, Douglas N.}, title = {Challenges and responses in mathematical research publishing}, booktitle = {Promoting Research Integrity in a Global Environment}, editor = {Mayer, Tony and Steneck, Nicholas}, pages = {301--304}, publisher = {World Scientific}, address = {Singapore}, year = {2012}, abstract = { This article discusses some of the challenges the mathematical community faces in maintaining ethical scholarly publishing, and some recent actions that have been made to confront them. }, keywords = {research integrity, publishing}, pdf = {http://umn.edu/~arnold/papers/integrityconf.pdf} }

@article{impact-factors, author = {Arnold, Douglas N. and Fowler, Kristine K.}, title = {Nefarious Numbers}, journal = {Notices Amer. Math. Soc.}, fjournal = {Notices of the American Mathematical Society}, volume = {58}, year = {2011}, number = {3}, pages = {434--437}, note = {Also appeared in Gazette of the Australian Mathematical Society 38(1):9--16, 2011; Newsletter of the European Mathematical Society 80:34--36, 2011; and, in Russian translation, in a volume on bibliometrics, Moscow, 2011.}, abstract = { We investigate the journal impact factor, focusing on the applied mathematics category. We discuss impact factor manipulation and demonstrate that the impact factor gives an inaccurate view of journal quality, which is poorly correlated with expert opinion. }, keywords = {impact factor}, pdf = {http://umn.edu/~arnold/papers/impact-factors.pdf} }

@article{serendipity, author = {Arnold, Douglas N. and Awanou, Gerard}, title = {The serendipity family of finite elements}, year = {2011}, journal = {Found. Comput. Math.}, fjournal = {Foundations of Computational Mathematics. The Journal of the Society for the Foundations of Computational Mathematics}, volume = {11}, pages = {337--344}, abstract = { We give a new, simple, dimension-independent definition of the serendipity finite element family. The shape functions are the span of all monomials which are linear in at least s-r of the variables where s is the degree of the monomial or, equivalently, whose superlinear degree (total degree with respect to variables entering at least quadratically) is at most r. The degrees of freedom are given by moments of degree at most r-2d on each face of dimension d. We establish unisolvence and a geometric decomposition of the space. }, keywords = {serendipity, finite element, unisolvence}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/serendipity.pdf}, doi = {10.1007/s10208-011-9087-3} }

@article{golf, author = {Arnold, Douglas N.}, title = {The science of a drive}, journal = {Notices Amer. Math. Soc.}, fjournal = {Notices of the American Mathematical Society}, volume = {57}, year = {2010}, number = {4}, pages = {498--501}, note = {Also appeared in the collection \emph{Mathematics and Sports}, J.~Gallian, ed., Mathematical Association of America Dolciani Mathematical Expositions \# 43, 2010, 149--156.}, abstract = { This article, prepared in conjunction with Mathematics Awareness Month 2010 with the theme of Mathematics and Sports, reviews several examples of how mathematics eluciates physical phenomena pertaining to the golf drive. Specifically it discusses the double-pendulum model of a golf swing, transfer of energy and momentum in the club head/ball impact, and drag and lift in the flight of the golf ball. }, keywords = {golf, mathematical modeling}, pdf = {http://umn.edu/~arnold/papers/golf.pdf} }

@article{bulletin, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Finite element exterior calculus: from {H}odge theory to numerical stability}, journal = {Bull. Amer. Math. Soc. (N.S.)}, fjournal = {American Mathematical Society. Bulletin. New Series}, volume = {47}, year = {2010}, pages = {281--354}, abstract = { This article reports on the confluence of two streams of research, one emanating from the fields of numerical analysis and scientific computation, the other from topology and geometry. In it we consider the numerical discretization of partial differential equations that are related to differential complexes so that de~Rham cohomology and Hodge theory are key tools for exploring the well-posedness of the continuous problem. The discretization methods we consider are finite element methods, in which a variational or weak formulation of the PDE problem is approximated by restricting the trial subspace to an appropriately constructed piecewise polynomial subspace. After a brief introduction to finite element methods, we develop an abstract Hilbert space framework for analyzing the stability and convergence of such discretizations. In this framework, the differential complex is represented by a complex of Hilbert spaces and stability is obtained by transferring Hodge theoretic structures that ensure well-posedness of the continuous problem from the continuous level to the discrete. We show stable discretization discretization is achieved if the finite element spaces satisfy two hypotheses: they can be arranged into a subcomplex of this Hilbert complex, and there exists a bounded cochain projection from that complex to the subcomplex. In the next part of the paper, we consider the most canonical example of the abstract theory, in which the Hilbert complex is the de~Rham complex of a domain in Euclidean space. We use the Koszul complex to construct two families of finite element differential forms, show that these can be arranged in subcomplexes of the de~Rham complex in numerous ways, and for each construct a bounded cochain projection. The abstract theory therefore applies to give the stability and convergence of finite element approximations of the Hodge Laplacian. Other applications are considered as well, especially the elasticity complex and its application to the equations of elasticity. Background material is included to make the presentation self-contained for a variety of readers. }, keywords = {finite element exterior calculus, exterior calculus, de Rham cohomology, Hodge theory, Hodge Laplacian, mixed finite elements}, subjclass = {65N30, 58A14}, pdf = {http://umn.edu/~arnold/papers/bulletin.pdf}, doi = {10.1090/S0273-0979-10-01278-4} }

@article{integrity, author = {Arnold, Douglas N.}, title = {Integrity under attack: the state of scholarly publishing}, journal = {SIAM News}, fjournal = {SIAM News}, volume = {42}, year = {2009}, number = {10}, pages = {2--3}, abstract = { An editorial on integrity in mathematical scholarly publishing}, keywords = {scholarly publishing, integrity, plagiarism, impact factor}, pdf = {http://umn.edu/~arnold/papers/integrity.pdf}, note = {Spanish translation in \emph{Gac.\ R.\ Soc.\ Mat.\ Esp.}, 13(1):21-25, 2010; Chinese translation in \emph{Mathematical Culture} 4, 2010} }

@article{lagrangemixed, author = {Arnold, Douglas N. and Rognes, Marie}, title = {Stability of {L}agrange elements for the mixed {L}aplacian}, journal = {Calcolo}, fjournal = {Calcolo. A Quarterly on Numerical Analysis and Theory of Computation}, volume = {46}, year = {2009}, pages = {245--260}, abstract = { The stability properties of simple element choices for the mixed formulation of the Laplacian are investigated numerically. The element choices studied use vector Lagrange elements, i.e., the space of continuous piecewise polynomials vector fields of degree at most r, for the vector variable, and divergence of this space, which consists of discontinuous piecewise polynomials of one degree lower, for the scalar variable. For polynomial degrees r equal 2 or 3, this pair of spaces was found to be stable for all mesh families tested. In particular, it is stable on diagonal mesh families, in contrast to its behaviour for the Stokes equations. For degree r equal 1, stability holds for some meshes, but not for others. Additionally, convergence was observed precisely for the methods that were observed to be stable. However, it seems that optimal order L2 estimates for the vector variable, known to hold for r>3, do not hold for lower degrees. }, keywords = {mixed finite elements, Lagrange finite elements, stability}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/lagrangemixed.pdf}, doi = {10.1007/s10092-009-0009-6} }

@article{decomp, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Geometric decompositions and local bases for spaces of finite element differential forms}, journal = {Comput. Methods Appl. Mech. Engrg.}, fjournal = {Computer Methods in Applied Mechanics and Engineering}, volume = {198}, year = {2009}, pages = {1660--1672}, abstract = { We study the two primary families of spaces of finite element differential forms with respect to a simplicial mesh in any number of space dimensions. These spaces are generalizations of the classical finite element spaces for vector fields, frequently referred to as Raviart--Thomas, Brezzi--Douglas--Marini, and N\'ed\'elec spaces. In the present paper, we derive geometric decompositions of these spaces which lead directly to explicit local bases for them, generalizing the Bernstein basis for ordinary Lagrange finite elements. The approach applies to both families of finite element spaces, for arbitrary polynomial degree, arbitrary order of the differential forms, and an arbitrary simplicial triangulation in any number of space dimensions. A prominent role in the construction is played by the notion of a consistent family of extension operators, which expresses in an abstract framework a sufficient condition for deriving a geometric decomposition of a finite element space leading to a local basis. }, keywords = {finite element exterior calculus, finite element bases, Bernstein bases}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/decomp.pdf}, doi = {10.1016/j.cma.2008.12.017} }

@article{moebius, author = {Arnold, Douglas N. and Rogness, Jonathan}, title = {M\"obius transformations revealed}, journal = {Notices Amer. Math. Soc.}, fjournal = {Notices of the American Mathematical Society}, volume = {55}, year = {2008}, pages = {1226--1231}, abstract = { The authors' video, also titled Moebius Transformations Revealed, presents a visualization of these complex functions. It has attracted the attention of the general public as well as a technical audience. Here, we explain the mathematics involved, as well as how the video was produced. }, keywords = {Moebius transformations, complex mappings, stereographic projection, visualization}, subjclass = {30-01}, pdf = {http://umn.edu/~arnold/papers/moebius.pdf}, note = {Chinese translation in \emph{Shu Xue Yi Lin}, 28(2):109--115, 2009} }

@article{hungry, author = {Arnold, Douglas N.}, title = {Is the public hungry for math?}, journal = {Notices Amer. Math. Soc.}, fjournal = {Notices of the American Mathematical Society}, volume = {55}, year = {2008}, pages = {1069}, abstract = { Opinion column concerning public interest in mathematics }, pdf = {http://umn.edu/~arnold/papers/hungry.pdf}, note = {Chinese translation in \emph{Shu Xue Yi Lin}, 28(2):171--173, 2009} }

@article{elas3dfamily, author = {Arnold, Douglas N. and Awanou, Gerard and Winther, Ragnar}, title = {Finite elements for symmetric tensors in three dimensions}, journal = {Math. Comput.}, fjournal = {Mathematics of Computation}, volume = {77}, year = {2008}, pages = {1229--1251}, abstract = { We construct finite element subspaces of the space of symmetric tensors with square-integrable divergence on a three-dimensional domain. These spaces can be used to approximate the stress field in the classical Hellinger--Reissner mixed formulation of the elasticty equations, when standard discontinous finite element spaces are used to approximate the displacement field. These finite element spaces are defined with respect to an arbitrary simplicial triangulation of the domain, and there is one for each positive value of the polynomial degree used for the displacements. For each degree, these provide a stable finite element discretization. The construction of the spaces is closely tied to discretizations of the elasticity complex, and can be viewed as the three-dimensional analogue of the triangular element family for plane elasticity previously proposed by Arnold and Winther. }, keywords = {finite element, elasticity, mixed method}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/elas3dfamily.pdf}, doi = {10.1090/S0025-5718-08-02071-1} }

@article{fedf, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Finite element differential forms}, journal = {Proc. Appl. Math. Mech.}, fjournal = {Proceedings in Applied Mathematics and Mechanics}, volume = {7}, year = {2007}, pages = {1021901--1021902}, abstract = { A differential form is a field which assigns to each point of a domain an alternating multilinear form on its tangent space. The exterior derivative operation, which maps differential forms to differential forms of the next higher order, unifies the basic first order differential operators of calculus, and is a building block for a great variety of differential equations. When discretizing such differential equations by finite element methods, stable discretization depends on the development of spaces of finite element differential forms. As revealed recently through the finite element exterior calculus, for each order of differential form, there are two natural families of finite element subspaces associated to a simplicial triangulation. In the case of forms of order zero, which are simply functions, these two families reduce to one, which is simply the well-known family of Lagrange finite element subspaces of the first order Sobolev space. For forms of degree $1$ and of degree $n-1$ (where $n$ is the space dimension), we obtain two natural families of finite element subspaces, unifying many of the known mixed finite element spaces developed over the last decades. }, keywords = {mixed method, finite element, differential form, exterior calculus}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/fedf.pdf}, doi = {10.1002/pamm.200700549} }

@article{mixedelas3d, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Mixed finite element methods for linear elasticity with weakly imposed symmetry}, journal = {Math. Comput.}, fjournal = {Mathematics of Computation}, volume = {76}, year = {2007}, pages = {1699--1723}, abstract = { In this paper, we construct new finite element methods for the approximation of the equations of linear elasticity in three space dimensions that produce direct approximations to both stresses and displacements. The methods are based on a modified form of the Hellinger--Reissner variational principle that only weakly imposes the symmetry condition on the stresses. Although this approach has been previously used by a number of authors, a key new ingredient here is a constructive derivation of the elasticity complex starting from the de~Rham complex. By mimicking this construction in the discrete case, we derive new mixed finite elements for elasticity in a systematic manner from known discretizations of the de Rham complex. These elements appear to be simpler than the ones previously derived. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field. }, keywords = {mixed method, finite element, elasticity}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/mixedelas3d.pdf}, doi = {10.1090/S0025-5718-07-01998-9} }

@article{ecbc, author = {Arnold, Douglas N. and Tarfulea, Nicolae}, title = {Boundary conditions for the {E}instein-{C}hristoffel formulation of {E}instein's equations}, journal = {Electron. J. Differential Equations}, fjournal = {Electronic Journal of Differential Equations}, volume = {Conf.~15}, year = {2007}, pages = {11--27}, abstract = { Specifying boundary conditions continues to be a challenge in numerical relativity in order to obtain a long time convergent numerical simulation of Einstein's equations in domains with artificial boundaries. In this paper, we address this problem for the Einstein--Christoffel (EC) symmetric hyperbolic formulation of Einstein's equations linearized around flat spacetime. First, we prescribe simple boundary conditions that make the problem well posed and preserve the constraints. Next, we indicate boundary conditions for a system that extends the linearized EC system by including the momentum constraints and whose solution solves Einstein's equations in a bounded domain. Finally, we extend our results to the case of inhomogeneous boundary conditions. }, keywords = {general relativity, Einstein equations, boundary condition}, subjclass = {35Q75, 35L50, 83C99}, pdf = {http://umn.edu/~arnold/papers/ecbc.pdf} }

@article{dgrm2, author = {Arnold, Douglas N. and Brezzi, Franco and Falk, Richard S. and Marini, Donatella}, title = {Locking-free {R}eissner--{M}indlin elements without reduced integration}, journal = {Comput. Methods Appl. Mech. Engrg.}, fjournal = {Computer Methods in Applied Mechanics and Engineering}, volume = {96}, year = {2007}, pages = {3660--3671}, abstract = { In a recent paper of Arnold, Brezzi, and Marini, the ideas of discontinuous Galerkin methods were used to obtain and analyze two new families of locking free finite element methods for the approximation of the Reissner--Mindlin plate problem. By following their basic approach, but making different choices of finite element spaces, we develop and analyze other families of locking free finite elements that eliminate the need for the introduction of a reduction operator, which has been a central feature of many locking-free methods. For k>1, all the methods use piecewise polynomials of degree k to approximate the transverse displacement and (possibly subsets) of piecewise polynomials of degree k-1 to approximate both the rotation and shear stress vectors. The approximation spaces for the rotation and the shear stress are always identical. The methods vary in the amount of interelement continuity required. In terms of smallest number of degrees of freedom, the simplest method approximates the transverse displacement with continuous, piecewise quadratics and both the rotation and shear stress with rotated linear Brezzi-Douglas-Marini elements. }, keywords = {discontinuous Galerkin, Reissner-Mindlin plate, locking}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/dgrm2.pdf}, doi = {10.1016/j.cma.2006.10.023} }

@article{acta, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Finite element exterior calculus, homological techniques, and applications}, journal = {Acta Numer.}, fjournal = {Acta Numerica}, volume = {15}, year = {2006}, pages = {1--155}, abstract = { Finite element exterior calculus is an approach to the design and understanding of finite element discretizations for a wide variety of systems of partial differential equations. This approach brings to bear tools from differential geometry, algebraic topology, and homological algebra to develop discretizations which are compatible with the geometric, topological, and algebraic structures which underlie well-posedness of the PDE problem being solved. In the finite element exterior calculus, many finite element spaces are revealed as spaces of piecewise polynomial differential forms. These connect to each other in discrete subcomplexes of elliptic differential complexes, and are also related to the continuous elliptic complex through projections which commute with the complex differential. Applications are made to the finite element discretization of a variety of problems, including the Hodge Laplacian, Maxwell's equations, the equations of elasticity, and elliptic eigenvalue problems, and also to preconditioners. }, keywords = {finite element, exterior calculus, mixed method, Hodge, de Rham}, subjclass = {65N30, 58A10, 58A12, 58A14, 65N12, 65N15, 65N25, 65N55, 74S05, 78M20}, pdf = {http://umn.edu/~arnold/papers/acta.pdf}, doi = {10.1017/S0962492906210018} }

@incollection{elascomplexes2d, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Differential complexes and stability of finite element methods {II}: The elasticity complex}, booktitle = {Compatible Spatial Discretizations}, editor = {Arnold, D. and Bochev, P. and Lehoucq, R. and Nicolaides, R. and Shashkov, M.}, series = {IMA Vol. Math. Appl.}, volume = {142}, pages = {47--68}, publisher = {Springer}, address = {Berlin}, year = {2006}, abstract = { A close connection between the ordinary de Rham complex and a corresponding elasticity complex is utilized to derive new mixed finite element methods for linear elasticity. For a formulation with weakly imposed symmetry, this approach leads to methods which are simpler than those previously obtained. For example, we construct stable discretizations which use only piecewise linear elements to approximate the stress field and piecewise constant functions to approximate the displacement field. We also discuss how the strongly symmetric methods proposed in [8] can be derived in the present framework. The method of construction works in both two and three space dimensions, but for simplicity the discussion here is limited to the two dimensional case. }, keywords = {mixed finite element method, Hellinger-Reissner principle, elasticity}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/elascomplexes2d.pdf}, doi = {10.1007/0-387-38034-5_3} }

@incollection{nacomplexes, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Differential complexes and stability of finite element methods {I}: The de {R}ham complex}, booktitle = {Compatible Spatial Discretizations}, editor = {Arnold, D. and Bochev, P. and Lehoucq, R. and Nicolaides, R. and Shashkov, M.}, series = {The IMA Volumes in Mathematics and its Applications}, volume = {142}, pages = {23--46}, publisher = {Springer}, address = {Berlin}, year = {2006}, abstract = { In this paper we explain the relation between certain piecewise polynomial subcomplexes of the de Rham complex and the stability of mixed finite element methods for elliptic problems. }, keywords = {mixed finite element method, de Rham complex, stability}, subjclass = {65N12}, pdf = {http://umn.edu/~arnold/papers/nacomplexes.pdf}, doi = {10.1007/0-387-38034-5_2} }

@article{mixedelasrect, author = {Arnold, Douglas N. and Awanou, Gerard}, title = {Rectangular mixed finite elements for elasticity}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, volume = {15}, year = {2005}, number = {9}, pages = {1417--1429}, issn = {0218-2025}, mrclass = {65N30 (74S05)}, mrnumber = {MR2166210 (2006f:65112)}, mrreviewer = {Carsten Carstensen}, abstract = { We present a family of stable rectangular mixed finite elements for plane elasticity. Each member of the family consists of a space of piecewise polynomials discretizing the space of symmetric tensors in which the stress field is sought, and another to discretize the space of vector fields in which the displacement is sought. These may be viewed as analogues in the case of rectangular meshes of mixed finite elements recently proposed for triangular meshes. As for the triangular case the elements are closely related to a discrete version of the elasticity differential complex. }, keywords = {mixed method, finite element, elasticity, rectangular}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/mixedelasrect.pdf}, doi = {10.1142/S0218202505000741} }

@article{dgrm, author = {Arnold, Douglas N. and Brezzi, Franco and Marini, L. Donatella}, title = {A family of discontinuous {G}alerkin finite elements for the {R}eissner--{M}indlin plate}, journal = {J. Sci. Comput.}, fjournal = {Journal of Scientific Computing}, volume = {22/23}, year = {2005}, pages = {25--45}, issn = {0885-7474}, coden = {JSCOEB}, mrclass = {65N30 (74K20 74S05)}, mrnumber = {MR2142189 (2006b:65160)}, mrreviewer = {Alessandro Veneziani}, abstract = { We develop a family of locking-free elements for the Reissner-Mindlin plate using Discontinuous Galerkin techniques, one for each odd degree, and prove optimal error estimates. A second family uses conforming elements for the rotations and nonconforming elements for the transverse displacement, generalizing the element of Arnold and Falk to higher degree. }, keywords = {discontinuous Galerkin, Reissner-Mindlin plate, locking}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/dgrm.pdf}, doi = {10.1007/s10915-004-4134-8} }

@article{vecquad, author = {Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S.}, title = {Quadrilateral {$H({\rm div})$} finite elements}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {42}, year = {2005}, number = {6}, pages = {2429--2451}, issn = {0036-1429}, mrclass = {65N30 (41A10 41A63)}, mrnumber = {MR2139400 (2006d:65129)}, mrreviewer = {Beny Neta}, abstract = { We consider the approximation properties of quadrilateral finite element spaces of vector fields defined by the Piola transform, extending results previously obtained for scalar approximation. The finite element spaces are constructed starting with a given finite dimensional space of vector fields on a square reference element, which is then transformed to a space of vector fields on each convex quadrilateral element via the Piola transform associated to a bilinear isomorphism of the square onto the element. For affine isomorphisms, a necessary and sufficientcondition for approximation of order r+1 in L2 is that each component ofthe given space of functions on the reference element contain all polynomial functions of total degree at most r. In the case of bilinear isomorphisms,the situation is more complicated and we give a precise characterization of what is needed for optimal order L2-approximation of the function and of its divergence. As applications, we demonstrate degradation of the convergence order on quadrilateral meshes as compared to rectangular meshes for some standard finite element approximations of H(div). We also derive new estimates for approximation by quadrilateral Raviart-Thomas elements (requiring less regularity) and propose a new quadrilateral finite element space which provides optimal order approximation in H(div). Finally, we demonstrate the theory with numerical computations of mixed and least squares finite element aproximations of the solution of Poisson's equation. }, keywords = {quadrilateral, finite element, approximation, mixed finite element}, subjclass = {65N30, 41A10, 41A25, 41A27, 41A63}, pdf = {http://umn.edu/~arnold/papers/vecquad.pdf}, doi = {10.1137/S0036142903431924} }

@article{einsteinhyper, author = {Alekseenko, Alexander M. and Arnold, Douglas N.}, title = {New first-order formulation for the {E}instein equations}, journal = {Phys. Rev. D (3)}, fjournal = {Physical Review. D. Third Series}, volume = {68}, year = {2003}, number = {6}, pages = {064013, 6}, issn = {0556-2821}, coden = {PRVDAQ}, mrclass = {83C05}, mrnumber = {MR2035939 (2004m:83012)}, mrreviewer = {Norbert Noutchegueme}, abstract = { We derive a new first-order formulation for Einstein's equations which involves fewer unknowns than other first-order formulations that have been proposed. The new formulation is based on the $3+1$ decomposition with arbitrary lapse and shift. In the reduction to first order form only 8 particular combinations of the 18 first derivatives of the spatial metric are introduced. In the case of linearization about Minkowski space, the new formulation consists of symmetric hyperbolic system in 14 unknowns, namely the components of the extrinsic curvature perturbation and the 8 new variables, from whose solution the metric perturbation can be computed by integration. }, keywords = {relativity, Einstein equations, symmetric hyperbolic}, subjclass = {PACS 04.20.Ex, 04.25.Dm}, pdf = {http://umn.edu/~arnold/papers/einsteinhyper.pdf}, doi = {10.1103/PhysRevD.68.064013} }

@article{hierarchical, author = {Arnold, Douglas N. and Madureira, Alexandre L.}, title = {Asymptotic estimates of hierarchical modeling}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, volume = {13}, year = {2003}, number = {9}, pages = {1325--1350}, issn = {0218-2025}, mrclass = {35J25 (35B20 74G10 74K20)}, mrnumber = {MR2005646 (2004j:35074)}, abstract = { In this paper we propose a way to analyze certain classes of dimension reduction models for elliptic problems in thin domains. We develop asymptotic expansions for the exact and model solutions, having the thickness as small parameter. The modeling error is then estimated by comparing the respective expansions, and the upper bounds obtained make clear the influence of the order of the model and the thickness on the convergence rates. The techniques developed here allows for estimates in several norms and semi-norms, and also interior estimates (which disregards boundary layers). }, keywords = {hierarchical modeling, dimension reduction, asymptotic estimates}, subjclass = {35C20}, pdf = {http://umn.edu/~arnold/papers/hierarchical.pdf}, doi = {10.1142/S0218202503002933} }

@incollection{icmsc, author = {Arnold, Douglas N. and Winther, Ragnar}, title = {Mixed finite elements for elasticity in the stress-displacement formulation}, booktitle = {Current trends in scientific computing (Xi'an, 2002)}, editor = {Chen, Z. and Glowinski, R. and Li, K.}, series = {Contemp. Math.}, volume = {329}, pages = {33--42}, publisher = {Amer. Math. Soc.}, address = {Providence, RI}, year = {2003}, mrclass = {65N30 (74S05)}, mrnumber = {MR2022629}, abstract = { We present a family of pairs of finite element spaces for the unaltered Hellinger--Reissner variational principle using polynomial shape functions on a single triangular mesh for stress and displacement. There is a member of the family for each polynomial degree, beginning with degree two for the stress and degree one for the displacement, and each is stable and affords optimal order approximation. The simplest element pair involves 24 local degrees of freedom for the stress and 6 for the displacement. We also construct a lower order element involving 21 stress degrees of freedom and 3 displacement degrees of freedom which is, we believe, likely to be the simplest possible conforming stable element pair with polynomial shape functions. For all these conforming elements the approximate stress not only belongs to H(div), but is also continuous at element vertices, which is more continuity than may be desired. We show that for conforming finite elements with polynomial shape functions, this additional continuity is unavoidable. To overcome this obstruction, we construct as well some non-conforming stable mixed finite elements, which we show converge with optimal order as well. The simplest of these involves only 12 stress and 6 displacement degrees of freedom on each triangle. }, keywords = {mixed method, finite element, elasticity}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/icmsc.pdf} }

@article{ncelas, author = {Arnold, Douglas N. and Winther, Ragnar}, title = {Nonconforming mixed elements for elasticity}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, volume = {13}, year = {2003}, number = {3}, pages = {295--307}, issn = {0218-2025}, mrclass = {65N30 (65N12 65N15 74G15 74S05)}, mrnumber = {MR1977627 (2004f:65176)}, mrreviewer = {Ji{\v{r}}{\'{\i}} Nedoma}, abstract = { We construct first order, stable, nonconforming mixed finite elements for plane elasticity and analyze their convergence. The mixed method is based on the Hellinger-Reissner variational formulation in which the stress and displacement fields are the primary unknowns. The stress elements use polynomial shape functions but do not involve vertex degrees of freedom. }, keywords = {mixed method, finite element, nonconforming, elasticity}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/ncelas.pdf}, doi = {10.1142/S0218202503002507} }

@inproceedings{icm2002, author = {Arnold, Douglas N.}, title = {Differential complexes and numerical stability}, booktitle = {Proceedings of the International Congress of Mathematicians, Vol. I (Beijing, 2002)}, editor = {Tatsien, L.}, pages = {137--157}, publisher = {Higher Ed. Press}, address = {Beijing}, year = {2002}, mrclass = {65N30 (65N12)}, mrnumber = {MR1989182 (2004h:65115)}, mrreviewer = {Andre Nicolet}, abstract = { Differential complexes such as the de Rham complex have recently come to play an important role in the design and analysis of numerical methods for partial differential equations. The design of stable discretizations of systems of partial differential equations often hinges on capturing subtle aspects of the structure of the system in the discretization. In many cases the differential geometric structure captured by a differential complex has proven to be a key element, and a discrete differential complex which is appropriately related to the original complex is essential. This new geometric viewpoint has provided a unifying understanding of a variety of innovative numerical methods developed over recent decades and pointed the way to stable discretizations of problems for which none were previously known, and it appears likely to play an important role in attacking some currently intractable problems in numerical PDE. }, keywords = {finite element, numerical stability, differential complex}, subjclass = {65N12}, pdf = {http://umn.edu/~arnold/papers/icm2002.pdf}, note = {Spanish translation in \emph{Gac.\ R.\ Soc.\ Mat.\ Esp.}, 8(2):335-360, 2005} }

@incollection{rmquad, author = {Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S.}, title = {Remarks on quadrilateral {R}eissner--{M}indlin plate elements}, booktitle = {WCCM V - Fifth World Congress on Computational Mechanics}, publisher = {Technical University of Vienna}, editor = {Mang, H. and Rammerstorfer, F. and Eberhardsteiner, J.}, pages = {137--157}, year = {2002}, abstract = { Over the last two decades, there has been an extensive effort to devise and analyze finite elements schemes for the approximation of the Reissner�Mindlin plate equations which avoid locking, numerical overstiffness resulting in a loss of accuracy when the plate is thin. There are now many triangular and rectangular finite elements, for which a mathematical analysis exists to certify them as free of locking. Generally speaking, the analysis for rectangular elements extends to the case of parallograms, which are defined by affine mappings of rectangles. However, for more general convex quadrilaterals, defined by bilinear mappings of rectangles, the analysis is more complicated. Recent results by the authors on the approximation properties of quadrilateral finite elements shed some light on the problems encountered. In particular, they show that for some finite element methods for the approximation of the Reissner-Mindlin plate, the obvious generalization of rectangular elements to general quadrilateral meshes produce methods which lose accuracy. In this paper, we present an overview of this situation. }, keywords = {Reissner-Mindlin plate, finite element, locking-free, isoparametric}, subjclass = {65N30,74K20}, pdf = {http://umn.edu/~arnold/papers/rmquad.pdf} }

@article{rmkl, author = {Arnold, Douglas N. and Madureira, Alexandre L. and Zhang, Sheng}, title = {On the range of applicability of the {R}eissner--{M}indlin and {K}irchhoff-{L}ove plate bending models}, journal = {J. Elasticity}, fjournal = {Journal of Elasticity. The Physical and Mathematical Science of Solids}, volume = {67}, year = {2002}, number = {3}, pages = {171--185 (2003)}, issn = {0374-3535}, coden = {JELSAY}, mrclass = {74K20 (74B05 74G99)}, mrnumber = {MR1997951 (2004e:74053)}, mrreviewer = {George Jaiani}, abstract = { We show that the Reissner-Mindlin plate bending model has a wider range of applicability than the Kirchhoff-Love model for the approximation of clamped linearly elastic plates. Under the assumption that the body force density is constant in the transverse direction, the Reissner-Mindlin model solution converges to the three-dimensional linear elasticity solution in the relative energy norm for the full range of surface loads. However, for loads with a significant transverse shear effect, the Kirchhoff-Love model fails. }, keywords = {plate, Reissner-Mindlin, Kirchhoff-Love}, subjclass = {73K10, 73C02}, pdf = {http://umn.edu/~arnold/papers/rmkl.pdf}, doi = {10.1023/A:1024986427134} }

@article{dgerr, author = {Arnold, Douglas N. and Brezzi, Franco and Cockburn, Bernardo and Marini, L. Donatella}, title = {Unified analysis of discontinuous {G}alerkin methods for elliptic problems}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {39}, year = {2002}, number = {5}, pages = {1749--1779}, issn = {0036-1429}, mrclass = {65N30}, mrnumber = {MR1885715 (2002k:65183)}, abstract = { We provide a framework for the analysis of a large class of discontinuous methods for second-order elliptic problems. It allows for the understanding and comparison of most of the discontinuous Galerkin methods that have been proposed for the numerical treatment of elliptic problems by diverse communities over three decades. }, keywords = {elliptic problems, discontinuous Galerkin, interior penalty}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/dgerr.pdf}, doi = {10.1137/S0036142901384162} }

@article{mixedelas, author = {Arnold, Douglas N. and Winther, Ragnar}, title = {Mixed finite elements for elasticity}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {92}, year = {2002}, number = {3}, pages = {401--419}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {65N30 (65N12 65N15 74G15 74S05)}, mrnumber = {MR1930384 (2003i:65103)}, mrreviewer = {Carsten Carstensen}, abstract = { There have been many efforts, dating back four decades, to develop stable mixed finite elements for the stress-displacement formulation of the plane elasticity system. This requires the development of a compatible pair of finite element spaces, one to discretize the space of symmetric tensors in which the stress field is sought, and one to discretize the space of vector fields in which the displacement is sought. Although there are number of well-known mixed finite element pairs known for the analogous problem involving vector fields and scalar fields, the symmetry of the stress field is a substantial additional difficulty, and the elements presented here are the first ones using polynomial shape functions which are known to be stable. We present a family of such pairs of finite element spaces, one for each polynomial degree, beginning with degree two for the stress and degree one for the displacement, and show stability and optimal order approximation. We also analyze some obstructions to the construction of such finite element spaces, which account for the paucity of elements available. }, keywords = {mixed method, finite element, elasticity}, subjclass = {65N30, 74S05}, pdf = {http://umn.edu/~arnold/papers/mixedelas.pdf}, doi = {10.1007/s002110100348} }

@article{quadapprox, author = {Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S.}, title = {Approximation by quadrilateral finite elements}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {71}, year = {2002}, number = {239}, pages = {909--922}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65N30}, mrnumber = {MR1898739 (2003c:65112)}, abstract = { We consider the approximation properties of finite element spaces on quadrilateral meshes. The finite element spaces are constructed starting with a given finite dimensional space of functions on a square reference element, which is then transformed to a space of functions on each convex quadrilateral element via a bilinear isomorphism of the square onto the element. It is known that for affine isomorphisms, a necessary and sufficient condition for approximation of order r+1 in L2 and order r in H1 is that the given space of functions on the reference element contain all polynomial functions of total degree at most r. In the case of bilinear isomorphisms, it is known that the same estimates hold if the function space contains all polynomial functions of separate degree r. We show, by means of a counterexample, that this latter condition is also necessary. As applications we demonstrate degradation of the convergence order on quadrilateral meshes as compared to rectangular meshes for serendipity finite elements and for various mixed and nonconforming finite elements. }, keywords = {quadrilateral, finite element, approximation, serendipity, mixed finite element}, subjclass = {65N30, 41A10, 41A25, 41A27, 41A63}, pdf = {http://umn.edu/~arnold/papers/quadapprox.pdf}, doi = {10.1090/S0025-5718-02-01439-4} }

@article{quadmeshes, author = {Arnold, Douglas N. and Boffi, Daniele and Falk, Richard S. and Gastaldi, Lucia}, title = {Finite element approximation on quadrilateral meshes}, journal = {Comm. Numer. Methods Engrg.}, fjournal = {Communications in Numerical Methods in Engineering}, volume = {17}, year = {2001}, number = {11}, pages = {805--812}, issn = {1069-8299}, coden = {CANMER}, mrclass = {65N30 (65N50)}, mrnumber = {MR1872639 (2002i:65121)}, abstract = { Quadrilateral finite elements are generally constructed by starting from a given finite dimensional space of polynomials V^ on the unit reference square K^. The elements of V^ are then transformed by using the bilinear isomorphisms F_K which map K^ to each convex quadrilateral element K. It has been recently proven that a necessary and sufficient condition for approximation of order r+1 in L^2 and r in H^1 is that V^ contains the space Q_r of all polynomial functions of degree r separately in each variable. In this paper several numerical experiments are presented which confirm the theory. The tests are taken from various examples of applications: the Laplace operator, the Stokes problem and an eigenvalue problem arising in fluid-structure interaction modeling. }, keywords = {quadrilateral, finite element, approximation, serendipity, mixed finite element}, subjclass = {65N30, 41A10, 41A25, 41A27, 41A63}, pdf = {http://umn.edu/~arnold/papers/quadmeshes.pdf}, doi = {10.1002/cnm.450} }

@incollection{enumath, author = {Arnold, Douglas N.}, title = {Numerical problems in general relativity}, booktitle = {Numerical Mathematics and Advanced Applications}, editor = {Neittaanm\"ki, P. and Tiihonen, T. and Tarvainen, P.}, publisher = {World Scientific}, year = {2000}, abstract = { The construction of gravitational wave observatories is one of the greatest scientific efforts of our time. As a result, there is presently a strong need to numerically simulate the emission of gravitation radiation from massive astronomical events such as black hole collisions. This entails the numerical solution of the Einstein field equations. We briefly describe the field equations in their natural setting, namely as statements about the geometry of space time. Next we describe the complicated system that arises when the field equations are recast as partial differential equations, and discuss procedures for deriving from them a more tractable system consisting of constraint equations to be satisfied by initial data and together with evolution equations. We present some applications of modern finite element technology to the solution of the constraint equations in order to find initial data relevant to black hole collisions. We conclude by enumerating some of the many computational challenges that remain. }, keywords = {general relativity, gravitational radiation, black holes, Einstein equations}, subjclass = {65N30, 83C05, 83C35, 83C57}, pdf = {http://umn.edu/~arnold/papers/enumath.pdf} }

@incollection{dg, author = {Arnold, Douglas N. and Brezzi, Franco and Cockburn, Bernardo and Marini, Donatella}, title = {Discontinuous {G}alerkin methods for elliptic problems}, booktitle = {Discontinuous Galerkin methods (Newport, RI, 1999)}, editor = {Cockburn, B. and Karniadakis, G. and Shu, C.}, series = {Lect. Notes Comput. Sci. Eng.}, volume = {11}, pages = {89--101}, publisher = {Springer}, address = {Berlin}, year = {2000}, mrclass = {65N30}, mrnumber = {MR1842165 (2002d:65118)}, abstract = { We provide a common framework for the understanding, comparison, and analysis of several discontinuous Galerkin methods that have been proposed for the numerical treatment of elliptic problems. This class includes the recently introduced methods of Bassi and Rebay (together with the variants proposed by Brezzi, Manzini, Marini, Pietra and Russo), the local discontinuous Galerkin methods of Cockburn and Shu, and the method of Baumann and Oden. It also includes the so-called interior penalty methods developed some time ago by Douglas and Dupont, Wheeler, Baker, and Arnold among others. }, pdf = {http://umn.edu/~arnold/papers/dg.pdf}, doi = {10.1007/978-3-642-59721-3_5} }

@article{hdivcurl, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Multigrid in {$H({\rm div})$} and {$H({\rm curl})$}}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {85}, year = {2000}, number = {2}, pages = {197--217}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {65N55}, mrnumber = {MR1754719 (2001d:65161)}, mrreviewer = {Do Young Kwak}, abstract = { We consider the solution of systems of linear algebraic equations which arise from the finite element discretization of variational problems posed in the Hilbert spaces H(div) and H(curl) in three dimensions. We show that if appropriate finite element spaces and appropriate additive or multiplicative Schwarz smoothers are used, then the multigrid V-cycle is an efficient solver and preconditioner for the discrete operator. All results are uniform with respect to the mesh size, the number of mesh levels, and weights on the two terms in the inner products. }, keywords = {multigrid, preconditioner, mixed method, finite element}, subjclass = {65N55, 65N30}, pdf = {http://umn.edu/~arnold/papers/hdivcurl.pdf}, doi = {10.1007/PL00005386} }

@article{bistet, author = {Arnold, Douglas N. and Mukherjee, Arup and Pouly, Luc}, title = {Locally adapted tetrahedral meshes using bisection}, journal = {SIAM J. Sci. Comput.}, fjournal = {SIAM Journal on Scientific Computing}, volume = {22}, year = {2000}, number = {2}, pages = {431--448}, issn = {1064-8275}, mrclass = {65N50 (65N30)}, mrnumber = {MR1780608 (2002h:65204)}, mrreviewer = {Zden{\v{e}}k Dost{\'a}l}, abstract = { We present an algorithm for the construction of locally adapted conformal tetrahedral meshes. The algorithm is based on bisection of tetrahedra. A new data structure is introduced, which simplifies both the selection of the refinement edge of a tetrahedron and the recursive refinement to conformity of a mesh once some tetrahedra have been bisected. We prove that repeated application of the algorithm leads to only finitely many tetrahedral shapes up to similarity, and bound the amount of additional refinement that is needed to achieve conformity. Numerical examples of the effectiveness of the algorithm are presented. }, keywords = {bisection, tetrahedral meshes, adaptive refinement, similarity classes, finite elements}, subjclass = {65N50}, pdf = {http://umn.edu/~arnold/papers/bistet.pdf}, doi = {10.1137/S1064827597323373} }

@incollection{bistetima, author = {Arnold, Douglas N. and Mukherjee, Arup}, title = {Tetrahedral bisection and adaptive finite elements}, booktitle = {Grid generation and adaptive algorithms (Minneapolis, MN, 1997)}, editor = {Bern, M. and Flahery, J. and Luskin, M.}, series = {IMA Vol. Math. Appl.}, volume = {113}, pages = {29--42}, publisher = {Springer}, address = {New York}, year = {1999}, mrclass = {65N30 (65N50)}, mrnumber = {MR1705561 (2000f:65131)}, abstract = { An adaptive finite element algorithm for elliptic boundary value problems in $\R^3$ is presented. The algorithm uses linear finite elements, a-posteriori error estimators, a mesh refinement scheme based on bisection of tetrahedra, and a multi-grid solver. We show that the repeated bisection of an arbitrary tetrahedron leads to only a finite number of dissimilar tetrahedra, and that the recursive algorithm ensuring conformity of the meshes produced terminates in a finite number of steps. A procedure for assigning numbers to tetrahedra in a mesh based on a-posteriori error estimates, indicating the degree of refinement of the tetrahedron, is also presented. Numerical examples illustrating the effectiveness of the algorithm are given. }, keywords = {finite elements, adaptive mesh refinement, error estimators, bisection of tetrahedra}, subjclass = {65N50}, pdf = {http://umn.edu/~arnold/papers/bistetima.pdf}, doi = {10.1007/978-1-4612-1556-1_2} }

@incollection{dundee, author = {Arnold, D. N. and Mukherjee, A. and Pouly, L.}, title = {Adaptive finite elements and colliding black holes}, booktitle = {Numerical analysis 1997 (Dundee)}, editor = {Griffiths, D. and Higham, D. and Watson, G.}, series = {Pitman Res. Notes Math. Ser.}, volume = {380}, pages = {1--15}, publisher = {Longman}, address = {Harlow}, year = {1998}, mrclass = {83C57 (65N55 83-08)}, mrnumber = {MR1636915 (99f:83039)}, abstract = { According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as the collision of black holes, should generate gravity waves of sufficient amplitude to detect on earth. The massive Laser Interferometer Gravitational-Wave Observatory, or LIGO, is now being constructed to detect gravity waves. Consequently there is great interest in the computer simulation of black hole collisions and similar events, based on the numerical solution of the Einstein field equations. In this note we introduce the scientific, mathematical, and computational problems and discuss the development of a computer code to solve the initial data problem for colliding black holes, a nonlinear elliptic boundary value problem posed in an unbounded three dimensional domain which is a key step in solving the full field equations. The code is based on finite elements, adaptive meshes, and a multigrid solution process. Here we will particularly emphasize the mathematical and algorithmic issues arising in the generation of adaptive tetrahedral meshes. }, keywords = {adaptivity, finite elements, black holes, Einstein equations}, subjclass = {65N50, 65N30, 83C05, 83C35, 83C57}, pdf = {http://umn.edu/~arnold/papers/dundee.pdf} }

@article{hdivn, author = {Arnold, D. N. and Falk, R. S. and Winther, R.}, title = {Multigrid preconditioning in {$H({\rm div})$} on non-convex polygons}, journal = {Comput. Appl. Math.}, fjournal = {Computational and Applied Mathematics}, volume = {17}, year = {1998}, number = {3}, pages = {303--315}, issn = {0101-8205}, mrclass = {65N55 (65F10)}, mrnumber = {MR1687885 (2000k:65224)}, mrreviewer = {Zden{\v{e}}k Dost{\'a}l}, abstract = { In an earlier paper we constructed and analyzed a multigrid preconditioner for the system of linear algebraic equations arising from the finite element discretization of boundary value problems associated to the differential operator I - grad div. In this paper we analyze the procedure without assuming that the underlying domain is convex and show that, also in this case, the preconditioner is spectrally equivalent to the inverse of the discrete operator. }, keywords = {preconditioner, finite element, multigrid, nonconvex domain}, subjclass = {65N55, 65N30}, pdf = {http://umn.edu/~arnold/papers/hdivn.pdf} }

@incollection{santiago, author = {Alessandrini, Stephan M. and Arnold, Douglas N. and Falk, Richard S. and Madureira, Alexandre L.}, title = {Dimensional reduction for plates based on mixed variational principles}, booktitle = {Shells (Santiago de Compostela, 1997)}, series = {Cursos Congr. Univ. Santiago de Compostela}, volume = {105}, pages = {25--28}, publisher = {Univ. Santiago de Compostela}, address = {Santiago de Compostela}, year = {1997}, mrclass = {74K20}, mrnumber = {MR1786907}, abstract = { We consider the derivation and rigorous justification of models for thin linearly elastic plates using mixed variational principles. }, keywords = {plate, dimensional reduction, mixed variational principle}, subjclass = {73K10, 73C02}, pdf = {http://umn.edu/~arnold/papers/santiago.pdf} }

@incollection{prermdd9, author = {Arnold, D. N. and Falk, R. S. and Winther, R.}, title = {Preconditioning discrete approximations of the {R}eissner--{M}indlin plate model}, booktitle = {Ninth International Conference on Domain Decomposition Methods}, editor = {Bj{\o}rstad, P. and Espedal, M. and Keyes, D.}, pages = {215--221}, publisher = {DDM.org}, year = {1996}, abstract = { We consider iterative methods for the solution of linear systems of equations arising from mixed finite element discretization of the Reissner-Mindlin plate model. We show how to construct symmetric positive definite block diagonal preconditioners for these indefinite systems such that the resulting systems have spectral condition numbers independent of both the mesh size h and the plate thickness t. }, keywords = {preconditioner, Reissner, Mindlin, plate, finite element}, subjclass = {65N30, 65N22, 65F10, 73V05}, pdf = {http://umn.edu/~arnold/papers/prermdd9.pdf} }

@incollection{hdivdd9, author = {Arnold, D. N. and Falk, R. S. and Winther, R.}, title = {Preconditioning in {$H$}(div) and applications}, booktitle = {Ninth International Conference on Domain Decomposition Methods}, editor = {Bj{\o}rstad, P. and Espedal, M. and Keyes, D.}, pages = {12--19}, publisher = {DDM.org}, year = {1996}, abstract = { Summarizing the work of [1], we show how to construct preconditioners using domain decomposition and multigrid techniques for the system of linear algebraic equations which arises from the finite element discretization of boundary value problems associated to the differential operator I - grad div. These preconditioners are shown to be spectrally equivalent to the inverse of the operator and thus may be used to precondition iterative methods so that any given error reduction may be achieved in a finite number of iterations independent of the mesh discretization. We describe applications of these results to the efficient solution of mixed and least squares finite element approximations of elliptic boundary value problems. }, keywords = {preconditioner, mixed method, least squares, finite element,}, subjclass = {65N55, 65N30}, pdf = {http://umn.edu/~arnold/papers/hdivdd9.pdf} }

@incollection{platederiv, author = {Alessandrini, Stephen M. and Arnold, Douglas N. and Falk, Richard S. and Madureira, Alexandre L.}, title = {Derivation and justification of plate models by variational methods}, booktitle = {Plates and shells (Qu\'ebec, QC, 1996)}, editor = {Fortin, M.}, series = {CRM Proc. Lecture Notes}, volume = {21}, pages = {1--20}, publisher = {Amer. Math. Soc.}, address = {Providence, RI}, year = {1999}, mrclass = {74K20}, mrnumber = {MR1696513 (2000j:74055)}, mrreviewer = {Petr E. Tovstik}, abstract = { We consider the derivation of two-dimensional models for the bending and stretching of a thin three-dimensional linearly elastic plate using variational methods. Specifically we consider restriction of the trial space in two different forms of the Hellinger-Reissner variational principle for 3-D elasticity to functions with a specified polynomial dependence in the transverse direction. Using this approach many different plate models are possible and we classify and investigate the most important. We study in detail a method which leads naturally not only to familiar plate models, but also to error bounds between the plate solution and the full 3-D solution. }, keywords = {plate, dimensional reduction, Reissner-Mindlin}, subjclass = {73K10, 73C02}, pdf = {http://umn.edu/~arnold/papers/platederiv.pdf} }

@incollection{cai, author = {Arnold, Douglas N.}, title = {Computer-aided instruction}, booktitle = {Encarta Encyclopedia}, publisher = {Microsoft}, year = {1997}, abstract = { Encyclopedia article on computer-aided instruction }, pdf = {http://umn.edu/~arnold/papers/cai.pdf} }

@article{prerm, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, Ragnar}, title = {Preconditioning discrete approximations of the {R}eissner--{M}indlin plate model}, journal = {RAIRO Mod\'el. Math. Anal. Num\'er.}, fjournal = {RAIRO Mod\'elisation Math\'ematique et Analyse Num\'erique}, volume = {31}, year = {1997}, number = {4}, pages = {517--557}, issn = {0764-583X}, mrclass = {73K10 (73V05)}, mrnumber = {MR1457459 (99b:73028)}, abstract = { We consider iterative methods for the solution of the linear system of equations arising from the mixed finite element discretization of the Reissner-Mindlin plate model. We show how to construct a symmetric positive definite block diagonal preconditioner such that the resulting linear system has spectral condition number independent of both the mesh size h and the plate thickness t. We further discuss how this preconditioner may be implemented and then apply it to efficiently solve this indefinite linear system. Although the mixed formulation of the Reissner-Mindlin problem has a saddle-point structure common to other mixed variational problems, the presence of the small parameter t and the fact that the matrix in the upper left corner of the partition is only positive semidefinite introduces new complications. }, keywords = {preconditioner, Reissner, Mindlin, plate, finite element}, subjclass = {65N30, 65N22, 65F10, 73V05}, pdf = {http://umn.edu/~arnold/papers/prerm.pdf}, doi = {10.1051/m2an/1997310405171} }

@article{hdiv, author = {Arnold, Douglas N. and Falk, Richard S. and Winther, R.}, title = {Preconditioning in {$H({\rm div})$} and applications}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {66}, year = {1997}, number = {219}, pages = {957--984}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65N55 (65N30)}, mrnumber = {MR1401938 (97i:65177)}, mrreviewer = {O. Widlund}, abstract = { We consider the solution of the system of linear algebraic equations which arises from the finite element discretization of boundary value problems associated to the differential operator I - grad div. The natural setting for such problems is in the Hilbert space H(div) and the variational formulation is based on the inner product in H(div). We show how to construct preconditioners for these equations using both domain decomposition and multigrid techniques. These preconditioners are shown to be spectrally equivalent to the inverse of the operator. As a consequence, they may be used to precondition iterative methods so that any given error reduction may be achieved in a finite number of iterations, with the number independent of the mesh discretization. We describe applications of these results to the efficient solution of mixed and least squares finite element approximations of elliptic boundary value problems. }, keywords = {preconditioner, mixed method, least squares, finite element, multigrid, domain decomposition}, subjclass = {65N55, 65N30}, pdf = {http://umn.edu/~arnold/papers/hdiv.pdf}, doi = {10.1090/S0025-5718-97-00826-0} }

@article{ozf, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {Analysis of a linear-linear finite element for the {R}eissner--{M}indlin plate model}, journal = {Math. Models Methods Appl. Sci.}, fjournal = {Mathematical Models \& Methods in Applied Sciences}, volume = {7}, year = {1997}, number = {2}, pages = {217--238}, issn = {0218-2025}, mrclass = {73V05 (65M60 73K10)}, mrnumber = {MR1440607 (98b:73034)}, mrreviewer = {Rolf Rannacher}, abstract = { An analysis is presented for a recently proposed finite element method for the Reissner-Mindlin plate problem. The method is based on the standard variational principle, uses nonconforming linear elements to approximate the rotations and conforming linear elements to approximate the transverse displacements, and avoids the usual "locking problem" by interpolating the shear stress into a rotated space of lowest order Raviart-Thomas elements. When the plate thickness t=O(h), it is proved that the method gives optimal order error estimates uniform in t. However, the analysis suggests and numerical calculations confirm that the method can produce poor approximations for moderate sized values of the plate thickness. Indeed, for t fixed, the method does not converge as the mesh size h tends to zero. }, keywords = {Reissner, Mindlin, plate, finite element, nonconforming}, subjclass = {65N30, 73K10, 73K25}, pdf = {http://umn.edu/~arnold/papers/ozf.pdf}, doi = {10.1142/S0218202597000141} }

@article{rmint, author = {Arnold, Douglas N. and Liu, Xiaobo}, title = {Interior estimates for a low order finite element method for the {R}eissner--{M}indlin plate model}, journal = {Adv. Comput. Math.}, fjournal = {Advances in Computational Mathematics}, volume = {7}, year = {1997}, number = {3}, pages = {337--360}, issn = {1019-7168}, mrclass = {73V05 (65N30 73K10)}, mrnumber = {MR1449685 (98i:73053)}, mrreviewer = {Juhani Pitk{\"a}ranta}, abstract = { Interior error estimates are obtained for a low order finite element introduced by Arnold and Falk for the Reissner-Mindlin plates. It is proved that the approximation error of the finite element solution in the interior domain is bounded above by two parts: one measures the local approximability of the exact solution by the finite element space and the other the global approximability of the finite element method. As an application, we show that for the soft simply supported plate, the Arnold-Falk element still achieves an almost optimal convergence rate in the energy norm away from the boundary layer, even though optimal order convergence cannot hold globally due to the boundary layer. Numerical results are given which support our conclusion. }, keywords = {Reissner-Mindlin plate, boundary layer, mixed finite element,}, subjclass = {65N30,73K10}, pdf = {http://umn.edu/~arnold/papers/rmint.pdf}, doi = {10.1023/A:1018907205385} }

@article{psri, author = {Arnold, Douglas N. and Brezzi, Franco}, title = {The partial selective reduced integration method and applications to shell problems}, journal = {Comput. Structures}, fjournal = {Computers \& Structures}, volume = {64}, year = {1997}, pages = {879--880}, abstract = { We briefly present the main idea of partial selective reduced integration as developed in other works of the authors. The idea is quite general and can be applied to a number of different situations, but we concentrate on the case of the Naghdi shell model. }, keywords = {partial selective reduced integration, Naghdi shell}, subjclass = {65N30, 73K15, 73V05}, pdf = {http://umn.edu/~arnold/papers/psri.pdf}, doi = {10.1016/S0045-7949(96)00151-4} }

@article{shellelt, author = {Arnold, Douglas N. and Brezzi, Franco}, title = {Locking-free finite element methods for shells}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {66}, year = {1997}, number = {217}, pages = {1--14}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {73V05 (65N30 73K15)}, mrnumber = {MR1370847 (97c:73061)}, mrreviewer = {Juhani Pitk{\"a}ranta}, abstract = { We propose a new family of finite element methods for the Naghdi shell model, one method associated with each nonnegative integer k. The methods are based on a nonstandard mixed formulation, and the kth method employs triangular Lagrange finite elements of degree k+2 augmented by bubble functions of degree k+3 for both the displacement and rotation variables, and discontinuous piecewise polynomials of degree k for the shear and membrane stresses. This method can be implemented in terms of the displacement and rotation variables alone, as the minimization of an altered energy functional over the space mentioned. The alteration consists of the introduction of a weighted local projection into part, but not all, of the shear and membrane energy terms of the usual Naghdi energy. The relative error in the method, measured in a norm which combines the H1 norm of the displacement and rotation fields and an appropriate norm of the shear and membrane stress fields, converges to zero with order k+1 uniformly with respect to the shell thickness for smooth solutions, at least under the assumption that certain geometrical coefficients in the Nagdhi model are replaced by piecewise constants. }, keywords = {shell, locking, finite element}, subjclass = {65N30,73K15,73V05}, pdf = {http://umn.edu/~arnold/papers/shellelt.pdf}, doi = {10.1090/S0025-5718-97-00785-0} }

@article{bdlayer2, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {Asymptotic analysis of the boundary layer for the {R}eissner--{M}indlin plate model}, journal = {SIAM J. Math. Anal.}, fjournal = {SIAM Journal on Mathematical Analysis}, volume = {27}, year = {1996}, number = {2}, pages = {486--514}, issn = {0036-1410}, coden = {SJMAAH}, mrclass = {73K10 (73C02)}, mrnumber = {MR1377485 (97i:73064)}, abstract = { We investigate the structure of the solution of the Reissner-Mindlin plate equations in its dependence on the plate thickness for various boundary conditions, developing asymptotic expansions in powers of the plate thickness for the main physical quantities. These expansions are uniform up to the boundary for the transverse displacement, but for other variables there is a boundary layer, whose strength depends on the boundary conditions. We give rigorous error bounds for the errors in the expansions in Sobolev norms and make various applications. }, keywords = {Reissner, Mindlin, plate, boundary layer}, subjclass = {73K10, 73K25}, pdf = {http://umn.edu/~arnold/papers/bdlayer2.pdf}, doi = {10.1137/S0036141093245276} }

@article{stokes, author = {Arnold, Douglas N. and Liu, Xiao Bo}, title = {Local error estimates for finite element discretizations of the {S}tokes equations}, journal = {RAIRO Mod\'el. Math. Anal. Num\'er.}, fjournal = {RAIRO Mod\'elisation Math\'ematique et Analyse Num\'erique}, volume = {29}, year = {1995}, number = {3}, pages = {367--389}, issn = {0764-583X}, mrclass = {76M10 (65N30 76D07)}, mrnumber = {MR1342712 (96d:76055)}, mrreviewer = {Yves Achdou}, abstract = { Local error estimates are derived which apply to most stable finite mixed finite element discretizations of the stationary Stokes equations. }, keywords = {Stokes equations, mixed finite element method, local}, subjclass = {65N30, 65N15, 76M10, 76D07}, pdf = {http://umn.edu/~arnold/papers/stokes.pdf}, doi = {10.1051/m2an/1995290303671} }

@article{ncstokes, author = {Arnold, Douglas N.}, title = {On nonconforming linear-constant elements for some variants of the {S}tokes equations}, journal = {Istit. Lombardo Accad. Sci. Lett. Rend. A}, fjournal = {Istituto Lombardo. Accademia di Scienze e Lettere. Rendiconti. Scienze Matematiche e Applicazioni. A}, volume = {127}, year = {1993}, number = {1}, pages = {83--93 (1994)}, issn = {0021-2504}, mrclass = {65N30 (65M60 73V05 76D07 76M10)}, mrnumber = {MR1284845 (95c:65175)}, abstract = { Nonconforming piecewise linear finite elements for the velocity field and piecewise constant elements for the pressure field give a simple stable, optimal order approximation to the Stokes equations, but are not stable for the equations of incompressible elasticty, which differ from the Stokes equations only in that the vector Laplace operator is replaced by the Lame operator. However, we show that if we replace the divergence by the rotation, then the nonconforming linear-constant element is stable both for the system involving the Laplacian and for that involving the Lame operator. Finally we discuss an application to the Reissner-Mindlin plate. }, keywords = {Stokes equations, mixed finite element}, subjclass = {65N30,65N12,76M10,76D07,73K10}, pdf = {http://umn.edu/~arnold/papers/ncstokes.pdf} }

@incollection{rmelts, author = {Arnold, Douglas N. and Brezzi, Franco}, title = {Some new elements for the {R}eissner--{M}indlin plate model}, booktitle = {Boundary value problems for partial differential equations and applications}, editor = {Baiocchi, C. and Lions, J.-L.}, series = {RMA Res. Notes Appl. Math.}, volume = {29}, pages = {287--292}, publisher = {Masson}, address = {Paris}, year = {1993}, mrclass = {73V05 (73K10)}, mrnumber = {MR1260452 (94k:73066)}, mrreviewer = {Luis Trabucho}, abstract = { We report on a new approach to obtaining stable locking free finite element discretizations for the Reissner-Mindlin plate. }, keywords = {Reissner-Mindlin plate, finite element}, subjclass = {65N30,73K10,73K25}, pdf = {http://umn.edu/~arnold/papers/rmelts.pdf} }

@incollection{imacs92, author = {Arnold, Douglas N. and Qin, Jinshui}, title = {Quadratic velocity/linear pressure Stokes elements}, booktitle = {Advances in Computer Methods for Partial Differential Equations--VII}, editor = {Vichnevetsky, R. and Knight, D. and Richter, G.}, publisher = {IMACS}, year = {1992}, pages = {28--34}, abstract = { We study the finite element approximation of the stationary Stokes equations in the velocity-pressure formulation using continuous piecewise quadratic functions for velocity and discontinuous piecewise linear functions for pressure. For some meshes this method is unstable, even after spurious pressure modes are removed. For other meshes there are spurious local pressure modes, but once they are removed the method is stable, and in particular, the velocity converges with optimal order. On yet other meshes there are no spurious pressure modes and the method is stable and optimally convergent for both pressure and velocity. }, keywords = {finite element, Stokes equations}, pdf = {http://umn.edu/~arnold/papers/imacs92.pdf} }

@article{plates, author = {Arnold, Douglas N.}, title = {Innovative finite element methods for plates}, journal = {Mat. Apl. Comput.}, fjournal = {Matem\'atica Aplicada e Computacional}, volume = {10}, year = {1991}, number = {2}, pages = {77--88}, issn = {0101-8205}, mrclass = {65N30 (73K10 73V05)}, mrnumber = {MR1172086 (93c:65136)}, abstract = { Finite element methods for the Reissner-Mindlin plate theory are discussed. Methods in which both the tranverse displacement and the rotation are approximated by finite elements of low degree mostly suffer from locking. However a number of related methods have been devised recently which avoid locking effects. Although the finite element spaces for both the rotation and transverse displacement contain little more than piecewise linear functions, optimal order convergence holds uniformly in the thickness. The main ideas leading to such methods are reviewed and the relationships between various methods are clarified. }, keywords = {Reissner, Mindlin, plate, finite element}, subjclass = {65N30, 73K10}, pdf = {http://umn.edu/~arnold/papers/plates.pdf} }

@article{mixed, author = {Arnold, Douglas N.}, title = {Mixed finite element methods for elliptic problems}, journal = {Comput. Methods Appl. Mech. Engrg.}, fjournal = {Computer Methods in Applied Mechanics and Engineering}, volume = {82}, year = {1990}, number = {1-3}, pages = {281--300}, issn = {0045-7825}, coden = {CMMECC}, mrclass = {65N30 (73V05 76M10)}, mrnumber = {MR1077658 (91h:65168)}, abstract = { This paper treats the basic ideas of mixed finite element methods at an introductory level. Although the viewpoint presented is that of a mathematician, the paper is aimed at practitioners and the mathematical prerequisites are kept to a minimum. A classification of variational principles and of the corresponding weak formulations and Galerkin methods--displacement, equilibrium, and mixed--is given and illustrated through four significant examples. The advantages and disadvantages of mixed methods are discussed. The concepts of convergence, approximability, and stability and their interrelations are developed, and a resume is given of the stability theory which governs the performance of mixed methods. The paper concludes with a survey of techniques that have been developed for the construction of stable mixed methods and numerous examples of such methods. }, keywords = {mixed method, finite element, variational principle}, subjclass = {65N30, 73C35, 73K25}, pdf = {http://umn.edu/~arnold/papers/mixed.pdf}, doi = {10.1016/0045-7825(90)90168-L} }

@article{bdlayer, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {The boundary layer for the {R}eissner--{M}indlin plate model}, journal = {SIAM J. Math. Anal.}, fjournal = {SIAM Journal on Mathematical Analysis}, volume = {21}, year = {1990}, number = {2}, pages = {281--312}, issn = {0036-1410}, coden = {SJMAAH}, mrclass = {73K10 (35B25)}, mrnumber = {MR1038893 (91c:73053)}, mrreviewer = {Fu Ru Jiang}, abstract = { The structure of the solution of the Reissner-Mindlin model of a clamped plate is investigated, emphasizing its dependence on the plate thickness. Asymptotic expansions in powers of the plate thickness are developed for the main physical quantities and the boundary layer is studied. Rigorous error bounds are given for the errors in the expansions in Sobolev norms. As applications, new regularity results for the solutions and new estimates for the difference between the Reissner-Mindlin solution and the solution to the biharmonic equation are derived. Boundary conditions for a clamped edge are considered for most of the paper, and the very similar case of a hard simply-supported plate is discussed briefly at the end. }, keywords = {Reissner, Mindlin, plate, boundary layer}, subjclass = {73K10, 73K25}, pdf = {http://umn.edu/~arnold/papers/bdlayer.pdf}, doi = {10.1137/S0036141093245276} }

@incollection{rmexamples, author = {Arnold, Douglas N. and Falk, Richard}, title = {Edge effects in the {R}eissner-{M}indlin plate theory}, booktitle = {Analytical and Computational Models for Shells}, editor = {Noor, A. and Belytschko, T. and Simo, J.}, publisher = {American Society of Mechanical Engineers}, pages = {71--90}, year = {1989}, abstract = { We investigate the structure of the solution of the Reissner-Mindlin plate equations in its dependence on the plate thickness for various boundary conditions, developing asymptotic expansions in powers of the plate thickness for the main physical quantities. These expansions are uniform up to the boundary for the transverse displacement, but for other variables there is a boundary layer, whose strength depends on the boundary conditions. We give rigorous error bounds for the errors in the expansions in Sobolev norms and make various applications. }, keywords = {Reissner, Mindlin, plate, boundary layer, edge effect}, subjclass = {73K10, 73K25}, pdf = {http://umn.edu/~arnold/papers/rmexamples.pdf} }

@article{heatnote, author = {Arnold, Douglas N. and Noon, Patrick J.}, title = {Coercivity of the single layer heat potential}, journal = {J. Comput. Math.}, fjournal = {Journal of Computational Mathematics. An International Journal on Numerical Methods, Analysis and Applications}, volume = {7}, year = {1989}, number = {2}, pages = {100--104}, issn = {0254-9409}, coden = {JCMMEB}, mrclass = {65P05 (35K05)}, mrnumber = {MR1016829 (90k:65186)}, abstract = { The single layer heat potential operator, K, arises in the solution of initial-boundary value problems for the heat equation using boundary integral methods. In this note we show that K maps a certain anisotropic Sobolev space isomorphically onto its dual, and, moreover, satisfies the coercivity inequality\ge c|q|^2. We thereby establish the well-posedness of the operator equation K q=f and provide a basis for the analysis of the discretizations. }, keywords = {coercivity, heat potential, anisotropic Sobolev space}, pdf = {http://umn.edu/~arnold/papers/heatnote.pdf} }

@article{mindlin, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {A uniformly accurate finite element method for the {R}eissner--{M}indlin plate}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {26}, year = {1989}, number = {6}, pages = {1276--1290}, issn = {0036-1429}, coden = {SJNAAM}, mrclass = {65N30 (73K10 73V05)}, mrnumber = {MR1025088 (91c:65068)}, mrreviewer = {Juhani Pitk{\"a}ranta}, abstract = { We present and analyze a simple finite element method for the Mindlin-Reissner plate model in the primitive variables. Our method uses nonconforming linear finite elements for the transverse displacement and conforming linear finite elements enriched by bubbles for the rotation, with the computation of the element stiffness matrix modified by the inclusion of a simple elementwise averaging. We prove that the method converges with optimal order uniformly with respect to thickness. }, keywords = {Reissner-Mindlin plate, finite element, nonconforming}, subjclass = {65N30,73K10,73K25}, pdf = {http://umn.edu/~arnold/papers/mindlin.pdf}, doi = {10.1137/0726074} }

@article{delta, author = {Cheng, R. S.-C. and Arnold, D. N.}, title = {The delta-trigonometric method using the single-layer potential representation}, journal = {J. Integral Equations Appl.}, fjournal = {Journal of Integral Equations and Applications}, volume = {1}, year = {1988}, number = {4}, pages = {517--547}, issn = {0897-3962}, mrclass = {65N30 (65N35)}, mrnumber = {MR1008405 (90j:65127)}, mrreviewer = {Gunther Schmidt}, abstract = { The Dirichlet problem for Laplace's equation is often solved by means of the single layer potential representation, leading to a Fredholm integral equation of the first kind with logarithmic kernel. We propose to solve this integral equation using a Petrov-Galerkin method with trigonometric polynomials as test functions and, as trial functions, a span of delta distributions centered at boundary points. The approximate solution to the boundary value problem thus computed converges exponentially away from the boundary and algebraically up to the boundary. We show that these convergence results hold even when the discretization matrices are computed via numerical quadratures. Finally, we discuss our implementation of this method using the fast Fourier transform to compute the discretization matrices, and present numerical experiments in order to confirm our theory and to examine the behavior of the method in cases where the theory doesn't apply due to lack of smoothness. }, subjclass = {65R20, 65N30, 65E05, 45L10}, pdf = {http://umn.edu/~arnold/papers/delta.pdf}, doi = {10.1216/JIE-1988-1-4-517} }

@article{corners, author = {Arnold, Douglas N. and Scott, L. Ridgway and Vogelius, Michael}, title = {Regular inversion of the divergence operator with {D}irichlet boundary conditions on a polygon}, journal = {Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4)}, fjournal = {Annali della Scuola Normale Superiore di Pisa. Classe di Scienze. Serie IV}, volume = {15}, year = {1988}, number = {2}, pages = {169--192}, issn = {0391-173X}, coden = {PSNAAI}, mrclass = {35F15 (35B65 35Qxx 47F05)}, mrnumber = {MR1007396 (91i:35043)}, mrreviewer = {Ramiz J. Aliev}, abstract = { We consider the existence of regular solutions to the boundary value problem div U = f on a plane polygonal domain with the Dirichlet boundary condition U=g. We formulate simultaneously necessary and sufficient conditions on f and g in order that a solution U exist in the Sobolev space W^{s+1}_p. In addition to the obvious regularity and integral conditions these consist of at most one compatibility condition at each vertex of the polygon. In the special case of homogeneous boundary data, it is necessary and sufficient that f belong to W^s_p, have mean value zero, and vanish at each vertex. (The latter condition only applies if s is large enough that the point values make sense.) We construct a solution operator which is independent of s and p. As intermediate results we obtain various new trace theorems for Sobolev spaces on polygons. }, keywords = {divergence, trace, Sobolev space}, subjclass = {35B65, 46E35}, pdf = {http://umn.edu/~arnold/papers/corners.pdf} }

@article{elas, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {A new mixed formulation for elasticity}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {53}, year = {1988}, number = {1-2}, pages = {13--30}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {73C60 (65N05 73-08 73C35 73K25)}, mrnumber = {MR946367 (89f:73020)}, mrreviewer = {N. Gass}, abstract = { We propose a new mixed variational formulation for the equations of linear elasticity. It does not require symmetric tensors and consequently is easy to discretize by adapting mixed finite elements developed for scalar second order elliptic equations. }, keywords = {mixed finite element method, elasticity}, subjclass = {65N30, 73C35, 73K25}, pdf = {http://umn.edu/~arnold/papers/elas.pdf}, doi = {10.1007/BF01395876} }

@incollection{heat, author = {Arnold, D. N. and Noon, P. J.}, title = {Boundary integral equations of the first kind for the heat equation}, booktitle = {Boundary elements IX, Vol.\ 3 (Stuttgart, 1987)}, editor = {Brebbia, C. and Wendland, W. and Kuhn, G.}, pages = {213--229}, publisher = {Comput. Mech.}, address = {Southampton}, year = {1987}, mrclass = {80A20 (45B05 80-08)}, mrnumber = {MR965345 (89j:80006)}, mrreviewer = {J. D. P. Donnelly}, abstract = { The solution of the heat equation with Dirichlet boundary conditions by the boundary integral method leads to an integral equation of the first kind to determine the boundary flux. We show that the linear operator so defined is an automorphism from a certain function to another defined on the boundary and is coercive, thereby establishing the well-posedness of the method. }, subjclass = {80A20}, keywords = {boundary integral method, heat equation}, pdf = {http://umn.edu/~arnold/papers/heat.pdf} }

@article{constrained, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {Well-posedness of the fundamental boundary value problems for constrained anisotropic elastic materials}, journal = {Arch. Rational Mech. Anal.}, fjournal = {Archive for Rational Mechanics and Analysis}, volume = {98}, year = {1987}, number = {2}, pages = {143--165}, issn = {0003-9527}, coden = {AVRMAW}, mrclass = {73C15 (73C30)}, mrnumber = {MR866719 (87m:73015)}, mrreviewer = {Stan Chiri{\c{t}}{\u{a}}}, abstract = { We consider the equations of linear homogeneous anisotropic elasticity admitting the possibility that the material is internally constrained, and formulate a simple necessary and sufficient condition for the fundamental boundary value problems to be well-posed. For materials fulfilling the condition, we establish continuous dependence of the displacement and stress on the elastic moduli and ellipticity of the elasticity system. As an application we determine the orthotropic materials for which the fundamental problems are well-posed in terms of their Young's moduli, shear moduli, and Poisson ratios. Finally, we derive a reformulation of the elasticity system that is valid for both constrained and unconstrained materials and involves only one scalar unknown in addition to the displacements. For a two-dimensional constrained material a further reduction to a single scalar equation is outlined. }, keywords = {elasticity, anisotropic, constraint, well-posed}, pdf = {http://umn.edu/~arnold/papers/constrained.pdf}, doi = {10.1007/BF00251231} }

@techreport{orthotropic, author = {Arnold, Douglas N. and Falk, Richard S.}, title = {Continuous dependence on the elastic coefficients for a class of anisotropic materials}, number = {165}, institution = {Institute for Mathematics \& its Applications}, year = {1985}, pages = {1--36}, abstract = { We prove apriori estimates and continuous dependence on the elastic moduli for the equations of homogeneous orthotropic elasticity. These results are uniform with respect to the three Poisson rations, Young's moduli, and shear moduli of the material for certain ranges of these constants. These ranges include the possibility that the compliance tensor is singular such as occurs for incompressible materials. }, keywords = {orthotropic elasticity, incompressible, constrained}, subjclass = {73C30; 73C35}, pdf = {http://umn.edu/~arnold/papers/orthotropic.pdf} }

@article{splinecol, author = {Arnold, Douglas N. and Wendland, Wolfgang L.}, title = {The convergence of spline collocation for strongly elliptic equations on curves}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {47}, year = {1985}, number = {3}, pages = {317--341}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {65N40 (65N30 65N35)}, mrnumber = {MR808553 (87f:65142)}, mrreviewer = {Michael Minkoff}, abstract = { Most boundary element methods for two-dimensional boundary value problems are based on point collocation on the boundary and the use of splines as trial functions. Here we present a unified asymptotic error analysis for even as well as for odd degree splines subordinate to uniform or smoothly graded meshes and prove asymptotic convergence of optimal order. The equations are collocated at the breakpoints for odd degree and the internodal midpoints for even degree splines. The crucial assumption for the generalized boundary integral and integro-differential operators is strong ellipticity. Our analysis is based on simple Fourier expansions. In particular, we extend results by J. Saranen and W.L. Wendland from constant to variable coefficient equations. Our results include the first convergence proof of midpoint collocation with piecewise constant functions, i.e., the panel method for solving systems of Cauchy singular integral equations. }, subjclass = {65R20, 65N99, 65N30, 65E05, 30C30, 73K30, 65N35}, pdf = {http://umn.edu/~arnold/papers/splinecol.pdf}, doi = {10.1007/BF01389582} }

@article{minc, author = {Arnold, D. N. and Brezzi, F.}, title = {Mixed and nonconforming finite element methods: implementation, postprocessing and error estimates}, journal = {RAIRO Mod\'el. Math. Anal. Num\'er.}, fjournal = {RAIRO Mod\'elisation Math\'ematique et Analyse Num\'erique}, volume = {19}, year = {1985}, number = {1}, pages = {7--32}, issn = {0764-583X}, mrclass = {65N30}, mrnumber = {MR813687 (87g:65126)}, abstract = { We discuss a technique of implementing certain mixed finite elements based on the use of Lagrange multipliers to impose interelement continuity. The matrices arising from this implementation are positive definite. Considering some well-known mixed methods, namely the Raviart-Thomas methods for second order elliptic problems and the Hellan-Hermann-Johnson method for biharmonic problems, we show that the computed Lagrange multipliers may be exploited in a simple postprocess to produce better approximation of the original variables. We further extablish an equivalence between the mixed methods and certain modified versions of well-known nonconforming methods, notably the Morley method in the case of the biharmonic problem. The equivalence is exploited to provide error estimates for both the mixed and nonconforming methods. }, keywords = {mixed finite element, Lagrange multiplier}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/minc.pdf}, doi = {10.1051/m2an/1985190100071} }

@article{mini, author = {Arnold, D. N. and Brezzi, F. and Fortin, M.}, title = {A stable finite element for the {S}tokes equations}, journal = {Calcolo}, fjournal = {Calcolo. A Quarterly on Numerical Analysis and Theory of Computation}, volume = {21}, year = {1984}, number = {4}, issn = {0008-0624}, pages = {337--344}, coden = {CDABAE}, mrclass = {65N30 (76-08 76D07)}, mrnumber = {MR799997 (86m:65136)}, mrreviewer = {Gheorghe Micula}, abstract = { We present in this paper a new velocity-pressure finite element for the computation of Stokes flow. We discretize the velocity field iwth continuous piecewise linear functions enriched by bubble functions, and the pressure by piecewise linear functions. We show that this element satisfies the usual inf-sup condition and converges with first order for both velocities and pressures. Finally we relate this element to families of higher order elements and to the popular Taylor-Hood element. }, keywords = {finite element, Stokes equations}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/mini.pdf}, doi = {10.1007/BF02576171} }

@incollection{newmixed, author = {Arnold, Douglas N.}, title = {A new mixed formulation for the numerical solution of elasticity problems}, booktitle = {Advances in Computer Methods for Partial Differential Equations--V}, editor = {Vichnevetsky, R. and Stepleman, R.}, publisher = {IMACS}, year = {1984}, pages = {353--356}, abstract = { A mixed formulation for boundary value problems in linear elastostatics is presented. This formulation differs slightly from the classical Hellinger-Reissner formulation. The unknown fields are the displacement and a tensor related but not equal to the stress. The tensors appearing in the formulation need not be symmetric, and consequently mixed finite elements developed for scalar second order elliptic problems may be applied directly. }, keywords = {elasticity, mixed finite element}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/newmixed.pdf} }

@incollection{singint, author = {Arnold, Douglas N.}, title = {The effect of the test functions on the convergence of spline projection methods for singular integral equations}, booktitle = {Numerical Solution of Singular Integral Equations}, editor = {Gerasoulis, A. and Vichnevetsky, R.}, publisher = {IMACS}, year = {1984}, pages = {1--4}, abstract = { We investigate the asymptotic convergence properties of a variety of methods for the numerical solution of the system of singular integral equations arising from the traction problem of plane elasticity. Various sorts of Galerkin methods and collocation methods are considered, all of which determine a spline approximation via paring with certain test functions; the test functions may be splines of the same degree as the trial functions (ordinary Galerkin methods), splines of different degree (Petrov-Galerkin methods), delta functions (collocation), or trigonometric polynomials (spline-trig methods). The choice of test functions is shown to have a significant influence on the convergence properties. }, keywords = {integral equations, Galerkin methods, collocation}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/singint.pdf} }

@article{peers, author = {Arnold, Douglas N. and Brezzi, Franco and Douglas, Jr., Jim}, title = {P{EERS}: a new mixed finite element for plane elasticity}, journal = {Japan J. Appl. Math.}, fjournal = {Japan Journal of Applied Mathematics}, volume = {1}, year = {1984}, number = {2}, pages = {347--367}, issn = {0910-2043}, mrclass = {65N30 (73-08)}, mrnumber = {MR840802 (87h:65189)}, mrreviewer = {C. Dierieck}, abstract = { A mixed finite element procedure for plane elasticity is introduced and analyzed. The symmetry of the stress tensor is enforced through the introduction of a Lagrange multiplier. An additional Lagrange multiplier is instroduced to simplify the algebraic system. Applications are made to incompressible elastic problems and to plasticity problems. }, keywords = {finite element methods, plane elasticity}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/peers.pdf}, doi = {10.1007/BF03167064} }

@article{hopes, author = {Arnold, Douglas N. and Douglas, Jr., Jim and Gupta, Chaitan P.}, title = {A family of higher order mixed finite element methods for plane elasticity}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {45}, year = {1984}, number = {1}, pages = {1--22}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {65N30 (73-08)}, mrnumber = {MR761879 (86a:65112)}, mrreviewer = {Reinhard Scholz}, abstract = { The Dirichler problem for the equations of plane elasticity is approximated by a mixed finite element method using a new family of composite finite elements having properties analogous to those possessed by the Raviart-Thomas mixed finite elements for a scalar, second-order elliptic equation. Estimates of optimal order and minimal regularity are derived for the errors in the displacement vector and the stress tensor in L^2 and optimal order negative norm estimates are obtained in (H^s)' for a range of s depending on the index of the finite element space. An optimal order estimate inL in L_infinity for the displacement error is given. Also, a quasioptimal estimate is derived in an appropriate space. All estimates are valid uniformly with respect to the compressibility and apply in the incompressible case. The formulation of the elements is presented in detail. }, keywords = {finite element methods, plane elasticity}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/hopes.pdf}, doi = {10.1007/BF01379659} }

@article{colpde, author = {Arnold, Douglas N. and Saranen, Jukka}, title = {On the asymptotic convergence of spline collocation methods for partial differential equations}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {21}, year = {1984}, number = {3}, pages = {459--472}, issn = {0036-1429}, coden = {SJNAAM}, mrclass = {65N35}, mrnumber = {MR744168 (86g:65221)}, mrreviewer = {Ian Gladwell}, abstract = { We examine the asymptotic accuracy of the method of collocation for the approximate solution of linear elliptic partial differential equations. Specifically we consider the nodal collocation of a second order equation in the plane with biperiodicity conditions using tensor product smooth splines of odd degree as trial functions. We prove optimal rates of convergence in L2 for partial derivatives of the approximate solution which are of order at least two in one variable, while the solution itself and its gradient converge in L2 at rates less than the optimal approximation theoretic results. }, pdf = {http://umn.edu/~arnold/papers/colpde.pdf}, doi = {10.1137/0721034} }

@article{principles, author = {Arnold, D. N. and Babu{\v{s}}ka, I. and Osborn, J.}, title = {Finite element methods: principles for their selection}, journal = {Comput. Methods Appl. Mech. Engrg.}, fjournal = {Computer Methods in Applied Mechanics and Engineering}, volume = {45}, year = {1984}, number = {1-3}, pages = {57--96}, issn = {0045-7825}, coden = {CMMECC}, mrclass = {65N30 (65M60 73-08 76-08)}, mrnumber = {MR759804 (86f:65187)}, abstract = { Principles for the selection of a finite element method for a particular problem are discusses. These principles are stated in terms of the notion of approximability, optimality, and stability. Several examples are discussed in details as illustrations. Conclusions regarding the selection of finite element methods are summarized in the final section of the paper. }, keywords = {finite elements, aproximability, stability, optimality}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/principles.pdf}, doi = {10.1016/0045-7825(84)90151-8} }

@article{collocation, author = {Arnold, Douglas N. and Wendland, Wolfgang L.}, title = {On the asymptotic convergence of collocation methods}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {41}, year = {1983}, number = {164}, pages = {349--381}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65N35 (65L10 65L60 65N30)}, mrnumber = {MR717691 (85h:65254)}, abstract = { We prove quasioptimal and optimal order estimates in various Sobolev norms for the approximation of linear strongly elliptic pseudodifferential equations in one independent variable by the method of nodal collocation by odd degree polynomial splines. The analysis pertains in particular to many of the boundary element methods used for numerical computation in engineering applications. Equations to which the analysis is applied include Fredholm integral equations of the second kind, certain first kind Fredholm equations, singular integral equations involving Cauchy kernels, a variety of integro-differential equations, and two-point boundary value problems for ordinary differential equations. The error analysis is based on an equivalence which we establish between the collocation methods and certain nonstandard Galerkin methods. We compare the collocation method with a standard Galerkin method using splines of the same degree, showing that the Galerkin method is quasioptimal in a Sobolev space of lower index and furnishes optimal order approximation for a range of Sobolev indices containing and extending below that for the collocation method, and so the standard Galerkin method achieves higher rates of convergence. }, keywords = {collocation, spline, integral equation}, subjclass = {65R20, 65R99, 65L10, 65N99, 45J05, 45L10, 45F15, 35S99, 30C30, 73K30, 31A30}, pdf = {http://umn.edu/~arnold/papers/collocation.pdf}, doi = {10.1090/S0025-5718-1983-0717691-6} }

@article{splinetrig, author = {Arnold, Douglas N.}, title = {A spline-trigonometric {G}alerkin method and an exponentially convergent boundary integral method}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {41}, year = {1983}, number = {164}, pages = {383--397}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65N15 (41A15 45L05 65D07 65R20)}, mrnumber = {MR717692 (84m:65117)}, mrreviewer = {G. T. Symm}, abstract = { We consider a Galerkin method for functional equations in one space variable which uses periodic cardinal splines as trial functions and trigonometric polynomials as test functions. We analyze the method applied to the integral equation of the first kind arising from a single layer potential formulation of the Dirichlet problem in the interior or exterior of an analytic plane curve. In contrast to ordinary spline Galerkin methods, we show that the method is stable, and so provides quasioptimal approximation, in a large family of Hilbert spaces including all the Sobolev spaces of negative order. As a consequence we prove that the approximate solution to the Dirichlet problem and all its derivatives converge pointwise with exponential rate. }, keywords = {spline, spline-trigonometric, Galerkin method, boundary integral}, subjclass = {65R20, 65N30, 65E05, 45L10, 45B05, 41A15}, pdf = {http://umn.edu/~arnold/papers/splinetrig.pdf}, doi = {10.1090/S0025-5718-1983-0717692-8} }

@incollection{selection, author = {Arnold, D. N. and Babu{\v{s}}ka, I. and Osborn, J.}, title = {Selection of finite element methods}, booktitle = {Hybrid and mixed finite element methods (Atlanta, Ga., 1981)}, editor = {Atluri, S. and Gallagher, R. and Zienkiewicz, O.}, series = {Wiley-Intersci. Pub.}, pages = {433--451}, publisher = {Wiley}, address = {New York}, year = {1983}, mrclass = {65N30 (73K25)}, mrnumber = {MR718351 (84m:65120)}, abstract = { The goal of engineering computations is to obtain quantitative information about engineering problems. This goal is usually achieved by the approximation solution of a mathematically formulated problem. Although a relevant mathematical formulation of the problem and its approximation solution are closely related, here we shall suppose that a mathematical formulation has already been determined and is amenable to an approximate treatment. We shall discuss a broad class of approaches based on variational methods of discretization which allow one to find the approximation solution within a desired range of accuracy. We discuss properties of these methods which enable us to distinguish among them and which aid in the selection or design of a method which is effective in achieving the given goals of the computation. }, keywords = {finite elements, variational methods}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/selection.pdf} }

@incollection{colvsgal, author = {Arnold, D. N. and Wendland, W. L.}, title = {Collocation versus {G}alerkin procedures for boundary integral methods}, booktitle = {Boundary element methods in engineering (Southampton, 1982)}, editor = {Brebbia, C.}, pages = {18--33}, publisher = {Springer}, address = {Berlin}, year = {1982}, mrclass = {73K30 (65N30 65N35)}, mrnumber = {MR737197 (85f:73074)}, abstract = { We compare the efficiency of the solution of two-dimensional elliptic boundary value problems via boundary integral methods using two different discretization procedures with comparable convergence rates: Galerkin procedures with numerical integration and collocation. }, keywords = {boundary integral methods, boundary element methods,}, pdf = {http://umn.edu/~arnold/papers/colvsgal.pdf} }

@article{intpen, author = {Arnold, Douglas N.}, title = {An interior penalty finite element method with discontinuous elements}, journal = {SIAM J. Numer. Anal.}, fjournal = {SIAM Journal on Numerical Analysis}, volume = {19}, year = {1982}, number = {4}, pages = {742--760}, issn = {0036-1429}, coden = {SJNAAM}, mrclass = {65N30}, mrnumber = {MR664882 (83f:65173)}, abstract = { A new semidiscrete finite element method for the solution of second order nonlinear parabolic boundary value problems is formulated and analyzed. The test and trial spaces consist of discontinuous piecewise polynomial functions over quite general meshes with interelement continuity enforced approximately by means of penalties. Optimal order error estimates in energy and L2-norms are stated in terms of locally expressed quantities. They are proved first for a model problem and then in general. }, pdf = {http://umn.edu/~arnold/papers/intpen.pdf}, doi = {10.1137/0719052} }

@article{kdv, author = {Arnold, Douglas N. and Winther, Ragnar}, title = {A superconvergent finite element method for the {K}orteweg-de~{V}ries equation}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {38}, year = {1982}, number = {157}, pages = {23--36}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65M60 (76A60 76B15)}, mrnumber = {MR637284 (82m:65087)}, abstract = { An unconditionally stable fully discrete finite element method for the Korteweg-de Vries equation is presented. In addition to satisfying optimal order global estimates, it is shown that this method is superconvergent at the nodes. The algorithm is derived from the conservative method proposed by the second author by the introduction of a small time-independent forcing term into the discrete equations. This term is a form of the quasiprojection which was first employed in the analysis of superconvergence phenomena for parabolic problems. However, in the present work, unlike in the parabolic case, the quasiprojection is used as perturbation of the discrete equations and does not affect the choice of initial values. }, keywords = {superconvergence, finite element, Korteweg-de Vries equation}, subjclass = {65N30, 65N15, 76B15}, pdf = {http://umn.edu/~arnold/papers/kdv.pdf}, doi = {10.1090/S0025-5718-1982-0637284-8} }

@article{paramdep, author = {Arnold, Douglas N.}, title = {Discretization by finite elements of a model parameter dependent problem}, journal = {Numer. Math.}, fjournal = {Numerische Mathematik}, volume = {37}, year = {1981}, number = {3}, pages = {405--421}, issn = {0029-599X}, coden = {NUMMA7}, mrclass = {65N30 (73K12)}, mrnumber = {MR627113 (82h:65077)}, abstract = { The discretization by finite elements of a model variational problem for a clamped loaded beam is studied with emphasis on the effect of the beam thickness, which appears as a parameter in the problem, on the accuracy. It is shown that the approximation achieved by a standard finite element method degenerates for thin beams. In contrast a large family of mixed finite element methods are shown to yield quasioptimal approximation independent of the thickness parameter. The most useful of these methods may be realized by replacing the integrals appearing in the stiffness matrix of the standard method by Gauss quadratures. }, keywords = {mixed finite element method, reduced integration, Timoshenko beam, parameter}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/paramdep.pdf}, doi = {10.1007/BF01400318} }

@incollection{robustness, author = {Arnold, Douglas N.}, title = {Robustness of finite element methods for a model parameter dependent problem}, booktitle = {Advances in Computer Methods for Partial Differential Equations--IV}, editor = {Vichnevetsky, R. and Stepleman, R.}, publisher = {IMACS}, year = {1981}, pages = {18--22}, abstract = { A convergence analysis is presented for standard and mixed finite element discretizations of a model system of equations for a transversely loaded beam. The equations depend parametrically on the beam thickness and the emphasis of the analysis is on the robustness of the methods with respect to this parameter. The mixed methods are shown to be far more robust than the standard methods employing elements of the same degree. Moveover they entail no additional computational expense. Computational results are included to illustrate the main results. }, keywords = {beam, mixed finite element, parameter-dependence}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/robustness.pdf} }

@article{soboleveqn, author = {Arnold, Douglas N. and Douglas, Jr., Jim and Thom{\'e}e, Vidar}, title = {Superconvergence of a finite element approximation to the solution of a {S}obolev equation in a single space variable}, journal = {Math. Comp.}, fjournal = {Mathematics of Computation}, volume = {36}, year = {1981}, number = {153}, pages = {53--63}, issn = {0025-5718}, coden = {MCMPAF}, mrclass = {65N30 (35K70)}, mrnumber = {MR595041 (82f:65108)}, mrreviewer = {J. D. Atkinson}, abstract = { A standard Galerkin method for a quasilinear equation of Sobolev type using continuous, piecewise-polynomial spaces is presented and analyzed. Optimal order error estimates are established in various norms, and nodal superconvergence is demonstrated. Discretization in time by explicit single-step methods is discussed. }, keywords = {superconvergence, finite element, Sobolev equation}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/soboleveqn.pdf}, doi = {10.1090/S0025-5718-1981-0595041-4} }

@article{supercon, author = {Arnold, D. N. and Douglas, Jr., J.}, title = {Superconvergence of the {G}alerkin approximation of a quasilinear parabolic equation in a single space variable}, journal = {Calcolo}, fjournal = {Calcolo}, volume = {16}, year = {1979}, number = {4}, pages = {345--369}, issn = {0008-0624}, coden = {CDABAE}, mrclass = {65N30}, mrnumber = {MR592476 (83b:65118)}, mrreviewer = {Rudolf Gorenflo}, abstract = { The asympotic expansion of the Galerkin solution of a parabolic equation by means of a sequence of elliptic projections that was introduced by Douglas, Dupont, and Wheeler is carried out for a quasilinear equation. This quasi-projection can be applied to establish knot superconvergence in the case of a single space variable. In addition, an optimal order error estimate in L-infinity(L-infinity) is derived for a single space variable. }, keywords = {finite element, parabolic equation, superconvergence}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/supercon.pdf}, doi = {10.1007/BF02576636} }

@phdthesis{thesis, author = {Arnold, D. N.}, title = {An interior penalty finite element method with discontinuous elements}, school = {University of Chicago}, year = {1979}, abstract = { A nonconforming finite element procedure for the solution of second order, nonlinear parabolic boundary value problems is formulated and analyzed. The finite element space consists of discontinuous piecewise polynomial functions over quite general meshes, with inter-element continuity being enforced approximately by means of penalties. Optimal order error estimates in energy and L2 norms are stated in terms of locally expressed quantities. They are proved first for a model problem and then in general. }, keywords = {finite element, interior penalty}, subjclass = {65N30}, pdf = {http://umn.edu/~arnold/papers/ArnoldThesis.pdf} }

*This file was generated by
bibtex2html 1.99.*