Engineering Papers⌕ Search

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.

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

At least 235 records · Page 13

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗

Some examples of image warping for low vision prosthesis

NASA has developed an image processor, the Programmable Remapper, for certain functions in machine vision. The Remapper performs a highly arbitrary geometric warping of an image at video rate. It might ultimately be shrunk to a size and cost that could allow its use in a low-vision prosthesis. Coordinate warpings have been developed for retinitis pigmentosa (tunnel vision) and for maculapathy (loss of central field) that are intended to make best use of the patient's remaining viable retina. The rationales and mathematics are presented for some warpings that we will try in clinical studies using the Remapper's prototype.

Juday, Richard D.↗

Effect of fiber orientation on the fracture toughness of brittle matrix composites

The effective fracture toughness of brittle matrix materials can be increased through the addition of short, poorly bonded fibers which bridge the growing crack. The orientation distribution of the fibers is likely to be biased, and not in an ideal random or aligned state. A micromechanical model is formulated for the postcracking behavior using the force-displacement relation for an arbitrary fiber bridging a crack, the fiber orientation density function, and the fiber location density function. This model is then used to determine an effective traction law for the bridging fibers, as well as the steady state bridging toughness increment. In most cases, the results may be placed in the form of a product of the aligned fiber results times a modifying integrated orientation factor. The frictional shear stress on fiber pull-out is allowed to vary during pull-out, modeling the effects of matrix breakdown, fiber surface smoothing or wear debris accumulation. Results are presented for a variety of representative planar and three-dimensional fiber orientation states.

Jain, L. K.↗

Development and Operation of a Material Identification and Discrimination Imaging Spectroradiometer

Many imaging applications require quantitative determination of a scene's spectral radiance. This paper describes a new system capable of real-time spectroradiometric imagery. Operating at a full-spectrum update rate of 30Hz, this imager is capable of collecting a 30 point spectrum from each of three imaging heads: the first operates from 400 nm to 950 nm, with a 2% bandwidth; the second operates from 1.5 micro-m to 5.5 micro-m with a 1.5% bandwidth; the third operates from 5 micro-m to 12 micro-m, also at a 1.5% bandwidth. Standard image format is 256 x 256, with 512 x 512 possible in the VIS/NIR head. Spectra of up to 256 points are available at proportionately lower frame rates. In order to make such a tremendous amount of data more manageable, internal processing electronics perform four important operations on the spectral imagery data in real-time. First, all data in the spatial/spectral cube of data is spectro-radiometrically calibrated as it is collected. Second, to allow the imager to simulate sensors with arbitrary spectral response, any set of three spectral response functions may be loaded into the imager including delta functions to allow single wavelength viewing; the instrument then evaluates the integral of the product of the scene spectral radiances and the response function. Third, more powerful exploitation of the gathered spectral radiances can be effected by application of various spectral-matched filtering algorithms to identify pixels whose relative spectral radiance distribution matches a sought-after spectral radiance distribution, allowing materials-based identification and discrimination. Fourth, the instrument allows determination of spectral reflectance, surface temperature, and spectral emissivity, also in real-time. The spectral imaging technique used in the instrument allows tailoring of the frame rate and/or the spectral bandwidth to suit the scene radiance levels, i.e., frame rate can be reduced, or bandwidth increased to improve SNR when viewing low radiance scenes. The unique challenges of design and calibration are described. Pixel readout rates of 160 MHz, for full frame readout rates of 1000 Hz (512 x 512 image) present the first challenge; processing rates of nearly 600 million integer operations per second for sensor emulation, or over 2 billion per second for matched filtering, present the second. Spatial and spectral calibration of 66,536 pixels (262,144 for the 512 x 512 version) and up to 1,000 spectral positions mandate novel decoupling methods to keep the required calibration memory to a reasonable size. Large radiometric dynamic range also requires care to maintain precision operation with minimum memory size.

Dombrowski, Mark↗

Application of Normal Mode Expansion to AE Waves in Finite Plates

Breckenridge et al. (1975), Hsu (1985) and Pao (1978) adapted approaches from seismology to calculate the response at the surface of an infinite half-space and an infinite plate. These approaches have found use in calibrating acoustic emission (AE) transducers. However, it is difficult to extend this theoretical approach to AE testing of practical structures. Weaver and Pao (1982) considered a normal mode solution to the Lamb equations. Hutchinson (1983) pointed out the potential relevance of Mindlin's plate theory (1951) to AE. Pao (1982) reviewed Medick s (1961) classical plate theory for a point source, but rejected it as useful for AE and no one seems to have investigated its relevance to AE any further. Herein, a normal mode solution to the classical plate bending equation was investigated for its applicability to AE. The same source-time function chosen by Weaver and Pao is considered. However, arbitrary source and receiver positions are chosen relative to the boundaries of the plate. This is another advantage of the plate theory treatment in addition to its simplicity. The source does not have to be at the center of the plate as in the axisymmetric treatment. The plate is allowed to remain finite and reflections are predicted. The importance of this theory to AE is that it can handle finite plates, realistic boundary conditions, and can be extended to composite materials.

Gorman, M. R.↗

Cyclotron Resonant Scattering Feature Simulations II. Description of the CRSF Simulation Process

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons o_ quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems.Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time.Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua.Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online (see link in footnote, page 1).

Schwarm, F.-W.↗

Optimal estimation of large structure model errors

In-flight estimation of large structure model errors is usually required as a means of detecting inevitable deficiencies in large structure controller/estimator models. The present paper deals with a least-squares formulation which seeks to minimize a quadratic functional of the model errors. The properties of these error estimates are analyzed. It is shown that an arbitrary model error can be decomposed as the sum of two components that are orthogonal in a suitably defined function space. Relations between true and estimated errors are defined. The estimates are found to be approximations that retain many of the significant dynamics of the true model errors. Current efforts are directed toward application of the analytical results to a reference large structure model.

Rodriguez, G.↗

Kinematics, influence functions and field quantities for disturbance propagation from moving disturbance sources

A unified method is presented for deriving the influence functions of moving singularities which determine the field quantities in aerodynamics and aeroacoustics. The moving singularities comprise volume and surface distributions having arbitrary orientations in space and to the trajectory. Hence one generally valid formula for the influence functions which reveal some universal relationships and remarkable properties in the disturbance fields. The derivations used are completely consistent with the physical processes in the propagation field, such that treatment renders new descriptions for some standard concepts. The treatment is uniformly valid for subsonic and supersonic Mach numbers.

Das, A.↗

Autoreturn Function for a Remotely Piloted Vehicle

An algorithm to maneuver an air vehicle to intercept and follow a pre-planned path while remaining within an arbitrary, closed boundary is outlined. The immediate application is for an autonomous lost-link return-to-runway function for a remotely piloted vehicle being developed by NASA, but other applications are hypothesized. Results of implementation in a flight simulator are given.

McMinn, J. D.↗

Derivation of a Tappered p-Version Beam Finite Element

A tapered p-version beam finite element suitable for dynamic applications is derived. The taper in the element is represented by allowing the area moments of inertia to vary as quartic polynomials along the length of the beam, and the cross-sectional area to vary as a quadratic polynomial. The p-version finite-element characteristics are implemented through a set of polynomial shape functions. The lower-order shape functions are identical to the classical cubic and linear shape functions normally associated with a beam element. The higher-order shape functions are a hierarchical set of polynomials that are integrals of orthogonal polynomials. Explicit expressions for the mass and stiffness matrices are presented for an arbitrary value of p. The element has been verified to be numerically stable using shape functions through 22nd order.

Hinnant, Howard E.↗

Estimation of Optimum Stimulus Amplitude for Balance Training using Electrical Stimulation of the Vestibular System

Sensorimotor changes such as posture and gait instabilities can affect the functional performance of astronauts after gravitational transitions. Sensorimotor Adaptability (SA) training can help alleviate decrements on exposure to novel sensorimotor environments based on the concept of 'learning to learn' by exposure to varying sensory challenges during posture and locomotion tasks (Bloomberg 2015). Supra-threshold Stochastic Vestibular Stimulation (SVS) can be used to provide one of many challenges by disrupting vestibular inputs. In this scenario, the central nervous system can be trained to utilize veridical information from other sensory inputs, such as vision and somatosensory inputs, for posture and locomotion control. The minimum amplitude of SVS to simulate the effect of deterioration in vestibular inputs for preflight training or for evaluating vestibular contribution in functional tests in general, however, has not yet been identified. Few studies (MacDougall 2006; Dilda 2014) have used arbitrary but fixed maximum current amplitudes from 3 to 5 mA in the medio-lateral (ML) direction to disrupt balance function in healthy adults. Giving this high level of current amplitude to all the individuals has a risk of invoking side effects such as nausea and discomfort. The goal of this study was to determine the minimum SVS level that yields an equivalently degraded balance performance. Thirteen subjects stood on a compliant foam surface with their eyes closed and were instructed to maintain a stable upright stance. Measures of stability of the head, trunk, and whole body were quantified in the ML direction. Duration of time they could stand on the foam surface was also measured. The minimum SVS dosage was defined to be that level which significantly degraded balance performance such that any further increase in stimulation level did not lead to further balance degradation. The minimum SVS level was determined by performing linear fits on the performance variable at different stimulation levels. Results from the balance task suggest that there are inter-individual differences and the minimum SVS amplitude was found to be in the range of 1 mA to 2.5 mA across subjects. SVS resulted in an average decrement of balance task performance in the range of 62%-73% across different measured variables at the minimum SVS amplitude in comparison to the control trial (no stimulus). Training using supra-threshold SVS stimulation is one of the sensory challenges used for preflight SA training designed to improve adaptability to novel gravitational environments. Inter-individual differences in response to SVS can help customize the SA training paradigms using minimal dosage required. Another application of using SVS is to simulate acute deterioration of vestibular sensory inputs in the evaluation of tests for assessing vestibular function.

Goel, R.↗

Finite difference grid generation by multivariate blending function interpolation

The General Interpolants Method (GIM) code which solves the multidimensional Navier-Stokes equations for arbitrary geometric domains is described. The geometry module in the GIM code generates two and three dimensional grids over specified flow regimes, establishes boundary condition information and computes finite difference analogs for use in the GIM code numerical solution module. The technique can be classified as an algebraic equation approach. The geometry package uses multivariate blending function interpolation of vector-values functions which define the shapes of the edges and surfaces bounding the flow domain. By employing blending functions which conform to the cardinality conditions the flow domain may be mapped onto a unit square (2-D) or unit cube (3-D), thus producing an intrinsic coordinate system for the region of interest. The intrinsic coordinate system facilitates grid spacing control to allow for optimum distribution of nodes in the flow domain.

Anderson, P. G.↗

Terminal Descent Sensor Simulation

Sulcata software simulates the operation of the Mars Science Laboratory (MSL) radar terminal descent sensor (TDS). The program models TDS radar antennas, RF hardware, and digital processing, as well as the physics of scattering from a coherent ground surface. This application is specific to this sensor and is flexible enough to handle end-to-end design validation. Sulcata is a high-fidelity simulation and is used for performance evaluation, anomaly resolution, and design validation. Within the trajectory frame, almost all internal vectors are represented in whatever coordinate system is used to represent platform position. The trajectory frame must be planet-fixed. The platform body frame is specified relative to arbitrary reference points relative to the platform (spacecraft or test vehicle). Its rotation is a function of time from the trajectory coordinate system specified via dynamics input (file for open loop, callback for closed loop). Orientation of the frame relative to the body is arbitrary, but constant over time. The TDS frame must have a constant rotation and translation from the platform body frame specified at run time. The DEM frame has an arbitrary, but time-constant, rotation and translation with respect to the simulation frame specified at run time. It has the same orientation as sigma0 frame, but is possibly translated. Surface sigma0 has the same arbitrary rotation and translation as DEM frame.

Chen, Curtis W.↗

Subsonic potential aerodynamics for complex configurations - A general theory

A general theory of subsonic potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the velocity potential is obtained for both supersonic and subsonic flow. Under the small perturbation assumption, the potential at any point in the field depends only upon the values of the potential and its normal derivative on the surface of the body. On the surface of the body, this representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The theory is applied to finite-thickness wings in subsonic steady and oscillatory flows.

Morino, L.↗

Analysis of acoustic radiation in a jet flow environment

The analysis of sound fields from arbitrary source distributions in terms of Legendre and spherical Hankel functions is well known. The purpose of this paper is to extend this classical method of analysis to environments such as jet flows where flow and flow gradients are inherently present. The wave-equation governing the radiation of sound in such an environment is derived. The steady state flow and flow gradients in the axial and transverse directions appear as coefficients in the terms of the wave-equation. A semi-numerical method is used to solve the wave-equation in terms of modified spherical harmonics yielding the phase velocities and the directivities of an infinite set of modes. The directivity of each mode is obtained in terms of modified Legendre functions by numerical integration. Some results of these directivity and phase-velocity calculations are presented for a limited number of frequency and flow parameters. Both convective and shear refraction are shown to be important.

Mungur, P.↗

A General Theory of Unsteady Compressible Potential Aerodynamics

The general theory of potential aerodynamic flow around a lifting body having arbitrary shape and motion is presented. By using the Green function method, an integral representation for the potential is obtained for both supersonic and subsonic flow. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface, sigma, of the body. Hence, if the point P approaches the surface of the body, the representation reduces to an integro-differential equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface sigma. For the important practical case of small harmonic oscillation around a rest position, the equation reduces to a two-dimensional Fredholm integral equation of second-type. It is shown that this equation reduces properly to the lifting surface theories as well as other classical mathematical formulas. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings.

Morino, L.↗

On the formulation of the aerodynamic characteristics in aircraft dynamics

The theory of functionals is used to reformulate the notions of aerodynamic indicial functions and superposition. Integral forms for the aerodynamic response to arbitrary motions are derived that are free of dependence on a linearity assumption. Simplifications of the integral forms lead to practicable nonlinear generalizations of the linear superpositions and stability derivative formulations. Applied to arbitrary nonplanar motions, the generalization yields a form for the aerodynamic response that can be compounded of the contributions from a limited number of well-defined characteristic motions, in principle reproducible in the wind tunnel. Further generalizations that would enable the consideration of random fluctuations and multivalued aerodynamic responses are indicated.

Tobak, M.↗

Transient thermal convection in microgravity

The unsteady two-dimensional thermal convection in a cylinder due to a transient acceleration solved for a step-function excitation. From this, the solution was obtained for an arbitrary time-dependent acceleration. The solutions are valied for sufficiently low Rayleigh numbers and therefore, relevant to microgravity fields. As an example, two graphs are presented for he convection resulting from the movement of an astronaut inside the Shuttle. The analysis can be applied to obtain any other convective flows such as those caused by g-jitter or variable rotation of the Shuttle.

Dressler, R. F.↗