Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASPEN”

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 55 records · Page 3

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal Integration of Advanced Nuclear Reactors with a Reference Refinery, Methanol Synthesis, and a Wood Pulp Plant (Rev.1)

The present report is intended to provide process flow diagrams (PFDs) and energy and mass balance data sheets for a U.S. industrial sector subset with which nuclear heat and power could be integrated—a subset that includes the oil refining, methanol and pulp and paper industries. Coupling options for integrating nuclear energy into these industries are quantitatively outlined for reference systems, and future work will extend this analysis in greater detail. Opportunities for integrating small modular nuclear reactors (SMNRs) were investigated for each of the industrial process configurations. Aspen HYSYS and Cycle-Tempo models for a high-temperature gas-cooled reactor were developed to evaluate the proposed integration. This introductory evaluation provides a general description and assessment of the operating principles, reactor coolant core outlet temperature, and reactor size to be integrated with industry. The industrial processes of oil refining and the production of methanol, pulp and paper were simulated by using Aspen HYSYS, Aspen Plus, and the PRELIM (Petroleum Refinery Life Cycle Inventory Model) tool to develop process details. Cycle-Tempo models then extend the process modeling results to obtain net energy demands (e.g., heat, steam, and electricity) when accounting for process steam and waste heat recovery. This information is intended to foster the analysis of integrating an SMNR to decarbonize industrial facilities. The SMNR would provide reliable, competitive, and sustainable clean energy while reducing carbon emissions and other environmental impacts, such as water withdrawals, consumption, and contamination. The refining industry, exhibited in Figure ES1, is a leading consumer of fossil -fuel-based heat, power, and hydrogen in the U.S. industrial sector, generating over 164 million metric tons (MMT) of CO 2 emissions in 2023. The overall mass and energy pertaining to a generalized complex refinery in the United States is reflected in Figure ES1, along with energy metrics regarding integration with a nuclear power plant (NPP). Data sheets were developed to indicate the energy requirements for the overall refinery and each refinery process. The data sheet for the overall refinery is shown in Table ES2.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Functional Imaging of Microbial Interactions With Tree Roots Using a Microfluidics Setup

Coupling microfluidics with microscopy has emerged as a powerful approach to study at cellular resolution the dynamics in plant physiology and root-microbe interactions (RMIs). Most devices have been designed to study the model plant Arabidopsis thaliana at higher throughput than conventional methods. However, there is a need for microfluidic devices which enable in vivo studies of root development and RMIs in woody plants. Here, we developed the RMI-chip, a simple microfluidic setup in which Populus tremuloides (aspen tree) seedlings can grow for over a month, allowing continuous microscopic observation of interactions between live roots and rhizobacteria. We find that the colonization of growing aspen roots by Pseudomonas fluorescens in the RMI-chip involves dynamic biofilm formation and dispersal, in keeping with previous observations in a different experimental set-up. Also, we find that whole-cell biosensors based on the rhizobacterium Bacillus subtilis can be used to monitor compositional changes in the rhizosphere but that the application of these biosensors is limited by their efficiency at colonizing aspen roots and persisting. These results indicate that functional imaging of dynamic root-bacteria interactions in the RMI-chip requires careful matching between the host plant and the bacterial root colonizer.

59 BASIC BIOLOGICAL SCIENCES↗

Development of an innovative process for post-combustion CO 2 capture to produce high-value NaHCO 3 nanomaterials

An innovative post-combustion process using aqueous sodium glycinates solutions (SGS) for CO 2 capture from a split flue gas stream emitted from the 600 MWe coal power plant, used in the Wolverine Clean Energy Venture (WCEV) project, was developed in Aspen Plus v.10. The flue gas flow rate used in the process was 12.43 kg/s (at 353.15 K and 101.33 kPa) and contained 0.0023 and 13.33 mol% of SO 2 and CO 2 , respectively. The overall process includes 5 main units designed to remove all SO 2 and capture more than 90 mol% of CO 2 in the flue gas stream, while producing high-value, salable sodium bicarbonate (NaHCO 3 ) nanomaterials to offset the total process costs. The hydraulics, mass transfer characteristics, and process performance obtained using Aspen Plus were discussed. Also, the capital expenditure (CAPEX), operating expenditure (OPEX), and Levelized cost of CO 2 capture (LCOC) were calculated to assess the feasibility of this process. Further, the hydraulics in the SO 2 washing and CO 2 capture units showed a pressure drop of 12 and 1 kPa, respectively, and the behaviors of the liquid holdup and normalized packing specific wetted surface area were similar in both units. The gas-side mass transfer coefficients were orders of magnitude greater than the liquid-side mass transfer coefficients. The process was able to capture 2.352 kg/s of CO 2 and produce 4.486 kg/s of valuable NaHCO 3 nanomaterials. Also, the calculated CAPEX, OPEX, and LCOC of the process for a 30-year plant lifetime were ($\$ $4,450,552), (233.00 $\$ $/h) and (35.49 $\$ $/ton of CO 2 captured), respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Predictive Model for Offgas Composition in Waste Glass Melters - Spring 2021 SULI General Audience Abstract

Approximately 56 million gallons of nuclear waste are currently stored in 177 underground storage tanks at the Hanford site in Washington state. A portion of the tank waste is targeted for immobilization by converting the nuclear waste to borosilicate glass through a process called vitrification. During vitrification, emissions are produced in addition to the intended glass product. These emissions, contained in a volume above the molten glass called the plenum, are composed of water vapor, air components, aerosols, particulates, and other gaseous species that result from evaporation, reactions, and air influx. The formation of these gasses causes the surface of the glass to foam. To reduce foaming, sucrose is added to the waste stream before the conversion to glass. Sucrose also reduces the formation of hazardous gasses in the plenum including nitrous oxides, NOx, and carboxides, COx. However, if excess sucrose is present, the reactions that produce the gasses in the plenum will be incomplete, and increased levels of NOx and COx can be observed. The Waste Treatment and Immobilization Plant at the Hanford site is equipped with the technology to treat the plenum emissions. Though many tests have been run to determine the presence of particulates in the plenum emissions, there is does not exist an established method to confidently predict the composition of the gaseous species in the plenum prior to treatment without physical testing. A model was created in the process modeling software Aspen Plus to predict the composition of NOx and COx in the plenum. The model employs a constant stirred tank reactor and user-defined reaction chemistry to determine the plenum composition. To verify the model, the offgas predicted by the model were compared to available offgas data for ten different feed compositions. It is demonstrated that the developed Aspen Plus model is capable of predicting the composition of NOx and COx in the plenum.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Expanding the horizon of bio-naphtha beyond gasoline blend: property characterization and conversion opportunity assessment through technoeconomic and life-cycle analyses

Bio-naphtha, a common by-product of biorefineries, is expected to experience substantial growth in supply due to increasing demands for renewable diesel and synthetic aviation fuel (SAF). However, demand for bio-naphtha itself as a gasoline blendstock is limited because of the electrification of light-duty vehicles. This work investigated valorization opportunities for bio-naphtha from catalytic fast pyrolysis, hydrothermal liquefaction, Fischer–Tropsch synthesis, and hydrotreated esters and fatty acids pathways. These opportunities include producing polymer-grade olefin via steam cracking, SAF via steam cracking followed by olefin oligomerization, and renewable aromatics benzene, toluene, and xylene (BTX) and hydrogen via catalytic reforming. Process models were developed in Aspen Plus V14 and Aspen HYSYS V14 to calculate the mass and energy balances for each conversion step. Technoeconomic assessment and life-cycle analysis were conducted to evaluate the minimum fuel/product selling price, conversion cost, and life-cycle CO2 equivalent (CO2e) emission reduction. Technoeconomic assessment results suggest a minimum fuel/product selling price as low as $1.9/kg of olefins, $6.30/gal of gasoline-equivalent SAF, and $1.2/kg of BTX without any incentives. For all pathways, these prices are dominated by bio-naphtha feedstock costs, which account for at least 76% of the total cost. Compared with petroleum baselines, bio-naphtha-derived SAF, olefins, and BTX can achieve significant CO2e emission reductions from the use of renewable carbon resources. The results of life-cycle analysis and subsequent technoeconomic assessment, incorporating carbon credits, indicate the economic viability of using bio-naphtha for polymer-grade olefin and BTX production, with product costs comparable to market prices.

Xu, Shuang↗

Economic and Environmental Analysis to Evaluate the Potential Value of Co-Optima Diesel Bioblendstocks to Petroleum Refiners

The U.S. petroleum refining sector is undergoing a period of historic transformation, catalyzed by the decarbonization of the U.S. economy. Diesel-boiling-range bioblendstocks have gained traction, owing to their superior fuel properties and environmental performance as compared to traditional petroleum fuels. This work couples refinery linear programming models with life cycle assessment to quantify the potential economic and environmental benefits, and trade-offs, of blending diesel-boiling-range bioblendstocks at petroleum refineries. Linear programming models were developed in Aspen Process Industry Modeling Systems (PIMS) for three representative petroleum refinery configurations of differing complexity. Seven diesel-boiling-range bioblendstocks: 4-butoxyheptane, 5-ethyl-4-propylnonane, soy biodiesel, sludge hydrothermal liquefaction diesel, polyoxymethylene ethers, renewable diesel, and hexyl hexanoate, were investigated to identify key fuel properties that influence refineries' economics and to track the effect of adding bioblendstocks on refinery-wide cradle-to-gate greenhouse gases (GHG) emissions. These analyses considered blending levels from 10 to 30 vol% and fuel demand projections over the period 2040 to 2050. This analysis determines that bioblendstock sulfur content and cetane number are the primary fuel attributes with the potential to provide value to refiners. Life cycle assessment results indicate that the use of diesel-boiling-range bioblendstocks can reduce cradle-to-gate refinery GHG emissions by up to ~ 40 % relative to conventional refinery operations when considering carbon uptake in the supply chain of the bioblendstock. Refinery-wide marginal GHG abatement costs range from 120 to 3,600 USD2016/metric tons carbon dioxide equivalent avoided across the scenarios evaluated. Reducing the price of bioblendstocks is identified as a key to their adoption.

ADVANCED PROPULSION SYSTEMS,BIOMASS FUELS,ENERGY P↗

Dual-map framework for noise characterization of quantum computers

In order to understand the capabilities and limitations of quantum computers, it is necessary to develop methods that efficiently characterize and benchmark error channels present on these devices. In this paper, we present a method that faithfully reconstructs a marginal (local) approximation of the effective noise (MATEN) channel, that acts as a single layer at the end of the circuit. We first introduce a dual-map framework that allows us to analytically derive expectation values of observables with respect to noisy circuits. These findings are supported by numerical simulations of the quantum approximate optimization algorithm (QAOA) that also justify the MATEN, even in the presence of nonlocal errors that occur during a circuit. Finally, we demonstrate the performance of the method on Rigetti's Aspen-11 quantum computer for QAOA circuits up to six qubits, successfully predicting the observed measurements on a majority of the qubits.

Sud, James↗

Direct pulse-level compilation of arbitrary quantum logic gates on superconducting qutrits

Advanced simulations and calculations on quantum computers require high-fidelity implementations of quantum operations. The universal gateset approach builds complex unitaries from a small set of primitive gates, often resulting in a long gate sequence, which is typically a leading factor in the total accumulated error. Compiling a complex unitary for processors with higher-dimensional logical elements, such as qutrits, exacerbates the accumulated error per unitary, since an even longer gate sequence is required. Optimal control methods promise time- and resource-efficient compact gate sequences and, therefore, higher fidelity. These methods generate pulses that can directly implement any complex unitary on a quantum device. In this work, we demonstrate that any arbitrary qubit and qutrit gate can be realized with high fidelity, which can significantly reduce the length of a gate sequence. We generate and test pulses for a large set of randomly selected arbitrary unitaries on several quantum processing units (QPUs): the Lawrence Livermore National Laboratory Quantum Device and Integration Testbed’s (QuDIT’s) standard QPU and three of Rigetti’s QPUs: Ankaa-2, Ankaa-9Q-1, and Aspen-M-3. On the QuDIT platform’s standard QPU, the average fidelity of random qutrit gates is 97.9 ± 0.5% measured with conventional QPT and 98.8 ± 0.6% from QPT with gate folding. Rigetti’s Ankaa-2 achieves random qubit gates with an average fidelity of 98.4 ± 0.5% (conventional QPT) and 99.7 ± 0.1% (QPT with gate folding). On Ankaa-9Q-1 and Aspen-M-3, the average fidelities with conventional qubit QPT measurements were higher than 99% (see Appendix). Here we show that optimal control gates are robust to drift for at least 3 h and that the same calibration parameters can be used for all implemented gates. Our work promises that the calibration overheads for optimal control gates can be made small enough to enable efficient quantum circuits based on this technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

CHESS 2025: Leaf Area Index (LAI) for meadow, shrub, tree, and understory vegetation

This dataset contains Leaf Area Index (LAI) measurements made as part of the Colorado Headwaters Ecological Spectroscopy Study (CHESS) during June and July of 2025. Data were collected in the Upper Gunnison Basin, Colorado, across three study domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). Field observations of LAI were collected within 72 hours of airborne data collection by the National Ecological Observatory Network’s Aerial Observation Platform (NEON AOP). The NEON AOP collected waveform LiDAR (Light Detection and Ranging) and imaging spectrometer data in 426 spectral bands from the visible to shortwave infrared. LAI measurements were collected using the LICOR LAI-2200C Plant Canopy Analyzer following protocols outlined in the instrument manual (LI-COR 2019). Sampling targeted four distinct vegetation types: meadows, shrubs, trees, and aspen forest understory. We have archived data separately by site type because different field methods were used for each. At meadow sites, measurements were made at the four corners of 1m x 1m plots, with the instrument moving inward toward the center of the plot. At shrub sites, we measured the canopies of individual shrubs. At tree sites, we made measurements within a 10m x 10m subplot centered around a focal tree, with 30 observations taken on a regular grid. At aspen understory sites, we measured overstory trees following the tree protocol and understory herbaceous vegetation following the meadow protocol. All measurements included above-canopy (A) and below-canopy (B) readings, with specific protocols for scattering correction measurements in direct-sun conditions. Data were processed using the R package `rlai` (Worsham 2025). This package includes functions to calculate LAI, gap fraction, apparent clumping factor (Ω), scattering correction, and other canopy metrics. Package contents: Full file descriptions appear in ‘flmd.csv’. Files named according to the convention ‘lai_*_summary_data_cleaned.csv’ contain summary values of LAI, apparent clumping factor (Ωapp), and scattering correction factors for each site. These are the analysis-ready products that most data users will work with. Files named ‘lai_*_metadata_cleaned.csv’ contain additional site-level observations made during field collection. We have also archived intermediate and supplementary data for users who wish to check our processing approach or apply alternative methods. ‘raw_lai_2200C.zip’ contains the raw files as read from the LI-COR instrument, with no processing applied, in TXT format. The zip archive contains subdirectories by site type, which are further subdivided by sampling area. Filenames correspond to the sampling site number. ‘intermediate_results.zip’ contains detailed output from the processing routines, in JSON format. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘scattering_correction_logs.zip’ contains logfiles from the implementation of Kobayashi et al.'s (2013) scattering correction algorithm. The logfiles report values of several parameters at each iteration of the algorithm, as the model converges toward a stable solution. They are intended for users who want to verify scattering correction performance. The zip archive contains subdirectories by site type; filenames correspond to the sampling site number. ‘spot_checks.csv’ reports LAI and other values for a small number of files processed with LI-COR FV2200 software (LI-COR 2013) using the same control parameters as in our R-based approach. Additional metadata are provided in a data dictionary describing column names and definitions (dd.csv), and in a file-level metadata file (flmd.csv). All zip files can be expanded with common archive utilities. TXT, CSV, and JSON files can be ingested into R or Python computing environments or read in common text editor utilities. Geospatial information: Geospatial data for mapping measurement site locations are in the files CHESS_polygons_lai_UTM.geojson, CHESS_polygons_shrub_UTM.geojson, and CHESS_polygons_meadow_UTM.geojson in the companion geospatial package for the 2025 CHESS campaign, ‘CHESS 2025: Location data for field observations and sampling’ (Henderson et al., 2026). CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. * Todorov and Worsham are co–first authors.

2018 NEON and 2025 CHESS Campaigns↗

AmeriFlux FLUXNET-1F US-UMd UMBS Disturbance

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-UMd UMBS Disturbance. This is the FLUXNET version of the carbon flux data for the site US-UMd UMBS Disturbance produced by applying the standard ONEFlux (1F) software. Site Description - The UMBS Disturbance site is an artificial disturbance site that has recently been created as part of the Forest Accelerate Succession ExperimenT (FASET). In Spring 2008, every aspen and birch tree (>6,700, ~35% canopy LAI), the dominant early successional trees, were girdled over 39 ha of the FASET treatment plot to stimulate a disturbance that will move the forest into a later successional stage, dominated by maples, oaks, and white pine. This treatment caused aspen and birch mortality within 2 - 3 years. As a result of the changed canopy structure, there is a divergence in net ecosystem exchange between the control plot (enhanced carbon uptake) and the treatment plot (reduced carbon uptake).

Classen, Aimee↗

AmeriFlux FLUXNET-1F CA-SF3 Saskatchewan - Western Boreal, forest burned in 1998

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site CA-SF3 Saskatchewan - Western Boreal, forest burned in 1998. This is the FLUXNET version of the carbon flux data for the site CA-SF3 Saskatchewan - Western Boreal, forest burned in 1998 produced by applying the standard ONEFlux (1F) software. Site Description - The 1998 burn site (F98) was in the east part of Prince Albert National Park, Saskatchewan, in the Waskesiu Fire, ignited by lightning that burned about 1700 ha in July 1998. The pre-fire forest consisted of jack pine and black spruce stands, with some intermixed aspen. The fire was severe, consuming much of the top layer of organic soil and killing all trees. In 2001, much of the regenerating vegetation consisted of aspen saplings about 1 m tall and shorter jack pine and black spruce seedlings. An overstory of dead, leafless jack pine trees dominated at a height of 18 m. Sparse grass and herbs, such as fireweed (Epilobium angustifolium L.) covered the ground. There were a large number of fallen dead trees, mostly perched above the ground and not decomposing quickly.

Amiro, Brian [University of Manitoba; Canadian For↗

Absorber Column CFD model validation against PNNL’s device-scale absorber column on the LCFS unit

Absorber column has been widely used for CO 2 capture in the coal-fired power plants. High-fidelity CFD models play an important role in absorber column design and solvents optimization, which help enhance the CO 2 capture efficiency and reduce the operation cost. This report provides a comprehensive description of the development of CFD absorber models at two different levels, the design and implementation of PNNL’s device-scale absorber column experiment, and the methodology to combine the CFD results, experiment data, and Aspen model for a better understanding of the interface area in packed column. A composite model is firstly proposed in the Discrete Element Method (DEM) packing process, which can model the complex geometry of various packing elements. This generates a realistic packing pattern and accurate packing porosity ε and specific area a p compared to the actual values for absorber column used in LCFS. Two level of CFD absorber models were developed, namely the full-size column model (FCM) to simulate the entire packed column with a focus on the wall/entrance effects, and the representative column model (RCM) to simulate a section of column with a focus on the sensitivity study of interface area. A Design of Experiment (DoE) plan was developed to guide the collection of 100 run CFD data and 12 experiment runs. The 100 CFD runs were carried out in the RCM with Pro-Pak packing and cover a wide operation range and solvent properties. Impact of influential factors on the interface area in packed column were investigated in details. The CFD interface area was then combined with the experimental data and Aspen model to infer some information of the effective contact angle in the column. The accuracy and uncertainties in the interface area and contact angle can be quantified.

01 COAL, LIGNITE, AND PEAT↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Quality Guidelines for Energy System Studies: Process Modeling Design Parameters

The National Energy Technology Laboratory (NETL) conducts systems analysis studies that require a large number of inputs, from ambient conditions to parameters for Aspen Plus ® (Aspen) process blocks. The sheer number of assumptions required makes it impractical to document all of them in each issued report. The purpose of the Quality Guidelines for Energy System Studies (QGESS) is to document the assumptions most commonly used in system analysis studies and the basis for those assumptions. In order to develop the systems analysis models presented in various NETL reports, significant vendor data have been obtained, and these data enhance the model outputs. Much of the vendor data obtained are considered proprietary and not suitable for public release or attribution to a specific vendor. As such, several sub-systems common in NETL reports and their process parameter data are not reported in this document to protect proprietary vendor information. The values and ranges of values presented in this report represent assumptions that have been made in previous studies.

97 MATHEMATICS AND COMPUTING↗

Process Design and Techno-Economic Analysis of the Modular Staged Pressurized Oxy-Combustion (SPOC) Power Plant for Biomass

This work describes the process design and techno-economic analysis (TEA) of the modular SPOC power plant for biomass firing and coal-biomass co-firing. Two Rankine cycles were considered: a supercritical steam cycle (242 bar, 593°C, 593°C) with 550 MWe net output and a subcritical cycle (166 bar, 566°C, 566°C) with 200 MWe net output. For both cases, 95% carbon capture was modeled, and hybrid poplar biomass was chosen to generate carbon-negative power. In addition, the supercritical 500 MWe case included a 25% biomass co-firing (carbon neutral) case. For both cycles, a 100% Powder River Basin coal firing case was used for comparison purposes. In the SPOC process, oxygen is produced via a cryogenic air separation unit (ASU) and the heat generated from the compression of air is integrated into the steam cycle and utilized for boiler feed water pre-heating. Unique to the SPOC process, the boilers are pressurized and arranged in a series-parallel configuration, with minimized flue gas recirculation. The flue gas is cooled and scrubbed in the direct-contact cooler (DCC) column, and the moisture in the flue gas is condensed, leaving the bottom of the DCC at a sufficiently high temperature such that it can be used for boiler feed water pre-heating, improving plant thermal efficiency. Following drying and purification, CO2 in the flue gas is at the purity required for storage or utilization. The performance data were obtained from process modelling via Aspen Plus®. The stream data from Aspen Plus® were used as an input for the AACE Class 5 cost study. Ultimately, the capital costs, Levelized Cost of Electricity (LCOE), and cost of CO2 captured and avoided were obtained. The HHV efficiency of the carbon negative 550 MWe supercritical SPOC case (34.8%) was clearly above those reported by NETL for the BECCS baseline cases of supercritical pulverized coal with capture (B12B, 31.5%) and the 49% biomass co-firing case with capture (PA3, 29.2%). The HHV efficiency of the carbon-negative subcritical plant is also higher than the subcritical baseline PC plant with capture (case B11B.95) presented by NETL (32% vs 29.7%). The LCOE for the SPOC 100% biomass case was similar to the LCOE for the BECCS 49% biomass with carbon capture case ($147/MWh), and the SPOC carbon neutral case LCOE was lower ($110/MWh) than the cost for the NETL baseline SC coal firing case with 90% carbon capture ($114/MWh).

Magalhaes, Duarte↗