Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Programming, Linear”

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 1,423 records · Page 79

A study of the sonic-boom characteristics of a blunt body at a Mach number of 6

An experimental and analytical study of the sonic boom static pressure signatures generated by a blunt body at mach 6 has shown that finite difference computer programs can be used to give reasonable estimates of the signatures. The calculated near field static pressure signature was extrapolated to the far field by a program using the method of characteristics. A comparison of this extrapolated signature with the signature predicted by far field sonic boom theory (linearized) shows that peak overpressures are about the same, at least up to mach, but the far field theory overestimates the length of the signature.

Ashby, G. C., Jr.↗

Wing Subsonic Aerodynamic Performance Estimates

Solution by iteration estimates performance of twisted and cambered wings of arbitrary planform. SUBAER based on linearized theory lifting-surface solution provides spanwise distribution of theoretical leading edge thrust, in addition to surface distribution of perturbation velocities. Program written in FORTRAN IV.

Carlson, H. W.↗

A three degree-of-freedom, typical section flutter analysis using harmonic transonic air forces

Typical section flutter solutions are obtained using generalized air forces calculated with OPTRAN2, a program based on the transonic small perturbation potential equation, which uses potentials from a nonlinear solution to the steady flow to obtain linear unsteady air forces. Flutter results are calculated with a conventional, linear V-g type solution. Results are presented for a NACA 64A010 airfoil oscillating in plunge and pitch, and with a quarter-chord control surface in the Mach number range of 0.80 to 0.90. The flutter boundary shows a significant 'transonic' bucket, together with changes in flutter mode as, with increasing Mach numbers, a shock appears on the airfoil and then moves aft. Single degree-of-freedom instabilities associated with both the pitch and the control surface motions are identified.

Weatherill, W. H.↗

A new method for determining acoustic-liner admittance in a rectangular duct with grazing flow from experimental data

A method is developed for determining acoustic liner admittance in a rectangular duct with grazing flow. The axial propagation constant, cross mode order, and mean flow profile is measured. These measured data are then input into an analytical program which determines the unknown admittance value. The analytical program is based upon a finite element discretization of the acoustic field and a reposing of the unknown admittance value as a linear eigenvalue problem on the admittance value. Gaussian elimination is employed to solve this eigenvalue problem. The method used is extendable to grazing flows with boundary layers in both transverse directions of an impedance tube (or duct). Predicted admittance values are compared both with exact values that can be obtained for uniform mean flow profiles and with those from a Runge Kutta integration technique for cases involving a one dimensional boundary layer.

Watson, W. R.↗

Automated problem scheduling and reduction of synchronization delay effects

It is anticipated that in order to make effective use of many future high performance architectures, programs will have to exhibit at least a medium grained parallelism. A framework is presented for partitioning very sparse triangular systems of linear equations that is designed to produce favorable preformance results in a wide variety of parallel architectures. Efficient methods for solving these systems are of interest because: (1) they provide a useful model problem for use in exploring heuristics for the aggregation, mapping and scheduling of relatively fine grained computations whose data dependencies are specified by directed acrylic graphs, and (2) because such efficient methods can find direct application in the development of parallel algorithms for scientific computation. Simple expressions are derived that describe how to schedule computational work with varying degrees of granularity. The Encore Multimax was used as a hardware simulator to investigate the performance effects of using the partitioning techniques presented in shared memory architectures with varying relative synchronization costs.

Saltz, Joel H.↗

An application of high authority/low authority control and positivity

Control Dynamics Company (CDy), in conjunction with NASA Marshall Space Flight Center (MSFC), has supported the U.S. Air Force Wright Aeronautical Laboratory (AFWAL) in conducting an investigation of the implementation of several DOD controls techniques. These techniques are to provide vibration suppression and precise attitude control for flexible space structures. AFWAL issued a contract to Control Dynamics to perform this work under the Active Control Technique Evaluation for Spacecraft (ACES) Program. The High Authority Control/Low Authority Control (HAC/LAC) and Positivity controls techniques, which were cultivated under the DARPA Active Control of Space Structures (ACOSS) Program, were applied to a structural model of the NASA/MSFC Ground Test Facility ACES configuration. The control systems design were accomplished and linear post-analyses of the closed-loop systems are provided. The control system designs take into account effects of sampling and delay in the control computer. Nonlinear simulation runs were used to verify the control system designs and implementations in the facility control computers. Finally, test results are given to verify operations of the control systems in the test facility.

Seltzer, S. M.↗

Probabilistic methods for rotordynamics analysis

This paper summarizes the development of the methods and a computer program to compute the probability of instability of dynamic systems that can be represented by a system of second-order ordinary linear differential equations. Two instability criteria based upon the eigenvalues or Routh-Hurwitz test functions are investigated. Computational methods based on a fast probability integration concept and an efficient adaptive importance sampling method are proposed to perform efficient probabilistic analysis. A numerical example is provided to demonstrate the methods.

Wu, Y.-T.↗

Software For Generating Graphs And Charts

Common Graphics Library (CGL) computer program designed to enable users to generate graphs and charts of quality sufficient for publication or optical projection. Quickly and easily generates linear, logarithmic, bar, pie, and composite charts. Features of linear charts include automatic scaling, increasing or decreasing numerical axes, and character axes. Features of pie charts include segment labels, exploded segments, and chart keys. One of general features, horizontal or vertical orientation of figures. Written in ANSI FORTRAN 77.

Hammond, Dana P.↗

Benchmarking Gas Path Diagnostic Methods: A Public Approach

Recent technology reviews have identified the need for objective assessments of engine health management (EHM) technology. The need is two-fold: technology developers require relevant data and problems to design and validate new algorithms and techniques while engine system integrators and operators need practical tools to direct development and then evaluate the effectiveness of proposed solutions. This paper presents a publicly available gas path diagnostic benchmark problem that has been developed by the Propulsion and Power Systems Panel of The Technical Cooperation Program (TTCP) to help address these needs. The problem is coded in MATLAB (The MathWorks, Inc.) and coupled with a non-linear turbofan engine simulation to produce "snap-shot" measurements, with relevant noise levels, as if collected from a fleet of engines over their lifetime of use. Each engine within the fleet will experience unique operating and deterioration profiles, and may encounter randomly occurring relevant gas path faults including sensor, actuator and component faults. The challenge to the EHM community is to develop gas path diagnostic algorithms to reliably perform fault detection and isolation. An example solution to the benchmark problem is provided along with associated evaluation metrics. A plan is presented to disseminate this benchmark problem to the engine health management technical community and invite technology solutions.

Simon, Donald L.↗

Mission-Maps For Outbound Cislunar Transfer Trajectories

This study quantifies the robustness and sensitivity of an outbound cislunar trajectory for a lunar lander in the form of mission-maps, or topological maps that allows either a computer program or mission designer to intuitively optimize the placement of critical outbound correction burns from the derived sensitivity data. The non-linear multi-body dynamics are applied to generate an outbound cislunar reference profile used by a linear covariance analysis (LinCov) tool to compute the expected Δv and trajectory dispersions due to the initial state uncertainty, sensor errors, maneuver execution errors, and disturbance accelerations along the outbound cislunar profile. The rapid performance analysis capabilities of LinCov are complimented with parallel processing techniques to evaluate hundreds and thousands of different translational burn locations, placements, and targeting constraints to identify the combination that minimizes the total Δv usage (nominal plus 3σ Δv) and trajectory dispersions at lunar orbit insertion. This study utilizes a generalized reference targeting algorithm to quickly assess the integrated closed-loop GN&C system performance due to different targeting configurations and constraints. The resulting mission maps provide an intuitive insight to ascertain each trajectory correction maneuver’s (TCM) sensitivity to different burn times along an outbound cislunar trajectory and quickly identify desirable engineering tradeoffs when performing analysis on the number and placement of these burns that nominally zero. Multiple mission maps are generated for a variety of different performance parameters that allow engineers to visually identify optimal solutions for trajectory correction maneuver placements, the number of correction burns, and the targeting constraints for each burn.

GN&C↗

A data-driven linear formulation of the optimal demand response scheduling problem for an industrial air separation unit

Demand response (DR) has become a key element in balancing the power grid as the contribution of time-varying renewable power generation increases. Chemical plants are appealing candidates for DR programs as they offer large, concentrated and flexible loads. DR participation calls for frequent production rate changes over time scales that overlap with the dominant dynamics of the plant. Production scheduling should therefore consider the process dynamics explicitly. Here we present a data-driven approach for modelling the scheduling-relevant dynamics based on historical closed-loop operating data using autoregressive with extra inputs (ARX) models. We introduce a new, linear scheduling problem formulation based on the ARX representation, and demonstrate its implementation on an industrial air separation unit.

42 ENGINEERING↗

Integrated Chemical and Microorganism Monitoring of Air Using Gas Chromatography/Ion Mobility Spectometry: Toward an Expanded-Use Volatile Organic Analyzer (VOA)

The work described in this research program originated with the choice by NASA of an ion mobility spectrometer for air quality monitoring on-board the international spacestation. Though the gas chromatograph-ion mobility spectrometer analyzer known as VOA met or exceeded expectations, limitations in the basic understanding of response and the utilization of foundational principles into usable technology was considered unacceptable. In this research program, a comprehensive model for the origins of mobility spectra was proposed, tested and verified. The principles considered responsible for the appearance of mobility spectra have now been elucidated through this project. This understanding has been applied in automated identification of mobility spectra using neural networks and routine procedures for this now exist. Finally, the limitation on linear range has been shown to be a technical limitation and not a fundamental limitation so that a hardware component was crafted to extend the linear range of a mobility spectrometer by 10X. This project has led to one Ph.D. dissertation and one MS thesis. In addition, over ten public presentations at professional meetings and six journal publications have resulted from this program of research. The findings are so plentiful that total analysis of the findings may require four to six years or more. The findings confirm that the decision to use VOA was sound and that the chemical and physical principles of mobility spectrometry are both understandable and predictable.

Eiceman, G. A.↗

Comparative kinetics and reciprocal inhibition of nitrate and nitrite uptake in roots of uninduced and induced barley (Hordeum vulgare L.) seedlings

Nitrate and NO2- transport by roots of 8-day-old uninduced and induced intact barley (Hordeum vulgare L. var CM 72) seedlings were compared to kinetic patterns, reciprocal inhibition of the transport systems, and the effect of the inhibitor, p-hydroxymercuribenzoate. Net uptake of NO3- and NO2- was measured by following the depletion of the ions from the uptake solutions. The roots of uninduced seedlings possessed a low concentration, saturable, low Km, possibly a constitutive uptake system, and a linear system for both NO3- and NO2-. The low Km system followed Michaelis-Menten kinetics and approached saturation between 40 and 100 micromolar, whereas the linear system was detected between 100 and 500 micromolar. In roots of induced seedlings, rates for both NO3- and NO2- uptake followed Michaelis-Menten kinetics and approached saturation at about 200 micromolar. In induced roots, two kinetically identifiable transport systems were resolved for each anion. At the lower substrate concentrations, less than 10 micromolar, the apparent low Kms of NO3- and NO2- uptake were 7 and 9 micromolar, respectively, and were similar to those of the low Km system in uninduced roots. At substrate concentrations between 10 and 200 micromolar, the apparent high Km values of NO3- uptake ranged from 34 to 36 micromolar and of NO2- uptake ranged from 41 to 49 micromolar. A linear system was also found in induced seedlings at concentrations above 500 micromolar. Double reciprocal plots indicated that NO3- and NO2- inhibited the uptake of each other competitively in both uninduced and induced seedlings; however, Ki values showed that NO3- was a more effective inhibitor than NO2-. Nitrate and NO2- transport by both the low and high Km systems were greatly inhibited by p-hydroxymercuribenzoate, whereas the linear system was only slightly inhibited.

NASA Program CELSS↗

Do nonlinearities play a significant role in short term, beat-to-beat variability?

Numerous studies of short-term beat-to-beat variability in cardiovascular signals have not resolved the debate about the completeness of linear analysis techniques. This aim of this paper is to evaluate further the role of nonlinearities in short-term, beat-to-beat variability. We compared linear autoregressive moving average (ARMA) and nonlinear neural network (NN) models for predicting instantaneous heart rate (HR) and mean arterial blood pressure (BP) from past HR and BP. To evaluate these models, we used HR and BP time series from the MIMIC database. Experimental results indicate that NN-based nonlinearities do not play a significant role and suggest that ARMA linear analysis techniques provide adequate characterization of the system dynamics responsible for generating short-term, beat-to-beat variability.

NASA Discipline Cardiopulmonary↗

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Convergence of Overlapping Schwarz Decomposition for Nonlinear Optimal Control

Here, we study the convergence properties of an overlapping Schwarz decomposition algorithm for solving nonlinear optimal control problems (OCPs). The algorithm decomposes the time domain into a set of overlapping subdomains, and solves all subproblems defined over subdomains in parallel. The convergence is attained by updating primal-dual information at the boundaries of overlapping subdomains. We show that the algorithm exhibits local linear convergence, and that the convergence rate improves exponentially with the overlap size. We also establish global convergence results for a general quadratic programming, which enables the application of the Schwarz scheme inside second-order optimization algorithms (e.g., sequential quadratic programming). The theoretical foundation of our convergence analysis is a sensitivity result of nonlinear OCPs, which we call "exponential decay of sensitivity" (EDS). Intuitively, EDS states that the impact of perturbations at domain boundaries (i.e., initial and terminal time) on the solution decays exponentially as one moves into the domain. Here, we expand a previous analysis available in the literature by showing that EDS holds for both primal and dual solutions of nonlinear OCPs, under uniform second-order sufficient condition, controllability condition, and boundedness condition. We conduct experiments with a quadrotor motion planning problem and a partial differential equations (PDE) control problem to validate our theory, and show that the approach is significantly more efficient than alternating direction method of multipliers and as efficient as the centralized interior-point solver.

42 ENGINEERING↗

Using the NASTRAN Thermal Analyzer to simulate a flight scientific instrument package

The NASTRAN Thermal Analyzer has proven to be a unique and useful tool for thermal analyses involving large and complex structures where small, thermally induced deformations are critical. Among its major advantages are direct grid point-to-grid point compatibility with large structural models; plots of the model that may be generated for both conduction and boundary elements; versatility of applying transient thermal loads especially to repeat orbital cycles; on-line printer plotting of temperatures and rate of temperature changes as a function of time; and direct matrix input to solve linear differential equations on-line. These features provide a flexibility far beyond that available in most finite-difference thermal analysis computer programs.

Lee, H.-P.↗