Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic derivatives”

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

Nonlinear aerodynamic modeling using multivariate orthogonal functions

A technique was developed for global modeling of nonlinear aerodynamic coefficients using multivariate orthogonal functions based on the data. Each orthogonal function retained in the model was decomposed into an expansion of ordinary polynomials in the independent variables, so that the final model could be interpreted as selectively retained terms from a multivariable power series expansion. A predicted squared-error metric was used to determine the orthogonal functions to be retained in the model; analytical derivatives were easily computed. The approach was demonstrated on the Z-body axis aerodynamic force coefficient (Cz) wind tunnel data for an F-18 research vehicle which came from a tabular wind tunnel and covered the entire subsonic flight envelope. For a realistic case, the analytical model predicted experimental values of Cz very well. The modeling technique is shown to be capable of generating a compact, global analytical representation of nonlinear aerodynamics. The polynomial model has good predictive capability, global validity, and analytical differentiability.

Morelli, Eugene A.↗

Linearized transfer between inclined circular orbits using low-thrust blow down propulsion system

Noncoplanar transfers between neighboring circular orbits are presented for spacecraft using their own low-thrust blow down propulsion system. It is assumed that the out-of-plane angle between the decaying thrust vector and the current orbit plane remains constant for each extended burn. Switching conditions are derived for the cutoff and relight of the propulsion system in order to carry out a given transfer with inclination change. Furthermore the location where the thrust acceleration is initially applied with respect to the line of nodes of the two orbits is uniquely determined. Finally an analytic derivation of the linearized coplanar motion for stationkeeping and terminal rendezvous studies is also presented and a scheme for deriving the second order correction shown.

Kechichian, J. A.↗

Examination of Spin-Orbit Resonance in Eccentric and Low Altitude Mars Orbits

Low orbit perturbations derive from mass concentrations and show high sensitivity to certain initial orbit parameters and spin-orbit resonance con-ditions. Evaluation of spin-orbit resonance and gravity perturbation effects are important in mission design applications. Here we examine low altitude and polar orbits of Mars. The method uses analytical derivations to identify regimes of interest and computational modeling and graphical design aids for detailed investigation.

Mars↗

Interactions between long and synoptic-scale waves. I - Instability of a nonzonal flow. II - Growth rate of long waves

The interactions between synoptic-scale and long waves were investigated analytically. First, the influence of the long wave on the synoptic-scale wave is examined. The structure of a synoptic-scale wave growing on a wavy basic scale was analyzed under an assumption that the synoptic-scale waves have the structure of the most unstable normal modes. The derived analytical solution, which is simple and is amenable to physical interpretation, can be understood in terms of eddies and their local growth rate. The analytical solution is then used to examine growth of a long wave in the presence of parameterized synoptic-scale waves. Two possibly unstable solutions were found; one is a modification of the linear long wave, and the other a strongly nonlinear solution. In both cases, the synoptic-scale wave increases the growth rate of the long wave.

Ebisuzaki, Wesley↗

Induced anisotropy of thermal conductivity of polymer solids under large strains

A simple analytical form of induced anisotropy of heat conductivity of initially isotropic polymer solids results from employing the simplified theory of the three-chain model of the non-Gaussian network. The analytical form appears to be valid up to a stretch ratio of 2.65, which is the limit of existing experimental data. The effect of induced anisotropy on the temperature distribution, due to the large deformations, is illustrated for a highly expanded spherical shell and a cylindrical tube under a steady-state heat flow using the derived analytical form of the strain-dependent heat conductivity.

Peng, S. T. J.↗

The Development of Cambered Airfoil Sections Having Favorable Lift Characteristics at Supercritical Mach Numbers

Several groups of new airfoil sections, designated as the NACA 8-series, are derived analytically to have lift characteristics at supercritical Mach numbers which are favorable in the sense that the abrupt loss of lift, characteristic of the usual airfoil section at Mach numbers above the critical, is avoided. Aerodynamic characteristics determined, from two-dimensional windtunnel tests at Mach numbers up to approximately 0.9 are presented for each of the derived airfoils. Comparisons are made between the characteristics of these airfoils and the corresponding characteristics of representative NPiCA 6-series airfoils. The experimental results confirm the design expectations in demonstrating for the NACA S-series airfoils either no variation, or an Increase from the low-speed design value, In the lift coefficient at a constant angle of attack with increasing Mach number above the critical. It was not found possible to improve the variation with Mach number of the slope of the lift curve for these airfoils above that for the NACA 6-series airfoils. The drag characteristics of the new airfoils are somewhat inferior to those of the NACA 6- series with respect to divergence with Mach number, but the pitching-moment characteristics are more favorable for the thinner new sections In demonstrating somewhat smaller variations of moment coefficient with both angle of attack and Mach number. The effect on the aero&ynamic characteristics at high Mach numbers of removing the cusp from the trailing-edge regions of two 10-percent-chord-thick NACA 6-series airfoils is determined to be negligible.

AIRFOILS - NACA 836B110 - CHAR , MOMENTS & ORDINAT↗

The development of cambered airfoil sections having favorable lift characteristics at supercritical Mach numbers

Several groups of new airfoil sections, designated as the NACA 8-series, are derived analytically to have lift characteristics at supercritical Mach numbers which are favorable in the sense that the abrupt loss of lift, characteristic of the usual airfoil section at Mach numbers above the critical, is avoided. Aerodynamic characteristics determined from two-dimensional wind-tunnel tests at Mach numbers up to approximately 0.9 are presented for each of the derived airfoils. Comparisons are made between the characteristics of these airfoils and the corresponding characteristics of representative NACA 6-series airfoils.

Graham, Donald J↗

Phase Function Determination in Support of Orbital Debris Size Estimation

To recover the size of a space debris object from photometric measurements, it is necessary to determine its albedo and basic shape: if the albedo is known, the reflective area can be calculated; and if the shape is known, the shape and area taken together can be used to estimate a characteristic dimension. Albedo is typically determined by inferring the object s material type from filter photometry or spectroscopy and is not the subject of the present study. Object shape, on the other hand, can be revealed from a time-history of the object s brightness response. The most data-rich presentation is a continuous light-curve that records the object s brightness for an entire sensor pass, which could last for tens of minutes to several hours: from this one can see both short-term periodic behavior as well as brightness variations with phase angle. Light-curve interpretation, however, is more art than science and does not lend itself easily to automation; and the collection method, which requires single-object telescope dedication for long periods of time, is not well suited to debris survey conditions. So one is led to investigate how easily an object s brightness phase function, which can be constructed from the more survey-friendly point photometry, can be used to recover object shape. Such a recovery is usually attempted by comparing a phase-function curve constructed from an object s empirical brightness measurements to analytically-derived curves for basic shapes or shape combinations. There are two ways to accomplish this: a simple averaged brightness-versus phase curve assembled from the empirical data, or a more elaborate approach in which one is essentially calculating a brightness PDF for each phase angle bin (a technique explored in unpublished AFRL/RV research and in Ojakangas 2011); in each case the empirical curve is compared to analytical results for shapes of interest. The latter technique promises more discrimination power but requires more data; the former can be assembled in its essentials from fewer measurements but will be less definitive in its assignments. The goal of the present study is to evaluate both techniques under debris survey conditions to determine their relative performance and, additionally, to learn precisely how a survey should be conducted in order to maximize their performance. Because the distendedness of objects has more of an effect than their precise shape in calculating a characteristic dimension, one is interested in the techniques discrimination ability to distinguish between an elongated rectangular prism and a short rectangular prism or cube, or an elongated cylinder from a squat cylinder or sphere. Sensitivity studies using simulated data will be conducted to determine discrimination power for both techniques as a function of amount of data collected and range (and specific region) of phase angles sampled. Empirical GEODSS photometry data for distended objects (dead payloads with solar panels, rocket bodies) and compact objects (cubesats, calibration spheres, squat payloads) will also be used to test this discrimination ability. The result will be a recommended technique and data collection paradigm for debris surveys in order to maximize this type of discrimination.

Hejduk, M. D.↗

Oxidation kinetics of a continuous carbon phase in a nonreactive matrix

Analytical solutions of and experimental results on the oxidation kinetics of carbon in a pore are presented. Reaction rate, reaction sequence, oxidant partial pressure, total system pressure, pore/crack dimensions, and temperature are analyzed with respect to the influence of each on an overall linear-parabolic rate relationship. Direct measurement of carbon recession is performed using two microcomposite model systems oxidized in the temperature range of 700 to 1200 C, and for times to 35 h. Experimental results are evaluated using the derived analytical solutions. Implications on the oxidation resistance of continuous-fiber-reinforced ceramic-matrix composites containing a carbon constituent are discussed.

Eckel, Andrew J.↗

Electron Beam Return-Current Losses in Solar Flares: Initial Comparison of Analytical and Numerical Results

Accelerated electrons play an important role in the energetics of solar flares. Understanding the process or processes that accelerate these electrons to high, nonthermal energies also depends on understanding the evolution of these electrons between the acceleration region and the region where they are observed through their hard X-ray or radio emission. Energy losses in the co-spatial electric field that drives the current-neutralizing return current can flatten the electron distribution toward low energies. This in turn flattens the corresponding bremsstrahlung hard X-ray spectrum toward low energies. The lost electron beam energy also enhances heating in the coronal part of the flare loop. Extending earlier work by Knight & Sturrock (1977), Emslie (1980), Diakonov & Somov (1988), and Litvinenko & Somov (1991), I have derived analytical and semi-analytical results for the nonthermal electron distribution function and the self-consistent electric field strength in the presence of a steady-state return-current. I review these results, presented previously at the 2009 SPD Meeting in Boulder, CO, and compare them and computed X-ray spectra with numerical results obtained by Zharkova & Gordovskii (2005, 2006). The phYSical significance of similarities and differences in the results will be emphasized. This work is supported by NASA's Heliophysics Guest Investigator Program and the RHESSI Project.

Holman, Gordon↗

Estimability of geodetic parameters from space VLBI observables

The feasibility of space very long base interferometry (VLBI) observables for geodesy and geodynamics is investigated. A brief review of space VLBI systems from the point of view of potential geodetic application is given. A selected notational convention is used to jointly treat the VLBI observables of different types of baselines within a combined ground/space VLBI network. The basic equations of the space VLBI observables appropriate for convariance analysis are derived and included. The corresponding equations for the ground-to-ground baseline VLBI observables are also given for a comparison. The simplified expression of the mathematical models for both space VLBI observables (time delay and delay rate) include the ground station coordinates, the satellite orbital elements, the earth rotation parameters, the radio source coordinates, and clock parameters. The observation equations with these parameters were examined in order to determine which of them are separable or nonseparable. Singularity problems arising from coordinate system definition and critical configuration are studied. Linear dependencies between partials are analytically derived. The mathematical models for ground-space baseline VLBI observables were tested with simulation data in the frame of some numerical experiments. Singularity due to datum defect is confirmed.

Adam, Jozsef↗

Estimation of dynamic rotor loads for the rotor systems research aircraft: Methodology development and validation

The Rotor Systems Research Aircraft uses load cells to isolate the rotor/transmission systm from the fuselage. A mathematical model relating applied rotor loads and inertial loads of the rotor/transmission system to the load cell response is required to allow the load cells to be used to estimate rotor loads from flight data. Such a model is derived analytically by applying a force and moment balance to the isolated rotor/transmission system. The model is tested by comparing its estimated values of applied rotor loads with measured values obtained from a ground based shake test. Discrepancies in the comparison are used to isolate sources of unmodeled external loads. Once the structure of the mathematical model has been validated by comparison with experimental data, the parameters must be identified. Since the parameters may vary with flight condition it is desirable to identify the parameters directly from the flight data. A Maximum Likelihood identification algorithm is derived for this purpose and tested using a computer simulation of load cell data. The identification is found to converge within 10 samples. The rapid convergence facilitates tracking of time varying parameters of the load cell model in flight.

Duval, R. W.↗

A Three-Dimensional Analytic Model for the Scattering of a Spherical Bush

Advanced climate models require a more realistic description of canopy radiation with reasonable computational efficiency. This paper develops the mathematics of scattering from a spherical object conceptualized to be a spherical bush to provide a building block that helps to address this need of climate models. It is composed of a homogeneous distribution of individual smaller objects that scatter isotropically. In the limit of small optical depth, incident radiation will scatter isotropically as the sum of that scattered by all the individual scatterers, but at large optical depth the radiation leaving the spherical bush in a given direction is reduced by mutual shadowing of the smaller objects. In the single scattering limit, the scattering phase function and so the albedo are obtained by simple but accurate analytic expressions derived from analytic integration and numerical evaluation. Except in the limit of thin canopies, the scattering and hence albedos are qualitatively and quantitatively different than those derived from 1-D modeling.

Dickinson, Robert E.↗

Low-Earth Orbit Trajectory Optimization in the Presence of Atmospheric Uncertainty

The previous 20 to 25 years have seen a tremendous increase in space exploration, and with that an increase in the level of logistics planning needed to ensure mission success. For spacecraft that are designed to be periodically re-supplied, a key logistics consumable is propellant, as it constitutes the greatest up-mass on re-supply vehicles. A trajectory design strategy is therefore desired that minimizes propellant usage in order to ease the demand for propellant re-supply missions. This thesis develops such a strategy in three stages, and uses the International Space Station (ISS) as its testbed, as no other LEO spacecraft is more challenging from a space logistics standpoint. First, the ISS trajectory planning problem is formulated as a constrained burn optimization problem assuming a deterministic atmosphere. The cost function is total ∆v, with constraints imposed on longitude of ascending viii node (LAN) and semi-major axis (SMA) altitude. Analytic derivatives are constructed for both the cost and constraints, which are necessary given the 6-week to 2-year time frames being considered. A gradient-based optimizer is then utilized to find locally-optimal solutions to real-world ISS trajectory planning problems. Second, atmospheric uncertainty is addressed by constructing a probabilistic model of space weather data using Gaussian Processes (GPs). Bayesian inference is performed using the GP model to generate mean and covariance estimates for space weather predictions, whose pedigree is assessed against test data. The predictions are then mapped into atmospheric density via the analytic Jacchia-Roberts density model, and the effect of space weather uncertainty on orbital lifetime is examined. Third, an ISS burn execution uncertainty model is developed. This model, along with the space weather uncertainty model, are deployed in a linear covariance analysis to ascertain their combined effect on LAN and SMA altitude dispersions. The deterministic constraints from the original problem are re-formulated as stochastic constraints, where now the constraint uncertainty interval is required to fall within specified bounds. An updated optimization framework is constructed using the original ∆v cost function along with the stochastic constraints to solve the trajectory optimization problem under atmospheric uncertainty. Finally, the complete architecture is summarized for deployment in an operational setting.

Trajectory Optimization↗

Probabilistic evaluation of on-line checks in fault-tolerant multiprocessor systems

The analysis of fault-tolerant multiprocessor systems that use concurrent error detection (CED) schemes is much more difficult than the analysis of conventional fault-tolerant architectures. Various analytical techniques have been proposed to evaluate CED schemes deterministically. However, these approaches are based on worst-case assumptions related to the failure of system components. Often, the evaluation results do not reflect the actual fault tolerance capabilities of the system. A probabilistic approach to evaluate the fault detecting and locating capabilities of on-line checks in a system is developed. The various probabilities associated with the checking schemes are identified and used in the framework of the matrix-based model. Based on these probabilistic matrices, estimates for the fault tolerance capabilities of various systems are derived analytically.

Nair, V. S. S.↗

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗

Spacecraft attitude control momentum requirements analysis

The relationship between attitude and angular momentum control requirements is derived for a fixed attitude, Earth orbiting spacecraft with large area articulating appendages. Environmental effects such as gravity gradient, solar radiation pressure, and aerodynamic forces arising from a dynamic, rotating atmosphere are examined. It is shown that, in general, each environmental effect contributes to both cyclic and secular momentum requirements both within and perpendicular to the orbit plane. The gyroscopic contribution to the angular momentum control requirements resulting from a rotating, Earth oriented spacecraft is also discussed. Special conditions are described where one or more components of the angular momentum can be made to vanish, or become purely cyclical. Computer generated plots for a candidate space station configuration are presented to supplement the analytically derived results.

Robertson, Brent P.↗

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗