Search NASA⌕ Search

SEARCH · Search NASA

Results for “sample average approximation”

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.

464 records · Page 26

ROSAT observations of Coma Cluster galaxies

The approximately 86 ks ROSAT Position Sensitive Proportional Counter (PSPC) image of the Coma Cluster is deeper than any previous X-ray observation of a galaxy cluster. We search for X-ray emission from 35 individual galaxies in a magnitude-limited sample, all of which lie within 20 arcmins of the optical axis in at least one of the four Coma pointings. We detect seven galaxies in the 0.4-2.4 keV band at a significance level exceeding 3 sigma, and a further four at above 2 sigma. Although we can set only upper limits on the individual flux from each of the other galaxies, we are able to measure their mean flux by stacking the observations. The X-ray luminosities of the seven detections range from 6.2 x 10(exp 40) to 1.5 x 10(exp 42) ergs/s (0.4-2.4 keV for H(sub 0) = 50 km/s/Mpc). For galaxies with a blue absolute magnitude of about -21 we find a mean X-ray luminosity of 1.3 x 10(exp 40) ergs/s. The ratio of X-ray to optical luminosity is substantially smaller for such subjects than for the brightest galaxies in the cluster. The X-ray luminosities of the four brightest galaxies are ill-defined, however, because of ambiguity in distinguishing galaxy emission from cluster emission. Each object appears to be related to significant structure in the diffuse intracluster medium. We also investigate emission in the softer 0.2-0.4 keV band where detections are less significant because of the higher background, and we discuss the properties of a number of interesting individual sources. The X-ray luminosities of the Coma galaxies are similar to those of galaxies in the Virgo Cluster and in other regions with relatively low galaxy density. We conclude that large-scale environmental effects do not significantly enhance or suppress the average X-ray emission from galaxies, but that individual objects vary in luminosity substantially in a way which may depend on the detailed history of their environment.

Dow, K. L.↗

Recession of FEP specimens from trays D11 and B7

We report work done at Boeing Defense and Space Group on analysis of silvered teflon specimens taken from selected locations of the Long Duration Exposure Facility under support from a contract provided by NASA LaRC. The samples discussed in this presentation were taken from the unexposed side of D11 and extended through the folded area of this blanket into the exposed area. Two similar areas were cut from blanket B7, one from the edge of the blanket near row six and one from the edge of the blanket near row eight and within a few centimeters of the copper grounding strap for B7. The specimens were each divided into three sections by cutting with a scapel. Two of the sections were mounted in a potting compound, which was cut and polished such that the cross-sectional thickness of each was exposed. One piece was mounted straight and the other was mounted in an attempt to configure the specimens such that it was bent with a radius of curvature similar to the on-orbit configuration. The third portion of each specimen was used for SEM images to help define the angle of exposure with respect to the ram at each location on the specimen. Photomicrographs were taken in cross section from the edge of the blanket through the curved transition region into the exposed area of the blanket. The thickness of the Fluorinated Ethylene Propylene (FEP) layer was determined at known distances from the edge of the blanket. SEM images were obtained at known distances to help define the angle with respect to ram and therefore establish the atomic oxygen fluence on each location and correlate this exposure with thickness. Thickness measurements made with a two to three centimeter distance minimized the uncertainty arising from variations in the as-manufactured thickness of each blanket. The nominal angle from ram of the exposed portion of each blanket, and the fact that the unexposed edge portions are approximately at right angles to the exposed portion were also used to help define the angles. Thickness measurements were taken at specified locations. An average thickness for the unexposed portion of the blanket was determined. Changes in thickness were then determined by difference.

Pippin, H. Gary↗

The granulite suite: Impact melts and metamorphic breccias of the early lunar crust

The granulite suite consists of two major types of rocks. One is coarse-grained and poikilitic with many euhedral crystals of olivine and plagioclase. These characteristics indicate crystallization from a melt; the poikilitic granulites are impact melt breccias. The other group is finer-grained and granoblastic, with numerous triple junctions; the granoblastic granulites are metamorphic rocks. Compositional groups identified by Lindstrom and Lindstrom contain both textural types. Two pyroxene thermometry indicates that both groups equilibrated at 1000 to 1150 C. Calculations suggest that the granoblastic group, which has an average grain size of about 80 microns, was annealed for less than 6 x 10 exp 4 y at 1000 C, and for less than 2500 y at 1150 C. Similar equilibration temperatures suggest that both groups were physically associated after impact events produced the poikilitic melts. Granulitic impactites hold important information about the pre-Nectarian bombardment history of the Moon, and the composition and thermal evolution of the early lunar crust. Granulitic impactites are widely considered to be an important rock type in the lunar crust, but how they formed is poorly understood. Metal compositions and elevated concentrations of meteoritic siderophile elements suggest that most lunar granulites are impact breccias. Their occurrence as clasts in approximately 3.9 Ga breccias, and Ar-(40-39) ages greater than or = 4.2 Ga for some granulites show that they represent a component of the lunar crust which formed prior to the Nectarian cataclysm. Petrographic characteristics of lunar granulites indicate at least two endmember textural variants which apparently formed in fundamentally different ways. One type has granoblastic textures consisting of equant, polygonal to rounded grains, and abundant triple junctions with small dispersions around 120 degrees indicating a close approach to textural equilibrium. As suggested by many authors, granoblastic granulites probably formed by subsolidus annealing and recrystallization of fragmental or glassy protoliths. Examples of this type include 15418, 78155, and 79215. The other textural type consists of poikilitic to poikiloblastic rocks with euhedral to subhedral plagioclase and olivine enclosed by interstitial pyroxene. In some cases, the texture resembles that of an orthocumulate. Examples of this type include 60035, 67955, and 77017. Rounding of grain edges is common in poikilitic granulites, but the regular crystal shapes and widely dispersed dihedral angles show they are far from textural equilibrium. The textures of poikilitic granulites are more consistent with the formation of these rocks by crystallization from a melt than by subsolidus metamorphism. A few samples have been recognized with textural characteristics transitional between those of the granoblastic and poikiloblastic endmembers (e.g., 72559, 78527). Pyroxene compositions taken from the literature and determined for this study by electron microprobe were used to calculate equilibration temperatures. The Kretz Ca transfer (solvus) thermometer and the Lindsley and Anderson graphical method both give similar temperatures, which range from approximately 1000 to 1150 C. There is no apparent temperature difference between granoblastic and poikilitic varieties, but there is a hint in these data that the more ferroan varieties equilibrated to lower temperatures. Additional studies are in progress to test this possibility.

Cushing, J. A.↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable 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) initiative. This dataset is 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. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current 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 simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

The Effect of Temperature on the Preservation of Volatile-Rich Lunar Samples

Introduction. The Moon’s south pole is a high-priority target for human exploration and scientific study. This interest is, in part, due to the presence of Permanently Shadowed Regions (PSRs), which could contain high concentrations of unique volatiles at cryogenic temperatures [1]. Returned samples from PSRs may include a unique combination of rocks, regolith, and volatile species, providing unprecedented insights into the history of the Solar System and the potential for resource utilization on the Moon. However, because PSR samples are cryogenic up-on collection, lunar polar sample return will eventually require cold stowage for the journey from the Moon to Earth. Without cold stowage, PSR sample return will likely result in phase changes and chemical reactions within the volatile component of the sample, which would negatively impact the resulting scientific studies of those samples. This abstract summarizes the initial results from an ongoing characterization of analog PSR samples at a range of temperatures, with the goal of defining the temperatures needed for a flight cold stowage freezer. Background. Based on remote sensing observations of the Moon [2], south polar PSRs range in temperature from ~120K for small and/or shallow PSRs to ~20K at the most extreme locations in large, deep PSRs. A range of volatiles have been hypothesized to exist at the surface or subsurface of the lunar poles [3-5 and others]. This hypothesis was verified when the LCROSS mission impacted the <50-K PSR in the crater Cabeus, detecting a range of volatiles from water to low condensation temperature species such as H2S and methane [6]. Species such as H2S and ammonia (also detected by LCROSS) are also highly reactive, and increase the likelihood of chemical reactions at elevated (non-cryogenic) temperatures. At the Johnson Space Center’s Planetary Exploration and Astromaterials Research Laboratory (JSC-PEARL), we have developed a volatile-bearing lunar simulant that incorporates several of the species detected by LCROSS [Table 1] mixed cryogenically with the USGS Lunar Highlands Type (LHT) regolith simulant. The new volatile-regolith simulant will be used to assess the degree of sample alteration at room temperature, -20°C, -80°C, and -196°C (liquid nitro-gen), over a two-week period. Room temperature samples represent those likely to be returned during initial missions without cold stowage, -20°C provides an analog to Apollo cold curated samples, -80°C is the temperature of multiple flight payload freezers (e.g., MELFI), and -196°C is analogous to lunar PSRs. Two weeks is an approximation of the time between sample collection and Earth return for initial Artemis missions. Over this period of time, sample head-space gases will be analyzed using a Universal Gas Analyzer (UGA, a type of mass spectrometer) coupled with a Baratron pressure sensor. After testing, the regolith component of the simulant will be purged of volatiles and preserved for future electron beam and/or FTIR analysis. Experimental Procedure. Volatile-regolith simulants will be produced as an initial homogenous batch; this batch will then be distributed into aliquots (gas chromatography/GC vials or cryo vials), ensuring that each sample has the same starting composition and conditions [Fig. 1]. In addition to the “full” simulant shown in Table 1, less complex simulant compositions will be used as baseline and control samples [Table 2]. Aliquots will be produced in triplicate for each simulant composition, storage temperature, and date of sampling. Headspace gases in all Day 0 samples will be analyzed by the UGA immediately. Cold storage samples for future analytical days will be placed in freezers appropriate to their target temperatures (-20°C, -80°C, -196°C). For ambient-temperature samples, regolith and regolith-water samples will be stored in a fume hood, while the full simulant will be stored in a sealed Parr vessel for safety; no other simulants (RWCM/ RWCM+) will be stored at ambient temperature for this test. Samples will be analyzed by UGA in this manner on each Analysis Day outlined in Table 2. Analytical Data. The UGA measures the partial pressures in a single sample aliquot over a set mass range of 0-105 atomic mass units (AMU) [Figure 2]; this set range was selected to slightly exceed the mass of the highest-mass expected reaction product (H2SO4). The Baratron complements the UGA by measuring the total pressure in the headspace of a sample vial. Coupled together, the quantitative abundances of gases will be monitored throughout the test. UGA analyses of the triplicate samples for each storage temperature, day, and simulant composition will be averaged, and standard deviations for each will be calculated. The compositions of starting species (shown in Table 1) will be characterized as a function of time, and the presence of any new compounds (reaction products) will be monitored as well. Total pressures will be recorded for each sample analysis, and any samples that show signs of leakage (e.g., a significant reduction in pressure or simulant volatiles) will be discarded. Anticipated Results. Testing is planned to begin in January 2022. The resulting data will allow compositional and phase changes in the volatile component of the simulants to be determined. Both the reduction in initial compounds and the addition of reaction products are expected to be observed. In addition, the relative efficacy of the different temperatures at pre-serving the initial composition of the simulants will be quantified. Finally, the regolith component of each sample will be argon-purged and stored in a controlled environment for future laboratory analysis. Compositional and morphological changes in the regolith are expected for samples above 0°C. This test will be repeated three times over the course of 2022. Understanding the effect of temperature on both the volatile and regolith components of analog lunar materials will allow requirements for a cold stowage freezer to be developed. The implementation of cold stowage for lunar polar missions will maximize the preservation of returned samples, enabling ground-breaking lunar and Solar System volatiles science for decades to come.

J L Mitchell↗

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗

Nano Icy Moons Propellant Harvester

As one of just a few bodies identified in the solar system with a liquid ocean, Europa has become a top priority in the search for life outside of Earth. However, cost estimates for exploring Europa have been prohibitively expensive, with estimates of a NASA Flagship class orbiter and lander approaching $5 billion. ExoTerra's NIMPH offers an affordable solution that can not only land, but return a sample from the surface to Earth. NIMPH combines solar electric propulsion (SEP) technologies being developed for the asteroid redirect mission and microsatellite electronics to reduce the cost of a full sample return mission below $500 million. A key to achieving this order-of-magnitude cost reduction is minimizing the initial mass of the system. The cost of any mission is directly proportional to its mass. By keeping the mission within the constraints of an Atlas V 551 launch vehicle versus an SLS, we can significantly reduce launch costs. To achieve this we reduce the landed mass of the sample return lander, which is the largest multiplier of mission mass, and shrink propellant mass through high efficiency SEP and gravity assists. The NIMPH projects first step in reducing landed mass focuses on development of a micro-In Situ Resource Utilization (micro-ISRU) system. ISRU allows us to minimize landed mass of a sample return mission by converting local ice into propellants. The project reduces the ISRU system to a CubeSat-scale package that weighs just 1.74 kg and consumes just 242 W of power. We estimate that use of this ISRU vs. an identical micro-lander without ISRU reduces fuel mass by 45 kg. As the dry mass of the lander grows for larger missions, these savings scale exponentially. Taking full advantage of the micro-ISRU system requires the development of a micro-liquid oxygen-liquid hydrogen engine. The micro-liquid oxygen-liquid hydrogen engine is tailored for the mission by scaling it to match the scale of the micro-lander and the low gravity of the target moon. We also tailor the engine for a near stoichiometric mixture ratio of 7.5. Most high-performance liquid oxygen-liquid hydrogen engines inject extra liquid hydrogen to lower the average molecular weight of the exhaust, which improves specific impulse. However, this extra liquid hydroden requires additional power and processing time on the surface for the ISRU to create. This increases mission cost, and on missions within high radiation environments such as Europa, increases radiation shielding mass. The resulting engine weighs just 1.36 kg and produces 71.5 newton of thrust at 364 s specific impulse. Finally, the mission reduces landed mass by taking advantage of the SEP modules solar power to beam energy to the surface using a collimated laser. This allows us to replace an 45 kg MMRTG with a 2.5 kg resonant array. By using the combination of ISRU, a liquid oxygen-liquid hydrogen engine, and beamed power, we reduce the initial mass of the lander to just 51.5 kg. When combined with an SEP module to ferry the lander to Europa the initial mission mass is just 6397 kg - low enough to be placed on an Earth escape trajectory using an Atlas V 551 launch vehicle. By comparison, we estimate a duplicate lander using an MMRTG and semi-storable propellants such as liquid oxygen-methane would result in an order of magnitude increase in initial lander mass to 445 kg. Attempting to perform the trajectory with a 450 s liquid oxygen-liquid hydrogen engine would increase initial mass to approximately 135,000 kg. Using an Atlas V 1 U.S. Dollar per kg rate to Earth escape value of $27.7k per kg, just the launch savings are over $3.5 billion.

liquid↗

The Nuclear Cryogenic Propulsion Stage

Nuclear Thermal Propulsion (NTP) development efforts in the United States have demonstrated the technical viability and performance potential of NTP systems. For example, Project Rover (1955 - 1973) completed 22 high power rocket reactor tests. Peak performances included operating at an average hydrogen exhaust temperature of 2550 K and a peak fuel power density of 5200 MW/m3 (Pewee test), operating at a thrust of 930 kN (Phoebus-2A test), and operating for 62.7 minutes in a single burn (NRX-A6 test). Results from Project Rover indicated that an NTP system with a high thrust-to-weight ratio and a specific impulse greater than 900 s would be feasible. Excellent results were also obtained by the former Soviet Union. Although historical programs had promising results, many factors would affect the development of a 21st century nuclear thermal rocket (NTR). Test facilities built in the US during Project Rover no longer exist. However, advances in analytical techniques, the ability to utilize or adapt existing facilities and infrastructure, and the ability to develop a limited number of new test facilities may enable affordable development, qualification, and utilization of a Nuclear Cryogenic Propulsion Stage (NCPS). Bead-loaded graphite fuel was utilized throughout the Rover/NERVA program, and coated graphite composite fuel (tested in the Nuclear Furnace) and cermet fuel both show potential for even higher performance than that demonstrated in the Rover/NERVA engine tests.. NASA's NCPS project was initiated in October, 2011, with the goal of assessing the affordability and viability of an NCPS. FY 2014 activities are focused on fabrication and test (non-nuclear) of both coated graphite composite fuel elements and cermet fuel elements. Additional activities include developing a pre-conceptual design of the NCPS stage and evaluating affordable strategies for NCPS development, qualification, and utilization. NCPS stage designs are focused on supporting human Mars missions. The NCPS is being designed to readily integrate with the Space Launch System (SLS). A wide range of strategies for enabling affordable NCPS development, qualification, and utilization should be considered. These include multiple test and demonstration strategies (both ground and in-space), multiple potential test sites, and multiple engine designs. Two potential NCPS fuels are currently under consideration - coated graphite composite fuel and tungsten cermet fuel. During 2014 a representative, partial length (approximately 16") coated graphite composite fuel element with prototypic depleted uranium loading is being fabricated at Oak Ridge National Laboratory (ORNL). In addition, a representative, partial length (approximately 16") cermet fuel element with prototypic depleted uranium loading is being fabricated at Marshall Space Flight Center (MSFC). During the development process small samples (approximately 3" length) will be tested in the Compact Fuel Element Environmental Tester (CFEET) at high temperature (approximately 2800 K) in a hydrogen environment to help ensure that basic fuel design and manufacturing process are adequate and have been performed correctly. Once designs and processes have been developed, longer fuel element segments will be fabricated and tested in the Nuclear Thermal Rocket Element Environmental Simulator (NTREE) at high temperature (approximately 2800 K) and in flowing hydrogen.

Houts, Michael G.↗

Retrieval of Daily Maps of PM2.5 Aerosol in the Problematic California Valleys: Bright, Speckled Reflectances, Thin AOT, but High Pollution

The San Joaquin Valley suffers from severe episodes of respirable aerosol (PM2.5) in wintertime.We provide maps of aerosol episodes using daily snapshots of PM2.5 and its changing features despite numerous difficulties inherent to sampling the region. Linear relationships relating aerosol optical thickness (AOT) to PM2.5 give an explained variance of approximately 3.The GEO-CAPE mission has as a goal the provision of relevant measures of respirable aerosol to the community,but has not formulated a science goal beyond the limited goal of retrieval of AOT, bringing the usefulness of GEO-CAPE into doubt.Our special focus was on the DISCOVER-AQ period, Jan-Feb 2013, which had many supporting measurements.Both high pollution and retrieval difficulties tend to occur in many Mediterranean agricultural regions like the San Joauin. One difficulty is the relatively bright surfaces with considerable exposed soil. NASAs MAIAC and MODIS Deep Blue retrieval techniques are shown to have considerable skill even at low aerosol optical thickness (AOT) values, as evaluated by concurrent AERONET sunphotometer measurements.More significantly, these AOT values can correspond to high daytime PM2.5 since aerosol mixed layer depth is thin and variable, 200m 600 m. The thin layers derive from typical subsidence of dry air between more stormy periods. This situation provides an advantage: water vapor column is also almost completely limited to a similar mixed layer depth, and can thus serve as a measure of aerosol dilution.Using the MAIAC Water Vapor Column:In order to make the maps below, we used the MAIAC data but subtracted partial water-vapor columns estimated from MERRA Reanalysis Data availabe from the GSFC GMAO using kriging. We did not use the mixed-layer estimates from MERRA, since such analyses were found problematic during our forecasting exercises for DISCOVER-AQ. Observations from the aircraft soundings suggested that this overlying moisture was mostly due to larger scale flows, not ML venting.However, the specific humidity at the surface and a nearly well-mixed ML was analyzed by kriging from the surface network (MesoWest, University of Utah). These were thought to be truer, uninfluenced by physical process modeling that combines with data observations. (TBD: How different are they?) Procedure: Subtract overlying partial water columns from MAIAC column water and divide this by a surface value of water vapor. (MAIAC column is expressed in cm of water, i.e., water vapor at surface conditions.This method appears to bring out useful details in the distribution of submicron particles in the very problematic Wintertime San Joaquin Valley, and allow analysis of pollution episodes throughout the valley, rather than long-term averages.

san joaquin valley↗

Comparative Soot Diagnostics: Preliminary Results

The motivation for the Comparative Soot Diagnostics (CSD) experiment lies in the broad practical importance of understanding combustion generated particulate. Depending upon the circumstances, particulate matter can affect the durability and performance of combustion equipment, can be a pollutant, can be used to detect fires and, in the form of soot, can be the dominant source of radiant energy from flames. The nonbuoyant structure of most flames of practical interest makes understanding of soot processes in low gravity flames important to our ability to predict fire behavior on earth. These studies also have direct applications to fire safety in human-crew spacecraft, since smoke is the indicator used for automated detection in current spacecraft. In the earliest missions (Mercury, Gemini and Apollo), the crew quarters were so cramped that it was considered reasonable that the astronauts would rapidly detect any fire. The Skylab module, however, included approximately 20 UV-sensing fire detectors. The Space Shuttle has 9 particle-ionization smoke detectors in the mid-deck and flight deck and Spacelab has six additional particle-ionization smoke detectors. The designated detectors for the ISS are laser-diode, forward-scattering, smoke or particulate detectors. Current plans for the ISS call for two detectors in the open area of the module, and detectors in racks that have both cooling air flow and electrical power. Due to the complete absence of data concerning the nature of particulate and radiant emission from incipient and fully developed low-g fires, all three of these detector systems were designed based upon l-g test data and experience. As planned mission durations and complexity increase and the volume of spacecraft increases, the need for and importance of effective, crew-independent, fire detection grows significantly. To provide this level of protection, more knowledge is needed concerning low-gravity fire phenomena and, in particular, how they might be detected and suppressed. Prior to CSD, no combustion-generated particulate samples had been collected near the flame zone for well-developed microgravity flames. All of the extant data either came from drop tower tests and therefore only corresponded to the early stages of a fire or were collected far from the flame zone. The fuel sources in the drop tower tests were restricted to laminar gas-jet diffusion flames and very rapidly overheated wire insulation. The gas-jet tests indicated, through thermophoretic sampling, (2) that soot primaries and aggregates (groups of primary particles) in low-gravity may be significantly larger than those in normal gravity (1-g). This raises new scientific questions about soot processes as well as practical issues for particulate size sensitivity and detection alarm threshold levels used in on-orbit smoke detectors. Preliminary tests in the 2.2 second drop tower suggest that particulate generated by overheated wire insulation may be larger in low-g than in 1-g. Transmission Electron Microscope (TEM) grids downstream of the fire region in the Wire Insulation Flammability experiment as well as visual observation of long string-like aggregates, further confirm this suggestion. The combined impact of these limited results and theoretical predictions is that, as opposed to extrapolation from l-g data, direct knowledge of low-g combustion particulate is needed for more confident design of smoke detectors for spacecraft. This paper describes the operation and preliminary results of the CSD, a project conceived and developed at NASA Lewis Research Center. The CSD flight experiment was conducted in the Middeck Glovebox Facility (MGBX) on USMP-3. The project is support by NASA Headquarters Microgravity Science and Applications Division and Code Q. The results presented here are from the microgravity portion of the experiment, including the temporal response of the detectors and average sizes of the primary and aggregate particles captured on the thermophoretic probes.

Urban, David L.↗

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗

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↗