Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed source”

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

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

At least 1,063 records · Page 59

Flight Testing of a Cryogenic Capillary Pumped Loop

Future space-based cryogenic systems will require enhanced integration flexibility, lower weight reduced parasitic penalties, better vibration isolation, and a variety of other improvements to meet performance goals. Additionally, there is an increasing need to locate cooling sources remotely from cooled components. In the past flexible conductive links were used and worked well in most cases. However, as the transport lengths increase, conductive couplings become heavier and less effective, and must be replaced by higher performance systems. One available option, which can meet many of these future requirements, is the cryogenic capillary pumped loop (CCPL). The development of the CCPL technology started in 1992, following the success of the room temperature CPLS. The extrapolation of CCPL technology to cryogenic temperatures offers many performance benefits, which are not currently within the reach of traditional heat pipes or conductive links. Specific advantages of the CCPL technology pertaining to cryocooler integration include: (1) greater capillary pumping pressure for improved ground testability; (2) improved mechanical isolation; (3) faster diode shutdown and lower reverse heat leaks; (4) tighter control of detector temperature; (5) variable or fixed conductance operation; and (6) ease of integration due to their flexibility. The applications of CCPL technology are numerous. Military and commercial applications include surveillance satellites, earth observing satellites, deep space observation systems, medical devices, and many other cryogenic systems. Over the past few years, several breadboard and prototype CCPLs have been built and ground tested. A prototype CCPL has demonstrated successful operation between 80K and 110K with heat loads between O.5W and 12W using nitrogen as the working fluid, and 35K and 40K with head loads of 0.25W to 3.5W using neon. In order to verify CCPL performance in a microgravity environment, a flight unit, CCPL-5, was tested onboard the Space Shuttle STS-95 in October of 1998 as part of the CRYOTSU Flight Experiment. This flight was the first in-space demonstration of the CCPL. The CCPL-5 utilized nitrogen as the working fluid and operated between 75K and 110K. Flight results indicated excellent performance of the CCPL-5 under zero-G environment The CCPL could start from a supercritical condition in all tests, and the loop operating temperature could be tightly controlled regardless of changes in the heat load and/or the sink temperature. In addition, the loop demonstrated successful operation with a heat load of 0.5W as well as with parasitic heat loads alone. There were no noticeable differences between zero-G and one-G operation.

Ku, Jentung↗

A Case Study Comparing Citizen Science Aurora Data with Global Auroral Boundaries Derived from Satellite Imagery and Empirical Models

Aurorasaurus is a citizen science project that offers a new, global data source consisting of ground-based reports of the aurora. For this case study, aurora data collected during the 17-18 March 2015 geomagnetic storm are examined to identify their conjunctions with Defense Meteorological Satellite Program (DMSP) satellite passes over the high latitude auroral regions. This unique set of aurora data can provide ground-truth validation of existing auroral precipitation models. Particularly, the solar wind driven, Oval Variation, Assessment, Tracking, Intensity, and Online Nowcasting (OVATION) Prime 2013 (OP-13) model and a Kp-dependent model of Zhang-Paxton (Z-P) are utilized for our boundary validation efforts. These two similar models are compared for the first time. Global equatorward auroral boundaries are derived from the OP 13 model and the DMSP Special Sensor Ultraviolet Spectrographic Imager (SSUSI) far ultraviolet (FUV) data using the Z-P model at a fixed flux level of 0.2 erg cm(exp -2)s(exp -1). These boundaries are then compared with citizen science reports as well as with each other. Even though there are some large differences between the global boundaries for a few cases, the average difference is about 1.5 deg in geomagnetic latitude, with OP-13 being equatorward of Z-P model. When these boundaries are compared with each other as a function of local time, no clear overall trend as a function of local time was observed. It is also found that the ground based reports are more consistent with the predictions of the OP-13 model.

Kosar, Burcu C.↗

A Computational and Experimental Study of Resonators in Three Dimensions

In a previous work by the present authors, a computational and experimental investigation of the acoustic properties of two-dimensional slit resonators was carried out. The present paper reports the results of a study extending the previous work to three dimensions. This investigation has two basic objectives. The first is to validate the computed results from direct numerical simulations of the flow and acoustic fields of slit resonators in three dimensions by comparing with experimental measurements in a normal incidence impedance tube. The second objective is to study the flow physics of resonant liners responsible for sound wave dissipation. Extensive comparisons are provided between computed and measured acoustic liner properties with both discrete frequency and broadband sound sources. Good agreements are found over a wide range of frequencies and sound pressure levels. Direct numerical simulation confirms the previous finding in two dimensions that vortex shedding is the dominant dissipation mechanism at high sound pressure intensity. However, it is observed that the behavior of the shed vortices in three dimensions is quite different from those of two dimensions. In three dimensions, the shed vortices tend to evolve into ring (circular in plan form) vortices, even though the slit resonator opening from which the vortices are shed has an aspect ratio of 2.5. Under the excitation of discrete frequency sound, the shed vortices align themselves into two regularly spaced vortex trains moving away from the resonator opening in opposite directions. This is different from the chaotic shedding of vortices found in two-dimensional simulations. The effect of slit aspect ratio at a fixed porosity is briefly studied. For the range of liners considered in this investigation, it is found that the absorption coefficient of a liner increases when the open area of the single slit is subdivided into multiple, smaller slits.

Tam, C. K. W.↗

Improved Detection System Description and New Method for Accurate Calibration of Micro-Channel Plate Based Instruments and Its Use in the Fast Plasma Investigation on NASA's Magnetospheric MultiScale Mission

The Fast Plasma Investigation (FPI) on NASAs Magnetospheric MultiScale (MMS) mission employs 16 Dual Electron Spectrometers (DESs) and 16 Dual Ion Spectrometers (DISs) with 4 of each type on each of 4 spacecraft to enable fast (30 ms for electrons; 150 ms for ions) and spatially differentiated measurements of the full 3D particle velocity distributions. This approach presents a new and challenging aspect to the calibration and operation of these instruments on ground and in flight. The response uniformity, the reliability of their calibration and the approach to handling any temporal evolution of these calibrated characteristics all assume enhanced importance in this application, where we attempt to understand the meaning of particle distributions within the ion and electron diffusion regions of magnetically reconnecting plasmas. Traditionally, the micro-channel plate (MCP) based detection systems for electrostatic particle spectrometers have been calibrated using the plateau curve technique. In this, a fixed detection threshold is set. The detection system count rate is then measured as a function of MCP voltage to determine the MCP voltage that ensures the count rate has reached a constant value independent of further variation in the MCP voltage. This is achieved when most of the MCP pulse height distribution (PHD) is located at higher values (larger pulses) than the detection system discrimination threshold. This method is adequate in single-channel detection systems and in multi-channel detection systems with very low crosstalk between channels. However, in dense multi-channel systems, it can be inadequate. Furthermore, it fails to fully describe the behavior of the detection system and individually characterize each of its fundamental parameters. To improve this situation, we have developed a detailed phenomenological description of the detection system, its behavior and its signal, crosstalk and noise sources. Based on this, we have devised a new detection system calibration method that enables accurate and repeatable measurement and calibration of MCP gain, MCP efficiency, signal loss due to variation in gain and efficiency, crosstalk from effects both above and below the MCP, noise margin, and stability margin in one single measurement. More precise calibration is highly desirable as the instruments will produce higher quality raw data that will require less post-acquisition data correction using results from in-flight pitch angle distribution measurements and ground calibration measurements. The detection system description and the fundamental concepts of this new calibration method, named threshold scan, will be presented. It will be shown how to derive all the individual detection system parameters and how to choose the optimum detection system operating point. This new method has been successfully applied to achieve a highly accurate calibration of the DESs and DISs of the MMS mission. The practical application of the method will be presented together with the achieved calibration results and their significance. Finally, it will be shown that, with further detailed modeling, this method can be extended for use in flight to achieve and maintain a highly accurate detection system calibration across a large number of instruments during the mission.

detection system↗

Pipe Parity Analysis of Seawater Desalination in the United States: Exploring Costs, Energy, and Reliability via Case Studies and Scenarios of Emerging Technology

As climate change-induced variables exacerbate water scarcities, the use of seawater reverse osmosis (SWRO) membrane desalination technology for treating seawater could potentially provide a long-term drought-proof source of drinking water. This study carried out a technoeconomic assessment (TEA) for three SWRO desalination plants in the U.S. and one in Israel to set baselines for the cost and energy consumption for seawater desalination. In addition, a breakeven curve for implementing SWRO desalination was estimated in relation to the cost of water conservation measures to meet drought-induced reduction of the traditional water supply. The results show that the cost of SWRO water production scales with the plant capacity while energy intensity is not dramatically different across the plants. The higher cost in some U.S. plants is due to high capital investment, including land acquisition and permitting. Variations in plant capacity utilization have the greatest impact on the levelized cost of water (LCOW) over the plant service life, suggesting the importance of reducing fouling and maintenance-related downtime. Scenario analysis of fixed labor cost reduction through process automation indicates investments in automation and sensing technology could result in long-term savings. Breakeven analysis shows the decision to adopt SWRO is highly dependent on the local cost associated with water conservation to meet water supply reduction. Moreover, a small reduction in SWRO cost can influence a shift toward the adoption of SWRO over water conservation measures. Incorporation of future water demand, water conservation potential, and water stress data around the nation indicates SWRO desalination could be an important contributor to the future municipal drinking water portfolio in the U.S.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) research platform. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗

Autonomous Cryogenic Load Operations: Knowledge-Based Autonomous Test Engineer

The Knowledge-Based Autonomous Test Engineer (KATE) program has a long history at KSC. Now a part of the Autonomous Cryogenic Load Operations (ACLO) mission, this software system has been sporadically developed over the past 20 years. Originally designed to provide health and status monitoring for a simple water-based fluid system, it was proven to be a capable autonomous test engineer for determining sources of failure in the system. As part of a new goal to provide this same anomaly-detection capability for a complicated cryogenic fluid system, software engineers, physicists, interns and KATE experts are working to upgrade the software capabilities and graphical user interface. Much progress was made during this effort to improve KATE. A display of the entire cryogenic system's graph, with nodes for components and edges for their connections, was added to the KATE software. A searching functionality was added to the new graph display, so that users could easily center their screen on specific components. The GUI was also modified so that it displayed information relevant to the new project goals. In addition, work began on adding new pneumatic and electronic subsystems into the KATE knowledge base, so that it could provide health and status monitoring for those systems. Finally, many fixes for bugs, memory leaks, and memory errors were implemented and the system was moved into a state in which it could be presented to stakeholders. Overall, the KATE system was improved and necessary additional features were added so that a presentation of the program and its functionality in the next few months would be a success.

Schrading, J. Nicolas↗

Modeling of Combined Lead Fast Reactor and Concentrating Solar Power Supercritical Carbon Dioxide Cycles to Demonstrate Feasibility, Efficiency Gains, and Cost Reductions

Solar power has innate issues with weather, grid demand and time of day, which can be mitigated through use of thermal energy storage for concentrating solar power (CSP). Nuclear reactors, including lead-cooled fast reactors (LFRs), can adjust power output according to demand; but with high fixed costs and low operating costs, there may not be sufficient economic incentive to make this worthwhile. We investigate potential synergies through coupling CSP and LFR together in a single supercritical CO 2 Brayton cycle and/or using the same thermal energy storage. Combining these cycles allows for the LFR to thermally charge the salt storage in the CSP cycle during low-demand periods to be dispatched when grid demand increases. The LFR/CSP coupling into one cycle is modeled to find the preferred location of the LFR heat exchanger, CSP heat exchanger, sCO 2 -to-salt heat exchanger (C2S), turbines, and recuperators within the supercritical CO 2 Brayton cycle. Three cycle configurations have been studied: two-cycle configuration, which uses CSP and LFR heat for dedicated turbocompressors, has the highest efficiencies but with less component synergies; a combined cycle with CSP and LFR heat sources in parallel is the simplest with the lowest efficiencies; and a combined cycle with separate high-temperature recuperators for both the CSP and LFR is a compromise between efficiency and component synergies. Additionally, four thermal energy storage charging techniques are studied: the turbine positioned before C2S, requiring a high LFR outlet temperature for viability; the turbine after the C2S, reducing turbine inlet temperature and therefore power; the turbine parallel to the C2S producing moderate efficiency; and a dedicated circulator loop. While all configurations have pros and cons, use of a single cycle offers component synergies with limited efficiency penalty. Using a turbine in parallel with the C2S heat exchanger is feasible but results in a low charging efficiency, while a dedicated circulator loop offers flexibility and near-perfect heat storage efficiency but increasing cost with additional cycle components.

14 SOLAR ENERGY↗

On the Unusually High Temperature of the Cluster of Galaxies 1E 0657-56

A recent X-ray observation of the cluster 1E 0657-56 (z = 0.296) with ASC,4 implied an unusually high temperature of approx. 17 keV. Such a high temperature would make it the hottest known cluster and severely constrain cosmological models since, in a Universe with critical density (Omega = 1) the probability of observing such a cluster is only approx. 4 x 10(exp -5). Here we test the robustness of this observational result since it has such important implications. We analysed the data using a variety of different data analysis methods and spectral analysis assumptions and find a temperature of approx. 11 - 12 keV in all cases, except for one class of spectral fits. These are fits in which the absorbing column density is fixed at the Galactic value. Using simulated data for a 12 keV cluster, we show that a high temperature of approx. 17 keV is artificially obtained if the true spectrum has a stronger low-energy cut-off than that for Galactic absorption only. The apparent extra absorption may be astrophysical in origin, (either intrinsic or line-of-sight), or it may be a problem with the low-energy CCD efficiency. Although significantly lower than previous measurements, this temperature of kT approx. 11 - 12 keV is still relatively high since only a few clusters have been found to have temperatures higher than 10 keV and the data therefore still present some difficulty for an Omega = 1 Universe. Our results will also be useful to anyone who wants to estimate the systematic errors involved in different methods of background subtraction of ASCA data for sources with similar signal-to-noise to that of the IE 0657-56 data reported here.

Yaqoob, Tahir↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗

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↗

Beamed Energy and Communications Optical Node (BEACON) Demonstrator

Due to long shadow periods (2 weeks or greater) on the south pole, concepts to raise solar arrays to a sufficient height at specific locations have shown the capability to provide power to surface users for much longer periods. Such a tower can also provide a 3rd Generation Partnership Project (3GPP) service for users up to 10 km away, dependent on terrain. An option to deliver power, albeit with low efficiency, using a laser beam coupled with the tower height could provide mobile and fixed users power during darkness, reducing their battery requirements. A demonstration of these technologies in the lunar environment is crucial to support future Artemis campaigns as well as potential emerging lunar infrastructures. A demonstrator design of a 15 m deployed boomon the south pole has been shown to enable both the gathering of kilowatts of power and the provision of 3GPP relay and backhaul given an appropriate lunar location. Using a laser to send power to a photovoltaic receiver has been proposed to transmit electrical power on the moon, particularly for applications such as powering a rover in near-polar permanently shadowed regions (PSR) where solar power is not available. In this work, the Compass team performed a conceptual engineering design study of a near-term laser surface-to-surface power beaming and relay station using a tower to simultaneously carry the power source (Vertical Solar Array Technologies (VSAT)[1]), the 3GPP relay antenna, and the laser telescope.

Deployable solar array tower laser beamed power 3G↗

Beamed Energy and Communications Optical Node (BEACON) Demonstrator

Due to long shadow periods (2 weeks or greater) on the south pole, concepts to raise solar arrays to a sufficient height at specific locations have shown the capability to provide power to surface users for much longer periods. Such a tower can also provide a 3rd Generation Partnership Project (3GPP) service for users up to 10 km away, dependent on terrain. An option to deliver power, albeit with low efficiency, using a laser beam coupled with the tower height could provide mobile and fixed users power during darkness, reducing their battery requirements. A demonstration of these technologies in the lunar environment is crucial to support future Artemis campaigns as well as potential emerging lunar infrastructures. A demonstrator design of a 15 m deployed boomon the south pole has been shown to enable both the gathering of kilowatts of power and the provision of 3GPP relay and backhaul given an appropriate lunar location. Using a laser to send power to a photovoltaic receiver has been proposed to transmit electrical power on the moon, particularly for applications such as powering a rover in near-polar permanently shadowed regions (PSR) where solar power is not available. In this work, the Compass team performed a conceptual engineering design study of a near-term laser surface-to-surface power beaming and relay station using a tower to simultaneously carry the power source (Vertical Solar Array Technologies (VSAT)[1]), the 3GPP relay antenna, and the laser telescope.

Lunar Relay↗

Machine Learning assisted optimization and parameter space exploration dataset of spin ice Hamiltonian

This repository contains both simulated and experimental structure factor data for the data challenge involving the inverse scattering problem. The simulated data were generated during a machine-learning-assisted optimization routine described in ref[1]. The experimental structure factor was measured on a rare-earth oxide, Dy2Ti2O7 using diffuse neutron scattering from time-of-flight techniques on the CORELLI instrument at the Spallation Neutron Source, Oak Ridge National Laboratory. A Metropolis Monte Carlo code implemented to run in a High-performance computing setting was used to calculated simulated structure factors for the spin-ice Hamiltonian at 680 mK, which is the same temperature as for the experimental data. The total size of all the files in this repository is 5.12 GB. A detailed description of the files is given below. ExperimentalData_630mK.dat – A linearized version of 3-dimensional experimental data of size 61×81×21. This data was processed to remove an estimation of non-magnetic background, including nuclear scattering signal and instrumentation background. Parameters.dat – 6700 samples were evaluated over the 4-dimensional parameter space (J_1, J_2, J_3 and J_(3^' )). There is an additional parameter, D in the spin Hamiltonian to mimic the dipolar interaction between magnetic ions. However, this parameter, D was fixed to a value determined by prior work. This file contains five columns for the parameters J_1, J_2, J_3, J_(3^' ) and D respectively. 3D_Simulation_Data.dat – The simulated structure factor, S(Q) data are included in this file. Each raw contains a linearized array of 3D volumes of S(Q) calculated for the parameter set given in the corresponding row of the file Parameters.dat. The size of the volume data was matched to the experimental data. Qx(h,-h,0).dat, Qy(k,k,-2k).dat, Qz(l,l,l).dat – These files contain the h, k, and l values along with the reciprocal vectors [h,-h,0], [k,k,-2k] and [l,l,l] respectively.

36 MATERIALS SCIENCE↗

Monitoring Global Sea Level: Eustatic Variations, Local Variations, and Solid Earth Effects

Project BIFROST (Baseline Inferences for Fennoscandian Rebound Observations. Sea-level and Tectonics) combines networks of continuously operating CPS receivers in Sweden and Finland to measure ongoing crustal deformation due to glacial isostatic adjustment, (CIA). We present an analysis of data collected in the years 1993-1998. We compare the CPS determinations of three-dimensional crustal motion to predictions calculated using the high resolution Fennoscandian deglaciation model recently proposed by Lambeck et al. We find that the the maximum observed uplift rate (approx. 10 mm/ yr) and the maximum predicted uplift rate agree to better than 1 mm/ yr. The patterns of uplift also agree quite well, although differences are discernible. The chi(exp 2) difference between predicted and GPS-observed radial rates is reduced by a factor of 5-6 compared to that for the "null" (no uplift) model, depending on whether a mean difference is first removed. The north components of velocity agree at about the same relative level. whereas the agreement for the east components is worse, a problem possibly related to the lack of bias fixing. We have also compared the values for the observed radial deformation rates to those based on sea-level rates from Baltic tide gauges. The weighted RMS difference between CPS and tide-gauge rates (after removing a mean) is 0.6 mm/ yr, giving an indication of the combined accuracy of the CPS and tide-gauge measurement systems. Spectral analysis of the time series of position estimates yields spectral indices in the range -1 to -2. An EOF analysis indicates, however, that much of this power is correlated among the sites. The correlation appears to be regional and falls off only slightly with distance. Some of this correlated noise is associated with snow accumulation on the antennas or, for those antennas with radomes, on the radomes. This problem has caused us to modify the radomes used several times, leading to one of our more significant sources of uncertainty.

Davis, James L.↗

Aerosol Midlatitude Cyclone Indirect Effects in Observations and High-Resolution Simulations

Aerosol-cloud interactions are a major source of uncertainty in inferring the climate sensitivity from the observational record of temperature. The adjustment of clouds to aerosol is a poorly constrained aspect of these aerosol-cloud interactions. Here, we examine the response of midlatitude cyclone cloud properties to a change in cloud droplet number concentration (CDNC). Idealized experiments in high-resolution, convection-permitting global aquaplanet simulations with constant CDNC are compared to 13 years of remote-sensing observations. Observations and idealized aquaplanet simulations agree that increased warm conveyor belt (WCB) moisture flux into cyclones is consistent with higher cyclone liquid water path (CLWP). When CDNC is increased a larger LWP is needed to give the same rain rate. The LWP adjusts to allow the rain rate to be equal to the moisture flux into the cyclone along the WCB. This results in an increased CLWP for higher CDNC at a fixed WCB moisture flux in both observations and simulations. If observed cyclones in the top and bottom tercile of CDNC are contrasted it is found that they have not only higher CLWP but also cloud cover and albedo. The difference in cyclone albedo between the cyclones in the top and bottom third of CDNC is observed by CERES to be between 0.018 and 0.032, which is consistent with a 4.6-8.3 Wm(exp -2) in-cyclone enhancement in upwelling shortwave when scaled by annual-mean insolation. Based on a regression model to observed cyclone properties, roughly 60% of the observed variability in CLWP can be explained by CDNC and WCB moisture flux.

Mccoy, Daniel T.↗

Strong Sensitivity of the Isotopic Composition of Methane to the Plausible Range of Tropospheric Chlorine

The C-13 isotopic ratio of methane, δC-13 of CH4, provides additional constraints on the CH4 budget to complement the constraints from CH4 observations. The interpretation of δC-13 observations is complicated, however, by uncertainties in the methane sink. The reaction of CH4 with Cl is highly fractionating, increasing the relative abundance of 13CH4, but there is currently no consensus on the strength of the tropospheric Cl sink. Global model simulations of halogen chemistry differ strongly from one another in terms of both the magnitude of tropospheric Cl and its geographic distribution. This study explores the impact of the inter-model diversity in Cl fields on the simulated δC-13 of CH4. We use a set of GEOS global model simulations with different predicted Cl fields to test the sensitivity of the δC-13 of CH4 to the diversity of Cl output from chemical transport models. We find that δC-13 is highly sensitive to both the amount and geographic distribution of Cl. Simulations with Cl providing 0.28% or 0.66% of the total CH4 loss bracket the δC-13 observations for a fixed set of emissions. Thus, even when Cl provides only a small fraction of the total CH4 loss and has a small impact on total CH4, it provides a strong lever on δC-13. Consequently, it is possible to achieve a good representation of total CH4 using widely different Cl concentrations, but the partitioning of CH4 loss between the OH and Cl reactions leads to strong differences in isotopic composition depending on which model’s Cl field is used. Comparing multiple simulations, we find that altering the tropospheric Cl field leads to approximately a 0.5‰ increase in δ(13)CH4for each percent increase in how much CH4 is oxidized by Cl. The geographic distribution and seasonal cycle of Cl also impacts the hemispheric gradient and seasonal cycle of δC-13. The large effect of Cl on δC-13 compared to total CH4 broadens the range of CH4 source mixtures that can be reconciled with δC-13 observations. Stronger constraints on tropospheric Cl are necessary to improve estimates of CH4 sources from δC-13 observations.

Global model simulations of halogen chemistry↗

Experimental Simulation of Active Control With On-line System Identification on Sound Transmission Through an Elastic Plate

An adaptive control algorithm with on-line system identification capability has been developed. One of the great advantages of this scheme is that an additional system identification mechanism such as an additional uncorrelated random signal generator as the source of system identification is not required. A time-varying plate-cavity system is used to demonstrate the control performance of this algorithm. The time-varying system consists of a stainless-steel plate which is bolted down on a rigid cavity opening where the cavity depth was changed with respect to time. For a given externally located harmonic sound excitation, the system identification and the control are simultaneously executed to minimize the transmitted sound in the cavity. The control performance of the algorithm is examined for two cases. First, all the water was drained, the external disturbance frequency is swept with 1 Hz/sec. The result shows an excellent frequency tracking capability with cavity internal sound suppression of 40 dB. For the second case, the water level is initially empty and then raised to 3/20 full in 60 seconds while the external sound excitation is fixed with a frequency. Hence, the cavity resonant frequency decreases and passes the external sound excitation frequency. The algorithm shows 40 dB transmitted noise suppression without compromising the system identification tracking capability.

Source record↗