Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Unobtrusive Software and System Health Management with R2U2 on a Parallel MIMD Coprocessor

Dynamic monitoring of software and system health of a complex cyber-physical system requires observers that continuously monitor variables of the embedded software in order to detect anomalies and reason about root causes. There exists a variety of techniques for code instrumentation, but instrumentation might change runtime behavior and could require costly software re-certification. In this paper, we present R2U2E, a novel realization of our real-time, Realizable, Responsive, and Unobtrusive Unit (R2U2). The R2U2E observers are executed in parallel on a dedicated 16-core EPIPHANY co-processor, thereby avoiding additional computational overhead to the system under observation. A DMA-based shared memory access architecture allows R2U2E to operate without any code instrumentation or program interference.

Schumann, Johann↗
NASA NTRS2018

The NANOGrav 11-year Data Set: High-Precision Timing of 45 Millisecond Pulsars

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, and six high-timing-precision pulsars were observed weekly. All were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and (for binary pulsars) orbital parameters; time-variable dispersion delays; and parameters that quantify pulse-profile evolution with frequency. The timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of significant orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. A companion paper uses these data to constrain the strength of the gravitational-wave background

Arzoumanian, Zaven↗
NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗
NASA NTRSDate not supplied

NSCOR for Evaluating Risk Factors and Biomarkers for Adaptation and Resilience to Spaceflight: Emotional Valence and Social Processes in ICC/ICE Environments

Space exploration class missions, such as a mission to Mars, will require optimization of human performance, adaptability, and resilience. This NASA Specialized Center of Research (NSCOR) utilizes the NIMH Research Domain Criteria (RDoC) framework to identify biological and behavioral markers of individual social adaptation and emotional resilience (as well as vulnerability) to spaceflight-relevant stressors such as living in extended isolation. The overarching goal of this NSCOR is to obtain novel information to help identify biomarkers of individuals who are resilient and/or adaptable to the stressors of isolated, confined, and controlled (ICC) and isolated, confined, and extreme (ICE) environments.A total of N=90 healthy adult astronaut surrogates are being studied in three spaceflight-analog environments: (1) n=40 healthy adults in the Isolation and Confinement Analog Research Unit (ICARUS), an ICC at the University of Pennsylvania, during 7-day missions, for a target total of 280 subject days; (2) n=32 healthy adult astronaut surrogates studied in NASA’s Human Exploration Research Analog (HERA), an ICC at Johnson Space Center during 45-day missions, for a target total of 2,112 subject days; and (3) n=18 healthy adults in the Alfred-Wegener-Institute’s Neumayer Station III, an ICE in Antarctica, during 14-month missions, for a target total of 7,560 subject days. Dr. Nindl’s Laboratory at the University of Pittsburgh is analyzing a priori selected protein biomarkers in blood, saliva, and urine. Complementary rodent models of exposure to early life stressors, confinement, and isolation are being evaluated at Dr. Hensch’s Laboratory to further validate the neurobehavioral and biological findings from the human studies.Given the inconsistency and varied definition of resilience/adaptation in the scientific literature, the NSCOR team developed a composite resilience/adaptation measure that reflects the most relevant outcomes to resilience/adaptation across psychosocial and neurobehavioral functions, as well as neurocognitive and spaceflight-relevant operational performance. To achieve this, group consensus was attained from subject matter experts to produce a rank-order of importance for each input variable. The final resilience/adaptation score included 36 variables that were collected across spaceflight analogs. As of 10/1/2021, the NSCOR project acquired data on n=27 subjects at ICARUS, n=16 at HERA, and n=18 at Neumayer. The COVID-19 pandemic delayed data acquisition at ICARUS and HERA.Among subjects studied to date, 99% of neural and neurobehavioral data (e.g., neuroimaging for structure and function, behavioral measures) as well as blood, saliva, and urine for biochemical assays have been acquired. For rodent models, Dr. Hensch’s laboratory has established biochemical and behavioral parameters reflecting confinement stress in social networks of mice for comparison to stress responses in the human spaceflight analog environments.Group social behaviors were measured with a Social Network Analysis (SNA) approach to define objective parameters associated with sociability and its plasticity by sex. This analytic approach may help identify a network of individuals who are more effective teammates or more likely to generate new social relationships. Data acquisition, biomarker assessment, and data quality control will continue through September 2022.

D F Dinges↗
NASA NTRS2021

Romie: A Domain-Independent Tool for Computer-Aided Robust Operations Management

Romie is a decision support tool based on AI's latest advances in the domain of robust scheduling. Unlike all its predecessors, the tool allows to (i) visually model the operational problem and context entirely (ii) optimize to find near-optimal schedules while taking uncertainty into account and deals with (iii) a combination of various {key performance indicators (KPIs). It comes with a web user interface. Part or all of the modelled activities may be associated to random variables describing their stochastic durations, in order to produce schedules that are robust w.r.t. temporal uncertainty. Hence, depending on the pursued KPIs, the schedules maximize a combination of the following terms: the probability of satisfying the problem constraints, the expected return/efficiency, the expected outcome quality, and even the operators' wellness by minimizing its expected extra-hours. Initially developed for spatial exploration and demonstration in the context of Mars analog missions, this versatile tool is here applied to operations management in both biotechnology manufacturing and robots parametrization in a cave exploration context.

Chien, Steve A.↗
DOE OSTI2019

MSE based drilling optimization using neural network simulaton

The method disclosed receives a data stream from an MWD system and determines the response of a specific energy (SE) relationship and a rate of penetration (ROP) relationship respectively to variables controllable by the operator, in order to enable operation at a lowest SE, or a highest Rate-of-Penetration (ROP) to SE ratio. The method utilizes artificial neural networks trained by MWD data to deduce a depth-of-cut and torque based on relationships manifesting between the various data points collected, and an SE equation and a predicted ROP is evaluated over a series of probable operating points. The method continuously gathers and analyzes MWD data during the drilling operation and allows an operator to manage the controllable parameters such that operation at the lowest SE or highest ROP or ROP to SE ratio can be achieved during the drilling operation.

Zhang, Wu↗
NASA NTRS1988

Adaptive servo control for umbilical mating

Robotic applications at Kennedy Space Center are unique and in many cases require the fime positioning of heavy loads in dynamic environments. Performing such operations is beyond the capabilities of an off-the-shelf industrial robot. Therefore Robotics Applications Development Laboratory at Kennedy Space Center has put together an integrated system that coordinates state of the art robotic system providing an excellent easy to use testbed for NASA sensor integration experiments. This paper reviews the ways of improving the dynamic response of the robot operating under force feedback with varying dynamic internal perturbations in order to provide continuous stable operations under variable load conditions. The goal is to improve the stability of the system with force feedback using the adaptive control feature of existing system over a wide range of random motions. The effect of load variations on the dynamics and the transfer function (order or values of the parameters) of the system has been investigated, more accurate models of the system have been determined and analyzed.

Zia, Omar↗
NASA NTRS1988

Spatial and temporal variability of VAS radiance measurements by structure and correlation analysis

The statistical structure function analysis presently applied to VISSR Atmospheric Sounder (VAS) measurements has been extended to include time, and yields structure plots in either two spatial dimensions or one spatial dimension and time that indicate three-dimensional measurement variability. The analyses that include time as a coordinate also yield an indication of the mean speed and direction of the analyzed data. Results for three-hourly VAS data indicate that sampling at a high, approximately 1-hour, frequency is required in order to correctly monitor VAS measurements' temporal variability in a way equivalent to high spatial resolution.

Hillger, Donald W.↗
DOE OSTI2019

A parallel-in-time algorithm for variable step multistep methods

This paper presents a multigrid reduction in time (MGRIT) algorithm for achieving time parallelism using multistep backward difference formula (BDF) methods on variably-spaced temporal grids. This MGRIT approach transforms the linear multistep methods into single step methods applied to groups of time steps. Stability considerations are addressed through lowering of the order on coarse grids. The methods are presented for fixed and variable time step formulations. Furthermore, numerical results show moderate speedups for the heat equation as well as two IEEE power grid test problems characterized by nonlinear differential algebraic equations (DAE) of index 1. The performance of MGRIT with BDF is compared to MGRIT with Runge-Kutta for both the fixed and variable step methods of the same order.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗
NASA NTRS1988

Viscous flow computations using a second-order upwind differencing scheme

In the present computations of a wide range of fluid flow problems by means of the primitive variables-incorporating Navier-Stokes equations, a mixed second-order upwinding scheme approximates the convective terms of the transport equations and the scheme's accuracy is verified for convection-dominated high Re number flow problems. An adaptive dissipation scheme is used as a monotonic supersonic shock flow capture mechanism. Many benchmark fluid flow problems, including the compressible and incompressible, laminar and turbulent, over a wide range of M and Re numbers, are presently studied to verify the accuracy and robustness of this numerical method.

Chen, Y. S.↗
NASA NTRS1991

Aerodynamic shape optimization using sensitivity analysis on third-order Euler equations

Two major advancements of an aerodynamic optimization method with two design variables are described. The first is represnted by an improved flow prediction method that uses the third-order accurate discretization of the Euler equations. This method makes it possible to predict the flowfield of a perturbed shape which generates shocks and other large gradients without intermediate CFD analysis. The second advancement is the use of every surface grid point as a design variable in the aerodynamic shape optimization problem. The improved algorithm is demonstrated by optimizing the ramp shape of a scramjet-afterbody configuration for maximum axial thrust. It is concluded that the improved flow-field prediction method eliminates most of the flow analysis and a priori guessing of all possible shapes from which the optimum is to be selected.

Baysal, Oktay↗
DOE OSTI2017

Arrhenius Rate: constant volume burn

A constant volume burn occurs for an idealized initial state in which a large volume of reactants at rest is suddenly raised to a high temperature and begins to burn. Due to the uniform spatial state, there is no fluid motion and no heat conduction. This reduces the time evolu tion to an ODE for the reaction progress variable. With an Arrhenius reaction rate, two characteristics of thermal ignition are illustrated: induction time and thermal runaway. The Frank-Kamenetskii approximation then leads to a simple expression for the adiabatic induction time. For a first order reaction, the analytic solution is derived and used to illustrate the effect of varying the activation temperature; in particular, on the induction time. In general, the ODE can be solved numerically. This is used to illustrate the effect of varying the reaction order. We note that for a first order reaction, the time evolution of the reaction progress variable has an exponential tail. In contrast, for a reaction order less than one, the reaction completes in a nite time. The reaction order also affects the induction time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Stochastic Models of Tropical Rain-Rate Fields

Because of the extreme variability of rain rate in space and time and the difficulties with remote sensing methods of measuring rain rates, accurate determination of rainfall over large areas and time periods has long been a problem for hydrologists, meteorogists, and climatologists. A number of statistical models of rain have been developed in order to investigate the impact of rain variability on satellite remote sensing methods, validation of satellite rain products, and generation of rain maps with accompanying error estimates. These models may be useful in examining 'sub-grid scale' issues in representing precipitation in numerical mdoels. A stochastic model will first be described which can generate time-dependent high-resolution spatial rain fields with space and time correlations similar to those seen in rain data, as well as representing the presence of areas with zero rain rate and log-normally distributed rain rates where there is rain. A simpler model derived from this, formulated in the spectral domain, seems to imply fractal-like rain statistics at small scales when fit to rain data.

Bell, Thomas L.↗
NASA NTRS1999

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗
DOE OSTI2024

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Gradients of galactic cosmic rays and anomalous components

Measurements of radial and latitudinal gradients of galactic cosmic rays and anomalous components now cover radii from 0.3 to 40 AU from the sun and latitudes up to 30 deg above the ecliptic plane for particle energies from approx. 10 MeV/n up to relativistic energies. The most accurate measurements cover the period 1972 through 1987, which includes more than one full 11 year cycle of solar activity. Radial gradients for glactic cosmic rays of all energies and species are small (similar to less than 10 percent AU), and variable in time, reaching a minimum of near 0 percent AU out to 30 AU for some species at solar maximum. Gradients for anomalous components are larger, of order 15 percent AU, may show similar time variability, and are relatively independent of particle species and energy. For the period 1985 through 1986 the intensity decreased away from the ecliptic for all species and energies. For galactic cosmic rays, the measured gradients are approx. 0.5 percent/degree near 20 AU, while for anomalous components the gradients are larger, ranging from 3 to 6 percent/degree. Comparison with a similar measurement for anomalous helium in 1975 through 1976 suggests that the latitude gradients for anomalous components have changed sign between 1975 and 1985. For galactic cosmic rays, the available evidence suggests no change in sign of the latitudinal gradient for relativistic particles.

Mckibben, R. B.↗
NASA NTRS1978

Interaction of a strong blast wave with a free surface

When a point source explosion is initiated at the ocean surface, the shock propagated into the water is reflected at the surface as a centered expansion wave. The solution in the neighborhood of the interaction point is obtained by writing the equations of motion in the appropriate similarity variables and then changing the independent variables to polar coordinates based at the interaction point. From the zero-order solution of the resulting equations the slopes of boundaries at the interaction point are obtained. A first-order perturbation of this solution provides more accurate representation of the flow variables and the curvature of the shock surface near the interaction point.

Falade, A.↗