BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Michael Holst

Michael Holst contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAmath.APmath-phmath.MPgr-qcmath.OCNumerical Analysisphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 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
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

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
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
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