Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “efficiency table”

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 685 records · Page 38

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation of a Perforated Plate Cooling Flow Part I

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hotor combusting gases. The design of more efficient film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained to assess current CFD capabilities and develop better modeling approaches for these turbulent flow fields where accurate calculation of turbulent heat flux is important. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean and root mean square (rms) velocities as well as the mean and rms temperatures in a film cooling flow field. In addition to off-body flow field measurements, infrared thermography (IR) and thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures (up to TR =2.7) and Mach numbers (up to Mach 0.9), across a 30.48 cm long plate equipped with three patches of 45 small (~1 mm) diameter cooling holes arranged in a staggered configuration. In addition, both streamwise 2-component PIV (along the plate centerline) and cross-stream 3-component Stereo PIV data at 14 axial stations were collected in the same flows. Only a subset of the data collected in the test program is included in this Part 1 report. The rest of the data will be published in a future report, Part 2, along with planned CFD predictions of the complex cooling film flow.The entire data set of Raman temperature data, PIV velocity data and IR camera data covering the Set Points 23 and 49 in the test matrix in Table 1 is available in an accompanying DVD (available online from www.sti.nasa.gov) for those interested in further analysis.

Mark P Wernet↗

Using the cFS Command and Data Dictionary (CCDD) to Automate Software Development on Habulous

Final paper is attached. The NASA developed Core Flight System (cFS) is a reusable software architecture that has been used on multiple spaceflight missions. By using this framework, missions are able to reuse code from other missions, as well as leverage deployment onto similar computer architectures (i.e. not "reinvent the wheel" on each new mission). The success in the cFS concept can be seen in the large number of projects using cFS at FSW-2018. The Habulous project is an Earth-based testbed, used for hardware and software that may one day be used on a future space habitat unit, with many participating groups from various NASA centers and aerospace organizations around the country. The distributed nature of the various teams mean that defining (and following) an interface definition is critical on the project. Additionally, since various groups use various types of computer hardware (32/64-bit, big/little endian, Linux/VxWorks/Windows) many additional complications exist in interfacing all the various components into a final integrated system. cFS is used on the majority the flight software (FSW) in running in Habulous. But some subsystems have elected to not use cFS, and use a software bridge (called SBN_lib) to interact with the other cFS nodes in Habulous. In order to most efficiently develop the FSW, a central database is used to define and store each message sent by cFS. A Command and Data Dictionary (CDD) is something nearly universal on spacecraft, but as a team we worked to develop the CDD before the SW development was complete, and not treat it like "as built" documentation. To manage the CDD, the cFS Command and Data Dictionary (CCDD) tool was chosen (available from NASA as open source software). The CCDD tool has successfully been used to automate/autocode a large amount of software used on Habulous, as we are hoping to use it to define even more items in the future (time-triggered Ethernet (TTE) network maps, CPU scheduling). Additionally, Habulous has been exploring the use of cFS on wildly heterogeneous CPUs, and how to coordinate all those various machines using/extending the software bus – network (SBN) application in cFS, as well as TTE to coordinate message passing between various synchronized machines. The major topics to be covered in the presentation are: (1) Updating to the CCSDS_v2 extended headers (and using CPU# as subsystem ID). (2) Managing all the message identification numbers for each cFS message sent/received on any of the various CPUs. (3) Using the CCDD information to automatically generate the C-header files that define the structure for all software bus (SB) commands/telemetry messages. (4) Using the CCDD to automatically generate XML Telemetry and Command Exchange (XTCE) files, which streams display production/integration/testing in a web based display architecture (5) Extending/customizing SBN to pass messages among computers on multiple networks. (6) Using "Protobetter" inside SBN to manage different endian-ness/architectures. (7) Using SBN_lib to allow non-cFS node to communicate with cFS nodes. (8) Developing TTE network and schedule tables for all the various CPUs to use.

Hirsh, Robert L.↗

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation of a Perforated Plate Cooling Flow: Part I

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data was obtained to assess current CFD capabilities and develop better modeling approaches for these turbulent flow fields where accurate calculation of turbulent heat flux is important. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean and root mean square (rms) velocities as well as the mean and rms temperatures in a film cooling flow field. In addition to off-body flow field measurements, infrared thermography (IR) and thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures (up to TR = 2.7) and Mach numbers (up to Mach 0.9), across a 30.48 cm long plate equipped with three patches of 45 small (~1 mm) diameter cooling holes arranged in a staggered configuration. In addition, both streamwise 2-component PIV (on the plate centerline) and cross-stream 3-component Stereo PIV data at 14 axial stations were collected in the same flows. Only a subset of the data collected in the test program is included in this Part I report. The rest of the data will be published in a future report, Part II, along with planned CFD predictions of the complex cooling film flow. The entire data set of Raman temperature data, PIV velocity data and IR camera data covering the Set Points 23 and 49 in the test matrix in Table 1 is available in an accompanying DVD (available online from www.sti.nasa.gov) for those interested in further analysis.

Nozzle↗

A Review of Value of Solar Studies In Theory and In Practice

This brief summarizes a collection of state- and utility-commissioned value-of-solar (VoS) studies and related literature, with a focus on who commissioned the study, which value and cost categories were discussed and/or quantified, and the methods used. Our objective is to compile information on prior VoS studies to inform state regulators and other stakeholders that may pursue related studies or integrate findings into rate design. The brief is organized into three parts: 1) an introduction to distributed solar photovoltaic (DPV) compensation; 2) a review of theoretical research on VoS; and 3) a review of VoS studies. The vast majority of VoS studies have served an informational role of quantifying the net benefits of PV. Three studies were commissioned in states or utility service territories that subsequently implemented VoS tariffs in California, New York, and Austin, Texas. When applied as a tariff, VoS aims to compensate PV output as efficiently as possible by doing so at rates that reflect the marginal benefits and costs of PV through value and cost categories that may vary temporally and/or geographically. This could lead to higher compensation in locations and times where more PV output is more valuable and consequently drive adoption in those locations to provide more societal benefits. Value and cost factors can be broadly grouped into five categories: generation, transmission, distribution, other utility, and other social categories. Those conducting VoS studies must weigh various tradeoffs when deciding which categories to include and quantify. Tradeoffs include prioritizing values based on their magnitude of value or cost impact, as well as taking into account the feasibility of data collection and accurate quantification. Values of higher magnitude and estimation feasibility are quantified in the majority of studies, including the earliest of studies conducted in the 2000s and 2010s. Additionally, some values of higher magnitude but low feasibility in the earliest of studies have become quantifiable in recent years. There are some values with low average system-wide levels but very high magnitude in specific locations or hours. The value magnitude in some cases can be tied to DPV penetration with low value in areas with little congestion and/or low penetration and vice versa. In these cases, values that are easier to quantify are often incorporated, while those that are more difficult are often addressed via a placeholder value. The placeholder value is paired with a discussion around data needs and methods to improve future estimates, as well as a conversation about when these value categories may increase in magnitude and necessitate more rigorous quantification. This brief summarizes findings from two meta-analyses of VoS studies that took place between 2005 and 2018, as well as findings from four additional studies published from 2018 to 2023. Table ES-1 summarizes the various value and cost categories included in each respective study and whether they were quantified, discussed, or omitted. Values such as avoided energy, capacity, transmission capacity, line losses, and avoided environmental costs are quantified in every study. Some categories were deemed harder to quantify and less impactful at the time of the study, so they were discussed but not quantified (e.g., ancillary services). Other categories, including many at the distribution level, were very locationally and/or temporally specific and dependent on high DPV penetration. These were sometimes quantified and at other times discussed. Notably, when it came to utility costs, integration costs were discussed in all cases, though they were deemed to have a small impact. Other utility costs were omitted for the most part; however, the utility-commissioned study (by NorthWestern Energy in Montana) included both lost utility revenue and programmatic/administrative cost categories. While there are some similarities across studies, each had fairly unique methods that are detailed in the body of this brief.

14 SOLAR ENERGY↗

Organic matter and soil structure in the Everglades Agricultural Area

This publication pertains to management of organic soils (Histosols) in the Everglades Agricultural Area (EAA). These former wetland soils are a major resource for efficient agricultural production and are important globally for their high organic matter content. Recognition of global warming has led to considerable interest in soils as a repository for carbon. Soils rich in organic matter essentially sequester or retain carbon in the profile and can contribute directly to keeping that sequestered carbon from entering the atmosphere. Identification and utilization of management practices that minimize the loss of carbon from organic soils to the atmosphere can minimize effects on global warming and increase the longevity of subsiding Histosols for agricultural use. Understanding and predicting how these muck soils will respond to current and changing land uses will help to manage soil carbon. The objectives of this document are to: a. Discuss organic soil oxidation relative to storing or releasing carbon and nitrogen b. Evaluate effects of cultivation (compare structure for sugarcane vs. uncultivated soil) Based upon the findings from the land-use comparison (sugarcane or uncultivated), organic carbon was higher with cultivation in the lower depths. There is considerable potential for minimum tillage and residue management to further enhance carbon sequestration in the sugarcane system. Carbon sequestration is improved and soil subsidence is slowed with sugarcane production, and both of these are positive outcomes. Taking action to increase or maintain carbon sequestration appears to be appropriate but may introduce some risk to farming operations. Additional management methods are needed to reduce this risk. For both the longevity of these organic soils and from a global perspective, slowing subsidence through BMP implementation makes sense. Since these BMPs also have considerable societal benefit, it remains to be seen if society will help to offset a part or all of the additional risk through some form of cost-sharing program or carbon credits trading. In general, the subsidence throughout the EAA has been slowed because of higher water table management and implementation of other selected BMPs. In addition, the comparison of soil with different land uses shows that the humification rate, conversion of organic matter from peat to humus, has changed. Another likely factor is a relative increase in the mineral content of soil as the organic constituents are lost through subsidence.

54 ENVIRONMENTAL SCIENCES↗

Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays

Single-slope analog-to-digital converters (ADCs) are particularly useful for onchip digitization in focal plane arrays (FPAs) because of their inherent monotonicity, relative simplicity, and efficiency for column-parallel applications, but they are comparatively slow. Squareroot encoding can allow the number of code values to be reduced without loss of signal-to-noise ratio (SNR) by keeping the quantization noise just below the signal shot noise. This encoding can be implemented directly by using a quadratic ramp. The reduction in the number of code values can substantially increase the quantization speed. However, in an FPA, the fixed pattern noise (FPN) limits the use of small quantization steps at low signal levels. If the zero-point is adjusted so that the lowest column is onscale, the other columns, including those at the center of the distribution, will be pushed up the ramp where the quantization noise is higher. Additionally, the finite frequency response of the ramp buffer amplifier and the comparator distort the shape of the ramp, so that the effective ramp value at the time the comparator trips differs from the intended value, resulting in errors. Allowing increased settling time decreases the quantization speed, while increasing the bandwidth increases the noise. The FPN problem is solved by breaking the ramp into two portions, with some fraction of the available code values allocated to a linear ramp and the remainder to a quadratic ramp. To avoid large transients, both the value and the slope of the linear and quadratic portions should be equal where they join. The span of the linear portion must cover the minimum offset, but not necessarily the maximum, since the fraction of the pixels above the upper limit will still be correctly quantized, albeit with increased quantization noise. The required linear span, maximum signal and ratio of quantization noise to shot noise at high signal, along with the continuity requirement, determines the number of code values that must be allocated to each portion. The distortion problem is solved by using a lookup table to convert captured code values back to signal levels. The values in this table will be similar to the intended ramp value, but with a correction for the finite bandwidth effects. Continuous-time comparators are used, and their bandwidth is set below the step rate, which smoothes the ramp and reduces the noise. No settling time is needed, as would be the case for clocked comparators, but the low bandwidth enhances the distortion of the non-linear portion. This is corrected by use of a return lookup table, which differs from the one used to generate the ramp. The return lookup table is obtained by calibrating against a stepped precision DC reference. This results in a residual non-linearity well below the quantization noise. This method can also compensate for differential non-linearity (DNL) in the DAC used to generate the ramp. The use of a ramp with a combination of linear and quadratic portions for a single-slope ADC is novel. The number of steps is minimized by keeping the step size just below the photon shot noise. This in turn maximizes the speed of the conversion. High resolution is maintained by keeping small quantization steps at low signals, and noise is minimized by allowing the lowest analog bandwidth, all without increasing the quantization noise. A calibrated return lookup table allows the system to maintain excellent linearity.

Wrigley, Chris J.↗

Metrology of Large Parts

As discussed in the first chapter of this book, there are many different methods to measure a part using optical technology. Chapter 2 discussed the use of machine vision to measure macroscopic features such as length and position, which was extended to the use of interferometry as a linear measurement tool in chapter 3, and laser or other trackers to find the relation of key points on large parts in chapter 4. This chapter looks at measuring large parts to optical tolerances in the sub-micron range using interferometry, ranging, and optical tools discussed in the previous chapters. The purpose of this chapter is not to discuss specific metrology tools (such as interferometers or gauges), but to describe a systems engineering approach to testing large parts. Issues such as material warpage and temperature drifts that may be insignificant when measuring a part to micron levels under a microscope, as will be discussed in later chapters, can prove to be very important when making the same measurement over a larger part. In this chapter, we will define a set of guiding principles for successfully overcoming these challenges and illustrate the application of these principles with real world examples. While these examples are drawn from specific large optical testing applications, they inform the problems associated with testing any large part to optical tolerances. Manufacturing today relies on micrometer level part performance. Fields such as energy and transportation are demanding higher tolerances to provide increased efficiencies and fuel savings. By looking at how the optics industry approaches sub-micrometer metrology, one can gain a better understanding of the metrology challenges for any larger part specified to micrometer tolerances. Testing large parts, whether optical components or precision structures, to optical tolerances is just like testing small parts, only harder. Identical with what one does for small parts, a metrologist tests large parts and optics in particular to quantify their mechanical properties (such as dimensions, mass, etc); their optical prescription or design (i.e. radius of curvature, conic constant, vertex location, size); and their full part shape. And, just as with small parts, a metrologist accomplishes these tests using distance measuring instruments such as tape measures, inside micrometers, coordinate measuring machines, distance measuring interferometers; angle measuring instruments such as theodolites, autocollimators; and surface measuring instruments including interferometers, stylus profilers, interference microscopes, photogrammetric cameras, or other tools. However, while the methodology may be similar, it is more difficult to test a large object for the simple reason that most metrologists do not have the necessary intuition. The skills used to test small parts or optics in a laboratory do not extrapolate to testing large parts in an industrial setting any more than a backyard gardener might successfully operate a farm. But first, what is a large part? A simple definition might be the part's size or diameter. For optics and diffuse surface parts alike, the driving constraint is ability to illuminate the part's surface. For reflective convex mirrors, large is typically anything greater than 1 meter. But, for refractive optics, flats or convex mirrors, large is typically greater than 0.5 meter. While a size definition is simple, it may be less than universal. A more nuanced definition might be that a large part is any component which cannot be easily tested in a standard laboratory environment, on a standard vibration isolated table using standard laboratory infrastructure. A micro-switch or a precision lens might be easily measured to nanometer levels under a microscope in a lab, but a power turbine spline or a larger telescope mirror will not fit under that microscope and may not even fit on the table.

Stahl, H. Philip↗

THF Co-Solvent Biomass Fractionation to Catalytic Fuel Precursors with High Yields (Final Report)

The overall goal of this project is to utilize the Co-Solvent Enhanced Lignocellulosic Fractionation (CELF) as a platform technology to realize transformation of real biomass feedstocks to fungible transportation fuels. CELF is a solvolysis method that seeks to break the cost barrier for biomass deconstruction so that high yield liquid intermediates from both sugars and lignin can be recovered for their subsequent catalytic conversion to fungible transportation fuels. It is important for the DOE EERE mission to focus on supporting platform technologies that enable greater end-to-end utilization of biomass whereby the co-development of both process and catalyst is streamlined. CELF is one of the most competitive platform technologies that uses multifunctional co-solvents THF and water to provide the cleanest intermediate product streams that can allow advanced catalysts to be developed on real feedstocks instead of model compounds. For this project, we developed a hybrid catalytic process that first applies CELF homogeneous catalysis to convert solid biomass to a liquid stream containing intermediate fuel precursors (FPs) such as furfural (FF) and 5-hydroxymethylfurfural (5-HMF). Afterwards, a second heterogeneous catalytic step was then applied to selectively hydrodeoxygenate (HDO) the FPs produced from biomass into gasoline blendstocks methylfuran (MF) and dimethylfuran (DMF). Concurrently, capitalizing on CELF’s efficient lignin extraction capabilities, we further explored production C9+ cyclic hydrocarbons from biomass lignin to improve overall process economics. For the first step, the CELF process was optimized to achieve yields of >95% for FF and >60% for HMF from the xylan and glucan fractions in poplar wood, respectively. For the second step, we created a patent-pending bifunctional Cu-Ni solid catalyst that is tolerant to real product streams from the CELF process and simultaneously produce MF and DMF blendstocks at >80% yields. To provide more market versatility, we further explored the production of 2-pentanol from MF and have achieved >70% yields. During this project’s three year period, we also developed a dynamically informed technoeconomic computer model that helped to guide key project decisions that allowed us to reach a final biofuel selling price (MFSP) of <$\$$3/GGE (gallon gasoline equivalent) and total biomass carbon utilization of 60%. We also determined that the MFSP in our model could be further reduced to below $\$$2/GGE if the value of the lignin-derived aromatic hydrocarbon fuels could be sold at diesel-range prices (>$\$$500/tonne). Our process model was designed to consume 288 tonnes per day of poplar wood chips with minor consumption of hydrogen to produce a range of fuel blendstocks suitable for use in the transportation sector. The project was successful in meeting DOE EERE program criteria for application of an all-catalytic biofuels process to achieve a biofuel selling price target of <$\$$3/GGE as well as to achieve a 60% reduction in greenhouse gas emissions (GHG) compared to conventional fossil petroleum processes. At project completion, the outlined catalytic process reported total utilization of up to 60% of the carbon found in biomass, rivaling more mature biological cellulosic ethanol technologies at a fraction of the fixed and operating costs. This project also resulted in 7 peer-reviewed publications, 2 provisional inventions, and 1 patent pending. NOTE: This public version of the document contains redacted text, tables, and figures as indicated within the document.

09 BIOMASS FUELS↗

Exploration of Near-Term Potential Routes and Procedures for Urban Air Mobility

Urban air mobility is gaining interest as the need for On Demand Mobility in today's congested traffic is becoming high in metropolitan areas. Urban Air Mobility (UAM) is envisioned as a concept to transport passengers and cargo safely and efficiently using innovative aircraft in the urban areas. It is expected to improve mobility for the general public, decongest road traffic, reduce transport time and reduce the strain on existing public transport networks. There exist several challenges to Urban Air Mobility (UAM) such as integration of procedures with airspace and the airport, noise levels that are acceptable to the general public, public safety, public acceptance, vehicle certification, and more. Most of the research in the United States and European skies (DLR - German Aerospace Center) related to urban areas has focused on small UAS (Unmanned Aircraft Systems) flights (NASA's UTM (UAS Traffic Management) research) and their integration with the airspace and building safe operations in densely populated areas. Previous studies on UAM have focused on fast time simulations of the routes that are separated via a separation service and network of routes. Similarly, research in Europe has focused on the approach profile for these innovative aircraft, vertiports and battery life among others. UAM as a part of the On-Demand Mobility effort has provided some guidelines for operations as shown below: Does not require additional ATC (Air Traffic Control) infrastructure; Does not impose additional workload on ATC; Does not restrict operations of traditional airspace users; Will meet appropriate safety thresholds and requirements; Will prioritize operational scalability; Will allow flexibility where possible and structure where necessary. This paper explores potential routes and procedures in a Human-In-The-Loop (HITL) experiment that could be applied in the near-term to allow integration of UAM flights into the airspace as well as a large airport. The airspace that was explored was Dallas Fort Worth (DFW) airspace managed by the DFW East Tower in South Flow only. In addition, Dallas Love Field (DAL) and Addison (ADS) airspace were also part of the testbed. The initial set of routes investigated in this study were published helicopter routes in the DFW area. Figure 1 shows class B airspace in DFW area and the origin/destination city pairs where UAM flights flew along with helicopter routes shown in blue. The research focused on exploring procedures for integrating UAM flights into Class Bravo and Class Delta airspace. Three different communication procedures, evaluated with three different levels of UAM traffic, are shown in Table 1. The current day routes were evaluated with current day communication procedures were explored as the first condition. The current day routes were also evaluated in the second condition with reduced communications, which was assumed due to the presence of a Letter Of Agreement (LOA). The purpose of the LOA was to reduce the verbiage associated with pilots getting clearance to Class B airspace from the controllers, pre-assigning beacons codes to the UAM flights, separate routes by assigning altitudes and speeds to flights going in any one direction. Flights were expected to automatically change frequency when exiting Class B airspace, thus transition points for entry and exit points were also specified in the LOA.

Urban Air Mobility↗

Accelerating Biomimetic Solar - Energy Harvesting: Mapping the Interaction Landscape of Plasmonic-Excitonic Hybrid Nanosystems (Final Report)

In general, excitonic and plasmonic nanoscale materials in close proximity show high potential for significant breakthroughs in energy related materials research. The interactions between these two kinds of materials result in coupled optical transitions (plexcitons), distinct from those of both the individual exciton and plasmon as well as from those of the sum of their constituents (synergistic effects). By linking together materials-research and physical-research approaches, this project contributes to a concerted approach on nanomaterials energy research. The project’s overall goal is to accelerate the development of well-defined plexcitonic model systems consisting of carefully engineered plasmonic and excitonic nanomaterial— essential for both gaining a fundamental understanding of plexcitonic nanomaterials and the development of novel design principles for biomimetic solar energy harvesting. During the 3-year project period and the terminal renewal with limited support for a 12-month period, we successfully synthesized and characterized (1) a robust excitonic nanomaterial and (2) a library of plasmonic nanoparticles as well as developed (3) a microfluidic platform for homogenous nanosynthesis as summarized below: (1) Robust Excitonic Nanomaterial. Supramolecular assemblies are Nature’s most successful material system for solar energy harvesting. However, photovoltaic devices based on artificial supramolecular assemblies continue to be stymied by disappointing efficiencies and poor stability. The conceptual failure may lie in current solar cell architectures, which rely on solidifying supramolecular assemblies as an ensemble into a solid matrix, neglecting the intrinsic fragility of the assemblies’ internal structure, thus disrupting or even destroying their delicate optoelectronic properties, that is, delicate Frenkel excitonic properties. Supramolecular assemblies may finally serve as usable light harvesting material systems for solar energy conversion technologies, only if they meet the following criteria: (a) Stability, that is, the fragile structure including its delicate Frenkel excitonic character needs to be stable, (b) Robustness, that is, resistant against elevated and fluctuating temperatures, and (c) Viability for device integration, that is, capable of being immobilized onto solid substrates. Here, by developing a nanocomposite via a tunable, cage-like scaffold design, we successfully provided stable supramolecular nanocomposites, that inhabit robust Frenkel excitons despite harming environmental conditions such as extreme heat stress. (2) Library of Plasmonic Nanoparticles. Naturally, current models describing plasmonic hybrid quantum states—plasmonic hybridizations—parallel those developed for molecular orbitals, equating individual plasmonic nanostructures with “atoms” and the plasmonic nanoassemblies with “molecules.” In analogy to organic synthesis, a suitably robust fabrication method would allow for “atom-like” manipulation of “molecule-like” plasmonic nanoassemblies; of high value for next-generation energy nanotechnologies. Despite this frequent comparison, current plasmonic nanoassembly fabrication methods favor top-down templating over wet-chemical synthesis, however, achieving precise control over nanostructure’s geometry and surface characteristics remain an art and a scientific challenge. The conceptual failure may lie in the current wet-chemical synthesis paradigm, as it relies on the accessibility of a multi-dimensional synthesis parameter space through limited, rather one-dimensional synthesis procedures by employing step-by-step approaches. Solution-based nanoarchitectonics for rational design of precisely built plasmonic nanoassemblies via solution-based fabrication may finally be possible only if multi-dimensional syntheses approaches are available that allow for comprehensive control over the plasmonic nanomaterials’ (a) Structural Properties and (b) Surface Properties. Here, by developing an innovative multidimensional 1,3-propanediol based polyol synthesis, we successfully provided control over the plasmonic building-block’s geometry (size and shape) together with its surface characteristics. Our results present a critical step toward the vision of a “periodic table-like” system for plasmonic materials based on straightforward wet-chemical syntheses for energy nanotechnologies. Developing deliberate modifications on this synthesis, we generated a library of plasmonic nanostructures covering the vast parameter space—opening the door for fundamental investigation of plexcitonic model systems. (3) Microfluidic Platform for Homogenous Nanosynthesis. Control over structural properties of plexcitonic nanocomposites remains a challenge due to current limitations in nanosynthesis techniques. Slight variations in nanostructure’s geometry impact their optoelectronic properties, demanding precise synthesis beyond the capabilities of solution-based (batch) synthesis processes. In contrast, the small, confined liquid volumes used in microfluidics—a reaction technique where the manipulation of fluids takes place in channels with dimensions of tens of micrometers—allows for homogenous synthesis conditions, providing excellent control of the reaction and, as a result, of the materials’ geopmetry and composition. However, thus far, the majority of plexcitonic systems has been developed via batch synthesis. Here, by successfully developing a two-channel microreactor, our microfluidic-supported synthesis approach combines the advantages of both microfluidics and batch platforms, allowing for precise spatio-temporal control over all synthesis parameters opening the possibility for homogenous nanosythnesis of well-defined plexcitonic model systems.

14 SOLAR ENERGY↗

Environmental performances of different configurations of a material recovery facility in a life cycle perspective

Highlights: • An attributional LCA assessed the environmental performances of an integrated MRF. • Different configurations of the MRF technological evolution have been analysed. • Reliable data provide a large LCI table for all the configurations and each sub-units. • A biological treatment in-situ improves the potentials of land occupation and global warming. • The inclusion of a small-scale gasifier could improve almost all the environmental potentials. - Abstract: The study evaluated the environmental performances of an integrated material recovery facility (MRF) able to treat 32 kt/y of unsorted mixed waste, made of residuals from household source separation and separate collection. The facility includes a mechanical sorting platform for the production of a solid recovered fuel (SRF) utilized in an external waste-to-energy plant, bio-cells for tunnel composting of organic fraction, and a sanitary landfill for the safe disposal of ultimate waste. All the MRF sub-units have been analysed in depth in order to acquire reliable data for a life cycle assessment study, focused on the environmental performances of different configurations of the facility. The study investigated a “past” configuration, including just mechanical sorting, landfilling and biogas combustion in a gas engine, and the “present” one, which includes also a composting unit. Two possible “future” configurations, having a gasifier inside the MRF battery limits, have been also analysed, assessing the performances of two fluidized bed reactors of different size, able to gasify only the residues generated by the sorting platform or the whole amount of produced SRF, respectively. The analysis evaluated the contributions of each unit in the different configurations and allowed a reliable assessment of the technological evolution of the facility. The results quantified the positive effect of the inclusion of an aerobic treatment of the waste organic fraction. The SRF gasification in situ appears to improve the MRF environmental performances in all the impact categories, with the exclusion of that of global warming.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Jointly Sponsored Research Program on Energy Related Research

Cooperative Agreements, DE-FC26-08NT43293, DOE-WRI Cooperative Research and Development Program for Fossil Energy-Related Resources began in June 2009. The goal of the Program was to develop, commercialize, and deploy technologies of value to the nation’s fossil and renewable energy industries. To ensure relevancy and early commercialization, the involvement of an industrial partner was encouraged. In that regard, the Program stipulated that a minimum of 20% cost share be achieved in a fiscal year. This allowed WRI to carry a diverse portfolio of technologies and projects at various development technology readiness levels. Depending upon the maturity of the research concept and technology, cost share for a given task ranged from none to as high as 67% (two-thirds). Over the course of the Program, a total of twenty six tasks were proposed for DOE approval. Over the period of performance of the Cooperative agreement, WRI has put in place projects utilizing a total of $\$$7,089,581 in USDOE funds. Against this funding, cosponsors have committed $\$$7,398,476 in private funds to produce a program valued at $\$$14,488,057. Tables 1 and 2 presented at the end of this section is a compilation of the funding for all the tasks conducted under the program. The goal of the Cooperative Research and Development Program for Fossil Energy-Related Resources was to through collaborative research with the industry, develop or assist in the development of innovative technology solutions that will: • Increase the production of United States energy resources – coal, natural gas, oil, and renewable energy resources; • Enhance the competitiveness of United States energy technologies in international markets and assist in technology transfer; • Reduce the nation's dependence on foreign energy supplies and strengthen both the United States and regional economies; and • Minimize environmental impacts of energy production and utilization. Success of the Program can be measured by several criteria. Using the deployment of the federal funding with industrial participation as a performance criterion, over the course of the program, the cosponsors contributed more dollars than the federal funds. As stated earlier, a little more than half of the funding for the Program was derived from industrial partners. The industrial partners also enthusiastically supported the research and development activities with cash contribution of $\$$4,710,372.67, nearly 65% of the required cost share. Work on all of the tasks proposed under the Cooperative Agreement has been completed. This report summarizes and highlights the results from the Program. Under the Cooperative Agreement Program, energy-related tasks emphasized petroleum processing, upgrading and characterization, coal and biomass beneficiation and upgrading, coal combustion systems development including oxy-combustion, emissions monitoring and abatement, coal gasification technologies including gas clean-up and conditioning, hydrogen and liquid fuels production, and the development of technologies for the utilization of renewable energy resources. Environmental-related activities emphasized cleaning contaminated soils using microbial fuel cells, development of processes and sorbents for emissions reduction and recovery of water from power plant flue gas, and biological carbon capture and reuse. Technology enhancement activities included resource characterization studies, development of improved methods, monitors and sensors. In general the goals of the tasks proposed were to enhance competitiveness of U.S. technology, increase production of domestic resources, and reduce environmental impacts associated with energy production and utilization. Technologies being brought to commercialization as a result of the funds provided by the Cooperative Agreement contribute to the overall goals of the USDOE and the nation. Each has broad applicability both within the United States and abroad, thereby helping to enhance the competitiveness of U.S. energy technologies in international markets and assisting in technology transfer. Under the Cooperative Agreement Program, WRI has furthered the development of two different coal upgrading technologies. River Basin Energy technology was scaled-up and demonstrated at a nominal 40 tpd size. Similarly, WRI’s patented mercury removal technology further developed into WRITE Coal technology which was then integrated into oxy-combustion and gasification systems for IGCC and fuels production. Integrated systems with WRITE Coal technology applied at the front end represent substantial environmental and efficiency gains. A variation of the RBE coal upgrading technology is being commercialized as a torrefaction technology for woody biomass. WRI worked with EPRI and NIST to develop and improve mercury calibration standards for emissions monitoring. Working with Chart Energy and Chemicals, WRI scaled-up compact reactor technology for the synthesis of fuels and chemicals from syngas. Compact reactor technology represents a five-fold increase in productivity over conventional reactors making smaller-scale distributed synthesis plants an economical viability. Similarly, WRI's patented mixed alcohol synthesis catalyst production is being scaled-up in collaboration with a commercial catalyst manufacturer.

99 GENERAL AND MISCELLANEOUS↗

Retrospective on Recent DOE-Funded Studies Concerning the Extraction of Rare Earth Elements & Lithium from Geothermal Brines (Final Report)

Rare earth elements (REE) and lithium are non-toxic metals that are considered critical materials due to their use in electronics, magnets, batteries, and a wide variety of industrial processes important for the economy and military preparedness. Demand for REE and lithium is increasing and these critical materials are imported, so identifying and exploiting domestic sources of REE and lithium is a national priority. The U.S. Department of Energy (DOE) Geothermal Technologies Office (GTO) has been in the forefront of sponsoring research investigating the potential recovery of REE, lithium, and other critical minerals from geothermal brines. It has been proposed that the future of geothermal energy should include “hybrid systems” that combine electricity generation with other revenue-generating activities, such as recovery of valuable and critical minerals, including REE and lithium. Two recent GTO funding opportunities have focused on the recovery of REE and other valuable minerals from geothermal brines. The research supported by the GTO’s mineral recovery program is focused on three areas: resource characterization, technology for the extraction of REE, and technology for the extraction of lithium (Tables 1 and 2). This report is a retrospective study examining the outcome of GTO’s two recent mineral recovery programs (DE-FOA-0001016 in FY 2014 and DE-FOA-0001376 in FY 2016). In this report, the knowledge, technology, and techniques that were developed by researchers funded by GTO are summarized and discussed. Four projects were funded to assess the concentrations and amounts of REE found in geothermal brines and oil field produced waters. The GTO-funded studies compiled publically available data on REE concentrations from brines and produced water from all over the USA. In addition, new samples were collected and characterized from major geothermal and hydrocarbon basins in the Western USA. The studies examined the relationship between lithology and REE concentrations and developed models examining the influence of geology on REE concentrations in produced brines. It was determined that REE are frequently found at higher concentrations in oil field produced water than geothermal brines, but that some geothermal areas had significant REE resources. Significant reservoirs of REE were identified in the Western USA. In some cases, concentrations of REE were more than 1000 times the concentrations found in seawater. Collectively, these studies represent a comprehensive picture of REE resources associated with geothermal and hydrocarbon systems in the USA. The studies did not examine lithium resources, but in some cases, lithium concentration data was collected. Data from these studies are housed in the Geothermal Data Repository (GDR) and represent a significant information resource and it is recommended that these data be further analyzed in a future study. Eight projects were funded to develop new technology for REE extraction from geothermal fluids. These projects investigated sorption as an approach for removal and recovery of REE from geothermal brines. The projects investigated cutting-edge technology for selective sorption of ions from complex solutions, including the application of metal-organic frameworks and biosorbent proteins. The REE sorption studies tested different combinations of metal-binding ligands and solid supports. The most promising metal-binding ligands for REE included phosphonic acid, thiol, and carboxylic acid functional groups. Ligands were attached or incorporated into a wide variety of solid supports. In most cases, attachment was via covalent bonding to organic resins, polymers, or silica-based supports. Most of the REE projects were conducted at a low technology readiness level (TRL) and showed promise, but direct comparison between technologies was not possible based on the available information. It is recommended that testing and reporting be standardized to the extent possible to facilitate comparisons between technologies. Two projects were directed at novel lithium extraction technology. Both projects investigated the use of inorganic sorbents, including manganese oxides. One study also examined the use of metal- ion imprinted polymers as selective ion-exchange resins for the separation of lithium and manganese from brines. Both approaches showed promise for the selective extraction of lithium from brines, including potentially geothermal brines. Results from these GTO studies indicated that selective REE and lithium extraction is possible, but interference from co-occurring solutes, such as calcium, magnesium, or heavy metals, will interfere with process efficiency and negatively impact process economics. Techno-economic analysis conducted as part of the resource and technology studies suggest extraction of REE from geothermal brines is unlikely to be economically viable, especially since non-geothermal produced waters frequently have higher REE concentrations. It is recommended that benchmarks for techno-economic analysis be established to the extent possible for future studies, to facilitate direct comparison of various technologies. Based on the collective results of this program, it appears that hybrid geothermal power would benefit more from recovery of lithium and other metals, rather than REE. It is recommended that future studies be conducted at a higher- TRL and that sorbents be tested against actual geothermal fluid samples. Prior higher-TRL efforts to extract metals from geothermal brines should be further evaluated for lessons learned.

36 MATERIALS SCIENCE↗

Distributed fiber sensing systems for 3D combustion temperature field monitoring in coal-fired boilers using optically generated acoustic waves (Final Report)

In this project, we have developed and tested three kinds of fiber optic sensing systems for real time monitoring of temperature variations within an industrial scale boiler furnace. The fiber optic sensing systems target spatial and temporal distributions of high temperature profiles in a boiler furnace in fossil power plants. The reconstructed temperature profile will provide critical input for the control mechanisms to optimize the combustion process. This temperature profile will address the essential problem for fossil power plants in achieving higher efficiency and fewer pollutant emissions. Acoustic pyrometer systems have been used to reconstruct temperature field of power plant boilers based on measuring TOF (times-of-flight) of sound waves along some straight paths in a 2D cross-section of the boiler. In this project, optically generated acoustic signals from a fiber optic sensing system have replaced the acoustic signals generated from an electrical transducer. A 3D reconstruction algorithm replaced the previous 2D model. In this project, three kinds of fiber optic sensing systems have been developed and tested. They are fiber optic sensing system I, fiber optic sensing system II (Distributed Sensing System I) and fiber optic sensing system III (Distributed Sensing System II). For fiber optic sensing system I, the fiber optic ultrasound generator acts as a signal generator. A microphone, hydrophone or other electronic devices serve as a signal receiver. In this system, there are one generator and one receiver. Distance test, water temperature test, air temperature test, air temperature reconstruction, and GE ISBF pilot test were performed by Fiber optic sensing system I. The fiber optic sensing system I successfully detected temperature in all these tests. We got 2D temperature reconstruction results by using the fiber optic sensing system I and it matched the reference data. The fiber optic sensing system I successfully survived in GE ISBF boiler environment (480 °F). For fiber optic sensing system II (Distributed Sensing System I), it is an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Fiber Bragg Grating (FBG) and Fabry-Perot (FP) sensor act as a signal receiver. In this system, there is one generator and one receiver. Aluminum plate temperature test, furnace high temperature test, and GE ISBF pilot test were performed by the fiber optic sensing system II. Fiber optic sensing system II successfully detected the temperature in all these tests. The fiber optic sensing system II successfully survived at up to 700 °C furnace environment and 320 °C GE ISBF boiler environment. For fiber optic sensing system III (Distributed Sensing System II), it is also an all optical ultrasound system. The fiber optic ultrasound generator acts as a signal generator. Multiple FP fiber sensors act as signal receivers. In this system, there is one generator and three receivers. Three GE ISBF pilot tests were performed by fiber optic sensing system III. The fiber optic sensing system III survived in the cold flow tests in GE’s ISBF pilot test facility. However, we didn’t get high temperature data by using this system since the nanosecond laser issues. During the period of the project, test trials, data simulation and algorithm optimization was performed successfully. For real time temperature field construction, the sampling rate must be fast enough to capture the field variations. The technology of Code-division multiple access (CDMA) is well studied which could allow parallel multiplexing, even if signals overlap in time or frequencies. Moreover, it has been known that extending the length of signal significantly improves SNR. For acoustic signals, these multiplexing techniques have also been widely used, mainly for sonar and acoustic communications. The CDMA modulation technique has been proposed and studied to guarantee high network throughput, low channel access delay and low energy consumption. We have studied the temperature field reconstruction using Gaussian Radial Basis Functions (GRBF)-based approximation approach. Reconstruction of 3D temperature field using Neural Networks with measured TOF and known propagation paths is feasible. 2D and 3D temperature field reconstruction simulation results are achieved. The milestone status is shown in Table 1. We finished milestone 1-8 and milestone 10. For milestone 9, we did three pilot tests by using the fiber optic sensing system III (Distributed Sensing System II) at GE Power. However, due to the failure of the ns laser, we did not get the temperature results. We conducted some additional tasks that were not originally proposed: 1) We fabricated a fiber optic sensing system I and did a pilot test based on this system. 2) In the proposal, we proposed two pilot tests at GE Power. In reality, we finished at least seven pilot tests at GE Power. GE Power has made a lot of efforts for supporting the pilot tests. 3) We got a simulation results based on CDMA. In summary, most of the tasks have been accomplished. The outcome of this project removed a few barriers that hinder the achievement of the final product of the distributed sensing systems. With the successful accomplishment of this project, a prototype of the fiber optic sensing system can be fabricated to attract more interests from companies and other funding agencies.

47 OTHER INSTRUMENTATION↗

NASA's Habitation Development Status: Current Concepts and ISRU Opportunities

Introduction: The National Aeronautics and Space Administration (NASA) is embarking on a bold journey to return humankind to the Moon and onward to Mars with innovative commercial, international, and academic partnerships [1]. Under the Artemis series of missions, NASA seeks to establish sustained human exploration of deep space through an objectives-based approach [2]. This approach drives the identification of needed system functionality and the current and future capabilities which will eventually allow humanity to sustainably live beyond Earth. Providing evolvable and scalable habitation is a cornerstone function that calls for the collection and integration of current, developmental, and future technologies that can meet near-term exploration needs while growing into long-term sustained presence. NASA is advancing in-space habitation through its Next Space Technologies for Exploration Partnerships (NextSTEP) model while designing lunar, Mars transit, and Mars surface habitat government reference concepts for Artemis missions. These efforts have unveiled possible near-term opportunities for the in-space resource utilization community if human habitation is considered a future customer of space resources. NextSTEP Habitation Development: NASA is closely working with commercial partners under its NextSTEP Appendix A model to advance habitation systems in the arena of inflatable and composite habitation structures. Such efforts promise efficiencies in volumetric packaging and overall spacecraft mass respectively. Recent testing by commercial partners have helped to quantify possible failure mechanisms for inflatable structures while advancing their technology towards eventual flight certification. The advancement of such Class II habitation structures, in which the habitat is only fully deployed once in-space or on a planetary surface [3], is critical to providing increased habitable volume for long-duration missions with no additional mass penalties. The progression of such technology is infused into NASA’s government reference concepts for notional deep space habitation concepts. Current Government Reference Concepts: To best inform the formulation of future collaborative solicitations, NASA employs the practice of internally developing reference concepts for future exploration elements. These concepts aid in identifying the functions and capabilities needed to complete NASA missions as well as feasible solutions within the timeframe needed. Government reference concepts for a lunar Surface Habitat, Mars Transit Habitat, and Mars Sur-face Habitat are continuously being developed and updated to better guide the Agency’s overall exploration architecture. Lunar Surface Habitat. As NASA returns to the Moon, it is evaluating possible lunar surface habitation concepts. The Surface Habitat (SH) reference concept entails a hybrid metallic-inflatable structure capable of initially housing two crew for surface stays of up to 30 days in duration [4]. While initial missions may span ~7 days in duration, consideration is being given to expanding SH’s capability to support a crew of four for up to 60 days over its 15-year design life [5]. Functionally, the SH serves as a ‘hub’ for all Artemis crewed surface operations, providing internal volume for maintenance, medical, logistics, science utilization, and extravehicular activity (EVA) support in addition to core habitation functionality such as environmental control and life support (ECLS) and power generation and distribution among many others. Additionally, NASA has entered a study agreement with the Italian Space Agency (Agenzia Spaziale Italiana – ASI) to investigation a possible Multi-Purpose Habitat (MPH) as an additional or augmenting habit-able element for the lunar surface [6]. Mars Transit Habitat. NASA’s current architectural concept for initial human missions to Mars entails the utilization of a transit habitat (TH) to transport a four-person crew to and from Mars orbit, departing from and returning to a lunar near-rectilinear halo orbit (NRHO), over the course of a ~1,200-day mission [5]. While holding a similar 15-year design lifetime, TH will also support a series of analog mission activities in NRHO to gradually test the systems and interaction with lunar surface elements, some of which may be adapted for Mars surface exploration. Holding similar functional capabilities as SH, TH is sized to support much longer durations in space and greater logistical independence. Mars Surface Habitat. NASA is still exploring the concept of operations for initial crewed missions to Mars. As such, the Mars Surface Habitat (MSH) concept is still in its infancy as options for mobile, pressurized habitation and stationary habitats are being explored. It is expected MSH will leverage heavily from the lunar SH and possible lunar pressurized rover, however the very different Martian environment will likely necessitate modifications. ISRU Opportunities: Despite the advancements under NASA’s NextSTEP habitation work and continually optimized reference concepts, NASA is facing near-term mass and power challenges that may create opportunities to the ISRU community by providing yet another possible customer for space resources. While NASA desires to use regenerative ECLS systems (ECLSS) for all habitation concepts, their operation comes with initial mass penalties and maintenance overheads when compared to simpler open-loop architectures. Because of this, NASA is currently pro-posing an open-loop, consumables-based architecture for its surface habitats to achieve initial launch and delivery lander mass targets while scarring for the in-corporation of regenerative ECLSS to meet longer mission durations and sustained presence. With such an architecture, oxygen and potable water are needed consumables, which initial ISRU systems may be able to provide in a pilot capacity. Although the TH is expecting to utilize regenerative ECLSS, advancements in lunar surface-based ISRU and possible re-supply of spacecraft in lunar or Mars orbit could significantly reduce the logistical need for ECLSS related spares on TH and possibly allow a similar open-loop and consumables-based architecture. Table 1 Notional ECLSS consumables per mission use for SH and TH [5]. SH (28-day) Open Loop: Water (kg) 288 | Oxygen (kg) 30 TH (1110-day) Closed-Loop: Water (kg) 182 | Oxygen (kg) 123 In addition to ECLSS mass concerns, the SH is facing challenges with energy storage to support operations over periods of darkness exceeding 100 hours. Both battery and fuel cell-based power architectures are being traded, opening options for external power generation and energy storage. One possibility is ISRU-produced H2 and O2 feeding external primary fuel cell systems which could supplement habitation power generation while reducing initial mass and volume until much more powerful fission power systems might be deployed. Conclusion: Significant advancements are being made in evolvable habitation concepts that span from near-term technologies, such as inflatable structures under NextSTEP, to potentially revolutionary capabilities like the lunar surface construction as funded through the Moon-to-Mars Planetary Autonomous Construction Technologies (MMPACT) project [7]. As NASA investigates both heritage capabilities and rap-idly advancing, disruptive technologies, there are likely many near-term opportunities for initial ISRU capabilities to significantly aid human habitation and increase self-sufficiency beyond Earth.

habitation↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗

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↗