Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Description”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Interlaboratory Study: Testing Reproducibility of Solid Biofuels Component Identification Using Reflected Light Microscopy

Considering global market trends and concerns about climate change and sustainability, increased biomass use for energy is expected to continue. As more diverse materials are being utilized to manufacture solid biomass fuels, it is critical to implement quality assessment methods to analyze these fuels thoroughly. One such method is reflected light microscopy (RLM), which has the potential to complement and enhance current standard testing, leading to improving fuel quality assessment and, ultimately, preventing avoidable air pollution. An interlaboratory study (ILS) was conducted to test the reproducibility of biomass fuels component identification using a reflected light microscopy technique. The exercise was conducted on thirty photomicrographs showing biomass and various undesired components (like plastics or mineral matter), which were purposely added (by the ILS organizers) to contaminate wood pellets and charcoal-based grilling fuels. Forty-six participants had various levels of difficulty identifying the marked components, and as a result, the percentage of correct answers ranged from 52.2 to 94.4%. Among the most difficult components to distinguish were petroleum products and inorganic matter. Various reasons led to the misidentification, including insufficient morphological descriptions of the components provided to participants, ambiguities of the nomenclature, limitations of the analytical and exercise method, and insufficient experience of the participants. Overall, the results indicate that RLM has the potential to enhance the quality assessment of biomass fuels. However, they also demonstrate that the petrographic classification used in this exercise requires further refinement before it can be standardized. While a new simplified classification of solid biomass fuels components was created as an outcome of this study, future research is necessary to refine the nomenclature, develop a microscopic morphological description of the components, and verify the accuracy of component identification with a follow-up ILS.

biomass↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape. The two variants of the hexagonal phases form a twin with respect to each other; (d) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

The three-dimensional (3D) quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam–plasma interactions. In this algorithm, the plasma response, i.e., plasma wave wake, to a non-evolving laser or particle beam is calculated using a set of Maxwell’s equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. Additionally, at the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the algorithmic complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r–z and a gridless description in . We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations is derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), a full 3D quasi-static code (QuickPIC), and the new quasi-static PIC code with azimuthal decomposition (QPAD) are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Efficacy of the symmetry-adapted basis for ab initio nucleon-nucleus interactions for light- and intermediate-mass nuclei

We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for 4 He, 16 O, and 20 Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasing model space sizes and symmetry-based selections on the SA-RGM norm and direct potential kernels, as well as on phase shifts, which are the input to calculations of cross sections. We demonstrate the efficacy of the SA basis and its scalability with particle numbers and model space dimensions, with a view toward ab initio descriptions of nucleon scattering and capture reactions up through the medium-mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, 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 degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.

03 NATURAL GAS↗
DOE OSTI2020

Gaseous jet through an outward opening injector: Details of mixing characteristic and turbulence scales

Direct injection (DI) strategy of natural gas (NG) into internal combustion engines (ICE) has led to higher thermal efficiency and lower exhaust emissions. In order to thoroughly understand the most relevant phenomena affecting the performances of such engines, computational fluid dynamics (CFD) plays a key role as an accurate description of the jet evolution and interaction within the combustion chamber is required to that aim. Accurate description of high-pressure gaseous jets is rather challenging at high Mach numbers, as the injected gas is strongly under-expanded once in the ambient, giving room to shocks due to compressibility effects. Also the interaction between shock waves and mixing layers needs to be carefully represented with a multi-dimensional model, calling for substantial computational resources requirements. In this paper a numerical investigation of the behavior of a gaseous jet (Argon) through an outward opening injector has been carried out. A Large Eddy Simulation (LES) approach has been used in order to track the structures derived by the interaction of the injected fuel with the surrounding ambient. Although already good results were obtained using a Reynolds Averaged Navier-Stokes (RANS) approach, the adoption of LES is required to characterize more accurately the jet properties in terms of vortex structures and mixing effectiveness. Finally, the effect of the Nozzle Pressure Ratio (NPR) on the jet evolution has been highlighted in the paper, showing how a higher NPR would give a faster injection process, compromising however the homogeneity of the mixture.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Critical evaluation and thermodynamic modeling of the Pd–Sn system

The single crystalline material PdSn 4 , a homologue of the Dirac nodal arc semimetal PtSn 4 , is a promising candidate to search for new topological states with fascinating quantum physical properties. In the present work, the thermochemistry and phase diagram information of the Pd–Sn system published in the literature was collected and critically reviewed, two thermodynamic reassessments of the Pd–Sn system were carried out in the frame of the CALPHAD approach to optimize the Gibbs free energy of each phase presented in the system. The liquid phase was described using respectively the Bragg-Williams random mixing model and the associate solution model. All the reliable thermodynamic property and phase equilibria data can be generally described using the presently obtained thermodynamic descriptions of the Pd–Sn system. Issues related to the further improvement of the thermodynamic descriptions of the Pd–Sn system are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2022

Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies

The standard Λ Cold Dark Matter (ΛCDM) cosmological model provides a good description of a wide range of astrophysical and cosmological data. However, there are a few big open questions that make the standard model look like an approximation to a more realistic scenario yet to be found. In this paper, we list a few important goals that need to be addressed in the next decade, taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant H 0 , the σ 8 – S 8 tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the 5.0 σ tension between the Planck CMB estimate of the Hubble constant H 0 and the SH0ES collaboration measurements. After showing the H 0 evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the Planck CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density Ω m , and the amplitude or rate of the growth of structure ( σ 8 , f σ 8 ). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the H 0 – S 8 tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals. Finally, we give an overview of upgraded experiments and next-generation space missions and facilities on Earth that will be of crucial importance to address all these open questions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Uranium Oxide Synthetic Pathway Discernment through Unsupervised Morphological Analysis

We present a novel unsupervised machine learning method for quantitative representation of scanning electron micrographs and its applications and performance for nuclear forensic analysis of uranium ore concentrates. The method uses a vector quantizing variational autoencoder followed by a histogram operation to encode a micrograph into a single dimensional representation, called the latent vector. The method requires no extant labeling of the data and can be applied over large datasets of micrographs with minimal human interaction. The representations generated are broadly descriptive of each micrograph and the microstructure of the material imaged. In the case of uranium ore concentrate analysis, the representations were amenable to processing reagent and ore concentrate species classification with accuracy of 81:8%, which is competitive with state-of-the-art supervised networks. The representations were also used to classify previously unseen processing routes, were able to classify imaging parameters such as magnification (to 76:0% accuracy), were able to classify fine grained process parameters such as calcining temperature (to 74:4% accuracy), and their informatic properties indicate that they are generally descriptive of the image represented. This method can be applied across microstructure analysis fields to perform quantitative analysis without the need for labor intensive and possibly biased human analysis.

Scanning Electron Microscopy, Vector Quantizing Va↗
DOE OSTI2022

Anisotropic temperature-dependent elastic constants and thermal conductivities of TRISO particle coatings

Tristructural isotropic (TRISO) particles show great promise as a candidate fuel for use in several next-generation high-temperature nuclear reactor designs due to their structural integrity and fuel performance at high temperatures and burnups. Computational codes exist that can simulate TRISO fuel performance characteristics and failure probabilities under extreme conditions which require knowledge of the TRISO coatings’ thermophysical properties. The thermophysical descriptions of the TRISO particle's layers (i.e., buffer, pyrolytic carbon, and silicon carbide) currently used in fuel performance codes, however, assume that many of these properties are constant with respect to temperature or texture. In this paper, we obtain the full elastic stiffness tensors of the carbon and silicon carbide layers, which have transversely isotropic symmetry. Additionally, the calculated elastic properties of each of the coating layers are in remarkable agreement with the current models at room temperature and correct orientations. Additionally, the textured 3C-SiC layer was found to exhibit novel auxetic behavior above 1500 °C. The anisotropic temperature-dependent thermal conductivities of the carbon layers were calculated using acoustical Grüneisen-Debye theory which are in excellent agreement with current models at room temperature and correct orientations. These texture- and temperature-dependent relationships can be incorporated into the thermophysical description of TRISO particles in order to more accurately model fuel performance and failure probabilities under extreme conditions in forthcoming high-fidelity computational simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The CLAS12 Data Acquisition System

The CLAS12 Data Acquisition System was designed and built as part of the CLAS12 detector project in Hall B at Jefferson Laboratory. This article contains a full description of the system, including requirements, design, hardware, and software descriptions, as well as the achieved performance. The associated computing, network, and slow controls systems are also described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The Dynamics of Ions on Phased Radio-frequency Carpets in High Pressure Gases and Application for Barium Tagging in Xenon Gas Time Projection Chambers

Radio-frequency (RF) carpets with ultra-fine pitches are examined for ion transport in gases at atmospheric pressures and above. We develop new analytic and computational methods for modeling RF ion transport at densities where dynamics are strongly influenced by buffer gas collisions. An analytic description of levitating and sweeping forces from phased arrays is obtained, then thermodynamic and kinetic principles are used to calculate ion loss rates in the presence of collisions. This methodology is validated against detailed microscopic SIMION simulations. We then explore a parameter space of special interest for neutrinoless double beta decay experiments: transport of barium ions in xenon at pressures from 1 to 10 bar. Our computations account for molecular ion formation and pressure dependent mobility as well as finite temperature effects. We discuss the challenges associated with achieving suitable operating conditions, which lie beyond the capabilities of existing devices, using presently available or near-future manufacturing techniques.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Burnable absorbers in nuclear reactors – A review

Burnable absorbers can benefit nuclear reactors of virtually any design by providing reactivity control for extended fuel cycles, tritium production, burning of long-lived radionuclides, and reactor safety. When selecting the ideal burnable absorber type and its design, one must consider the resulting impact on the reactor’s fuel cycle design and cost, reactivity, thermal hydraulics, manufacturing, and radiation response. These selection criteria, as well as neutronic and thermophysical material property requirements, may be vastly different depending on whether the burnable absorber is intended for use in a commercial water-cooled reactor, a research reactor, or a next-generation advanced reactor system. A recent integration, synthesis, description of past and present technologies, and identification of existing gaps and areas of future research is lacking on these important topics. Here, this paper includes a fundamental description of the use of burnable absorbers and their impacts on reactivity, absorber depletion, self-shielding, basic thermophysical properties, and the use of burnable absorbers in next-generation nuclear applications..

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Multiscale modeling and cinematic visualization of photosynthetic energy conversion processes from electronic to cell scales

We report conversion of sunlight into chemical energy, namely photosynthesis, is the primary energy source of life on Earth. A visualization depicting this process, based on multiscale computational models from electronic to cell scales, is presented in the form of an excerpt from the fulldome show Birth of Planet Earth. This accessible visual narrative shows a lay audience, including children, how the energy of sunlight is captured, converted, and stored through a chain of proteins to power living cells. The visualization is the result of a multi-year collaboration among biophysicists, visualization scientists, and artists, which, in turn, is based on a decade-long experimental-computational collaboration on structural and functional modeling that produced an atomic detail description of a bacterial bioenergetic organelle, the chromatophore. Software advancements necessitated by this project have led to significant performance and feature advances, including hardware-accelerated cinematic ray tracing and instanced visualizations for efficient cell-scale modeling. The energy conversion steps depicted feature an integration of function from electronic to cell levels, spanning nearly 12 orders of magnitude in time scales. This atomic detail description uniquely enables a modern retelling of one of humanity’s earliest stories—the interplay between light and life.

97 MATHEMATICS AND COMPUTING↗