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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Cooling histories of Apollo 15 quartz-normative basalts

Controlled cooling rate studies of vitrophyre 15597 are used to calculate the physical properties of a quartz-normative magma as well as cooling histories of Apollo 15 quartz-normative basalts (QNB). Pyroxene phenocrysts of a QNB magma do not settle appreciably during the cooling of the flow unit; this behavior is consistent with the small degree of chemical differentiation that occurred between members of the suite. Two types of cooling histories for QNB are discerned. The vitrophyres experienced initial slow cooling followed by rapid cooling, and the coarse-grained microgabbros underwent slow nearly linear cooling histories. Textural evidence suggests that lunar vitrophyres crystallized near the top of the flow. It is suggested that the pattern of cooling may have been caused by the chilling of liquids that inherited early crystals in the interior of the flow.

Grove, T. L.↗

NASA NTRS1984

Lubrication of Machine Elements

The understanding of hydrodynamic lubrication began with the classical experiments of Tower and Petrov. Reynolds used a reduced form of the Navier-Stokes equations and the continuity equation to generate a second order differential equation for the pressure in the narrow, converging gap of a bearing contact. Such a pressure enables a load to be transmitted between the surfaces with very low friction since the surfaces are completely separated by a film of fluid. In such a situation it is the physical properties of the lubricant, notably the dynamic viscosity, that dictate the behavior of the contact. The understanding of boundary lubrication is normally attributed to Hardy and Doubleday. In boundary lubrication it is the physical and chemical properties of thin films of molecular proportions and the surfaces to which they are attached that determine contact behavior. The lubricant viscosity is not an influential parameter. Research is devoted to a better understanding and more precise definition of other lubrication regimes between these extremes. One such regime, elastohydrodynamic lubrication, occurs in nonconformal contacts, where the pressures are high and the bearing surfaces deform elastically. In this situation the viscosity of the lubricant may raise considerably, and this further assists the formation of an effective fluid film. The science of these three lubrication regimes (hydrodynamic, elastohydrodynamic, and boundary) are described and the manner in which this science is used in the design of machine elements is examined.

Hamrock, Bernard J.↗

NASA NTRS1991

Exotic containers for capillary surfaces

This paper discusses 'exotic' rotationally symmetric containers that admit an entire continuum of distinct equilibrium capillary free surfaces. The paper extends earlier work to a larger class of parameters and clarifies and simplifies the governing differential equations, while expressing them in a parametric form appropriate for numerical integration. A unified presentation suitable for both zero and nonzero gravity is given. Solutions for the container shapes are depicted graphically along with members of the free-surface continuum, and comments are given concerning possible physical experiments.

Concus, Paul↗

NASA NTRS1993

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

NASA NTRS1993

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗

NASA NTRS1974

Efficient stream distributions in radiative transfer theory

This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.

Whitney, C.↗

NASA NTRS1979

Axial-flow compressor turning angle and loss by inviscid-viscous interaction blade-to-blade computation

A method for computation of the flow field around an arbitrary airfoil cascade on an axially symmetric blade-to-blade surface was developed which takes into account the development and separation of the blade surface boundary layers and mixing in the wake. The method predicts the overall fluid turning and total pressure loss in the context of an inviscid-viscous interaction scheme. The inviscid flow solution is obtained from a compressible flow matrix method. The viscous flow is obtained from a differential boundary layer method which calculates laminar, transitional and turbulent boundary layers. Provisions for the calculation of laminar and turbulent separation regions were added to the viscous scheme. The combined inviscid-viscous interaction scheme described yields results which are quantitatively consistent with experimental data. This suggests that the physical basis for the interactive system is correct and justifies continued exploration and use of the method.

Hansen, E. C.↗

NASA NTRS1991

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations: High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

NASA NTRS1991

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations - High Reynolds number flows

Goodrich, John W.↗

NASA NTRS1992

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗

NASA NTRS1975

Special analysis of community annoyance with aircraft noise reported by residents in the vicinity of JFK Airport, 1972

During the summer of 1972, about 1500 residents were interviewed twice in 11 communities near JFK airport. Detailed aircraft operations reports were also collected for this period, and an effort has been made to analyze recorded human response data in relation to a number of physical exposure parameters. A series of exposure indexes, based on an arithmetic integration of aircraft operations, were correlated with summated aircraft noise annoyance responses. None of these correlations were as good as the CNR index which assumes a logrithmetic integration of numbers of aircraft exposures and includes a day-night differential weighting of 10:1. There were substantial variations in average annoyance responses among communities with similar CNR exposures, substantiating previous findings that attitudinal and other personal variables also play an important role in determining annoyance differences.

Borsky, P. N.↗

NASA NTRS1993

Modeling bristle lift-off in idealized brush seal configurations

We attempt in this paper to develop a model for the flow through brush seals and determine their elastic behavior in order to predict the dependence of brush/journal clearance on geometry and operating conditions. Several idealizations regarding brush seal configuration, flow conditions, and elastic behavior are made in the analysis in order to determine closed form parametric dependence. This formulation assumes that there is no initial interference between the bristle tip and the rotor. Also, interbristle, bristle-backing plate, and bristle-rotor friction is neglected. The bristle bundle or the brush seal as it is alternately called is assumed homogeneous and isotropic on a macroscopic scale so that a physical property like permeability is uniform. The fluid is assumed to be homogeneous, incompressible, viscous, and flowing under steady conditions. A schematic of a brush seal is shown. If the nominal bristle-shaft interference is absent then under static conditions the bristles may deflect axially due to the imposed pressure differential. This axial deflection may create a clearance permitting leakage flow in excess of that which occurs through the porous matrix formed by the bristle bundles. Under dynamic conditions the Couette flow created by shaft motion could be strong enough to cause bristle deflection and once again a clearance may develop. The paper proposes a means to determine this clearance (or at least describe its parametric dependence on geometry and operating conditions) under static as well as dynamic conditions.

Modi, Vijay↗

NASA NTRS1998

A Novel Sensor for Attitude Determination Using Global Positioning System Signals

An entirely new sensor approach for attitude determination using Global Positioning System (GPS) signals is developed. The concept involves the use of multiple GPS antenna elements arrayed on a single sensor head to provide maximum GPS space vehicle availability. A number of sensor element configurations are discussed. In addition to the navigation function, the array is used to find which GPS space vehicles are within the field-of-view of each antenna element. Attitude determination is performed by considering the sightline vectors of the found GPS space vehicles together with the fixed boresight vectors of the individual antenna elements. This approach has clear advantages over the standard differential carrier-phase approach. First, errors induced by multipath effects can be significantly reduced or eliminated altogether. Also, integer ambiguity resolution is not required, nor do line biases need to be determined through costly and cumbersome self-surveys. Furthermore, the new sensor does not require individual antennas to be physically separated to form interferometric baselines to determine attitude. Finally, development potential of the new sensor is limited only by antenna and receiver technology development unlike the physical limitations of the current interferometric attitude determination scheme. Simulation results indicate that accuracies of about 1 degree (3 omega) are possible.

Crassidis, John L.↗

NASA NTRS2003

Development of Algorithms for Control of Humidity in Plant Growth Chambers

Algorithms were developed to control humidity in plant growth chambers used for research on bioregenerative life support at Kennedy Space Center. The algorithms used the computed water vapor pressure (based on measured air temperature and relative humidity) as the process variable, with time-proportioned outputs to operate the humidifier and de-humidifier. Algorithms were based upon proportional-integral-differential (PID) and Fuzzy Logic schemes and were implemented using I/O Control software (OPTO-22) to define and download the control logic to an autonomous programmable logic controller (PLC, ultimate ethernet brain and assorted input-output modules, OPTO-22), which performed the monitoring and control logic processing, as well the physical control of the devices that effected the targeted environment in the chamber. During limited testing, the PLC's successfully implemented the intended control schemes and attained a control resolution for humidity of less than 1%. The algorithms have potential to be used not only with autonomous PLC's but could also be implemented within network-based supervisory control programs. This report documents unique control features that were implemented within the OPTO-22 framework and makes recommendations regarding future uses of the hardware and software for biological research by NASA.

Costello, Thomas A.↗

NASA NTRS2010

Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple

The development of more efficient thermoelectric couple technology capable of operating with high-grade heat sources up to 1,275 K is key to improving the performance of radioisotope thermoelectric generators. Lanthanum telluride La3-xTe4 and 14-1-11 Zintls (Yb14MnSb11) have been identified as very promising materials. The fabrication of advanced high-temperature thermoelectric couples requires the joining of several dissimilar materials, typically including a number of diffusion bonding and brazing steps, to achieve a device capable of operating at elevated temperatures across a large temperature differential (up to 900 K). A thermoelectric couple typically comprises a heat collector/ exchanger, metallic interconnects on both hot and cold sides, n-type and ptype conductivity thermoelectric elements, and cold-side hardware to connect to the cold-side heat rejection and provide electrical connections. Differences in the physical, mechanical, and chemical properties of the materials that make up the thermoelectric couple, especially differences in the coefficients of thermal expansion (CTE), result in undesirable interfacial stresses that can lead to mechanical failure of the device. The problem is further complicated by the fact that the thermoelectric materials under consideration have large CTE values, are brittle, and cracks can propagate through them with minimal resistance. The inherent challenge of bonding brittle, high-thermal-expansion thermoelectric materials to a hot shoe material that is thick enough to carry the requisite electrical current was overcome. A critical advantage over prior art is that this device was constructed using all diffusion bonds and a minimum number of assembly steps. The fabrication process and the materials used are described in the following steps: (1) Applying a thin refractory metal foil to both sides of lanthanum telluride. To fabricate the n-type leg of the advanced thermoelectric couple, the pre-synthesized lanthanum telluride coupon was diffusion bonded to the metal foil using a thin adhesion layer. (2) Repeating a similar process for the 14-1-11 Zintl p-type leg of the advanced thermoelectric couple. (3) Bonding thick CTE-matched metal plates on the metallized lanthanum telluride and Yb14MnSb11 to form the hot and cold sides of the thermoelectric couple. The calculated conversion efficiency of such an advanced couple would be about 10.5 percent, about 35 percent better than heritage radioisotope thermoelectric technology that relies on Si-Ge alloys. In addition, unlike Si-Ge alloys, these materials can be combined with many other thermoelectric materials optimized for operation at lower temperatures to achieve conversion efficiency in excess of 15 percent (a factor of 2 increase over heritage technology).

Ravi, Vilupanur A.↗

NASA NTRS2009

Optimization of Layer Densities for Spacecraft Multilayered Insulation Systems

Numerous tests of various multilayer insulation systems have indicated that there are optimal densities for these systems. However, the only method of calculating this optimal density was by a complex physics based algorithm developed by McIntosh. In the 1970's much data were collected on the performance of these insulation systems with many different variables analyzed. All formulas generated included number of layers and layer density as geometric variables in solving for the heat flux, none of them was in a differentiable form for a single geometric variable. It was recently discovered that by converting the equations from heat flux to thermal conductivity using Fourier's Law, the equations became functions of layer density, temperatures, and material properties only. The thickness and number of layers of the blanket were merged into a layer density. These equations were then differentiated with respect to layer density. By setting the first derivative equal to zero, and solving for the layer density, the critical layer density was determined. Taking a second derivative showed that the critical layer density is a minimum in the function and thus the optimum density for minimal heat leak, this is confirmed by plotting the original function. This method was checked and validated using test data from the Multipurpose Hydrogen Testbed which was designed using McIntosh's algorithm.

Johnson, W. L.↗

NASA NTRS1998

Mars Volatiles and Climate Surveyor (MVACS) Integrated Payload for the Mars Polar Lander Mission

The Mars Volatiles and Climate Surveyor (MVACS) integrated payload for the Mars Polar Lander will be launched in January 1999, with a scheduled landing on Mars' south-polar layered deposits in December 1999. Over the course of its 90-day nominal mission during the martian southern spring and summer seasons, it will make in situ measurements that will provide new insights into the behavior and distribution of martian volatiles, MVACS consists of four major instrument systems: a surface stereo imager (SSI), which will acquire multispectral stereo images of the surface and atmosphere; a 2-m robotic arm (RA), which will dig a O.5-m deep trench and acquire surface and subsurface samples that will be imaged by a focusable robotic arm camera (RAC), which will take close-up images of surface and subsurface samples at a spatial resolution of 21 micron; a meteorology package (MET), which will make the first measurements of surface pressure, temperature, and winds in Mars' southern hemisphere and employ a tunable diode laser (TDL) spectrometer to measure the water-vapor concentration and isotopic composition of CO2 in the martian atmosphere; and a thermal and evolved gas analyzer (TEGA), which will use differential scanning calorimetry and TDL-evolved gas analysis to determine the concentrations of ices, adsorbed volatiles, and volatile-bearing minerals in surface and subsurface soil samples. The unique in situ measurements made by MVACS at its high-latitude landing site will define a number of important aspects of the physical, isotopic, and chemical nature of the martian near-surface and subsurface environment that will be valuable in better understanding Mars meteorites and returned samples, as well as in the search for martian resources that could be utilized by humans.

Paige, D. A.↗

NASA NTRS1998

Differential Absorption Lidar (DIAL) Measurements of Atmospheric Water Vapor Utilizing Robotic Aircraft

A new unpiloted air vehicle (UAV) based water vapor DIAL system will be described. This system is expected to offer lower operating costs, longer test duration and severe weather capabilities. A new high-efficiency, compact, light weight, diode-pumped, tunable Cr:LiSAF laser will be developed to meet the UAV payload weight and size limitations and its constraints in cooling capacity, physical size and payload. Similarly, a new receiver system using a single mirror telescope and an avalanche photo diode (APD) will be developed. Projected UAV parameters are expected to allow operation at altitudes up to 20 km, endurance of 24 hrs and speed of 400 km/hr. At these conditions measurements of water vapor at an uncertainty of 2-10% with a vertical resolution of 200 m and horizontal resolution of 10 km will be possible.

Hoang, Ngoc↗