SEARCH · Engineering Papers

Results for “Differentiable physics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Elliptic Grid Generation of Spiral-Bevel Pinion Gear Typical of OH-58 Helicopter Transmission

This paper discusses the source term treatment in the numerical solution of elliptic partial differential equations for an interior grid generation problem in generalized curvilinear coordinates. The geometry considered is that of a planar cross-section of a generic spiral-bevel gear tooth typical of a pinion in the OH-58 helicopter transmission. The source terms used are appropriate for an interior grid domain where all the boundaries are prescribed via a combination of Dirichlet and Neumann boundary conditions. New constraints based on the Green's Theorem are derived which uniquely determine the coefficients in the source terms. These constraints are designed for boundary clustered grids where gradients in physical quantities need to be resolved adequately. However, it is seen that the present formulation works satisfactorily for mild clustering also. Thus, a fully automated elliptic grid generation technique is made possible where there is no need for a parametric study of these parameters since the new relations fix these free parameters uniquely.

Kaul, Upender K.↗

NASA NTRS2010

Defining an Abrasion Index for Lunar Surface Systems as a Function of Dust Interaction Modes and Variable Concentration Zones

Unexpected issues were encountered during the Apollo era of lunar exploration due to detrimental abrasion of materials upon exposure to the fine-grained, irregular shaped dust on the surface of the Moon. For critical design features involving contact with the lunar surface and for astronaut safety concerns, operational concepts and dust tolerance must be considered in the early phases of mission planning. To systematically define material selection criteria, dust interaction can be characterized by two-body or three-body abrasion testing, and subcategorically by physical interactions of compression, rolling, sliding and bending representing specific applications within the system. Two-body abrasion occurs when a single particle or asperity slides across a given surface removing or displacing material. Three-body abrasion occurs when multiple particles interact with a solid surface, or in between two surfaces, allowing the abrasives to freely rotate and interact with the material(s), leading to removal or displacement of mass. Different modes of interaction are described in this paper along with corresponding types of tests that can be utilized to evaluate each configuration. In addition to differential modes of abrasion, variable concentrations of dust in different zones can also be considered for a given system design and operational protocol. These zones include: (1) outside the habitat where extensive dust exposure occurs, (2) in a transitional zone such as an airlock or suitport, and (3) inside the habitat or spacesuit with a low particle count. These zones can be used to help define dust interaction frequencies, and corresponding risks to the systems and/or crew can be addressed by appropriate mitigation strategies. An abrasion index is introduced that includes the level of risk, R, the hardness of the mineralogy, H, the severity of the abrasion mode, S, and the frequency of particle interactions, F.

Kobrick, Ryan L.↗

NASA NTRS2019

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

NASA NTRS1974

On the thermal evolution of the terrestrial planets

Physical and chemical constraints for such different planetary objects as the earth, the moon and meteorite parent bodies can best be satisfied by thermal history models having high initial temperatures. On the basis of thermal calculations it is suggested that the evolution of the other terrestrial planets (Mars, Venus and Mercury) was also characterized by high initial temperatures. Under these conditions, melting and, consequently, fractionation would set in at an early stage. Because of the resulting redistribution of the long-lived radioactive heat sources and the concentration of these elements in the surface layers, large-scale differentiation could be achieved by partial melting.

Fricker, P. E.↗

NASA NTRS1991

A Large-Grain Mapping Approach for Multiprocessor Systems Through Data Flow Model Ph.D. Thesis

A large-grain level mapping method is presented of numerical oriented applications onto multiprocessor systems. The method is based on the large-grain data flow representation of the input application and it assumes a general interconnection topology of the multiprocessor system. The large-grain data flow model was used because such representation best exhibits inherited parallelism in many important applications, e.g., CFD models based on partial differential equations can be presented in large-grain data flow format, very effectively. A generalized interconnection topology of the multiprocessor architecture is considered, including such architectural issues as interprocessor communication cost, with the aim to identify the 'best matching' between the application and the multiprocessor structure. The objective is to minimize the total execution time of the input algorithm running on the target system. The mapping strategy consists of the following: (1) large-grain data flow graph generation from the input application using compilation techniques; (2) data flow graph partitioning into basic computation blocks; and (3) physical mapping onto the target multiprocessor using a priority allocation scheme for the computation blocks.

Kim, Hwa-Soo↗

NASA NTRS2002

Spacewire on Earth orbiting scatterometers

The need for a high speed, reliable and easy to implement communication link has led to the development of a space flight oriented version of IEEE 1355 called SpaceWire. SpaceWire is based on high-speed (200 Mbps) serial point-to-point links using Low Voltage Differential Signaling (LVDS). SpaceWIre has provisions for routing messages between a large network of processors, using wormhole routing for low overhead and latency. {additionally, there are available space qualified hybrids, which provide the Link layer to the user's bus}. A test bed of multiple digital signal processor breadboards, demonstrating the ability to meet signal processing requirements for an orbiting scatterometer has been implemented using three Astrium MCM-DSPs, each breadboard consists of a Multi Chip Module (MCM) that combines a space qualified Digital Signal Processor and peripherals, including IEEE-1355 links. With the addition of appropriate physical layer interfaces and software on the DSP, the SpaceWire link is used to communicate between processors on the test bed, e.g. sending timing references, commands, status, and science data among the processors. Results are presented on development issues surrounding the use of SpaceWire in this environment, from physical layer implementation (cables, connectors, LVDS drivers) to diagnostic tools, driver firmware, and development methodology. The tools, methods, and hardware, software challenges and preliminary performance are investigated and discussed.

serial↗

NASA NTRS2005

Diffraction-based optical correlator

Method and system for wavelength-based processing of a light beam. A light beam, produced at a chemical or physical reaction site and having at least first and second wavelengths, ?1 and ?2, is received and diffracted at a first diffraction grating to provide first and second diffracted beams, which are received and analyzed in terms of wavelength and/or time at two spaced apart light detectors. In a second embodiment, light from first and second sources is diffracted and compared in terms of wavelength and/or time to determine if the two beams arise from the same source. In a third embodiment, a light beam is split and diffracted and passed through first and second environments to study differential effects. In a fourth embodiment, diffracted light beam components, having first and second wavelengths, are received sequentially at a reaction site to determine whether a specified reaction is promoted, based on order of receipt of the beams. In a fifth embodiment, a cylindrically shaped diffraction grating (uniform or chirped) is rotated and translated to provide a sequence of diffracted beams with different wavelengths. In a sixth embodiment, incident light, representing one or more symbols, is successively diffracted from first and second diffraction gratings and is received at different light detectors, depending upon the wavelengths present in the incident light.

Spremo, Stevan M.↗

NASA NTRS1973

A qualitative study of the complete set of solutions of the differential equation of motion of a test particle in the equatorial plane of the Kerr gravitational field

A study is made of the mathematical solution of the differential equation of motion of a test particle in the equatorial plane of the Kerr gravitational field, using S (Schwarzschild-like) coordinates. A qualitative solution of this equation leads to the conclusion that there can only be 25 different types of orbits. For each value of a, the results are presented in a master diagram for which h and e are the parameters. A master diagram divides the h, e parameter space into regions such that at each point within one of these regions the types of admissible orbits are qualitatively the same. A pictorial representation of the physical orbits in the r, phi plane is also given.

Montgomery, H. E.↗

NASA NTRS2005

Foundations of Tensor Analysis for Students of Physics and Engineering With an Introduction to the Theory of Relativity

Tensor analysis is one of the more abstruse, even if one of the more useful, higher math subjects enjoined by students of physics and engineering. It is abstruse because of the intellectual gap that exists between where most physics and engineering mathematics leave off and where tensor analysis traditionally begins. It is useful because of its great generality, computational power, and compact, easy to use, notation. This paper bridges the intellectual gap. It is divided into three parts: algebra, calculus, and relativity. Algebra: In tensor analysis, coordinate independent quantities are sought for applications in physics and engineering. Coordinate independence means that the quantities have such coordinate transformations as to leave them invariant relative to a particular observer s coordinate system. Calculus: Non-zero base vector derivatives contribute terms to dynamical equations that correspond to pseudoaccelerations in accelerated coordinate systems and to curvature or gravity in relativity. These derivatives have a specific general form in tensor analysis. Relativity: Spacetime has an intrinsic geometry. Light is the tool for investigating that geometry. Since the observed geometry of spacetime cannot be made to match the classical geometry of Euclid, Einstein applied another more general geometry differential geometry. The merger of differential geometry and cosmology was accomplished in the theory of relativity. In relativity, gravity is equivalent to curvature.

Kolecki, Joseph C.↗

NASA NTRS2013

An Extension of the Time-Spectral Method to Overset Solvers

Relative motion in the Cartesian or overset framework causes certain spatial nodes to move in and out of the physical domain as they are dynamically blanked by moving solid bodies. This poses a problem for the conventional Time-Spectral approach, which expands the solution at every spatial node into a Fourier series spanning the period of motion. The proposed extension to the Time-Spectral method treats unblanked nodes in the conventional manner but expands the solution at dynamically blanked nodes in a basis of barycentric rational polynomials spanning partitions of contiguously defined temporal intervals. Rational polynomials avoid Runge's phenomenon on the equidistant time samples of these sub-periodic intervals. Fourier- and rational polynomial-based differentiation operators are used in tandem to provide a consistent hybrid Time-Spectral overset scheme capable of handling relative motion. The hybrid scheme is tested with a linear model problem and implemented within NASA's OVERFLOW Reynolds-averaged Navier- Stokes (RANS) solver. The hybrid Time-Spectral solver is then applied to inviscid and turbulent RANS cases of plunging and pitching airfoils and compared to time-accurate and experimental data. A limiter was applied in the turbulent case to avoid undershoots in the undamped turbulent eddy viscosity while maintaining accuracy. The hybrid scheme matches the performance of the conventional Time-Spectral method and converges to the time-accurate results with increased temporal resolution.

Leffell, Joshua Isaac↗

NASA NTRSDate not supplied

A Verification Framework for Runtime Assurance of Autonomous UAS

Runtime Assurance (RTA) is a design-time architecture for safety-critical systems where an internal monitor acts upon detecting a violation of a property. The simplex architecture is an instance of RTA, where the action taken is to hand control of the overall system to a trusted controller when an untrusted one violates a safety property. Simplex RTA is emerging as a method for allowing AI/ML and other unverified software to be integrated into safety-critical applications like aircraft. To this end, the American Society for Testing and Materials (ASTM) and NASA have each published guidelines on the use of RTA in such systems. In the simplex RTA framework, a system has an advanced controller (AC) and a reversionary controller (RC). The system is allowed to operate with the AC until a runtime monitor detects that some property has been violated and then the RC takes over. Assuming that the sample rate of the monitor will detect improper functioning with enough time for the RC to correct the impending problem, and that the RC is trusted, the system will operate as intended. This use of the simplex RTA framework can allow for the integration of untrusted, but possibly more performant, controllers in a safe way. This paper presents a formalization of a simplex RTA framework in the Prototype Verification System (PVS) theorem prover using an embedding of differential dynamic logic (DDL) called Plaidypvs. A novel feature of this framework is that it can be instantiated at different levels of abstraction. This feature allows for the formal verification of a system with an untrusted black box component, such as an AI/ML controller. This paper does not address the many difficulties in deploying RTA in an industrial-level system. Instead, the focus is on the formal verification of the simplex RTA framework in the language of hybrid programs. Hybrid programs are programs that include both discrete and continuous dynamics and can be used to model complex cyber-physical systems. Plaidypvs is a tool that enables formalization of hybrid programs in the PVS theorem prover. Plaidypvs enables the verification of the general simplex RTA framework and then, by specializing some components of the hybrid program, verifying instances of the framework while treating the untrusted component as a black box. A selection of Unmanned Aircraft Systems (UAS) operations are shown as instances of the general RTA framework in PVS. This offers the benefit of design time verification of relevant safety properties to the system, and it also gives requirements on the sample rate of sensors that determine the time interval in which the ‘switch’ property of the RTA framework is checked.

PVS↗

NASA NTRS2011

Effects of Hygrothermal Cycling on the Chemical, Thermal, and Mechanical Properties of 862/W Epoxy Resin

The hygrothermal aging characteristics of an epoxy resin were characterized over 1 year, which included 908 temperature and humidity cycles. The epoxy resin quickly showed evidence of aging through color change and increased brittleness. The influence of aging on the material s glass transition temperature (Tg) was evaluated by Differential Scanning Calorimetry (DSC) and Dynamic Mechanical Analysis (DMA). The Tg remained relatively constant throughout the year long cyclic aging profile. The chemical composition was monitored by Fourier Transform Infrared Spectroscopy (FTIR) where evidence of chemical aging and advancement of cure was noted. The tensile strength of the resin was tested as it aged. This property was severely affected by the aging process in the form of reduced ductility and embrittlement. Detailed chemical evaluation suggests many aging mechanisms are taking place during exposure to hygrothermal conditions. This paper details the influence of processes such as: advancement of cure, chemical degradation, and physical aging on the chemical and physical properties of the epoxy resin.

Miller, Sandi G.↗

NASA NTRS1988

Charged-particle transport in one dimension

A numerical solution to high-energy charged-particle transport is found by evaluation of the integral equations obtained by inverting the Boltzmann differential operator. An algorithm has been written for either continuous or discrete spectra at the boundary, which allows efficient and accurate evaluation of this rather complicated problem. Present results are compared with analytic solutions based on the perturbation theory.

NASA Discipline Radiation Health↗

NASA NTRS2022

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗

NASA NTRSDate not supplied

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

trajectory prediction↗

NASA NTRS1987

Selection principles and pattern formation in fluid mechanics and nonlinear shell theory

Research accomplishments are summarized and publications generated under the contract are listed. The general purpose of the research was to investigate various symmetry breaking problems in fluid mechanics by the use of structure parameters and selection principles. Although all of the nonlinear problems studied involved systems of partial differential equations, many of these problems led to the study of a single nonlinear operator equation of the form F(w, lambda, gamma) = 0, (w is an element of H), (lambda is an element of R1), (gamma is an element of R1). Instead of varying only the load parameter lambda, as is often done in the study of such equations, one of the main ideas used was to vary the structure parameter gamma in such a way that stable solutions were obtained. In this way one determines detailed stability results by making use of the structure of the model equations and the known physical parameters of the problem. The approach was carried out successfully for Benard-type convection problems, Taylor-like problems for short cylinders, rotating Couette-Poiseuille channel flows, and plane Couette flows. The main focus of the research was on wave theory of vortex breakdown in a tube. A number of preliminary results for inviscid axisymmetric flows were obtained.

Sather, Duane P.↗

NASA NTRS1992

Major issues and outstanding questions

The physical, chemical, geologic, and biological issues underlying current efforts to understand Mars are introduced and summarized. Attention is given not only to the present status and past history of Mars, and also to Mars as one member of the class of terrestrial planets and as one product of the evolution from the primordial accretion disk. The issue of the origin and evolution of Mars' inventory of volatile elements is treated in detail, as are the SNC meteorites, both for what these putative Martian samples indicate about volatile distributions and because of their profound impact on assessments on bulk chemical composition and the chronology of planetary differentiation and late-stage volcanism. Overviews of the major problems arising from multidisciplinary investigations of the planet are also presented.

Pepin, Robert O.↗

NASA NTRS1996

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗