SEARCH · Engineering Papers

Results for “Operator approximation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Investigation of Radiological Contamination at Parks in the Gateway National Recreation Area - 20350

In 2005, the National Park Service was notified that elevated radiation levels were detected at Great Kills Park during a baseline radiological mapping effort conducted by the Counter Terrorism Bureau of the New York City Police Department and the United States Department of Energy. Great Kills Park is in the Staten Island Unit of the Gateway National Recreation Area. A records search determined that approximately 280 acres of the park was established over an area in which the New York City Department of Sanitation placed approximately 15.1 million cubic yards of waste fill from 1944 to 1948. Since 2005, a series of response actions and investigations were conducted on behalf of the National Park Service to mitigate risk to human health and the environment and to determine the nature of radiological contamination at Great Kills Park. The findings support that radiological contamination is related to certain components of waste fill, including incinerator residue from both on and offsite incinerators, coal ash from the historical use of coal-fired boilers for heating, and incidental radiological artifacts in the waste fill due to their historical use and subsequent disposal. A majority of the radiological artifacts were determined to contain radium-226 and included radioluminescent markers, radium impregnated switch board buttons, radium paint containers, and brachytherapy devices. The National Park Service has several other Parks in the Gateway National Recreation Area that were established in areas with a similar history to Great Kills Park, and overlie waste fill historically placed during landfilling operations, including Spring Creek Park and Dead Horse Bay. Over the last year, the National Park Service has performed additional radiological investigations at Great Kills Park as part of the ongoing Remedial Investigation under the Comprehensive Environmental Response, Compensation, and Liability Act. In addition, preliminary radiological investigations were performed at Spring Creek Park and Dead Horse Bay Park in the Jamaica Bay Unit of Gateway National Recreation Area. These investigation activities identified gamma radiation levels above the ambient level, several of which were attributable to radiological artifacts similar in nature to those recovered at Great Kills Park. Ongoing investigations at Parks in the Gateway National Recreation Area support an association between radiological contamination and historical landfilling practices that were common prior to the regulation and control of radioactive materials. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

AquaPV: Regulatory and Environmental Considerations for Floating Photovoltaic Projects Located on Federally Controlled Reservoirs in the United States

To meet the nation's decarbonization goals, the U.S. Department of Energy's Solar Futures study forecasts that installed solar photovoltaic (PV) capacity must increase nearly tenfold, from 80 gigawatts (GW) in 2020 to approximately 760 GW cumulative installed capacity by 2035. Ground-mounted PV is expected to dominate future solar deployment and will require more than 3.5 million acres of land to meet annual demand projections (of nearly 45 GW) by 2030. However, various competing demands for land (e.g., agricultural production, conservation) and high land acquisition costs in specific locations could be challenges to meeting future PV demand solely with ground-mounted PV deployment. Floating photovoltaics (FPV) may be an alternative in locations where ground-mounted PV is not feasible and aid in reaching the nation's PV deployment and decarbonization goals. FPV is a newer siting approach in which a PV array is affixed to a floating apparatus and sited on a water body like a reservoir behind a dam. FPV systems may be stand-alone or co-located at new or existing hydroelectric facilities or pumped storage hydropower (PSH) facility reservoirs. Co-located FPV systems may or may not be operationally paired and work in tandem with the hydroelectric or PSH facility. This report provides novel analysis to understand the opportunities and challenges associated with developing stand-alone and co-located FPV projects on reservoirs in the United States. Specifically, the report explores potential environmental and energy benefits and environmental impacts associated with the siting, construction, and operation of FPV projects. The report also identifies and analyzes U.S. federal- and state-issued permits and authorizations required by federal laws to understand the licensing pathways and regulatory requirements for FPV projects sited on reservoirs licensed by the Federal Energy Regulatory Commission and on powered and non-powered reservoirs owned by the Bureau of Reclamation or U.S. Army Corps of Engineers.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗

DOE OSTI2024

AquaPV: Regulatory and Environmental Considerations for Floating Photovoltaic Projects Located on Federally Controlled Reservoirs in the United States

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗

DOE OSTI2025

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Comparative study on the formation of Cr and Ti ohmic contacts to (001) β -Ga 2 O 3

Here, a comparative study of Cr/Au and Ti/Au ohmic contacts on (001) β-Ga 2 O 3 was conducted. The electrical behavior from current-voltage measurements and the interfacial composition and microstructure as determined from high-resolution transmission electron microscopy (TEM) with energy dispersive x-ray analysis were compared for the different contacts at selected points in an annealing series (300–700 °C, 1 min. anneals in N 2 ). Cr/Au contacts became ohmic at temperatures (300–350 °C) approximately 50–100 °C lower than Ti/Au contacts (400–450 °C). Cr/Au and Ti/Au contacts demonstrated optimal ohmic behavior (lowest resistance) when annealed to 450–500 and 500–600 °C, respectively, with Ti/Au contacts yielding a lower total resistance than Cr/Au. Cross-sectional TEM images of Cr/Au contacts annealed at 450 °C revealed the presence of Au nanoclusters at the Ga 2 O 3 interface and CrO x layers at both the top of the contact and the Ga 2 O 3 interface. Whereas TiO x also formed at the top and bottom interfaces in 450 °C-annealed Ti/Au contacts, the TiO x surface layer appeared to be variable in thickness and/or discontinuous, unlike the CrO x surface layer. Au nanoclusters were not detected at the interface in the Ti/Au contacts. The interdiffusion and oxidation observed in both contact metallizations point to the need for diffusion barriers that may allow these contacts to be used in future Ga 2 O 3 -based devices that operate at elevated temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Validation Testing for Molten Chloride Reactor Experiment Equipment Removal and Disposal Techniques

The Molten Chloride Reactor Experiment (MCRE) will be the first reactor featuring a fast-spectrum molten chloride circulating nuclear fuel system in the world. Planning for equipment removal and disposal (ERD) of MCRE has identified several technology gaps due to the unique environment of this nuclear experiment. Some of the gaps arise from the application of existing disassembly and/or sizing methods to novel material forms or in novel configurations. Others arise from unknown material behavior. This paper summarizes proposed test plans for ERD validation experiments to address these complicated or unknown equipment removal procedures. At the Waste Management Symposia in 2024, the Idaho National Laboratory (INL) MCRE ERD team presented the challenges associated with hosting multiple nuclear experiments in series with only brief transition periods between systems. Such difficulties include higher dose rates, the presence of radioisotopes infrequently encountered in reactor decommissioning and radioactive waste management, lack of intrinsic remote-operations infrastructure in the test bed, space constraints in the test bed, and contamination minimization requirements. To address these challenges, remote or semi-remote technologies are planned to be implemented in a non-hot cell environment with limited space availability. The team also discussed how a systems engineering approach is being used for conceptual development and design of equipment removal systems to address these challenges. For example, to reduce constraints for the removal of more difficult components, non-activated, noncontaminated elements are planned to be taken out first where possible. Still, there are complexities associated with the remaining components. In this work, the operational framework for MCRE ERD was reviewed for technical gaps and open questions, and test plans were drafted to address these areas. The tests plans were written for the following categories: vision systems, pipe cutting, drill/grout/filler, flush salt, and miscellaneous, with the miscellaneous group consisting of tests like techniques for removing bearings and reflector bricks. The test plans explore material, infrastructure, and staffing requirements needed for test execution. The test plans additionally focus on the evaluation of success. Determining the outcome of a test is imperative - as these explorative actions have the potential to rearrange or re-scope planned ERD activities. Success criteria identified thus far include required tool output, required area(s), debris production and mitigation, and repeatability. Test plans are an essential aspect of the systems engineering approach to MCRE ERD. They are used as the beginning steps in defining use cases for the ERD system. Performance of the validation tests is expected to begin in the summer of 2025 and will take approximately 9 to 12 months to complete. Execution of these plans will be expedited by specifying test needs ahead of time, facilitating efficient interactions with any subcontractors tasked with running the requested tests. Evaluating the outcomes of these tests will inform MCRE ERD procedures and timing and will also identify additional technical constraints for the MCRE ERD System. This upfront process optimization effort will help the project save time and resources at the end of the experiment.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2024

FEED Study of CarbonCapture Inc DAC and CarbonCure Utilization Technologies Using United States Steel’s Gary Works Plant Waste Heat (Final Report)

The University of Illinois at Urbana-Champaign (UIUC) led this project to produce a front-end engineering design (FEED) study of an advanced Direct Air Capture and Utilization System (DACUS) system that can remove a minimum of 5,000 tonnes/yr net of carbon dioxide from air (based on cradle-to-gate LCA) and utilizing the CO 2 to produce low carbon intensity ready mix concrete. The designed system, if built, would be larger than any currently existing Direct Air Capture (DAC) collector in the U.S. Such carbon capture technologies are critical to meeting the goals of the DOE’s program to accelerate climate-critical technology. In addition to the power sector, industrial facilities for the manufacture of steel and cement/concrete are among the major sources of anthropogenic CO 2 . DAC is a promising new technology for reducing CO 2 , a potent greenhouse gas, in the atmosphere but is expensive, in part due to the energy required to adsorb and desorb captured CO 2 during cycles. By integrating CarbonCapture Inc. (CCI) DAC modules at United States Steel's Gary Works (USS) and utilizing the site's waste heat, energy, and location this project evaluates the feasibility of utilizing the captured CO 2 and the logistics of transportation. CarbonCapture Inc. has developed an innovative DAC system using novel adsorbents to cost-effectively capture CO 2 . The captured, liquified gas will be trucked to ready-mix concrete plants within the region, the closest of which is approximately 3.5 miles away, where CarbonCure will inject it into concrete during the mixing process at the facilities. The carbon dioxide reacts with concrete, mineralizing into calcium carbonate (CaCO 3 ), permanently locking the greenhouse gas into the matrix of the building material. This FEED study demonstrated a full CO 2 value chain for DACUS from industrial facilities. It also provided a means for Visage Energy Corp. (Visage) to assess the impact of this holistic approach on job creation, regional economic impact, and environmental justice. The project team also included Sargent & Lundy (S&L) to provide the constructability review and costing of the integration of the DAC with the steel plant. Ecotek Engineering USA, LLC designed the outside battery limit (OSBL) infrastructure to connect the DAC and the plant. Activities performed during the project included: (1) Project Management Plan; (2) Technology Maturation Plan (TMP); (3) Initial Workforce Readiness Plan; (4) Workforce Readiness Plan; (5) Front-End Engineering Design (FEED) Study; (6) Project Design Basis; (7) Hazards and Operability (HAZOP) Study; (8) Constructability Review; (9) Project Cost Assessment; (10) Logistics Analysis of CO 2 Transportation to the Utilization Site; (11) Business Case Analysis; (12) Life Cycle Analysis (LCA); (13) Environmental Health and Safety (EH&S) Analysis; (14) Environmental Justice Analysis; and (15) Economic Revitalization and Job Creation Outcomes Analysis. This report provides a summary of the outcomes and results of the project, which was performed between Oct. 1, 2022, through Sept. 30, 2024.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Production of Fischer-Tropsch Synfuels at Nuclear Plants

A case study analysis was performed to evaluate nuclear-powered synthetic fuel production in the midwestern United States (U.S.). A Fischer-Tropsch (FT) fuel synthesis plant design was used as the basis for the analysis. The FT plant design was configured to produce a product slate consisting of diesel fuel, jet fuel, and motor gasoline blend stocks from carbon dioxide (CO 2 ) and hydrogen (H 2 ) feedstocks. The CO 2 feedstock for the FT plant was assumed to be sourced from biorefineries in the region around a Midwest light water reactor (LWR) nuclear power plant (NPP). The analysis specifies that power from the LWR is used to produce H 2 via high-temperature steam electrolysis and to operate the FT synfuel production plant. Capital costs were estimated for the FT plant while capital costs for the electrolysis plant were based on previous Idaho National Laboratory (INL) studies. In addition to labor and maintenance costs for the FT and electrolysis plants, operating costs also include the costs for CO 2 feedstock transport. An analysis was performed to determine the cost of transporting CO 2 from the distributed biorefinery sources to the centralized fuel synthesis plant as a function of the synfuel plant capacity and corresponding CO 2 demand. The primary revenue streams are associated with sales of the synthetic fuel products. The synthetic fuel products will likely follow the same market trends as the conventional fuel products. The synfuel price data was thus based on projections made by the U.S. Energy Information Administration (EIA) 2021 Annual Energy Outlook (AEO) for conventional fuel products minus federal and state taxes, as well as marketing and distribution costs. The economic analysis also considered cases that included and excluded revenues from the 2022 Inflation Reduction Act (IRA) clean hydrogen production tax credit (PTC) of $\$ $3.00/kg for the first ten years of operation. The economic analysis calculated the net present value (NPV) for cases involving steady-state synfuel production for comparison with the NPV for a business-as-usual case in which NPP continues to sell only electric power to the grid. A synfuel production “Reference Case” was considered in addition to sensitivity cases in which the plant capacity, electricity price, and synthetic fuel product prices were perturbed. The synfuel production Reference Case considered a scenario in which the electrolysis and synfuel plants utilized a combined electrical load of 1000 megawatt electrical (MWe) from the LWR with the balance of the LWR power output being sold to the electric grid. The economic analysis suggests that the synfuel production Reference Case evaluated in this analysis would lead to considerable economic potential for near-term deployment of a nuclear-based synfuel production plant. Specifically, the economic analysis suggests that the deployment of a 1000 megawatt (MW) nuclear-powered synfuel plant could result in a NPV increase of approximately $\$ $1.7 billion for a case with no clean synfuel price premium relative to conventional petroleum fuels when accounting for the additional revenues from the 2022 IRA clean hydrogen PTCs of $\$ $3/kg. Sensitivity analysis was performed to evaluate the effect of perturbation of selected model input parameters on the NPV for the synfuel production Reference Case. The sensitivity analysis indicates that the plant capacity has the largest impact on the differential NPV, with a smaller synfuel production capacity resulting in a decrease in revenue when a larger fraction of the power from the NPP is sold to the grid and a smaller fraction of the power is used to produce synthetic fuel products. The synfuel product pricing has the next largest impact on the differential NPV, with lower synfuel prices resulting in decreased NPV from decreased synfuel sales revenue while higher synfuel prices result in increased NPV from increased synfuel sales revenue. Electricity pricing has a smaller effect on the NPV than the fuel sales price since, in the Reference Case, most of the energy from the NPP is used for synfuel production and a smaller amount of the system revenues are associated with electrical power sales. However, the electricity price sensitivity does indicate that the Synfuel Integrated Energy System (IES) would have a greater NPV than the business-as-usual case (e.g., grid power sales only) when electricity market prices are low, suggesting that synfuel production could provide a strategy for decreasing the economic risks to NPPs posed by a loss of revenues attributed to falling electricity market prices.

10 SYNTHETIC FUELS↗

DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Torrefaction of almond shell as a renewable reinforcing agent for plastics: techno-economic analyses and comparison to bioethanol process

Abstract In 2016, the US state of California alone produced nearly 3.5 billion kilograms of almonds, accounting for approximately 84% of the world’s almond production. This generated about 2.58 million metric tons (MTs) of almond residues. Almond shells are currently either burned to generate power or disposed of in landfill. Valorizing almond shells and hulls provides an opportunity to replace petroleum-derived products and divert organic material from landfill. Here we demonstrate a detailed techno-economic analysis (TEA) of an almond shell torrefaction process capable of utilizing the 520 000 MTs of almond shells produced annually in California. Our process also includes preprocessing the torrefied biomass to exploit it as a reinforcing agent for plastics. We further compared the revenue generated from the torrefied biomass and bioethanol derived from the same quantity of almond shells. We considered three different torrefaction facility scales to evaluate trade-offs between economies of scale at the facility and trucking costs to deliver almond shells. A facility that takes in 200 000 MT yr –1 of almond shells results in lower per-unit-output basis capital and operating cost relative to other smaller-scale torrefaction facilities, including 10 000 MT yr –1 and 50 000 MT yr –1 , considered for analysis in this study. The large-sale facility results in a minimum selling price (MSP) of the torrefied biomass of $311.4 MT –1 . An analogous TEA on converting almond residues into bioethanol is also investigated. The MSP of almond shell derived ethanol ($1.71 kg −1 ) is higher than that of corn ($0.48 kg −1 ) or cellulosic biomass ($0.88 kg −1 ) derived ethanol. Compared with the bioethanol route, the torrefied almond shells result in three times more revenue if utilized as a reinforcing agent for plastics.

Liu, Dupeng (ORCID:0000000269514064)↗

DOE OSTI2020

Machine Learning Algorithms for Automated NIF Capsule Mandrel Selection

Small shells, approximately 2 mm in diameter, made from Poly(α-methylstyrene) (PAMS) are used as mandrels in the production of glow discharge polymer capsules located at the center of inertial confinement fusion experiments. The visual inspection process of microscope images of these shell mandrels, including detection of micron-sized defects on the shell surface as well as the handling and sorting, is a very labor-intensive, repetitive, and highly subjective process that stands to benefit greatly from automation. As part of an effort to decrease the number of labor hours spent in capsule handling, inspection, and metrology, the development of robotic systems was presented in a paper by Carlson et al., “Automation in Target Fabrication” [Fusion Sci. Technol., Vol. 70, p. 274 (2016)]. The current work expands the automated image acquisition systems developed previously and adds the use of convolutional neural networks to select capsules best suited for use in the downstream production process. Through the use of these machine learning algorithms, the selection process becomes robust, repeatable, and operator independent. Finally, as an added benefit the system developed as part of this work is able to provide defect statistics on entire shell batches and feed this information upstream to the production team.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Emergency department overcrowding and its associated factors at HARME medical emergency center in Eastern Ethiopia

Introduction: Emergency department (ED) overcrowding has become a significant concern as it can lead to compromised patient care in emergency settings. Various tools have been used to evaluate overcrowding in ED. However, there is a lack of data regarding this issue in resource-limited countries, including Ethiopia. This study aimed to validate NEDOCS, assess level of ED overcrowding and identify associated factors at HARME Medical Emergency Center, located in Hiwot Fana Comprehensive Specialized Hospital, Harar, Ethiopia. Methods: A cross-sectional study was conducted at the HARME Medical Emergency Center, Hiwot Fana Comprehensive Specialized Hospital, involving a total of 899 patients during 120 sampling intervals. The area under the receiver operating characteristic curves (AUC) was calculated to evaluate the agreement between objective and subjective assessments of ED overcrowding. A multivariable logistic regression analysis was employed to identify factors associated with ED overcrowding and statistically significant association was declared using 95% confidence level and a p-value < 0.05. Results: The interrater agreement showed a strong correlation with a Cohen's kappa (κ) of 0.80. The National Emergency Department Overcrowding Study Score demonstrated a strong association with subjective assessments from residents and case team nurses, with an AUC of 0.81 and 0.79, respectively. According to residents' perceptions, ED were considered overcrowded 65.8% of the time. Factors significantly associated with ED overcrowding included waiting time for triage (AOR: 2.24; 95% CI: 1.54–3.27), working time (AOR: 2.23; 95% CI: 1.52–3.26), length of stay (AOR: 2.40; 95% CI: 1.27–4.54), saturation level (AOR: 2.35; 95% CI: 1.31–4.20), chronic illness (AOR: 2.19; 95% CI: 1.37–3.53), and abnormal pulse rate (AOR: 1.52; 95% CI: 1.06–2.16). Conclusion: The study revealed that ED were overcrowded approximately two-thirds of the time.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Nested active learning for efficient model contextualization and parameterization: pathway to generating simulated populations using multi-scale computational models

There is increasing interest in the use of mechanism-based multi-scale computational models (such as agent-based models (ABMs)) to generate simulated clinical populations in order to discover and evaluate potential diagnostic and therapeutic modalities. The description of the environment in which a biomedical simulation operates (model context) and parameterization of internal model rules (model content) requires the optimization of a large number of free parameters. In this work, we utilize a nested active learning (AL) workflow to efficiently parameterize and contextualize an ABM of systemic inflammation used to examine sepsis. Contextual parameter space was examined using four parameters external to the model’s rule set. The model’s internal parameterization, which represents gene expression and associated cellular behaviors, was explored through the augmentation or inhibition of signaling pathways for 12 signaling mediators associated with inflammation and wound healing. We have implemented a nested AL approach in which the clinically relevant (CR) model environment space for a given internal model parameterization is mapped using a small Artificial Neural Network (ANN). The outer AL level workflow is a larger ANN that uses AL to efficiently regress the volume and centroid location of the CR space given by a single internal parameterization. We have reduced the number of simulations required to efficiently map the CR parameter space of this model by approximately 99%. In addition, we have shown that more complex models with a larger number of variables may expect further improvements in efficiency.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Development of a framework for sequential Bayesian design of experiments: Application to a pilot-scale solvent-based CO 2 capture process

In this paper, a methodology is developed for sequential design of experiments (SDoE) for process systems and applied to a solvent-based CO 2 capture system. In this approach, the prior knowledge of the system is used to prioritize process data collection at specific operating conditions. These data are then incorporated into a Bayesian inference methodology for updating a stochastic model by refining estimations of its underlying parameters, and the updated model is then used to generate the next set of test runs. Thus, the new knowledge obtained from the data is used to guide subsequent iterations of the experimental runs, ensuring that the overall data collection is maximally informative given that most experimental campaigns, especially at pilot or higher-scale plants, are costly, time-consuming, and resource-limited. The test run objective for this work was to minimize the maximum model prediction uncertainty for key output variables, but the methodology is generic and can be readily applied to other test run objectives. This methodology is applied to an aqueous monoethanolamine (MEA) pilot plant campaign at the National Carbon Capture Center (NCCC) in Wilsonville, Alabama, USA. The SDoE framework was utilized for two iterations, while collecting 18 sets of data representing different process conditions, and this resulted in an overall average reduction in uncertainty of approximately 50% in the prediction of CO 2 capture percentage. Moreover, 11 additional data sets were obtained with variation of absorber packing height for further model validation. This work shows the capability of the SDoE framework to maximize learning given limited resources, allowing for the reduction of model uncertainty, which is of great importance for many applications including reduction of technical risk associated with scale-up and economic analysis.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2020

Online Energy-optimal Routing for Electric Vehicles with Combinatorial Multi-arm Semi-Bandit

We have seen a rapid growth in the adoption of electric vehicle in today’s commercial mobility service market, from carrying passengers to the delivery of goods. One cardinal issue concerning the operation of EVs is the optimal routing of EV fleets with limited battery capacity. In this study, we investigate the energy-optimal online routing problem for the fleet of EVs, which focuses on identifying real-time minimum electricity consumption paths (MECP) for multiple OD pairs with limited information. We develop a multi-OD combinatorial multi-arm semi-bandit model (MCMAB) that uses the fleet of EVs as sensors in the transportation network and promotes the utilization of common information shared by different OD pairs. We further enrich the model with the path elimination policy to obtain MECP of high confidence while significantly reducing the number of learning iterations and the number of explorations needed. We demonstrate the effectiveness of the MCMAB and the efficiency of the path elimination policy with comprehensive numerical experiments in Manhattan, NYC. The results show that the proposed online routing algorithms can achieve near-optimal MECPs efficiently, and the quality of the solutions is significantly better than using the shortest travel time paths as approximate MECPs.

Chen, Xiaowei↗