Source author record

Benjamin Seibold

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

11 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.

preprint2016arXiv

Relations between WENO3 and Third-order Limiting in Finite Volume Methods

Weighted essentially non-oscillatory (WENO) and finite volume (FV) methods employ different philosophies in their way to perform limiting. We show that a generalized view on limiter functions, which considers a two-dimensional, rather than a one-dimensional dependence on the slopes in neighboring cells, allows to write WENO3 and $3^\text{rd}$-order FV schemes in the same fashion. Within this framework, it becomes apparent that the classical approach of FV limiters to only consider ratios of the slopes in neighboring cells, is overly restrictive. The hope of this new perspective is to establish new connections between WENO3 and FV limiter functions, which may give rise to improvements for the limiting behavior in both approaches.

preprint2013arXiv

A comparison of data-fitted first order traffic models and their second order generalizations via trajectory and sensor data

The Aw-Rascle-Zhang (ARZ) model can be interpreted as a generalization of the first order Lighthill-Whitham-Richards (LWR) model, possessing a family of fundamental diagram curves, rather than a single one. We investigate to which extent this generalization increases the predictive accuracy of the models. To that end, a systematic comparison of two types of data-fitted LWR models and their second order ARZ counterparts is conducted, via a version of the three-detector problem test. The parameter functions of the models are constructed using historic fundamental diagram data. The model comparisons are then carried out using time-dependent data, of two very different types: vehicle trajectory data, and single-loop sensor data. The study of these PDE models is carried out in a macroscopic sense, i.e., continuous field quantities are constructed from the discrete data, and discretization effects are kept negligibly small.

preprint2013arXiv

Effect of the choice of stagnation density in data-fitted first- and second-order traffic models

For a class of data-fitted macroscopic traffic models, the influence of the choice of the stagnation density on the model accuracy is investigated. This work builds on an established framework of data-fitted first-order Lighthill-Whitham-Richards (LWR) models and their second-order Aw-Rascle-Zhang (ARZ) generalizations. These models are systematically fitted to historic fundamental diagram data, and then their predictive accuracy is quantified via a version of the three-detector problem test, considering vehicle trajectory data and single-loop sensor data. The key outcome of this study is that with commonly suggested stagnation densities of 120 vehicles/km/lane and above, information travels backwards too slowly. It is then demonstrated that the reduction of the stagnation density to 90-100 vehicles/km/lane addresses this problem and results in a significant improvement of the predictive accuracy of the considered models.

preprint2012arXiv

The sound of an evolving floating sculpture

Commissioned by MIT's in-house artist Jane Philbrick, we evolve an abstract 2D surface (resembling Marta Pan's 1961 "Sculpture Flottante I") under mean curvature, all the while calculating the eigenmodes and eigenvalues of the Laplace-Beltrami operator on the resulting shapes. These are then synthesized into a sound-wave embodying the "swan song" of the surfaces as the evolve to points and vanish. The surface is approximated by a triangulation, and we present a robust approach to approximate the normal directions and the mean curvature. The resulting video and sound-track were parts in the Jane Philbrick's exhibition "Everything Trembles" in Lund, Sweden, 2009.

preprint2009arXiv

A gradient-augmented level set method with an optimally local, coherent advection scheme

The level set approach represents surfaces implicitly, and advects them by evolving a level set function, which is numerically defined on an Eulerian grid. Here we present an approach that augments the level set function values by gradient information, and evolves both quantities in a fully coupled fashion. This maintains the coherence between function values and derivatives, while exploiting the extra information carried by the derivatives. The method is of comparable quality to WENO schemes, but with optimally local stencils (performing updates in time by using information from only a single adjacent grid cell). In addition, structures smaller than the grid size can be located and tracked, and the extra derivative information can be employed to obtain simple and accurate approximations to the curvature. We analyze the accuracy and the stability of the new scheme, and perform benchmark tests.

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.

preprint2009arXiv

An exactly conservative particle method for one dimensional scalar conservation laws

A particle scheme for scalar conservation laws in one space dimension is presented. Particles representing the solution are moved according to their characteristic velocities. Particle interaction is resolved locally, satisfying exact conservation of area. Shocks stay sharp and propagate at correct speeds, while rarefaction waves are created where appropriate. The method is variation diminishing, entropy decreasing, exactly conservative, and has no numerical dissipation away from shocks. Solutions, including the location of shocks, are approximated with second order accuracy. Source terms can be included. The method is compared to CLAWPACK in various examples, and found to yield a comparable or better accuracy for similar resolutions.

preprint2008arXiv

Self-sustained nonlinear waves in traffic flow

In analogy to gas-dynamical detonation waves, which consist of a shock with an attached exothermic reaction zone, we consider herein nonlinear traveling wave solutions, termed "jamitons," to the hyperbolic ("inviscid") continuum traffic equations. Generic existence criteria are examined in the context of the Lax entropy conditions. Our analysis naturally precludes traveling wave solutions for which the shocks travel downstream more rapidly than individual vehicles. Consistent with recent experimental observations from a periodic roadway (Sugiyama et al., New Journal of Physics, 10, 2008), our numerical calculations show that, under appropriate road conditions, jamitons are attracting solutions, with the time evolution of the system converging towards a jamiton-dominated configuration. Jamitons are characterized by a sharp increase in density over a relatively compact section of the roadway. Applications of our analysis to traffic modeling and control are examined by way of a detailed example.