Researcher profile

Joshua C. Dolence

Joshua C. Dolence contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
5topics
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)

preprint2022arXiv

Not-Quite Transcendental Functions and their Applications

Transcendental functions, such as exponentials and logarithms, appear in a broad array of computational domains: from simulations in curvilinear coordinates, to interpolation, to machine learning. Unfortunately they are typically expensive to compute accurately. In this note, we argue that in many cases, the properties of the function matters more than the exact functional form. We present new functions, which are not transcendental, that can be used as drop-in replacements for the exponential and logarithm in many settings for a significant performance boost. We show that for certain applications using these functions result in no drop in the accuracy at all, as they are perfectly accurate representations of themselves, if not the original transcendental functions.

preprint2022arXiv

PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion

The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.

preprint2019arXiv

Convection-Aided Explosions in One-Dimensional Core-Collapse Supernova Simulations I: Technique and Validation

Most one-dimensional core-collapse simulations fail to explode, yet multi-dimensional simulations often explode. A dominant multi-dimensional effect aiding explosion is neutrino-driven convection. We incorporate a convection model in approximate one-dimensional core-collapse supernova (CCSN) simulations. This is the 1D+ method. This convection model lowers the neutrino luminosity required for explosion by 30%, similar to the reduction observed in multi-dimensional simulations. The model is based upon the global turbulence model of Mabanta & Murphy (2018) and models the mean-field turbulent flow of neutrino-driven convection. In this preliminary investigation, we use simple neutrino heating and cooling algorithms to compare the critical condition in the 1D+ simulations with the critical condition observed in two-dimensional simulations. Qualitatively, the critical conditions in the 1D+ and the two-dimensional simulations are similar. The assumptions in the convection model affect the radial profiles of density, entropy, and temperature, and comparisons with the profiles of three dimensional simulations will help to calibrate these assumptions. These 1D+ simulations are consistent with the profiles and explosion conditions of equivalent two-dimensional CCSN simulations but are ~100 times faster, and the 1D+ prescription has the potential to be ~100,000 faster than three-dimensional CCSN simulations. The 1D+ technique will be ideally suited to test the explodability of thousands of progenitor models.

preprint2019arXiv

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

We present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress-energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Radiation four-fources are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way, with Compton scattering further approximated by the Kompaneets equation. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that MOCMC converges at least as $\sim N^{-1}$, rather than the canonical Monte Carlo $N^{-1/2}$, where $N$ is the number of samples per zone. On several problems we demonstrate substantial improvement over Eddington and M1 closures and gray opacities.