Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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 757 records · Page 42

Use of Digital Real-Time Simulation and Optimization to Identify Maximum Real Power Injection on the Banshee Distribution Network: Preprint

This study investigates the hosting capacity of the Banshee Distribution Network by optimizing the real power injection at carefully selected Distributed Energy Resource (DER) locations. The analysis is conducted within the framework of power system operational constraints, including bus voltage ranges, thermal line ratings, and transformer loading limits. A Python-based Genetic Algorithm (GA), implemented using the PyGAD library, is employed to iteratively identify the optimal power injection configuration that maximizes network utilization while preserving system reliability. The methodology integrates a real-time simulation environment using the Real-Time Digital Simulator (RTDS), allowing high-fidelity evaluation of power flow and voltage behavior under each proposed injection scenario. By coupling the optimization algorithm with real-time simulation feedback, this approach ensures that both static and dynamic constraints are enforced during the evaluation process. The GA leverages evolutionary operators such as selection, crossover, and mutation to navigate the nonlinear search space efficiently. The results of the study delineate the feasible hosting capacity at three targeted buses, reflecting maximum real power levels that can be injected without causing voltage violations, transformer overloading, or line congestion. These findings provide a decision- support tool for distribution planners and utilities aiming to integrate higher penetrations of DERs in existing infrastructure. Additionally, the work lays the foundation for extending such optimization techniques to multi-objective formulations, including economic dispatch and reactive power coordination, in future studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

Performance of a CO 2 sorbent for indoor air cleaning applications: Effects of environmental conditions, sorbent aging, and adsorption of co-occurring formaldehyde

Indoor air cleaning systems that incorporate CO 2 sorbent materials enable HVAC load shifting and efficiency improvements. This study developed a bench-scale experimental system to evaluate the performance of a sorbent under controlled operation conditions. A thermostatic holder containing 3.15 g sorbent was connected to a manifold that delivered CO 2 -enriched air at a known temperature and relative humidity (RH). The air stream was also enriched with 0.8-2.1 ppm formaldehyde. The CO 2 concentration was monitored in real-time upstream and downstream of the sorbent, and integrated formaldehyde samples were collected at different times using DNPH-coated silica cartridges. Sorbent regeneration was carried out by circulating clean air in countercurrent. Almost 200 loading/regeneration cycles were performed in the span of 17 months, from which 104 were carried out at reference test conditions defined by loading with air at 25°C, 38% RH, and 1000 ppm CO 2 , and regenerating with air at 80°C, 3% RH and 400 ppm CO 2 . The working capacity decreased slightly from 43-44 mg CO 2 per g sorbent to 39-40 mg per g over the 17 months. The capacity increased with lower loading temperature (in the range 15-35°C) and higher regeneration temperature, between 40 and 80°C. The CO 2 capacity was not sensitive to the moisture content in the range 6-9 g/m 3 , and decreased slightly when dry air was used. Loading isothermal breakthrough curves were fitted to three simple adsorption models, verifying that pseudo-first-order kinetics appropriately describes the adsorption process. The model predicted that equilibrium capacities decreased with increasing temperature from 15 to 35°C, while adsorption rate constants slightly increased. The formaldehyde adsorption efficiency was 80%-99% in different cycles, corresponding to an average capacity of 86 ± 36 µg/g. Formaldehyde was not quantitatively released during regeneration, but its accumulation on the sorbent did not affect CO 2 adsorption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗

Extinction Monitoring of Pulsed Proton Beams Using FPGA-Based Peak Detection

The Mu2e experiment at Fermilab imposes stringent requirements on the elimination of out-of-time beam in its pulsed proton beam - a requirement known as "extinction". We present a method to measure the out-of-time particle rates to calculate the level of extinction in the inter-pulse gaps. The proposed method utilizes an array of quartz Cherenkov radiators and photomultiplier tubes to detect particles scattered from a vacuum chamber in the M4 transfer beamline at Fermilab.The measurement will employ a new μTCA-based FPGA system for data acquisition and signal processing, utilizing real-time peak detection algorithms to count scattered beam particles. By integrating data over many transfers, the time profile of the out-of-time beam will be resolved to fractional levels relative to that of the in-time beam. These results are compared with G4beamline simulations to validate models of beam transport, dynamics, and extinction, providing critical input for optimizing beam delivery to Mu2e.

Hensley, Ryan [UC, Davis]↗

First Search for Neutral Current Coherent Single-Photon Production in MicroBooNE

This article presents the first search for neutrino-induced neutral current coherent single-photon production (NC coherent 1$\gamma$). The search makes use of data from the MicroBooNE 85-tonne active volume liquid argon time projection chamber detector, situated in the Fermilab Booster Neutrino Beam (BNB), with an average neutrino energy of $\langle E_{\nu}\rangle \sim 0.8$ GeV. A targeted selection of candidate neutrino interactions with a single photon-like electromagnetic shower in the final state and no visible vertex activity was developed to search for the NC coherent 1$\gamma$ process, along with two auxiliary selections used to constrain the dominant background from NC$\pi^0$ production. With an integrated exposure of $6.87 \times 10^{20}$ protons on target delivered by the BNB, we set the world's first limit for this rare process, corresponding to an upper limit on the flux-averaged cross section of $\sigma<1.49 \times 10^{-41}\text{cm}^2$ at 90% C.L.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cyber-Physical Simulation of the Cold Startup of Solid Oxide Fuel Cell – Gas Turbine (SOFC-GT) Hybrid Systems

This work introduces experimental studies for the cold startup process (CPS) of the SOFC-GT hybrid system using the cyber-physical simulation approach. The physical gas turbine is coupled with a cyber-physical SOFC stack, which is represented using the integration of a real time dynamic SOFC model with physical components (e.g., pressure chamber, natural gas burner, etc.). Different ramp rates of the turbine speed were tested out during the startup processes. Bypass valves were also used to manipulate the airflow during SOFC-GT hybrid system start-up process. Different ramp rates enable the rapid start-up of the turbine to avoid surge and stall, meanwhile enable acceptable warm rate of the fuel cell stack without damaging the cell material. CPS can enable dynamic characterizations of highly integrated systems at lower cost.

Zhou, Nana↗

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗

A Novel Modeling Framework for Computationally Efficient and Accurate Real-Time Ensemble Flood Forecasting With Uncertainty Quantification

A novel modeling framework that simultaneously improves accuracy, predictability, and computational efficiency is presented. It embraces the benefits of three modeling techniques integrated together for the first time: surrogate modeling, parameter inference, and data assimilation. The use of polynomial chaos expansion (PCE) surrogates significantly decreases computational time. Parameter inference allows for model faster convergence, reduced uncertainty, and superior accuracy of simulated results. Ensemble Kalman filters assimilate errors that occur during forecasting. To examine the applicability and effectiveness of the integrated framework, we developed 18 approaches according to how surrogate models are constructed, what type of parameter distributions are used as model inputs, and whether model parameters are updated during the data assimilation procedure. We conclude that (1) PCE must be built over various forcing and flow conditions, and in contrast to previous studies, it does not need to be rebuilt at each time step; (2) model parameter specification that relies on constrained, posterior information of parameters (so-called Selected specification) can significantly improve forecasting performance and reduce uncertainty bounds compared to Random specification using prior information of parameters; and (3) no substantial differences in results exist between single and dual ensemble Kalman filters, but the latter better simulates flood peaks. The use of PCE effectively compensates for the computational load added by the parameter inference and data assimilation (up to ~80 times faster). Therefore, the presented approach contributes to a shift in modeling paradigm arguing that complex, high-fidelity hydrologic and hydraulic models should be increasingly adopted for real-time and ensemble flood forecasting.

54 ENVIRONMENTAL SCIENCES↗

Recent and future developments in pultrusion technology with consideration for curved geometries: A review

Herein this paper examined the current state and future developments in pultrusion with particular emphasis on its application in curved part manufacturing. The relationship between factors such as resin chemistry, fiber characteristics, and die geometry that influences the properties of pultruded product were highlighted. Moreover, the specific challenges associated with pultruding curved parts such as the complexities in achieving uniformity and structural integrity in such geometries were discussed. The review emphasized mold design, process improvement, adaptive control systems for precise resin impregnation and material selection to address these challenges. Additionally, the paper suggests the integration of real-time monitoring and data analytics as ways to enhance quality control during curved parts pultrusion. These advancements will help to overcome challenges specific to curved pultrusion and make the process more efficient. Other manufacturing techniques such as filament winding, thermoforming, pulforming were mentioned as alternatives to curved parts pultrusion. The review also explores pultruded variable curvature processes, highlighting some notable patents and article related to this subject matter. Production of pultruded variable curvature parts was seen as a key driver that can shape the future of pultrusion. Finally, the paper anticipates future trends, with sustainability, customization, integration of advanced materials, and development of techniques for pultrusion of composites parts.

42 ENGINEERING↗

First observations from the Kr multi-monochromatic x-ray imager for time and spatially resolved diagnosis of hot implosion cores

This paper presents initial findings from the recently deployed Kr multi-monochromatic x-ray imager (Kr MMI) at the Omega laser facility. The experiment focuses on exploring implosion dynamics in exploding pusher capsules at three distinct initial gas fill densities. Utilizing time-gated and spatially integrated measurements, core size, electron temperature (Te), and electron densities (ne) are extracted through the analysis of the spectral region encompassing the Kr He α and its satellite lines. A comprehensive spectral database, incorporating atomic kinetics, spectroscopic-quality radiation transport, and Stark broadened line shapes, has been developed for rigorous data analysis. These measurements underscore the utility of the new Kr MMI instrument, which, combined with sophisticated analysis techniques, enables the diagnosis of plasma conditions at Te>2000 eV, thereby extending the capabilities beyond the prior Ar MMI design. This is an important stepping stone for achieving time-gated and space-resolved diagnostics of electron temperature, electron density, and heat transport in high temperature implosion cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nanosecond Gated CMOS Camera (NSGCC) ICD (Rev. 2.1)

The Ultra-Fast X-ray Imager (UXI) program is an ongoing effort at Sandia National Laboratories to create high speed, multi-frame, time-gated Read Out Integrated Circuits (ROICs), and a corresponding suite of photodetectors to image a wide variety of High Energy Density (HED) physics experiments on both Sandia’s Z-Machine and LLNL’s National Ignition Facility (NIF). Several cameras have been designed over the length of the program; one of the most recent is the Icarus, which is an improvement on past imagers (Furi and Hippogriff). A second sensor that can be connected is the Daedalus sensor. The Icarus is a 1024 × 512-pixel array with either 25 μm or 8 µm spatial resolution containing four frames of storage per pixel and has improved timing generation and distribution components while achieved 2 ns time gating. The Daedalus sensor is also a 1024 x 512-pixel array with 25 µm special resolution containing three frames of storage per pixel and has an increased set of features for a wider variety of applications from interlacing of rows in each frame to configurability of all shutters. See Section 14 for details regarding the Icarus implementation of the firmware and Section 15 for details regarding the Daedalus implementation. Due to the unique test environments UXI sensors are targeted for, full custom hardware was required to physically mount an Icarus or Daedalus sensor, manage its various functions, and read out pixel data for transfer to a host computer. Beyond experimental functionality, the hardware also needed to accommodate sensor characterization requirements. Lawrence Livermore National Laboratory’s ‘Version 4.0 Board’ was the result of these efforts. It mounts all the components required to fully utilize the Icarus and Daedalus sensors including analog to digital converters to convert pixel data and various system voltages to digital form for readout and analysis, DAC channels for remote configuration of critical bias voltages, static random-access memories to buffer pixel data, RS422 and Gigabit Ethernet communications for remote access, and an FPGA to tie these components together. This document describes the FPGA electrical interfaces in detail to allow the reader a greater understanding of the device, and to facilitate implementation of custom software to control and manage it. The Version 4.0 Board is a continuation of the Nano-second Gated CMOS hardware design that retains much of the functionality of the Version 1.0 Board while adding features including a DAC instead of digital potentiometers, as well as sensors for pressure and radiation. The Version 4.0 board is intended for applications requiring tight form-factor enclosures. It is composed of two stacking boards; one holds the FPGA and regulators to power the various components of the board while the other contains the mating connector to the sensor, image-readoff ADCs, the DAC, and other components.

42 ENGINEERING↗

Ectomycorrhizal symbiosis prepares its host locally and systemically for abiotic cue signaling

Tree growth and survival are dependent on their ability to perceive signals, integrate them, and trigger timely and fitted molecular and growth responses. While ectomycorrhizal symbiosis is a predominant tree-microbe interaction in forest ecosystems, little is known about how and to what extent it helps trees cope with environmental changes. We hypothesized that the presence of Laccaria bicolor influences abiotic cue perception by Populus trichocarpa and the ensuing signaling cascade. We submitted ectomycorrhizal or non-ectomycorrhizal P. trichocarpa cuttings to short-term cessation of watering or ozone fumigation to focus on signaling networks before the onset of any physiological damage. Poplar gene expression, metabolite levels, and hormone levels were measured in several organs (roots, leaves, mycorrhizas) and integrated into networks. We discriminated the signal responses modified or maintained by ectomycorrhization. Ectomycorrhizas buffered hormonal changes in response to short-term environmental variations systemically prepared the root system for further fungal colonization and alleviated part of the root abscisic acid (ABA) signaling. Furthermore, the presence of ectomycorrhizas in the roots also modified the leaf multi-omics landscape and ozone responses, most likely through rewiring of the molecular drivers of photosynthesis and the calcium signaling pathway. In conclusion, P. trichocarpa-L. bicolor symbiosis results in a systemic remodeling of the host's signaling networks in response to abiotic changes. In addition, ectomycorrhizal, hormonal, metabolic, and transcriptomic blueprints are maintained in response to abiotic cues, suggesting that ectomycorrhizas are less responsive than non-mycorrhizal roots to abiotic challenges.

60 APPLIED LIFE SCIENCES↗

Implementation of McMurchie–Davidson Algorithm for Gaussian AO Integrals Suited for SIMD Processors

We report an implementation of the McMurchie− Davidson evaluation scheme for 1- and 2-particle Gaussian AO integrals designed for processors with Single Instruction Multiple Data (SIMD) instruction sets. Like in our recent MD implementation for graphical processing units (GPUs) [Asadchev, A.; Valeev, E. F.. J. Chem. Phys. 2024, 160, 244109.], variable-sized batches of shellsets of integrals are evaluated at a time. By optimizing for the floating point instruction throughput rather than minimizing the number of operations, this approach achieves up to 50% of the theoretical hardware peak FP64 performance for many common SIMD-equipped platforms (AVX2, AVX512, NEON), which translates to speedups of up to 30 over the state-of-the-art one-shellset-at-a-time implementation of Obara−Saika-type schemes in Libint for a variety of primitive and contracted integrals. As with our previous work, we rely on the standard C++ programming language such as the std::simd standard library feature to be included in the 2026 ISO C++ standard without any explicit code generation to keep the code base small and portable. The implementation is part of the open source LibintX library freely available at https://github.com/ValeevGroup/libintx.

Basis sets↗

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tracing U.S. fuel life-cycle greenhouse gas emissions in a multi-sector dynamics model using LC-GCAM

Model-based analysis of fuel pathways is essential for informing energy and environmental policy. Two major model types are typically used: multi-sector dynamics models, which capture the broader energy-economy, such as GCAM (Global Change Analysis Model), and life cycle assessment models, such as GREET (Greenhouse Gases, Regulated Emissions, and Energy Use in Transportation). Each has distinct strengths and limitations, and recent studies increasingly adopt hybrid approaches to harness the advantages of both. However, such integration is often time-consuming and complicated by inconsistencies in system boundaries and technology definitions. We present LC-GCAM, a new tool that enables estimation of life-cycle greenhouse gas emissions and primary energy use for any fuel pathway represented in GCAM. We apply LC-GCAM to 300 scenarios designed to explore key uncertainties affecting the life-cycle performance of future fuel options in the U.S. freight sector. To evaluate LC-GCAM, we compare its results with those from GREET for nine fuel types in a 2030 reference scenario. When input assumptions are modestly aligned, LC-GCAM and GREET estimates typically agree within 10% (absolute sum-based mean absolute percentage error). LC-GCAM offers a flexible and efficient approach to generating life-cycle metrics within an integrated modeling framework, supporting robust policy analysis across a wide range of interacting energy system uncertainties.

Wolfram, Paul↗