Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equation discovery”

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 91 records · Page 5

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of very strong gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Recently this situation has changed dramatically, with a series of amazing breakthroughs. This talk will take you on this quest for the holy grail of numerical relativity, showing how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LIGO and LISA.

Centrella, Joan↗

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of extreme gravity, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. This situation has changed dramatically in the past year, with a series of amazing breakthroughs. This talk will take you on this quest for the holy grail of numerical relativity, showing how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LISA and LIGO.

Centrella, Joan↗

Crossing Mars: Past and Future Missions to a Cold, Dry Desert

Dr. Geoffrey A. Landis of the Photovoltaics and Space Environment Effects Branch presented an overview of recent discoveries about the environment of Mars. He covered missions from the 1966 Mariner IV that returned those first grainy close-up pictures of Mars showing an ancient cratered terrain to the Mars Odyssey mission with its tantalizing evidence of recent water flows on Mars. Mars is one of the most interesting planets in the solar system, featuring enormous canyons, giant volcanoes, and indications that, early in its history, it might have had rivers and perhaps even oceans. Five years ago, in July of 1997, the Pathfinder mission landed on Mars, bringing with it the microwave-oven sized Sojourner rover to wander around on the surface and analyze rocks. Pathfinder is only the first of an armada of spacecraft that will examine Mars from the pole to the equator in the next decade, culminating (someday, we hope!) with a mission to bring humans to Mars.

Landis, Geoffrey A.↗

Complexity for Survival of Living Systems

A logical connection between the survivability of living systems and the complexity of their behavior (equivalently, mental complexity) has been established. This connection is an important intermediate result of continuing research on mathematical models that could constitute a unified representation of the evolution of both living and non-living systems. Earlier results of this research were reported in several prior NASA Tech Briefs articles, the two most relevant being Characteristics of Dynamics of Intelligent Systems (NPO- 21037), NASA Tech Briefs, Vol. 26, No. 12 (December 2002), page 48; and Self-Supervised Dynamical Systems (NPO- 30634) NASA Tech Briefs, Vol. 27, No. 3 (March 2003), page 72. As used here, living systems is synonymous with active systems and intelligent systems. The quoted terms can signify artificial agents (e.g., suitably programmed computers) or natural biological systems ranging from single-cell organisms at one extreme to the whole of human society at the other extreme. One of the requirements that must be satisfied in mathematical modeling of living systems is reconciliation of evolution of life with the second law of thermodynamics. In the approach followed in this research, this reconciliation is effected by means of a model, inspired partly by quantum mechanics, in which the quantum potential is replaced with an information potential. The model captures the most fundamental property of life - the ability to evolve from disorder to order without any external interference. The model incorporates the equations of classical dynamics, including Newton s equations of motion and equations for random components caused by uncertainties in initial conditions and by Langevin forces. The equations of classical dynamics are coupled with corresponding Liouville or Fokker-Planck equations that describe the evolutions of probability densities that represent the uncertainties. The coupling is effected by fictitious information-based forces that are gradients of the information potential, which, in turn, is a function of the probability densities. The probability densities are associated with mental images both self-image and nonself images (images of external objects that can include other agents). The evolution of the probability densities represents mental dynamics. Then the interaction between the physical and metal aspects of behavior is implemented by feedback from mental to motor dynamics, as represented by the aforementioned fictitious forces. The interaction of a system with its self and nonself images affords unlimited capacity for increase of complexity. There is a biological basis for this model of mental dynamics in the discovery of mirror neurons that learn by imitation. The levels of complexity attained by use of this model match those observed in living systems. To establish a mechanism for increasing the complexity of dynamics of an active system, the model enables exploitation of a chain of reflections exemplified by questions of the form, "What do you think that I think that you think...?" Mathematically, each level of reflection is represented in the form of an attractor performing the corresponding level of abstraction with more details removed from higher levels. The model can be used to describe the behaviors, not only of biological systems, but also of ecological, social, and economics ones.

Zak, Michail↗

The Scintillation Prediction Observations Research Task (SPORT): an International Science Mission Using a Cubesat

The Scintillation Prediction Observations Research Task (SPORT) is a 6U CubeSat mission to address the compelling but difficult problem of understanding the preconditions leading to equatorial plasma bubbles. The scientific literature describes the preconditions in both the plasma drifts and the density profiles related to bubble formations that occur several hours later in the evening. Most of the scientific discovery has resulted from observations at a single site, within a single longitude sector, from Jicamarca, Peru. SPORT will provide a systematic study of the state of the pre-bubble conditions at all longitudes sectors to enhance understanding between geography and magnetic geometry. SPORT is an international partnership between National Aeronautics and Space Administration (NASA), the Brazilian National Institute for Space Research (INPE), and the Technical Aeronautics Institute under the Brazilian Air Force Command Department (DCTA/ITA), and encouraged by U.S. Southern Command. This talk will present an overview of the SPORT mission, observation strategy, and science objectives to improve predictions of ionospheric disturbances that affect radio propagation of telecommunication signals. The science goals will be accomplished by a unique combination of satellite observations from a nearly circular middle inclination orbit and the extensive operation of ground based observations from South America near the magnetic equator.

Scintillation Prediction Observations↗

Hypersonic Navier-Stokes Comparisons to Orbiter Flight Data

During the STS-119 flight of Space Shuttle Discovery, two sets of surface temperature measurements were made. Under the HYTHIRM program3 quantitative thermal images of the windward side of the Orbiter with a were taken. In addition, the Boundary Layer Transition Flight Experiment 4 made thermocouple measurements at discrete locations on the Orbiter wind side. Most of these measurements were made downstream of a surface protuberance designed to trip the boundary layer to turbulent flow. In this paper, we use the US3D computational fluid dynamics code to simulate the Orbiter flow field at conditions corresponding to the STS-119 re-entry. We employ a standard two-temperature, five-species finite-rate model for high-temperature air, and the surface catalysis model of Stewart.1 This work is similar to the analysis of Wood et al . 2 except that we use a different approach for modeling turbulent flow. We use the one-equation Spalart-Allmaras turbulence model8 with compressibility corrections 9 and an approach for tripping the boundary layer at discrete locations. In general, the comparison between the simulations and flight data is remarkably good

Candler, Graham V.↗

Cratering Rates in the Outer Solar System

We use several independent constraints on the number of ecliptic comets (aka JFCs) to determine impact cratering rates from Jupiter to Pluto. Long period comets and asteroids are currently unimportant on most worlds at most sizes. The size- number distribution of comets smaller than 20 km is inferred from size-number distributions of impact craters on Europa, Ganymede, and Triton; while the size- number distribution of comets bigger than 50 km is equated to the size-number distribution of Kuiper Belt Objects. The gap is bridged by interpolation. It is notable that small craters on Jupiter's moons indicate a pronounced paucity of small impactors, while small craters on Triton imply a collisional population rich in small bodies. However it is unclear whether the craters on Triton are of heliocentric or planetocentric origin. We therefore consider two cases for Saturn and beyond: a Case A in which the size-number distribution is like that inferred at Jupiter, and a Case B in which small objects obey a more nearly collisional distribution. Known craters on Saturnian and Uranian satellites are consistent with either Case, although surface ages are much younger in Case B, especially at Saturn and Uranus. At Neptune and especially at Saturn our cratering rates are much higher than rates estimated by Shoemaker and colleagues, presumably because Shoemaker's estimates mostly predate discovery of the Kuiper Belt. We also estimate collisional disruption rates of moons and compare these to estimates in the literature .

Zahnle, K.↗

Positron Interactions with Atoms and Ions

Dirac, in 1928, combining the ideas of quantum mechanics and the ideas of relativity invented the well-known relativistic wave equation. In his formulation, he predicted an antiparticle of the electron of spin n-bar/2. He thought that this particle must be a proton. Dirac published his interpretation in a paper 'A theory of electrons and protons.' It was shown later by the mathematician Hermann Weyl that the Dirac theory was completely symmetric between negative and positive particles and the positive particle must have the same mass as that of the electron. In his J. Robert Oppenheimer Memorial Prize Acceptance Speech, Dirac notes that 'Blackett was really the first person to obtain hard evidence for the existence of a positron but he was afraid to publish it. He wanted confirmation, he was really over cautious.' Positron, produced by the collision of cosmic rays in a cloud chamber, was detected experimentally by Anderson in 1932. His paper was published in Physical Review in 1933. The concept of the positron and its detection were the important discoveries of the 20th century. I have tried to discuss various processes involving interactions of positrons with atoms and ions. This includes scattering, bound states and resonances. It has not been possible to include the enormous work which has been carried out during the last 40 or 50 years in theory and measurements.

Bhatia, Anand K.↗

A Study of the Adequacy of Quasi-geostrophic Dynamics for Modeling the Effect of Frontal Cyclones on the Larger Scale Flow

The major objectives of this study are to test the validity of quasi-geostrophic (QG) dynamics, compared to primitive equation (PE) dynamics for modeling the effect of cyclone waves on the larger scale flow, and to study the formation of frontal cyclones and the dynamics of occluded frontogenesis. In order to extend the realism of the initial conditions over those used in earlier work of this nature and in order to include strong surface frontal zones in the initial conditions, a horizontal resolution on the order of 100 km between gridpoints is needed. In order to allow for integrations of up to two weeks simulated time and to keep the logistics and cost of the runs feasible, the horizontal domain of the models must be reduced. Based on work concerning polar lows, the tropospheric static stability was reduced and was also allowed to vary in different tropospheric layers. Previous work used a 'standard atmosphere' temperature lapse rate. Stratospheric stability remains unchanged. The use of one of these 'reduced stability' cases, with no surface frontal zone present resulted in the discovery of normal mode perturbations.

Mudrick, S.↗

Kepler and Ground-Based Transits of the exo-Neptune HAT-P-11b

We analyze 26 archival Kepler transits of the exo-Neptune HAT-P-11b, supplemented by ground-based transits observed in the blue (B band) and near-IR (J band). Both the planet and host star are smaller than previously believed; our analysis yields Rp = 4.31 R xor 0.06 R xor and Rs = 0.683 R solar mass 0.009 R solar mass, both about 3 sigma smaller than the discovery values. Our ground-based transit data at wavelengths bracketing the Kepler bandpass serve to check the wavelength dependence of stellar limb darkening, and the J-band transit provides a precise and independent constraint on the transit duration. Both the limb darkening and transit duration from our ground-based data are consistent with the new Kepler values for the system parameters. Our smaller radius for the planet implies that its gaseous envelope can be less extensive than previously believed, being very similar to the H-He envelope of GJ 436b and Kepler-4b. HAT-P-11 is an active star, and signatures of star spot crossings are ubiquitous in the Kepler transit data. We develop and apply a methodology to correct the planetary radius for the presence of both crossed and uncrossed star spots. Star spot crossings are concentrated at phases 0.002 and +0.006. This is consistent with inferences from Rossiter-McLaughlin measurements that the planet transits nearly perpendicular to the stellar equator. We identify the dominant phases of star spot crossings with active latitudes on the star, and infer that the stellar rotational pole is inclined at about 12 deg 5 deg to the plane of the sky. We point out that precise transit measurements over long durations could in principle allow us to construct a stellar Butterfly diagram to probe the cyclic evolution of magnetic activity on this active K-dwarf star.

Deming, Drake↗

Future Visions for Scientific Human Exploration

Today, humans explore deep-space locations such as Mars, asteroids, and beyond, vicariously here on Earth, with noteworthy success. However, to achieve the revolutionary breakthroughs that have punctuated the history of science since the dawn of the Space Age has always required humans as "the discoverers," as Daniel Boorstin contends in this book of the same name. During Apollo 17, human explorers on the lunar surface discovered the "genesis rock," orange glass, and humans in space revamped the optically crippled Hubble Space Telescope to enable some of the greatest astronomical discoveries of all time. Science-driven human exploration is about developing the opportunities for such events, perhaps associated with challenging problems such as whether we can identify life beyond Earth within the universe. At issue, however, is how to safely insert humans and the spaceflight systems required to allow humans to operate as they do best in the hostile environment of deep space. The first issue is minimizing the problems associated with human adaptation to the most challenging aspects of deep space space radiation and microgravity (or non-Earth gravity). One solution path is to develop technologies that allow for minimization of the exposure time of people to deep space, as was accomplished in Apollo. For a mission to the planet Mars, this might entail new technological solutions for in-space propulsion that would make possible time-minimized transfers to and from Mars. The problem of rapid, reliable in-space transportation is challenged by the celestial mechanics of moving in space and the so-called "rocket equation." To travel to Mars from Earth in less than the time fuel-minimizing trajectories allow (i.e., Hohmann transfers) requires an exponential increase in the amount of fuel. Thus, month-long transits would require a mass of fuel as large as the dry mass of the ISS, assuming the existence of continuous acceleration engines. This raises the largest technological stumbling block to moving humans on site as deep-space explorers, delivering the masses required for human spaceflight systems to LEO or other Earth orbital vantage points using the existing or projected fleet of Earth-to-orbit (ETO) launch vehicles. Without a return to Saturn V-class boosters or an alternate path, one cannot imagine emplacing the masses that would be required for any deep-space voyage without a prohibitive number of Shuttle-class launches. One futurist solution might involve mass launch systems that could be used to move the consumables, including fuel, water, food, and building materials, to LEO in pieces rather than launching integrated systems. This approach would necessitate the development of robotic assembly and fuel-storage systems in Earth orbit, but could provide for a natural separation of low-value cargo (e.g., fuel, water).

Garvin, James↗

OSIRIS-REx Encounters Bennu: Initial Assessment from the Approach Phase

The OSIRIS-REx spacecraft launched on September 8, 2016, on a seven-year journey to return samples from asteroid (101955) Bennu. This presentation summarizes the scientific results from the Approach and Preliminary Survey phases. Bennu observations are set to begin on August 17, 2018,when the asteroid is bright enough for detection by the PolyCam. PolyCam and MapCam collect data to survey the asteroid environment for any hazards and characterize the asteroid point-source photometric properties. Resolved images acquired during final approach, starting in late October 2018, allow the creation of a shape model using stereophotoclinometry (SPC), needed by both the navigation team and science planners. The OVIRS and OTES spectrometers characterize the point- source spectral properties over a full rotation period, providing a first look at any features and thermophysical properties. TAGSAM is released from the launch container and deployed into the sampling configuration then returned to the stow position.Preliminary Survey follows the Approach Phase in early December 2018. This phase consists of a series of hyperbolic trajectories that cross over the North and South poles and the equator of Bennu at a close-approach distance of 7 km. Images from these Preliminary Survey passes provide data to complete the 75-cm resolution SPC global shape model and solve for the rotation state. Once the shape model is complete, the asteroid coordinate system is defined for co-registration of all data products. These higher-resolution images also constrain the photometric properties and allow for an initial assessment of the geology. In Preliminary Survey the team also obtains the first OLA data, providing a measure of the surface topography. OVIRS and OTES collect data as "ride-along" instruments, with the spacecraft pointing driven by imaging constraints. These data provide a first look at the spectral variation across the surface of Bennu. Radio science measurements, combined with altimetry and imagery, determine Bennu's mass, a prerequisite to placing the spacecraft into orbit in late December 2018. Together, data from the Approach and Preliminary Survey phases set the stage for the extensive mapping planned for 2019. These dates are the baseline plan. Any contingency or unexpected discovery may change this mission profile.

Lauretta, Dante S.↗

Neutron Stars and Black Holes Seen with the Rossi X-Ray Timing Explorer (RXTE)

Astrophysical X-rays bring information about location, energy, time, and polarization. X-rays from compact objects were seen in the first explorations to vary in time. Eclipses and pulsations have simple explanations that identified the importance of X-ray binaries and magnetic neutron stars in the first decade of X-ray astronomy. The dynamics of accretion onto stellar and supermassive black holes and onto neutron stars with relatively low magnetic fields shows up as more complex variations, quasi-periodic oscillations, noise with characteristic frequency spectra, broad-band changes in the energy spectra. To study these variations, RXTE instruments needed to have large area and operational flexibility to find transient activity and observe when it was present. Proportional counters and Phoswich scintillators provided it in a modest mission that has made textbook level contributions to understanding of compact objects. The first seen, and the brightest known, X-ray binary, Sco X-1 is one of a class of neutron stars with low mass companions. Before RXTE, none of these had been seen to show pulsations, though they were hypothesized to be the precursors of radio pulsars with millisecond periods and low magnetic fields. RXTE's large area led to identifying coherent millisecond pulsars in a subset which are relatively faint transients. It also led to identifying short episodes of pulsation during thermonuclear bursts, in sources where a steady signal is not seen. The X-ray stage verifies the evolution that produces millisecond radio pulsars.Masses and radii of neutron stars are being determined by various techniques, constraining the equation of state of matter at nuclear densities. Accretion should lead to a range of neutron star masses. An early stage of superstrong magnetic field neutron stars is now known to produce X-ray and gamma-ray bursts in crust quakes and magnetic field reconnection releases of energy. Soft Gamma Repeaters, Anomolous X-ray Pulsars, and high magnetic field rotation-powered pulsars are all now called magnetars, because they have pulse periods indicating they are slowing down as they would with magnetic dipole radiation for a surface field above 5x1013 gauss. The accretion disk has been connected to the launching of radio jets from black holes, and even from neutron stars. Estimates of the angular momenta of black holes are being made from different approaches, modelling a high frequency oscillation that may be related to how close the inner part of the accretion disk is to the black hole, modelling the continua spectra of the X-ray emission, and modeling the emission of red-shifted iron that may be emitted from the accretion disk. These investigations require early discovery of the black hole transient with the All Sky Monitor on RXTE or other monitoring information, frequent extended observations, and coordinated observations with missions that give higher energy resolution, or radio and infrared information.

Swank, Jean↗

Discovery of the Neutron Star Spin Frequency in EXO 0748-676

We report the results of a search for burst oscillations during thermonuclear X-ray bursts from the low mass X-ray binary (LMXB) EXO 0748-676. With the proportional counter array (PCA) onboard the Rossi X-ray Timing Explorer (RXTE) we have detected a 45 Hz oscillation in the average power spectrum of 38 thermonuclear X-ray bursts from this source. We computed power spectra with 1 Hz frequency resolution for both the rising and decaying portions of 38 X-ray bursts from the public RXTE archive. We averaged the 1 Hz power spectra and detected a significant signal at 45 Hz in the decaying phases of the bursts. The signal is detected at a significance level of 4 x 10 (exp -8) similar signal was detected in the rising intervals. The oscillation peak is unresolved at 1 Hz frequency resolution, indicating an oscillation quality factor, Q = nu (sub 0)/Delta nu (sub fwhm) greater than 45, and the average signal amplitude is approximately equal to 3% (rms) The detection of 45 Hz burst oscillations from EXO 0748-676 provides compelling evidence that this is the neutron star spin frequency in this system. We use the inferred spin frequency to model the widths of absorption lines from the neutron star surface and show that the widths of the absorption lines from EXO 0748-676 recently reported by Cottam et al. are consistent with a 45 Hz spin frequency as long as the neutron star radius is in the range from about 9.5 - 15 km. With a known spin frequency, precise modelling of the line profiles from EXO 0748-676 holds great promise for constraining the dense matter equation of state.

Villarreal, Adam R.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Autonomous Ocean World Exploration: Advancement of a Virtual Testbed

The search for life (extinct or extant) and potentially habitable bodies in our solar system and beyond is one of the 12 priority science questions outlined in the National Acadamies’ 2022 decadal survey [5]. Extraterrestrial destinations containing liquid water present an opportunity to search for life as we know it, and in recent years an increasing number of such locations have been discovered within our solar system. Several Jovian moons—Europa, Ganymede, and Callisto [10]—and the Saturnian moons Enceladus [8] and Titan [9] are known or suspected to harbor massive subsurface oceans. Of these "ocean worlds", Europa is the focus of at least one planned NASA orbiter mission, Europa Clipper [4], and an early lander mission concept, the Europa Lander [2, 3]. Whereas most robotic missions to the Moon and Mars (e.g. orbiters, rovers, landers) to date have had ground controllers on Earth tightly involved in mission operations, missions to more distant worlds will require a high degree of onboard autonomy due to long communication lags and blackouts, harsh environments (radiation, cold), and more limited battery and hardware life. The past decade has seen great advances in both AI technologies and computing scalability and performance that offer promising solutions for spacecraft autonomy and motivate the software system and research programs described in this paper. The Ocean Worlds Autonomy Testbed for Exploration, Research, and Simulation (OceanWATERS) [1], which has been in development at the NASA Ames Research Center since 2018, is a virtual environment for testing lander autonomy solutions. It is built on the Robot Operating System (ROS), runs on consumer-grade Linux workstations, and was released as open source in 2020. OceanWATERS provides a physical and visual simulation of a prototypical lander in a Europa-like environment (Figure 1). The lander was modeled after requirements and specifications made in JPL’s Europa Lander Study of 2016 [3]. Simulated lander systems include stereo cameras and spotlights mounted on an antenna mast that pans and tilts, a 6 degrees of freedom (DoF) robotic arm with a force-torque sensor and two interchangeable end effectors, and a battery pack power system. The environment consists of multiple terrain models including a highly detailed model sourced from the FROST dataset [11], simulation of surrounding planetary bodies based on an ephemeris model, and lighting from the sun with associated surface illumination, reflectance, and shadows. Operations supported by OceanWATERS include panoramic and directed imaging of the environment and lander workspace, Cartesian and joint-level arm commanding, grinding of the terrain surface (e.g. digging a trench), and scooping of ground material (Figure 2) which can be discarded or collected as science samples in a receptacle that can be emptied (science operations themselves are not simulated). These operations are realized as ROS Actions and are complimented by a wide selection of telemetry that is continually produced by each lander subsystem. The power system model is driven by the open-source Generic Software Architecture for Prognostics (GSAP) [11] that predicts the battery’s remaining useful life and other characteristics. As a testbed for high-level autonomy, OceanWATERS provides an execution framework based on PLEXIL [12], an open-source plan specification language and execution engine developed largely at Ames. NASA's initial development of OceanWATERS, as well the Ocean Worlds Lander Autonomy Testbed (OWLAT) [6], a complimentary physical testbed developed at JPL, was the first step in a plan for realizing candidate onboard autonomy solutions for such planetary landers. In 2020 NASA solicited applications for its Autonomous Robotics Research for Ocean Worlds (ARROW) program, and in 2021 the similar Concepts for Ocean worlds Life Detection Technology (COLDTech) program. Collectively six research teams, based in universities and companies across the United States, were awarded grants to develop and demonstrate autonomy solutions on OceanWATERS and OWLAT. These 1–2-year projects have now finished or are nearing completion, and a wide variety of autonomy challenges in ocean world surface missions were addressed. Prototyped and demonstrated solutions have included autonomous discovery, response and adaptation to system faults and unexpected environmental events, world model synthesis through perception, plan synthesis using learned models, methods to optimize sample target selection and prioritize science data transmission, extension of PLEXIL for stochastic decision-making, and an integration of a model of JPL’s mission-ready COLDArm [7]. Technologies used in these projects include many forms of machine learning, causal reasoning, automated planning, Markov decision processes, formal methods, and other advanced techniques. A more detailed summary of the ARROW and COLDTech projects is given herein. OceanWATERS has had significant enhancements since its open-source release in 2020. Many of its new features were driven or shaped by feedback from the ARROW and COLDTech teams and requirements of their projects. In support of enabling autonomous adaptation to spacecraft faults (a specific capability solicited by both programs), a fault injection and detection framework was developed that supports a wide and growing range of fault types such as locked joints, image loss, and battery failures. The power system model was completed and integrated into the simulator, starting as a single-cell battery model and later upgraded to a multi-cell model with associated faults such as cell disconnection. Arm/terrain interaction was improved by adding a force-torque sensor and associated faults, and an analytic dig force model based on the Balovnev bucket force equations. Environment fidelity was increased by modeling terrain deformation resulting from digging and scooping; visual improvements were made in textures, lighting, and shadows. To facilitate interoperation with OWLAT, a unified command and telemetry interface between the testbeds was developed at the ROS level, along with a PLEXIL interface. The number of lander operations was greatly expanded (e.g. with Cartesian-based arm and antenna movement), and a framework was designed for users to build their own lander actions. A GUI for PLEXIL plan selection was created (Figure 3), and an expansive set of plans were added, such as those that illustrate patterns for fault handling. This paper provides a self-contained high-level description of OceanWATERS, focusing on more detailed coverage of the aforementioned enhancements. It provides a high-level summary of the projects undertaken by participants in the ARROW and COLDTech programs and how these efforts have helped shape OceanWATERS. Finally, potential future work and directions for the testbed are listed, as likely informed by the recent planetary science decadal survey [5].

K Michael Dalal↗

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗

Recent Progress in Understanding the Sun's Magnetic Dynamo

100 years ago we thought that the Sun and stars shone as a result of slow gravitational contraction over a few tens of millions of years - putting astronomers at odds with geologists who claimed that the Earth was much, much older. That mystery was solved in the 1920s and 30s with the discovery of nuclear energy (proving that the geologists had it right all along). Other scientific mysteries concerning the Sun have come and gone but three major mysteries remain: 1) How does the Sun produce sunspots with an 11-year cycle? 2) What produces the huge explosions that result in solar flares, prominence eruptions, and coronal mass ejections? and 3) Why is the Sun's outer atmosphere, the corona, so darned hot? Recent progress in solar astronomy reveals a single key to understanding all three of these mysteries.The 11-year time scale for the sunspot cycle indicates the presence of a magnetic dynamo within the Sun. For decades this dynamo was though to operate within the Sun's convection zone - the outmost 30% of the Sun where convective currents transport heat and advect magnetic lines of force. The two leading theories for the dynamo had very different models for the dynamics of the convection zone. Actual measurements of the dynamics using the techniques of helioseismology showed that both of these models had to be wrong some 20 years ago. A thin layer of strongly sheared flow at the base of the convection zone (now called the tachocline) was then taken to be the seat of the dynamo. Over the last 10 years it has become apparent that a weak meridional circulation within the convection zone also plays a key role in the dynamo. This meridional circulation has plasma rising up from the tachocline in the equatorial regions, spreading out toward the poles at a top speed of about 10-20 m/s at the surface, sinking back down to the tachocline in the polar regions, and then flowing back toward the equator at a top speed of about 1-2 m/s in the tachocline itself. Recent dynamo models that include this meridional flow now appear to have some power for predicting the size of future sunspot cycles.

Hathaway, David. H.↗