Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite precision”

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

Advanced Intensity-Modulation Continuous-Wave Lidar Techniques for ASCENDS O2 Column Measurements

Global atmospheric carbon dioxide (CO2) measurements for the NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) space mission are critical for improving our understanding of global CO2 sources and sinks. Advanced Intensity- Modulated Continuous-Wave (IM-CW) lidar techniques are investigated as a means of facilitating CO2 measurements from space to meet the ASCENDS measurement requirements. In recent numerical, laboratory and flight experiments we have successfully used the Binary Phase Shift Keying (BPSK) modulation technique to uniquely discriminate surface lidar returns from intermediate aerosol and cloud contamination. We demonstrate the utility of BPSK to eliminate sidelobes in the range profile as a means of making Integrated Path Differential Absorption (IPDA) column CO2 measurements in the presence of optically thin clouds, thereby eliminating the need to correct for sidelobe bias errors caused by the clouds. Furthermore, high accuracy and precision ranging to the surface as well as to the top of intermediate cloud layers, which is a requirement for the inversion of column CO2 number density measurements to column CO2 mixing ratios, has been demonstrated using new hyperfine interpolation techniques that takes advantage of the periodicity of the modulation waveforms. This approach works well for both BPSK and linear swept-frequency modulation techniques. The BPSK technique under investigation has excellent auto-correlation properties while possessing a finite bandwidth. A comparison of BPSK and linear swept-frequency is also discussed in this paper. These results are extended to include Richardson-Lucy deconvolution techniques to extend the resolution of the lidar beyond that implied by limit of the bandwidth of the modulation, where it is shown useful for making tree canopy measurements.

Campbell, Joel F.↗

Support trusses for large precision segmented reflectors - Preliminary design and analysis

The Precision Segmented Reflector (PSR) primary structures plan is outlined. Geometries and design considerations for erectable and deployable reflector support structures are discussed. Support truss requirements and goals for the PSR are given, and the results of static and dynamic analyses of a prototype four meter diameter structure are presented. In addition, similar results are presented for two 20-meter diameter support trusses. Implications of the analyses for the PSR program are considered and the formulation and limitations of current PSR finite element models are discussed. It is shown that if the secondary optical system is supported by a simple tripod design, the first six vibration modes are likely to be dominated by the secondary system. The 20-meter diameter support trusses are found to be quite stiff for structures of such large size.

Collins, Timothy J.↗

'Backscratching' alternative to PRSI pattern recognition

The Backscratching optical correlation algorithm has been proposed for four degree of freedom tracking. In an alternating Cartesian and log-polar implementation, the tracked parameters are scale, rotation, and two-dimenisonal translation. The algorithm has a finite capture radius in the four-dimensional tracking space. The capture radius is dependent on the tracked object, the correlator architecture, and the method of filter computation. Some methods of extending the capture radius are discussed. One is a modification of matched filters, another is a careful consideration of log-polar transform center, and another is an operational method. Some simulations of the filter construction method, in which a larger capture radius is gained at the expense of precision in determining the four parameters are presented.

Monroe, Stanley E., Jr.↗

NASTRAN finite element idealization study

The investigation of the effects of variations of mesh refinement and mesh pattern were conducted using a basic rectangular mesh pattern. When employing the constant strain TRMEM element, the basic rectangular pattern was subdivided into triangles. This subdivision employs two different triangular patterns and allows results to be obtained which demonstrate the effect of modelling bias. Errors in tip deflection, direct stress, and shearing stress as a function of mesh size and element aspect ratio were obtained as well as mid-span stress distributions. All problems were solved on an IBM 360/95 computer using MacNeal-Schwendler Version MSC-38 Rigid Format-1 of the NASTRAN computer program. While NASTRAN uses double precision arithmetic for the solution of the global equations for displacements, subsequent computations to obtain element stresses are carried out in single precision. This suggests that some improvement in stress recovery might be expected when using the higher precision DCD machines.

Case, W. R.↗

Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays

Single-slope analog-to-digital converters (ADCs) are particularly useful for onchip digitization in focal plane arrays (FPAs) because of their inherent monotonicity, relative simplicity, and efficiency for column-parallel applications, but they are comparatively slow. Squareroot encoding can allow the number of code values to be reduced without loss of signal-to-noise ratio (SNR) by keeping the quantization noise just below the signal shot noise. This encoding can be implemented directly by using a quadratic ramp. The reduction in the number of code values can substantially increase the quantization speed. However, in an FPA, the fixed pattern noise (FPN) limits the use of small quantization steps at low signal levels. If the zero-point is adjusted so that the lowest column is onscale, the other columns, including those at the center of the distribution, will be pushed up the ramp where the quantization noise is higher. Additionally, the finite frequency response of the ramp buffer amplifier and the comparator distort the shape of the ramp, so that the effective ramp value at the time the comparator trips differs from the intended value, resulting in errors. Allowing increased settling time decreases the quantization speed, while increasing the bandwidth increases the noise. The FPN problem is solved by breaking the ramp into two portions, with some fraction of the available code values allocated to a linear ramp and the remainder to a quadratic ramp. To avoid large transients, both the value and the slope of the linear and quadratic portions should be equal where they join. The span of the linear portion must cover the minimum offset, but not necessarily the maximum, since the fraction of the pixels above the upper limit will still be correctly quantized, albeit with increased quantization noise. The required linear span, maximum signal and ratio of quantization noise to shot noise at high signal, along with the continuity requirement, determines the number of code values that must be allocated to each portion. The distortion problem is solved by using a lookup table to convert captured code values back to signal levels. The values in this table will be similar to the intended ramp value, but with a correction for the finite bandwidth effects. Continuous-time comparators are used, and their bandwidth is set below the step rate, which smoothes the ramp and reduces the noise. No settling time is needed, as would be the case for clocked comparators, but the low bandwidth enhances the distortion of the non-linear portion. This is corrected by use of a return lookup table, which differs from the one used to generate the ramp. The return lookup table is obtained by calibrating against a stepped precision DC reference. This results in a residual non-linearity well below the quantization noise. This method can also compensate for differential non-linearity (DNL) in the DAC used to generate the ramp. The use of a ramp with a combination of linear and quadratic portions for a single-slope ADC is novel. The number of steps is minimized by keeping the step size just below the photon shot noise. This in turn maximizes the speed of the conversion. High resolution is maintained by keeping small quantization steps at low signals, and noise is minimized by allowing the lowest analog bandwidth, all without increasing the quantization noise. A calibrated return lookup table allows the system to maintain excellent linearity.

Wrigley, Chris J.↗

Characterization of acoustic disturbances in linearly sheared flows

The equation describing the plane wave propagation, the stability, or the rectangular duct mode characteristics in a compressible inviscid linearly sheared parallel, but otherwise homogeneous flow, is shown to be reducible to Whittaker's equation. The resulting solutions, which are real, viewed as functions of two variables, depend on a parameter and an argument, the values of which have precise physical meanings depending on the problem. The exact solutions in terms of Whittaker functions are used to obtain a number of known results of plane wave propagation and stability in linearly sheared flows as limiting cases in which the speed of sound goes to infinity (incompressible limit) or the shear layer thickness, or wave number, goes to zero (vortex sheet limit). The usefulness of the exact solutions is then discussed in connection with the problems of plane wave propagation and the stability of a finite thickness layer with a linear velocity profile.

Koutsoyannis, S. P.↗

Modification to Scott's method for direct integration of gaseous transmission to improve speed and accuracy

Two limitations to the method of Scott for rapid direct integration of gaseous transmission in a finite spectral interval over an inhomogeneous optical path are discussed. One of these is the inability to use a precomputed form of the Voigt line shape, which can be parameterized via the linewidth alone. This is rectified using the Voigt shape approximation of Kielkopf, which involves a linear combination of three functions which can be pre-computed: the Lorentz shape, the Doppler shape and a correction function. The other limitation involves situations where lines must be treated individually, either line positions must be known with great precision or a uniform value for the widths of all lines in the spectral interval does not describe the transmission adequately. In this case, the strongest lines are treated individually, while the remainder are treated assuming uniform line widths and line positions coincident with positions on the frequency grid, i.e. following Scott's method directly.

Orton, G. S.↗

Calculation of area-averaged vertical profiles of the horizontal wind velocity from volume-imaging lidar data

Area-averaged horizontal wind measurements are derived from the motion of spatial inhomogeneities in aerosol backscattering observed with a volume-imaging lidar. Spatial averaging provides high precision, reducing sample variations of wind measurements well below the level of turbulent fluctuations, even under conditions of very light mean winds and strong convection or under the difficult conditions represented by roll convection. Wind velocities are measured using the two-dimensional spatial cross correlation computed between successive horizontal plane maps of aerosol backscattering, assembled from three-dimensional lidar scans. Prior to calculation of the correlation function, three crucial steps are performed: (1) the scans are corrected for image distortion by the wind during a finite scan time; (2) a temporal high pass median filtering is applied to eliminate structures that do not move with the wind; and (3) a histogram equalization is employed to reduce biases to the brightest features.

Schols, J. L.↗

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗

Mechanical behavior of thermal barrier coatings for gas turbine blades

Plasma-sprayed thermal barrier coatings (TBCs) will enable turbine components to operate at higher temperatures and lower cooling gas flow rates; thereby improving their efficiency. Future developments are limited by precise knowledge of the material properties and failure mechanisms of the coating system. Details of this nature are needed for realistic modeling of the coating system which will, in turn, promote advancements in coating technology. Complementary experiments and analytical modeling which were undertaken in order to define and measure the important failure processes for plasma-sprayed coatings are presented. The experimental portion includes two different tests which were developed to measure coating properties. These are termed tensile adhesion and acoustic emission tests. The analytical modeling section details a finite element method which was used to calculate the stress distribution in the coating system. Some preliminary results are presented.

Berndt, C. C.↗

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, M. G.↗

Analysis of artificial viscosity effects on reacting flows using a spectral multidomain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, Michele G.↗

A proof of the Woodward-Lawson sampling method for a finite linear array

An extension of the continuous aperture Woodward-Lawson sampling theorem has been developed for a finite linear array of equidistant identical elements with arbitrary excitations. It is shown that by sampling the array factor at a finite number of specified points in the far field, the exact array factor over all space can be efficiently reconstructed in closed form. The specified sample points lie in real space and hence are measurable provided that the interelement spacing is greater than approximately one half of a wavelength. This paper provides insight as to why the length parameter used in the sampling formulas for discrete arrays is larger than the physical span of the lattice points in contrast with the continuous aperture case where the length parameter is precisely the physical aperture length.

Somers, Gary A.↗

Chaos and Beyond in a Water Filled Ultrasonic Resonance System

Finite amplitude ultrasonic wave resonances in a one-dimensional liquid-filled cavity, formed by a narrow band transducer and a plane reflector, are reported. The resonances are observed to include not only the expected harmonic and subharmonic signals (1,2) but chaotic signals as well. The generation mechanism requires attaining a threshold value of the driving amplitude that the liquid-filled cavity system becomes sufficiently nonlinear in response. The nonlinear features of the system were recently investigated via the construction of an ultrasonic interferometer having optical precision. The transducers were compressional, undamped quartz and lithium niobate crystals having the frequency range 1-10 MHz, driven by a high power amplifier. Both an optical diffraction system to characterize the diffraction pattern of laser light normally incident to the cavity and a receiving transducer attached to an aligned reflector with lapped flat and parallel surfaces were used to assess the generated resonance response in the cavity. At least 5 regions of excitation are identified.

Lazlo, Adler↗

Eigenvalues of the Laplacian of a graph

Let G be a finite undirected graph with no loops or multiple edges. The Laplacian matrix of G, Delta(G), is defined by Delta sub ii = degree of vertex i and Delta sub ij = -1 if there is an edge between vertex i and vertex j. The structure of the graph G is related to the eigenvalues of Delta(G); in particular, it is proved that all the eigenvalues of Delta(G) are nonnegative, less than or equal to the number of vertices, and less than or equal to twice the maximum vertex degree. Precise conditions for equality are given.

Anderson, W. N., Jr.↗

Thermal calibration and analysis of a Bridgman flight furnace

A combination of computer modeling and experimental measurements is proposed which permits the calibration of the temperature profile for a spacecraft flight furnace. The Fluent code based on the finite difference method uses the furnace set-point temperatures as boundary conditions to model a Bridgman directional solidification furnace. Experimental measurements are employed as a comparison for the numerical analysis, and the radial gradients and wall temperatures are selected when the calculated values are sufficiently close to those from measurements. The calculations are found to agree with experimental data, and some results of the modeling are given which relate to specific furnace parameters. The CFD approach facilitates the transposition of furnace calibration to microgravitational operation in addition to providing precise initial radial temperature gradients and furnace-wall temperatures.

Knuteson, D. J.↗

Bi-directional reflectance of finite and infinite clouds

The transfer of solar irradiance in plane parallel and broken cloud fields is modeled using a Monte Carlo method. The angular distribution pattern of radiances exiting the cloud layer is studied with varying cloud geometries, optical thicknesses, cloudiness and solar zenith angles. A rather large anisotropy of the reflected flux is found, usually increasing with solar zenith angle and with patterns that strongly depend on cloud geometry. The main features are: a local maximum of reflected intensity in the forward direction for all cases; a limb darkening for the plane parallel case; and a limb brightening and a local maximum of reflected intensity in the backward direction for broken clouds. The method is used to study the cloud reflectance sensitivity to various parameters. A more precise description of cloud field internal and external structure is needed in order to obtain accurate bi directional reflectance diagrams.

Breon, Francois-Marie↗

Advances in Non-contact Measurement of Creep Properties

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility, and heated with a laser. The laser is aligned off-center so that the absorbed photons transfer their momentum to the sample, causing it to rotate at up to 250,000+ RPM. The rapid rotation loads the sample through centripetal acceleration, causing it to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the noncontact method exploits stress gradients within the sample to determine the stress exponent in a single test. This method was validated in collaboration with the University of Tennessee for niobium at 1985 C, with agreement within the uncertainty of the conventional measurements. A similar method is being employed on Ultra-High-Temperature ZrB2- SiC composites, which may see application in rocket nozzles and sharp leading edges for hypersonic vehicles.

Hyers, Robert W.↗