Source author record

Gregory Herschlag

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

ResearcherUnclaimed source record

math.PR physics.flu-dyn Computation math.NA Numerical Analysis physics.chem-ph

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

Multi-Scale Merge-Split Markov Chain Monte Carlo for Redistricting

We develop a Multi-Scale Merge-Split Markov chain on redistricting plans. The chain is designed to be usable as the proposal in a Markov Chain Monte Carlo (MCMC) algorithm. Sampling the space of plans amounts to dividing a graph into a partition with a specified number of elements which each correspond to a different district. The districts satisfy a collection of hard constraints and the measure may be weighted with regard to a number of other criteria. The multi-scale algorithm is similar to our previously developed Merge-Split proposal, however, this algorithm provides improved scaling properties and may also be used to preserve nested communities of interest such as counties and precincts. Both works use a proposal which extends the ReCom algorithm which leveraged spanning trees merge and split districts. In this work we extend the state space so that each district is defined by a hierarchy of trees. In this sense, the proposal step in both algorithms can be seen as a "Forest ReCom." We also expand the state space to include edges that link specified districts, which further improves the computational efficiency of our algorithm. The collection of plans sampled by the MCMC algorithm can serve as a baseline against which a particular plan of interest is compared. If a given plan has different racial or partisan qualities than what is typical of the collection of plans, the given plan may have been gerrymandered and is labeled as an outlier.

preprint2020arXiv

Non-reversible Markov chain Monte Carlo for sampling of districting maps

Evaluating the degree of partisan districting (Gerrymandering) in a statistical framework typically requires an ensemble of districting plans which are drawn from a prescribed probability distribution that adheres to a realistic and non-partisan criteria. In this article we introduce novel non-reversible Markov chain Monte-Carlo (MCMC) methods for the sampling of such districting plans which have improved mixing properties in comparison to previously used (reversible) MCMC algorithms. In doing so we extend the current framework for construction of non-reversible Markov chains on discrete sampling spaces by considering a generalization of skew detailed balance. We provide a detailed description of the proposed algorithms and evaluate their performance in numerical experiments.

preprint2015arXiv

A consistent hierarchy of generalized kinetic equation approximations to the chemical master equation applied to surface catalysis

We develop a hierarchy of approximations to the master equation for systems that exhibit translational invariance and finite-range spatial correlation. Each approximation within the hierarchy is a set of ordinary differential equations that considers spatial correlations of varying lattice distance; the assumption is that the full system will have finite spatial correlations and thus the behavior of the models within the hierarchy will approach that of the full system. We provide evidence of this convergence in the context of one- and two-dimensional numerical examples. Lower levels within the hierarchy that consider shorter spatial correlations, are shown to be up to three orders of magnitude faster than traditional kinetic Monte Carlo methods (KMC) for one-dimensional systems, while predicting similar system dynamics and steady states as KMC methods. We then test the hierarchy on a two-dimensional model for the oxidation of CO on RuO2(110), showing that low-order truncations of the hierarchy efficiently capture the essential system dynamics. By considering sequences of models in the hierarchy that account for longer spatial correlations, successive model predictions may be used to establish empirical approximation of error estimates. The hierarchy may be thought of as a class of generalized phenomenological kinetic models since each element of the hierarchy approximates the master equation and the lowest level in the hierarchy is identical to a simple existing phenomenological kinetic models.

preprint2015arXiv

Optimal reservoir conditions for fluid extraction through permeable walls in the viscous limit

In biological transport mechanisms such as insect respiration and renal filtration, fluid travels along a leaky channel allowing exchange with systems exterior the the channel. The channels in these systems may undergo peristaltic pumping which is thought to enhance the material exchange. To date, little analytic work has been done to study the effect of pumping on material extraction across the channel walls. In this paper, we examine a fluid extraction model in which fluid flowing through a leaky channel is exchanged with fluid in a reservoir. The channel walls are allowed to contract and expand uniformly, simulating a pumping mechanism. In order to efficiently determine solutions of the model, we derive a formal power series solution for the Stokes equations in a finite channel with uniformly contracting/expanding permeable walls. This flow has been well studied in the case of weakly permeable channel walls in which the normal velocity at the channel walls is proportional to the wall velocity. In contrast we do not assume weakly driven flow, but flow driven by hydrostatic pressure, and we use Dacry's law to close our system for normal wall velocity. We use our flow solution to examine flux across the channel-reservoir barrier and demonstrate that pumping can either enhance or impede fluid extraction across channel walls. We find that associated with each set of physical flow and pumping parameters, there are optimal reservoir conditions that maximizes the amount of material flowing from the channel into the reservoir.

preprint2010arXiv

Reynolds number limits for jet propulsion: A numerical study of simplified jellyfish

The Scallop Theorem states that reciprocal methods of locomotion, such as jet propulsion or paddling, will not work in Stokes flow (Reynolds number = 0). In nature the effective limit of jet propulsion is still in the range where inertial forces are significant. It appears that almost all animals that use jet propulsion swim at Reynolds numbers (Re) of about 5 or more. Juvenile squid and octopods hatch from the egg already swimming in this inertial regime. The limitations of jet propulsion at intermediate Re is explored here using the immersed boundary method to solve the two-dimensional Navier Stokes equations coupled to the motion of a simplified jellyfish. The contraction and expansion kinematics are prescribed, but the forward and backward swimming motions of the idealized jellyfish are emergent properties determined by the resulting fluid dynamics. Simulations are performed for both an oblate bell shape using a paddling mode of swimming and a prolate bell shape using jet propulsion. Average forward velocities and work put into the system are calculated for Reynolds numbers between 1 and 320. The results show that forward velocities rapidly decay with decreasing Re for all bell shapes when Re < 10. Similarly, the work required to generate the pulsing motion increases significantly for Re < 10. When compared actual organisms, the swimming velocities and vortex separation patterns for the model prolate agree with those observed in Nemopsis bachei. The forward swimming velocities of the model oblate jellyfish after two pulse cycles are comparable to those reported for Aurelia aurita, but discrepancies are observed in the vortex dynamics between when the 2D model oblate jellyfish and the organism.