books.bib

@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.},
}

publist.bib

@article{landscapelaw0,
  author = {Arnold, Douglas N. and Filoche, Marcel and Mayboroda, Svitlana and Wang, Wei and Zhang, Shiwen},
  title = {The landscape law for tight binding {H}amiltonians},
  journal = {Comm. Math. Phys.},
  fjournal = {Communications in Mathematical Physics},
  volume = {396 },
  number = {},
  pages = {1339-1391},
  year = {2022},
  abstract = {The present paper extends the landscape theory pioneered in Filoche and
Mayboroda (Proc. Natl. Acad. Sci. USA 109(37):14761-14766, 2012), Arnold et al.
(Commun. Partial Differ. Equ. 44(11):1186-1216, 2019) and David et al. (Adv Math
390:107946, 2021) to the tight-binding Schr\"odinger operator on ${\mathbb{Z}^d}$. In particular, we
establish upper and lower bounds for the integrated density of states in terms of the
counting function based upon the localization landscape.
},
  keywords = {Integrated Density of States, Anderson models, tight-binding model, Landscape Law},
  skipsubjclass = {35J10, 35P20, 82B44},
  pdf = {https://umn.edu/~arnold/papers/landscapelaw0.pdf},
  doi = {10.1007/s00220-022-04494-8}
}
@article{spectral-functions,
  author = {Pelletier, Pierre and Delande, Dominique and Josse, Vincent and Aspect, Alain and Mayboroda, Svitlana and Arnold, Douglas N. and Filoche, Marcel},
  title = {Spectral functions and localization-landscape theory in speckle potentials},
  journal = {Phys. Rev. A},
  fjournal = {Physical Review A},
  volume = {105},
  number = {},
  pages = {023314},
  year = {2022},
  skipnote = {To appear. arXiv preprint 2111.13155},
  abstract = {Spectral function is a key tool for understanding the behavior
of Bose-Einstein condensates of cold atoms in random potentials generated
by a laser speckle. In this paper we introduce a new method for computing
the spectral functions in disordered potentials. Using a combination of the
Wigner-Weyl approach with the landscape theory, we build an approximation
for the Wigner distributions of the eigenstates in the phase space and
show its accuracy in all regimes, from the deep quantum regime to the
intermediate and semiclassical. Based on this approximation, we devise a
method to compute the spectral functions using only the landscape-based
effective potential. The paper demonstrates the efficiency of the proposed
approach for disordered potentials with various statistical properties
without requiring any adjustable parameters.
},
  keywords = {spectral function, localization landscape, speckle potential},
  skipsubjclass = {},
  pdf = {https://umn.edu/~arnold/papers/spectral-functions.pdf},
  doi = {10.1103/PhysRevA.105.023314}
}
@article{l2bddproj,
  author = {Arnold, Douglas N. and G\'uzman, Johnny},
  title = {Local {$L^2$-bounded} commuting projections in {FEEC}},
  volume = {55},
  number = {5},
  pages = {2169-2184},
  year = {2021},
  journal = {ESAIM: M2AN},
  fjournal = {ESAIM: Mathematical Modelling and Numerical Analysis},
  skipnote = {arXiv preprint 2104.00184},
  abstract = {We construct local projections into canonical finite element spaces that appear
in the finite element exterior calculus.  These projections are bounded in $L^2$ and commute
with the exterior derivative.
},
  keywords = {cochain projection, commuting projection, finite element exterior calculus},
  skipsubjclass = {65N30},
  pdf = {https://umn.edu/~arnold/papers/l2bddproj.pdf},
  doi = {10.1051/m2an/2021054}
}
@article{bgg,
  author = {Arnold, Douglas N. and Hu, Kaibo},
  title = {Complexes from complexes},
  year = {2021},
  journal = {Found. Comput. Math.},
  volume = {21},
  pages = {1739-1774},
  skipnote = {arXiv preprint 2005.12437},
  abstract = {This paper is concerned with the derivation and properties of
differential complexes arising from a variety of problems in differential
equations, with applications in continuum mechanics, relativity, and
other fields.  We present a systematic procedure which, starting from
well-understood differential complexes such as the de Rham complex,
derives new complexes and deduces the properties of the new complexes
from the old.  We relate the cohomology of the derived complex to that of
the input complexes and show that the new complex has closed ranges, and,
consequently, satisfies a Hodge decomposition, Poincar\'e type inequalities,
well-posed Hodge-Laplacian boundary value problems, regular decomposition,
and compactness properties on general Lipschitz domains.
},
  keywords = {Differential complex, Hilbert complex, de Rham complex, elasticity complex, BGG resolution, finite element exterior calculus},
  subjclass = {35J10, 65N30},
  pdf = {https://umn.edu/~arnold/papers/bgg.pdf},
  doi = {10.1007/s10208-021-09498-9}
}
@article{landscapelaw1,
  author = {Desforges, Perceval and Mayboroda, Svitlana and Zhang, Shiwen and David, Guy and Arnold, Douglas N. and Wang, Wei and Filoche, Marcel},
  title = {Sharp estimates for the integrated density of states in {A}nderson tight-binding models},
  year = {2021},
  journal = {Phys. Rev. A},
  fjournal = {Physical Review A},
  volume = {104},
  number = {1},
  pages = {012207},
  numpages = {11},
  abstract = {Recent work
has proved the existence of bounds from above and below for the integrated
density of states (IDOS) of the Schr\"odinger operator throughout the
spectrum, called the landscape law. These bounds involve dimensional
constants whose optimal values are yet to be determined. Here, we
investigate the accuracy of the landscape law in 1D and 2D tight-binding
Anderson models, with binary or uniform random distributions. We show, in
particular, that in 1D, the IDOS can be approximated with high accuracy
through a single formula involving a remarkably simple multiplicative
energy shift. In 2D, the same idea applies but the prefactor has to be
changed between the bottom and top parts of the spectrum.
},
  keywords = {Integrated Density of States, Anderson models, tight-binding model, Landscape Law},
  skipsubjclass = {35J10, 65N30},
  pdf = {https://umn.edu/~arnold/papers/landscapelaw1.pdf},
  doi = {10.1103/PhysRevA.104.012207}
}
@article{localization-reply,
  author = {Filoche, M. and Arnold, D. and David, G. and Jerison, D. and Mayboroda, S.},
  title = {Reply to comment on ``{E}ffective Confining Potential of Quantum States in Disordered Media''},
  journal = {Phys. Rev. Lett.},
  volume = {124},
  pages = {219702},
  numpages = {2},
  year = {2020},
  month = {May},
  doi = {10.1103/PhysRevLett.124.219702},
  pdf = {https://umn.edu/~arnold/papers/localization-reply.pdf}
}
@article{curved_hhj,
  author = {Arnold, Douglas N. and Walker, Shawn W.},
  title = {The {Hellan-Herrmann-Johnson} method with curved elements},
  journal = {SIAM J. Numer. Anal.},
  fjournal = {SIAM Journal on Numerical Analysis},
  volume = {58},
  number = {5},
  pages = {2829-2855},
  year = {2020},
  skipnote = {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 = {https://umn.edu/~arnold/papers/curved_hhj.pdf},
  doi = {10.1137/19M1288723}
}
@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 = {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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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},
  number = {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 = {https://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},
  number = {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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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},
  number = {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 = {https://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},
  number = {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 = {https://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},
  number = {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 = {https://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 finite element spaces of differential forms which satisfy
the crucial assumptions of the finite 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 define and understand the
spaces. These tools allow us to treat a wide variety of situations, which are often treated separately,
in a unified fashion.},
  keywords = {finite element differential form, finite element exterior calculus},
  subjclass = {65N30},
  pdf = {https://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 = {https://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 finite 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 finite element method, in which
the rotation of the solution is introduced as a second unknown, is
advantageous, and appropriate choices of mixed finite 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 finite 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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://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 = {https://umn.edu/~arnold/papers/ArnoldThesis.pdf}
}

This file was generated by bibtex2html 1.99.