Search NASASearch

SEARCH · Search NASA

Results for “path integral formulation”

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 19 records

Coherent laser radar performance for general atmospheric refractive turbulence

A general theory for the signal-to-noise ratio (SNR) of a coherent detection laser radar is developed using the path-integral formulation (Fresnel approximation), which is valid for any typical path-integrated atmospheric refractive turbulence. The principal effects of refractive turbulence are discussed, and analytical expressions are presented for the case of untruncated Gaussians for the transmitted field, local oscillator field, and transmitter/receiver optics. The physical mechanisms that reduce heterodyne efficiency are identified.

Frehlich, Rod G.

Description of quantum noise by a Langevin equation

General features of the quantum noise problem expressed as the equations of motion for a particle coupled to a set of oscillators are investigated analytically. Account is taken of the properties of the companion oscillators by formulating quantum statistical correlation Langevin equations (QSLE). The frequency of the oscillators is then retained as a natural cut-off for the quantum noise. The QSLE is further extended to encompass the particle trajectory and is bounded by initial and final states of the oscillator. The states are expressed as the probability of existence at the moment of particle collision that takes the oscillator into a final state. Two noise sources then exist: a statistical uncertainty of the initial state and the quantum dynamical uncertainty associated with a transition from the initial to final state. Feynman's path-integral formulation is used to characterize the functional of the particle trajectory, which slows the particle. It is shown that the energy loss may be attributed to friction, which satisfies energy conservation laws.

Metiu, H.

Physics-Based Methods of Failure Analysis and Diagnostics in Human Space Flight

The Integrated Health Management (IHM) for the future aerospace systems requires to interface models of multiple subsystems in an efficient and accurate information environment at the earlier stages of system design. The complexity of modern aeronautic and aircraft systems (including e.g. the power distribution, flight control, solid and liquid motors) dictates employment of hybrid models and high-level reasoners for analysing mixed continuous and discrete information flow involving multiple modes of operation in uncertain environments, unknown state variables, heterogeneous software and hardware components. To provide the information link between key design/performance parameters and high-level reasoners we rely on development of multi-physics performance models, distributed sensors networks, and fault diagnostic and prognostic (FD&P) technologies in close collaboration with system designers. The main challenges of our research are related to the in-flight assessment of the structural stability, engine performance, and trajectory control. The main goal is to develop an intelligent IHM that not only enhances components and system reliability, but also provides a post-flight feedback helping to optimize design of the next generation of aerospace systems. Our efforts are concentrated on several directions of the research. One of the key components of our strategy is an innovative approach to the diagnostics/prognostics based on the real time dynamical inference (DI) technologies extended to encompass hybrid systems with hidden state trajectories. The major investments are into the multiphysics performance modelling that provides an access of the FD&P technologies to the main performance parameters of e.g. solid and liquid rocket motors and composite materials of the nozzle and case. Some of the recent results of our research are discussed in this chapter. We begin by introducing the problem of dynamical inference of stochastic nonlinear models and reviewing earlier results. Next, we present our analytical approach to the solution of this problem based on the path integral formulation. The resulting algorithm does not require an extensive global search for the model parameters, provides optimal compensation for the effects of dynamical noise, and is robust for a broad range of dynamical models. In the following Section the strengths of the algorithm are illustrated illustrated by inferring the parameters of the stochastic Lorenz system and comparing the results with those of earlier research. Next, we discuss a number of recent results in application to the development of the IHM for aerospace system. Firstly, we apply dynamical inference approach to a solution of classical three tank problems with mixed unknown continuous and binary parameters. The problem is considered in the context of ground support system for filling fuel tanks of liquid rocket motors. It is shown that the DI algorithm is well suited for successful solution of a hybrid version of this benchmark problem even in the presence of additional periodic and stochastic perturbation of unknown strength. Secondly, we illustrate our approach by its application to an analysis of the nozzle fault in a solid rocket motor (SRM). The internal ballistics of the SRM is modelled as a set of one-dimensional partial differential equations coupled to the dynamics of the propellant regression. In this example we are specifically focussed on the inference of discrete and continuous parameters of the nozzle blocking fault and on the possibility of an application of the DI algorithm to reducing the probability of "misses" of an on-board FD&P for SRM. In the next section re-contact problem caused by first stage/upper stage separation failure is discussed. The reaction forces imposed on the nozzle of the upper stage during the re-contact and their connection to the nozzle damage and to the thrust vector control (TVC) signal are obtained. It is shown that transient impact induced torquean be modelled as a response of an effective damped oscillator. A possible application of the DI algorithm to the inference of damage parameters and predicting fault dynamics ahead of time using the actuator signal is discussed. Finally, we formulate Bayesian inferential framework for development of the IHM system for in-flight structural health monitoring (SHM) of composite materials. We consider the signal generated by piezoelectric actuator mounted on composite structure generating elastic waves in it. The signal received by the sensor is than compared with the baseline signal. The possibility of damage inference is discussed in the context of development of the SHM.

Smelyanskiy, Vadim N.

Path Planning: Differential Dynamic Programming and Model Predictive Path Integral Control on VTOL Aircraft

This paper explores two optimal control approaches, widely used in robotics, to establish their viability as real-time trajectory planners for vehicle configurations envisioned for the emerging aviation sector of Urban Air Mobility (UAM). Differential Dynamic Programming (DDP) enables planning over highly nonlinear dynamics using second-order approximations along a nominal trajectory, and displays quadratic convergence to a local solution. Model Predictive Path Integral (MPPI) is a stochastic sampling-based algorithm that can optimize for general cost criteria, including potentially highly nonlinear formulations, and supports parallel computation through the use of modern GPU hardware. In this work, DDP and MPPI were implemented using model predictive control (MPC), and the results indicate they are able to successfully transition the aircraft over different flight envelopes and generate trajectories unique to UAM vehicles.

Differential Dynamic Programming

New Millennium EO3

The New Millennium Program (NMP) is a NASA technology program that focuses on the validation of advanced spacecraft and instrumentation technologies in space. This program specifically seeks technologies that could significantly benefit future space and Earth science missions by enabling new science capabilities and reducing life cycle costs. These technologies must also require a validation in space to mitigate risks to the first science users, and provide cross-cutting benefits to both NASA's Earth and Space Science enterprises. The NASA Office of Earth Science (OES) directed the NMP to focus the third Earth Observing mission, E03, on innovative measurement concepts that would facilitate remote sensing observations from orbits beyond conventional low-Earth orbit (LEO). These orbits include geosynchronous orbits, highly elliptical orbits, mid-Earth and high-Earth orbits, and other unique vantage points such as L1 and L2. To maximize the input from the Earth science community, a NASA Research Announcement (NRA) was released to solicit innovative measurement concepts for this NMP flight. Because the NMP is a technology validation program, rather than a conventional science program, the NRA required that these measurement concepts employ revolutionary technologies and/or measurement strategies that will enable future science missions from orbits beyond LEO. Another requirement was that a validation in space was needed to reduce real or perceived risks of this concept to future science users. The proposals submitted in response to this NRA were peer reviewed by the NASA OES. The measurement concepts selected through this process will be summarized in this presentation. The E03 measurement concept NRA did not solicit complete mission concepts or flight hardware. Instead, the selected investigators will join integrated project formulation teams to define the mission for the demonstration of the measurement technique and participate in mission design trades and implementation planning. The teams will define and document the measurement approach, the required technologies, the validation plan, the scalability of the design to future science missions, and the infusion path for the completely validated instrument. During this project formulation phase, the specific technologies required to implement each measurement concept will be solicited through an open technology solicitation, and the technology providers will be added to the integrated teams. Once a complete candidate mission concept has been formulated for each measurement concept, these concepts will be subjected to a comprehensive review. The Associate Administrator of the OES will use the results of this review to select a single concept for the E03 flight opportunity. That mission concept will enter its implementation phase, which is expected to last for less than 36 months before launch. During this phase, and after launch, the implementation team will work with the NMP to document and disseminate the information gathered during the flight qualification testing and operations of the measurement concept and its component technologies. This information will be archived and published, where appropriate, to encourage the rapid infusion of these validated technologies into future space and Earth science missions. I

Li, Fuk K.

Elevated Temperature Crack Propagation

This paper is a summary of two NASA contracts on high temperature fatigue crack propagation in metals. The first evaluated the ability of fairly simple nonlinear fracture parameters to correlate crack propagation. Hastelloy-X specimens were tested under isothermal and thermomechanical cycling at temperatures up to 980 degrees C (1800 degrees F). The most successful correlating parameter was the crack tip opening displacement derived from the J-integral. The second evaluated the ability of several path-independent integrals to correlate crack propagation behavior. Inconel 718 specimens were tested under isothermal, thermomechanical, temperature gradient, and creep conditions at temperatures up to 650 degrees C (1200 degrees F). The integrals formulated by Blackburn and by Kishimoto correlated the data reasonably well under all test conditions.

Orange, Thomas W.

Fast algorithms for visualizing fluid motion in steady flow on unstructured grids

The plotting of streamlines is an effective way of visualizing fluid motion in steady flows. Additional information about the flowfield, such as local rotation and expansion, can be shown by drawing in the form of a ribbon or tube. In this paper, we present efficient algorithms for the construction of streamlines, streamribbons and streamtubes on unstructured grids. A specialized version of the Runge-Kutta method has been developed to speed up the integration of particle paths. We have also derived closed-form solutions for calculating angular rotation rate and radius to construct streamribbons and streamtubes, respectively. According to our analysis and test results, these formulations are two to four times better in performance than previous numerical methods. As a large number of traces are calculated, the improved performance could be significant.

Ueng, S. K.

Vector formulation for interferogram surface fitting

Interferometry is an optical testing technique that quantifies the optical path difference (OPD) between a reference wave front and a test wave front based on the interference of light. Fringes are formed when the OPD is an integral multiple of the illuminating wavelength. The resultant two-dimensional pattern is called an interferogram. The function of any interferogram analysis program is to extract this OPD and to produce a representation of the test wave front (or surface). This is accomplished through a three-step process of sampling, ordering, and fitting. We develop a generalized linear-algebra vector-notation model of the interferogram sampling and fitting process.

Fischer, David J.

A Fully Associative, Non-Linear Kinematic, Unified Viscoplastic Model for Titanium Based Matrices

Specific forms for both the Gibb's and complementary dissipation potentials are chosen such that a complete (i.e., fully associative) potential based multiaxial unified viscoplastic model is obtained. This model possesses one tensorial internal state variable that is associated with dislocation substructure, with an evolutionary law that has nonlinear kinematic hardening and both thermal and strain induced recovery mechanisms. A unique aspect of the present model is the inclusion of non-linear hardening through the use of a compliance operator, derived from the Gibb's potential, in the evolution law for the back stress. This non-linear tensorial operator is significant in that it allows both the flow and evolutionary laws to be fully associative (and therefore easily integrated) and greatly influences the multiaxial response under non-proportional loading paths. In addition to this nonlinear compliance operator, a new consistent, potential preserving, internal strain unloading criterion has been introduced to prevent abnormalities in the predicted stress-strain curves, which are present with nonlinear hardening formulations, during unloading and reversed loading of the external variables. Specification of an experimental program for the complete determination of the material functions and parameters for characterizing a metallic matrix, e.g., TIMETAL 21S, is given. The experiments utilized are tensile, creep, and step creep tests. Finally, a comparison of this model and a commonly used Bodner-Partom model is made on the basis of predictive accuracy and numerical efficiency.

Arnold, S. M.

Models for infrared atmospheric radiation

Line and band models for infrared spectral absorption are discussed. Radiative transmittance and integrated absorptance of Lorentz, Doppler, and voigt line profiles were compared for a range of parameters. It was found that, for the intermediate path lengths, the combined Lorentz-Doppler (Voigt) profile is essential in calculating the atmospheric transmittance. Narrow band model relations for absorptance were used to develop exact formulations for total absorption by four wide band models. Several continuous correlations for the absorption of a wide band model were compared with the numerical solutions of the wide band models. By employing the line-by-line and quasi-random band model formulations, computational procedures were developed for evaluating transmittance and upwelling atmospheric radiance. Homogeneous path transmittances were calculated for selected bands of CO, CO2, and N2O and compared with experimental measurements. The upwelling radiance and signal change in the wave number interval of the CO fundamental band were also calculated.

Tiwari, S. N.

PLANETARY DEFENSE: FINDINGS AND RECOMMENDATIONS FROM THE NATIONAL ACADEMIES PLANETARY SCIENCE AND ASTROBIOLOGY DECADAL SURVEY 2023 - 2032

The National Academies of Sciences, Engineering, and Medicine was chartered by the NASA Science Mission Directorate (SMD), with support from the National Science Foundation, to conduct a new decadal survey to develop a comprehensive science and mission strategy for planetary science (1). The survey was tasked to determine the current state of knowledge and to identify the most important scientific questions to be addressed in astrobiology and planetary science during the upcoming decade of 2023 - 2032. For the first time, the decadal survey also included and addressed the topic of planetary defense, now recognized as a key element of NASA’s planetary science activities under the direction of NASA’s SMD Planetary Defense Coordination Office. The topic of planetary defense was included in the decadal survey as a dedicated chapter, divided into five sections based on the framework from the U.S. National NEO Preparedness Strategy and Action Plan (2). These sections discussed: NEO Detection, Tracking, and Characterization NEO Modelling, Prediction, and Integration NEO Deflection and Disruption Missions International Cooperation on NEO Preparation NEO Impact Emergency Procedures and Action Protocols Content of these sections included identifying the activities, ground- and space-based assets, facilities, resources, technologies, and missions that should be considered to enable NEO preparedness over the next decade. Supporting inputs to this chapter involved contributions from submitted white papers authored by many members of the NEO science and planetary defense communities. In addition, a planetary defense rapid mission architecture study commissioned by the decadal survey small bodies panel and led by the Jet Propulsion Laboratory examined several different types of in-situ characterization and mitigation demonstration missions. All of this input contributed to the formulation of 42 findings and 11 recommendations that suggest a path forward for increasing our understanding of the NEO population and their associated physical characteristics, modelling the impact probabilities and hazard potential of these objects, and implementing plans to develop and test NEO impact mitigation capabilities. Including planetary defense in the decadal survey has helped sharpen the focus on this international issue and provided key inputs for stakeholders and decision makers to consider when developing future planetary defense-related activities. Such information is vital for developing a long-term integrated planetary defense capability able to protect humanity in the decades to come.

NEO Detection and Characterization

Thermal Testing and Integration: Magnetospheric MultiScale (MMS) Observatories with Digital 1-Wire Sensors

Thermocouples require two thin wires to be routed out of the spacecraft to connect to the ground support equipment used to monitor and record the temperature data. This large number of wires that exit the observatory complicates integration and creates an undesirable heat path during testing. These wires exiting the spacecraft need to be characterized as a thermal short that will not exist during flight. To minimize complexity and reduce thermal variables from these ground support equipment (GSE) wires, MMS pursued a hybrid path for temperature monitoring, utilizing thermocouples and digital 1-wire temperature sensors. Digital 1-wire sensors can greatly reduce harness mass, length and complexity as they can be spliced together. For MMS, 350 digital 1-wire sensors were installed on the spacecraft with only 18 wires exiting as opposed to a potential 700 thermocouple wires. Digital 1-wire sensors had not been used in such a large scale at NASAGSFC prior to the MMS mission. During the MMS thermal vacuum testing a lessons learned matrix was formulated that will assist future integration of 1-wires into thermal testing and one day into flight.

Sensors

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.

Dynamic Flow Management Problems in Air Transportation

In 1995, over six hundred thousand licensed pilots flew nearly thirty-five million flights into over eighteen thousand U.S. airports, logging more than 519 billion passenger miles. Since demand for air travel has increased by more than 50% in the last decade while capacity has stagnated, congestion is a problem of undeniable practical significance. In this thesis, we will develop optimization techniques that reduce the impact of congestion on the national airspace. We start by determining the optimal release times for flights into the airspace and the optimal speed adjustment while airborne taking into account the capacitated airspace. This is called the Air Traffic Flow Management Problem (TFMP). We address the complexity, showing that it is NP-hard. We build an integer programming formulation that is quite strong as some of the proposed inequalities are facet defining for the convex hull of solutions. For practical problems, the solutions of the LP relaxation of the TFMP are very often integral. In essence, we reduce the problem to efficiently solving large scale linear programming problems. Thus, the computation times are reasonably small for large scale, practical problems involving thousands of flights. Next, we address the problem of determining how to reroute aircraft in the airspace system when faced with dynamically changing weather conditions. This is called the Air Traffic Flow Management Rerouting Problem (TFMRP) We present an integrated mathematical programming approach for the TFMRP, which utilizes several methodologies, in order to minimize delay costs. In order to address the high dimensionality, we present an aggregate model, in which we formulate the TFMRP as a multicommodity, integer, dynamic network flow problem with certain side constraints. Using Lagrangian relaxation, we generate aggregate flows that are decomposed into a collection of flight paths using a randomized rounding heuristic. This collection of paths is used in a packing integer programming formulation, the solution of which generates feasible and near-optimal routes for individual flights. The algorithm, termed the Lagrangian Generation Algorithm, is used to solve practical problems in the southwestern portion of United States in which the solutions are within 1% of the corresponding lower bounds.

Patterson, Sarah Stock

A mixed-mode crack analysis of rectilinear anisotropic solids using conservation laws of elasticity

A very simple and convenient method of analysis for studying two-dimensional mixed-mode crack problems in rectilinear anisotropic solids is presented. The analysis is formulated on the basis of conservation laws of anisotropic elasticity and of fundamental relationships in anisotropic fracture mechanics. The problem is reduced to a system of linear algebraic equations in mixed-mode stress intensity factors. One of the salient features of the present approach is that it can determine directly the mixed-mode stress intensity solutions from the conservation integrals evaluated along a path removed from the crack-tip region without the need of solving the corresponding complex near-field boundary value problem. Several examples with solutions available in the literature are solved to ensure the accuracy of the current analysis. This method is further demonstrated to be superior to other approaches in its numerical simplicity and computational efficiency. Solutions of more complicated and practical engineering problems dealing with the crack emanating from a circular hole in composites are presented also to illustrate the capacity of this method.

Wang, S. S.

A fully associative, nonisothermal, nonlinear kinematic, unified viscoplastic model for titanium alloys

Specific forms for both the Gibb's and complementary dissipation potentials are chosen such that a complete (i.e., fully associative) potential base multiaxial, nonisothermal unified viscoplastic model is obtained. This model possesses one tensorial internal state variable (that is, associated with dislocation substructure) and an evolutionary law that has nonlinear kinematic hardening and both thermal and strain induced recovery mechanisms. A unique aspect of the present model is the inclusion of nonlinear hardening through the use of a compliance operator, derived from the Gibb's potential, in the evolution law for the back stress. This nonlinear tensorial operator is significant in that it allows both the flow and evolutionary laws to be fully associative (and therefore easily integrated), greatly influences the multiaxial response under non-proportional loading paths, and in the case of nonisothermal histories, introduces an instantaneous thermal softening mechanism proportional to the rate of change in temperature. In addition to this nonlinear compliance operator, a new consistent, potential preserving, internal strain unloading criterion has been introduced to prevent abnormalities in the predicted stress-strain curves, which are present with nonlinear hardening formulations, during unloading and reversed loading of the external variables. The specific model proposed is characterized for a representative titanium alloy commonly used as the matrix material in SiC fiber reinforced composites, i.e., TIMETAL 21S. Verification of the proposed model is shown using 'specialized' non-standard isothermal and thermomechanical deformation tests.

Arnold, S. M.

Hypersonic vehicle control law development using H infinity and mu-synthesis

Applicability and effectiveness of robust control techniques to a single-stage-to-orbit (SSTO) airbreathing hypersonic vehicle on an ascent accelerating path and their effectiveness are explored in this paper. An SSTO control system design problem, requiring high accuracy tracking of velocity and altitude commands while limiting angle of attack oscillations, minimizing control power usage and stabilizing the vehicle all in the presence of atmospheric turbulence and uncertainty in the system, was formulated to compare results of the control designs using H infinity and mu-synthesis procedures. The math model, an integrated flight/propulsion dynamic model of a conical accelerator class vehicle, was linearized as the vehicle accelerated through Mach 8. Controller analysis was conducted using the singular value technique and the mu-analysis approach. Analysis results were obtained in both the frequency and the time domains. The results clearly demonstrate the inherent advantages of the structured singular value framework for this class of problems. Since payload performance margins are so critical for the SSTO mission, it is crucial that adequate stability margins be provided without sacrificing any payload mass.

Gregory, Irene M.

Methods of Attenuation Correction for Dual-Wavelength and Dual-Polarization Weather Radar Data

In writing the integral equations for the median mass diameter and number concentration, or comparable parameters of the raindrop size distribution, it is apparent that the forms of the equations for dual-polarization and dual-wavelength radar data are identical when attenuation effects are included. The differential backscattering and extinction coefficients appear in both sets of equations: for the dual-polarization equations, the differences are taken with respect to polarization at a fixed frequency while for the dual-wavelength equations, the differences are taken with respect to frequency at a fixed polarization. An alternative to the integral equation formulation is that based on the k-Z (attenuation coefficient-radar reflectivity factor) parameterization. This-technique was originally developed for attenuating single-wavelength radars, a variation of which has been applied to the TRMM Precipitation Radar data (PR). Extensions of this method have also been applied to dual-polarization data. In fact, it is not difficult to show that nearly identical equations are applicable as well to dualwavelength radar data. In this case, the equations for median mass diameter and number concentration take the form of coupled, but non-integral equations. Differences between this and the integral equation formulation are a consequence of the different ways in which attenuation correction is performed under the two formulations. For both techniques, the equations can be solved either forward from the radar outward or backward from the final range gate toward the radar. Although the forward-going solutions tend to be unstable as the attenuation out to the range of interest becomes large in some sense, an independent estimate of path attenuation is not required. This is analogous to the case of an attenuating single-wavelength radar where the forward solution to the Hitschfeld-Bordan equation becomes unstable as the attenuation increases. To circumvent this problem, the equations can be expressed in the form of a final-value problem so that the recursion begins at the far range gate and proceeds inward towards the radar. Solving the problem in this way traditionally requires estimates of path attenuation to the final gate: in the case of orthogonal linear polarizations, the attenuations at horizontal and vertical polarizations (same frequency) are required while in the dual-wavelength case, attenuations at the two frequencies (same polarization) are required.

Meneghini, R.