Search NASA⌕ Search

SEARCH · Search NASA

Results for “outer 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.

At least 703 records · Page 39

Crater Topography on Titan: Implications for Landscape Evolution

We present a comprehensive review of available crater topography measurements for Saturn's moon Titan. In general, the depths of Titan's craters are within the range of depths observed for similarly sized fresh craters on Ganymede, but several hundreds of meters shallower than Ganymede's average depth vs. diameter trend. Depth-to-diameter ratios are between 0.0012 +/- 0.0003 (for the largest crater studied, Menrva, D approximately 425 km) and 0.017 +/- 0.004 (for the smallest crater studied, Ksa, D approximately 39 km). When we evaluate the Anderson-Darling goodness-of-fit parameter, we find that there is less than a 10% probability that Titan's craters have a current depth distribution that is consistent with the depth distribution of fresh craters on Ganymede. There is, however, a much higher probability that the relative depths are uniformly distributed between 0 (fresh) and 1 (completely infilled). This distribution is consistent with an infilling process that is relatively constant with time, such as aeolian deposition. Assuming that Ganymede represents a close 'airless' analogue to Titan, the difference in depths represents the first quantitative measure of the amount of modification that has shaped Titan's surface, the only body in the outer Solar System with extensive surface-atmosphere exchange.

Implications↗

Development of a Model for the Small-Particle Orbital Debris Population Based on the STS Impact Record

In preparation for the release of the Orbital Debris Engineering Model (ORDEM) version 3.1, the NASA Orbital Debris Program Office (ODPO) revisited how orbiting debris populations of characteristic sizes smaller than 1 cm were modeled. The primary contributor to the population of sub-centimeter debris particles is the surface deterioration or erosion of spacecraft materials exposed to the outer-space environment. Because small particulates are not directly trackable by remote sensing, the primary means of detection is via historical counts of small impact features on flown radiator and window surfaces of the U.S. Space Transportation System (STS, also known as the Space Shuttle) from 1995-2011. Historic NASA studies of high-velocity impact tests have related impact-feature size to particle mass and velocity for certain STS surfaces, so that a corresponding particle size may be inferred from each small-impact feature observed. Micro-debris populations are then estimated by modeling the path and orientation of an STS mission through a simulated debris environment, and the densities of this simulated environment are rescaled to approximate the number of observed STS impact features. Monte-Carlo methods are further employed to gauge the estimation uncertainty of the rescaled environment. A description of the chosen methodologies for estimating and adjusting the micro-debris population model, and the results, are presented.

Seago, J..↗

Hot proton anisotropies and cool proton temperatures in the outer magnetosphere

The hot protons of the outer magnetosphere typically exhibit a temperature anistropy such that T(sub perp)/T(sub parallel) greater than 1, where perpendicular and parallel symbols denote directions relative to the background magnetic field. If this anisotrpy is sufficiently large, the electomagneitc proton cyclotron anistropy instability may be excited. This instability is studied using linear Vlasov theory and one-dimensional hybrid simulations for a homogeneous plasma model representative of conditions in the outer magnetosphere with a hot anisotropic proton component (denoted by subscript h) and a cool, initially isotropic proton component (subscript c). Linear theory yields an instability threshold condition on the hot proton temperature anistropy where as the simulations imply an upper bound on T(sub perp h)/T(sub parallel h); both the threshold and the upper bound have similar scaling with the maximum growth rate gamma (sub m), the parallel beta of the hot component, beta(sub parallel h), and the relative density of the hot component n(sub h)/n(sub e). An anlysis of plasma observations from the Los Alamos magnetospheric plasma analyzer (MPA) in geosynchronous orbits finds that the maximum value of the hot proton temperature anisotropy approximately satisfies the predicted scaling with beta(sub parallel h) and nu(sub h)/n(sub e) and yields the proportionality factor that quantifies this upper bound. The simulations are also used to examine the heating of the cool proton cyclotron instability. The simulations yield a scaling for the dimensionless late-time cool proton average temperature T(sub c)/T(sub parallel h) as (n(sub h)/n(sub e))/beta(sub parallel h exp 0.5). Analysis of MPA data shows that the observed values of T(sub c)/T(sub parallel h) have similar scaling and again yield the proportionality factor which quantifies this relationship.

Gary, S. Peter↗

Electric Sail Propulsion for Exploring Nearby Interstellar Space

An Electric Sail is a revolutionary propellant-less propulsion system that is ideal for deep space missions to the outer planets, the Heliopause, and beyond. It is revolutionary in that it uses momentum exchange with the hypersonic solar wind to propel a spacecraft within the heliosphere. The momentum exchange is affected by the deflection of charged solar wind particles by an array of electrically biased wires that extend outward up to 30 km from a slowly rotating spacecraft. A high-voltage, positive bias on the wires, which are oriented normal to the solar wind flow, deflects the streaming protons, resulting in a reaction force on the wires that is also directed radially away from the sun. Over a period of months, this small force can accelerate the spacecraft to enormous speeds-on the order of 100-150 km/s (approximately 20 to 30 AU/yr). Unlike solar sails, Electric Sails do not rely on a fixed area to produce thrust. In fact, as they move away from the Sun and solar wind pressure decreases, the area for solar proton momentum transfer becomes larger, increasing system efficiency. As a result, thrust decreases at ≈1/r**(7/6) instead of the ≈1/r**2 rate typical for solar sails. The net effect is that an increased radial range of operation, together with increased thrust, both contribute to higher velocities and shorter total trip times to distant destinations. The MSFC Advanced Concepts Office (ACO) was awarded a Phase II NASA Innovative Advanced Concepts (NIAC) study to mature the technology for possible future demonstration and implementation. Preliminary results indicate that the physics of the system is viable and that a spacecraft propelled by an Electric Sail could reach the Heliopause in less than 15 years - and could be developed within a decade.

Johnson, Les↗

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗

Realistic NLTE Radiative Transfer for Modeling Stellar Winds

This NASA grant supported the development of codes to solve the non-LTE multi-level spherical radiative transfer problem in the presence of velocity fields. Much of this work was done in collaboration with Graham Harper (CASA, University of Colorado). These codes were developed for application to the cool, low-velocity winds of evolved late-type stars. Particular emphasis was placed on modeling the wind of lambda Velorum (K4 lb), the brightest K supergiant in the sky, based on extensive observations of the ultraviolet spectrum with the HST/GHRS from GO program 5307. Several solution techniques were examined, including the Eddington factor Approach described in detail by Bennett & Harper (1997). An Eddington factor variant of Harper's S-MULTI code (Harper 1994) for stationary atmospheres was developed and implemented, although full convergence was not realized. The ratio of wind terminal velocity to turbulent velocity is large (approx. 0.3-0.5) in these cool star winds so this assumption of stationarity provides reasonable starting models. Final models, incorporating specified wind laws, were converged using the comoving CRD S-MULTI code. Details of the solution procedure were published by Bennett & Harper (1997). Our analysis of the wind of lambda Vel, based on wind absorption superimposed on chromospheric emission lines in the ultraviolet, can be found in Carpenter et al. (1999). In this paper, we compare observed wind absorption features to an exact CRD calculation in the comoving frame, and also to a much quicker, but approximate, method using the SEI (Sobolev with Exact Integration) code of Lamers, Cerruti-Sola, & Perinotto (1987). Carpenter et al. (1999) provide detailed comparisons of the exact CRD and approximate SEI results and discuss when SEI is adequate to use for computing wind line profiles. Unfortunately, the observational material is insufficient to unambiguously determine the wind acceleration law for lambda Vel. Relatively few unblended Fe II lines of optical depth sensitive to the wind acceleration region are present in the existing HST/GHRS data set. Most of the Fe II lines are either too optically thick (resulting in a board, black wind absorption profile) or too optically thin (resulting in no wind absorption feature present). Also, most of the ultraviolet spectra obtained from HST GO-5307 was at medium resolution (R approx. 40,000, corresponding to a velocity resolution of 7.5 km/s). This is simply inadequate to resolve the turbulence in the outer wind; a key parameter in theoretical wind models. We can now say that an unambiguous determination of the wind velocity law in lambda Vel will require complete coverage of the ultraviolet spectrum at high dispersion (R approx. 10(exp 5), or 3 km/s). This is now feasible usin, the STIS echelle spectrograph on-board HST.

Bennett, Philip D.↗

Performance and Application of Parallel OVERFLOW Codes on Distributed and Shared Memory Platforms

The presentation discusses recent studies on the performance of the two parallel versions of the aerodynamics CFD code, OVERFLOW_MPI and _MLP. Developed at NASA Ames, the serial version, OVERFLOW, is a multidimensional Navier-Stokes flow solver based on overset (Chimera) grid technology. The code has recently been parallelized in two ways. One is based on the explicit message-passing interface (MPI) across processors and uses the _MPI communication package. This approach is primarily suited for distributed memory systems and workstation clusters. The second, termed the multi-level parallel (MLP) method, is simple and uses shared memory for all communications. The _MLP code is suitable on distributed-shared memory systems. For both methods, the message passing takes place across the processors or processes at the advancement of each time step. This procedure is, in effect, the Chimera boundary conditions update, which is done in an explicit "Jacobi" style. In contrast, the update in the serial code is done in more of the "Gauss-Sidel" fashion. The programming efforts for the _MPI code is more complicated than for the _MLP code; the former requires modification of the outer and some inner shells of the serial code, whereas the latter focuses only on the outer shell of the code. The _MPI version offers a great deal of flexibility in distributing grid zones across a specified number of processors in order to achieve load balancing. The approach is capable of partitioning zones across multiple processors or sending each zone and/or cluster of several zones into a single processor. The message passing across the processors consists of Chimera boundary and/or an overlap of "halo" boundary points for each partitioned zone. The MLP version is a new coarse-grain parallel concept at the zonal and intra-zonal levels. A grouping strategy is used to distribute zones into several groups forming sub-processes which will run in parallel. The total volume of grid points in each group are approximately balanced. A proper number of threads are initially allocated to each group, and in subsequent iterations during the run-time, the number of threads are adjusted to achieve load balancing across the processes. Each process exploits the multitasking directives already established in Overflow.

Djomehri, M. Jahed↗

Decadal Polar Motion of the Earth Excited by the Convective Outer Core From Geodynamo Simulations

Long time geodetic observation records show that the orientation of the Earth's rotation axis with respect to the terrestrial reference frame, or polar motion, changes on a broad range of timescales. Apart from external torques from the luni-solar tides, these changes are excited by interactions among different components of the Earth system. The convective fluid outer core has long been conjectured a likely contributor to the observed polar motion on timescales upward of decades, such as the 30 year Markowitz wobble. We investigated the electromagnetic coupling scenario across the core-mantle boundary via numerical geodynamo simulation for different geodynamo parameters (Rayleigh numbers and magnetic Rossby numbers). Our simulated polar motion varies strongly with the dynamo parameters, while its excitation on decadal timescales appear to converge asymptotically within the adopted range of numerical Rossby numbers. Three strongest asymptotic modes emerge from numerical results, with periods around 30, 40, and 60 years for the prograde excitation and around 24, 30, and 60 years for the retrograde excitation. Their amplitudes are all larger than 5 times 10 (sup minus 8), or approximately 10 milliseconds of arc.The results suggest that the electromagnetic core-mantle coupling could explain a substantial portion, if not all, of the observed decadal polar motion. In particular, the predicted 60 year polar motion deserves special attention for future observations and studies.

Kuang, W.↗

Studies of asteroids, comets, and Jupiter's outer satellites

The work comprises observational, theoretical, and computational research on asteroids, together with a smaller effort concerning the astrometry of comets and Jupiter's satellites JVI through JXIII. Two principal areas of research, centering on astrometry and photometry, are interrelated in their aim to study the overall structure of the asteroid belt and the physical and orbital properties of individual asteroids. About 2000 accurate photographic positions of asteroids and comets, including a number from the Lowell, Palomar, and Goethe-Link archival plate collections, the last of which was donated to us last winter by Indiana University were measured and published. Charge coupled device (CCD) astrometry of 36 faint targets was undertaken, including 4 comets; JVI, JVII, JVIII, JLX, JXI, and JXII; and 26 asteroids, most of which are Earth-approachers. A deep, bias-correctable asteroid survey (LUKAS), the aim of which is to determine the true spatial distribution of asteroids down to subkilometer diameters was started. A series of eight plates at the UK Schmidt telescope that contain images of asteroids as faint as V approximately 22 mag was obtained. Analysis of microdensitometric scans of two plates has shown that about 98 percent of the asteroid images could be identified completely automatically.

Bowell, Edward↗

High-luminosity single carbon stars in stellar and galactic evolution

In the solar neighborhood, approximately half of all intermediate mass main sequence stars with initially between 1 solar mass and about 5 solar masses become carbon stars with luminosities near 10,000 lunar luminosities for typically less than 1 million years. These high luminosity carbon stars lose mass at rates nearly always in excess of 10 to the -7th solar mass/yr and sometimes in excess of 0.00001 solar mass/yr. Locally, close to half of the mass returned into the interstellar medium by intermediate mass stars before they become white dwarfs is during the carbon star phase. A much greater fraction of lower metallicity stars become carbon-rich before they evolve into planetary nebulae, than do higher-metallicity stars; therefore, carbon stars are much more important in the outer than in the inner Galaxy.

Jura, M.↗

Large-Aperture [O I] 6300 A Photometry of Comet Hale-Bopp: Implications for the Photochemistry of OH

Large-aperture photometric observations of comet Hale-Bopp (C/1995 O1) in the forbidden red line of neutral oxygen ([O I] 6300 angstroms) with the 150 mm dual-etalon Fabry-Perot spectrometer that comprises the Wisconsin H-alpha Mapper and a 50 mm dual-etalon Fabry-Perot spectrometer at the McMath-Pierce main telescope from 1997 late February to mid April yield a total metastable O((sup 1)D) production rate of (2.3-5.9) x 10(exp 30)/s. Applying the standard H2O and OH photodissociation branching ratios, we derive a water production rate, Q(H2O), of (2.6-6.1) x 10(exp 31)/s, which disagrees with Q(H2O = 1x10(exp 31)/s determined by independent H2O, OH, and H measurements. Furthermore, our own [O I] 6300 observations of the inner coma (< 30,000 km) using the 3.5 m Wisconsin-Indiana-Yale-NOAO telescope Hydra and Densepak multi-object spectrographs yield Q(H2O) = 1 x 10(exp 31)/s. Using our [O I] 6300 data, which cover spatial scales ranging from 2,000 to 1x10(exp 6) km, and a complementary set of wide-field ground-based OH images, we can constrain the sources of the apparent excess O((sup 1)D) emission to the outer coma, where photodissociation of OH is assumed to be the dominant O((sup 1)D) production mechanism. From production rates of other oxygen-bearing volatiles (e.g., CO and CO2), we can account for at most 30% of the observed excess O((sup 1)D) emission. Since even less O((sup 1)D) should be coming from other sources (e.g., electron excitation of neutral O and distributed nonnuclear sources of H2O), we hypothesize that the bulk of the excess O((sup 1)D) is likely coming from photodissociating OH. Using the experimental OH photo-dissociation cross section of Nee and Lee at Ly-alpha as a guide in modifying the theoretical OH cross sections of van Dishoeck and Dalgarno, we can account for approximately 60% of the observed O((sup 1)D) excess without requiring major modifications to the other OH branching ratios or the total OH photodissociation lifetime.

Morgenthaler, Jeffrey P.↗

Outer magnetosphere near midnight at quiet and disturbed times.

Ogo 5 magnetic-field and energetic-electron (E greater than 50 keV) data are used to study both the quiet-time, steady-state configuration of the outer magnetosphere or near tail region near midnight and the disturbed time changes of this configuration. The nighttime cusp is found to be a distinct feature within the plasma sheet at quiet times but indistinguishable from the plasma sheet at disturbed times. The sequence of thinning and expansion of the plasma sheet in this region in association with the substorms is studied. The response of the plasma sheet in the near tail at about 10 earth radii is found to be similar to that in the more distant tail at more than 20 earth radii. Finally, the nature of field-aligned currents flowing on the plasma-sheet boundary is investigated. Assuming infinite current sheets, the sheet current density at Ogo 5 is found to be approximately .01 A/m.

Aubry, M. P.↗

Helium abundances on the moon: Assumptions and estimates

Nuclear energy is a highly desirable source of energy, and He-3 is the most prized of the fusion reactants. As the Wisconsin Group has emphasized, He-3 may be the only true economic ore on the Moon. The lack of a shielding atmosphere on the Moon permits solar-wind alpha particles to impinge upon the lunar regolith and become implanted into the various solid components. In particular, large quantities of helium (5 to 50 ppm) are presented. The measured parameter of I(sub s)/FeO, a direct indicator of maturity and exposure age, can be used as a first approximation to predict the abundances of many solar-wind components in the soils. However, because ilmenite has a much higher retentivity for helium than the other phases, the TiO2 contents of the soils are better indicators of helium contents (Taylor, Space 90). High-Ti mare bassalt regions, such as at the Apollo 17 locale, appear to be the best areas for He mining (15 to 50 ppm He(sub T)), versus 3 to 9 ppm in the Highlands. However, the relationships between I(sub s)/FeO, TiO2 and He-3 contents are complicated - e.g., many of the most He-rich soils are immature to submature. The amount of He-3 in the regolith of the moon is estimated at 220,000 tons in the outer 2 m of the Maria.

Taylor, Lawrence A.↗

Optical Design of the Origins Space Telescope

The optical design is presented for the Origins Space Telescope, one of four large missions under study in preparation for the 2020 Decadal Survey in Astronomy and Astrophysics is presented. Sensitive to the mid- and far-infrared spectrum (between 2.8 and 588 µm) and set to launch in approximately 2035, Origins sets out to answer a number of important scientific questions by addressing NASA's three key science goals in astrophysics: How does the universe work?, How did we get here, and Are we alone? To do so, Origins seeks to answer three main questions: 1) How do galaxies form stars, build up metals, and grow their central supermassive black holes from reionization to today? 2) How do the conditions for habitability develop during the process of planet formation? 3) Do planets orbiting M-dwarf stars support life?The Origins telescope has a 5.9 m diameter primary mirror and operates at f/14. The large on-axis primary consists of 18 'keystone' segments of two different prescriptions arranged in two annuli (six inner and twelve outer segments) that together form a circular aperture in the goal of achieving a symmetric point spread function. To simplify the design, the telescope will be launched with both the primary mirror segments and the secondary mirror already deployed. To accommodate the 46 x 15 arcminute full field of view of the telescope at the design wavelength of = 30 µm, a three-mirror anastigmat configuration is used. The design is diffraction-limited across its instruments' fields of view. A brief discussion of each of the three baselined instruments within the Instrument Accommodation Module (IAM) is presented: 1) Origins Survey Spectrometer (OSS), 2) Mid-infrared Spectrometer, Camera (MISC) transit spectrometer channel, and 3) Far-Infrared Polarimeter/Imager (FIP). In addition, the upscope options for the observatory are laid out as well, including a fourth instrument: the Heterodyne Receiver for Origins (HERO).

Howard, Joseph↗

Feasibility test for a V-slit star mapper for pioneer spacecraft terminal navigation

A laboratory demonstration of the feasibility of using a V-slit star mapper to meet the sensitivity and accuracy of on-board navigational requirements for future Pioneer Missions to the outer planets was conducted by the Control and Sensors Laboratory of TRW. The breadboard was extremely simple in configuration, consisting of an end-on photomultiplier tube and a V-slit reticle located at the focal plane of the objective lens. In addition, a plano-convex lens was used between the reticle and the PMT in a Fabry-Perot configuration. The analytical effort indicated that the sensor should easily meet the requirements. The Pioneer SRA test set was examined to determine its basic accuracy and modify it where necessary to bring its accuracy into the 1-3 arc second range. The test results show that it is feasible to use this type of star mapper in the 10 arc second accuracy range. The test equipment accuracy (approximately 5 arc Sec) was sufficient to bound the sensor errors at less than 10 arc seconds.

Gates, R. F.↗

The Galileo Probe: How it Has Changed Our Understanding of Jupiter

The Galileo Mission to Jupiter, which arrived in December of 1995, provided the first study by an orbiter, and the first in-situ sampling via an entry probe, of an outer planet atmosphere. The rationale for an entry probe is that, even from an orbiter, remote sensing of the jovian atmosphere could not adequately retrieve the information desired. This paper provides a current summary of the most significant aspects of the data returned from the Galileo entry probe. As a result of the probe measurements, there has been a reassessment of our understanding of outer planet formation and evolution of the solar system. The primary scientific objective of the Galileo probe was to determine the composition of the jovian atmosphere, which from remote sensing remained either very uncertain, or completely unknown, with respect to several key elements. The probe found that the global He mass fraction is. significantly above the value reported from the Voyager Jupiter flybys but is slightly below the protosolar value, implying that there has been some settling of He to the deep jovian interior. The probe He measurements have also led to a reevaluation of the Voyager He mass fraction for Saturn, which is now determined to be much closer to that of Jupiter. The elements C, N, S, Ar, Kr, Xe were all found to have global abundances approximately 3 times their respective solar abundances. This result has raised a number of fundamental issues with regard to properties of planetesimals and the solar nebula at the time of giant planet formation. Ne, on the other hand, was found to be highly depleted, probably as the result of it being carried along with helium as helium settles towards the deep interior. The global abundance of O was not obtained by the probe because of the influence of local processes at the probe entry site (PES), processes which depleted condensible species, in this case H2O, well below condensation levels. Other condensible species, namely NH3 and H2S, were similarly affected but attained their deep equilibrium mixing ratios before the maximum depth sampled by the probe. Processes that might be capable of producing such effects on the condensibles are still under investigation. Measured isotopic ratios of noble gases and other heavy elements are solar, and (D + (Sup 3)He)/H is the same to within measurement uncertainties as in the local interstellar medium. No thick clouds were detected, and in particular no significant water cloud, but the PES location clearly affected the probe measurements of clouds. In fact, the probe data must be understood in the context of the location of the PES, which was within what is termed a 5 micron hot spot, a local clearing in the clouds that is bright near the 5 microns spectral region. The thermal structure at the PES was determined from approximately 1000 km above the 1 bar pressure level (10(exp -9 bars)) to 132 km 1 bar (22bars). The probe showed the atmosphere to have a generally sub-adiabatic temperature gradient (static stability) of = 0.1 K/km to as deep as the probe made measurements. In the upper atmosphere the probe derived a maximum positive vertical temperature gradient of approximately 5 K/km, and maximum temperature of = 900 K. The energy sources producing the warm upper atmosphere have yet to be completely identified. At first glance, Doppler tracking of the probe indicates that the long observed cloud level zonal winds extend to levels at least as deep as the probe made measurements. Zonal wind increases from = 80 m/s at pressures less than a bar to about 180 m/s near 5 bars, and remains approximately constant with depth thereafter. However, there is a question as to whether the winds measured from probe tracking are representative of the general wind field, or are considerably influenced by localized winds associated with the PES.

Young, Richard E.↗

A Resolved Near-Infrared Image of the Inner Cavity in the GM Aur Transitional Disk

We present high-contrast H-band polarized intensity (PI) images of the transitional disk around the young solar like star GM Aur. The near-infrared direct imaging of the disk was derived by polarimetric differential imaging using the Subaru 8.2 m Telescope and HiCIAO. An angular resolution and an inner working angle of 0 07 and radius approximately 0 05, respectively, were obtained. We clearly resolved a large inner cavity, with a measured radius of 18+/ 2 au, which is smaller than that of a submillimeter interferometric image (28 au). This discrepancy in the cavity radii at near-infrared and submillimeter wavelengths may be caused by a 34M(sub Jup) planet about 20 au away from the star, near the edge of the cavity. The presence of a near-infrared inner cavity is a strong constraint on hypotheses for inner cavity formation in a transitional disk. A dust filtration mechanism has been proposed to explain the large cavity in the submillimeter image, but our results suggest that this mechanism must be combined with an additional process. We found that the PI slope of the outer disk is significantly different from the intensity slope obtained from HSTNICMOS, and this difference may indicate the grain growth process in the disk.

high-contrast H-band polarized intensity (PI)↗

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↗