SEARCH · Engineering Papers
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.
Software to evaluate functions of one or two arguments is too often tested or certified, or the accuracy assessed, by ad hoc methods. Systematic methods to test, certify or assess the accuracy of single- or double-precision functions of one or two arguments, and supporting software, are described here.
Rapid generation of synthetic aperture radar (SAR) based flood extent and flood depth maps provide valuable data in disaster response efforts. We present a simple but powerful method using dual polarimetric SAR imagery. A RGB false-color map is generated using pre- and post-flooding imagery, allowing operators to distinguish between existing standing water in pre-flooding data and recently flooded areas. This method works very well in areas of standing water, while large omission errors can be seen in urban areas due to the double-bounce effect. A flood depth map is also estimated by using an external DEM. Compared with FEMA flood product, flood water depth from the proposed method showed low bias with small dispersion. This automatic flood mapping system will contribute to the rapid assessment for disaster relief efforts.
Rapid generation of synthetic aperture radar (SAR) based flood extent and flood depth maps provide valuable data in disaster response efforts. We present a simple but powerful method using dual-polarimetric SAR imagery. A RGB false-color map is generated using pre- and post-flooding imagery, allowing operators to distinguish between existing standing water in pre-flooding data and recently flooded areas. This method works very well in areas of standing water, while large omission errors can be seen in urban areas due to the double-bounce effect. A flood depth map is also estimated by using an external DEM. Compared with FEMA flood product, flood water depth from the proposed method showed low bias with small dispersion. This automatic flood mapping system will contribute to the rapid assessment for disaster relief efforts.
A versatile and efficient technique for the solution of the resonance line scattering problem with frequency redistribution in planetary atmospheres is introduced. Similar to the doubling approach commonly used in monochromatic scattering problems, the technique has been extended to include the frequency dependence of the radiation field. Methods for solving problems with external or internal sources and coupled spectral lines are presented, along with comparison of some sample calculations with results from Monte Carlo and Feautrier techniques. The doubling technique has also been applied to the solution of resonance line scattering problems where the R-parallel redistribution function is appropriate, both neglecting and including polarization as developed by Yelle and Wallace (1989). With the constraint that the atmosphere is illuminated from the zenith, the only difficulty of consequence is that of performing precise frequency integrations over the line profiles. With that problem solved, it is no longer necessary to use the Monte Carlo method to solve this class of problem.
A complete theoretical electric dipole moment function mu(r) is computed for OH(X 2Pi) at the full configuration-interaction (FCI) level for DZP + diffuse Gaussian basis set. The energy and dipole moment are also determined with 2d functions on oxygen at seven r values. The following methods are compared to the FCI for both the oxygen 1d and 2d basis sets: single-reference singles-plus-doubles configuration interaction (SDCI), SDCI with Davidson's correction, the coupled pair functional approach (CPF) and a modified form of CPF, CASSCF, and CASSCF-MRSDCI, and the external contracted CI(CCI) method. The dipole moments are evaluated both as an expectation value and by finite-field methods. The results support the superiority of evaluating the dipole moment as an energy derivative rather than as an expectation value.
Air pollutants are identified, and the needs for their measurement from satellites and aircraft are discussed. An assessment is made of the properties of these pollutants and of the normal atmosphere, including interactions with light of various wavelengths and the resulting effects on transmission and scattering of optical signals. The possible methods for active remote measurement are described; the relative performance capabilities of double-ended and single-ended systems are compared qualitatively; and the capabilities of the several single-ended or backscattering techniques are compared quantitatively. The differential-absorption lidar (DIAL) technique is shown to be superior to the other backscattering techniques. The lidar system parameters and their relationships to the environmental factors and the properties of pollutants are examined in detail. A computer program that models both the atmosphere (including pollutants) and the lidar system is described. The performance capabilities of present and future lidar components are assessed, and projections are made of prospective measurement capabilities for future lidar systems. Following a discussion of some important operational factors that affect both the design and measurement capabilities of airborne and satellite-based lidar systems, the extensive analytical results obtained through more than 1000 individual cases analyzed with the aid of the computer program are summarized and discussed. The conclusions are presented. Recommendations are also made for additional studies to investigate cases that could not be explored adequately during this study.
Two large, complex deep-focus earthquakes are analyzed, using body waves, to determine the relative locations and equivalent forces of the component events. Details of the source mechanism determination are given, with special emphasis on the means of comparing results for double-couple and more general models. The formalism includes a description of a method for calculating the confidence ellipses around the principal axes of the moment tensor. An interpretation of the results suggests that the mechanism of deep-focus earthquakes is controlled by the buoyancy of the slab and the change in viscosity at the 650 km discontinuity.
This program was undertaken to initiate the development of a test method for testing adhesive joints in metal-adhesive-composite systems. The uniform double cantilever beam (UDCB) and the width tapered beam (WTB) specimen geometries were evaluated for measuring Mode I fracture toughness in these systems. The WTB specimen is the preferred geometry in spite of the fact that it is more costly to machine than the UDCB specimen. The use of loading tabs attached to thin sheets of composites proved to be experimentally unsatisfactory. Consequently, a new system was developed to load thin sheets of adherends. This system allows for the direct measurement of displacement along the load line. In well made joints separation occurred between the plies rather than in the adhesive.
Measurement does not affect joint parameters. Verification of analytical and strength-prediction methods for bolted composite joints based generally on data obtained experimentally from double-lap-joint specimens. In mechanically fastened joints, stresses maximal at fastener holes. Ability to measure accurately hole elongations without affecting joint parameters provides better understanding of elastic and plastic behavior of joint material leading to failure mechanisms in mechanically fastened joints required for design of more-efficient, lightweight composite joints.
Several International Solar Terrestrial Physics (ISTP) missions require the design of complex gravity assisted trajectories in order to investigate the interaction of the solar wind with the Earth's magnetic field. These trajectories present a formidable trajectory design and optimization problem. The philosophy and methodology that enable an analyst to design and analyse such trajectories are discussed. The so called 'floating end point' targeting, which allows the inherently nonlinear multiple body problem to be solved with simple linear techniques, is described. The combination of floating end point targeting with analytic approximations with a Newton method targeter to achieve trajectory design goals quickly, even for the very sensitive double lunar swingby trajectories used by the ISTP missions, is demonstrated. A multiconic orbit integration scheme allows fast and accurate orbit propagation. A prototype software tool, Swingby, built for trajectory design and launch window analysis, is described.
Several International Solar-Terrestrial Physics (ISTP) missions require the design of complex gravity-assisted trajectories in order to investigate the interaction of the solar wind with the earth's magnetic field. These trajectories present a formidable trajectory design and optimization problem. This paper discusses the philosophy and methodology that enable an analyst to design and analyze such trajectories. The authors describe what is called 'floating end-point' targeting, which allows the inherently nonlinear multiple body problem to be solved with simple linear techniques. This paper demonstrates how floating end-point targeting combines analytic approximations with a Newton method targeter to achieve trajectory design goals quickly, even for the very sensitive double lunar swingby trajectories used by the ISTP missions. A Multiconic orbit integration scheme allows fast and accurate orbit propagation. This paper also describes a prototype software tool, Swingby, which has been built for trajectory design and launch window analysis.
The energy absorption response and crushing characteristics of geometrically scaled graphite-Kevlar epoxy composite plates were investigated. Three different trigger mechanisms including chamfer, notch, and steeple geometries were incorporated into the plate specimens to initiate crushing. Sustained crushing was achieved with a simple test fixture which provided lateral support to prevent global buckling. Values of specific sustained crushing stress (SSCS) were obtained which were comparable to values reported for tube specimens from previously published data. Two sizes of hybrid plates were fabricated; a baseline or model plate, and a full-scale plate with in-plane dimensions scaled by a factor of two. The thickness dimension of the full-scale plates was increased using two different techniques; the ply-level method in which each ply orientation in the baseline laminate stacking sequence is doubled, and the sublaminate technique in which the baseline laminate stacking sequence is repeated as a group. Results indicated that the SSCS is independent of trigger mechanism geometry. However, a reduction in the SSCS of 10-25 percent was observed for the full-scale plates as compared with the baseline specimens, indicating a scaling effect in the crushing response.
A shock capturing finite difference method for systems of hyperbolic conservation laws is presented which avoids the need to solve Riemann problems while being competitive in performance with other current methods. A staggered spatial mesh is employed, so that complicated nonlinear waves generated at cell interfaces are averaged over cell interiors at the next time level. The full method combines to form a conservative version of the modified method of characteristics. The advantages of the method are discussed, and numerical results are presented for the two-dimensional double ellipse problem.
The energy absorption response and crushing characteristics of geometrically scaled graphite-Kevlar epoxy composite plates were investigated. Two different trigger mechanisms including notch, and steeple geometries were incorporated into the plate specimens to initiate crushing. Sustained crushing was achieved with a new test fixture which provided lateral support to prevent global buckling. Values of specific sustained crushing stress (SSCS) were obtained which were lower than values reported for tube specimens from previously published data. Two sizes of hybrid plates were fabricated; a baseline or model plate, and a full-scale plate with inplane dimensions scaled by a factor of two. The thickness dimension of the full-scale plates was increased using two different techniques: the ply-level method in which each ply orientation in the baseline laminate stacking sequence is doubled, and the sublaminate technique in which the baseline laminate stacking sequence is repeated as a group. Results indicated that the SSCS has a small dependence on trigger mechanism geometry. However, a reduction in the SSCS of 10-25% was observed for the full-scale plates as compared with the baseline specimens, indicating a scaling effect in the crushing response.
This paper presents the results of a numerical study of shock interactions resulting from Mach 10 flow about sharp double cones. Computations are made by using the direct simulation Monte Carlo (DSMC) method of Bird. The sensitivity and characteristics of the interactions are examined by varying flow conditions, model size, and configuration. The range of conditions investigated includes those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel.
The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.
Incident shock-wave velocities were measured in the Langley 6-inch expansion tube, operated as a shock tube, with air, argon, carbon dioxide, and helium as test gases. Unheated helium was used as the driver gas and most data were obtained at pressures of approximately 34 and 54 MN/sq m. A range of pressure ratio across the diaphragm was obtained by varying the quiescent test-gas pressure, for a given driver pressure, from 0.0276 to 34.5 kN/sq m. Single- and double-diaphragm modes of operation were employed and diaphragms of various materials tested. Shock velocity was determined from microwave interferometer measurements, response of pressure transducers positioned along interferometer measurements, response of pressure transducers positioned along the driven section (time-of-arrival gages), and to a lesser extent, measured tube-wall pressure. Velocities obtained from these methods are compared and limitations of the methods discussed. The present results are compared with theory and the effects of diaphragm mode (single or double diaphragm), diaphragm material, heating of the driver gas upon pressurization of the driver section, diaphragm opening time, interface mixing, and two-dimensional (nonplanar) flow are discussed.
An improved panel method for the solution of three dimensional flow and wing and wing-body combinations with leading edge vortex separation is presented. The method employs a three dimensional inviscid flow model in which the configuration, the rolled-up vortex sheets, and the wake are represented by quadratic doublet distributions. The strength of the singularity distribution as well as shape and position of the vortex spirals are computed in an iterative fashion starting with an assumed initial sheet geometry. The method calculates forces and moments as well as detail surface pressure distributions. Improvements include the implementation of improved panel numerics for the purpose of elimination the highly nonlinear effects of ring vortices around double panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A complete description of the method is included. A variety of cases generated by the computer program implementing the method are presented which verify the mathematical assumptions of the method and which compare computed results with experimental data to verify the underlying physical assumptions made by the method.