Search NASA⌕ Search

SEARCH · Search NASA

Results for “deterministic”

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 667 records · Page 37

Position-Dependency of Fuel Assembly Homogenization in a Pressurized Water Reactor

2-step nuclear reactor core analysis has long been established as the conventional calculational procedure for the pressurized water reactor (PWR) due to its efficiency. This is because the procedure significantly avoids the potentially huge computational burden arising from an explicit geometrical modeling and consideration of over ∼ 500 energy groups. It instead approaches the problem with spatial homogenization and group condensation: a 2-D lattice transport calculation to homogenize the fuel assembly cross-sections and a 3-D whole core nodal diffusion. One notes that the few-group constant generation using deterministic method could be complicated due to its limited flexibilities in geometrical modeling and inherent drawbacks of the multi-group approximation. Furthermore, accuracy of the deterministic few-group constant generation can especially be compromised during the resonance treatment calculation. As such, an alternative lattice calculation using the Monte Carlo (MC) approach, in place of the deterministic fuel assembly homogenization (few-group constant generation), has previously been proposed. During the assembly homogenization, directions and magnitudes of the current at boundaries are unknowns. Therefore, zero net current boundary conditions (infinite lattice) are instead used as the assembly node. The assumption is quite reasonable considering heterogeneities inside the node are more dominant than heterogeneities due to nodal location in the core. In fact, this approximation has successfully been applied to many LWR core analyses to reasonable accuracy. It should nonetheless be noted that the net current zero boundary assumption effectively ignores position-dependency of the fuel assembly homogenization in the general 2-step method. As such, there are a number of well-known researches on the extended assembly homogenization techniques, which introduce correction factors to compensate effects of the assembly environment, such as Generalized Equivalence Theory (GET), Super-homogenization method (SPH), and color set calculation method. Even though the assembly environment effect has long been acknowledged, its effect has not properly been investigated due to the limited computing power in the past and due to the reasonable accuracy of the extended assembly homogenization techniques. Nonetheless, quantification of position-dependency of the assembly homogenization in the core could still potentially offer additional value in terms of calculation accuracy to the 2-step reactor analysis. It is upon this motivation that the paper is prepared. In this research, position-dependent few-group constants were calculated and compared against those of infinite lattices. All calculations were completed using the MC Serpent2 code with ENDF/B7.0 nuclear data library. In this study, position-dependency of the fuel assembly homogenized cross sections in a PWR core has been investigated through comparison of the conventional 2-step reactor analysis with the reference high fidelity Monte Carlo simulation. It was found that the homogenized two-group constants are clearly position-dependent as the discrepancies between infinite lattice and whole-core calculations are quite significant especially in a rodded assembly. As such, the corresponding correction required due to the assembly environment effect is not negligible compared to the correction required due to assembly heterogeneity itself. It is therefore quite clear that position-dependency of the group cross sections should be considered in the 2-step procedure for a more accurate PWR core analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

EBR-II Static Neutronic Calculations by PHISICS / MCNP6 codes

The International Atomic Energy Agency (IAEA) launched a Coordinated Research Project (CRP) on the Shutdown Heat Removal Tests (SHRT) performed in the 80's at the Experimental fast Breeder Reactor EBR-II, USA. The scope of the CRP is to improve and validate the simulation tools for the study and the design of the liquid metal cooled fast reactors. Moreover, training of the next generation of fast reactor analysts is being also considered the other scope of the CRP. In this framework, a static neutronic model was developed, using state-of-the art neutron transport codes like SCALE/PHISICS (deterministic solution) and MCNP6 (stochastic solution). Preliminary criticality calculations have been performed using both models (deterministic and stochastic). Results show the good agreement between different models. It should also be noted that deterministic solution could still be further improved by implementing super-homogenization factors and by increasing the transport calculation order.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Comparison of the FW-CADIS and MR-CADIS Variance Reduction Methods

The Consistent Adjoint Driven Importance Sampling (CADIS) method, which uses an adjoint deterministic estimate for making variance reduction parameters that optimize Monte Carlo (MC) shielding calculations, was introduced in 1998. This method was designed to optimize a single tally, typically an integrated response such as a detector count rate or a dose rate. Because both deterministic and stochastic methods are employed, CADIS is referred to as a hybrid method. In 2007, the forward-weighted CADIS (FW-CADIS) method was introduced to optimize several tallies at once, giving more uniform relative uncertainties among the tallies for both high-flux and low-flux areas of the problem. This method can be used for several small detectors, for semi-global problems (mesh tally over a large fraction of the problem), or for global problems (mesh tally over the entire problem). Optimization can be done for an integrated response such as a dose rate D(r) as a function of position, the total flux {phi}(r) as a function of position, or the space/energy flux {phi}(r,E). Both the CADIS and FW-CADIS methods have been in the MAVRIC sequence of SCALE for three versions (6.0 in 2009, 6.1 in 2011, and 6.2 in 2016). The ADVANTG package released in 2015 uses the CADIS and FW-CADIS methods to accelerate MCNP neutron/photon shielding calculations. In 2014, the multi-response CADIS (MR-CADIS) method was presented as a way to obtain even more uniform relative uncertainties for problems with multiple tallies than those obtained with FW-CADIS. The method does not use a forward deterministic calculation but instead uses one adjoint calculation for every tally to be optimized in the Monte Carlo calculation. This work compares the above methods.

42 ENGINEERING↗

Implementation and Validation of Reference Sensitivity Profile Calculations in TRIPOLI 4 - 14247

In current neutron transport calculations, particularly on continuous-energy Monte-Carlo codes, the major source of uncertainty is coming from nuclear data. In order to improve nuclear data evaluations, integral experiments are used. Feedback is then precisely obtained by re-estimating cross-section levels and correlations according to their sensitivity. k{sub eff}-sensitivity coefficients are traditionally obtained using deterministic codes since hands-on methods were not yet implemented in continuous energy Monte Carlo codes. In this paper, an Iterated-Fission-Probability-based method newly developed in the continuous energy Monte Carlo code TRIPOLI-4 will be used to calculate sensitivity coefficient breakdown on isotopes and reactions and results will be compared to those of ERANOS deterministic code. Objective is twofold: to show the reliability of the new TRIPOLI-4 method and to identify where and why profiles are slightly different. Sensitivity profiles, forward and adjoint fluxes are thus meticulously compared using a fine 1968-group energy structure. An OECD international benchmark (Uncertainty Analysis for Criticality Safety Assessment - UACSA - Phase III) is taken as working example to enforce validation and to compare to independent results (e.g. MCNP-6). This benchmark is constituted of two bare spheres filled with homogeneous mixture of uranium and polyethylene (respectively 2% and 50% U-235 enrichments) and gives an overview of issues encountered in thermal reactor cores. Several energy-group structures (1968G, 172G, 33G) are also tested in ERANOS to identify multi-group and implicit effects. Lessons learnt are that continuous Monte Carlo methods seem to be mature and compare well with deterministic codes. Only damper the U-238 elastic - and on certain conditions the U-238 capture - sensitivity coefficient for which difference can reach 30% between TRIPOLI-4 and ERANOS. Indeed concerning U-238, implicit effects are emphasized and multi-group theory remains limited by hypothesis; even a 1968-group structure is not refined enough to properly calculate U-238 elastic scattering k{sub eff}-sensitivity. This reference 3D continuous-energy sensitivity capability in TRIPOLI4-IFP will be extensively used in the uncertainty automated tool RIB of the French Criticality-Safety package CRISTAL. (authors)

42 ENGINEERING↗

Incorporating realistic conservatism into criticality alarm placement and emergency response: a regulatory perspective - 14425

Review of historical process plant criticality accidents shows that their consequences vary widely and are extremely difficult to predict accurately. Because of the wide variability and inherent uncertainty in estimating consequences, criticality detector placement analysis and consequence estimation for emergency response typically utilize very conservative assumptions with regard to the total fissions, time-dependent fission rate, and intervening shielding. For detector placement, a 'minimum accident of concern' is typically defined, whereas emergency consequence evaluations typically define a maximum fission source term. These require different mindsets because what is conservative for one type of analysis may be non-conservative for the other. While some conservatism is commonly practiced in both cases, excessive conservatism can result in a net increase in risk, such as resulting from inadvertent alarm activations, that can be costly to operations and counterproductive to safety. Conversely, reliance on very detailed calculational models to reduce conservatism could increase the number and complexity of controls, such as shielding for dose mitigation. Deterministic regulatory criteria have also been applied for determining where criticality detectors are to be deployed and what areas must be evacuated in the event of an accident. This paper explores several issues associated with criticality alarm placement and emergency response from a regulatory perspective. This paper will discuss insights gained from historical fuel facility accidents and significant events regarding the inherent risk of criticality, the most likely locations for criticality to occur, and the most probable characteristics of a critical excursion. This will include discussion of the previously disputed relevance of the United States Nuclear Regulatory Commission's (USNRC's) longstanding model for a critical excursion. While USNRC regulations deterministically require dual alarm coverage in all areas having greater than a 'critical mass' of fissionable material, a process is provided to grant exemptions when a risk-informed argument is technically justified. Similar risk-informed considerations have been used to justify compensatory measures allowing limited operation during periods of alarm system outage due to component failure or maintenance. This paper will explore the application of risk information and performance history to address requirements for criticality alarm system coverage, granting exemptions from deterministic criteria, allowing compensatory measures during system outages, determining evacuation areas and routes, and assessing the significance associated with alarm system issues. Finally, the impact of shielding on the USNRC's regulatory risk basis and approved subcritical margin will be reviewed. (authors)

42 ENGINEERING↗

A Methodology for defining intermediate level waste package characteristics that ensure criticality safety in a geological disposal facility - 14599

Radioactive Waste Management (RWM) is responsible for implementing geological disposal of the UK's higher activity radioactive wastes. These wastes include fissile nuclides and, therefore, demonstration of the criticality safety of the wastes under disposal conditions forms an important component of RWM's disposal system safety case. The basis for ensuring criticality safety of intermediate level waste (ILW) packages is to specify controls on the fissile material content of individual packages through the application of deterministic criticality safety assessments (CSAs). The approach is defined and controlled through RWM's Disposability Assessment process, which aims to ensure that: criticality is prevented during waste package transport to a GDF; the risk of criticality is tolerable and as low as reasonably practicable (ALARP) during GDF operations; and the likelihood and consequences of criticality arising from the evolution of waste packages after GDF closure are low. The most restrictive phase determines the package fissile material limit and, until now, in many cases the GDF post-closure phase has provided the most restrictive limit. Currently, the GDF post-closure limits are derived deterministically using highly stylised and conservative assumptions about fissile material accumulations occurring on very long timescales into the future. RWM, waste producers and regulators do not currently consider this approach to be yielding an appropriate balance of risk. Therefore, RWM is modifying the process for setting fissile material limits on ILW packages by making use of the outcomes of its recent research projects to establish the likelihood of post-closure criticality in a GDF and the consequences of hypothetical post-closure criticality events (presented in companion papers). This paper presents the revised methodology, which is based on the requirement that criticality in a GDF can be demonstrated to be of low and tolerable consequences with regard the performance of a GDF. A probabilistic 'low-likelihood package envelope' is proposed that establishes the packaging and disposal facility conditions under which post-closure criticality is considered to be of low likelihood and low consequence. The bounding conditions for establishing the extent of the package envelope have been evaluated through application of the methodologies and models developed in RWM's Likelihood of Criticality research project. It is expected that the vast majority of ILW packages would meet the package envelope requirements. The methodology also allows the use of other options, such as deterministic CSAs, for waste packages that fall outside the low-likelihood package envelope. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Stochastic level-set method for shape optimisation

We present a new method for stochastic shape optimisation of engineering structures. The method generalises an existing deterministic scheme, in which the structure is represented and evolved by a level-set method coupled with mathematical programming. The stochastic element of the algorithm is built on the methods of statistical mechanics and is designed so that the system explores a Boltzmann–Gibbs distribution of structures. In non-convex optimisation problems, the deterministic algorithm can get trapped in local optima: the stochastic generalisation enables sampling of multiple local optima, which aids the search for the globally-optimal structure. The method is demonstrated for several simple geometrical problems, and a proof-of-principle calculation is shown for a simple engineering structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimal sizing study of hybrid wind/PV/diesel power generation unit

In this paper, a methodology of sizing optimization of a stand-alone hybrid wind/PV/diesel energy system is presented. This approach makes use of a deterministic algorithm to suggest, among a list of commercially available system devices, the optimal number and type of units ensuring that the total cost of the system is minimized while guaranteeing the availability of the energy. The collection of 6 months of data of wind speed, solar radiation and ambient temperature recorded for every hour of the day were used. The mathematical modeling of the main elements of the hybrid wind/PV/diesel system is exposed showing the more relevant sizing variables. A deterministic algorithm is used to minimize the total cost of the system while guaranteeing the satisfaction of the load demand. A comparison between the total cost of the hybrid wind/PV/diesel energy system with batteries and the hybrid wind/PV/diesel energy system without batteries is presented. The reached results demonstrate the practical utility of the used sizing methodology and show the influence of the battery storage on the total cost of the hybrid system. (author)

14 SOLAR ENERGY↗

Approximating Nash Equilibrium in Day-ahead Electricity Market Bidding with Multi-agent Deep Reinforcement Learning

In this paper, a day-ahead electricity market bidding problem with multiple strategic generation company (GEN-CO) bidders is studied. The problem is formulated as a Markov game model, where GENCO bidders interact with each other todevelop their optimal day-ahead bidding strategies. Considering unobservable information in the problem, a model-free and data-driven approach, known as multi-agent deep deterministic policy gradient (MADDPG), is applied for approximating the Nash equilibrium (NE) in the above Markov game. The MADDPG algorithm has the advantage of generalization due to the automatic feature extraction ability of the deep neural networks. The algorithm is tested on an IEEE 30-bus system with three competitive GENCO bidders in both an uncongested caseand a congested case. Comparisons with a truthful bidding strategy and state-of-the-art deep reinforcement learning methods including deep Q network and deep deterministic policy gradient (DDPG) demonstrate that the applied MADDPG algorithm can find a superior bidding strategy for all the market participants with increased profit gains. In addition, the comparison with a conventional model-based method shows that the MADDPG algorithm has higher computational efficiency, which is feasible for real-world applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quasi-Stochastic Approximation and Off-Policy Reinforcement Learning: Preprint

The Robbins-Monro stochastic approximation algorithm is a foundation of many algorithmic frameworks for reinforcement learning and often an efficient approach to solving (or approximating the solution to) complex optimal control problems. In many cases, however, practitioners are unable to apply these techniques because of an inherent high variance. This paper aims to provide a general foundation for 'quasi-stochastic approximation,' in which all of the processes under consideration are deterministic, much like quasi-Monte Carlo for variance reduction in simulation. The variance reduction can be substantial, subject to tuning of pertinent parameters in the algorithm. This paper introduces a new coupling argument to establish the optimal rate of convergence provided the gain is sufficiently large. These results are established for linear models and tested also in nonideal settings. A major application of these general results is a new class of reinforcement learning algorithms for deterministic state space models. In this setting, the main contribution is a class of algorithms for approximating the value function for a given policy, using a different policy designed to introduce exploration.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Chernobyl, the HBO miniseries: Fact and fiction (Part III)

With episode four of the series, we moved even further from reality than in prior episodes. First, we need background. Exposure to ionizing radiations causes two types of medical effects: deterministic and stochastic. Deterministic effects are predictable, dose-dependent, and occur in everyone exposed to the same dose. For example, everyone exposed to an acute whole-body dose of 5 gray (5,000 millisieverts) will have a marked immediate decrease in blood granulocytes.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-Performance Mesostructured Organic Hybrid Pseudocapacitor Electrodes

The electrodes of a hybrid electrochemical capacitor which utilize the quinone (Q)-hydroquinone (QH 2 ) couple, a prototypical organic redox system known to provide fast and reversible proton-coupled electron-transfer reactions, are deterministically mesostructured via a colloidal templating strategy to provide good ion and electron transport pathways, enabling a high rate performance. Specifically, a conducting polymer, polypyrrole (PPy), is functionalized with a pseudocapacitive material, a Q/QH 2 -containing catechol derivative, by noncovalent interactions. The mesostructure of this hybrid material is formed into an ordered 3D porous structure by a polystyrene colloidal crystal template-assisted electrosynthesis. The catechol derivative is sufficiently bound to the PPy through noncovalent interactions to provide a volumetric capacitance as high as ≈130 F cm –3 and a capacitance retention of ≈75% over 10 000 charging/discharging cycles. Furthermore, when compared with a randomly structured electrode, the deterministically structured electrode exhibits an improved rate performance due to the mesostructure facilitated electron and ion transport.

36 MATERIALS SCIENCE↗

Unraveling the Main Chain and Side Chain Effects on Thin Film Morphology and Charge Transport in Quinoidal Conjugated Polymers

Three series of low-bandgap polymers based on a novel quinoidal para -azaquinodimethane ( p -AQM) unit are devised and synthesized, enabling an in-depth study of the impact of structural factors such as polymer main chain, branching point of the side chain, and the length of the branch chains on the thin film morphologies and charge transport properties. Morphological studies reveal that the polymers composed of larger repeating units exhibit a stronger tendency to form edge-on lamella. In contrast, altering the side chain structures of polymers with the same main chain configuration indicates that the branching point position has a more deterministic impact than the branch chain length on the interchain interactions and the crystallite orientation. These results demonstrate a compound odd-even effect of the branching point on the chain packing and morphology, which correlates well with the corresponding field effect transistor performances. The polymer with the branching point at the fourth carbon displays the highest charge carrier mobility over 1.0 cm 2 V-Three series of low-bandgap polymers based on a novel quinoidal para -azaquinodimethane ( p -AQM) unit are devised and synthesized, enabling an in-depth study of the impact of structural factors such as polymer main chain, branching point of the side chain, and the length of the branch chains on the thin film morphologies and charge transport properties. Morphological studies reveal that the polymers composed of larger repeating units exhibit a stronger tendency to form edge-on lamella. On the other hand, altering the side chain structures of polymers with the same main chain configuration indicates that the branching point position has a more deterministic impact than the branch chain length on the interchain interactions and the crystallite orientation. These results reflect a compound odd-even effect of the branching point on the chain packing and morphology, which correlates well with the corresponding field effect transistor performances. The polymer with the branching point at the fourth carbon displays the highest charge carrier mobility over 1.0 cm 2 V -1 s -1 , concurrent with a bimodal texture. This study provides a comprehensive description of the correlations between polymer structures, thin film morphology, and device performances, providing a clear path to desirable bimodal thin film texture for charge transport.

36 MATERIALS SCIENCE↗

Mutual Ferromagnetic-Ferroelectric Coupling in Multiferroic Copper-Doped ZnO

A mutual ferromagnetic and ferroelectric coupling (multiferroic behavior) in Cu-doped ZnO is demonstrated via deterministic control of Cu doping and defect engineering. The coexistence of multivalence Cu ions and oxygen vacancies is important to multiferroic behaviors in ZnO:Cu. The samples show clear ferroelectric and ferromagnetic domain patterns. These domain structures may be written reversibly via electric and magnetic bias.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Conversion of Classical Light Emission from a Nanoparticle‐Strained WSe 2 Monolayer into Quantum Light Emission via Electron Beam Irradiation

Abstract Solid‐state single photon emitters (SPEs) within atomically thin transition metal dichalcogenides (TMDs) have recently attracted interest as scalable quantum light sources for quantum photonic technologies. Among TMDs, WSe 2 monolayers (MLs) are promising for the deterministic fabrication and engineering of SPEs using local strain fields. The ability to reliably produce isolatable SPEs in WSe 2 is currently impeded by the presence of numerous spectrally overlapping states that occur at strained locations. Here nanoparticle (NP) arrays with precisely defined positions and sizes are employed to deterministically create strain fields in WSe 2 MLs, thus enabling the systematic investigation and control of SPE formation. Using this platform, electron beam irradiation at NP‐strained locations transforms spectrally overlapped sub‐bandgap emission states into isolatable, anti‐bunched quantum emitters. The dependence of the emission spectra of WSe 2 MLs as a function of strain magnitude and exposure time to electron beam irradiation is quantified and provides insight into the mechanism for SPE production. Excitons selectively funnel through strongly coupled sub‐bandgap states introduced by electron beam irradiation, which suppresses spectrally overlapping emission pathways and leads to measurable anti‐bunched behavior. The findings provide a strategy to generate isolatable SPEs in 2D materials with a well‐defined energy range.

2D material↗

Hierarchical Composites Patterned via 3D Printed Cellular Fluidics

Additive manufacturing of freeform structures containing multiple materials with deterministic spatial arrangement and interactions remains a challenge for most 3D printing processes, due to complex fabrication tool requirements and limitations in printability of some material classes. Here, in this paper, a versatile method is reported to produce architected composites using the concept of cellular fluidics, in which lattices of unit cells are used as templating scaffolds to guide flowable infill materials in a programmed spatial pattern, upon which they are cured in place to produce a deterministically ordered multimaterial solid. The lattice design relies on the unit cell size, type, strut diameter, surface wetting, and distribution of cellular structures to control liquid flow and retention. Individual unit cells are tuned to achieve reliable infilling and combined into higher-order architectures to achieve multiscale composite materials with disparate mechanical properties, including those considered non-printable. Lattice design considerations for leveraging capillary phenomena and demonstrate several methods of patterning polymers in 3D-printed cellular fluidic structures are presented. The concept of tuning the compressive response of an architected composite using a flexible-elastomer as the lattice and a stiff-epoxy as the infill material is illustrated.

36 MATERIALS SCIENCE↗

Quantitative Electron Beam‐Single Atom Interactions Enabled by Sub‐20‐pm Precision Targeting

The ability to probe and control matter at the picometer scale is essential for advancing quantum and energy technologies. Scanning transmission electron microscopy offers powerful capabilities for materials analysis and modification, but sample damage, drift, and scan distortions hinder single atom analysis and deterministic manipulation. Materials analysis and modification via electron–solid interactions can be transformed by precise delivery of electrons to a specified atomic location, maintaining the beam position despite drift, and minimizing collateral dose. Here a fast, low-dose, sub-20-pm precision electron beam positioning technique is developed, “atomic lock-on,” (ALO), which offers the ability to position the beam on a specific atomic column without previously irradiating that column. This technique is used to lock onto a single selected atomic location to repeatedly measure its weak electron energy loss signal despite sample drift. Moreover, electron beam-matter interactions in single atomic events are measured with μ s time resolution. This enables observation of single-atom dynamics, such as atomic bistability, revealing partially bonded atomic configurations and recapture phenomena. This opens prospects for using electron microscopy for high-precision measurements and deterministic control of matter for quantum technologies.

2D materials↗

Nth-plant supply: corn stover supplies and costs in a fleet of biorefineries

Feedstock cost and cost variability is expected to increase with the number of biorefineries. To quantify this effect, this spatial-economic analysis simulates feedstock cost and cost variability of an industry based on corn stover as a function of the number of biorefineries. Results are reported for nine scenarios (a base case and sensitivity analysis of four variables – harvest efficiency, sustainability constraints, opportunity cost, and corn grain yield) under deterministic and stochastic simulations, assuming biorefineries using 658 000 Mg (725 000 tons) year- 1 of corn stover in 2019. The resulting supply curves are highly elastic (i.e. little change in cost) for the first 50 of the 121 biorefineries, with price increases in subsequent biorefineries depending on scenario. In the base-case deterministic scenario, weighted-average stover costs are $66 Mg -1 ($60 ton- 1 ), $69 Mg -1 ($62 ton -1 ), and $156 Mg -1 ($142 ton -1 ), at the first, 60th, and 121st biorefineries, respectively. The stochastic simulations, subject to observed 30-year corn yield variability, follow a similar pattern, with price distributions that vary by scenario. The base-case stochastic simulations illustrate minimal cost variability for the first 60 biorefineries, but rapid increases in cost variability in the second half of potential biorefineries, with similar patterns observed in the other scenarios. Of the four variables explored, price was most sensitive to harvest efficiency, followed by sustainability constraints, corn yield, and opportunity cost. Results suggest that, under conventional logistics, about half of the US corn stover resource is reliably available with minimum cost increase and variability. Interactive visualization is available at https://doi.org/10.11578/1828779.

09 BIOMASS FUELS↗