Search NASA⌕ Search

SEARCH · Search NASA

Results for “generator parameter estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,495 records · Page 83

Validation of the Archived CERES Surface and Atmosphere Radiation Budget (SARB) at SGP

The CERES Surface and Atmosphere Radiation Budget (SARB) product (Charlock et al, 2002) includes the vertical profile of broadband SW, broadband LW, and 8-12 micron window (WN) fluxes; upwelling and downwelling at TOA, 70 hPa, 200 hPa, 500 hPa, and the surface; and for all-sky and clear-sky conditions. We test the archived CERES TRMM record of SARB for January-August 1998 and focus on discrepancies with ground-based measurements at SGP. The CERES SARB is generated by a highly modified Fu-Liou radiative transfer code (Fu and Liou, 1993). The most critical inputs for this application are cloud optical properties (fractional area, optical depth, particle size and phase, height of top, and estimate of geometrical thickness Minnis et al., 2002) from the narrowband VIRS imager. Numerous VIRS pixels (approx. 2km resolution at nadir) are matched to each of the large (approx. 20km) CERES broadband footprints (Wielicki et al, 1996). Other inputs include temperature and humidity from ECMWF (Rabier et al, 1998) , NCEP ozone profiles from SBUV and TOVS (Yang et al, 2001), aerosol optical thickness (AOT) from the Model for Atmospheric Transport and Chemistry (MATCH) aerosol assimilation (Collins et al., 2001) or alternately from the VIRS imager (Ignatov and Stowe, 2000). VIRS AOT is available for clear and partly cloudy ocean footprints during daylight; and only when viewing geometry renders a contribution from sunglint as unlikely. For other footprints, AOT is taken from MATCH. AOT is apportioned into fractions of dust (Tegan and Lacis, 1996), sea salt, sulfate, dust, soluble organic, insoluble organic, and soot (Hess et al., 1996) using the 6-hourly MATCH output. Tuned fluxes are retrieved by adjusting inputs to nudge computed TOA fluxes toward CERES observations (Rose et al, 1997). In clear conditions, the fields of humidity, surface skin temperature, surface albedo and AOT are adjusted to produce a closer match of computed and observed fluxes at TOA. When CERES footprints have clouds, the cloud optical thickness, fractional area within the footprint, and temperature of cloud top are adjusted by the tuning algorithm. Both tuned and untuned fluxes are archived, as are the respective adjustments to any parameters at the surface or within the atmosphere.

Charlock, Thomas P.↗

Virtual Sensors: Using Data Mining to Efficiently Estimate Spectra

Detecting clouds within a satellite image is essential for retrieving surface geophysical parameters, such as albedo and temperature, from optical and thermal imagery because the retrieval methods tend to be valid for clear skies only. Thus, routine satellite data processing requires reliable automated cloud detection algorithms that are applicable to many surface types. Unfortunately, cloud detection over snow and ice is difficult due to the lack of spectral contrast between clouds and snow. Snow and clouds are both highly reflective in the visible wavelen,ats and often show little contrast in the thermal Infrared. However, at 1.6 microns, the spectral signatures of snow and clouds differ enough to allow improved snow/ice/cloud discrimination. The recent Terra and Aqua Moderate Resolution Imaging Spectro-Radiometer (MODIS) sensors have a channel (channel 6) at 1.6 microns. Presently the most comprehensive, long-term information on surface albedo and temperature over snow- and ice-covered surfaces comes from the Advanced Very High Resolution Radiometer ( AVHRR) sensor that has been providing imagery since July 1981. The earlier AVHRR sensors (e.g. AVHRR/2) did not however have a channel designed for discriminating clouds from snow, such as the 1.6 micron channel available on the more recent AVHRR/3 or the MODIS sensors. In the absence of the 1.6 micron channel, the AVHRR Polar Pathfinder (APP) product performs cloud detection using a combination of time-series analysis and multispectral threshold tests based on the satellite's measuring channels to produce a cloud mask. The method has been found to work reasonably well over sea ice, but not so well over the ice sheets. Thus, improving the cloud mask in the APP dataset would be extremely helpful toward increasing the accuracy of the albedo and temperature retrievals, as well as extending the time-series of albedo and temperature retrievals from the more recent sensors to the historical ones. In this work, we use data mining methods to construct a model of MODIS channel 6 as a function of other channels that are common to both MODIS and AVHRR. The idea is to use the model to generate the equivalent of MODIS channel 6 for AVHRR as a function of the AVHRR equivalents to MODIS channels. We call this a Virtual Sensor because it predicts unmeasured spectra. The goal is to use this virtual channel 6. to yield a cloud mask superior to what is currently used in APP . Our results show that several data mining methods such as multilayer perceptrons (MLPs), ensemble methods (e.g., bagging), and kernel methods (e.g., support vector machines) generate channel 6 for unseen MODIS images with high accuracy. Because the true channel 6 is not available for AVHRR images, we qualitatively assess the virtual channel 6 for several AVHRR images.

Srivastava, Ashok↗

Identification of NPP transients using artificial intelligence - 335

If the accidents happen in nuclear power plants (NPPs), the NPP operators identify the transients through observing the behaviors of the major parameters. However, the integrity of the instrumentation signals is not ensured under the accident circumstances and the major parameters essential in identifying the transients can rapidly change although the instruments generate accurate signals. Thus, to successfully control the accidents by taking necessary measures relatively faster, it is very important to provide the operators with the information that can identify the transients and accidents accurately. In this study, the transients of severe accidents in NPPs are identified and classified using artificial intelligence methodologies. Several transients are identified using support vector classification (SVC) and probabilistic neural network (PNN) as artificial intelligence methodologies. The used data for the transients are the simulation data using the modular accident analysis program (MAAP) code for the optimized power reactors (OPR1000s). The time-integrated values of chosen simulated sensor signals are used as the input variables for identifying the transients. As a result of this study, the proposed models identified the transients quite precisely. Therefore, the quite accurate identification of the transients indicates the excellence of these artificial intelligence methods and they will produce useful supporting information for the operators even under extreme transient circumstances in NPPs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Vulnerability of OFDR-based distributed sensors to radiations

Silica-based optical fibers have recently attracted much interest for their use in harsh environments such as the ones encountered in space, military or high energy physics applications. Small size, fast response, light weight and immunity to electromagnetic fields are favorable advantages that often become decisive for fiber sensing to be chosen over other conventional sensing technologies. As an important and representative example, Fukushima's accident highlighted weaknesses in the safety of nuclear power plants. Since, one of the strategic research axis of the nuclear industry is devoted to the development of novel technologies and sensors to enhance and reinforce the safety in nuclear power plants, especially in the case of accidental conditions associated with a strong increase of the constraints applied to the fiber-based system. The objective of this research field is to develop classes of distributed fiber-based sensors using scattering-based techniques, powerful solutions for various measurands measurement. Optical fiber properties, indeed, depend on several external parameters such as temperature, strain and therefore the fiber itself can be used as the sensitive element. Different classes of fiber-based sensing techniques have been recently investigated such as Fiber Bragg Gratings (FBGs) for discrete measurements and Brillouin, Raman and Rayleigh [8,9] scattering based techniques for distributed measurements of various environmental parameters. Whereas Brillouin and Raman sensor resolutions remain in the range of one meter, the advantage of Rayleigh scattering based technique is that it offers very high spatial resolution from 1 cm down to few μm over several hundred meters of fiber length down to few meters respectively. For nuclear industry, integrating fibers-based sensors has to improve the performances (resolution, operating range,...) of security systems in current nuclear power plants (NPPs) and offers new alternative technologies that may overcome the issues identified for next generation of NPPs. Such integration will only be possible if the OFDR based systems are able to resist to the constraints associated with industrial environments, one of the most constraining being the presence of high level of radiations. In this work, we carry out a systematic study to highlight the OFDR interest and sensitivity to probe the optical samples at high irradiation dose levels. The responses of five optical fibers types, from radiation hardened to radiation sensitive ones, are investigated to explore the influence of both the material compositions and the γ-irradiation on the ODFR sensors. Using these samples, we should highlight the influence of the core dopant concentration on the observed radiation-induced changes as well as the difference observed when the cladding is either radiation resistant or radiation sensitive. Our samples were irradiated using a {sup 60}Co source facility reaching total doses varying from 1 MGy up to a maximum of 10 MGy. All the measurements are performed after diverse months from irradiation to study permanents effects induced from these high γ-rays doses. We'll present at the conference all the experimental results acquired and use them to estimate the potential of OFDR-based systems for operation in radiation environments. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of an All-Aqueous Thermally Regenerative Redox Flow Battery to Support Fossil Fuel Assets

Low-temperature thermal energy, a largely untapped resource, holds significant promise for large-scale electrical power generation globally. Various stationary sources, including industrial entities and thermal power plants, emit considerable low-temperature heat that currently remains unutilized. This energy is often overlooked because its low temperature renders it unsuitable for efficient power generation using conventional methods. However, current research is exploring diverse technologies capable of converting low-temperature heat into grid-scale power, aiming to enhance grid efficiency, further decarbonization initiatives, and facilitate a shift toward more decentralized power systems. One such innovative technology is the thermally regenerative battery (TRB), noted for its high power and energy densities compared to similar technologies, positioning it as a potential game-changer in power generation. TRBs integrate two scalable and well-established unit operations: a redox flow battery and a distillation column. This integration suggests that once an effective TRB chemistry is established, the pathway to commercialization could be expedited. The copper-based thermally regenerative ammonia battery (Cu aq -TRAB) stands out as the first TRB that circumvents electrodeposition/dissolution reactions, stabilizing Cu(I) and Cu(II) within the electrolyte and maintaining stability of all electroactive species in an aqueous phase. This stabilization has led to improvements in coulombic efficiency, open circuit potential, and copper solubility, thereby enhancing power density, energy density, and overall energy efficiency. Preliminary tests were conducted to determine the effects of various electrolyte species on the performance metrics of the battery, both theoretically and experimentally. These tests revealed that the solubility of copper in the Cu aq -TRAB electrolyte was constrained by the Cu(I)-NH 3 complex. Adjusting the background electrolyte to 5 M NH4Br and the ligand concentration to 4 M NH 3 enabled the copper concentration to reach a maximum of 0.6 M. This modification led to an estimated theoretical maximum energy density of 9.5 Wh L -1 for the Cu aq -TRAB. Additionally, full cell testing indicated a tradeoff between peak power and energy density with varying copper and ammonia concentrations. Increasing the applied current density during discharge linearly raised the average power output, with a minimal reduction in energy density due to a balance between higher ohmic overpotential and reduced time for undesirable ammonia crossover. Furthermore, a comprehensive numerical sensitivity analysis of the complete Cu aq -TRAB system was performed. This analysis aimed to assess how the battery and the distillation column responded to changes in system input parameters, providing insights into optimal research directions for enhancing system performance. The analysis revealed that at room temperature, battery power was significantly more sensitive to ohmic losses than to mass transfer, with reaction rates having minimal impact. This trend continued even at higher temperatures. Also, the thermal energy required for ammonia separation was studied, showing that increased temperatures generally reduced energy requirements, except in low-pressure scenarios above 65 °C. An investigation into membrane performance in the Cu aq -TRAB was undertaken, given the significant impact of ammonia transport control and ohmic losses on system performance. Various membranes were evaluated to identify key performance metrics. Among the tested membranes, Selemion CMVN exhibited the highest performance, with a peak power density of 84 mW cm -2 and average values of 26 ± 6.8 mW cm -2 for power density and 2.9 Wh L -1 for energy density at an applied current density of 50 mA cm -2 . An economic assessment indicated a levelized cost of storage at $410 per MWh under optimal conditions, highlighting the commercial potential of the Cu aq -TRAB when utilizing cost-effective, readily available materials.

25 ENERGY STORAGE↗

Initial Sensorimotor and Cardiovascular Data Acquired from Soyuz Landings: Establishing a Functional Performance Recovery Time Constant

INTRODUCTION Testing of crew responses following long-duration flights has not been previously possible until a minimum of more than 24 hours after landing. As a result, it has not been possible to determine the trend of the early recovery process, nor has it been possible to accurately assess the full impact of the decrements associated with long-duration flight. To overcome these limitations, both the Russian and U.S. programs have implemented joint testing at the Soyuz landing site. This International Space Station research effort has been identified as the functional Field Test, and represents data collect on NASA, Russian, European Space Agency, and Japanese Aerospace Exploration Agency crews. RESEARCH The primary goal of this research is to determine functional abilities associated with long-duration space flight crews beginning as soon after landing as possible on the day of landing (typically within 1 to 1.5 hours). This goal has both sensorimotor and cardiovascular elements. To date, a total of 15 subjects have participated in a 'pilot' version of the full 'field test'. The full version of the 'field test' will assess functional sensorimotor measurements included hand/eye coordination, standing from a seated position (sit-to-stand), walking normally without falling, measurement of dynamic visual acuity, discriminating different forces generated with the hands (both strength and ability to judge just noticeable differences of force), standing from a prone position, coordinated walking involving tandem heel-to-toe placement (tested with eyes both closed and open), walking normally while avoiding obstacles of differing heights, and determining postural ataxia while standing (measurement of quiet stance). Sensorimotor performance has been obtained using video records, and data from body worn inertial sensors. The cardiovascular portion of the investigation has measured blood pressure and heart rate during a timed stand test in conjunction with postural ataxia testing (quiet stance sway) as well as cardiovascular responses during sensorimotor testing on all of the above measures. We have also collected motion sickness data associated with each of the postflight tests. When possible rudimentary cerebellar assessment was undertaken. In addition to the immediate post-landing collection of data, postflight data has been acquired twice more within 24 hours after landing and measurements continue until sensorimotor and cardiovascular responses have returned to preflight normative values (approximately 60 days postflight). SUMMARY The level of functional deficit observed in the crew tested to date is more severe than expected, clearly triggered by the acquisition of gravity loads immediately after landing when the demands for crew intervention in response to emergency operations will be greatest. Measureable performance parameters such as ability to perform a seat egress, recover from a fall or the ability to see clearly when walking, and related physiologic data (orthostatic responses) are required to provide an evidence base for characterizing programmatic risks and the degree of variability among crewmembers for exploration missions where the crew will be unassisted after landing. Overall, these early functional and related physiologic measurements will allow the estimation of nonlinear sensorimotor and cardiovascular recovery trends that have not been previously captured.

Reschke, M. F.↗

Toward Right-Fidelity Rotorcraft Conceptual Design

The aviation Advanced Design Office (ADO) of the US Army Aeroflightdynamics Directorate (AMRDEC) performs conceptual design of advanced Vertical Takeoff and Landing (VTOL) concepts in support of the Army's development and acquisition of new aviation systems. In particular, ADO engages in system synthesis to assess the impact of new technologies and their application to satisfy emerging warfighter needs and requirements. Fundamental to ADO being successful in accomplishing its role; is the ability to evaluate a wide array of proposed air vehicle concepts, and independently synthesize new concepts to inform Army and DoD decision makers about the tradespace in which decisions will be made (Figure 1). ADO utilizes a conceptual design (CD) process in the execution of its role. Benefiting from colocation with NASA rotorcraft researchers at the Ames Research Center, ADO and NASA have engaged in a survey of the current rotorcraft PD practices and begun the process of improving those capabilities to enable effective design and development of the next generation of VTOL systems. A unique aspect of CD in ADO is the fact that actual designs developed in-house are not intended to move forward in the development process. Rather, they are used as reference points in discussions about requirements development and technology impact. The ultimate products of ADO CD efforts are technology impact assessments and specifications which guide industry design activity. The fact that both the requirement and design are variables in the tradespace adds to the complexity of the CD process. A frequent need is ability to assess the relative "cost" of variations in requirement for a diverse set of VTOL configurations. Each of these configurations may have fundamentally different response characteristics to this requirement variation, and such insight into how different requirements drive different designs is a critical insight ADO attempts to provide decision makers. The processes and tools utilized are driven by the timeline in which questions must be answered. This can range from quick "back-of-the-envelope" assessments of a configuration made in an afternoon, to more detailed tradespace explorations that can take upwards of a year to complete. A variety of spreadsheet based tools and conceptual design codes are currently in use. The in-house developed conceptual sizing code RC (Rotorcraft) has been the preferred tool of choice for CD activity for a number of years. Figure 2 illustrates the long standing coupling between RC and solid modeling tools for layout, as well as a number of ad-hoc interfaces with external analyses. RC contains a sizing routine that is built around the use of momentum theory for rotors, classic finite wing theory, a referred parameter engine model, and semi-emperical weight estimation techniques. These methods lend themselves to rapid solutions, measured in seconds and minutes. The successful use of RC, however requires careful consideration of model input parameters and judicious comparison with existing aircraft to avoid unjustified extrapolation of results. RC is in fact a legacy of a series of codes whose development started in the early 1970s, and is best suited to the study of conventional helicopters and XV-15 style tiltrotors. Other concepts have been analyzed with RC, but typically it became necessary to modify the source code and methods for each unique configuration. Recent activity has lead to the development of a new code, NASA Design and Analysis of Rotorcraft (NDARC). NDARC uses a similar level of analytical fidelity as RC, but is built on a new framework intended to improve modularity and ability to rapidly model a wider array of concepts. Critical to achieving this capability is the decomposition of the aircraft system into a series of fundamental components which can then be assembled to form a wide-array of configurations. The paper will provide an overview of NDARC and its capabilities.

Sinsay, Jeffrey D.↗

An Integrated Approach To Offshore Wind Energy Assessment: Great Lakes 3D Wind Experiment

This grant supported fundamental research into the characterization of flow parameters of relevance to the wind energy industry focused on offshore and the coastal zone. A major focus of the project was application of the latest generation of remote sensing instrumentation and also integration of measurements and numerical modeling to optimize characterization of time-evolving atmospheric flow parameters in 3-D. Our research developed a new data-constrained Wind Atlas for the Great Lakes, and developed new insights into flow parameters in heterogeneous environments. Four experiments were conducted during the project: At a large operating onshore wind farm in May 2012; At the National Renewable Energy Laboratory National Wind Technology Center (NREL NWTC) during February 2013; At the shoreline of Lake Erie in May 2013; and At the Wind Energy Institute of Canada on Prince Edward Island in May 2015. The experiment we conducted in the coastal zone of Lake Erie indicated very complex flow fields and the frequent presence of upward momentum fluxes and resulting distortion of the wind speed profile at turbine relevant heights due to swells in the Great Lakes. Additionally, our data (and modeling) indicate the frequent presence of low level jets at 600 m height over the Lake and occasions when the wind speed profile across the rotor plane may be impacted by this phenomenon. Experimental data and modeling of the fourth experiment on Prince Edward Island showed that at 10-14 m escarpment adjacent to long-overseas fetch the zone of wind speed decrease before the terrain feature and the increase at (and slightly downwind of) the escarpment is ~3–5% at turbine hub-heights. Additionally, our measurements were used to improve methods to compute the uncertainty in lidar-derived flow properties and to optimize lidar-scanning strategies. For example, on the basis of the experimental data we collected plus those from one of our research partners we advanced a new methodology to estimate a priori the uncertainty in wind speed retrievals from arc scans based on site characteristics such as wind velocity, turbulence intensity and proposed scan geometry. Insights regarding use of remote sensing technologies deriving from project experiments were used to compile a best practice document http://doi.org/10.7298/X4QV3JGF for measuring wind speeds and turbulence offshore through in-situ and remote sensing technologies. A project-specific web-site was developed and is available at: http://www.geo.cornell.edu/eas/PeoplePlaces/Faculty/spryor/DoE_AIATOWEA/index.html

17 WIND ENERGY↗

Failure Mode Classification for Life Prediction Modeling of Solid-State Lighting

Since the passing of the Energy Independence and Security Act of 2007, the U.S. government has mandated greater energy independence which has acted as a catalyst for accelerating and facilitating research efforts toward the development and deployment of market-driven solutions for energy-saving homes, buildings and manufacturing, as well as sustainable transportation and renewable electricity generation. As part of this effort, an emphasis toward advancing solid-state lighting technology through research, development, demonstration, and commercial applications is assisting in the phase out of the common incandescent light bulb, as well as developing a more economical lighting source that is less toxic than compact fluorescent lighting. This has led lighting manufacturers to pursue SSL technologies for a wide range of consumer lighting applications. An SSL luminaire’s lifetime can be characterized in terms of lumen maintenance life. Lumen maintenance or lumen depreciation is the percentage decrease in the relative luminous flux from that of the original, pristine luminous flux value. Lumen maintenance life is the estimated operating time, in hours, when the desired failure threshold is projected to be reached at normal operating conditions. One accepted failure threshold of SSL luminaires is lumen maintenance of 70% -- a 30% reduction in the light output of the luminaire. Currently, the only approved lighting standard that puts forth a recommendation for long-term luminous flux maintenance projections towards a specified failure threshold of an SSL luminaire is the IES TM-28-14 (TM28) standard. iii TM28 was derived as a means to compare luminaires that have been tested at different facilities, research labs or companies. TM28 recommends the use of the Arrhenius equation to determine SSL device specific reaction rates from thermally driven failure mechanisms used to characterize a single failure mode – the relative change in the luminous flux output or “light power” of the SSL luminaire. The use of the Arrhenius equation necessitates two different temperature conditions, 25°C and 45°C are suggested by TM28, to determine the SSL lamp specific activation energy. One principal issue with TM28 is the lack of additional stresses or parameters needed to characterize non-temperature dependent failure mechanisms. Another principal issue with TM28 is the assumption that lumen maintenance or lumen depreciation gives an adequate comparison between SSL luminaires. Additionally, TM28 has no process for the determination of acceleration factors or lifetime estimations. Currently, a literature gap exists for established accelerated test methods for SSL devices to assess quality, reliability and durability before being introduced into the marketplace. Furthermore, there is a need for Physics-of-Failure based approaches to understand the processes and mechanisms that induce failure for the assessment of SSL reliability in order to develop generalized acceleration factors that better represent SSL product lifetime. This and the deficiencies in TM28 validate the need behind the development of acceleration techniques to quantify SSL reliability under a variety of environmental conditions. The ability to assess damage accrual and investigate reliability of SSL components and systems is essential to understanding the life time of the SSL device itself. The methodologies developed in this work increases the understanding of SSL devices iv through the investigation of component and device reliability under a variety of accelerated test conditions. The approaches for suitable lifetime predictions through the development of novel generalized acceleration factors, as well as a prognostics and health management framework, will greatly reduce the time and effort needed to produce SSL acceleration factors for the development of lifetime predictions.

Sakalaukus, Peter Joseph↗

Wake Vortex Inverse Model User's Guide

NorthWest Research Associates (NWRA) has developed an inverse model for inverting landing aircraft vortex data. The data used for the inversion are the time evolution of the lateral transport position and vertical position of both the port and starboard vortices. The inverse model performs iterative forward model runs using various estimates of vortex parameters, vertical crosswind profiles, and vortex circulation as a function of wake age. Forward model predictions of lateral transport and altitude are then compared with the observed data. Differences between the data and model predictions guide the choice of vortex parameter values, crosswind profile and circulation evolution in the next iteration. Iterations are performed until a user-defined criterion is satisfied. Currently, the inverse model is set to stop when the improvement in the rms deviation between the data and model predictions is less than 1 percent for two consecutive iterations. The forward model used in this inverse model is a modified version of the Shear-APA model. A detailed description of this forward model, the inverse model, and its validation are presented in a different report (Lai, Mellman, Robins, and Delisi, 2007). This document is a User's Guide for the Wake Vortex Inverse Model. Section 2 presents an overview of the inverse model program. Execution of the inverse model is described in Section 3. When executing the inverse model, a user is requested to provide the name of an input file which contains the inverse model parameters, the various datasets, and directories needed for the inversion. A detailed description of the list of parameters in the inversion input file is presented in Section 4. A user has an option to save the inversion results of each lidar track in a mat-file (a condensed data file in Matlab format). These saved mat-files can be used for post-inversion analysis. A description of the contents of the saved files is given in Section 5. An example of an inversion input file, with preferred parameters values, is given in Appendix A. An example of the plot generated at a normal completion of the inversion is shown in Appendix B.

Lai, David↗

ROBUST ONLINE MONITORING FOR CALIBRATION ASSESSMENT OF TRANSMITTERS AND INSTRUMENTATION

Robust online monitoring (OLM) technologies are expected to enable the extension or elimination of periodic sensor calibration intervals in operating and new reactors. Specifically, the next generation of OLM technology is expected to include newly developed advanced algorithms that improve monitoring of sensor/system performance and enable the use of plant data to derive information that currently cannot be measured. These advances in OLM technologies will improve the safety and reliability of current and planned nuclear power systems through improved accuracy and increased reliability of sensors used to monitor key parameters. In this paper, we discuss an overview of research being performed within the Nuclear Energy Enabling Technologies (NEET)/Advanced Sensors and Instrumentation (ASI) program, for the development of OLM algorithms to use sensor outputs and, in combination with other available information, 1) determine whether one or more sensors are out of calibration or failing and 2) replace a failing sensor with reliable, accurate sensor outputs. Algorithm development is focused on the following OLM functions: • Signal validation – fault detection and selection of acceptance criteria • Virtual sensing – signal value prediction and acceptance criteria • Response-time assessment – fault detection and acceptance criteria selection A GP-based uncertainty quantification (UQ) method previously developed for UQ in OLM, was adapted for use in sensor-fault detection and virtual sensing. For signal validation, the various components to the OLM residual (which is computed using an AAKR model) were explicitly defined and modeled using a GP. Evaluation was conducted using flow loop data from multiple sources. Results using experimental data from laboratory-scale flow loops indicate that the approach, while capable of detecting sensor drift, may be incapable of discriminating between sensor drift and model inadequacy. This may be due to a simplification applied in the initial modeling, where the sensor degradation is assumed to be stationary. In the case of virtual sensors, the GP model was used in a predictive mode to estimate the correct sensor reading for sensors that may have failed. Results have indicated the viability of using this approach for virtual sensing. However, the GP model has proven to be computationally expensive, and so alternative algorithms for virtual sensing are being evaluated. Finally, automated approaches to performing noise analysis for extracting sensor response time were developed. Evaluation of this technique using laboratory-scale data indicates that it compares well with manual techniques previously used for noise analysis. Moreover, the automated and manual approaches for noise analysis also compare well with the current “gold standard”, hydraulic ramp testing, for response time monitoring. Ongoing research in this project is focused on further evaluation of the algorithms, optimization for accuracy and computational efficiency, and integration into a suite of tools for robust OLM that are applicable to monitoring sensor calibration state in nuclear power plants.

calibration interval extension↗

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗

Acoustic sensors for fission gas characterization: R and D skills devoted to innovative instrumentation in MTR, non-destructive devices in hot lab facilities and specific transducers for measurements of LWR rods in nuclear plants

First of all, we will present the main principle of the method. A piezoelectric transducer, driven by a pulse generator, generates the acoustic waves in a cavity that may be the fuel rod or a chamber connected to an instrumented rod. The composition determination consists in measuring the time of flight of the acoustic signal emitted. The pressure can be estimated by a calibration process, above the measurement of the amplitude of the signal. Two projects will then be detailed. The first project consists in the development of advanced instrumentation for in-pile experiments in Material Testing Reactor. It constitutes a main goal for the improvement of the nuclear fuel behavior knowledge. This acoustic method was tested with success during a first experiment called REMORA 3, and the results were used to differentiate helium and fission gas release kinetics under transient operating conditions. This experiment was lead at OSIRIS reactor (CEA Saclay, France). As a first step of the development program, we performed in-pile tests on the most sensitive component, i.e., the piezoelectric transducer. For this purpose, the active part of this sensor has been qualified on gamma and neutron radiations and at high temperature. Various industrial piezo-ceramics were exposed to a high activity Cobalt source for few days. The cumulated dose was ranged from 50 kGy up to 2 MGy. Next, these devices were placed inside a Material Test Reactor to investigate their reliability towards neutron fluence. The final fluence after 150 days of irradiation was up to 1.6.10{sup 21}n/cm{sup 2} (for thermal neutron). Irreversible variations have been measured. Next, a specific sensor has been implemented on an instrumented fuel rod and tested in the frame of a REMORA 3 Irradiation test. It was the first experiment under high mixed, temperature neutron and gamma flux. A first irradiation phase took place in March 2010 in the OSIRIS reactor and in November 2010 for the second step of the irradiation. The instrumented fuel rod incorporating the ultrasonic gas composition sensor was finally irradiated during 2 weeks in nominal conditions. Neutronics calculation will be performed in order to calculate the thermal and fast neutron fluence and the gamma dose absorbed by acoustic sensor. A first evaluation gives a thermal fluence about 4,5.10{sup 19} n/cm{sup 2}, a fast neutrons fluence about 4,5.1018 n/cm{sup 2} and a total gamma dose up to 0,25 MGy The maximal temperature during the irradiation test was about 150 C. Although the ultrasonic sensor appears to be damaged, the optimization of the electrical attack parameters and the development of a new signal processing maintain the measurement feasibility up the end of the irradiation campaign. It was the first time that the composition of fission gas has been monitored all along an irradiation experiment in a MTR, giving access to the gas release kinetics. New researches involve thick film transducers produced by screen-printing process in order to propose piezoelectric structures for harsh temperature and irradiation measurements. The second project consists in the development of a non-destructive device that can be directly applied on a LWR fuel rod. The problem to be solved relates to the measurement of the fission gas pressure and composition in a fuel rod using a non-destructive method. Fuel rod internal pressure is one of the safety criteria applied in nuclear power analyses. This criterion must be verified in order to avoid any fuel-cladding gap reopening risk and therefore any local clad ballooning. Apart from the safety implications, this parameter is also a fuel behaviour indicator and reflects the overall fuel performance in operation, but also during shipping and long-term storage. Rod internal pressure is one criterion amongst others, like cladding corrosion, against which the acceptable fuel burn-up limit is set. A sensor has been achieved in 2007. A full-scale hot cell test of the internal gas pressure and composition measurement by an acoustic sensor was conducted successfully between 2008 and 2010 on 5 high burn-up MOX fuel rods and 2 very high burn-up UO{sub 2} fuel rods in LECA Facility at Cadarache Centre. An improvement of this sensor has been proposed, allowing us to divide by two the uncertainty on the pressure measurement. In the case of hot-cell measurements, viscous liquid can be used to couple the sensor with the rod. For gas content with a pressure exceeding 15 bars and a 10% Xe/Kr ratio, such coupling may reduce relative acoustic method accuracy by ±7% for pressure measurement result and ±0.25 % for the assessment of gas composition. These results make it possible to demonstrate the feasibility of the technique on LWR fuel rods. The transducer and the associated methodology are now operational for non-destructive measurements in hot lab facilities and allow characterising the fission gas without puncturing the fuel rods. Up to now, any other non-destructive method can be proposed. A next step will be the development of an industrial application in a fuel storage pool in order to perform a large number of measurements on a fuel assembly in nuclear plants.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

The Acoustic Emission signal acquired by the microphones placed in the CABRI test device along the fourteen last R.I.A. experiments: an example of reproducible research in nuclear science

The CABRI facility is an experimental nuclear reactor of the French Atomic Energy Commission. It is located at the Cadarache Research Centre in southern France and it is designed to act as a support to the French nuclear infrastructure. The purpose of the new testing programme termed, 'CABRI International Programme' (CIP) is to study the behaviour of PWR-type fuel rods at high burnup, equipped with an 'advanced' cladding, under Reactivity Initiated Accident (RIA) conditions (such as the scenario of a control rod ejection). Within the framework of this programme, piloted and funded by the French Institute of Nuclear Radioprotection and Safety (IRSN), ten tests are to be conducted with a frequency of two tests per year. The LPRE laboratory of the CEA which is in charge of the Preparation, realisation and breakdown of the test results studies the possibility to set up a new test analysis based on the processing of signals coming from sensors placed within the test equipment. During the experimental phase, the behaviour of the fuel element generates acoustic waves which can be detected by two microphones placed upstream and downstream of the test device. Studies showed the interest to realize temporal and spectral analyses on these signals by showing the existence of signatures which can be correlated with physical phenomena as the rod failure or the test shutdown (i.e. scram). The work presented in this article results from the will to consolidate these studies. Since the main phenomenon to be tracked is the fuel rod failure, the aim would be to highlight specific events which would have been precursors of the rod failure in order to use in the future these signals for further interpretation. Such an antecedent works resumption leads to a better understanding of the experimental needs and constitutes a good initial state to prepare the new very fast digital data acquisition systems. Although all the raw data are accessible in the form of text files, analyses and graphics representations were not straightforward to reproduce from the ancient studies since that, on one hand, people who were in charge of the original work left the laboratory and on the other hand because it is not easy when the time passes, even with our own work, to be able to remember the steps of data manipulations and the exact setup: - During the ancient experiments the use of analog data acquisition systems required to digitize tapes to be able to realize computer treatments. That had had for consequence to lose the initial dating. This one must be correctly edited to do temporal comparisons. - Analyses require functions for calculations whose parameters has to be well-known to reach the same results. We thus wished to manage our workflow in the idea that it can be easily reproducible on all the experiments. The object of the work presented in this article was to put in practice this strong bind between the data, treatments and generation of the document in order not to hesitate to do the iteration principle in action. We do not have to be afraid by the data driven analyses. According to the philosophy of the literate programming, the text of the technical document is woven with the computer code that produces all the printed output as tables, graphs for the study eliminating hence the unrealistic cut and paste. This difficulty is not specific to the nuclear domain. For many years, researchers have been worked out solutions to this mundane issue. And, presently, new technologies and high-level programming languages offer us actual answers. We will firstly present the tools applied in our laboratory to implement this workflow, then we will describe the global perception carried out to continue the study of the Acoustic Emission signals recorded by the two microphones during the fourteen last CABRI R.I.A. test.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An Integrated Supercritical System for Efficient Produced Water Treatment and Power Generation

The main goal of this project was to evaluate the feasibility of an integrated supercritical (SC) water desalination and power generation system for treatment of produced waters from CO 2 sequestration, oil production, and coal-bed methane recovery. This system may offer a transformative approach to generating power from coal or natural gas (NG) and purifying saline water or produced water with high total dissolved solids (TDS) in a zero-liquid discharge plant. The closed-loop steam cycle of conventional power plants was replaced with an open-loop steam cycle that heated the pretreated produced water above the SC point of water (374 °C and 221 bar) to precipitate and separate the salts dissolved in the water. After separation of the precipitated salt, the SC steam required further treatment to reduce the salt content from ~100 mg/L (ppm) to <1 ppm to make it suitable for use in turbines for power generation. A SC membrane distillation unit was used to polish the SC steam, which could be sent through a series of turbines for power generation, and the steam exiting the low-pressure turbine could be condensed by cooling water to produce pure water. Project tasks included a process simulation and a techno-economic evaluation of the integrated system; design and assembly of SC salt precipitation and membrane distillation systems; development and characterization of advanced carbon membranes for SC membrane distillation; and pretreatment, desalination, and purification of different produced water samples with salt concentrations of 30,000–200,000 ppm to a high-purity (<1 ppm) level. Process simulations were performed for 550 MW SC coal-fired or NG baseline cases and various scenarios for the proposed integrated SC system using produced water as the water source and coal or NG as the fuel source. Simulation results indicated that the net plant efficiency loss of the proposed SC cogeneration process would be in the range of <1% to ~4%, in which NG and lower TDS produced water cases would have higher efficiencies. The proposed process for coal-fired cases could provide all water needs for the power generation process and generate additional pure water for beneficial use. A techno-economic analysis was performed to estimate the cost of produced water treatment (dollars /kgal of treated water) by conventional treatment processes compared with the proposed cogeneration process. The cost estimation for the proposed cogeneration system included the cost of pretreatment, the cost of power loss to integrate the water desalination process into the power generation system, and the cost of plant modification or retrofitting. The cost of produced water treatment was calculated by assuming three different scenarios in which converting a baseline power plant to the proposed cogeneration plant would increase the cost of the boiler and feed water system by 10%, 50%, and 100%. Results indicated that the cost of produced water treatment, even with a 100% cost increase, was in the range of $8.6–11.6 dollars per kgal, which is about half of the treatment cost by conventional methods. The cost of produced water treatment by the proposed cogeneration system could be further reduced or a considerable income could be generated if the values of saltwater and freshwater products are considered. Overall, the techno-economic analysis results indicated that the proposed integrated SC system could potentially treat produced water at a cost significantly lower than the cost of produced water treatment by conventional technologies. To evaluate the feasibility of the proposed treatment system at a laboratory scale, a SC salt precipitation and membrane distillation system was designed and fabricated to treat a simulated brine or pretreated produced water by precipitating salts and further polishing the SC steam containing ~100 ppm salts to a high-purity steam. Through baseline tests with NaCl solutions, with concentrations ranging from 30,000 to 200,000 ppm, it was demonstrated that concentrated feed solutions could be treated by precipitating salt at SC conditions to a TDS level of 130 ppm. Produced water samples were collected from a potential CO 2 sequestration site, a coal-bed methane site, and an oilfield in Illinois. Representative samples were characterized for various water quality parameters, including pH; turbidity; conductivity; and concentrations of selected anions, cations, alkalinity, dissolved organic carbon, total suspended solids, TDS, and total petroleum hydrocarbons. The produced water samples were first treated by conventional pretreatment methods to remove residual oil, dissolved organic matter, and scale-forming species. Pretreatment methods included coagulation, sedimentation, filtration, activated carbon adsorption, and ion-exchange adsorption. Pretreated produced waters were treated at SC conditions in which salts were precipitated. Analysis of treated samples revealed that the concentration of Na (the main cation) in produced waters was reduced dramatically from 14,000–47,400 ppm in the feed to 21–80 ppm in treated waters. Concentrations of other main cations (Ca, Mg, Sr, K) also decreased significantly from initial values of ~1,000–22,000 ppm to ~1 ppm or less. However, a significant increase in the concentrations of Fe, Ni, Cu, and Co was observed, mainly caused by leaching of these elements from the system (vessels, filters, tubing), which is manufactured from alloys containing Fe, Ni, Cu, and Co. To polish the SC steam containing a salt content of ~100 ppm to a high-purity level of <1 ppm, a SC membrane distillation treatment using carbon membranes was employed. Different types of self-standing or supported carbon membranes were developed from graphite, graphene, carbon nanotubes, carbon nanofibers, and pyrolytic carbon. Screening tests revealed that only a limited number of developed carbon membranes were suitable for the purification of steam at SC conditions. Membranes that were either fragile and lacked mechanical strength or had large pores and showed limited salt rejection (<20%) were found to be unsuitable for the purification of SC steam. Promising carbon membranes that were successfully tested for purifying SC steam included two graphite membranes (i.e., AF-5, a porous graphite material, and FGS, a flexible graphite sheet) and a carbon nanotube membrane grown on a Hastelloy substrate and then further modified using a proprietary method (i.e., membrane C). The AF-5 membrane exhibited a stable water flux of ~50 kg/m 2 ∙h and ~40% salt rejection. The FGS membrane had a high salt rejection of 95%–98% but a lower water flux of ~15 kg/m 2 ∙h. Membrane C showed an extremely high flux of 200-800 kg/m 2 ∙h, which is one order of magnitude higher than the water flux of conventional membranes. Membrane C also exhibited 87%–93% salt rejection. Overall, a two-stage or multistage polishing treatment by these membranes could lower the salt content of the SC steam from initial values of ~100–200 ppm to the steam quality of <1 ppm or to parts per billion levels required for the power cycle.

20 FOSSIL-FUELED POWER PLANTS↗

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↗