Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input estimation”

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 1,171 records · Page 65

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

Improved Bayesian regularization of inverse problems in vibrations and acoustics using noise-only measurements

Here, this paper studies Tikhonov regularization (ridge regression) parameter selection for problems in vibrations and acoustics. The selection method is based on a popular Bayesian method, but it incorporates measurements of sensor noise. The regularization parameter is closely related to the ratio of system input energy to noise energy, so noise measurements inform the inference procedure and improve parameter identification. In cases where standard Bayesian regularization identifies zero as the optimal regularization parameter, noise measurements guarantee a unique nonzero optimum. Sufficient theoretical criteria are developed for this guarantee. The method is verified in even-determined and under-determined configurations in an acoustic source localization simulation and a vibration load identification experiment. It is shown to yield significant improvements over existing empirical Bayesian regularization. Improvements are larger in the even-determined case and smaller in the under-determined case, wherein the inverse solution is less sensitive to the regularization parameter.

42 ENGINEERING↗

Cross-national analysis of food security drivers: comparing results based on the Food Insecurity Experience Scale and Global Food Security Index

Abstract The second UN Sustainable Development Goal establishes food security as a priority for governments, multilateral organizations, and NGOs. These institutions track national-level food security performance with an array of metrics and weigh intervention options considering the leverage of many possible drivers. We studied the relationships between several candidate drivers and two response variables based on prominent measures of national food security: the 2019 Global Food Security Index (GFSI) and the Food Insecurity Experience Scale’s (FIES) estimate of the percentage of a nation’s population experiencing food security or mild food insecurity (FI ). We compared the contributions of explanatory variables in regressions predicting both response variables, and we further tested the stability of our results to changes in explanatory variable selection and in the countries included in regression model training and testing. At the cross-national level, the quantity and quality of a nation’s agricultural land were not predictive of either food security metric. We found mixed evidence that per-capita cereal production, per-hectare cereal yield, an aggregate governance metric, logistics performance, and extent of paid employment work were predictive of national food security. Household spending as measured by per-capita final consumption expenditure (HFCE) was consistently the strongest driver among those studied, alone explaining a median of 92% and 70% of variation (based on out-of-sample R 2 ) in GFSI and FI , respectively. The relative strength of HFCE as a predictor was observed for both response variables and was independent of the countries used for model training, the transformations applied to the explanatory variables prior to model training, and the variable selection technique used to specify multivariate regressions. The results of this cross-national analysis reinforce previous research supportive of a causal mechanism where, in the absence of exceptional local factors, an increase in income drives increase in food security. However, the strength of this effect varies depending on the countries included in regression model fitting. We demonstrate that using multiple response metrics, repeated random sampling of input data, and iterative variable selection facilitates a convergence of evidence approach to analyzing food security drivers.

42 ENGINEERING↗

Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.21 (20200601)

The continental flood inundation mapping (CFIM) framework is a high-performance computing (HPC)-based computational framework for the Height Above Nearest Drainage (HAND)-based inundation mapping methodology. Using the 10m Digital Elevation Model (DEM) data produced by U.S. Geological Survey (USGS) 3DEP (the 3-D Elevation Program) and the NHDPlus hydrography dataset produced by USGS and the U.S. Environmental Protection Agency (EPA), a hydrological terrain raster called HAND is computed for HUC6 units in the conterminous U.S. (CONUS). The value of each raster cell in HAND is an approximation of the relative elevation between the cell and its nearest water stream. Derived from HAND, a hydraulic property table is established to calculate river geometry properties for each of the 2.7 million river reaches covered by NHDPlus (5.5 million kilometers in total length). This table is a lookup table for water depth given an input stream flow value. Such lookup is available between water depth 0m and 25m at 1-foot interval. The flood inundation map is then computed by using HAND and this lookup table based on the near real-time water forecast from the National Water Model (NWM) at the National Oceanic and Atmospheric Administration (NOAA). HAND and the Hydraulic Property Table version v0.21 is an upgrade to v0.2. Version 0.21 applies NHD HR etching to DEMs for all HUC6 units before running the HAND workflow. This new version also fixed the issue of void areas along HUC boundary. This dataset comprises 331 HUC6 units for CONUS (excluding the five great lakes units), each is a downloadable zip file. Version 0.1 was computed at the National Center for Supercomputing Applications (NCSA) at the University of Illinois at Urban-Champaign in 2016 and is currently hosted at the Texas Advanced Computing Center (TACC). Please cite this data DOI and the following publications and data DOIs when you use HAND and the hydraulic property table: Liu, Yan Y. Height Above Nearest Drainage (HAND) and Hydraulic Property Table for CONUS - Version 0.2. (20200301). Oak Ridge Leadership Computing Facility. (2020). DOI: 10.13139/ORNLNCCS/1608331 Liu, Yan Y., David R. Maidment, David G. Tarboton, Xing Zheng, and Shaowen Wang. 'A CyberGIS integration and computation framework for high-resolution continental-scale flood inundation mapping.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 770-784. DOI: 10.1111/1752-1688.12660 Zheng, Xing, David G. Tarboton, David R. Maidment, Yan Y. Liu, and Paola Passalacqua. 'River channel geometry and rating curve estimation using height above the nearest drainage.' JAWRA Journal of the American Water Resources Association 54, no. 4 (2018): 785-806. DOI: 10.1111/1752-1688.12661 Liu, Y. (2018). Height Above Nearest Drainage (HAND) for CONUS, HydroShare, https://doi.org/10.4211/hs.69f7d237675c4c73938481904358c789 LICENSE FOR USE -- MAPS AND DATA DISCLAIMER This resource is shared under the Creative Commons Attribution CC BY, http://creativecommons.org/licenses/by/4.0/ MAPS AND DATA DISCLAIMER The Oak Ridge National Laboratory (ORNL) shall not be held liable for improper or incorrect use of the data described or information contained on this map or associated series of maps. The data and related map graphics are not legal, land survey or engineering documents and are not intended to be used as such. ORNL gives no warranty, express or implied, as to the accuracy, reliability, utility or completeness of this information. The user of these maps and data assumes all responsibility and risk for the use of the maps and data. ORNL disclaims all warranties, representations or endorsements either express or implied, with regard to the information contained in this map product, including, but not limited to, all implied warranties of merchantability, fitness for a particular purpose or non-infringement. This preliminary map product is for research and review purposes only. It is not intended to be used for emergency management operational or life safety decisions at the local or regional governmental level or by the general public. Users requiring information regarding hazardous conditions or meteorological conditions for specific geographic areas should consult directly with their city or county emergency management office.

54 ENVIRONMENTAL SCIENCES↗

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↗

National Ignition Campaign Data Analysis by an ICF Outsider

A simple DT burn model, with no alpha-particle warming feedback, using a combination of ICF inputs and the measured ion temperature, T, volume compression, C, ablator velocity, v, and burn width time, Δt, can reproduce the national-ignition campaign (NIC) fusion yields with a rms percent difference of ~40%. Analysis suggestion the fusion-yield reproducibility of NIC events with nominally identically measured T, C, v, and Δt is <30%. The highest fusion yields, at a fixed laser energy are obtained with depleted-uranium (DU) hohlraums and high-density carbon (HDC) ablators. Small diameter fit tubes (≤5 μm) enhance fusion yields by an average of 26±15% relative to 10-μm diameter fit tubes. A dependence of fusion yield on the diameter of the capsule relative to the diameter of the hohlraum is inferred for fixed T, C, v, and Δt. No evidence of conditions close to ignition, or alpha-particle warming is found in the existing NIC data with LEH< 2MJ. To obtain a predictive model for extrapolation proposes we map the T, C, and v into power laws as a function of LEH. The above mentioned extrapolatable model has a rms percent difference from the yields of the second half of the NIC of (56 ± 4)%. Higher (and significantly lower) DT temperatures are obtained with HDC ablators (and low-foot drives). The highest compressions are obtained via low-foot drives, but this does not compensate for the corresponding lower temperatures. An extrapolation of the presented model suggests significant alpha-particle heating at an LEH of ~5 MJ. It is hoped that others more familiar it the conditions needed for ignition will use to presented model to estimate the LEH required for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗

Optimizing the Heisenberg Vortex Tube for Hydrogen Cooling (Final Technical Report)

Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.

08 HYDROGEN↗

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Pilot Testing of a Highly Efficient Pre-combustion Sorbent-based Carbon Capture System

TDA developed and demonstrated a highly efficient pre-combustion carbon capture system. The overall objective of this work was to develop a new sorbent-based pre-combustion carbon capture technology for Integrated Gasification Combined Cycle (IGCC) power plants. In this project our goal was to demonstrate the techno-economic viability of the new technology by 1) demonstrating it in large-scale slipstream tests, and 2) carrying out a high fidelity engineering and cost analysis. TDA’s process used an advanced physical adsorbent that selectively removes CO 2 from coal-derived synthesis gas above the dew point of the gas at temperatures as high as 300°C. The sorbent consists of a mesoporous carbon whose surface was grafted with functional groups that remove CO 2 via a well-known acid-base interaction. As documented in bench-scale experiments and field tests with actual coal gas, the sorbent achieved a very high capacity for CO 2 at temperatures as high as 300°C. The sorbent bound CO 2 more strongly than common physical adsorbents, providing the chemical potential needed for the high temperature operation. However, because CO 2 does not form a true covalent bond with the surface sites (as is the case with chemical absorbents), the sorbent regeneration could be carried out with only a very small energy input. The heat input to regenerate our sorbent was only 4.9 kcal per mol of CO 2 , which is much lower than that for chemical absorbents (e.g., 29.9 kcal/mol CO 2 for sodium carbonate) and was similar to the requirements of physical solvents (e.g., 4 kcal/mol CO 2 for Selexol TM ). Because the sorbent operates above the dew point of the synthesis gas (unlike the Selexol TM process), a higher power cycle efficiency can be achieved. With previous DOE/NETL funding (Contract No. DE-FE-0000469), we demonstrated the techno-economic viability of the technology in bench-scale tests and slipstream demonstrations at the National Carbon Capture Center (NCCC), Wilsonville, Alabama and Wabash River IGCC plant in Terra Haute, Indiana. We demonstrated a stable working CO 2 capacity for over 11,650 cycles with simulated synthesis gas. We also evaluated its performance with actual synthesis gas in two test campaigns at the Wabash River IGCC Plant, Terre Haute, IN and the National Carbon Capture Center (NCCC), Wilsonville, AL. The slipstream tests clearly showed that the actual coal gas constituents and potential contaminants (e.g. trace metals, halides, tars) had no effect on the sorbent’s ability to remove CO 2 (the same sorbent beds were used in both field tests with no sign of deactivation for 2,000 cycles with over 26,750 SCF of gas treated). As expected, due to the high temperature CO 2 removal capability and low energy needed to regenerate the sorbent, the power cycle efficiency with our process was greater than 34% on a higher heating value (HHV) basis; in comparison, the same IGCC plant equipped with the Selexol TM solvent for carbon capture can only achieve 31.4% HHV efficiency. The capital cost for an IGCC system with TDA’s process is estimated as $2,417/kW e , which is 12% lower than that of the IGCC/ Selexol TM process. The levelized cost for electricity including the transport, storage and monitoring (TS&M) cost for CO 2 was calculated as $\$ $92.9/MWh (lowest reported to our knowledge), which is much better than the $105.2/MWh estimated for the IGCC/ Selexol TM process. In this project (DE-FE0013105), TDA Research, in collaboration with our partners Gas Technology Institute (GTI), Illinois Clean Coal Institute (ICCI), University of California, Irvine (UCI), University of Alberta (UOA), Siemens, NCCC and Sinopec advanced the technical maturity of the technology; scaling it up by a factor of 100. We optimized the reactor design using computational fluid dynamics (CFD); using adsorption modeling we improved the pressure swing adsorption (PSA) cycle sequence. We carried out two field test campaigns with a fully-equipped 0.1 MW e prototype unit (for a total of 844 hours) using actual synthesis gas to prove the viability of the new technology. A successful 30 day (707 hrs) evaluation was completed at NCCC under air blown gasification conditions. We demonstrated 97.3% carbon capture at 1,500 SLPM, 93% carbon capture at 1,800 SLPM, and 90% carbon capture at 2,100 SLPM in the NCCC tests. We also demonstrated the system for 137 hours at a Sinopec petrochemical plant under oxygen blown gasification, demonstrating 86% carbon capture at 2,660 SLPM. In collaboration with University of California, Irvine (UCI), we completed a techno-economic analysis (TEA) for TDA’s warm gas cleanup technology integrated to IGCC power plant. The net plant efficiencies (on a coal HHV basis) for the warm gas cleanup cases were estimated to be 34.0% for E-GasTM gasifier, 34.4% for GE gasifier, 33.4 for the Shell gasifier and 34.2 for the TRIG TM gasifier (Cases 2, 4, 6 and 8 in this study) with a catalytic combustor for CO 2 purification, which are significantly higher than those for the Cold Gas Case, or an increase of as much as 12% in the heat rate for Case 2, 6% for Case 4, 9% for Case 6, and 9% for Case 8. The 1st year cost of electricity with the transport, storage and monitoring (TS&M) costs for the CO 2 included was $\$ $129.2/MWh for the E-GasTM gasifier Warm Gas Cleanup Case, $\$ $131.9/MWh for the GE gasifier Warm Gas Cleanup Case, $\$ $146.8/MWh for the Shell Gasifier Warm Gas Cleanup Case, and $\$ $129.9/MWh for the TRIG TM gasifier Warm Gas Cleanup Case. For comparison, the costs for the baseline Cold Gas CO 2 removal with Selexol for the different gasifiers were: $\$ $146.6/MWh for the E-Gas TM gasifier, $\$ $142.2/MWh for the GE gasifier, $\$ $159.0/MWh for the Shell gasifier and $\$ $144.3/MWh for the TRIG TM gasifier. In summary, the costs for our system were 7 to 12% lower than the corresponding Cold Gas Cleanup cases. The results of this techno-economic analysis suggested that TDA’s high temperature PSA-based Warm Gas Clean-up Technology can make a substantial improvement in the IGCC plant thermal performance for achieving near zero CO 2 emissions for E-Gas TM , GE, Shell and TRIG TM gasifier based IGCC power plants. The capital expenses were estimated to be lower than that of Selexol’s™. Taken together, the higher net plant efficiency and lower capital and operating costs resulted in substantial reduction in the cost of carbon capture for the IGCC plant equipped with TDA’s high temperature PSA-based carbon capture system. Finally, in collaboration with Gas Technology Institute (GTI) we completed the environmental health and safety assessment for TDA’s warm gas carbon capture technology.

01 COAL, LIGNITE, AND PEAT↗

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Catalytic Bio-crude Production in a Novel, Short-Contact Time Reactor

RTI has developed a catalytic fast pyrolysis (CFP) process to produce biomass-based hydrocarbon replacement fuels. From a technology perspective, this advanced biofuels process produces liquid transportation fuels that can leverage capital expenditures in the existing petroleum distribution infrastructure. From a business perspective, this technology does not face potential market limitations caused by oxygenate blend limitations for light duty vehicle applications. A promising catalyst has been developed with excellent potential for deoxygenation to produce a hydrocarbon-rich biocrude that is more suitable for upgrading with traditional hydroprocessing technology. The role of the catalyst is to control the chemistry during biomass pyrolysis to produce a biocrude in a single step that has lower oxygen content and is more thermally stable than conventional biomass fast pyrolysis oil. Catalyst properties are optimized to minimize gas and coke production and improve catalytic deoxygenation and biocrude yields. Results from this project have demonstrated that biocrudes from catalytic fast pyrolysis contain <20 wt% oxygen, with the organic fraction as low as 12% oxygen. These oxygen numbers are > 50% lower compared to traditional fast pyrolysis oils. The biocrude produced from our CFP process can be upgraded to gasoline range hydrocarbons. A very conservative approach for upgrading included a preliminary stabilization step followed by a hydroprocessing step. The oxygen content in the upgraded product was decreased by 97% (to less than 3wt %) from the oxygen content of the biocrude sample. Additionally, the acid content was reduced by 99% from the acid content of the biocrude sample. In an optimized integrated process, the CFP step and the hydroprocessing (upgrading) step need to be carefully balanced to maximize the input biomass energy that is recovered in the finished biofuels. Hydrogen utilization and carbon recovery during both catalytic biomass pyrolysis and biocrude upgrading need to be optimized to maximize energy recovery in the finished biofuel. A nominal 1 ton per day (TPD) catalytic biomass pyrolysis system was designed based on the laboratory experimental data. The system was designed to achieve the short residence times (0.5–2 sec) and high heat transfer rates required for maximum liquid biocrude yields while optimizing process integration to maintain catalyst activity by continuous regeneration. The catalytic biomass pyrolysis reactor is a continuously circulating single-loop transport reactor design that is flexible enough to allow sensitivity studies around temperature, residence time, biomass feed rate, and catalyst-to-biomass ratio (i.e., catalyst circulation rate) for process optimization. For a nominal 100 lb/hr feed rate of dry biomass, the products include ~25 lb/hr char and coke, ~20 lb/hr permanent gases, and ~55 lb/hr condensable vapors. Pyrolysis temperatures are 350–600°C, and the system is operated close to ambient pressure. Production of the catalyst developed within the project was scaled up to provide enough material (200kg) for bench-scale testing. The 1TPD bench-scale unit was fabricated, installed and commissioned with this project. Biocrude samples from wood were produced in this unit with the CFP catalyst with material balance > 90%. The intermediate biocrude yields and chemical composition compare well with previous laboratory results and a hydroprocessing strategy for upgrading the resulting biocrude was developed. The experimental results obtained within this project were used to develop a process model for a 2000 TPD commercial facility. Detailed flow sheets were developed for all major unit operations in the process and provided a basis for an economic analysis to develop capital and operating costs. Estimates of the biofuel production costs were used to develop internal rates of return for the yields and compositions measured for the current state of technology at the end of this project and for the yields anticipated for an optimized (nth plant) scenario as additional process improvements are made.

09 BIOMASS FUELS↗

Cyote-attack Chain Estimator

Attack Chain Estimator (ACE) Application Overview The Attack Chain Estimator (ACE) Application is a sophisticated tool designed for the ingestion, classification, sequencing, and enrichment of cybersecurity threat reports. This application leverages advanced machine learning models and extensive historical data to provide comprehensive insights into cyber threats, specifically targeting Industrial Control Systems (ICS). Purpose The primary functions of the ACE Application include: Ingestion of Cybersecurity Threat Reporting: Capable of ingesting text-based threat reports in markdown or text file format. Supports ingestion of structured data from other sources in STIX/JSON format. Classification of Report’s Text-Based Events: Utilizes a DeBERTa classifier, specifically trained on cybersecurity data, to map the events to MITRE ATT&CK for ICS Tactics and Techniques. Classification is performed using multiple Jupyter notebooks and machine learning workflows hosted as FastAPI microservices: regex_data deberta_base_35_train_hft_classifier_mlflow.ipynb hft_regex_classifier_mlflow.ipynb param_train_hft_classifier_mlflow.ipynb regex_tactic_tech.ipynb Ordering of Tactics, Techniques, and Observable Events: Sequences the identified tactics, techniques, and events to form a coherent attack chain. Enrichment with Historical Attack Chain Details: Enhances the attack chain with details from historical attacks using a Markov model developed from CyOTE Precursor Analysis Report data. The Markov model is available as a FastAPI endpoint for seamless integration. Enrichment with Adversary Emulation Capabilities Data: Integrates adversary emulation capabilities data using MITRE Caldera for OT adversary abilities UUIDs. Export of Output Files: Provides options to export the enriched attack chain in JSON or CSV formats. Routing of Output to Other Applications: Facilitates routing of output to various platforms and applications, including: Threat Intelligence Platforms COREII Scout for Threat Intelligence Analysis COREII Modeling and Simulation for Adversary Emulation Technical Description The ACE Application is an advanced cybersecurity tool designed to provide detailed threat analysis and sequence generation. It is built on a robust architecture that integrates natural language processing, machine learning, and historical data modeling. Key Components: Data Ingestion Module: Handles the input of threat reports and data from various formats, ensuring flexibility in data sources. Classification Engine: Employs DeBERTa-based classifiers hosted as FastAPI microservices to analyze and classify threat report events in accordance with the MITRE ATT&CK framework for ICS. Sequence Generator: Orders the classified events into a logical attack chain, providing clear insight into the sequence of tactics and techniques used in the threat. Enrichment Engine: Integrates historical data and adversary emulation capabilities to enhance the attack chain with valuable context and additional details. The historical data enrichment is powered by a Markov model, which is available as a FastAPI endpoint. Export and Routing Module: Facilitates the export of the enriched attack chain in multiple formats and routes the output to designated applications for further analysis or emulation.

Paul, Tony [Idaho National Laboratory (INL), Idaho↗

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Process with Decoupled Absorber Kinetics and Solvent Regeneration through Membrane Dewatering and In-Column Heat Transfer (Final Report)

This report summarizes the work conducted on project DE-FE0031604 where University of Kentucky Center for Applied Energy (UK CAER) has validated its intensified CO 2 capture process through substantial enhancements to the kinetics of the absorption process and energy reductions by absorber temperature profile modification, dewatering and heat integration technologies for achieving significant capital and operating cost reductions. To address DOE’s objective of improving post-combustion CO 2 capture technology and reducing associated cost, UK CAER employed an intensified process which combined three key aspects targeted at overcoming inherent limitations or barriers in the conventional CO 2 capture and desorption process. The process designed to be independent of the type of solvent used, included (1) the use of 3-D printed two-channel structured packing material to control the temperature profile and increase the CO 2 absorption rate in the absorber, (2) a zeolite membrane dewatering unit for dewatering of the carbon-rich solvent to decouple solvent concentration needs for CO 2 absorption and desorption, and (3) a rich-split feed with two-phase flow heat transfer prior to the stripper that provided a secondary point of vapor generation to provide energy savings in steam extraction and solvent regeneration. The project was executed over two budget periods. This involved testing of individual process components which included the advanced heat transfer packing and the dewatering membrane on UK CAER’s 30 liter per minute (L/min) CO 2 (3” Column) capture bench unit with simulated flue gas in the first budget period. Subsequent scaled-up testing of these components together with the split-feed configuration were also tested in UK CAER’s 0.1 MWth CO 2 capture unit with coal-derived flue gas in the second-budget period. Long term studies were done during this period to assess process and solvent performance over extended duration. Project partners Lawrence Livermore National Laboratory (LLNL) and Media & Process Technology (MPT) led the development efforts for the advanced packing material and dewatering membrane respectively. Data from the long term testing was used as input for an environmental, health and safety (EH&S) assessment for the process and scaled technology performed by ALL4 LLC. Trimeric Corporation also completed a techno-economic analysis (TEA) for the UK CAER technology which was compared to the DOE reference Case B12B. Tests on the 3” column capture unit showed that the advanced heat transfer packing could be used to lower the bulge temperature in the absorber, and this was also proven in the scaled testing in the 0.1 MWth CO 2 capture unit. The bulge temperature could be lowered by >10 °C, changing the temperature profile in the absorber, and showed potential to enhance absorption with the ability to tailor the profile to provide conditions suitable for a solvent’s properties and kinetics. Conditions for short term evaluation of a 19” zeolite dewatering membrane on the 3” column capture unit yielded desirable fluxes and sustained rejection rates of >80%. However, for the scaled testing of six membrane modules consisting of 21 parts of 31-inch-long membrane tubes in each bundle (surface area 0.3 m 2 ), over a more extended duration, similar rejection rates could not be achieved. With the split-feed of the rich stream to the stripper, improved heat recovery minimized waste heat exiting the top of the stripper. The stripper exhaust temperatures could be reduced by >10 °C; reducing the amount of water vaporization contributed to lowering the regeneration energy by ~ 15%. The energy benefit could be sustained from the long term monitoring of the solvent performance. The solvent properties were not significantly impacted over the long-term operations. The benefits of the UK CAER process demonstrated experimentally were mostly validated from the TEA comparing a commercial scale application of the technology to DOE reference Case B12B. The cost of CO 2 capture for the UK CAER technology was estimated to be ~$34.97/tonne of CO 2 captured; a reduction of 23% compared to Case B12B. The increase in cost of electricity was also shown to be 16% lower than that of Case B12B. The total parasitic demand was also shown to be 11% lower. The key drivers for the benefits are a result of the process intensification approaches employed in the UK CAER technology for enhanced solvent performance, effective heat recovery and improved energy performance. The EH&S assessment did not find any major environmental concerns or barriers to the full scale implementation of the technology.

20 FOSSIL-FUELED POWER PLANTS↗

Modeling dust mineralogical composition: sensitivity to soil mineralogy atlases and their expected climate impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ~ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ~ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

54 ENVIRONMENTAL SCIENCES↗

Covariance Shaping Over Riemannian Manifolds for Massive MIMO Communication

Acquiring accurate instantaneous channel state information (CSI) is a challenging aspect of massive multi-input multi-output (MIMO) communication. Utilizing statistical information, such as channel covariance matrix, to design statistical beamforming vectors is robust when compared to instantaneous CSI. In this paper, we propose a novel MIMO covariance shaping scheme over Riemannian manifolds. It serves as an effective statistical beamforming solution to a number of close proximity user equipment (UE) that are undergoing substantial channel correlation. Proposed algorithm exploits the Hermitian positive definite nature of covariance matrices lying over Riemannian manifold. We introduce Wasserstein distance function as a Riemannian metric to measure distances between channel covariance matrices. Furthermore, K-means clustering technique is utilized to effectively identify the optimal shape of effective optimal covariance matrices. Our findings suggest that maximizing the geodesic distance between covariance matrices ultimately leads to a corresponding increase in the network throughput, as determined by the beamforming vector used to shape the covariance matrices. Simulation results validate that the proposed solution converges faster than Euclidean-based state-of-the-art, while maintaining the same computational complexity. Finally, the sum rate performance asymptotically achieves full capacity for two-UE case and more than 96% of the upper bound exhaustive search benchmark for multi-UE scenario.

42 ENGINEERING↗