Researcher profile

Ari Stern

Ari Stern contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
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

6 published item(s)

preprint2025arXiv

Finite element exterior calculus for time-dependent Hamiltonian partial differential equations

The success of symplectic integrators for Hamiltonian ODEs has led to a decades-long program of research seeking analogously structure-preserving numerical methods for Hamiltonian PDEs. In this paper, we construct a large class of such methods by combining finite element exterior calculus (FEEC) for spatial semidiscretization with symplectic integrators for time discretization. The resulting methods satisfy a local multisymplectic conservation law in space and time, which generalizes the symplectic conservation law of Hamiltonian ODEs, and which carries finer information about Hamiltonian structure than other approaches based on global function spaces. We give particular attention to conforming FEEC methods and hybridizable discontinuous Galerkin (HDG) methods. The theory and methods are illustrated by application to the semilinear Hodge wave equation.

preprint2020arXiv

Structure-preserving numerical integrators for Hodgkin-Huxley-type systems

Motivated by the Hodgkin-Huxley model of neuronal dynamics, we study explicit numerical integrators for "conditionally linear" systems of ordinary differential equations. We show that splitting and composition methods, when applied to the Van der Pol oscillator and to the Hodgkin-Huxley model, do a better job of preserving limit cycles of these systems for large time steps, compared with the "Euler-type" methods (including Euler's method, exponential Euler, and semi-implicit Euler) commonly used in computational neuroscience, with no increase in computational cost. These limit cycles are important to preserve, due to their role in neuronal spiking. Splitting methods even compare favorably to the explicit exponential midpoint method, which is twice as expensive per step. The second-order Strang splitting method is seen to perform especially well across a range of non-stiff and stiff dynamics.

preprint2019arXiv

Invariant connections, Lie algebra actions, and foundations of numerical integration on manifolds

Motivated by numerical integration on manifolds, we relate the algebraic properties of invariant connections to their geometric properties. Using this perspective, we generalize some classical results of Cartan and Nomizu to invariant connections on algebroids. This has fundamental consequences for the theory of numerical integrators, giving a characterization of the spaces on which Butcher and Lie-Butcher series methods, which generalize Runge-Kutta methods, may be applied.

preprint2019arXiv

Multisymplecticity of hybridizable discontinuous Galerkin methods

In this paper, we prove necessary and sufficient conditions for a hybridizable discontinuous Galerkin (HDG) method to satisfy a multisymplectic conservation law, when applied to a canonical Hamiltonian system of partial differential equations. We show that these conditions are satisfied by the "hybridized" versions of several of the most commonly-used finite element methods, including mixed, nonconforming, and discontinuous Galerkin methods. (Interestingly, for the continuous Galerkin method in dimension greater than one, we show that multisymplecticity only holds in a weaker sense.) Consequently, these general-purpose finite element methods may be used for structure-preserving discretization (or semidiscretization) of canonical Hamiltonian systems of ODEs or PDEs. This establishes multisymplecticity for a large class of arbitrarily-high-order methods on unstructured meshes.

preprint2015arXiv

Explicit formulas for reaction probability in reaction-diffusion experiments

A computational procedure is developed for determining the conversion probability for reaction-diffusion systems in which a first-order catalytic reaction is performed over active particles. We apply this general method to systems on metric graphs, which may be viewed as 1-dimensional approximations of 3-dimensional systems, and obtain explicit formulas for conversion. We then study numerically a class of 3-dimensional systems and test how accurately they are described by model formulas obtained for metric graphs. The optimal arrangement of active particles in a 1-dimensional multiparticle system is found, which is shown to depend on the level of catalytic activity: conversion is maximized for low catalytic activity when all particles are bunched together close to the point of gas injection, and for high catalytic activity when the particles are evenly spaced.

preprint2010arXiv

Discrete Hamilton-Pontryagin mechanics and generating functions on Lie groupoids

We present a discrete analog of the recently introduced Hamilton-Pontryagin variational principle in Lagrangian mechanics. This unifies two, previously disparate approaches to discrete Lagrangian mechanics: either using the discrete Lagrangian to define a finite version of Hamilton's action principle, or treating it as a symplectic generating function. This is demonstrated for a discrete Lagrangian defined on an arbitrary Lie groupoid; the often encountered special case of the pair groupoid (or Cartesian square) is also given as a worked example.