Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate dynamic programming”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

A GPS-Based Pitot-Static Calibration Method Using Global Output-Error Optimization

Pressure-based airspeed and altitude measurements for aircraft typically require calibration of the installed system to account for pressure sensing errors such as those due to local flow field effects. In some cases, calibration is used to meet requirements such as those specified in Federal Aviation Regulation Part 25. Several methods are used for in-flight pitot-static calibration including tower fly-by, pacer aircraft, and trailing cone methods. In the 1990 s, the introduction of satellite-based positioning systems to the civilian market enabled new inflight calibration methods based on accurate ground speed measurements provided by Global Positioning Systems (GPS). Use of GPS for airspeed calibration has many advantages such as accuracy, ease of portability (e.g. hand-held) and the flexibility of operating in airspace without the limitations of test range boundaries or ground telemetry support. The current research was motivated by the need for a rapid and statistically accurate method for in-flight calibration of pitot-static systems for remotely piloted, dynamically-scaled research aircraft. Current calibration methods were deemed not practical for this application because of confined test range size and limited flight time available for each sortie. A method was developed that uses high data rate measurements of static and total pressure, and GPSbased ground speed measurements to compute the pressure errors over a range of airspeed. The novel application of this approach is the use of system identification methods that rapidly compute optimal pressure error models with defined confidence intervals in nearreal time. This method has been demonstrated in flight tests and has shown 2- bounds of approximately 0.2 kts with an order of magnitude reduction in test time over other methods. As part of this experiment, a unique database of wind measurements was acquired concurrently with the flight experiments, for the purpose of experimental validation of the optimization method. This paper describes the GPS-based pitot-static calibration method developed for the AirSTAR research test-bed operated as part of the Integrated Resilient Aircraft Controls (IRAC) project in the NASA Aviation Safety Program (AvSP). A description of the method will be provided and results from recent flight tests will be shown to illustrate the performance and advantages of this approach. Discussion of maneuver requirements and data reduction will be included as well as potential applications.

Foster, John V.↗
DOE OSTI2021

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

An investigation of turbulence models

The accuracy to which a turbulent boundary layer or wake can be predicted numerically depends on the validity of the turbulence closure model used. The modeling of turbulence physics is one of the most difficult problems in computational fluid dynamics (CFD). In fact, it is one of the pacing factors in the development of CFD. In general, there are three main approaches to the description of trubulence physics. First is turbulence modeling in which the Reynolds averaged Navier-Stokes equations are used and some closure approximation is made for the Reynolds stresses. A second approach to turbulence is large eddy simulation (LES) in which the computational mesh is taken to be fine enough that the large scale structure of the turbulence can be calculated directly. An empirical assumption must be made for the small scale sub-grid turbulence. The third approach is direct simulation. In this technique the Navier-Strokes equations are solved directly on a mesh which if fine enough to resolve the smallest length scale of the turbulence. The Reynolds averaged equations are not used and no closure assumption is required. These last two approaches require extensive computer resources and as such are not engineering tools. The purpose of the work was to investigate the various engineering turbulence models for accuracy and ease of programming. This involved comparison of the models with each other and with experimental data.

Doria, Michael L.↗
NASA NTRS1997

The Solar-Stellar Connection (NAG5-6124: SOHO Guest Investigator Program)

The following is a final report from the SOHO Guest Investigator program to use the SUMER far-UV spectrometer to obtain imaging spectroscopy in support of the goals of the so-called "solar-stellar connection." In particular, a major emphasis was utilization of the long-slit time-resolved maps of the solar surface in bright far-UV emission lines to deduce how particular aspects of the temporally and spatially averaged line profiles trace back to individual structural features of the magnetically disturbed outer atmosphere; to help interpret the unresolved line profiles from high quality stellar observations (say, with the Space Telescope Imaging Spectrograph). The researchers served two tours of duty in the SOHO Operations Center as SUMER planners, during which time we conducted an extensive series of observing programs. These can be divided into three general categories: surface mapping, translimb spectroscopy, and active region diagnostics. We have analyzed some of the large volumes of data to the point where we have presented them in poster papers, and in invited papers at national and international meetings. Listed below are the titles of the preliminary publications we have written, including brief abstracts to indicate the main results. (1) Chromospheric structure and Dynamics-- Observations -- The chromosphere is a highly structured dynamic 'layer' of the solar outer atmosphere. Here, not only are the effects of mechanical heating first evident (moving upward in altitude from the deep photosphere), but also the amount of nonradiative energy deposited is far greater than in the albeit much hotter overlying transition region and corona. Further, the chromosphere is by far the thickest zone of the solar atmosphere with respect to the pressure scale height. A major goal of stellar astrophysics is to understand how the chromosphere is heated and why it adopts its peculiar structure. A cursory examination of solar filtergrams and high-resolution movies demonstrates that much of the chromospheric "action" must be occurring on fine spatial scales and short times; particularly in the cell interior transient brightenings, but also in the longer-lived network fragments. That regime of investigation is far removed from what one usually associates with "synoptic" measurements. Nevertheless, synoptic observations of chromospheric indices, filtergrams, and globally-averaged profile parameters not only can provide important insight concerning the crucial role of the cycle-variable part of the solar magnetic field; but they also can forge a key link with analogous measurements of the stars, where often the phenomena can be significantly exaggerated from the solar case, but high spatial resolution reconnaissance is not even a remote possibility. In addition to discussing the synoptic aspects of chromospheric structure and dynamics, I summarized new insights into the general problem of the solar chromosphere that have been obtained recently with the SUMER far-ultraviolet spectrometer on SOHO. (2) CO and the Temperature Structure of the Solar Atmosphere -- The surface layers of the Sun provide a crucial boundary condition for many of the processes that occur in the deep interior. The stratification of the outer solar atmosphere once was thought to be well understood. However, studies of thermally sensitive molecular absorptions in the infrared revealed puzzling anomalies. Strong lines of the CO fundamental vibration-rotation bands near 5 microns showed very cool temperatures at the extreme limb, and remarkable off-limb emissions extending well into the supposedly hot chromosphere. The conflicting pictures of the photosphere/chromosphere interface, from the widely separated wavelength regimes, has raised suspicions that those "layers" of the atmosphere are much more inhomogeneous than previously suspected. One proposal is that the low chromosphere is dominated by cool gas, the "COmosphere," which is threaded by a network of persistent small-scale hot magnetic filaments and occasionally disrupted by localized acoustic disturbances. The COmosphere is capped by the merged fields of the network elements in the chromospheric "canopy." (3) Translimb Spectroscopy with SOHO/SUMER -- We have used SUMER to obtain deep exposures of the 1300-1400 A spectrum, at the extreme limb and off-limb. Previous "translimb" studies in the thermal infrared had revealed remarkable extensions of cold material (T is approximately equal to 3000 K)-traced by carbon monoxide emission lines-into the heart of the hot chromosphere. A main objective of our program was to search for corresponding far-UV signatures of the "thermally-bifurcated" low chromosphere; for example, radiatively fluoresced emissions of the CO A-X 4th-positive system (collisional excitation would be negligible in cold gas). We conducted two separate observing programs with SUMER. Both made use of the 1 inch-diameter circular aperture, translated across the limb in the minimum motor step increments of 0."375, along the central meridian in the Northern polar coronal hole. The first program executed for nine hours beginning 19UT 25 Oct 1996. The full wavelength range was 1340-1400 A. It was recorded in two overlapping segments, placing key regions of the spectrum alternately on the KBr and bare parts of the detector, to help isolate 2nd-order features. Each segment was integrated for 500 s, and 32 pairs were obtained to span a 12" swath centered on the optical limb. The second program was conducted 00-09UT 01 Dec 1996. It consisted of a single wavelength setting (1300-1340 A) with exposure time 500 s, but twice the spatial coverage of the earlier series: 64 steps, for a total displacement of 24". The strong chromospheric resonance lines of atomic oxygen (1302-1305 A) and ionized carbon (1334-1335 A) were observed on the bare part of the MCP camera.

Ayres, Thomas R.↗
NASA NTRS2006

Using Apex To Construct CPM-GOMS Models

process for automatically generating computational models of human/computer interactions as well as graphical and textual representations of the models has been built on the conceptual foundation of a method known in the art as CPM-GOMS. This method is so named because it combines (1) the task decomposition of analysis according to an underlying method known in the art as the goals, operators, methods, and selection (GOMS) method with (2) a model of human resource usage at the level of cognitive, perceptual, and motor (CPM) operations. CPM-GOMS models have made accurate predictions about behaviors of skilled computer users in routine tasks, but heretofore, such models have been generated in a tedious, error-prone manual process. In the present process, CPM-GOMS models are generated automatically from a hierarchical task decomposition expressed by use of a computer program, known as Apex, designed previously to be used to model human behavior in complex, dynamic tasks. An inherent capability of Apex for scheduling of resources automates the difficult task of interleaving the cognitive, perceptual, and motor resources that underlie common task operators (e.g., move and click mouse). The user interface of Apex automatically generates Program Evaluation Review Technique (PERT) charts, which enable modelers to visualize the complex parallel behavior represented by a model. Because interleaving and the generation of displays to aid visualization are automated, it is now feasible to construct arbitrarily long sequences of behaviors. The process was tested by using Apex to create a CPM-GOMS model of a relatively simple human/computer-interaction task and comparing the time predictions of the model and measurements of the times taken by human users in performing the various steps of the task. The task was to withdraw $80 in cash from an automated teller machine (ATM). For the test, a Visual Basic mockup of an ATM was created, with a provision for input from (and measurement of the performance of) the user via a mouse. The times predicted by the automatically generated model turned out to approximate the measured times fairly well (see figure). While these results are promising, there is need for further development of the process. Moreover, it will also be necessary to test other, more complex models: The actions required of the user in the ATM task are too sequential to involve substantial parallelism and interleaving and, hence, do not serve as an adequate test of the unique strength of CPM-GOMS models to accommodate parallelism and interleaving.

John, Bonnie↗
NASA NTRS2013

Mixed Integer Programming and Heuristic Scheduling for Space Communication

Optimal planning and scheduling for a communication network was created where the nodes within the network are communicating at the highest possible rates while meeting the mission requirements and operational constraints. The planning and scheduling problem was formulated in the framework of Mixed Integer Programming (MIP) to introduce a special penalty function to convert the MIP problem into a continuous optimization problem, and to solve the constrained optimization problem using heuristic optimization. The communication network consists of space and ground assets with the link dynamics between any two assets varying with respect to time, distance, and telecom configurations. One asset could be communicating with another at very high data rates at one time, and at other times, communication is impossible, as the asset could be inaccessible from the network due to planetary occultation. Based on the network's geometric dynamics and link capabilities, the start time, end time, and link configuration of each view period are selected to maximize the communication efficiency within the network. Mathematical formulations for the constrained mixed integer optimization problem were derived, and efficient analytical and numerical techniques were developed to find the optimal solution. By setting up the problem using MIP, the search space for the optimization problem is reduced significantly, thereby speeding up the solution process. The ratio of the dimension of the traditional method over the proposed formulation is approximately an order N (single) to 2*N (arraying), where N is the number of receiving antennas of a node. By introducing a special penalty function, the MIP problem with non-differentiable cost function and nonlinear constraints can be converted into a continuous variable problem, whose solution is possible.

Lee, Charles H.↗
NASA NTRS1968

Flight Test of a 40-Foot Nominal Diameter Disk-Gap-Band Parachute Deployed at a Mach Number of 2.72 and a Dynamic Pressure of 9.7 Pounds per Square Foot

A 40-foot-nominal-diameter (12.2 meter) disk-gap-band parachute was flight tested as part of the NASA Supersonic Planetary Entry Decelerator (SPED-I) Program. The test parachute was deployed from an instrumented payload by means of a deployment mortar when the payload was at an altitude of 158,500 feet (48.2 kilometers), a Mach number of 2.72, and a free-stream dynamic pressure of 9.7 pounds per foot(exp 2) (465 newtons per meter(exp 2)). Suspension line stretch occurred 0.46 second after mortar firing and the resulting snatch force loading was -8.lg. The maximum acceleration experienced by the payload due to parachute opening was -27.2g at 0.50 second after the snatch force peak for a total elapsed time from mortar firing of 0.96 second. Canopy-shape variations occurred during the higher Mach number portion of the flight test (M greater than 1.4) and the payload was subjected to large amplitude oscillatory loads. A calculated average nominal axial-force coefficient ranged from about 0.25 immediately after the first canopy opening to about 0.50 as the canopy attained a steady inflated shape. One gore of the test parachute was damaged when the deployment bag with mortar lid passed through it from behind approximately 2 seconds after deployment was initiated. Although the canopy damage caused by the deployment bag penetration had no apparent effect on the functional capability of the test parachute, it may have affected parachute performance since the average effective drag coefficient of 0.48 was 9 percent less than that of a previously tested parachute of the same configuration.

EXPERIMENTS↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
NASA NTRS2006

Flight Test Determined Aerodynamics Force and Moment Characteristics of the X-43A Research Vehicle at Mach 7.0

The second flight of the HYPER-X Program afforded a unique opportunity to determine the aerodynamic force and moment characteristics of an airframe integrated scramjet powered aircraft in hypersonic flight. These data were gathered via a repeated series of pitch, yaw, and roll doublets, frequency sweeps, and pull-up/push-over maneuvers performed throughout the X-43A cowl-closed descent phase. The subject flight research maneuvers were conducted in a Mach number range of 6.8 to 0.95 at altitudes from 92,000 ft to sea level. In this flight regime, the dynamic pressure varied from 1300 psf to 400 psf with angle-of-attack ranging from 0 deg to 14 deg. The flight-extracted aerodynamics were compared with pre-flight predictions based on wind tunnel test data. The X-43A flight-derived axial force was found to be 10 to 15 percent higher than prediction. Under-predictions of similar magnitude were observed for the normal force. For Mach numbers greater than 4, the X-43A flight-derived stability and control characteristics resulted in larger than predicted static margins, with the largest discrepancy approximately 5-inches forward along the X(CG) at Mach 6. This would result in less static margin in pitch. The X-43A predicted lateral-directional stability and control characteristics matched well with flight data when allowance was made for the high uncertainty in angle-of-sideslip.

Davis, Mark C.↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code: Preprint

A passive yaw implementation is developed, validated, and explored for the Wave Energy Converter Simulator (WEC-Sim), an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01 was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than or equal to 1, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
NASA NTRS1999

PARC Analysis of the NASA/GE 2D NRA Mixer/Ejector Nozzle

Interest in developing a new generation supersonic transport has increased in the past several years. Current projections indicate this aircraft would cruise at approximately Mach 2.4, have a range of 5000 nautical miles and carry at least 250 passengers. A large market for such an aircraft will exist in the next century due to a predicted doubling of the demand for long range air transportation by the end of the century and the growing influence of the Pacific Rim nations. Such a proposed aircraft could more than halve the flying time from Los Angeles to Tokyo. However, before a new economically feasible supersonic transport can be built, many key technologies must be developed. Among these technologies is noise suppression. Propulsion systems for a supersonic transport using current technology would exceed acceptable noise levels. All new aircraft must satisfy FAR 36 Stage III noise regulations. The largest area of concern is the noise generated during takeoff. A concerted effort under NASA's High Speed Research (HSR) program has begun to address the problem of noise suppression. One of the most promising concepts being studied in the area of noise suppression is the mixer/ejector nozzle. This study analyzes a typical noise suppressing mixer ejector nozzle at take off conditions, using a Full Navier-Stokes (FNS) computational fluid dynamics (CFD) code.

DeBonis, J. R.↗
NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasizes the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M qi) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in turns of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, Derck↗
NASA NTRS2002

The Winds of B Supergiants

We present the most suitable data sets available in the International Ultraviolet Explorer (IUE) archive for the study of time-dependent stellar winds in early B supergiants. The UV line profile variability in 11 B0 to B3 stars is analyzed, compared and discussed, based on 16 separate data sets comprising over 600 homogeneously reduced high-resolution spectrograms. The targets include 'normal' stars with moderate rotation rates and examples of rapid rotators. A gallery of grey-scale images (dynamic spectra) is presented, which demonstrates the richness and range of wind variability and highlights different structures in the winds of these stars. This work emphasises the suitability of B supergiants for wind studies, under-pinned by the fact that they exhibit unsaturated wind lines for a wide range of ionization. The wind activity of B supergiants is substantial and has highly varied characteristics. The variability evident in individual stars is classified and described in terms of discrete absorption components, spontaneous absorption, bowed structures, recurrence, and ionization variability and stratification. Similar structures can occur in stars of different fundamental parameters, but also different structures may occur in the same star at a given epoch. We discuss the physical phenomena that may be associated with the spectral signatures, and highlight the challenges that these phenomena present to theoretical studies of time-dependent outflows in massive stars. In addition, SEI line-synthesis modelling of the UV wind lines is used to provide further information about the state of the winds in our program stars. Typically the range, implied by the line profile variability, in the product of mass-loss rate and ion fraction (M (dot) q(sub i)) is a factor of approximately 1.5, when integrated between 0.2 and 0.9 v infinity; it can however be several times larger over localized velocity regions. At a given effective temperature the mean relative ion ratios can differ by a factor of 5. The general excess in predicted (forward-scattered) emission in the low velocity regime is discussed in terms of structured outflows. Mean ion fractions are estimated over the B0 to B1 spectral classes, and trends in the ionic ratios as a function of wind velocity are described. The low values obtained for the ion fractions of UV resonance lines may reflect the role of clumping in the wind.

Massa, D.↗
NASA NTRS1986

Dynamic posture analysis of Spacelab-1 crew members

Dynamic posture testing was conducted on the science crew of the Spacelab-1 mission on a single axis linear motion platform. Tests took place in pre- and post-flight sessions lasting approximately 20 min each. The pre-flight tests were widely spaced over the several months prior to the mission while the post-flight tests were conducted over the first, second, fourth, and sixth days after landing. Two of the crew members were also tested on the day of landing. Consistent with previous postural testing conducted on flight crews, these crew members were able to complete simple postural tasks to an acceptable level even in the first few hours after landing. Our tests were designed to induce dynamic postural responses using a variety of stimuli and from these responses, evaluate subtle changes in the postural control system which had occurred over the duration of the flight. Periodic sampling post-flight allowed us to observe the time course of readaptation to terrestrial life. Our observations of hip and shoulder position, when subjected to careful analysis, indicated modification of the postural response from pre- to post-flight and that demonstrable adjustments in the dynamic control of their postural systems were taking place in the first few days after flight. For transient stimuli where the platform on which they were asked to stand quickly moved a few centimeters fore or aft then stopped, ballistic or open loop 'programs' would closely characterize the response. During these responses the desired target position was not always achieved and of equal importance not always properly corrected some 15 seconds after the platform ceased to move. The persistent observation was that the subjects had a much stronger dependence on visual stabilization post-flight than pre-flight. This was best illustrated by a slow or only partial recovery to an upward posture after a transient base-of-support movement with eyes open. Postural responses to persistent wideband pseudorandom base-of-support translation were modeled as time invarient linear systems arrived at by Kalman adaptive filter techniques. Derived model parameters such as damping factor and fundamental frequency of the closed loop system showed significant modification between pre- and post-flight. This phenomenon is best characterized by movement of the poles toward increasing stability. While pre-flight data tended to show shoulders and hips moving in phase with each other, post-flight data showed a more disjoint behavior.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
NASA NTRS1991

Evolutionary dynamics of plants and animals: a comparative approach

Patterns of longevity and rate of appearance of taxa in the fossil record indicate a different evolutionary dynamic between land plants and marine invertebrates. Among marine invertebrates, rates of taxonomic turnover declined through the Phanerozoic, with increasingly extinction-resistant, long-lived, clades coming to dominate. Among terrestrial vascular plants, rates of turnover increased through the Phanerozoic, with short-lived, extinction-prone clades coming to dominate from the Devonian to the present. Terrestrial vertebrates appear to approximate the marine invertebrate pattern more closely than the plant record. We identify two features which individually or jointly may have influenced this distinction. First, land plants continuously invaded stressful environments during their evolution, while marine invertebrates and terrestrial vertebrates did not. Second, the relative structural simplicity and indeterminate mode of plant growth vs. the relative structural complexity and determinate mode of animal growth may have influenced the timing of major clade origin in the two groups.

Review↗
NASA NTRS2001

CFD Simulations in Support of Shuttle Orbiter Contingency Abort Aerodynamic Database Enhancement

Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20deg-60deg, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). Presented here are details of the methodology and comparisons of computed aerodynamic coefficients against the values in the current Orbiter Operational Aerodynamic Data Book (OADB). While approximately 40 cases have been computed, only a sampling of the results is provided here. The computed results, in general, are in good agreement with the OADB data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects. The aerodynamic coefficients and detailed surface pressure distributions of the present simulations are being used by the Shuttle Program in the evaluation of the capabilities of the Orbiter in contingency abort scenarios.

Papadopoulos, Periklis E.↗
DOE OSTI2024

Semiglobal Safety-Filtered Extremum Seeking With Unknown CBFs

We introduce a safe extremum-seeking (Safe ES) algorithm which achieves the minimization of an unknown objective function while ensuring that an unknown, yet measured, control barrier function (CBF) remains above an arbitrarily small negative value for all time. In other words, “practical safety” is maintained during the entire period of convergence to the constrained extremum. Our design is based on quadratic program (QP) CBF style filters for safety, which is applied in an average and estimated sense. Using nonsmooth analysis tools, we guarantee semiglobal practical asymptotic (SPA) stability of the global constrained optimum, practical convergence to the safe set if starting in a condition violating the CBF, and practical safety for all time—semiglobally—if starting in safe set. The safety result of the paper is analogous with modern notions of SPA stability, guaranteeing that, for any small violation of safety, there exist design coefficients which guarantee that such a small violation is not exceeded. The paper outlines a set of sufficient conditions on the barrier and objective functions, and by way of a Lyapunov argument, we demonstrate that nonconvex constrained optimization problems can be solved. We present these results in the setting of a static map and a dynamical system. A simulation example illustrates the results.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

New York City Panel on Climate Change 2019 Report Chapter 9: Perspectives on a City in a Changing Climate 2008-2018

Cities experience multiple environmental shifts, stresses, and shocks—such as air and water pollution—and a variety of extreme events simultaneously and continuously. Current urban programs have focused on limiting the impacts of these conditions through a portfolio of multifaceted strategies, such as regulations and codes, management and restoration projects, and citizen engagement. Global climate change represents a new environmental dynamic to which cities now have to respond. While global climate change by definition has impacts world wide, residents and managers of cities, like New York, typically perceive changes in their own local environments. In most cities, temperature is warming with increasingly hotter and longer heatwaves, and heavier downpours are leading to more frequent inland flooding. In coastal cities, sea levels are rising, exacerbating coastal flooding. Analyzing and understanding the impacts of climate change on cities is important because of the dramatic growth in urban populations throughout the world. An estimated nearly 4.0 billion people reside in urban areas, accounting for 52% of the world’s population (UN, 2017). That percentage will increase dramatically in the coming decades as almost all of the growth to take place up to 2050 will be in urban areas (UN, 2017). The New York City metropolitan region (NYMR)—the five boroughs (equivalent to counties) of New York City and the adjacent 26 counties in the states of New York, New Jersey, and Connecticut—is an ideal model of an urban agglomeration. Approximately 8.6 million people live in the five boroughs and more than 15 million people live in the neighboring smaller cities, towns, and villages (City of New York, 2018a; US Census, 2017). The population of the five boroughs is projected to add 1 million people by 2030, while the total region is projected to reach 26.1 million (NYTC, 2015).

Solecki, William↗