Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “joint state 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 91 records · Page 5

Overview of Preliminary Concepts for a Port of Alaska Resilient Microgrid

The Port of Alaska in Anchorage enables the economic vitality of the Municipality of Anchorage and State of Alaska. It also provides significant support to defense activities across Alaska, especially to the Joint Base Elmendorf-Richardson (JBER) that is immediately adjacent to the Port. For this reason, stakeholders are interested in the resilience of the Ports operations. This report documents a preliminary feasibility analysis for developing an energy system that increases electric supply resilience for the Port and for a specific location inside JBER. The project concept emerged from prior work led by the Municipality of Anchorage and consultation with Port stakeholders. The project consists of a microgrid with PV, storage and diesel generation, capable of supplying electricity to loads at the Port a specific JBER location during utility outages, while also delivering economic value during blue-sky conditions. The study aims to estimate the size, configuration and concept of operations based on existing infrastructure and limited demand data. It also explores potential project benefits and challenges. The report goal is to inform further stakeholder consultation and next steps.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cosmological parameter estimation with a joint-likelihood analysis of the cosmic microwave background and big bang nucleosynthesis

Here, we present a joint-likelihood analysis of big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) data, consistently combining likelihoods and taking into account uncertainties in nuclear reaction rates for the first time. Bayesian inference is performed on the baryon abundance and the effective number of neutrino species, 𝑁 eff , using a CMB Boltzmann solver in combination with LINX , a new flexible and efficient BBN code. We marginalize over Planck nuisance parameters and nuclear rates to find 𝑁 eff =3.0⁢8$^{+0.15}_{−0.14}$, 2.9⁢4$^{+0.16}_{−0.15}$, or 2.96$^{+0.13}_{−0.14}$, for three separate reaction networks. This framework enables robust testing of the lambda cold dark matter paradigm and its variants with CMB and BBN data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tail Dependence as a Measure of Teleconnected Warm and Cold Extremes of North American Wintertime Temperatures

Current models for spatial extremes are concerned with the joint upper (or lower) tail of the distribution at two or more locations. Such models cannot account for teleconnection patterns of 2-m surface air temperature ( T 2m ) in North America, where very low temperatures in the contiguous United States may coincide with very high temperatures in Alaska in the wintertime. This dependence between warm and cold extremes motivates the need for a model with opposite-tail dependence in spatial extremes. This work develops a statistical modeling framework that has flexible behavior in all four pairings of high and low extremes at pairs of locations. In particular, we use a mixture of rotations of common Archimedean copulas to capture various combinations of four-corner tail dependence. We study teleconnected T 2m extremes using ERA5 of daily average 2-m temperature during the boreal winter. Further, the estimated mixture model quantifies the strength of opposite-tail dependence between warm temperatures in Alaska and cold temperatures in the midlatitudes of North America, as well as the reverse pattern. These dependence patterns are shown to correspond to blocked and zonal patterns of midtropospheric flow. This analysis extends the classical notion of correlation-based teleconnections to considering dependence in higher quantiles.

54 ENVIRONMENTAL SCIENCES↗

A new method for diagnosing effective radiative forcing from aerosol–cloud interactions in climate models

Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.

54 ENVIRONMENTAL SCIENCES↗

Valuing Resilience for Microgrids: Challenges, Innovative Approaches, and State Needs

The United States depends on the delivery of reliable, affordable, clean, and safe electricity. Electric utilities invest billions of dollars each year in generation, transmission, and distribution assets to meet this need. However, experiences with recent natural disasters of increasing frequency and duration demonstrate the shortcomings of this approach in the face of modern threats. Further, as customers rely on electricity for a broader range of important needs, such as transportation, as well as critical life-saving services and mission critical facilities such as water treatment, medical care, shelters, telecommunications, and more, the need to minimize the likelihood and impacts of outages grows. Against this backdrop, resilience has emerged as a key consideration to guide electricity spending, whether from utilities, customers, or taxpayers. Although reliability has been defined and measured for decades with broadly accepted metrics that measure how many customers lose power and at what frequency and duration, resilience considers the electricity system’s response to a disruption and its subsequent impacts on customers. Developing tools and methods to accurately assess the costs and benefits of resilience investments is a critical step toward the goal of mitigating the impacts of outages on customers and society. Today, electric system resilience is largely treated as an externality due to challenges estimating the costs of long-duration outages, impacts of outages on society, and increasing reliance on electricity for a growing set of interdependent services. These interdependencies include the water, wastewater, telecommunications, natural gas, and health sectors. Without knowing how much a given resilience investment will benefit customers or society more broadly, investors, policymakers, and regulators are less likely to make or approve such investments, and less able to prioritize those investments. State Energy Offices and public utility commissions (PUCs) lead the development of state-level energy policy and utility regulation, respectively, and each has interests in encouraging appropriate public and private investments in resilience. To this end, the National Association of State Energy Officials (NASEO) and National Association of Regulatory Utility Commissioners (NARUC), with the support of the U.S. Department of Energy (DOE) Office of Electricity (OE), formed a joint Microgrids State Working Group to explore the costs and benefits of microgrids, barriers to broader deployment of microgrids to meet resilience and other objectives, and policy and regulatory strategies to optimize investments in resilience, including but not limited to microgrids. Although no universally accepted valuation tool for resilience exists, National Laboratories, utilities, researchers, and state and federal agencies have collaborated to develop, apply, and improve a number of approaches to quantify resilience, several of which are still in progress at the time this report is published. This report seeks to share these important advances by discussing current definitions of resilience (Section 1), how microgrids are defined and used to meet resilience objectives (Section 2), new approaches to valuing resilience (Section 3), steps State Energy Offices and PUCs have taken to further resilience valuation efforts (Section 4), and finally, considerations and suggested next steps for State Energy Offices and PUCs (Section 5). Relevant examples of specific microgrid projects and resilience valuation efforts are included throughout the report. While this report is written specifically for NASEO and NARUC members, it may be useful for utilities, local governments, and individual customers interested in improving the way public and private dollars are spent to achieve resilience outcomes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗

Draft genome of the switchgrass head smut pathogen Tilletia maclaganii

Tilletia maclaganii is a smut fungal pathogen that causes significant biomass reduction of switchgrass ( Panicum virgatum ) used for animal forage and biofuel production. Here we present the annotated genome of T. maclaganii , strain Tm001-NY21, estimated at 42.79 Mb in size, in 53 assembled contigs and encoding 10,235 predicted genes. This genome will be important for future comparative studies of Ustilaginales across its geographic and host range.

PacBio↗

Energy Northwest - Horn Rapids Solar and Storage: An Assessment of Battery Technical Performance

Chartered in 1957 as a joint action agency of the state, Energy Northwest (ENW) is a consortium of 27 public utility districts and municipalities across Washington state. ENW takes advantage of economies of scale and shared services to help utilities run their operations more efficiently and at lower cost, to the benefit of more than 1.5 million customers. ENW develops, owns, and operates a diverse mix of electricity generating resources, including hydro, solar, and wind projects – and the Northwest’s only active nuclear energy facility. These projects provide enough reliable, affordable, and environmentally responsible energy to power more than a million homes each year, and that carbon-free electricity is provided at the cost of generation. The agency continually explores new generation projects to meet its members’ needs. In 2017, as part of the second round of funding from the Washington state Clean Energy Fund, the Washington State Department of Commerce granted up to $3 million in matching funds to develop an estimated $6.5 million project that deployed a 4-MW, 20-acre solar generating array of photovoltaic (PV) panels coupled with a 1 MW/5.5 MWh lithium-iron-phosphate battery energy storage system (BESS) in Richland, Washington. The combination of PV and BESS will provide a predictable, renewable generating source and will also serve as a training ground for solar and battery technicians throughout the nation. The City of Richland will purchase the power from the project and utilize the benefits of the energy storage. The project provides Washington state with its first opportunity to integrate a large-scale solar and storage facility into its clean mix of hydro, nuclear, and wind resources. This first-of-its-kind facility combines solar generation with battery storage and technician training. In 2019, Pacific Northwest National Laboratory (PNNL) worked with ENW to assess the integrated PV and BESS in representative use cases that could benefit the City of Richland. Between March and May 2022, extensive testing was conducted, and the results were used to assess the technical performance of the BESS subjected to actual field operations. Both reference performance and use case tests were performed: (A) Reference performance tests assess the general technical capabilities of the BESS, such as energy capacity, round-trip efficiency (RTE), ramp rate, and signal tracking capability. These are the first tests performed (baseline), and they are repeated after use case tests (post cycle). A standardized U.S. Department of Energy (DOE) energy storage performance protocol was used to characterize the BESS, including representative duty cycle profiles, test procedure guidance, and calculation guidance for determining key characteristics. (B) Use case tests examine the performance of the BESS for specific use cases using duty cycles developed by PNNL in collaboration with ENW. Five use cases were selected for testing: 1) demand charge reduction, 2) load shaping, 3) transmission charge reduction, 4) Volt-VAR service, and 5) outage mitigation. The use case duty cycles were developed based on utility and site-specific characteristics in addition to the technical characteristics and physical capabilities of the BESS. Use case tests were performed between the baseline and post cycle tests. This report describes the BESS and its components, presents testing and performance analysis results, and shares key insights and lessons learned from this project. Outcomes of the tests and analyses will help ENW understand the performance of the Horn Rapids BESS in its current state and design appropriate operational strategies for this and other BESSs over the long term.

14 SOLAR ENERGY↗

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ~5-50 cm -1 . This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effect of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal/mol -1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inverse rendering of fusion plasmas: Inferring plasma composition from imaging systems

In this work, we develop a differentiable rendering pipeline for visualising plasma emission within tokamaks, and estimating the gradients of the emission and estimating other physical quantities. Unlike prior work, we are able to leverage arbitrary representations of plasma quantities and easily incorporate them into a non-linear optimisation framework. The efficiency of our method enables not only estimation of a physically plausible image of plasma, but also recovery of the neutral Deuterium distribution from imaging and midplane measurements alone. We demonstrate our method with three different levels of complexity showing first that a poloidal neutrals density distribution can be recovered from imaging alone, second that the distributions of neutral Deuterium, electron density and electron temperature can be recovered jointly, and finally, that this can be done in the presence of realistic imaging systems that incorporate sensor cropping and quantisation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗

Tree-level carbon stock estimations across diverse species using multi-source remote sensing integration

Forests are critical carbon sinks, and remote sensing has been increasingly widely used for forest monitoring and biomass estimations. However, species-specific tree-level studies remain limited. In this study, we demonstrated the feasibility of integrating UAV-based LiDAR with high-resolution optical satellite imagery (0.5 m) to estimate biomass for individual trees across different species. The proposed method accurately estimated biomass for 53 trees (R² = 0.82, rRMSE = 0.44), with species-specific datasets, showing an average 25.2% increase in R² and a 14.8% reduction in rRMSE. A novel vegetation index combining forest structure parameters with vegetation indices (VIs) was developed using high-resolution multispectral satellite data (3 m) to explore its relationship with individual tree biomass. Combining forest structural parameters with VIs further improved estimation accuracy, achieving an R²of 0.89 and an rRMSE of 0.34. Species-specific datasets show an 11.6% increase in R²compared to methods without VIs, and a 22.2% improvement over methods using only VIs. SHapley Additive exPlanations (SHAP) analysis shows that the volume feature played a key role in model performance and remained stable throughout the training process. Altogether, the proposed approach enhances individual tree biomass and carbon sink estimations, showing great potential for large-scale precise forest carbon monitoring using multi-source remote sensing data.

59 BASIC BIOLOGICAL SCIENCES↗

Constraints on Local Primordial Non-Gaussianity with 3D Velocity Reconstruction from the Kinetic Sunyaev-Zeldovich Effect

The cosmic velocity field is an unbiased probe of the total matter distribution but is challenging to measure directly at intermediate and high redshifts. The large-scale velocity field imprints a signal in the cosmic microwave background (CMB) through the kinetic Sunyaev-Zeldovich (kSZ) effect. We perform the first 3D reconstruction of the large-scale velocity field from the kSZ effect by applying a quadratic estimator to CMB temperature maps and the 3D positions of galaxies. We do so by combining CMB data from the fifth data release of the Atacama Cosmology Telescope (in combination with Planck) and a spectroscopic galaxy sample from the Sloan Digital Sky Survey. We then measure the galaxy-velocity cross-power spectrum and detect the presence of the kSZ signal at a signal-to-noise ratio of 7.2⁢𝜎. Using this galaxy-velocity cross-correlation alone, we constrain the amplitude of local primordial non-Gaussianity finding 𝑓 NL =−9⁢0$^{+210}_{−350}$. In conclusion, this pathfinder measurement sets the stage for joint galaxy-CMB kSZ constraints to significantly enhance the 𝑓 NL information obtained from galaxy surveys through sample variance cancellation.

79 ASTRONOMY AND ASTROPHYSICS↗

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Source Characterization at Low Yields from Physics Experiment One A

The detonation of a well-recorded coupled explosion at the well-characterized Nevada National Security Site, as part of Physics Experiment One, has provided us with an opportunity to compare and contrast methods and results of yield estimation and other event characterization methods. Yield estimates are made using several methods, including local/regional magnitudes, waveform envelopes, acoustic amplitudes, joint seismoacoustic yield, and moment tensors. We include in our analysis a discussion and comparison of applicability, calibrations, assumptions, and uncertainties involved in each. With the exception of magnitude-based methods, all the techniques were more or less successful in recovering the known yield. It is, however, useful to be cognizant of requirements and calibrations, such as station coverage, emplacement conditions, propagation characterization, and uncertainties, to apply the methods to other areas with less complete information. We demonstrate that we can be successful at monitoring low-yield explosions using a variety of methods assuming that the signal-to-noise ratio is adequate at recording stations, and that the material properties are taken into account. Depth estimates are also important, particularly if one allows the possibility of a surface explosion.

Pasyanos, Michael E. [Lawrence Livermore National ↗

Disaggregating Future Retail Electricity Rate Growth [Slides]

Recent Berkeley Lab research found that modest retail rate increases over the past 10 years were mostly driven by large increases in capital expenditures (CapEx) that were offset in part by substantial wholesale price reductions. Decision-makers are increasingly concerned about the potential future rate impacts of a number of policies and industry trends that support rapid decarbonization, electrification, and grid modernization. Using historical FERC Form 1 data and the existing literature on policies and industry trends that are likely to affect utility-incurred costs and retail sales, Berkeley Lab researchers developed ranges of forecasted growth rates for cost-related rate drivers (i.e., fuel and purchased power; transmission, distribution, generation, and other categories of both non-fuel operations & maintenance and CapEx) and non-cost related rate drivers (i.e., retail sales, peak demand, and customers). These were then used as inputs to a pro-forma utility financial model (FINDER) that estimated the growth in retail electric rates between 2020 and 2030 for a prototypical vertically-integrated investor-owned utility in the United States. The analysis produced the following results: 1. Assuming average growth rates in all rate drivers, future retail rate growth is driven by sizable increases in all CapEx costs, where fuel and purchased power costs are replaced by generation CapEx as the largest rate component between 2020 and 2030. 2. Growth in sales/peak demand/customers, generation CapEx costs, and fuel and purchased power (FPP) costs, in isolation, produce the most uncertainty in rate growth. Specifically, a 1% increase in the compound annual growth rate (CAGR) of retail sales, coincident peak demand (CP), and customers (Sales-CP-Cust) results in a 0.88-0.93% decrease in the CAGR of rates, in isolation. However, a 1% increase in the CAGR of generation CapEx budgets results in a 0.07-0.14% increase in the CAGR of rates, while a 1% increase in the CAGR of FPP costs causes a 0.10-0.14% increase in the CAGR of rates, all else being equal. 3. Taking into account the correlation and variability of the growth in all rate drivers jointly, generation CapEx is expected to be both the largest and most uncertain rate component by 2030 (20-25% share of the retail rate). Transmission and distribution CapEx, along with fuel and purchased power costs are each expected to comprise between 12% and 17% of retail rates.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗

Efficient Active Learning for Gaussian Process Classification by Error Reduction

Active learning sequentially selects the best instance for labeling by optimizing an acquisition function to enhance data/label efficiency. The selection can be either from a discrete instance set (pool-based scenario) or a continuous instance space (query synthesis scenario). In this work, we study both active learning scenarios for Gaussian Process Classification (GPC). The existing active learning strategies that maximize the Estimated Error Reduction (EER) aim at reducing the classification error after training with the new acquired instance in a onestep-look-ahead manner. The computation of EER-based acquisition functions is typically prohibitive as it requires retraining the GPC with every new query. Moreover, as the EER is not smooth, it can not be combined with gradient-based optimization techniques to efficiently explore the continuous instance space for query synthesis. To overcome these critical limitations, we develop computationally efficient algorithms for EER-based active learning with GPC. Further, we derive the joint predictive distribution of label pairs as a one-dimensional integral, as a result of which the computation of the acquisition function avoids retraining the GPC for each query, remarkably reducing the computational overhead. We also derive the gradient chain rule to efficiently calculate the gradient of the acquisition function, which leads to the first query synthesis active learning algorithm implementing EER-based strategies. Our experiments clearly demonstrate the computational efficiency of the proposed algorithms. We also benchmark our algorithms on both synthetic and real-world datasets, which show superior performance in terms of sampling efficiency compared to the existing state-of-the-art algorithms.

97 MATHEMATICS AND COMPUTING↗