Researcher profile

Henry O. Jacobs

Henry O. Jacobs contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
11topics
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

10 published item(s)

preprint2016arXiv

Qualitatively accurate spectral schemes for advection and transport

The transport and continuum equations exhibit a number of conservation laws. For example, scalar multiplication is conserved by the transport equation, while positivity of probabilities is conserved by the continuum equation. Certain discretization techniques, such as particle based methods, conserve these properties, but converge slower than spectral discretization methods on smooth data. Standard spectral discretization methods, on the other hand, do not conserve the invariants of the transport equation and the continuum equation. This article constructs a novel spectral discretization technique that conserves these important invariants while simultaneously preserving spectral convergence rates. The performance of this proposed method is illustrated on several numerical experiments.

preprint2016arXiv

The role of symmetry and dissipation in biolocomotion

In this paper we illustrate the potential role which relative limit cycles may play in biolocomotion. We do this by describing, in great detail, an elementary example of reduction of a lightly dissipative system modeling crawling-type locomotion in 3D. The symmetry group SE(2) is the set of rigid transformations of the horizontal (ground) plane. Given a time-periodic perturbation, the system will admit a relative limit cycle whereupon each period is related to the previous by a fixed translation and rotation along the ground. This toy model identifies how symmetry reduction and dissipation can conspire to create robust behavior in crawling, and possibly walking, locomotion.

preprint2015arXiv

Higher-order Spatial Accuracy in Diffeomorphic Image Registration

We discretize a cost functional for image registration problems by deriving Taylor expansions for the matching term. Minima of the discretized cost functionals can be computed with no spatial discretization error, and the optimal solutions are equivalent to minimal energy curves in the space of $k$-jets. We show that the solutions convergence to optimal solutions of the original cost functional as the number of particles increases with a convergence rate of $O(h^{d+k})$ where $h$ is a resolution parameter. The effect of this approach over traditional particle methods is illustrated on synthetic examples and real images.

preprint2014arXiv

How to stare at the higher-order n-dimensional chain rule without losing your marbles

Given two real functions on the real line f and g, the Faa di Bruno provides the higher order derivative of the composition of f and g, as a summation over the lower order derivatives of f and g individually. The corresponding multi-dimensional generalization is substantially more difficult due to the complicated combinatorial considerations one must take into account when dealing with standard multi-indices. In this note we provide a simple statement and derivation of the multi-dimensional Faa di Bruno formula by resorting to notions from multiset theory.

preprint2014arXiv

Lagrangian mechanics on centered semi-direct product

There exists two types of semi-direct products between a Lie group $G$ and a vector space $V$. The left semi-direct product, $G \ltimes V$, can be constructed when $G$ is equipped with a left action on $V$. Similarly, the right semi-direct product, $G \rtimes V$, can be constructed when $G$ is equipped with a right action on $V$. In this paper, we will construct a new type of semi-direct product, $G \Join V$, which can be seen as the `sum' of a right and left semi-direct product. We then parallel existing semi-direct product Euler-Poincaré theory. We find that the group multiplication, the Lie bracket, and the diamond operator can each be seen as a sum of the associated concepts in right and left semi-direct product theory. Finally, we conclude with a toy example and the group of 2-jets of diffeomorphisms above a fixed point. This final example has potential use in the creation of particle methods for problems on diffeomorphism groups.

preprint2014arXiv

Symmetry in Image Registration and Deformation Modeling

We survey the role of symmetry in diffeomorphic registration of landmarks, curves, surfaces, images and higher-order data. The infinite dimensional problem of finding correspondences between objects can for a range of concrete data types be reduced resulting in compact representations of shape and spatial structure. This reduction is possible because the available data is incomplete in encoding the full deformation model. Using reduction by symmetry, we describe the reduced models in a common theoretical framework that draws on links between the registration problem and geometric mechanics. Symmetry also arises in reduction to the Lie algebra using particle relabeling symmetry allowing the equations of motion to be written purely in terms of Eulerian velocity field. Reduction by symmetry has recently been applied for jet-matching and higher-order discrete approximations of the image matching problem. We outline these constructions and further cases where reduction by symmetry promises new approaches to registration of complex data types.

preprint2013arXiv

Fluid-structure interaction in the Lagrange-Poincare formalism: the Navier-Stokes and inviscid regimes

In this paper, we derive the equations of motion for an elastic body interacting with a perfect fluid via the framework of Lagrange-Poincare reduction. We model the combined fluid-structure system as a geodesic curve on the total space of a principal bundle on which a diffeomorphism group acts. After reduction by the diffeomorphism group we obtain the fluid-structure interactions where the fluid evolves by the inviscid fluid equations. Along the way, we describe various geometric structures appearing in fluid-structure interactions: principal connections, Lie groupoids, Lie algebroids, etc. We finish by introducing viscosity in our framework as an external force and adding the no-slip boundary condition. The result is a description of an elastic body immersed in a Navier-Stokes fluid as an externally forced Lagrange-Poincare equation. Expressing fluid-structure interactions with Lagrange-Poincare theory provides an alternative to the traditional description of the Navier-Stokes equations on an evolving domain.

preprint2013arXiv

On the coupling between an ideal fluid and immersed particles

In this paper we use Lagrange-Poincare reduction to understand the coupling between a fluid and a set of Lagrangian particles that are supposed to simulate it. In particular, we reinterpret the work of Cendra et al. by substituting velocity interpolation from particle velocities for their principal connection. The consequence of writing evolution equations in terms of interpolation is two-fold. First, it gives estimates on the error incurred when interpolation is used to derive the evolution of the system. Second, this form of the equations of motion can inspire a family of particle and hybrid particle-spectral methods where the error analysis is "built-in". We also discuss the influence of other parameters attached to the particles, such as shape, orientation, or higher-order deformations, and how they can help with conservation of momenta in the sense of Kelvin's circulation theorem.

preprint2013arXiv

The role of SE(d)-reduction for swimming in Stokes and Navier-Stokes fluids

Steady swimming appears both periodic and stable. These characteristics are the very definition of limit cycles, and so we ask "Can we view swimming as a limit cycle?" In this paper we will not be able to answer this question in full. However, we shall find that reduction by SE(d)-symmetry brings us closer. Upon performing reduction by symmetry, we will find a stable fixed point which corresponds to a motionless body in stagnant water. We will then speculate on the existence of periodic orbits which are "approximately" limit cycles in the reduced system. When we lift these periodic orbits from the reduced phase space, we obtain dynamically robust relatively periodic orbits wherein each period is related to the previous by an SE(d)-phase. Clearly, an SE(d) phase consisting of nonzero translation and identity rotation means directional swimming, while non-trivial rotations correspond to turning with a constant turning radius.

preprint2012arXiv

Swimming as a limit cycle

Steady swimming can be characterized as both periodic and stable. These characteristics are the very definition of limit cycles, and so we ask "Can we view swimming as a limit cycle?" In this paper we will find that the answer is "yes". We will define a class of dissipative systems which correspond to the passive dynamics of a body immersed in a Navier-Stokes fluid (i.e. the dynamics of a dead fish). Upon performing reduction by symmetry we will find a hyperbolically stable fixed point which corresponds to the stability of a dead fish in stagnant water. Given a periodic force on the shape of the body we will invoke the persistence theorem to assert the existence of a loop which approximately satisfies the exact equations of motion. If we lift this loop with a phase reconstruction formula we will find that the lifted loops are not loops, but stable trajectories which represent regular periodic motion reminiscent of swimming.