Search NASA⌕ Search

SEARCH · Search NASA

Results for “commercial technology regional development program”

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.

84 records · Page 5

Quantification of Energy Savings and Demand Reduction for a Heat Pump Integrated with Thermal Energy Storage (Final Report)

To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries or installing new peaking power plants. Stakeholders and policy makers across the United States have expressed interests in promoting TES, as demonstrated by the US Department of Energy’s Grid-Interactive Efficient Buildings program and the efforts of various state legislatures. However, the cost value provided by TES are unclear. If reliable cost benefits were determined, stakeholders would have a clearer picture of the financial returns that can be gained from their investment in TES. The study in this report is conducted by ORNL with collaboration with Emerson the Helix Innovation Center. In the first part of this report, EnergyPlus was used to perform whole-building simulations for two residential buildings in Indianapolis and Atlanta. The HVAC system in both buildings were equipped with phase change material TES. The TES tank was charged in off-peak hours and discharged in peak hours to perform load shifting. First, the economic value implied by existing time-of-use (TOU) rates offered by utility companies was analyzed via whole-building simulation. Second, existing demand reduction (DR) incentives sourced from 3 different electrical grid administrators (i.e., California, Texas, and New England region) were surveyed to determine their implied value. The study suggests that the traditional value analysis that focuses on ROI for the building owner significantly undervalues TES technology making economic viability difficult. A more comprehensive value analysis that includes peak demand management and deferred capital for peaking power plants shows that TES should be economically viable but here the value is greater for the utility and requires large market penetration and aggregation to fully realize the benefits. Therefore, to facilitate commercialization, new business models are needed that include a broader range of stakeholders and distribute the value of TES proportionally. In the 2 nd part of this project, the benefits of a novel phase change material (PCM) integrated heat pump configuration were evaluated via detailed component based simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation. To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Flexible Stellarator Physics Facility

We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). The FESAC Long Range Plan recognized the quasi-symmetric stellarator as “the leading US approach to developing disruption-free, low-recirculating-power fusion configurations”. To deliver the ambitious Decadal Vision for Commercial Fusion Energy, we must establish a persuasive stellarator program in parallel to the tokamak one: the stellarator will prove to be a better path to a reactor if theoretical predictions are confirmed and novel optimization techniques and strategies work as desired. Since the release of the Long Range Plan in 2020, stellarators have arguably made the most significant advances of all fusion concepts. Groundbreaking results from W7-X demonstrated low neoclassical transport and the successful operation of the island divertor. Advances in theory and modeling now allow us to minimize turbulent transport, to achieve equilibria with precise quasisymmetry, to reduce neoclassical transport and fast ion loss to levels far below what has been previously achieved, and to minimize the effect of coil manufacturing errors. If realized, these advances will lead to cost-effective stellarator designs with confinement comparable to tokamaks but without the fundamental challenges of disruptions and current drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Use of NASA Satellite Data in Aiding Mississippi Barrier Island Restoration Projects

This presentation discusses a NASA Stennis Space Center project in which NASA-supported satellite and aerial data is being used to aid state and federal agencies in restoring the Mississippi barrier islands. Led by the Applied Science and Technology Project Office (ASTPO), this project will produce geospatial information products from multiple NASA-supported data sources, including Landsat, ASTER, and MODIS satellite data as well as ATLAS multispectral, CAMS multispectral, AVIRIS hyperspectral, EAARL, and other aerial data. Project objectives include the development and testing of a regional sediment transport model and the monitoring of barrier island restoration efforts through remote sensing. Barrier islands provide invaluable benefits to the State of Mississippi, including buffering the mainland from storm surge impacts, providing habitats for valuable wildlife and fisheries habitat, offering accessible recreational opportunities, and preserving natural environments for educating the public about coastal ecosystems and cultural resources. Unfortunately, these highly valued natural areas are prone to damage from hurricanes. For example, Hurricane Camille in 1969 split Ship Island into East and West Ship Island. Hurricane Georges in 1998 caused additional land loss for the two Ship Islands. More recently, Hurricanes Ivan, Katrina, Rita, Gustav, and Ike impacted the Mississippi barrier islands. In particular, Hurricane Katrina caused major damage to island vegetation and landforms, killing island forest overstories, overwashing entire islands, and causing widespread erosion. In response, multiple state and federal agencies are working to restore damaged components of these barrier islands. Much of this work is being implemented through federally funded Coastal Impact Assessment and Mississippi Coastal Improvement programs. One restoration component involves the reestablishment of the island footprints to that in 1969. Our project will employ NASA remote sensing data and products to support these federally funded efforts on multiple fronts. Landsat and ASTER data is being analyzed to assess changes in barrier island land cover over the last 35 years. ASTER, SRTM, and EAARL terrain products and other NASA airborne imagery are being applied in assessing changes in barrier island geomorphology and geospatial extent. MODIS data is being examined as a tool for sediment transport modeling by supplying geospatial data that quantifies in-water sediment concentrations. MODIS satellite data is being assessed for monitoring changes in the spatial extent of individual barrier islands. Results thus far indicate that NASA data products are useful in assessing barrier island conditions and changes. This value is enhanced with additional historical geospatial data, commercial high resolution satellite data, other non-NASA aerial imagery, and field survey data. The project s products are relevant to the Gulf of Mexico Alliance priority issues, including coastal habitat conservation, restoration and coastal community resilience. Such products will be available to state and federal agencies involved with coastal restoration. Potential end-users of these products include the National Park Service, U.S. Geological Survey, U.S. Army Corps of Engineers, Environmental Protection Agency, Mississippi Department of Environmental Quality, and Mississippi Department of Marine Resources.

Giardino, Marco↗

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection Description This dataset contains input and output data for the manuscript Mongird, K. et al. (under review) titled "Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection". Input data corresponds to gridded spatial siting attributes that are necessary to conduct a random forest machine learning analysis of siting feature importance. Output data includes SHAP feature analysis outputs, and classification report values. For data on power plant siting results referred to in the manuscript, please refer to the CERF: IM3 Projected Western US Power Plant Locations data download page. The downloadable data includes values for eight different future scenarios for the Western US. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Technical Information The dataset includes two sets of data files: (1) CERF gridded siting parameters and (2) Feature analysis outputs and classification reports. All downloadable data is in csv file format. Files with x/y coordinate information use the Albers Equal Area Conic projection (ESRI:102003). 1. CERF Gridded Siting Parameters This directory provides a balanced sample of gridded CERF siting parameters data for eight different scenarios for the Western US through 2055, seven different technologies, and eight timesteps. This data serves as input to the feature analysis. It contains the following parameters. region_name - name of region (i.e., state) sited - binary value representing whether the grid cell received a siting of that technology type (1=True) rcp - binary value representing scenario resource concentration pathway (0 = RCP4.5, 1 = RCP8.5) ssp - binary value representing scenario shared socioeconomic pathway (0 = SSP3, 1 = SSP5) climate - binary value representing cooler (0) or hotter (1) GCM forcing tech_name - generation technology name sited_year - year that values correspond to transmission_cost - cost of transmission interconnection pipeline_cost - cost of natural gas pipeline interconnection interconnection_cost - total interconnection cost (sum of transmission cost and gas pipeline cost) lmp - associated locational marginal value ($/MWh) associated with the grid cell, timestep, scenario, and technology xcoord - x-coordinate of location ycoord - y-coordinate of location 2a. Feature Analysis Output The dataset includes the feature analysis shap output for locational marginal price and interconnection cost. It contains the following parameters. technology - generator technology name scenario - name of scenario feature - name of feature, either locational_marginal_price or interconnection_cost value - the mean of absolute value of SHAP values for given feature 2b. Feature Analysis Classification Report This download includes the classification report associated with each random forest model. The dataset contains the following parameters. technology - generation technology name scenario - name of scenario test - one of precision (the proportion of predicted positives that are actually correct), recall (the proportion of actual positives that were correctly identified), f1-score (the harmonic mean of precision and recall) 0.0 - value of test for classification of 0 (grid cell not chosen for siting) 1.0 - value of test for classification of 1 (grid cell chosen for siting) accuracy - accuracy of model (i.e., fraction of all predictions that were right) macro avg - Simple average of test values for all classes weighted avg - Weighted average of test values for all classes, weighted based on Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program. License This data is made available under a CCBY4 License Disclaimer This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor the Contractor, nor any or their employees, nor any jurisdiction or organization that has cooperated in the development of these materials, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness or any information, apparatus, product, software, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORYoperated byBATTELLEfor theUNITED STATES DEPARTMENT OF ENERGYunder Contract DE-AC05-76RL01830

Mongird, Kendall [Pacific Northwest National Labor↗

Building an EPA Class VI Permit Application

Summary To accelerate the commercialization of carbon capture and storage (CCS), the US Department of Energy (US DOE) is building on decades of characterization efforts and pilot-scale projects through their CarbonSAFE program. Administered through their National Energy Technology Laboratory, this program seeks to bring fully integrated projects to the sector that can store more than 50 million tonnes of CO2 over a 30-year period. The program, which was enacted before the enhancement of Internal Revenue Code Section 45Q, is in the capture assessment, characterization, and permitting phase. The objectives of this paper are to discuss (a) the injection permitting requirements of the CarbonSAFE projects; (b) information gathering in support of the permit; (c) the timelines of field development and permit-related activities; (d) the major technical components of the field development plan; and (e) early feedback from the regulators toward acceptance of the permit. In Mississippi, more than 30,000 acres have been characterized by six deep characterization wells, a deep groundwater well, and 92 line miles of 2D seismic as part of the CarbonSAFE Project ECO2S. During the acquisition of seismic data, all receiver lines were live, which resulted in the generation of a pseudo-3D seismic design. The incorporation of a 3D seismic survey was not included as part of this project due to logistical difficulties presented by the undulating, wooded surface terrain. A suite of openhole geophysical logs was taken from each well, allowing for a detailed interpretation of prospective storage reservoirs and confining intervals to complement the analysis carried out on the 290 ft of a whole core that was cut through the prospective confining zone and storage reservoir. The detailed geologic and reservoir data were assembled and entered into a 3D model to assess the injection capacity and the area of review (AoR). This information fed into the detailed corrective action, monitoring, testing, and postinjection site care (PISC) modeling. The results have been exceptional. The geologic assessment has revealed three primary storage targets, ranging in depth from 3,500 ft to 6,000 ft. These storage reservoirs net 1,300 ft of sandstone, with mean porosity and permeability of 29% and 3.6 darcies, respectively. Together, these reservoirs have storage capacities that may exceed 20 million tonnes per square mile, making this a gigatonne prospect. Forward modeling of the project resulted in an AoR of 16 sq miles, injecting about 8000 t/d, for 30 years, via two deep injection wells. The excellent confining characteristics of the caprock, relatively simple geologic structure, and lack of historical well drilling activity in this area provide excellent containment of the injected CO2. Based on this work, the project has proposed 20 years of PISC. To date, only two US CO2 injection permits have been granted. These projects relied on a singular capture point feeding a singular sequestration point (source to sink), and considerations have not been made to garner CO2 emissions from other industrial sources. The Kemper County Storage Complex is a first-of-its-kind storage hub concept that looks to develop an area capable of storing significant quantities of CO2 from the region. Also, this work will show how characterization efforts, geological and numerical modeling efforts, and plan development were constructed in support of permit and incentives acceptance.

Energy & Fuels↗

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Residential Building Energy Efficiency Field Studies: Low-Rise Multifamily

In recent years, the U.S. Department of Energy (DOE) has conducted a series of research studies to validate energy efficient building technologies in the field. Much of the work has focused on single-family construction, and some has also addressed commercial energy codes. The work detailed in this DOE-funded study (EE0007616) focuses on low-rise multifamily buildings (three stories or fewer above grade) in various regions of the United States, and reports on how state-level building codes are being implemented, both in terms of observed characteristics and also in terms of estimated energy impacts. Nearly 100 buildings across four states—Illinois, Minnesota, Oregon, and Washington—were sampled, which represent a range of climate types from mild temperature to very cold continental. Both common entry and outdoor entry buildings were included, and a parallel research project evaluated envelope air tightness and current still-evolving air tightness testing methods. Finally, a set of structured interviews of building designers and other relevant professionals was carried to out to gain more insight into this market. To the greatest extent possible, the methodology developed under the project for low-rise multifamily buildings mirrored the approach established by Pacific Northwest National Laboratory (PNNL) for single-family residential buildings (https://www.energy.gov/eere/buildings/downloads/residential-building-energy-code-field-study). This included the general approach to sampling, recruitment, and data collection, as well as data analysis and presentation. The range of permitting dates for the sites encompassed two energy code cycles in most regions. All states in the study had adopted a variation of the International Energy Conservation Code (IECC) for the structure of their state code. The low-rise multifamily occupancy presents a hybrid building type: most of the building’s conditioned floor area was covered by the residential chapter of the code while portions of the building (such as corridors and common spaces) fell under the commercial code chapter. The key items assessed in this work were: Building Shell—exterior wall insulation, ceiling insulation, foundation insulation, windows. Common Areas—HVAC and lighting. Living Units—lighting, ventilation. A few items were not assessed in detail, given their relative paucity in this occupancy type; these included duct leakage, pipe insulation, and hot water circulation controls. Building characteristics were collected via a combination of architectural, mechanical, electrical, and plumbing plan reviews and field inspections, and entered into a spreadsheet-based tool that was later queried to build a database. Data went through quality control both upon arrival and via a later semi-automated review and assurance process. Most of the data are presented graphically so that the reader can quickly assess compliance with the applicable energy codes (both by state and by code year). As a final step, EnergyPlus™ simulations were created for all buildings in the study to estimate both the as-found energy use intensity (EUI) and the energy and CO 2 that could be saved if features that were found to not meet code minimums were brought up to code. The savings estimates were tabulated for each of the four states in the study. The research team found that the single-family approach was largely applicable to low-rise multifamily buildings. This applies to both the data collection and the prototype EUI analysis. Most of the occupied space is living units and falls under residential energy codes, and many characteristics use similar envelope construction and relatively straightforward mechanical systems and lighting. One of the most challenging aspects of this work was to build an effective spreadsheet-based data collection instrument that could allow efficient collection of both building plan and field data. The research team is of the view that other methods could be equally effective if the work is done carefully with diligent quality control. The primary findings for the work center around the thermal envelope and mechanical systems and lighting at the sites: For thermal envelope components, the majority of buildings met or were better than the prescriptive code.This suggests that building designers and builders are aware of code requirements. In some cases, surveyed buildings were designed to qualify for energy efficiency certification programs. These buildings made up at least 20% of sampled buildings in each state. Almost all buildings met mechanical system efficiency requirements (for both living units and common areas). In some cases, sites employed systems that were considerably more efficient than required by the applicable energy code. Dwelling units had a majority of high-efficacy lighting, often in excess of the state’s residential code requirements. While high-efficacy fixtures were also typical in common areas (corridors and stairwells), lighting power densities (LPDs) in these areas were sometimes higher than levels dictated by the applicable part of the state commercial energy code. The simulation models run on a series of low-rise multifamily prototypes, informed by a composite of the field data collected, calculated annual EUIs of between 20 and 50 kBtu/ft2-yr, with the range representing the effects of both building characteristics and building location (climate zone). A detailed process (based on simulations of prototype buildings) was used to estimate the amount of avoided energy use that would occur if 100% adherence to energy codes were attained. The results indicated modest savings are attainable for items such as window thermal performance and common area lighting. The result is overall only a modest potential for additional energy savings, averaging about 10% of EUI.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Thermocouple Testing in Support of the AGR-5/6/7 Experiment

This report documents thermocouple testing performed in IRC Lab C-15 over a period of seven years. This testing supported selection and characterization of the thermocouple set used in the AGR-5/6/7 experiment. The following summary was taken directly from the report. Temperature measurement is a challenging aspect of very high temperature irradiation experiments because commonly used high-temperature commercial thermocouples such as platinum-rhodium (Types S, R, and B) and tungsten-rhenium (Type C), suffer dramatic drift because of neutron-induced transmutation. As a result, these types of thermocouples, which are used routinely for industrial temperature measurements outside of reactors, are used only in very special circumstances for reactor experiments. Conversely, because of their low neutron cross-sections, Type N thermocouples are affected to only a limited extent by neutron irradiation. However, the use of these nickel-based thermocouples is limited when the temperature exceeds 1050°C due to drift arising from minor alloying elements migrating from the thermocouple's metal sheath to the thermoelements. This change in the composition of the thermo-elements results in significant decalibration of the signal. The issues described above were recognized during the early planning stages of the final AGR experiment (designated AGR-5/6/7), and a thermocouple furnace testing program was performed over a seven-year period (2014-2019, 2021) to first select and then characterize the best thermocouple set for the high temperature regions of the AGR-5/6/7 irradiation experiment. The calculated temperature range of the AGR-5/6/7 experiment was 600–1500°C. For temperatures below 1000°C standard Type N thermocouples were deemed adequate. The furnace testing campaign identified two thermocouple types suitable for measuring temperatures above 1000°C, a Mo/Nb thermocouple developed at INL called HTIR-TC, and a Type N thermocouple developed by Cambridge University (called herein Cambridge Type N), which featured a custom high nickel alloy sheath. One of the original goals of the furnace testing program was to identify a thermocouple capable of low drift operation near the peak temperature expected in AGR-5/6/7, i.e., about 1400°C. The HTIR-TC design appeared promising in this regard, however a manufacturing difficulty proved to be a barrier and instead the furnace testing focused on drift performance at 1250°C. The manufacturing difficulty was that the Nb sheaths of the HTIR-TCs experienced extreme embrittlement when heat treated at 1600°C or greater. Heat treatment is needed to stabilize the emf output of this TC type, and the higher the heat treatment temperature the higher the peak temperature of stable operation. Because of the sheath embrittlement the heat treatment temperature had to be lowered to 1450°C resulting in a stable operating temperature of about 1250°C. One of the successes of the furnace testing program was identification of a shortcoming in the heat treatment procedure that had been traditionally used in the production of HTIR-TCs. The shortcoming was that the entire heated length of the HTIR-TC sensor was not being heat treated, but rather only the part of the sensor expected to experience temperatures above 1000°C. The problem manifested itself when the thermocouples were removed from the heat treat furnace and placed in another furnace with a different geometry, their indicated temperatures would be widely scattered, but mostly in the negative direction. The solution was to heat treat the entire heated length of the sensor. Since the deepest immersion depth in the AGR-5/6/7 experiment was about 40 inches, a heat treatment length of 48 inches was used. After this change was implemented, thermocouples which were moved into a new environment with a different temperature profile (i.e., a different furnace), produced accurate temperature measurements. Although assembly of the AGR-5/6/7 experiment was completed in September of 2017 (and irradiation begun in 2018), furnace testing of thermocouples continued in 2018 and 2019. The main purpose of this testing was to establish very long-term drift characteristics of the HTIR and Cambridge Type N thermocouples installed in the experiment. Representative thermocouples from the same lots as those installed in the AGR-5/6/7 experiment were used. Additionally, thermocouples of different designs, (particularly variations on the HTIR-TC design) were "piggy-backed" on this testing program to provide insights for instrumenting future very high temperature irradiation experiments. This two-year testing program demonstrated that HTIR-TCs and Cambridge Type N TCs could operate at 1250°C for up to 10,000 hrs (and in some cases longer) while experiencing negative drifts on the order of 2-4°C/1000 hrs. This performance was considered acceptable given the extreme operating environment the sensors faced.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Austin Sustainable and Holistic Integration of Energy Storage and Solar PV [Austin SHINES]. Final Report, Version 2

The Austin SHINES project and solution is a software management platform, for an electric grid with a high penetration of dispersed photovoltaic (PV) solar generation sites, which maintains the traditional power quality and reliability associated with grid service. This project developed and deployed the platform as a Distributed Energy Resource Management System (DERMS), engaging multiple advanced controls, to evaluate operation and optimization of a fleet of diverse DER assets, installed at several locations among Austin Energy’s customers and distribution system. The project also produced a methodology to create a replicable DERMS template, adaptable to other regions and market structures. Last, Austin SHINES aimed to demonstrate the solution’s methodology would enable the DER grid ecosystem to serve load at a technical cost (System Levelized Cost of Electricity, or System LCOE) of less than the U.S. Department of Energy SHINES program metric of $0.14/kWh, in a defined boundary, while enabling a high penetration of distributed PV. Research was categorized in 6 reports (Final Deliverables = FD) listed below, with titles and descriptions indicating which area of understanding was investigated: FD-1: System Levelized Cost of Electricity (System LCOE) Methodology The creation and use of the System LCOE to Serve Load metric that encompasses the holistic, system-level costs and benefits of all resources, and enables them to be evaluated based on their ability to support an efficient and low-cost integrated grid ecosystem. FD-2: Software Platform Product Description The creation of new DER control methodologies deployable within a utility-grade software platform that enable DER's to maximize their benefit within a grid, that is capable of serving load enabling a high penetration of distributed PV generation. FD-3: Optimal Design Methodology Optimal design methodologies for individual DER installations that enable utilities to determine the optimal combinations and sizing for individual DER sites. FD-4: Austin SHINES Ownership and Operation Models for DER System Performance A comparison of multiple DER aggregation and ownership methodologies including direct utility control, third-party aggregator, and autonomous. FD-5: Economic Modeling & Optimization A comparison of multiple DER technology mixes and configurations within the distribution system, providing insight into an optimal blend of technologies that best enable the distribution system to serve load at the lowest cost at high penetrations of solar. FD-6: Fielded Assets Deployed DER assets within the Austin Energy SHINES circuits. Austin SHINES provided an opening for state-of-the-art technology products to be deployed, providing a rich opportunity for improving how each of the products perform as stand-alone products, and in concert with other complementary products. The Austin SHINES project comprised of two key metrics for System LCOE: SystemLCOE_SHINES<$0.14/kWh Modeled ΔSystemLCOE_SHINES/ΔSystemLCOE_Base≥20% at same solar penetration The System LCOE calculation uses the costs of the utility-owned infrastructure as it exists today, the cost of the DERs that exist in the system today, and the cost of the purchase of energy from ERCOT wholesale markets over the course of the calendar year. All costs are on an annualized basis. The capital and operating costs are derived from the rate case, which produces a yearly cost. The net cost of energy and services imported to the system is integrated over the test year, as is the load served and solar penetration. The first metric was easily achieved by every scenario considered. The goal was set when the Department of Energy’s SHINES Funding Opportunity Announcement was written in 2015 and was a more difficult target at the time. Due mostly to rapidly declining costs for DERs and the significant decrease in the Electric Reliability Council of Texas (ERCOT) energy market prices, which results in lower net cost of energy purchases, the System LCOE is well below this target for all scenarios considered. A fleet of DERs can assume different mixtures, each of which serves the load at a different LCOE. The optimal mixture of DERs serves load at the smallest System LCOE. The second metric (hereinafter %delta metric) asks that the holistic DERMS controls reduce the incremental cost above the baseline of going to a high solar penetration future by at least 20% as compared to the case of a DER deployment with no sophisticated controls (autonomous). Many comparison sets were created throughout this project. Physical technology was installed for informing utility engineering and testing several types of operational control schemes, through the DERMS. The types of operational control which were compared for valuation of the System LCOE Metric were: Holistic control = using the full suite of the DERMS platform to decide and optimize how/why the systems operate depending on weather, market, and reliability signal input. Autonomous control = a local mode at the asset site, wherein a schedule operates the asset, with visibility into performance only No control = the baseline for comparing value against the other two types of control The types of ownership control included: Direct Utility control = the utility dispatches a signal to each asset Third-Party Aggregator = a third party aggregates a fleet of assets and the utility dispatches one signal for all Autonomous = a local mode is set for operation at the asset site, wherein a schedule operates the asset, with visibility into performance only The types of control methodologies deployable within a utility-grade software platform included: Utility Peak Load Reduction = Lower transmission cost obligation Day-Ahead Energy Arbitrage = Realize economic value through price differential Real-Time Price Dispatch = Realize economic value from real-time price spikes Voltage support = Reduce losses and increase solar generation Distribution Congestion Management = Increase local grid reliability Demand Charge Reduction = Lower customer bills and realize system benefit The fielded assets deployed for the project were: Utility Scale Kingsbery Energy Storage System: 1.5 MW / 3 MWh Li-Ion battery storage Mueller Energy Storage System: 1.75 MW / 3.2 MWh Li-Ion battery storage, 7 Energy Storage Units (250 kW each) La Loma Community Solar: 2.6 MW Commercial Scale Aggregated storage installations at 3 sites, with existing solar (300+ kW): One 18 kW / 36 kWh Li-Ion battery storage Two 72 kW / 144 kWh Li-Ion battery storage Residential Scale Aggregated storage installations: -Six stationary battery storage systems (10 kWh each) at homes with existing solar -One Electric Vehicle installed as Vehicle-to-Grid (V2G) Utility-Controlled Solar via Smart Inverters at 12 homes Autonomously-Controlled Smart Inverters at 6 homes Over the course of the project, Austin SHINES undertook installing more than 3 MW of distributed battery energy storage, smart PV inverters, a DER control platform, and other enabling technologies utilizing customer and utility locations and aggregation models. All of these resources were to be integrated and optimized at the utility level. DER assets and control methodologies were designed to achieve a credible pathway to a System LCOE for energy delivered to load of $0.14//kWh or less by 2020, while maximizing distributed solar generation and maintaining acceptable standards of power quality. The project also established a template for other regions to follow, to maximize the adoption of distributed solar PV in support of an economic and efficient grid. In total, the Austin SHINES project added value to the DER subject area in each layer of integration. From utility, to commercial to residential scales, the sheer hierarchy of communication and coordination was a significant accomplishment in addition to learnings from what these communications revealed was unique to each. Economically, the most effective method demonstrated was the criticality of planning phases. Contingencies and multiple projection scenarios helped guide the project to deploy optimal design as close as feasible, in real world conditions. The project and reports will serve public benefit by outlining specific areas of DER strategy and installation where many stakeholders and needs can be addressed with improved efficiency. Overall, communities and utilities should use the results to guide the increasing options available for powering the grid with DER, renewables, and carbon considerate energy.

14 SOLAR ENERGY↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗