-
Filtering After Shading With Stochastic Texture Filtering
Authors:
Matt Pharr,
Bartlomiej Wronski,
Marco Salvi,
Marcos Fajardo
Abstract:
2D texture maps and 3D voxel arrays are widely used to add rich detail to the surfaces and volumes of rendered scenes, and filtered texture lookups are integral to producing high-quality imagery. We show that applying the texture filter after evaluating shading generally gives more accurate imagery than filtering textures before BSDF evaluation, as is current practice. These benefits are not merel…
▽ More
2D texture maps and 3D voxel arrays are widely used to add rich detail to the surfaces and volumes of rendered scenes, and filtered texture lookups are integral to producing high-quality imagery. We show that applying the texture filter after evaluating shading generally gives more accurate imagery than filtering textures before BSDF evaluation, as is current practice. These benefits are not merely theoretical, but are apparent in common cases. We demonstrate that practical and efficient filtering after shading is possible through the use of stochastic sampling of texture filters.
Stochastic texture filtering offers additional benefits, including efficient implementation of high-quality texture filters and efficient filtering of textures stored in compressed and sparse data structures, including neural representations. We demonstrate applications in both real-time and offline rendering and show that the additional error from stochastic filtering is minimal. We find that this error is handled well by either spatiotemporal denoising or moderate pixel sampling rates.
△ Less
Submitted 14 May, 2024;
originally announced July 2024.
-
Error Field Predictability and Consequences for ITER
Authors:
Matthew Pharr,
Nikolas Logan,
Carlos Paz-Soldan,
Jong-Kyu Park,
Chistopher Hansen
Abstract:
ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils…
▽ More
ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scaling laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set.
△ Less
Submitted 9 June, 2024; v1 submitted 3 June, 2024;
originally announced June 2024.
-
MANTA: A Negative-Triangularity NASEM-Compliant Fusion Pilot Plant
Authors:
MANTA Collaboration,
G. Rutherford,
H. S. Wilson,
A. Saltzman,
D. Arnold,
J. L. Ball,
S. Benjamin,
R. Bielajew,
N. de Boucaud,
M. Calvo-Carrera,
R. Chandra,
H. Choudhury,
C. Cummings,
L. Corsaro,
N. DaSilva,
R. Diab,
A. R. Devitre,
S. Ferry,
S. J. Frank,
C. J. Hansen,
J. Jerkins,
J. D. Johnson,
P. Lunia,
J. van de Lindt,
S. Mackie
, et al. (16 additional authors not shown)
Abstract:
The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicin…
▽ More
The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Two-Shot Optimization of Compositionally Complex Refractory Alloys
Authors:
James D. Paramore,
Brady G. Butler,
Michael T. Hurst,
Trevor Hastings,
Daniel O. Lewis,
Eli Norris,
Benjamin Barkai,
Joshua Cline,
Braden Miller,
Jose Cortes,
Ibrahim Karaman,
George M. Pharr,
Raymundo Arroyave
Abstract:
In this paper, a synergistic computational/experimental approach is presented for the rapid discovery and characterization of novel alloys within the compositionally complex (i.e., "medium/high entropy") refractory alloy space of Ti-V-Nb-Mo-Hf-Ta-W. This was demonstrated via a material design cycle aimed at simultaneously maximizing the objective properties of high specific hardness (hardness norm…
▽ More
In this paper, a synergistic computational/experimental approach is presented for the rapid discovery and characterization of novel alloys within the compositionally complex (i.e., "medium/high entropy") refractory alloy space of Ti-V-Nb-Mo-Hf-Ta-W. This was demonstrated via a material design cycle aimed at simultaneously maximizing the objective properties of high specific hardness (hardness normalized by density) and high specific elastic modulus (elastic modulus normalized by density). This framework utilizes high-throughput computational thermodynamics and intelligent filtering to first reduce the untenably large alloy space to a feasible size, followed by an iterative design cycle comprised of high-throughput synthesis, processing, and characterization in batch sizes of 24 alloys. After the first iteration, Bayesian optimization was utilized to inform selection of the next batch of 24 alloys. This paper demonstrates the benefit of using batch Bayesian optimization (BBO) in material design, as significant gains in the objective properties were observed after only two iterations or "shots" of the design cycle without using any prior knowledge or physical models of how the objective properties relate to the design inputs (i.e., composition). Specifically, the hypervolume of the Pareto front increased by 54% between the first and second iterations. Furthermore, 10 of the 24 alloys in the second iteration dominated all alloys from the first iteration.
△ Less
Submitted 11 May, 2024;
originally announced May 2024.
-
ART-Owen Scrambling
Authors:
Abdalla G. M. Ahmed,
Matt Pharr,
Peter Wonka
Abstract:
We present a novel algorithm for implementing Owen-scrambling, combining the generation and distribution of the scrambling bits in a single self-contained compact process. We employ a context-free grammar to build a binary tree of symbols, and equip each symbol with a scrambling code that affects all descendant nodes. We nominate the grammar of adaptive regular tiles (ART) derived from the repetit…
▽ More
We present a novel algorithm for implementing Owen-scrambling, combining the generation and distribution of the scrambling bits in a single self-contained compact process. We employ a context-free grammar to build a binary tree of symbols, and equip each symbol with a scrambling code that affects all descendant nodes. We nominate the grammar of adaptive regular tiles (ART) derived from the repetition-avoiding Thue-Morse word, and we discuss its potential advantages and shortcomings. Our algorithm has many advantages, including random access to samples, fixed time complexity, GPU friendliness, and scalability to any memory budget. Further, it provides two unique features over known methods: it admits optimization, and it is invertible, enabling screen-space scrambling of the high-dimensional Sobol sampler.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices
Authors:
C. Hansen,
I. G. Stewart,
D. Burgess,
M. Pharr,
S. Guizzo,
F. Logak,
A. O. Nelson,
C. Paz-Soldan
Abstract:
In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like description…
▽ More
In this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial orders for analytic and cross-code test cases.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
Rapid Quantification of Dynamic and Spall Strength of Metals Across Strain Rates
Authors:
Suhas Eswarappa Prameela,
Christopher C. Walker,
Christopher S. DiMarco,
Debjoy D. Mallick,
Xingsheng Sun,
Stephanie Hernandez,
Taisuke Sasaki,
Justin W. Wilkerson,
K. T. Ramesh,
George M. Pharr,
Timothy P. Weihs
Abstract:
The response of metals and their microstructures under extreme dynamic conditions can be markedly different from that under quasistatic conditions. Traditionally, high strain rates and shock stresses are measured using cumbersome and expensive methods such as the Kolsky bar or large spall experiments. These methods are low throughput and do not facilitate high-fidelity microstructure-property link…
▽ More
The response of metals and their microstructures under extreme dynamic conditions can be markedly different from that under quasistatic conditions. Traditionally, high strain rates and shock stresses are measured using cumbersome and expensive methods such as the Kolsky bar or large spall experiments. These methods are low throughput and do not facilitate high-fidelity microstructure-property linkages. In this work, we combine two powerful small-scale testing methods, custom nanoindentation, and laser-driven micro-flyer shock, to measure the dynamic and spall strength of metals. The nanoindentation system is configured to test samples from quasistatic to dynamic strain rate regimes (10$^{-3}$ s$^{-1}$ to 10$^{+4}$ s$^{-1}$). The laser-driven micro-flyer shock system can test samples through impact loading between 10$^{+5}$ s$^{-1}$ to 10$^{+7}$ s$^{-1}$ strain rates, triggering spall failure. The model material used for testing is Magnesium alloys, which are lightweight, possess high-specific strengths and have historically been challenging to design and strengthen due to their mechanical anisotropy. Here, we modulate their microstructure by adding or removing precipitates to demonstrate interesting upticks in strain rate sensitivity and evolution of dynamic strength. At high shock loading rates, we unravel an interesting paradigm where the spall strength of these materials converges, but the failure mechanisms are markedly different. Peak aging, considered to be a standard method to strengthen metallic alloys, causes catastrophic failure, faring much worse than solutionized alloys. Our high throughput testing framework not only quantifies strength but also teases out unexplored failure mechanisms at extreme strain rates, providing valuable insights for the rapid design and improvement of metals for extreme environments.
△ Less
Submitted 25 September, 2023;
originally announced September 2023.
-
Simultaneous access to high normalized current, pressure, density, and confinement in strongly-shaped diverted negative triangularity plasmas
Authors:
C. Paz-Soldan,
C. Chrystal,
P. Lunia,
A. O. Nelson,
K. E. Thome,
M. E. Austin,
T. B. Cote,
A. W. Hyatt,
A. Marinoni,
T. H. Osborne,
M. Pharr,
O. Sauter,
F. Scotti,
T. M. Wilks,
H. S. Wilson
Abstract:
Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized current, pressure, density, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangu…
▽ More
Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized current, pressure, density, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the operating space, edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are simultaneously observed, all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
Stochastic Texture Filtering
Authors:
Marcos Fajardo,
Bartlomiej Wronski,
Marco Salvi,
Matt Pharr
Abstract:
2D texture maps and 3D voxel arrays are widely used to add rich detail to the surfaces and volumes of rendered scenes, and filtered texture lookups are integral to producing high-quality imagery. We show that filtering textures after evaluating lighting, rather than before BSDF evaluation as is current practice, gives a more accurate solution to the rendering equation. These benefits are not merel…
▽ More
2D texture maps and 3D voxel arrays are widely used to add rich detail to the surfaces and volumes of rendered scenes, and filtered texture lookups are integral to producing high-quality imagery. We show that filtering textures after evaluating lighting, rather than before BSDF evaluation as is current practice, gives a more accurate solution to the rendering equation. These benefits are not merely theoretical, but are apparent in common cases. We further show that stochastically sampling texture filters is crucial for enabling this approach, which has not been possible previously except in limited cases. Stochastic texture filtering offers additional benefits, including efficient implementation of high-quality texture filters and efficient filtering of textures stored in compressed and sparse data structures, including neural representations. We demonstrate applications in both real-time and offline rendering and show that the additional stochastic error is minimal. Furthermore, this error is handled well by either spatiotemporal denoising or moderate pixel sampling rates.
△ Less
Submitted 15 May, 2023; v1 submitted 9 May, 2023;
originally announced May 2023.
-
A Generalized Ray Formulation For Wave-Optics Rendering
Authors:
Shlomi Steinberg,
Ravi Ramamoorthi,
Benedikt Bitterli,
Eugene d'Eon,
Ling-Qi Yan,
Matt Pharr
Abstract:
Under ray-optical light transport, the classical ray serves as a linear and local "point query" of light's behaviour. Linearity and locality are crucial to the formulation of sophisticated path tracing and sampling techniques, that enable efficient solutions to light transport problems in complex, real-world settings and environments. However, such formulations are firmly confined to the realm of…
▽ More
Under ray-optical light transport, the classical ray serves as a linear and local "point query" of light's behaviour. Linearity and locality are crucial to the formulation of sophisticated path tracing and sampling techniques, that enable efficient solutions to light transport problems in complex, real-world settings and environments. However, such formulations are firmly confined to the realm of ray optics, while many applications of interest -- in computer graphics and computational optics -- demand a more precise understanding of light: as waves. We rigorously formulate the generalized ray, which enables linear and weakly-local queries of arbitrary wave-optical distributions of light. Generalized rays arise from photodetection states, and therefore allow performing backward (sensor-to-source) wave-optical light transport. Our formulations are accurate and highly general: they facilitate the application of modern path tracing techniques for wave-optical rendering, with light of any state of coherence and any spectral properties. We improve upon the state-of-the-art in terms of the generality and accuracy of the formalism, ease of application, as well as performance. As a consequence, we are able to render large, complex scenes, as in Fig. 1, and even do interactive wave-optical light transport, none of which is possible with any existing method. We numerically validate our formalism, and make connection to partially-coherent light transport.
△ Less
Submitted 7 January, 2024; v1 submitted 28 March, 2023;
originally announced March 2023.
-
Decorrelating ReSTIR Samplers via MCMC Mutations
Authors:
Rohan Sawhney,
Daqi Lin,
Markus Kettunen,
Benedikt Bitterli,
Ravi Ramamoorthi,
Chris Wyman,
Matt Pharr
Abstract:
Monte Carlo rendering algorithms often utilize correlations between pixels to improve efficiency and enhance image quality. For real-time applications in particular, repeated reservoir resampling offers a powerful framework to reuse samples both spatially in an image and temporally across multiple frames. While such techniques achieve equal-error up to 100 times faster for real-time direct lightin…
▽ More
Monte Carlo rendering algorithms often utilize correlations between pixels to improve efficiency and enhance image quality. For real-time applications in particular, repeated reservoir resampling offers a powerful framework to reuse samples both spatially in an image and temporally across multiple frames. While such techniques achieve equal-error up to 100 times faster for real-time direct lighting and global illumination, they are still far from optimal. For instance, unchecked spatiotemporal resampling often introduces noticeable correlation artifacts, while reservoirs holding more than one sample suffer from impoverishment in the form of duplicate samples. We demonstrate how interleaving Markov Chain Monte Carlo (MCMC) mutations with reservoir resampling helps alleviate these issues, especially in scenes with glossy materials and difficult-to-sample lighting. Moreover, our approach does not introduce any bias, and in practice we find considerable improvement in image quality with just a single mutation per reservoir sample in each frame.
△ Less
Submitted 31 October, 2022;
originally announced November 2022.
-
A non-local magneto-curvature instability in a differentially rotating disk
Authors:
Fatima Ebrahimi,
Matthew Pharr
Abstract:
A global mode is shown to be unstable to non-axisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic field where local sta…
▽ More
A global mode is shown to be unstable to non-axisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic field where local standard MRI becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén continuum driven unstable modes. As its mode structure and relative dominance over MRI is inherently determined by the global spatial curvature as well as the flow shear in the presence of magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from a MRI-driven turbulence to a state dominated by global non-axisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at higher magnetic field than predicted by local models.
△ Less
Submitted 10 August, 2022; v1 submitted 28 June, 2022;
originally announced June 2022.
-
The contribution of plastic sink-in to the static friction of single asperity microscopic contacts
Authors:
Owen Brazil,
John B. Pethica,
George M. Pharr
Abstract:
We report microscale friction experiments for diamond/metal and diamond/silica contacts under gigapascal contact pressures. Using a new nanoprobe technique which has sufficient dynamic range of force and stiffness, we demonstrate the processes involved in the transition from purely interface sliding at the nanoscale to the situation where at least one of the sliding bodies undergoes some plastic d…
▽ More
We report microscale friction experiments for diamond/metal and diamond/silica contacts under gigapascal contact pressures. Using a new nanoprobe technique which has sufficient dynamic range of force and stiffness, we demonstrate the processes involved in the transition from purely interface sliding at the nanoscale to the situation where at least one of the sliding bodies undergoes some plastic deformation. For sliding of micrometer-scale diamond spherical tips on metallic substrates, additional local plastic yielding of the substrate resulting from tangential tractions causes the tip to sink into the surface, increasing the contact area in the direction of loading and resulting in a static friction coefficient higher than the kinetic during ploughing. This sink-in is largely absent in fused silica, and no friction drop is observed, along with lower friction in general. The transition from sinking in within the static friction regime to ploughing in the sliding friction regime is mediated by failure of the contact interface, indicated by a sharp increase in energy dissipation. At lower contact pressures, the elastic interfacial sliding behaviour characteristic of scanning probe or surface force apparatus experiments is recovered, bridging the gap between the exotic realm of nanotribology and plasticity dominated macroscale friction. We delineate the material and geometric factors which determine the transition. We also find that some unexpected light is cast on the origins of the difference between static and dynamic friction.
△ Less
Submitted 22 June, 2021;
originally announced June 2021.
-
A System for Acquiring, Processing, and Rendering Panoramic Light Field Stills for Virtual Reality
Authors:
Ryan S. Overbeck,
Daniel Erickson,
Daniel Evangelakos,
Matt Pharr,
Paul Debevec
Abstract:
We present a system for acquiring, processing, and rendering panoramic light field still photography for display in Virtual Reality (VR). We acquire spherical light field datasets with two novel light field camera rigs designed for portable and efficient light field acquisition. We introduce a novel real-time light field reconstruction algorithm that uses a per-view geometry and a disk-based blend…
▽ More
We present a system for acquiring, processing, and rendering panoramic light field still photography for display in Virtual Reality (VR). We acquire spherical light field datasets with two novel light field camera rigs designed for portable and efficient light field acquisition. We introduce a novel real-time light field reconstruction algorithm that uses a per-view geometry and a disk-based blending field. We also demonstrate how to use a light field prefiltering operation to project from a high-quality offline reconstruction model into our real-time model while suppressing artifacts. We introduce a practical approach for compressing light fields by modifying the VP9 video codec to provide high quality compression with real-time, random access decompression.
We combine these components into a complete light field system offering convenient acquisition, compact file size, and high-quality rendering while generating stereo views at 90Hz on commodity VR hardware. Using our system, we built a freely available light field experience application called Welcome to Light Fields featuring a library of panoramic light field stills for consumer VR which has been downloaded over 15,000 times.
△ Less
Submitted 20 October, 2018;
originally announced October 2018.