Researcher profile

Benjamin Seibold

Benjamin Seibold contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
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

4 published item(s)

preprint2020arXiv

Massively Parallel Stencil Strategies for Radiation Transport Moment Model Simulations

The radiation transport equation is a mesoscopic equation in high dimensional phase space. Moment methods approximate it via a system of partial differential equations in traditional space-time. One challenge is the high computational intensity due to large vector sizes (1600 components for P39) in each spatial grid point. In this work, we extend the calculable domain size in 3D simulations considerably, by implementing the StaRMAP methodology within the massively parallel HPC framework NAStJA, which is designed to use current supercomputers efficiently. We apply several optimization techniques, including a new memory layout and explicit SIMD vectorization. We showcase a simulation with 200 billion degrees of freedom, and argue how the implementations can be extended and used in many scientific domains.

preprint2020arXiv

Quantitative Relationship between Population Mobility and COVID-19 Growth Rate based on 14 Countries

This study develops a framework for quantification of the impact of changes in population mobility due to social distancing on the COVID-19 infection growth rate. Using the Susceptible-Infected-Recovered (SIR) epidemiological model we establish that under some mild assumptions the growth rate of COVID-19 deaths is a time-delayed approximation of the growth rate of COVID-19 infections. We then hypothesize that the growth rate of COVID-19 infections is a function of population mobility, which leads to a statistical model that predicts the growth rate of COVID-19 deaths as a delayed function of population mobility. The parameters of the statistical model directly reveal the growth rate of infections, the mobility-dependent transmission rate, the mobility-independent recovery rate, and the critical mobility, below which COVID-19 growth rate becomes negative. We fitted the proposed statistical model on publicly available data from 14 countries where daily death counts exceeded 100 for more than 3 days as of May 6th, 2020. The publicly available Google Mobility Index (GMI) was used as a measure of population mobility at the country level. Our results show that the growth rate of COVID-19 deaths can be accurately estimated 20 days ahead as a quadratic function of the transit category of GMI (adjusted R-squared = 0.784). The estimated 95% confidence interval for the critical mobility is in the range between 36.1% and 47.6% of the pre-COVID-19 mobility. This result indicates that a significant reduction in population mobility is needed to reverse the growth of COVID-19 epidemic. Moreover, the quantitative relationship established herein suggests that a readily available, population-level metric such as GMI can be a useful indicator of the course of COVID-19 epidemic.

preprint2019arXiv

Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts

For a wide range of phenomena, current computational ability does not always allow for fully atomistic simulations of high-dimensional molecular systems to reach time scales of interest. Coarse-graining (CG) is an established approach to alleviate the impact of computational limits while retaining the same algorithms used in atomistic simulations. It is of importance to understand how algorithms such as Langevin integrators perform on non-trivial CG molecular systems, and in particular how large of an integration time step can be used without introducing unacceptable amounts of error into averaged quantities of interest. To investigate this, we examined three different Langevin integrators on a CG polymer melt: the recently developed BAOAB method by Leimkuhler and Matthews, the Gronbech-Jensen and Farago method, or G-JF, and the frequently used Brunger-Brooks-Karplus integrator, also known as BBK. We compute and analyze key statistical properties for each. Our results indicate that the three integrators perform similarly when using a small friction parameter; however, outside of this regime the use of large integration steps produces significant deviations from the predicted diffusivity and steady-state distributions for all integration methods examined with the exception of G-JF.

preprint2009arXiv

A rarefaction-tracking method for hyperbolic conservation laws

We present a numerical method for scalar conservation laws in one space dimension. The solution is approximated by local similarity solutions. While many commonly used approaches are based on shocks, the presented method uses rarefaction and compression waves. The solution is represented by particles that carry function values and move according to the method of characteristics. Between two neighboring particles, an interpolation is defined by an analytical similarity solution of the conservation law. An interaction of particles represents a collision of characteristics. The resulting shock is resolved by merging particles so that the total area under the function is conserved. The method is variation diminishing, nevertheless, it has no numerical dissipation away from shocks. Although shocks are not explicitly tracked, they can be located accurately. We present numerical examples, and outline specific applications and extensions of the approach.