Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normalization”

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,549 records · Page 86

Metal Side Reflectors for Trapping Light in QWIPs

Focal-plane arrays of quantum-well infrared photodetectors (QWIPs) equipped with both light-coupling diffraction gratings and metal side reflectors have been proposed, and prototypes are expected to be fabricated soon. The purpose served by the metal side reflectors is to increase quantum efficiency by helping to trap light in the photosensitive material of each pixel. The reasons for using diffraction gratings were discussed in several prior NASA Tech Briefs articles. To recapitulate: In an array of QWIPs, the quantum-well layers are typically oriented parallel to the focal plane and therefore perpendicular or nearly perpendicular to the direction of incidence of infrared light. By virtue of the applicable quantum selection rules, light polarized parallel to the focal plane (as normally incident light is) cannot excite charge carriers and, hence, cannot be detected. Diffraction gratings scatter normally or nearly normally incident light into directions more nearly parallel to the focal plane, so that a significant portion of the light attains a component of polarization normal to the focal plane and, hence, can excite charge carriers. Unfortunately, light scattered in directions parallel or nearly parallel to the focal plane can escape sideways from the QWIP of a given pixel, as illustrated in Figure 1. The escaped light has made only a single pass through the interior photosensitive volume of the QWIP. The quantum efficiency of the QWIP would be increased by trapping light so that it makes multiple passes through the photosensitive volume. As shown in Figure 2, the sides of the QWIP of each pixel would be coated with gold to reflect escaping light back into the interior.

Gunapala, Sarath↗

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Evaluation of Delamination Growth Characterization Methods Under Mode I Fatigue Loading

Reliable delamination characterization data for laminated composites are needed for input to analytical models of structures to predict delamination. The double-cantilevered beam (DCB) specimen is used with laminated composites to measure fracture toughness, G(sub Ic), delamination onset strain energy release rate, and growth rate data under cyclic loading. In the current study, DCB specimens of IM7/8552 graphite/epoxy supplied by two different manufacturers were tested in static and fatigue to compare the measured characterization data from the two sources, and to evaluate a proposed ASTM standard for generating Paris Law equations. Static results were used to generate compliance calibration constants for the fatigue data, and a delamination resistance curve, G(sub IR), which was used to determine the effects of fiber-bridging on delamination growth. Specimens were tested in fatigue at a cyclic G(sub Imax) level equal to 50, 40 or 30% of G(sub Ic), to determine a delamination onset curve and delamination growth rate. The delamination onset curve equations had similar exponents and the same trends. Delamination growth rate was calculated by fitting a Paris Law to the da/dN versus G(sub Imax) data. Both a 2-point and a 7-point data reduction method were used and the Paris Law equations were compared. To determine the effects of fiber-bridging, growth rate results were normalized by the delamination resistance curve for each material and compared to the non-normalized results. Paris Law exponents were found to decrease by 31% to 37% due to normalizing the growth data. Normalizing the data also greatly reduced the amount of scatter between the different specimens. Visual data records from the fatigue testing were used to calculate individual compliance calibration constants from the fatigue data for some of the specimens. The resulting da/dN versus G(sub Imax) plots showed much improved repeatability between specimens. Gretchen

Murri, Gretchen B.↗

Quasi One-Dimensional Unsteady Modeling of External Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of an axisymmetric external compression supersonic inlet is being developed. The model utilizes compressible flow computational fluid dynamics to model the internal inlet segment as well as the external inlet portion between the cowl lip and normal shock, and compressible flow relations with flow propagation delay to model the oblique shocks upstream of the normal shock. The external compression portion between the cowl-lip and the normal shock is also modeled with leaking fluxes crossing the sonic boundary, with a moving CFD domain at the normal shock boundary. This model has been verified in steady state against tunnel inlet test data and it s a first attempt towards developing a more comprehensive model for inlet dynamics.

Kopasakis, George↗

Effects of Radiation and Dietary Iron on Expression of Genes and Proteins Involved in Drug Metabolism

Liver function, especially the rate of metabolic enzyme activities, determines the concentration of circulating drugs and the duration of their efficacy. Most pharmaceuticals are metabolized by the liver, and clinically-used medication doses are given with normal liver function in mind. A drug overdose can result in the case of a liver that is damaged and removing pharmaceuticals from the circulation at a rate slower than normal. Alternatively, if liver function is elevated and removing drugs from the system more quickly than usual, it would be as if too little drug had been given for effective treatment. Because of the importance of the liver in drug metabolism, we want to understand any effects of spaceflight on the enzymes of the liver. Dietary factors and exposure to radiation are aspects of spaceflight that are potential oxidative stressors and both can be modeled in ground experiments. In this experiment, we examined the effects of high dietary iron and low dose gamma radiation (individually and combined) on the gene expression of enzymes involved in drug metabolism, redox homeostasis, and DNA repair. METHODS All procedures were approved by the JSC Animal Care and Use Committee. Male Sprague-Dawley rats were divided into 4 groups (n=8); control, high Fe diet (650 mg iron/kg), radiation (fractionated 3 Gy exposure from a Cs- 137 source) and combined high Fe diet + radiation exposure. Animals were euthanized 24h after the last treatment of radiation; livers were removed immediately and flash -frozen in liquid nitrogen. Expression of genes thought to be involved in redox homeostasis, drug metabolism and DNA damage repair was measured by RT-qPCR. Where possible, protein expression of the same genes was measured by western blotting. All data are expressed as % change in expression normalized to reference gene expression; comparisons were then made of each treatment group to the sham exposed/ normal diet control group. Data was considered significant at p< 0.5. RESULTS Among the redox homeostasis genes examined, metallothionein showed a significant down regulation in the radiation treated group (-3.85 fold) and a trend toward down regulation in the high Fe + rad group. Metallothionein is involved in the regulation of physiological metals and also has antioxidant activities. Among the drug metabolism genes examined, ATP binding cassette subfamily B (Abcb1b) gene expression increased more than 10-fold in both groups that received radiation treatments. This increased expression was also seen at the protein level. This ABC transporter carries many different compounds across cell membranes, including administered medications. The cytochrome P450 2E1 enzyme, a mixed-function oxidase that deactivates some medications and activates others, showed about a 2-fold increase in gene expression in both radiation-treated groups, with a trend toward increased expression at the protein level. Expression of epoxide hydrolase, which detoxifies polycyclic aromatic hydrocarbons, showed similar 2-fold increases. Among the DNA repair genes examined, expression of RAD51 was significantly down regulated (1.5 fold) in the radiation treated group. RAD51 is involved in repair of double-stranded DNA breaks. CONCLUSION This experiment used 2 different sources of physiological oxidative stress, administered separately and together, and examined their impacts on liver gene and protein expression. It is clear that significant changes occurred in expression of several genes and proteins in the radiation-treated animals. If the results from this ground analog of portions of the spaceflight environment hold true for the spaceflight environment itself, the physiological roles of the affected enzymes (drug transport and metabolism, redox homeostasis) could mean consequences in redox homeostasis or the pharmacokinetics of administered medications

Faust, K. M.↗

Vertical Structure of Ice Cloud Layers From CloudSat and CALIPSO Measurements and Comparison to NICAM Simulations

The shape of the vertical profile of ice cloud layers is examined using 4 months of CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) global measurements taken on January, April, July, and October 2007. Ice clouds are selected using temperature profiles when the cloud base is located above the 253K temperature level. The obtained ice water content (IWC), effective radius, or extinction coefficient profiles are normalized by their layer mean values and are expressed in the normalized vertical coordinate, which is defined as 0 and 1 at the cloud base and top heights, respectively. Both CloudSat and CALIPSO observations show that the maximum in the IWC and extinction profiles shifts toward the cloud bottom, as the cloud depth increases. In addition, clouds with a base reaching the surface in a high-latitude region show that the maximum peak of the IWC and extinction profiles occurs near the surface, which is presumably due to snow precipitation. CloudSat measurements show that the seasonal difference in normalized cloud vertical profiles is not significant, whereas the normalized cloud vertical profile significantly varies depending on the cloud type and the presence of precipitation. It is further examined if the 7 day Nonhydrostatic Icosahedral Atmospheric Model (NICAM) simulation results from 25 December 2006 to 1 January 2007 generate similar cloud profile shapes. NICAM IWC profiles also show maximum peaks near the cloud bottom for thick cloud layers and maximum peaks at the cloud bottom for low-level clouds near the surface. It is inferred that oversized snow particles in the NICAM cloud scheme produce a more vertically inhomogeneous IWC profile than observations due to quick sedimentation.

Ham, Seung-Hee↗

Influence of Mixed Mode I-Mode II Loading on Fatigue Delamination Growth Characteristics of a Graphite Epoxy Tape Laminate

Mixed mode I-mode II interlaminar tests were conducted on IM7/8552 tape laminates using the mixed-mode bending test. Three mixed mode ratios, G(sub II)/G(sub T) = 0.2, 0.5, and 0.8, were considered. Tests were performed at all three mixed-mode ratios under quasi-static and cyclic loading conditions, where the former static tests were used to determine initial loading levels for the latter fatigue tests. Fatigue tests at each mixed-mode ratio were performed at four loading levels, Gmax, equal to 0.5G(sub c), 0.4G(sub c), 0.3G(sub c), and 0.2G(sub c), where G(sub c) is the interlaminar fracture toughness of the corresponding mixed-mode ratio at which a test was performed. All fatigue tests were performed using constant-amplitude load control and delamination growth was automatically documented using compliance solutions obtained from the corresponding quasi-static tests. Static fracture toughness data yielded a mixed-mode delamination criterion that exhibited monotonic increase in Gc with mixed-mode ratio, G(sub II)/G(sub T). Fatigue delamination onset parameters varied monotonically with G(sub II)/G(sub T), which was expected based on the fracture toughness data. Analysis of non-normalized data yielded a monotonic change in Paris law exponent with mode ratio. This was not the case when normalized data were analyzed. Fatigue data normalized by the static R-curve were most affected in specimens tested at G(sub II)/G(sub T)=0.2 (this process has little influence on the other data). In this case, the normalized data yielded a higher delamination growth rate compared to the raw data for a given loading level. Overall, fiber bridging appeared to be the dominant mechanism, affecting delamination growth rates in specimens tested at different load levels and differing mixed-mode ratios.

Ratcliffe, James G.↗

The Diversity of Anomalous HEDs: Isotopic Constraints on the Connection of EET 92023, GRA 98098, and Dhofar 700 With Vesta

The possibility for multiple parent bodies, instead of a common parent body of Vesta, for eucrites has been suggested based on the variable oxygen isotopic composition observed in some eucrites.. Recently, we added an extra dimension to the discussion based on the (epsilon)54Cr composition of the same eucrites with known (delta)17O to compare with the normal eucrites. The combined (delta)17O and (epsilon)54Cr isotope systematics for Pasamonte, PCA 91007, A-881394, and Ibitira indicate their likely origin from multiple different parent bodies than the normal eucrites. Often the qualifier anomalous is used to identify HEDs with (delta)17O values that deviate significantly (>3(sigma)) from the mean HED (delta)17O. However, variations in eucrites and diogenites also include unique geochemical characteristics such as bulk composition, trace element abundances, or volatile concentrations, in addition to (delta)17O. Here, we investigate three such geochemically anomalous HEDs: Elephant Moraine (EET) 92023, Graves Nunataks (GRA) 98098, and Dhofar 700. In addition, to verify the homogeneity of (epsilon)54Cr observed for normal HEDs thus far, a set of seven eucrites and diogenites considered normal samples were also investigated.

Sanborn, M. E.↗

Mixed Convection Blowoff Limits as a Function of Oxygen Concentration and Upward Forced Stretch Rate for Burning Pmma Rods of Various Sizes

Normal gravity flame blowoff limits in an axisymmetric pmma rod geometry in upward axial stagnation flow are compared with microgravity Burning and Suppression of Solids II (BASS-II) results recently obtained aboard the International Space Station. This testing utilized the same BASS-II concurrent rod geometry, but with the addition of normal gravity buoyant flow. Cast polymethylmethacrylate (pmma) rods of diameters ranging from 0.635 cm to 3.81 cm were burned at oxygen concentrations ranging from 14 to 18 by volume. The forced flow velocity where blowoff occurred was determined for each rod size and oxygen concentration. These blowoff limits compare favorably with the BASS-II results when the buoyant stretch is included and the flow is corrected by considering the blockage factor of the fuel. From these results, the normal gravity blowoff boundary for this axisymmetric rod geometry is determined to be linear, with oxygen concentration directly proportional to flow speed. We describe a new normal gravity upward flame spread test method which extrapolates the linear blowoff boundary to the zero stretch limit to resolve microgravity flammability limits, something current methods cannot do. This new test method can improve spacecraft fire safety for future exploration missions by providing a tractable way to obtain good estimates of material flammability in low gravity.

blowoff↗

Quantifying Pilot Contribution to Flight Safety During Dual Generator Failure

Accident statistics cite flight crew error in over 60% of accidents involving transport category aircraft. Yet, a well-trained and well-qualified pilot is acknowledged as the critical center point of aircraft systems safety and an integral safety component of the entire commercial aviation system. No data currently exists that quantifies the contribution of the flight crew in this role. Neither does data exist for how often the flight crew handles non-normal procedures or system failures on a daily basis in the National Airspace System. A pilot-in-the-loop high fidelity motion simulation study was conducted by the NASA Langley Research Center in partnership with the Federal Aviation Administration (FAA) to evaluate the pilot's contribution to flight safety during normal flight and in response to aircraft system failures. Eighteen crews flew various normal and non-normal procedures over a two-day period and their actions were recorded in response to failures. To quantify the human's contribution, crew complement was used as the experiment independent variable in a between-subjects design. Pilot actions and performance when one of the flight crew was unavailable were also recorded for comparison against the nominal two-crew operations. This paper details diversion decisions, perceived safety of flight, workload, time to complete pertinent checklists, and approach and landing results while dealing with a complete loss of electrical generators. Loss of electrical power requires pilots to complete the flight without automation support of autopilots, flight directors, or auto throttles. For reduced crew complements, the additional workload and perceived safety of flight was considered unacceptable.

Etherington, Timothy J.↗

Unusual Black Hole Binary LMC X-3: A Transient High-Mass X-Ray Binary That Is Almost Always On?

We have analyzed a rich, multimission, multiwavelength data set from the black hole X-ray binary (BHXB) LMC X-3, covering a new anomalous low state (ALS), during which the source flux falls to an unprecedentedly low and barely detectable level, and a more normal low state. Simultaneous X-ray and UV/optical monitoring data from Swift are combined with pointed observations from the Rossi X-ray Timing Explorer (RXTE) and X-ray Multi- Mirror Mission (XMM-Newton) and light curves from the Monitor of All-Sky X-ray Image (MAXI) instrument to compare the source characteristics during the ALS with those seen during the normal low state. An XMM-Newton spectrum obtained during the ALS can be modeled using an absorbed power law with Gamma = 1.41‚+/- 0.65 and a luminosity of 7.97 x 10(exp 33) erg/s (0.6-5 keV). The Swift X-ray and UV light curves indicate an X-ray lag of approx. 8 days as LMC X-3 abruptly exits the ALS, suggesting that changes in the mass accretion rate from the donor drive the X-ray lag. The normal low state displays an asymmetric profile in which the exit occurs more quickly than the entry, with minimum X-ray flux a factor of approx. 4300 brighter than during the ALS. The UV brightness of LMC X-3 in the ALS is also fainter and less variable than during normal low states. The existence of repeated ALSs in LMC X-3, as well as a comparison with other BHXBs, implies that it is very close to the transient/persistent X-ray source dividing line. We conclude that LMC X-3 is a transient source that is almost always "on."

Simultaneous X-ray and UV/optical monitoring dat↗

Load Schedule Design Improvements for the Calibration of NASA's Semi-span Balances

Improvements of the load schedule design for the calibration of NASA's semi-span balances are discussed. These improvements became possible after a minor hardware change was made to Calspan's Large Load Rig. This change moves normal and axial force load points 12 inches closer to the balance moment center. Resulting data quality improvements are significant because Calspan's Large Load Rig cannot currently apply forces at the balance moment center of a semi-span balance. Now, the supported set of rolling and yawing moment arms is 28 inches, 38 inches, and 48 inches. Consequently, the new load schedule design better supports expected normal and axial force locations on semi-span models that are tested in the NASA Ames 11-ft Transonic Wind Tunnel. In addition, estimates of the primary sensitivities of the normal force, axial force, rolling moment, and yawing moment have become more accurate as the applied normal and axial forces have moved closer to the balance moment center. Fundamental moment arm range requirements and Calspan's hardware updates are discussed in detail. In addition, calibration analysis results of NASA's ARC30K semi-span balance are presented to illustrate benefits of the hardware changes.

wind tunnel test↗

Climate Conditions During a Rift Valley Fever Post-epizootic Period in Free State, South Africa, 2014–2019

Rift Valley fever virus (RVFV) activity in Southern Africa tends to occur during periods of sustained elevated rainfall, cooler than normal conditions, and abundant vegetation cover creating ideal conditions for the increase and propagation of populations of RVFV mosquito vectors. These climatic and ecological conditions are modulated by large-scale tropical-wide El Niño–Southern Oscillation (ENSO) phenomena. The aim of this 5-year study was to investigate climatic conditions during Rift Valley fever “post-epizootic” period in Free State province of the Republic of South Africa, which historically experienced the largest RVF outbreaks in this country. We collected satellite-derived rainfall, land surface temperature (LST), and normalized difference vegetation index (NDVI) data since 2014 to understand broad environmental conditions in the years following a period of sustained and widespread large RVF outbreaks (2008–2011) in the region. We found this post-epizootic/interepizootic period to be characterized by below-normal rainfall (~-500 mm), above LSTs (~+12°C), depressed NDVI (60% below normal), and severe drought as manifested particularly during the 2015–2016 growing season. Such conditions reduce the patchwork of appropriate habitats available for emergence of RVFV vectors and diminish chances of RVFV activity. However, the 2016–2017 growing season saw a marked return to somewhat wetter conditions without any reported RVFV transmission. In general, the aggregate vector collections during this 5-year period follow patterns observed in climate measurements. During the 2017–2018 growing season, late and seasonally above average rainfall resulted in a focal RVF outbreak in one location in the study region. This unanticipated event is an indicator of cryptic RVF activity during post-epizootic period and may be a harbinger of RVFV activity in the coming years.

climate variability↗

Load Schedule Design Improvements for the Calibration of NASA's Semi-span Balances

Improvements of the load schedule design for the calibration of NASA's semi-span balances are discussed. These improvements became possible after a minor hardware change was made to Calspan's Large Load Rig. This change moves normal and axial force load points 12 inches closer to the balance moment center. Resulting data quality improvements are significant because Calspan's Large Load Rig cannot currently apply forces at the balance moment center of a semi-span balance. Now, the supported set of rolling and yawing moment arms is 28, 38, and 48 inches. Consequently, the new load schedule design better supports expected normal and axial force locations on semi-span models that are tested in the NASA Ames 11-ft Transonic Wind Tunnel. In addition, estimates of the primary sensitivities of the normal force, axial force, rolling moment, and yawing moment have become more accurate as the applied normal and axial forces have moved closer to the balance moment center. Fundamental moment arm range requirements and Calspan's hardware updates are discussed. In addition, calibration and check load analysis results of NASA's ARC30K semi-span balance are presented to illustrate benefits of the hardware changes.

wind tunnel test↗

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↗

Formaldehyde evolution in U.S. wildfire plumes during FIREX-AQ

Formaldehyde (HCHO) is one of the most abundant non-methane volatile organic compounds (VOCs) emitted by fires. HCHO also undergoes chemical production and loss as a fire plume ages, and it can be an important oxidant precursor. In this study, we disentangle the processes controlling HCHO by examining its evolution in wildfire plumes sampled by the NASA DC-8 during the FIREX-AQ field campaign. In nine of the twelve analyzed plumes, dilution-normalized HCHO increases with physical age (range 1 – 6 h). The balance of HCHO loss (mainly via photolysis) and production (via OH-initiated VOC oxidation) seems to control the sign and magnitude of this trend. Plume-average OH concentrations, calculated from VOC decays, range from –0.5 (± 0.5) × 10 6 to 5.3 (±0.7) × 10 6 cm -3 . The production and loss rates of dilution-normalized HCHO seem to decrease with plume age. Plume-to-plume variability in dilution-normalized secondary HCHO production correlates with OH abundance rather than normalized OH reactivity, suggesting that OH is the main driver of fire-to-fire variability in HCHO secondary production. Analysis suggests an effective HCHO yield of 0.33 (± 0.05) per VOC molecule oxidized for the 12 wildfire plumes. This finding can help connect space-based HCHO observations to the oxidizing capacity of the atmosphere and to VOC emissions.

Air quality↗

Optimizing Landsat Next Shortwave Infrared Bands for Crop Residue Characterization

This study focused on optimizing the placement of shortwave infrared (SWIR) bands for pixel-level estimation of fractional crop residue cover (f R ) for the upcoming Landsat Next mission. We applied an iterative wavelength shift approach to a database of crop residue field spectra collected in Beltsville, Maryland, USA (n = 916) and computed generalized two- and three-band spectral indices for all wavelength combinations between 2000 and 2350 nm, then used these indices to model field-measured f R . A subset of the full dataset with a Normalized Difference Vegetation Index (NDVI) < 0.3 threshold (n = 643) was generated to evaluate green vegetation impacts on f R estimation. For the two-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized normalized difference using 2226 nm and 2263 nm bands produced the top f R estimation performance (R 2 = 0.8222; RMSE = 0.1296). These findings were similar to the established two-band Shortwave Infrared Normalized Difference Residue Index (SINDRI) (R 2 = 0.8145; RMSE = 0.1324). Performance of the two-band generalized normalized difference and SINDRI decreased for the full-NDVI dataset (R 2 = 0.5865 and 0.4144, respectively). For the three-band wavelength shift analyses applied to the NDVI < 0.3 dataset, a generalized ratio-based index with a 2031–2085–2216 nm band combination, closely matching established Cellulose Absorption Index (CAI) bands, was top performing (R 2 = 0.8397; RMSE = 0.1231). Three-band indices with CAI-type wavelengths maintained top f R estimation performance for the full-NDVI dataset with a 2036–2111–2217 nm band combination (R 2 = 0.7581; RMSE = 0.1548). The 2036–2111–2217 nm band combination was also top performing in f R estimation (R 2 = 0.8690; RMSE = 0.0970) for an additional analysis assessing combined green vegetation cover and surface moisture effects. Our results indicate that a three-band configuration with band centers and wavelength tolerances of 2036 nm (±5 nm), 2097 nm (±14 nm), and 2214 (±11 nm) would optimize Landsat Next SWIR bands for f R estimation.

Landsat Next↗

X-ray analysis of fatigue damage in copper

A Warren-Averbach x-ray line profile analysis was applied to, broadened x-ray diffraction peaks from copper deformed in fatigue. The copper specimens were fatigued by four point bending at various strain amplitudes, and measurements were made at various fractions of the total fatigue life. The analysis results in an estimation of 1) an average coherently diffracting domain size in a direction normal to the diffracting planes; and 2) a strain distribution function. The strain distribution function is in terms of the root mean square strain, <ϵ^(2)> where ϵ is averaged over a given distance at all points in the diffracting crystals. This rms strain is expressed as a function of averaging distance for a direction normal to the diffracting planes. The rms strains were found to be a function of the fatigue strain amplitude level, increasing with amplitude and reaching a maximum limiting value during the first 2% of the fatigue life. The strains were normally larger in the 4lO0 direction compared to the <111> direction and the absolute values were quite small. These results may be explained by assuming the dislocations are present in the form of dipoles because the lattice strain around a dipole is much less than around single dislocations. Conclusions arrived at in this analysis are in agreement with current theories of fatigue damage in metals. The result that maximum hardening is achieved during the initial few percent of the fatigue life is confirmed by many other investigators using entirely different analytical techniques. It is felt that the structure of the material is changing during this initial hardening stage to an eventual steady state characteristic of the testing conditions. A unique substructure is developed which thereafter persists, and the remainder of the fatigue test is taken up by crack propagation through this substructure until ultimate failure occurs by fracture. The formation of this steady state structure is dependent on the maximum fatigue strain amplitude and on the frequency and temperature of cycling for a particular material.

COPPER↗