SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

NASA NTRS1989

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

NASA NTRS2000

Combustion Response of AP Composite Propellants

The Cohen & Strand model for AP composite propellants is applied as boundary conditions, one for AP and one for binder, in solving the heat conduction equation in each to compute linear and non-linear combustion response properties for each and for the aggregate propellant. Iterations couple AP and binder through the quasi-steady flame processes. Illustrative results for linear response functions (pressure-coupled and velocity-coupled) are presented for a monomodal AP propellant varying AP size, pressure and crossflow speed, and examples of non-linear responses to arbitrary waveforms are shown. A quantitative comparison with response function data is limited to one well-characterized research formulation. Mechanistic implications are discussed, including recommendations for future modeling work.

Shusser, Michael↗

NASA NTRSDate not supplied

Towards Energy Scale Calibration and Drift Correction of TES Detectors for Athena X-IFU

The Athena X-Ray Integral Field Unit (X-IFU) comprises a 2376-pixel array of transition edge sensors (TES) read out with time-division multiplexing (TDM). X-IFU will provide spatially resolved, high-resolution spectroscopy (2.5 eV full-width-half-maximum up to 7 keV) over the energy range 0.2 to 12 keV, with an absolute energy scale accuracy of 0.4 eV. The energy scale function maps the optimally filtered pulse height, in arbitrary engineering units, to real calibrated energy. Uncertainties in the calibration can result from imperfect fitting of the energy scale between the known calibration points. Furthermore, temporal changes in the TES operating environment, such as heat-sink temperature, magnetic field and bias voltage, can cause significant variations in the detector gain function over time. If not properly corrected, this can result in degradation of the energy resolution, and systematic errors in the absolute energy scale. The non-linear nature of TES detectors, coupled with the possibility of multiple simultaneously occurring sources of drift, can make effective corrections over the full bandpass of the instrument extremely challenging. Athena X-IFU will employ an on-board calibration source that provides known reference x-ray lines. This provides real-time monitoring of the gain stability of the detector system and information that can be used to correct for gain drifts. For X-IFU the baseline approach is to measure a series of calibration curves under different environmental conditions, which bound the expected drifts the instrument is predicted to see over the course of the mission. Using the information from the in-flight calibration source, these energy scale functions can be interpolated to generate a new corrected energy scale as a function of time. In this paper we discuss progress towards demonstrating that the X-IFU energy scale requirements can be met. We present measurements on ~ 200 pixels in a prototype X-IFU array read out with 8-column x 32-row TDM. We use a rotating target source containing 12 fluorescent targets to generate x-ray lines covering the energy range 4 keV (Sc-Kα) to 12 keV (Br-Kα). We present measurements of the non-linear energy scale function and show how variations in heat-sink temperature, TES bias voltage and magnetic field affect the shape of TES energy scale differently and introduce different residual gain errors over the bandpass. We explore different drift correction algorithms that use either a single or multiple referential lines to track and correct the gain from these various sources of drift. In addition to the pulse-height, the DC ‘baseline’ level of the TES can contain information about its bias conditions. Thus, we test a multi-parameter gain correction algorithm that attempts to incorporate both the pulse height and the additional baseline information into the algorithm.

Stephen J Smith↗

NASA NTRS1984

Generalized serial search code acquisition - The equivalent circular state diagram approach

A transform-domain method for deriving the generating function of the acquisition process resulting from an arbitrary serial search strategy is presented. The method relies on equivalent circular state diagrams, uses Mason's formula from flow-graph theory, and employs a minimum number of required parameters. The transform-domain approach is briefly described and the concept of equivalent circular state diagrams is introduced and exploited to derive the generating function and resulting mean acquisition time for three particular cases of interest, the continuous/center Z search, the broken/center Z search, and the expanding window search. An optimization of the latter technique is performed whereby the number of partial windows which minimizes the mean acquisition time is determined. The numerical results satisfy certain intuitive predictions and provide useful design guidelines for such systems.

Polydoros, A.↗

NASA NTRS1969

On the bound of first excursion probability

Method has been developed to improve the lower bound of the first excursion probability that can apply to the problem with either constant or time-dependent barriers. The method requires knowledge of the joint density function of the random process at two arbitrary instants.

Yang, J. N.↗

NASA NTRS1972

Computer program draws three-dimensional surfaces

Computer plotting program PLOT 3D draws views of surface forms z = f(x,y). Surface thus defined by program may be drawn after arbitrary rotations. Program portrays behavior of various functions involving two variables in many engineering, physics, and mathematical relationships.

Canright, R. B., Jr.↗

NASA NTRS1971

Calculation of compressible nonadiabatic boundary layers in laminar, transitional and turbulent flow by the method of integral relations

A computer program was developed to do the calculations for two-dimensional or axisymmetric configurations from low speeds to hypersonic speeds with arbitrary streamwise pressure, temperature, and Mach number distributions. Options are provided for obtaining initial conditions either from experimental information or from a theoretical similarity solution. The transition region can be described either by an arbitrary distribution of intermittency or by a function based on Emmons' probability theory. Correlations were developed for use in estimating the parameters of the theoretical intermittency function. Correlations obtained from other sources are used for estimating the transition point. Comparisons were made between calculated and measured boundary layer quantities for laminar, transitional, and turbulent flows on flat plates, cones, cone flares, and a waisted body of revolution. Excellent agreement was obtained between the present theory and two other theories based on the method of finite differences. The intermittency required to reproduce some experimental heat transfer results in hypersonic flow was found to be quite different from the theoretical function. It is suggested that the simple probability theory of Emmons may not be valid for representing the intermittency of hypersonic transitional boundary layers and that the program could be useful as a tool for detailed study of the intermittency of the transition region.

Kuhn, G. D.↗

NASA NTRS1971

Shell splitting due to electric fields.

Drift shells of trapped particles starting from a single field line in the earth's magnetic field may become separated owing to an electric field. By using the formalism of the shell-splitting function, this separation is calculated for an arbitrary curl-free time-independent electric field satisfying E.B = 0, with no limitation on the particle's longitudinal excursion. The effect of nondipole field components is also taken into account to the lowest order.

Stern, D. P.↗

NASA NTRS1972

Pattern Recognition and Linear Discriminant Function

The properties of a linear discriminant function are discussed for the case of arbitrary distributions with equal covariance matrices. Using two examples, a comparison is made showing how the difference of the means relates to the covariance matrices.

Finch, L. H.↗

NASA NTRS1976

The prediction of zenith range refraction from surface measurements of meteorological parameters

In the last two decades, increasingly sophisticated deep space missions have placed correspondingly stringent requirements on navigational accuracy. As part of the effort to increase navigational accuracy, and hence the quality of radiometric data, much effort has been expended in an attempt to understand and compute the tropospheric effect on range (and hence range rate) data. The general approach adopted has been that of computing a zenith range refraction, and then mapping this refraction to any arbitrary elevation angle via an empirically derived function of elevation. The prediction of zenith range refraction derived from surface measurements of meteorological parameters is presented. Refractivity is separated into wet (water vapor pressure) and dry (atmospheric pressure) components. The integration of dry refractivity is shown to be exact. Attempts to integrate wet refractivity directly prove ineffective; however, several empirical models developed by the author and other researchers at JPL are discussed. The best current wet refraction model is here considered to be a separate day/night model, which is proportional to surface water vapor pressure and inversely proportional to surface temperature. Methods are suggested that might improve the accuracy of the wet range refraction model.

Berman, A. L.↗

NASA NTRS1987

Building validation tools for knowledge-based systems

The Expert Systems Validation Associate (EVA), a validation system under development at the Lockheed Artificial Intelligence Center for more than a year, provides a wide range of validation tools to check the correctness, consistency and completeness of a knowledge-based system. A declarative meta-language (higher-order language), is used to create a generic version of EVA to validate applications written in arbitrary expert system shells. The architecture and functionality of EVA are presented. The functionality includes Structure Check, Logic Check, Extended Structure Check (using semantic information), Extended Logic Check, Semantic Check, Omission Check, Rule Refinement, Control Check, Test Case Generation, Error Localization, and Behavior Verification.

Stachowitz, R. A.↗

NASA NTRS1999

Remote Sensing of Rain

The first problem addressed concerns passive-microwave rain retrievals. Most current approaches start by building off-line a cloud-model-derived database. Given data, the retrieval algorithms search the database for the microwave temperatures "closest" to the observed data, then after some fine-tuning (performed in different ways by different implementations) the rain is estimated to be that which corresponds to the selected (and fine-tuned) set of database temperatures. These approaches have three drawbacks: they cannot properly take into account the ambiguities which arise from the fact that several rain scenarios can produce the same observed temperatures; they are quite inefficient since they require manipulating a large database along with often complex "fine-tuning" procedures; and they cannot refine their estimates if additional data is available. This past year we have derived closed formulae relating observed microwave brightness temperatures, T(sub b), and the underlying rain rates, R: average T(sub b) =f (rain) and average rain = g (T(sub b)), along with the corresponding covariance matrices. These results are sufficient to describe the conditional probabilities p(R/T(sub b)) and p(T(sub b)/R) to second order. Progress has also been made towards deriving a robust description of the rain drop size distribution (DSD). The widespread approach consisting in parameterizing the DSD as a gamma-distribution in terms of the drop diameter D suffers from the facts that, in reality, the DSD is not a smooth function of D and that the largely arbitrary Gamma model imposes unintended behavior, which has implications on any quantities derived from the DSD model. We have therefore developed a non-parametric yet practical description of the DSD, which is particularly well-suited for use in remote-sensing applications. The diagram on the left shows a comparison between an actual DSD sample and the truncated non-parametric representation. One figure shows the relation between radar reflectivity and rain rate derived using this representation. Validation of the Tropical Rainfall Measuring Mission (TRMM) radar-radiometer combined R and DSD algorithm is underway. This algorithm was designed to make optimal use of the instantaneous reflectivity profiles measured by the TRMM radar and the microwave brightness temperatures measured by the TRMM passive radiometer. So far, it appears to be the most reliable TRMM rain algorithm.

Haddad, Ziad S.↗

NASA NTRS2003

A Model for the Oxidation of C/SiC Composite Structures

A mathematical theory and an accompanying numerical scheme have been developed for predicting the oxidation behavior of C/SiC composite structures. The theory is derived from the mechanics of the flow of ideal gases through a porous solid. Within the mathematical formulation, two diffusion mechanisms are possible: (1) the relative diffusion of one species with respect to the mixture, which is concentration gradient driven and (2) the diffusion associated with the average velocity of the gas mixture, which is total gas pressure gradient driven. The result of the theoretical formulation is a set of two coupled nonlinear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations must be solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of space and time. The local rate of carbon oxidation is determined as a function of space and time using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The nonlinear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual method, allowing for the solution of the differential equations numerically. The end result is a numerical scheme capable of determining the variation of the local carbon oxidation rates as a function of space and time for any arbitrary C/SiC composite structures.

Sullivan, Roy M.↗

NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗

NASA NTRS1994

The effects of gravity darkening on the ultraviolet continuum polarization produced by circumstellar disks

We investigate the effects of gravity darkening on the UV continuum polarization produced by an axisymmetric disk that surrounds a rapidly rotating star. Although the model is a single scattering approximation, we do include the effects of attenuation (electron scattering plus hydrogen bound-free absorption) by the disk, using an approach similar to that of Sobolev (1963). Because of the gravity darkening of the star and the attenuation within the disk, the radiation field is not axially symmetric about the radius vector. This implies that the polarization source functions are no longer provided by the finite disk depolarization factors of Cassinelli, Nordsieck, & Murison (1987), which are functions of the intensity moments in a spherically symmetric atmosphere. We reformulate the polarization source functions using generalized intensity moment tensors (J, H(sub i), K(sub ij)) that are valid for an arbitrary radiation field and envelope geometry. We find that the polarization source functions are simplest when using intensity moments in the observer's reference frame. On the other hand, the intensity moments are most easily evaluated in the stellar reference frame. Using the rotation transformation properties of the generalized intensity moments, we relate the observer's moments to those evaluated in the stellar reference frame. Our procedure for determining the polarization source functions thus merely involves choosing a set of Euler angles for the coordinate rotations, and then evaluating the associated rotation matrix. The geometrical complications of polarization transfer are thus reduced to obtaining a coordinate rotation matrix.

Bjorkman, J. E.↗

NASA NTRS1980

A modified Barrett-Lampard expansion and its application to bandpass nonlinearities with both AM-AM and AM-PM conversion

A two-dimensional Barrett-Lampard expansion is developed to handle two-dimensional nonlinear systems with random Gaussian inputs. The theory is applied to the problem of evaluating the mean signal and the inphase and quadrature noise correlation functions at the output of a TWT nonlinearity exhibiting arbitrary AM-AM and AM-PM characteristics, when the input consists of a narrow-band signal and Gaussian noise. As a specialized application, the classical problem of determining the coherent and incoherent noise correlation functions at the output of a power-law bandpass nonlinearity is also considered. The results are of interest in assessing the performance of coherent satellite communication channels. The theory developed is accompanied by numerical examples of practical interest.

Chie, C. M.↗

NASA NTRS1985

Response functions of free mass gravitational wave antennas

The work of Gursel, Linsay, Spero, Saulson, Whitcomb and Weiss (1984) on the response of a free-mass interferometric antenna is extended. Starting from first principles, the earlier work derived the response of a 2-arm gravitational wave antenna to plane polarized gravitational waves. Equivalent formulas (generalized slightly to allow for arbitrary elliptical polarization) are obtained by a simple differencing of the '3-pulse' Doppler response functions of two 1-arm antennas. A '4-pulse' response function is found, with quite complicated angular dependences for arbitrary incident polarization. The differencing method can as readily be used to write exact response functions ('3n+1 pulse') for antennas having multiple passes or more arms.

Estabrook, F. B.↗