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

Evan S. Gawlik

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

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.comp-phmath-phmath.DGmath.MP

Catalog footprint

What is connected

8works
6topics
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

8 published item(s)

preprint2023arXiv

Finite element approximation of scalar curvature in arbitrary dimension

We analyze finite element discretizations of scalar curvature in dimension $N \ge 2$. Our analysis focuses on piecewise polynomial interpolants of a smooth Riemannian metric $g$ on a simplicial triangulation of a polyhedral domain $Ω\subset \mathbb{R}^N$ having maximum element diameter $h$. We show that if such an interpolant $g_h$ has polynomial degree $r \ge 0$ and possesses single-valued tangential-tangential components on codimension-1 simplices, then it admits a natural notion of (densitized) scalar curvature that converges in the $H^{-2}(Ω)$-norm to the (densitized) scalar curvature of $g$ at a rate of $O(h^{r+1})$ as $h \to 0$, provided that either $N = 2$ or $r \ge 1$. As a special case, our result implies the convergence in $H^{-2}(Ω)$ of the widely used "angle defect" approximation of Gaussian curvature on two-dimensional triangulations, without stringent assumptions on the interpolated metric $g_h$. We present numerical experiments that indicate that our analytical estimates are sharp.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Conservative Finite Element Method for the Incompressible Euler Equations with Variable Density

We construct a finite element discretization and time-stepping scheme for the incompressible Euler equations with variable density that exactly preserves total mass, total squared density, total energy, and pointwise incompressibility. The method uses Raviart-Thomas or Brezzi-Douglas-Marini finite elements to approximate the velocity and discontinuous polynomials to approximate the density and pressure. To achieve exact preservation of the aforementioned conserved quantities, we exploit a seldom-used weak formulation of the momentum equation and a second-order time-stepping scheme that is similar, but not identical, to the midpoint rule. We also describe and prove stability of an upwinded version of the method. We present numerical examples that demonstrate the order of convergence of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing the Fréchet Derivative of the Polar Decomposition

We derive iterative methods for computing the Fréchet derivative of the map which sends a full-rank matrix $A$ to the factor $U$ in its polar decomposition $A=UH$, where $U$ has orthonormal columns and $H$ is Hermitian positive definite. The methods apply to square matrices as well as rectangular matrices having more rows than columns. Our derivation relies on a novel identity that relates the Fréchet derivative of the polar decomposition to the matrix sign function $\mathrm{sign}(X) = X (X^2)^{-1/2}$ applied to a certain block matrix $X$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Embedding-Based Interpolation on the Special Orthogonal Group

We study schemes for interpolating functions that take values in the special orthogonal group $SO(n)$. Our focus is on interpolation schemes obtained by embedding $SO(n)$ in a linear space, interpolating in the linear space, and mapping the result onto $SO(n)$ via the closest point projection. The resulting interpolants inherit both the order of accuracy and the regularity of the underlying interpolants on the linear space. The values and derivatives of the interpolants admit efficient evaluation via either explicit formulas or iterative algorithms, which we detail for two choices of embeddings: the embedding of $SO(n)$ in the space of $n \times n$ matrices and, when $n=3$, the identification of $SO(3)$ with the set of unit quaternions. Along the way, we point out a connection between these interpolation schemes and geodesic finite elements. We illustrate the utility of these interpolation schemes by numerically computing minimum acceleration curves on $SO(n)$, a task which is handled naturally with $SO(n)$-valued finite elements having $C^1$-continuity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Universal Meshes for the Simulation of Brittle Fracture and Moving Boundary Problems

Universal meshes have recently appeared in the literature as a compu- tationally efficient and robust paradigm for the generation of conforming simpli- cial meshes for domains with evolving boundaries. The main idea behind a univer- sal mesh is to immerse the moving boundary in a background mesh (the universal mesh), and to produce a mesh that conforms to the moving boundary at any given time by adjusting a few of elements of the background mesh. In this manuscript we present the application of universal meshes to the simulation of brittle fracturing. To this extent, we provide a high level description of a crack propagation algorithm and showcase its capabilities. Alongside universal meshes for the simulation of brit- tle fracture, we provide other examples for which universal meshes prove to be a powerful tool, namely fluid flow past moving obstacles. Lastly, we conclude the manuscript with some remarks on the current state of universal meshes and future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Order Finite Element Methods for Moving Boundary Problems with Prescribed Boundary Evolution

We introduce a framework for the design of finite element methods for two-dimensional moving boundary problems with prescribed boundary evolution that have arbitrarily high order of accuracy, both in space and in time. At the core of our approach is the use of a universal mesh: a stationary background mesh containing the domain of interest for all times that adapts to the geometry of the immersed domain by adjusting a small number of mesh elements in the neighborhood of the moving boundary. The resulting method maintains an exact representation of the (prescribed) moving boundary at the discrete level, or an approximation of the appropriate order, yet is immune to large distortions of the mesh under large deformations of the domain. The framework is general, making it possible to achieve any desired order of accuracy in space and time by selecting a preferred and suitable finite-element space on the universal mesh for the problem at hand, and a preferred and suitable time integrator for ordinary differential equations. We illustrate our approach by constructing a particular class of methods, and apply them to a prescribed-boundary variant of the Stefan problem. We present numerical evidence for the order of accuracy of our schemes in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Supercloseness of Orthogonal Projections onto Nearby Finite Element Spaces

We derive upper bounds on the difference between the orthogonal projections of a smooth function $u$ onto two finite element spaces that are nearby, in the sense that the support of every shape function belonging to one but not both of the spaces is contained in a common region whose measure tends to zero under mesh refinement. The bounds apply, in particular, to the setting in which the two finite element spaces consist of continuous functions that are elementwise polynomials over shape-regular, quasi-uniform meshes that coincide except on a region of measure $O(h^γ)$, where $γ$ is a nonnegative scalar and $h$ is the mesh spacing. The projector may be, for example, the orthogonal projector with respect to the $L^2$- or $H^1$-inner product. In these and other circumstances, the bounds are superconvergent under a few mild regularity assumptions. That is, under mesh refinement, the two projections differ in norm by an amount that decays to zero at a faster rate than the amounts by which each projection differs from $u$. We present numerical examples to illustrate these superconvergent estimates and verify the necessity of the regularity assumptions on $u$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometric, Variational Discretization of Continuum Theories

This study derives geometric, variational discretizations of continuum theories arising in fluid dynamics, magnetohydrodynamics (MHD), and the dynamics of complex fluids. A central role in these discretizations is played by the geometric formulation of fluid dynamics, which views solutions to the governing equations for perfect fluid flow as geodesics on the group of volume-preserving diffeomorphisms of the fluid domain. Inspired by this framework, we construct a finite-dimensional approximation to the diffeomorphism group and its Lie algebra, thereby permitting a variational temporal discretization of geodesics on the spatially discretized diffeomorphism group. The extension to MHD and complex fluid flow is then made through an appeal to the theory of Euler-Poincaré systems with advection, which provides a generalization of the variational formulation of ideal fluid flow to fluids with one or more advected parameters. Upon deriving a family of structured integrators for these systems, we test their performance via a numerical implementation of the update schemes on a cartesian grid. Among the hallmarks of these new numerical methods are exact preservation of momenta arising from symmetries, automatic satisfaction of solenoidal constraints on vector fields, good long-term energy behavior, robustness with respect to the spatial and temporal resolution of the discretization, and applicability to irregular meshes.

0 citations0 reviewsNo code linkedOpen access