Search NASA⌕ Search

SEARCH · Search NASA

Results for “Logging”

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 613 records · Page 34

Model-Invariant Hybrid RANS-LES Computations on Unstructured Meshes

Hybrid RANS-LES computations combine the bene ts of RANS and LES so that LES is used in regions where the accuracy of RANS deteriorates. The numerous hybrid approaches are limited by the speci cation of the LES-RANS interface, which can cause nonphysical results such as log-layer mismatch and low shear stress. The hybrid RANS-LES approach based on the concept of model invariance, mitigates these problems, enabling seamless blending of the RANS and LES regions while forming the basis for interpreting the results in the interface region. This hybrid formulation was implemented in the NASA FUN3D unstructured code and computations for ow in a channel at Reynolds number of 3300 (based on the channel half width h and the bulk in ow velocity u(infinity) were carried out. An isotropic stochastic turbulence generator was implemented to generate in ow turbulence. The present approach was able to eliminate the log-layer mismatch and predict the shear stress fairly well. Thus, the model-invariant hybrid formulation coupled with the isotropic turbulence in ow generation serves as a physically meaningful way of performing hybrid RANS-LES computations.

Sreekanth, R.↗

Safe2Ditch Steer-To-Clear Development and Flight Testing

This paper describes a series of small unmanned aerial system (sUAS) flights performed at NASA Langley Research Center in April and May of 2019 to test a newly added Steer-to-Clear feature for the Safe2Ditch (S2D) prototype system. S2D is an autonomous crash management system for sUAS. Its function is to detect the onset of an emergency for an autonomous vehicle, and to enable that vehicle in distress to execute safe landings to avoid injuring people on the ground or damaging property. Flight tests were conducted at the City Environment Range for Testing Autonomous Integrated Navigation (CERTAIN) range at NASA Langley. Prior testing of S2D focused on rerouting to an alternate ditch site when an occupant was detected in the primary ditch site. For Steer-to-Clear testing, S2D was limited to a single ditch site option to force engagement of the Steer-to-Clear mode. The implementation of Steer-to-Clear for the flight prototype used a simple method to divide the target ditch site into four quadrants. An RC car was driven in circles in one quadrant to simulate an occupant in that ditch site. A simple implementation of Steer-to- Clear was programmed to land in the opposite quadrant to maximize distance to the occupant’s quadrant. A successful mission was tallied when this occurred. Out of nineteen flights, thirteen resulted in successful missions. Data logs from the flight vehicle and the RC car indicated that unsuccessful missions were due to geolocation error between the actual location of the RC car and the derived location of it by the Vision Assisted Landing component of S2D on the flight vehicle. Video data indicated that while the Vision Assisted Landing component reliably identified the location of the ditch site occupant in the image frame, the conversion of the occupant’s location to earth coordinates was sometimes adversely impacted by errors in sensor data needed to perform the transformation. Logged sensor data was analyzed to attempt to identify the primary error sources and their impact on the geolocation accuracy. Three trends were observed in the data evaluation phase. In one trend, errors in geolocation were relatively large at the flight vehicle’s cruise altitude, but reduced as the vehicle descended. This was the expected behavior and was attributed to sensor errors of the inertial measurement unit (IMU). The second trend showed distinct sinusoidal error for the entire descent that did not always reduce with altitude. The third trend showed high scatter in the data, which did not correlate well with altitude. Possible sources of observed error and compensation techniques are discussed.

Petty, Bryan J.↗

Atmospheric Characterization and Further Orbital Modeling of κ Andromeda b

We present κ Andromeda b's photometry and astrometry taken with Subaru/SCExAO+HiCIAO and Keck/NIRC2, combined with recently published SCExAO/CHARIS low-resolution spectroscopy and published thermal infrared photometry to further constrain the companion's atmospheric properties and orbit. The Y/Y−K colors of κ And b are redder than field dwarfs, consistent with its youth and lower gravity. Empirical comparisons of its Y-band photometry and CHARIS spectrum to a large spectral library of isolated field dwarfs reaffirm the conclusion from Currie et al. that it likely has a low gravity but admit a wider range of most plausible spectral types (L0–L2). Our gravitational classification also suggests that the best-fit objects for κ And b may have lower gravity than those previously reported. Atmospheric models lacking dust/clouds fail to reproduce its entire 1–4.7 μm spectral energy distribution (SED), and cloudy atmosphere models with temperatures of ~1700–2000 K better match κ And b data. Most well-fitting model comparisons favor 1700–1900 K, a surface gravity of log(g) ~ 4–4.5, and a radius of 1.3–1.6 R (sub Jup); the best-fit model (Drift-Phoenix) yields the coolest and lowest-gravity values: T (sub eff) = 1700 K and log g = 4.0. An update to κ And b's orbit with ExoSOFT using new astrometry spanning 7 yr reaffirms its high eccentricity (0.77 ± 0.08). We consider a scenario where unseen companions are responsible for scattering κ And b to a wide separation and high eccentricity. If three planets, including κ And b, were born with coplanar orbits, and one of them was ejected by gravitational scattering, a potential inner companion with mass ≲10 M (sub Jup) could be located at ≲25 au.

Exoplanet astronomy↗

The Warm Gas in the MW: A Kinematical Model

We develop a kinematical model for the Milky Way Si IV-bearing gas to determine its density distribution and kinematics. This model is constrained by a column density line-shape sample extracted from the Hubble Space Telescope/Cosmic Origins Spectrograph archival data, which contains 186 active galactic nucleus sight lines. We find that the Si IV ion density distribution is dominated by an extended disk along the z-direction (above or below the midplane), i.e., n(z) = n0 exp(-(z z/0 ) 0.82 , where z0 is the scale height of - 6.3+1.5 -1.6 kpc (northern hemisphere) and 3.6+0.9 -1.0 kpc (southern hemisphere). The density distribution of the disk in the radial direction shows a sharp edge at 15–20 kpc given by, n(rXY) = n0 exp(-(rXY/r0) 3.36 , where r0≈12.5±0.6 kpc. The difference of density distributions over rXY and z directions indicates that the warm gas traced by Si IV is mainly associated with disk processes (e.g., feedback or cycling gas) rather than accretion. We estimate the mass of the warm gas (within 50 kpc) is log(M(50kpc) M(sun)) ≈ 8.1 (assuming Z≈0.5 Z(sun)), and a 3σ upper limit of log(M(250kpc) M(sun)) ≈ 9.1 (excluding the Magellanic system). Kinematically, the warm gas disk is nearly co-rotating with the stellar disk at vrot = 215 ± 3 km s-1, which lags the midplane rotation by about 8 km s-1 kpc-1 (within 5 kpc). Meanwhile, we note that the warm gas in the northern hemisphere has significant accretion with vacc of 69±7 km s-1 at 10 kpc (an accretion rate of -0.60+0.13- 0.11 M(sun) yr-1), while in the southern hemisphere, there is no measurable accretion, with an upper limit of 0.4M(sun) yr−1.

Zhijie Qu↗

Resolving the Soft X-Ray Ultrafast Outflow in PDS 456

Past X-ray observations of the nearby luminous quasar PDS 456 (at z = 0.184) have revealed a wide-angle accretion disk wind with an outflow velocity of ∼−0.25c, as observed through observations of its blueshifted iron K-shell absorption line profile. Here we present three new XMM-Newton observations of PDS 456: one in 2018 September where the quasar was bright and featureless and two in 2019 September, 22 days apart, occurring when the quasar was five times fainter and where strong blueshifted lines from the wind were present. During the second 2019 September observation, three broad (σ = 3000 km s(exp -1)) absorption lines were resolved in the high-resolution Reflection Grating Spectrometer spectrum that are identified with blueshifted O VIII Lyα, Ne IX Heα, and Ne X Lyα. The outflow velocity of this soft X-ray absorber was found to be v/c = -0.258 ± 0.003, fully consistent with an iron K absorber with v/c = -0.261 ± 0.007. The ionization parameter and column density of the soft X-ray component (log ξ = 3.4, N(sub H) = 2 × 10(exp 21) cm(exp -2)) outflow was lower by about 2 orders of magnitude when compared to the high-ionization wind at iron K (log ξ = 5, N(sub H) = 7 × 10(exp 23) cm(exp -2). Substantial variability was seen in the soft X-ray absorber between the 2019 observations, declining from N(sub H) = 10(exp 23) to 10(exp 21) cm(exp -2) over 20 days, while the iron K component was remarkably stable. We conclude that the soft X-ray wind may originate from an inhomogeneous wind streamline passing across the line of sight that, due to its lower ionization, is located further from the black hole, on parsec scales, than the innermost disk wind

Active galactic nuclei↗

Analysis of CO2 Composite Spray Cleaning System Results

In a previous study, the cleaning efficiency of a CO2 composite cleaning system for particulate removal was tested. The study covered particles from spores to fluorescent particles of different sizes as well as a variety of substrate surfaces, including aluminum, titanium, stainless steel, and nitinol. Particles were deposited using aerosol (dry) and droplet (wet) deposition. Results from the previous study show that the CO2 composite spray system is capable of cleaning to sterility for aerosol deposited spores and is capable of cleaning a minimum of a 4-log reduction for droplet deposited spores. This minimum 4-log reduction matches current Planetary Protection dry heat microbial reduction requirements. In this paper we will present new data to further correlate the cleaning efficiency with how contamination was introduced to the surface, the surface roughness, and particle size. Possible causes for such correlations will be discussed.

Chung, Shirley↗

Spectral Range Within Global aCDOM(440) Algorithms for Oceanic, Coastal, and Inland Waters with Application to Airborne Measurements

The optically active component of dissolved organic material in aquatic ecosystems, or colored dissolved organic matter (CDOM), is represented by the coefficient of absorption due to the dissolved aquatic constituents at 440 nm, aCDOM(440). Remote sensing of aCDOM(440) enables characterization of ecosystem processes and aids in retrieval of chlorophyll a, a proxy for phytoplankton biomass. Spectrally adjacent band-ratio domains, e.g., blue to green, have previously been applied for remote sensing of aCDOM(440) in coastal and oceanic waters with similar results compared to more complex semi-analytical algorithms. Estimation of aCDOM(440) from ratios of the most spectrally separated ocean color wavebands (end members), e.g., ultraviolet (UV) to near-infrared (NIR), termed end-member analysis (EMA), has previously been shown to increase the accuracy of global aCDOM(440) retrievals from in-water observations of diffuse attenuation and to enable a unified algorithmic perspective without requiring regional adjustment of internal bio-optical parameters. EMA of above-water observations is evaluated herein, with a focus on coastal and inland waters in which increasing optical complexity and likelihood of bottom reflectance challenge the oceanic algorithms developed for deep and optically simple (case-1) waters. Analysis herein of three independent, in situ, bio-optical datasets indicates significant correlation between aCDOM(440) and end-member band ratios (next-generation 320 and 780 nm or legacy 412 and 670 nm ratios) with a coefficient of determination, R2, of 0.87 (log-scale) or higher based on a dataset spanning the dynamic range of global, conservative water bodies. For applicable wavelengths, EMA algorithms are shown to agree with case-1 relationships and to produce consistent log-scale uncertainties across more than three orders of magnitude in aCDOM(440) values (0.001–2.305 m􀀀1). EMA using UV and NIR wavelengths (320 and 780 nm) is applied to low-altitude airborne observations and satisfies 25% uncertainty based on unbiased percent differences (UPDs) within each of three dissimilar match-up sites ranging in aCDOM(440) from 0.02–0.57 m􀀀1. Results demonstrate that EMA is a useful and robust approach for the remote sensing of aCDOM(440) in coastal and inland waters, which are generally shallower, contain more optically complex environments, and span a greater range in aCDOM(440) than oceanic waters.

CDOM↗

Estimating the Ultraviolet Emission of M Dwarfs with Exoplanets from Ca ii and Hα

M dwarf stars are excellent candidates around which to search for exoplanets, including temperate, Earth-sized planets. To evaluate the photochemistry of the planetary atmosphere, it is essential to characterize the UV spectral energy distribution of the planet's host star. This wavelength regime is important because molecules in the planetary atmosphere such as oxygen and ozone have highly wavelength-dependent absorption cross sections that peak in the UV (900–3200 Å). We seek to provide a broadly applicable method of estimating the UV emission of an M dwarf, without direct UV data, by identifying a relationship between noncontemporaneous optical and UV observations. Our work uses the largest sample of M dwarf star far- and near-UV observations yet assembled. We evaluate three commonly observed optical chromospheric activity indices—Hα equivalent widths and log(10) L(Hα)/L(bol), and the Mount Wilson Ca ii H&K S and R'(HK) indices—using optical spectra from the HARPS, UVES, and HIRES archives and new HIRES spectra. Archival and new Hubble Space Telescope COS and STIS spectra are used to measure line fluxes for the brightest chromospheric and transition region emission lines between 1200 and 2800 Å. Our results show a correlation between UV emission-line luminosity normalized to the stellar bolometric luminosity and Ca ii R'(HK) with standard deviations of 0.31–0.61 dex (factors of ∼2–4) about the best-fit lines. We also find correlations between normalized UV line luminosity and Hα log(10) L(Hα)/L(bol) and the S index. These relationships allow one to estimate the average UV emission from M0 to M9 dwarfs when UV data are not available.

Katherine Melbourne↗

WISE 2150-7520AB: A Very Low-mass, Wide Comoving Brown Dwarf System Discovered through the Citizen Science Project Backyard Worlds: Planet 9

We report the discovery of WISE 2150-7520AB (W2150AB): a widely separated (∼341 au) very low-mass L1 + T8 co-moving system. The system consists of the previously known L1 primary 2MASS J21501592-7520367 and a newly discovered T8 secondary found at position 21:50:18.99–75:20:54.6 (MJD = 57947) using Wide-field Infrared Survey Explorer data via the Backyard Worlds: Planet 9 citizen science project. We present Spitzer ch1 and ch2 photometry (ch1–ch2 = 1.41 ± 0.04 mag) of the secondary and Folded-port InfraRed Echellette prism spectra of both components. The sources show no peculiar spectral or photometric signatures, indicating that each component is likely field age. Using all observed data and the Gaia DR2 parallax of 41.3593 ± 0.2799 mas for W2150A we deduce fundamental parameters of log(L(bol)/L(☉)) = −3.69 ± 0.01, T(eff) = 2118 ± 62 K, and an estimated mass = 72 ± 12 M(Jup) for the L1 and log(L(bol)/L(☉)) = −5.64 ± 0.02, T(eff) = 719 ± 61 K, and an estimated mass = 34 ± 22 M(Jup) for the T8. At a physical separation of ∼341 au this system has E(bin) = 1041 erg, making it the lowest binding energy system of any pair with M(tot) < 0.1 M(Sun) not associated with a young cluster. It is equivalent in estimated mass ratio, E(bin), and physical separation to the ∼2 Myr M7.25 + M8.25 binary brown dwarf 2MASS J11011926-7732383AB (2M1101AB) found in the Chameleon star-forming region. W2150AB is the widest companion system yet observed in the field where the primary is an L dwarf or later.

Jacqueline Faherty↗

Thermodynamics of the Lu2O3 – SiO2 system and comparison to other Rare Earth Silicates

Environmental barrier coatings are necessary to protect SiC based ceramics and composites from water vapor degradation in harsh engine environments. Currently, rare earth (RE) silicates are the most promising systems to protect SiC based ceramics and composites. This protection is largely based on reduced silica activity in these rare earth silicates which results in a lowered reactivity with water vapor. To that end, previous Knudsen effusion mass spectrometry (KEMS) studies have explored RE = Yb, Y silicates to measure the reduced silica activity and subsequent water vapor reactivity. Similarly, this work employs the KEMS technique to measure the SiO(g) vapor pressure in Lu containing RE silicates to calculate the activity of silica within the monosilicate [log(a(SiO2)) = -2351.1*1/T-1.6731] and the disilicate (log(a(SiO2)) = -4884.0*1/T + 2.2208) as a function of temperature. The enthalpies of formation for Lu monosilicate from the oxides and the elements were calculated to be -45±3 kJ/mol at 1550 K and -2831.1±12 kJ/mol at 298 K, respectively. The measured enthalpy of formation and those found in literature are compared to modeled values from density functional theory and those estimated using electronegativity.

Knudesen Effusion Mass Spectrometry↗

Comparison of U.S. and Russian Slow Crack Growth Data and Models

Both the US and Russian space programs use similar predictive models for design of fused silica windows on the International space station. The Russian model can be derived from the power expression for slow crack growth (SCG) or “static fatigue.” The US uses both power and exponential models. Despite the similarity of models and data fitting approach (linear regression and right censoring), different SCG parameters have been derived by US and Russian parties for the same material (Russian fused silica) tested in a similar manner. The difference appears to be related to the use of short-term strength data along with the longer-term static fatigue data, with the power law parameter n being very sensitive to the conversion of strength data into equivalent static data. This hybrid approach is feasible if strength data is measured with a constant stress rate and is appropriately converted to a static equivalent. More research into the approach is needed. However, because of the nonlinear behavior of fused silica in log(v) – log(K(I) ) space and the sensitivity of parameter estimation to fit range, the exponential model is a better choice regardless of test method. Functions are given to convert parameters from the Russian model to those in the US model. However, accurate conversion is hampered by the lack of inert strength data. When the same test technique is used, US and Russian materials exhibit very similar parameters.

Silica, strength, impact, crack growth, windows, I↗

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

Aerial Object Trajectory Classification by Training on Flight Controller Data and Testing on RADAR Generated Tracks

Onboard collision avoidance is needed to enable safe, autonomous flight operations for NASA projects such as Advanced Air Mobility (AAM), as well as many commercial applications. Real-time aerial object classification will improve onboard collision avoidance algorithm decision making and may reduce unnecessary activation of avoidance systems. This work trains an aircraft trajectory classifier using trajectories from flight controller logs and tests the classifier using RADAR collected trajectories during air to air experiments and ground to air experiments. In contrast to RADAR data, these flight controller logs are relatively abundant, which makes the possibility of substituting flight data for RADAR data an attractive, cost-effective option. The SVM model developed in this work achieved a 79.7% classification accuracy on the first second of radar trajectories of GA, multirotor sUAS, and fixed wing sUAS. Findings from this work indicate that it is feasible to classify sensor collected trajectories using a classifier trained on flight controller data.

Henry Holbrook↗

DataSet: Aerial Object Trajectory Classification by Training on Flight Controller Data and Testing on RADAR Generated

Onboard collision avoidance is needed to enable safe, autonomous flight operations for NASA projects such as Advanced Air Mobility (AAM), as well as many commercial applications. Real-time aerial object classification will improve onboard collision avoidance algorithm decision making and may reduce unnecessary activation of avoidance systems. This work trains an aircraft trajectory classifier using trajectories from flight controller logs and tests the classifier using RADAR collected trajectories during air to air experiments and ground to air experiments. In contrast to RADAR data, these flight controller logs are relatively abundant, which makes the possibility of substituting flight data for RADAR data an attractive, cost-effective option. The SVM model developed in this work achieved a 79.7% classification accuracy on the first second of radar trajectories of GA, multirotor sUAS, and fixed wing sUAS. Findings from this e that it is feasible to classify sensor collected trajectories using a classifier trained flight controller data.

Chester V Dolph↗

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↗

AdvoCATE User Guide

This tutorial gives an overview of AdvoCATE - the assurance case automation toolset - which provides a suite of features for safety risk management, including hazard analysis, the linking of hazards to requirement logs, risk modeling using bow tie diagrams, structured arguments, and evidence logs. In this tutorial we show how these related assurance artifacts can be combined with each other in an integrated assurance case.

Safety assurance↗

Role of RpoS in regulating stationary phase Salmonella Typhimurium pathogenesis-related stress responses under physiological low fluid shear force conditions

The discovery that biomechanical forces regulate microbial virulence was established with the finding that physiological low fluid shear (LFS) forces altered gene expression, stress responses, and virulence of the enteric pathogen Salmonella enterica serovar Typhimurium during log phase. These log phase LFS-induced phenotypes were independent of the master stress response regulator, RpoS (S). Given the central importance of RpoS in regulating stationary phase stress responses of S. Typhimurium cultured under conventional shake flask and static conditions, we examined its role in stationary phase cultures grown under physiological LFS. We constructed an isogenic rpoS mutant derivative of wild type S. Typhimurium and compared the ability of these strains to survive in vitro pathogenesis-related stresses that mimic those encountered in the infected host and environment. We also compared the ability of these strains to colonize (adhere, invade, and survive within) human intestinal epithelial cell cultures. Unexpectedly, LFS-induced resistance of stationary phase S. Typhimurium cultures to acid and bile salts stresses did not rely on RpoS. Likewise, RpoS was dispensable for stationary phase LFS cultures to adhere to and survive within intestinal epithelial cells. In contrast, the resistance of these cultures to challenge with oxidative and thermal stresses, and their invasion into intestinal epithelial cells was influenced by RpoS. These findings expand our mechanistic understanding of how physiological fluid shear forces modulate stationary phase S. Typhimurium physiology in unexpected ways and provide clues into microbial mechanobiology and nuances of Salmonella responses to microenvironmental niches in the infected host.

Karla Franco Meléndez↗

Quantum-Accelerated Distributed Algorithms for Approximate Steiner Trees and Directed Minimum Spanning Trees

We present two algorithms in the Quantum CONGEST-CLIQUE model of distributed computation that succeed with high probability; one for producing an approximately optimal Steiner Tree, and one for producing an exact spanning arborescence of minimum weight, the analog of a Minimum Spanning Tree in a directed graph, each of which uses O~(n^(1/4)) rounds of communication and O~(n^(9/4)) messages, achieving a lower round and message complexity than any known algorithms in the classical CONGEST-CLIQUE model. The CONGEST distributed computational model allows limited-sized messages to be transmitted within a network described by a communication graph of size n in a series of rounds to address a computational problem. The size limitation for such messages isO(log(n)) bits at each edge of the communication graph per round. The communication graph in the CONGEST-CLIQUE model is fully connected. In the Quantum CONGEST-CLIQUE model, at most O(log(n)) classical and quantum bits (qubits) can be communicated across each edge of the communication graph per round. At a high level, we achieve these results by combining classical algorithms with fast quantum subroutines. These speedups further contribute to understanding what problems can be solved more efficiently when we allow quantum communication in this CONGEST-CLIQUE model of distributed computation.

quantum distributed algorithms↗