BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michael Holst

Michael Holst appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APgr-qcmath-phmath.MPmath.DGmath.FAmath.OCphysics.comp-phNumerical Analysisphysics.data-anphysics.hist-phphysics.ins-det

Catalog footprint

What is connected

38works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2025arXiv

Towards Real Time Control of Water Engineering with Nonlinear Hyperbolic Partial Differential Equations

This paper examines aspirational requirements for software addressing mixed-integer optimization problems constrained by the nonlinear Shallow Water partial differential equations (PDEs), motivated by applications such as river-flow management in hydropower cascades. Realistic deployment of such software would require the simultaneous treatment of nonlinear and potentially non-smooth PDE dynamics, limited theoretical guarantees on the existence and regularity of control-to-state mappings under varying boundary conditions, and computational performance compatible with operational decision-making. In addition, practical settings motivate consideration of uncertainty arising from forecasts of demand, inflows, and environmental conditions. At present, the theoretical foundations, numerical optimization methods, and large-scale scientific computing tools required to address these challenges in a unified and tractable manner remain the subject of ongoing research across the associated research communities. Rather than proposing a complete solution, this work uses the problem as a case study to identify and organize the mathematical, algorithmic, and computational components that would be necessary for its realization. The resulting framework highlights open challenges and intermediate research directions, and may inform both more circumscribed related problems and the design of future large-scale collaborative efforts aimed at addressing such objectives.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Note On Determining Projections for Non-Homogeneous Incompressible Fluids

In this note, we consider a viscous incompressible fluid in a finite domain in both two and three dimensions, and examine the question of determining degrees of freedom (projections, functionals, and nodes). Our particular interest is the case of non-constant viscosity, representing either a fluid with viscosity that changes over time (such as an oil that loses viscosity as it degrades), or a fluid with viscosity varying spatially (as in the case of two-phase or multi-phase fluid models). Our goal is to apply the determining projection framework developed by the second author in previous work for weak solutions to the Navier-Stokes equations, in order to establish bounds on the number of determining functionals for this case, or equivalently, the dimension of a determining set, based on the approximation properties of an underlying determining projection. The results for the case of time-varying viscosity mirror those for weak solutions established in earlier work for constant viscosity. The case of space-varying viscosity, treated within a single-fluid Navier-Stokes model, is quite challenging to analyze, but we explore some preliminary ideas for understanding this case.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Symmetry Breaking in Density Functional Theory due to Dirac Exchange for a Hydrogen Molecule

We study symmetry breaking in the mean field solutions to the 2 electron hydrogen molecule within Kohn Sham (KS) local spin density function theory with Dirac exchange (the XLDA model). This simplified model shows behavior related to that of the (KS) spin density functional theory (SDFT) predictions in condensed and molecular systems. The Kohn Sham solutions to the constrained SDFT variation problem undergo spontaneous symmetry breaking as the relative strength of the non-convex exchange term increases. This results in the change of the molecular ground state from a paramagnetic state to an antiferromagnetic ground states and a stationary symmetric delocalized 1st excited state. We further characterize the limiting behavior of the minimizer when the strength of the exchange term goes to infinity. This leads to further bifurcations and highly localized states with varying character. The stability of the various solution classes is demonstrated by Hessian analysis. Finite element numerical results provide support for the formal conjectures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite Element Exterior Calculus for Evolution Problems

Arnold, Falk, and Winther [Bull. Amer. Math. Soc. 47 (2010), 281--354] showed that mixed variational problems, and their numerical approximation by mixed methods, could be most completely understood using the ideas and tools of Hilbert complexes. This led to the development of the Finite Element Exterior Calculus (FEEC) for a large class of linear elliptic problems. More recently, Holst and Stern [Found. Comp. Math. 12:3 (2012), 263--293 and 363--387] extended the FEEC framework to semi-linear problems, and to problems containing variational crimes, allowing for the analysis and numerical approximation of linear and nonlinear geometric elliptic partial differential equations on Riemannian manifolds of arbitrary spatial dimension, generalizing surface finite element approximation theory. In this article, we develop another distinct extension to the FEEC, namely to parabolic and hyperbolic evolution systems, allowing for the treatment of geometric and other evolution problems. Our approach is to combine the recent work on the FEEC for elliptic problems with a classical approach to solving evolution problems via semi-discrete finite element methods, by viewing solutions to the evolution problem as lying in time-parameterized Hilbert spaces (or Bochner spaces). Building on classical approaches by Thomee for parabolic problems and Geveci for hyperbolic problems, we establish a priori error estimates for Galerkin FEM approximation in the natural parametrized Hilbert space norms. In particular, we recover the results of Thomee and Geveci for two-dimensional domains and lowest-order mixed methods as special cases, effectively extending their results to arbitrary spatial dimension and to an entire family of mixed methods. We also show how the Holst and Stern framework allows for extensions of these results to certain semi-linear evolution problems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Emergence of Gravitational Wave Science: 100 Years of Development of Mathematical Theory, Detectors, Numerical Algorithms, and Data Analysis Tools

On September 14, 2015, the newly upgraded Laser Interferometer Gravitational-wave Observatory (LIGO) recorded a loud gravitational-wave (GW) signal, emitted a billion light-years away by a coalescing binary of two stellar-mass black holes. The detection was announced in February 2016, in time for the hundredth anniversary of Einstein's prediction of GWs within the theory of general relativity (GR). The signal represents the first direct detection of GWs, the first observation of a black-hole binary, and the first test of GR in its strong-field, high-velocity, nonlinear regime. In the remainder of its first observing run, LIGO observed two more signals from black-hole binaries, one moderately loud, another at the boundary of statistical significance. The detections mark the end of a decades-long quest, and the beginning of GW astronomy: finally, we are able to probe the unseen, electromagnetically dark Universe by listening to it. In this article, we present a short historical overview of GW science: this young discipline combines GR, arguably the crowning achievement of classical physics, with record-setting, ultra-low-noise laser interferometry, and with some of the most powerful developments in the theory of differential geometry, partial differential equations, high-performance computation, numerical analysis, signal processing, statistical inference, and data science. Our emphasis is on the synergy between these disciplines, and how mathematics, broadly understood, has historically played, and continues to play, a crucial role in the development of GW science. We focus on black holes, which are very pure mathematical solutions of Einstein's gravitational-field equations that are nevertheless realized in Nature, and that provided the first observed signals.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite Element Exterior Calculus for Parabolic Evolution Problems On Riemannian Hypersurfaces

Over the last ten years, the Finite Element Exterior Calculus (FEEC) has been developed as a general framework for linear mixed variational problems, their numerical approximation by mixed methods, and their error analysis. The basic approach in FEEC, pioneered by Arnold, Falk, and Winther in two seminal articles in 2006 and 2010, interprets these problems in the setting of Hilbert complexes, leading to a more general and complete understanding. Over the last five years, the FEEC framework has been extended to a broader set of problems. One such extension, due to Holst and Stern in 2012, was to problems with variational crimes, allowing for the analysis and numerical approximation of linear and geometric elliptic partial differential equations on Riemannian manifolds of arbitrary spatial dimension. Their results substantially generalize the existing surface finite element approximation theory in several respects. In 2014, Gillette, Holst, and Zhu extended the FEEC in another direction, namely to parabolic and hyperbolic evolution systems by combining the FEEC framework for elliptic operators with classical approaches for parabolic and hyperbolic operators, by viewing solutions to the evolution problem as lying in Bochner spaces (spaces of Banach-space valued parametrized curves). Related work on developing an FEEC theory for parabolic evolution problems has also been done independently by Arnold and Chen. In this article, we extend the work of Gillette-Holst-Zhu and Arnold-Chen to evolution problems on Riemannian manifolds, through the use of framework developed by Holst and Stern for analyzing variational crimes. We establish a priori error estimates that reduce to the results from earlier work in the flat (non-criminal) setting. Some numerical examples are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence and Optimality of Adaptive Mixed Methods on Surfaces

In a 1988 article, Dziuk introduced a nodal finite element method for the Laplace-Beltrami equation on 2-surfaces approximated by a piecewise-linear triangulation, initiating a line of research into surface finite element methods (SFEM). Demlow and Dziuk built on the original results, introducing an adaptive method for problems on 2-surfaces, and Demlow later extended the a priori theory to 3-surfaces and higher order elements. In a separate line of research, the Finite Element Exterior Calculus (FEEC) framework has been developed over the last decade by Arnold, Falk and Winther and others as a way to exploit the observation that mixed variational problems can be posed on a Hilbert complex, and Galerkin-type mixed methods can be obtained by solving finite dimensional subproblems. In 2011, Holst and Stern merged these two lines of research by developing a framework for variational crimes in abstract Hilbert complexes, allowing for application of the FEEC framework to problems that violate the subcomplex assumption of Arnold, Falk and Winther. When applied to Euclidean hypersurfaces, this new framework recovers the original a priori results and extends the theory to problems posed on surfaces of arbitrary dimensions. In yet another seemingly distinct line of research, Holst, Mihalik and Szypowski developed a convergence theory for a specific class of adaptive problems in the FEEC framework. Here, we bring these ideas together, showing convergence and optimality of an adaptive finite element method for the mixed formulation of the Hodge Laplacian on hypersurfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence of Goal-Oriented Adaptive Finite Element Methods for Semilinear Problems

In this article we develop a convergence theory for goal-oriented adaptive finite element algorithms designed for a class of second-order semilinear elliptic equations. We briefly discuss the target problem class, and introduce several related approximate dual problems that are crucial to both the analysis as well as to the development of a practical numerical method. We then review some standard facts concerning conforming finite element discretization and error-estimate-driven adaptive finite element methods (AFEM). We include a brief summary of a priori estimates for this class of semilinear problems, and then describe some goal-oriented variations of the standard approach to AFEM (GOAFEM). Following the recent approach of Mommer-Stevenson and Holst-Pollock for increasingly general linear problems, we first establish a quasi-error contraction result for the primal problem. We then develop some additional estimates that make it possible to establish contraction of the combined primal-dual quasi-error, and subsequently show convergence with respect to the quantity of interest. Finally, a sequence of numerical experiments are then carefully examined. It is observed that the behavior of the implementation follows the predictions of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-CMC Solutions to the Einstein Constraint Equations on Asymptotically Euclidean Manifolds with Apparent Horizon Boundaries

In this article we further develop the solution theory for the Einstein constraint equations on an n-dimensional, asymptotically Euclidean manifold M with interior boundary S. Building on recent results for both the asymptotically Euclidean and compact with boundary settings, we show existence of far-from-CMC and near-CMC solutions to the conformal formulation of the Einstein constraints when nonlinear Robin boundary conditions are imposed on S, similar to those analyzed previously by Dain (2004), by Maxwell (2004, 2005), and by Holst and Tsogtgerel (2013) as a model of black holes in various CMC settings, and by Holst, Meier, and Tsogtgerel (2013) in the setting of far-from-CMC solutions on compact manifolds with boundary. These "marginally trapped surface" Robin conditions ensure that the expansion scalars along null geodesics perpendicular to the boundary region S are non-positive, which is considered the correct mathematical model for black holes in the context of the Einstein constraint equations. Assuming a suitable form of weak cosmic censorship, the results presented in this article guarantee the existence of initial data that will evolve into a space-time containing an arbitrary number of black holes. A particularly important feature of our results are the minimal restrictions we place on the mean curvature, giving both near- and far-from-CMC results that are new.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Alternative Between Non-unique and Negative Yamabe Solutions to the Conformal Formulation of the Einstein Constraint Equations

The conformal method has been effective for parametrizing solutions to the Einstein constraint equations on closed 3-manifolds. However, it is still not well-understood; for example, existence of solutions to the conformal equations for zero or negative Yamabe metrics is still unknown without the so-called ``CMC'' or ``near-CMC'' assumptions. The first existence results without such assumptions, termed the ``far-from-CMC'' case, were obtained by Holst, Nagy, and Tsogtgerel in 2008 for positive Yamabe metrics. However, their results are based on topological arguments, and as a result solution uniqueness is not known. Indeed, Maxwell gave evidence in 2011 that far-from-CMC solutions are not unique in certain cases. In this article, we provide further insight by establishing a type of alternative theorem for general far-from-CMC solutions. For a given manifold M that admits a metric of positive scalar curvature and scalar flat metric g(0) with no conformal Killing fields, we first prove existence of an analytic, one-parameter family of metrics g(z) through g(0) such that R(g(z)) = z. Using this family of metrics and given data (tau,sigma,rho,j), we form a one-parameter family of operators F((phi,w),z) whose zeros satisfy the conformal equations. Applying Liapnuov-Schmidt reduction, we determine an analytic solution curve for F((phi,w),z) = 0 through a critical point where the linearization of F((phi,w),z) vanishes. The regularity of this curve, the definition of F((phi,w),z), and the earlier far-from-CMC results of Holst et al. allow us to then prove the following alternative theorem for far-from-CMC solutions: either (1) there exists a z_1 >0 such that (positive Yamabe) solutions to the z_1-parameterized conformal equations are non-unique; or (2) there exists z_2 < 0 such that (negative Yamabe) solutions to the z_2-parameterized conformal equations exist.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convergence of Goal-Oriented Adaptive Finite Element Methods for Nonsymmetric Problems

In this article we develop convergence theory for a class of goal-oriented adaptive finite element algorithms for second order nonsymmetric linear elliptic equations. In particular, we establish contraction results for a method of this type for Dirichlet problems involving the elliptic operator L u = div (A grad u) - (b,grad u) - cu, with A Lipschitz, almost-everywhere symmetric positive definite, with b divergence-free, and with c >= 0. We first describe the problem class and review some standard facts concerning conforming finite element discretization and error-estimate-driven adaptive finite element methods (AFEM). We then describe a goal-oriented variation of standard AFEM (GOAFEM). Following the recent work of Mommer and Stevenson for symmetric problems, we establish contraction of GOAFEM and convergence in the sense of the goal function. Our analysis approach is signficantly different from that of Mommer and Stevenson, combining the recent contraction frameworks developed by Cascon, Kreuzer, Nochetto and Siebert; by Nochetto, Siebert and Veeser; and by Holst, Tsogtgerel and Zhu. We include numerical results demonstrating performance of our method with standard goal-oriented strategies on a convection problem.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalized Solutions to Semilinear Elliptic PDE with Applications to the Lichnerowicz Equation

In this article we investigate the existence of a solution to a semilinear, elliptic, partial differential equation with distributional coefficients and data. The problem we consider is a generalization of the Lichnerowicz equation that one encounters in studying the constraint equations in general relativity. Our method for solving this problem consists of solving a net of regularized, semilinear problems with data obtained by smoothing the original, distributional coefficients. In order to solve these regularized problems, we develop a priori pointwise bounds and sub- and super-solutions and then apply a fixed-point argument for order-preserving maps. We then show that the net of solutions obtained through this process satisfies certain decay estimates by determining estimates for the sub- and super-solutions and by utilizing classical, a priori elliptic estimates. The estimates for this net of solutions allow us to regard this collection of functions as a solution in a Colombeau-type algebra. We motivate this Colombeau algebra framework by first solving an ill-posed critical exponent problem. To solve this ill-posed problem, we use a collection of smooth, "approximating" problems and then use the resulting sequence of solutions and a compactness argument to obtain a solution to the original problem. This approach is modeled after the more general Colombeau framework that we develop, and it conveys the potential that solutions in these abstract spaces have for obtaining classical solutions to ill-posed nonlinear problems with irregular data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Lichnerowicz equation on compact manifolds with boundary

In this article we initiate a systematic study of the well-posedness theory of the Einstein constraint equations on compact manifolds with boundary. This is an important problem in general relativity, and it is particularly important in numerical relativity, as it arises in models of Cauchy surfaces containing asymptotically flat ends and/or trapped surfaces. Moreover, a number of technical obstacles that appear when developing the solution theory for open, asymptotically Euclidean manifolds have analogues on compact manifolds with boundary. As a first step, here we restrict ourselves to the Lichnerowicz equation, also called the Hamiltonian constraint equation, which is the main source of nonlinearity in the constraint system. The focus is on low regularity data and on the interaction between different types of boundary conditions, which has not been carefully analyzed before. In order to develop a well-posedness theory that mirrors the existing theory for the case of closed manifolds, we first generalize the Yamabe classification to nonsmooth metrics on compact manifolds with boundary. We then extend a result on conformal invariance to manifolds with boundary, and prove a uniqueness theorem. Finally, by using the method of sub- and super-solutions (order-preserving map iteration), we establish several existence results for a large class of problems covering a broad parameter regime, which includes most of the cases relevant in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local Multilevel Preconditioners for Elliptic Equations with Jump Coefficients on Bisection Grids

The goal of this paper is to design optimal multilevel solvers for the finite element approximation of second order linear elliptic problems with piecewise constant coefficients on bisection grids. Local multigrid and BPX preconditioners are constructed based on local smoothing only at the newest vertices and their immediate neighbors. The analysis of eigenvalue distributions for these local multilevel preconditioned systems shows that there are only a fixed number of eigenvalues which are deteriorated by the large jump. The remaining eigenvalues are bounded uniformly with respect to the coefficients and the meshsize. Therefore, the resulting preconditioned conjugate gradient algorithm will converge with an asymptotic rate independent of the coefficients and logarithmically with respect to the meshsize. As a result, the overall computational complexity is nearly optimal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multilevel Preconditioners for Discontinuous Galerkin Approximations of Elliptic Problems with Jump Coefficients

We introduce and analyze two-level and multi-level preconditioners for a family of Interior Penalty (IP) discontinuous Galerkin (DG) discretizations of second order elliptic problems with large jumps in the diffusion coefficient. Our approach to IPDG-type methods is based on a splitting of the DG space into two components that are orthogonal in the energy inner product naturally induced by the methods. As a result, the methods and their analysis depend in a crucial way on the diffusion coefficient of the problem. The analysis of the proposed preconditioners is presented for both symmetric and non-symmetric IP schemes; dealing simultaneously with the jump in the diffusion coefficient and the non-nested character of the relevant discrete spaces presents extra difficulties in the analysis which precludes a simple extension of existing results. However, we are able to establish robustness (with respect to the diffusion coefficient) and nearly-optimality (up to a logarithmic term depending on the mesh size) for both two-level and BPX-type preconditioners. Following the analysis, we present a sequence of detailed numerical results which verify the theory and illustrate the performance of the methods. The paper includes an Appendix with a collection of proofs of several technical results required for the analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-uniqueness of Solutions to the Conformal Formulation

It is well-known that solutions to the conformal formulation of the Einstein constraint equations are unique in the cases of constant mean curvature (CMC) and near constant mean curvature (near-CMC). However, the new far-from-constant mean curvature (far-from-CMC) existence results due to Holst, Nagy, and Tsogtgerel in 2008, to Maxwell in 2009, and to Dahl, Gicquaud and Humbert in 2010, are based on degree theory rather than on the (uniqueness-providing) contraction arguments that had been used for all non-CMC existence results prior to 2008. In fact, Maxwell demonstrated in 2011 that solutions are non-unique in the far-from-CMC case for certain types of low-regularity mean curvature. In this article, we investigate uniqueness properties of solutions to the Einstein constraint equations on closed manifolds using tools from bifurcation theory. For positive, constant scalar curvature and constant mean curvature, we first demonstrate existence of a critical energy density for the Hamiltonian constraint with unscaled matter sources. We then show that for this choice of energy density, the linearization of the elliptic system develops a one-dimensional kernel in both the CMC and non-CMC (near and far) cases. Using Liapunov-Schmidt reduction and standard tools from nonlinear analysis, we demonstrate that solutions to the conformal formulation with unscaled data are non-unique by determining an explicit solution curve, and by analyzing its behavior in the neighborhood of a particular solution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Two-Grid Methods for Semilinear Interface Problems

In this article we consider two-grid finite element methods for solving semilinear interface problems in d space dimensions, for d=2 or d=3. We first describe in some detail the target problem class with discontinuous diffusion coefficients, which includes problems containing sub-critical, critical, and supercritical nonlinearities. We then establish basic quasi-optimal a priori error estimate for Galerkin approximations. In the critical and subcritical cases, we follow our recent approach to controling the nonlinearity using only pointwise control of the continuous solution and a local Lipschitz property, rather than through pointwise control of the discrete solution; this eliminates the requirement that the discrete solution satisfy a discrete form of the maximum principle, hence eliminating the need for restrictive angle conditions in the underlying mesh. The supercritical case continues to require such mesh conditions in order to control the nonlinearity. We then design a two-grid algorithm consisting of a coarse grid solver for the original nonlinear problem, and a fine grid solver for a linearized problem. We analyze the quality of approximations generated by the algorithm, and show that the coarse grid may be taken to have much larger elements than the fine grid, and yet one can still obtain approximation quality that is asymptotically as good as solving the original nonlinear problem on the fine mesh. The algorithm we describe, and its analysis in this article, combines four sets of tools: the work of Xu and Zhou on two-grid algorithms for semilinear problems; the recent results for linear interface problems due to Li, Melenk, Wohlmuth, and Zou; recent work on the Poisson-Boltzmann equation; and recent work on a priori estimates for semilinear problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adaptive Finite Element Methods with Inexact Solvers for the Nonlinear Poisson-Boltzmann Equation

In this article we study adaptive finite element methods (AFEM) with inexact solvers for a class of semilinear elliptic interface problems. We are particularly interested in nonlinear problems with discontinuous diffusion coefficients, such as the nonlinear Poisson-Boltzmann equation and its regularizations. The algorithm we study consists of the standard SOLVE-ESTIMATE-MARK-REFINE procedure common to many adaptive finite element algorithms, but where the SOLVE step involves only a full solve on the coarsest level, and the remaining levels involve only single Newton updates to the previous approximate solution. We summarize a recently developed AFEM convergence theory for inexact solvers, and present a sequence of numerical experiments that give evidence that the theory does in fact predict the contraction properties of AFEM with inexact solvers. The various routines used are all designed to maintain a linear-time computational complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Barrier methods for critical exponent problems in geometric analysis and mathematical physics

We consider the design and analysis of numerical methods for approximating positive solutions to nonlinear geometric elliptic partial differential equations containing critical exponents. This class of problems includes the Yamabe problem and the Einstein constraint equations, which simultaneously contain several challenging features: high spatial dimension n >= 3, varying (potentially non-smooth) coefficients, critical (even super-critical) nonlinearity, non-monotone nonlinearity (arising from a non-convex energy), and spatial domains that are typically Riemannian manifolds rather than simply open sets in Rn. These problems may exhibit multiple solutions, although only positive solutions typically have meaning. This creates additional complexities in both the theory and numerical treatment of such problems, as this feature introduces both non-uniqueness as well as the need to incorporate an inequality constraint into the formulation. In this work, we consider numerical methods based on Galerkin-type discretization, covering any standard bases construction (finite element, spectral, or wavelet), and the combination of a barrier method for nonconvex optimization and global inexact Newton-type methods for dealing with nonconvexity and the presence of inequality constraints. We first give an overview of barrier methods in non-convex optimization, and then develop and analyze both a primal barrier energy method for this class of problems. We then consider a sequence of numerical experiments using this type of barrier method, based on a particular Galerkin method, namely the piecewise linear finite element method, leverage the FETK modeling package. We illustrate the behavior of the primal barrier energy method for several examples, including the Yamabe problem and the Hamiltonian constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite Element Error Estimates for Critical Growth Semilinear Problems without Angle Conditions

In this article we consider a priori error and pointwise estimates for finite element approximations of solutions to semilinear elliptic boundary value problems in d>=2 space dimensions, with nonlinearities satisfying critical growth conditions. It is well-understood how mesh geometry impacts finite element interpolant quality, and leads to the reasonable notion of shape regular simplex meshes. It is also well-known how to perform both mesh generation and simplex subdivision, in arbitrary space dimension, so as to guarantee the entire hierarchy of nested simplex meshes produced through subdivision continue to satisfy shape regularity. However, much more restrictive angle conditions are needed for basic a priori quasi-optimal error estimates, as well as for a priori pointwise estimates. These angle conditions, which are particularly difficult to satisfy in three dimensions in any type of unstructured or adaptive setting, are needed to gain pointwise control of the nonlinearity through discrete maximum principles. This represents a major gap in finite element approximation theory for nonlinear problems on unstructured meshes, and in particular for adaptive methods. In this article, we close this gap in the case of semilinear problems with critical or sub-critical nonlinear growth, by deriving a priori estimates directly, without requiring the discrete maximum principle, and hence eliminating the need for restrictive angle conditions. Our main result is a type of local Lipschitz property that relies only on the continuous maximum principle, together with the growth condition. We also show that under some additional smoothness assumptions, the a priori error estimate itself is enough to give pointwise control the discrete solution, without the need for restrictive angle conditions. Numerical experiments confirm our theoretical conclusions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Geometric variational crimes: Hilbert complexes, finite element exterior calculus, and problems on hypersurfaces

A recent paper of Arnold, Falk, and Winther [Bull AMS, 47 (2010)] showed that a large class of mixed finite element methods can be formulated naturally on Hilbert complexes, where using a Galerkin-like approach, one solves a variational problem on a finite-dimensional subcomplex. In a seemingly unrelated research direction, Dziuk [Lect Notes in Math, vol 1357 (1988)] analyzed a class of nodal finite elements for the Laplace-Beltrami equation on smooth 2-surfaces approximated by a piecewise-linear triangulation; Demlow later extended this analysis [SIAM J Numer Anal, 47 (2009)] to 3-surfaces, as well as to higher-order surface approximation. In this article, we bring these lines of research together, first developing a framework for the analysis of variational crimes in abstract Hilbert complexes, and then applying this abstract framework to the setting of finite element exterior calculus on hypersurfaces. Our framework extends the work of Arnold, Falk, and Winther to problems that violate their subcomplex assumption, allowing for the extension of finite element exterior calculus to approximate domains, most notably the Hodge-de Rham complex on approximate manifolds. As an application of the latter, we recover Dziuk's and Demlow's a priori estimates for 2- and 3-surfaces, demonstrating that surface finite element methods can be analyzed completely within this abstract framework. Moreover, our results generalize these earlier estimates dramatically, extending them from nodal finite elements for Laplace-Beltrami to mixed finite elements for the Hodge Laplacian, and from 2- and 3-dimensional hypersurfaces to those of arbitrary dimension. By developing this analytical framework using a combination of general tools from differential geometry and functional analysis, we are led to a more geometric analysis of surface finite element methods, whereby the main results become more transparent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Goal-Oriented Adaptivity and Multilevel Preconditioning for the Poisson-Boltzmann Equation

In this article, we develop goal-oriented error indicators to drive adaptive refinement algorithms for the Poisson-Boltzmann equation. Empirical results for the solvation free energy linear functional demonstrate that goal-oriented indicators are not sufficient on their own to lead to a superior refinement algorithm. To remedy this, we propose a problem-specific marking strategy using the solvation free energy computed from the solution of the linear regularized Poisson-Boltzmann equation. The convergence of the solvation free energy using this marking strategy, combined with goal-oriented refinement, compares favorably to adaptive methods using an energy-based error indicator. Due to the use of adaptive mesh refinement, it is critical to use multilevel preconditioning in order to maintain optimal computational complexity. We use variants of the classical multigrid method, which can be viewed as generalizations of the hierarchical basis multigrid and Bramble-Pasciak-Xu (BPX) preconditioners.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients

In this paper, we present a multigrid preconditioner for solving the linear system arising from the piecewise linear nonconforming Crouzeix-Raviart discretization of second order elliptic problems with jump coefficients. The preconditioner uses the standard conforming subspaces as coarse spaces. Numerical tests show both robustness with respect to the jump in the coefficient and near-optimality with respect to the number of degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Semilinear mixed problems on Hilbert complexes and their numerical approximation

Arnold, Falk, and Winther recently showed [Bull. Amer. Math. Soc. 47 (2010), 281-354] that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article [arXiv:1005.4455], we extended the Arnold-Falk-Winther framework by analyzing variational crimes (a la Strang) on Hilbert complexes. In particular, this gave a treatment of finite element exterior calculus on manifolds, generalizing techniques from surface finite element methods and recovering earlier a priori estimates for the Laplace-Beltrami operator on 2- and 3-surfaces, due to Dziuk [Lecture Notes in Math., vol. 1357 (1988), 142-155] and later Demlow [SIAM J. Numer. Anal., 47 (2009), 805-827], as special cases. In the present article, we extend the Hilbert complex framework in a second distinct direction: to the study of semilinear mixed problems. We do this, first, by introducing an operator-theoretic reformulation of the linear mixed problem, so that the semilinear problem can be expressed as an abstract Hammerstein equation. This allows us to obtain, for semilinear problems, a priori solution estimates and error estimates that reduce to the Arnold-Falk-Winther results in the linear case. We also consider the impact of variational crimes, extending the results of our previous article to these semilinear problems. As an immediate application, this new framework allows for mixed finite element methods to be applied to semilinear problems on surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Nonlinear Elasticity Model of Macromolecular Conformational Change Induced by Electrostatic Forces

In this paper we propose a nonlinear elasticity model of macromolecular conformational change (deformation) induced by electrostatic forces generated by an implicit solvation model. The Poisson-Boltzmann equation for the electrostatic potential is analyzed in a domain varying with the elastic deformation of molecules, and a new continuous model of the electrostatic forces is developed to ensure solvability of the nonlinear elasticity equations. We derive the estimates of electrostatic forces corresponding to four types of perturbations to an electrostatic potential field, and establish the existance of an equilibrium configuration using a fixed-point argument, under the assumption that the change in the ionic strength and charges due to the additional molecules causing the deformation are sufficiently small. The results are valid for elastic models with arbitrarily complex dielectric interfaces and cavities, and can be generalized to large elastic deformation caused by high ionic strength, large charges, and strong external fields by using continuation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Adaptive Finite Element Modeling Techniques for the Poisson-Boltzmann Equation

We develop an efficient and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the regularization technique of Chen, Holst, and Xu; this technique made possible the first a priori pointwise estimates and the first complete solution and approximation theory for the Poisson-Boltzmann equation. It also made possible the first provably convergent discretization of the PBE, and allowed for the development of a provably convergent AFEM for the PBE. However, in practice the regularization turns out to be numerically ill-conditioned. In this article, we examine a second regularization, and establish a number of basic results to ensure that the new approach produces the same mathematical advantages of the original regularization, without the ill-conditioning property. We then design an AFEM scheme based on the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori pointwise estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. The stability advantages of the new regularization are demonstrated using an FETK-based implementation, through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Adaptive Numerical Treatment of Elliptic Systems on Manifolds

Adaptive multilevel finite element methods are developed and analyzed for certain elliptic systems arising in geometric analysis and general relativity. This class of nonlinear elliptic systems of tensor equations on manifolds is first reviewed, and then adaptive multilevel finite element methods for approximating solutions to this class of problems are considered in some detail. Two a posteriori error indicators are derived, based on local residuals and on global linearized adjoint or dual problems. The design of Manifold Code (MC) is then discussed; MC is an adaptive multilevel finite element software package for 2- and 3-manifolds. It employs a posteriori error estimation, adaptive simplex subdivision, unstructured algebraic multilevel methods, global inexact Newton methods, and numerical continuation methods for the numerical solution of nonlinear covariant elliptic systems on 2- and 3-manifolds. Some of the more interesting features of MC are described in detail, including some new ideas for topology and geometry representation in simplex meshes, and an unusual partition of unity-based method for exploiting parallel computers. A short example is then given which involves the Hamiltonian and momentum constraints in the Einstein equations, a representative nonlinear 4-component covariant elliptic system on a Riemannian 3-manifold which arises in general relativity. A number of operator properties and solvability results recently established are first summarized, making possible two quasi-optimal a priori error estimates for Galerkin approximations which are then derived. These two results complete the theoretical framework for effective use of adaptive multilevel finite element methods. A sample calculation using the MC software is then presented.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Applications of Domain Decomposition and Partition of Unity Methods in Physics and Geometry

We consider a class of adaptive multilevel domain decomposition-like algorithms, built from a combination of adaptive multilevel finite element, domain decomposition, and partition of unity methods. These algorithms have several interesting features such as very low communication requirements, and they inherit a simple and elegant approximation theory framework from partition of unity methods. They are also very easy to use with highly complex sequential adaptive finite element packages, requiring little or no modification of the underlying sequential finite element software. The parallel algorithm can be implemented as a simple loop which starts off a sequential local adaptive solve on a collection of processors simultaneously. We first review the Partition of Unity Method (PUM) of Babuvska and Melenk, and outline the PUM approximation theory framework. We then describe a variant we refer to here as the Parallel Partition of Unity Method (PPUM), which is a combination of the Partition of Unity Method with the parallel adaptive algorithm of Bank and Holst. We then derive two global error estimates for PPUM, by exploiting the PUM analysis framework it inherits, and by employing some recent local estimates of Xu and Zhou. We then discuss a duality-based variant of PPUM which is more appropriate for certain applications, and we derive a suitable variant of the PPUM approximation theory framework. Our implementation of PPUM-type algorithms using the FETK and MC software packages is described. We then present a short numerical example involving the Einstein constraints arising in gravitational wave models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence and Optimality of Adaptive Mixed Finite Element Methods

The convergence and optimality of adaptive mixed finite element methods for the Poisson equation are established in this paper. The main difficulty for mixed finite element methods is the lack of minimization principle and thus the failure of orthogonality. A quasi-orthogonality property is proved using the fact that the error is orthogonal to the divergence free subspace, while the part of the error that is not divergence free can be bounded by the data oscillation using a discrete stability result. This discrete stability result is also used to get a localized discrete upper bound which is crucial for the proof of the optimality of the adaptive approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Determining Projections and Functionals for Weak Solutions of the Navier-Stokes Equations

In this paper we prove that an operator which projects weak solutions of the two- or three-dimensional Navier-Stokes equations onto a finite-dimensional space is determining if it annihilates the difference of two "nearby" weak solutions asymptotically, and if it satisfies a single appoximation inequality. We then apply this result to show that the long-time behavior of weak solutions to the Navier-Stokes equations, in both two- and three-dimensions, is determined by the long-time behavior of a finite set of bounded linear functionals. These functionals are constructed by local surface averages of solutions over certain simplex volume elements, and are therefore well-defined for weak solutions. Moreover, these functionals define a projection operator which satisfies the necessary approximation inequality for our theory. We use the general theory to establish lower bounds on the simplex diameters in both two- and three-dimensions. Furthermore, in the three dimensional case we make a connection between their diameters and the Kolmogoroff dissipation small scale in turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local Convergence of Adaptive Methods for Nonlinear Partial Differential Equations

In this article we develop convergence theory for a general class of adaptive approximation algorithms for abstract nonlinear operator equations on Banach spaces, and use the theory to obtain convergence results for practical adaptive finite element methods (AFEM) applied to several classes of nonlinear elliptic equations. In the first part of the paper, we develop a weak-* convergence framework for nonlinear operators, whose Gateaux derivatives are locally Lipschitz and satisfy a local inf-sup condition. The framework can be viewed as extending the recent convergence results for linear problems of Morin, Siebert and Veeser to a general nonlinear setting. We formulate an abstract adaptive approximation algorithm for nonlinear operator equations in Banach spaces with local structure. The weak-* convergence framework is then applied to this class of abstract locally adaptive algorithms, giving a general convergence result. The convergence result is then applied to a standard AFEM algorithm in the case of several semilinear and quasi-linear scalar elliptic equations and elliptic systems, including: a semilinear problem with subcritical nonlinearity, the steady Navier-Stokes equations, and a quasilinear problem with nonlinear diffusion. This yields several new AFEM convergence results for these nonlinear problems. In the second part of the paper we develop a second abstract convergence framework based on strong contraction, extending the recent contraction results for linear problems of Cascon, Kreuzer, Nochetto, and Siebert and of Mekchay and Nochetto to abstract nonlinear problems. The contraction result is then applied to a standard AFEM algorithm for semilinear problems with sub- and super-critical nonlinearities and for the Hamiltonian constraint in general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local Refinement and Multilevel Preconditioning: Implementation and Numerical Experiments

In this paper, we examine a number of additive and multiplicative multilevel iterative methods and preconditioners in the setting of two-dimensional local mesh refinement. While standard multilevel methods are effective for uniform refinement-based discretizations of elliptic equations, they tend to be less effective for algebraic systems which arise from discretizations on locally refined meshes, losing their optimal behavior in both storage and computational complexity. Our primary focus here is on BPX-style additive and multiplicative multilevel preconditioners, and on various stabilizations of the additive and multiplicative hierarchical basis method (HB), and their use in the local mesh refinement setting. In this article, we describe in detail the implementation of these types of algorithms, including detailed discussions of the datastructures and traversal algorithms we employ for obtaining optimal storage and computational complexity in our implementations. We show how each of the algorithms can be implemented using standard datatypes available in languages such as C and FORTRAN, so that the resulting algorithms have optimal (linear) storage requirements, and so that the resulting multilevel method or preconditioner can be applied with optimal (linear) computational costs. Our implementations are performed in both C and MATLAB using the Finite Element ToolKit (FETK), an open source finite element software package. We finish the paper with a sequence of numerical experiments illustrating the effectiveness of a number of BPX and stabilized HB variants for several examples requiring local refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimality of multilevel preconditioners for local mesh refinement in three dimensions

In this article, we establish optimality of the Bramble-Pasciak-Xu (BPX) norm equivalence and optimality of the wavelet modified (or stabilized) hierarchical basis (WHB) preconditioner in the setting of local 3D mesh refinement. In the analysis of WHB methods, a critical first step is to establish the optimality of BPX norm equivalence for the refinement procedures under consideration. While the available optimality results for the BPX norm have been constructed primarily in the setting of uniformly refined meshes, a notable exception is the local 2D red-green result due to Dahmen and Kunoth. The purpose of this article is to extend this original 2D optimality result to the local 3D red-green refinement procedure introduced by Bornemann-Erdmann-Kornhuber (BEK), and then to use this result to extend the WHB optimality results from the quasiuniform setting to local 2D and 3D red-green refinement scenarios. The BPX extension is reduced to establishing that locally enriched finite element subspaces allow for the construction of a scaled basis which is formally Riesz stable. It is possible to show that the number of degrees of freedom used for smoothing is bounded by a constant times the number of degrees of freedom introduced at that level of refinement, indicating that a practical implementable version of the resulting BPX preconditioner for the BEK refinement setting has provably optimal (linear) computational complexity per iteration. An interesting implication of the optimality of the WHB preconditioner is the a priori H1-stability of the L2-projection. The theoretical framework employed supports arbitrary spatial dimension d >= 1 and requires no coefficient smoothness assumptions beyond those required for well-posedness in H1.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Schwarz Methods: To Symmetrize or Not to Symmetrize

A preconditioning theory is presented which establishes sufficient conditions for multiplicative and additive Schwarz algorithms to yield self-adjoint positive definite preconditioners. It allows for the analysis and use of non-variational and non-convergent linear methods as preconditioners for conjugate gradient methods, and it is applied to domain decomposition and multigrid. It is illustrated why symmetrizing may be a bad idea for linear methods. It is conjectured that enforcing minimal symmetry achieves the best results when combined with conjugate gradient acceleration. Also, it is shown that absence of symmetry in the linear preconditioner is advantageous when the linear method is accelerated by using the Bi-CGstab method. Numerical examples are presented for two test problems which illustrate the theory and conjectures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Finite Element Approximation of the Nonlinear Poisson-Boltzmann Equation

A widely used electrostatics model in the biomolecular modeling community, the nonlinear Poisson-Boltzmann equation, along with its finite element approximation, are analyzed in this paper. A regularized Poisson-Boltzmann equation is introduced as an auxiliary problem, making it possible to study the original nonlinear equation with delta distribution sources. A priori error estimates for the finite element approximation are obtained for the regularized Poisson-Boltzmann equation based on certain quasi-uniform grids in two and three dimensions. Adaptive finite element approximation through local refinement driven by an a posteriori error estimate is shown to converge. The Poisson-Boltzmann equation does not appear to have been previously studied in detail theoretically, and it is hoped that this paper will help provide molecular modelers with a better foundation for their analytical and computational work with the Poisson-Boltzmann equation. Note that this article apparently gives the first rigorous convergence result for a numerical discretization technique for the nonlinear Poisson-Boltzmann equation with delta distribution sources, and it also introduces the first provably convergent adaptive method for the equation. This last result is currently one of only a handful of existing convergence results of this type for nonlinear problems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Analysis of a General Family of Regularized Navier-Stokes and MHD Models

We consider a general family of regularized Navier-Stokes and Magnetohydrodynamics (MHD) models on n-dimensional smooth compact Riemannian manifolds with or without boundary, with n greater than or equal to 2. This family captures most of the specific regularized models that have been proposed and analyzed in the literature, including the Navier-Stokes equations, the Navier-Stokes-alpha model, the Leray-alpha model, the Modified Leray-alpha model, the Simplified Bardina model, the Navier-Stokes-Voight model, the Navier-Stokes-alpha-like models, and certain MHD models, in addition to representing a larger 3-parameter family of models not previously analyzed. We give a unified analysis of the entire three-parameter family using only abstract mapping properties of the principle dissipation and smoothing operators, and then use specific parameterizations to obtain the sharpest results. We first establish existence and regularity results, and under appropriate assumptions show uniqueness and stability. We then establish results for singular perturbations, including the inviscid and alpha limits. Next we show existence of a global attractor for the general model, and give estimates for its dimension. We finish by establishing some results on determining operators for subfamilies of dissipative and non-dissipative models. In addition to establishing a number of results for all models in this general family, the framework recovers most of the previous results on existence, regularity, uniqueness, stability, attractor existence and dimension, and determining operators for well-known members of this family.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Solving the Einstein constraint equations on multi-block triangulations using finite element methods

In order to generate initial data for nonlinear relativistic simulations, one needs to solve the Einstein constraints, which can be cast into a coupled set of nonlinear elliptic equations. Here we present an approach for solving these equations on three-dimensional multi-block domains using finite element methods. We illustrate our approach on a simple example of Brill wave initial data, with the constraints reducing to a single linear elliptic equation for the conformal factor $ψ$. We use quadratic Lagrange elements on semi-structured simplicial meshes, obtained by triangulation of multi-block grids. In the case of uniform refinement the scheme is superconvergent at most mesh vertices, due to local symmetry of the finite element basis with respect to local spatial inversions. We show that in the superconvergent case subsequent unstructured mesh refinements do not improve the quality of our initial data. As proof of concept that this approach is feasible for generating multi-block initial data in three dimensions, after constructing the initial data we evolve them in time using a high order finite-differencing multi-block approach and extract the gravitational waves from the numerical solution.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Rough solutions of the Einstein constraints on closed manifolds without near-CMC conditions

We consider the conformal decomposition of Einstein's constraint equations introduced by Lichnerowicz and York, on a closed manifold. We establish existence of non-CMC weak solutions using a combination of a priori estimates for the individual Hamiltonian and momentum constraints, barrier constructions for the Hamiltonian constraint, and topological fixed-point arguments. An important new feature of these results is the absense of the near-CMC assumption when the rescaled background metric is in the positive Yamabe class, if the freely specifiable part of the data given by the matter fields (if present) and the traceless-transverse part of the rescaled extrinsic curvature are taken to be sufficiently small. In this case, the mean extrinsic curvature can be taken to be an arbitrary smooth function without restrictions on the size of its spatial derivatives, giving what are apparently the first non-CMC existence results without the near-CMC assumption. Standard bootstrapping arguments to increase the regularity of the conformal factor are blocked by the use of a weak background metric. In the CMC case, we recover Maxwell's rough solution results as a special case. Our results extend the 1996 non-CMC result of Isenberg and Moncrief in three ways: (1) the near-CMC assumption is removed in the case of the positive Yamabe class; (2) regularity is extended down to the maximum allowed by the background metric and the matter; and (3) the result holds for all three Yamabe classes. This last extension was also accomplished recently by Allen, Clausen and Isenberg, although their result is restricted to the near-CMC case and to smoother background metrics and data.

0 citations0 reviewsNo code linkedOpen access