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.
Algorithms are developed for computing stresses in a semi-infinite body when loaded by a uniform pressure acting over a circular area. The algorithm allows easy determination of any stress component in a semi-infinite body having a known Poisson's ratio. Example curves are plotted for Portland cement grout and metal representative values.
The measurements of wind speed and rain rate over a large dynamic range has been demonstrated by using passive microwave remote sensing measurements by a C-band radiometer flying above a hurricane. The measurements were obtained by using a unique variable frequency radiometer, the Stepped Frequency Microwave Radiometer (SFMR), which is capable of near simultaneous measurements at several widely spaced frequencies in a near octave bandwidth in the C-band microwave frequency region. In addition, an algorithm was developed which computes the surface wind speed and rain rate from the SFMR measured radiometric brightness temperature at two frequencies, and this algorithm was tested using experimental data collected during flights into Hurricane Allen during 1980. Results showed a reasonable agreement between the SFMR measured data and measurements obtained by weather radar.
The atmospheric backscatter coefficient, beta, measured with an airborne CO Laser Doppler Velocimeter (LDV) system operating in a continuous wave, focussed model is discussed. The Single Particle Mode (SPM) algorithm, was developed from concept through analysis of an extensive amount of data obtained with the system on board a NASA aircraft. The SPM algorithm is intended to be employed in situations where one particle at a time appears in the sensitive volume of the LDV. In addition to giving the backscatter coefficient, the SPM algorithm also produces as intermediate results the aerosol density and the aerosol backscatter cross section distribution. A second method, which measures only the atmospheric backscatter coefficient, is called the Volume Mode (VM) and was simultaneously employed. The results of these two methods differed by slightly less than an order of magnitude. The measurement uncertainties or other errors in the results of the two methods are examined.
A conventional prime factor discrete Fourier transform (DFT) algorithm is used to realize a discrete Fourier-like transform on the finite field, GF(q sub n). A pipeline structure is used to implement this prime factor DFT over GF(q sub n). This algorithm is developed to compute cyclic convolutions of complex numbers and to decode Reed-Solomon codes. Such a pipeline fast prime factor DFT algorithm over GF(q sub n) is regular, simple, expandable, and naturally suitable for VLSI implementation. An example illustrating the pipeline aspect of a 30-point transform over GF(q sub n) is presented.
The objective of the advanced detection, isolation, and accommodation (ADIA) program is to improve the overall demonstrated reliability of digital electronic control systems for turbine engines. For this purpose, algorithms were developed which detect, isolate, and accommodate sensor failures using analytical redundancy. Preliminary results of a full scale engine demonstration of the ADIA algorithm are presented. Minimum detectable levels of sensor failures for an F100 turbofan engine control system are determined and compared to those obtained during a previous evaluation of this algorithm using a real-time hybrid computer simulation of the engine.
The objective of the advanced detection, isolation, and accommodation (ADIA) program is to improve the overall demonstrated reliability of digital electronic control systems for turbine engines. For this purpose, algorithms were developed which detect, isolate, and accommodate sensor failures using analytical redundancy. Preliminary results of a full scale engine demonstration of the ADIA algorithm are presented. Minimum detectable levels of sensor failures for an F100 turbofan engine control system are determined and compared to those obtained during a previous evaluation of this algorithm using a real-time hybrid computer simulation of the engine.
A conventional prime factor discrete Fourier transform (DFT) algorithm is used to realize a discrete Fourier-like transform on the finite field, GF(q sub n). This algorithm is developed to compute cyclic convolutions of complex numbers and to decode Reed-Solomon codes. Such a pipeline fast prime factor DFT algorithm over GF(q sub n) is regular, simple, expandable, and naturally suitable for VLSI implementation. An example illustrating the pipeline aspect of a 30-point transform over GF(q sub n) is presented.
A large portion of the earth is covered by thin semi-transparent cirrus cloud. The cirrus results from the natural injection of moisture into the upper troposphere by deep convection (i.e., anvils) and from man-made moisture injected into the upper troposphere by jet aircraft. Although most cirrus clouds are semi-transparent to infrared wavelengths, their heights, thicknesses, and spectral absorption properties must be known in order to retrieve atmospheric temperature and moisture profiles from the data. An algorithm is developed for accounting for the radiative properties of semi-transparent cloud in the retrieval of vertical temperature and moisture profiles. The algorithm is to be applied to the NASA ER2 HIS data collected during the FIRE cirrus field program.
Two algorithms were developed which utilized multi-sensor analysis techniques to complement the current Space Shuttle Main Engine (SSME) safety monitoring system. The first algorithm analyzed the accumulative error between actual and predicted values of the engine parameter set, while the second algorithm combined these error terms into a response pattern and correlated each pattern with a standard pattern. These algorithms were applied to twelve SSME anomalous test firings and were found to produce improved failure detection times in eight of those twelve compared to the current engine safety monitoring system. Of the eight detected anomalous test firings, the first algorithm detected all eight, while the second algorithm detected seven of the eight. No false alarms were indicated by either algorithm for twelve nominal test firings. An initial parametric study of these algorithms for optimized parameter selection is presented and algorithm robustness to sensor failure is demonstrated.
A method has been developed to determine the oceanic daily mean mixed layer depth from satellite observations and a mixed layer thermal inertia (MLTI) model. The algorithms were developed to use remotely-sensed values of sea surface temperature, albedo, and surface wind speeds to calculate the thermal inertia and to predict changes in subsurface diurnal mixed layer depth. The MLTI model, based on a mixed layer model of the upper ocean, has been used to simulate the diurnal mixing process and thermal inertia distribution in the Sargasso Sea around 34 deg N, 70 deg W. Sea surface temperature and albedo have been obtained from the NOAA7-AVHRR images. Surface wind speeds have been derived from the Scanning Multichannel Microwave Radiometer (SMMR) aboard Nimbus 7. Image processing was performed for images gathered between June and July 1982. The daily mean mixed layer depths predicted by the MLTI model agree well with data gathered at the LOTUS mooring located in the Sargasso Sea. This suggests that vertical mixing is the dominant physical process that controls the thermal inertia distribution in the midocean, far from major current systems, and that remote sensing is a promising tool to study such upper ocean processes.
The application of probability density function (pdf) turbulence models is addressed. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional computational fluid dynamic (CFD) flow solver with the Monte Carlo simulation of the pdf evolution equation was developed. The algorithm was validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames was carried out using this algorithm. Numerical results compared favorably with experimental data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.
Near infrared spectroscopy (NIRS) is commonly used to measure muscle oxygenation during exercise and recovery. Current NIRS algorithms do not account for variation in water content and optical pathlength during exercise. The current effort attempts to validate a newly developed NIRS algorithm during rhythmic handgrip exercise and recovery. Six female subjects, aver age 28 +/- 6 yrs, participated in the study. A venous catheter was placed in the retrograde direction in the antecubital space. A NIRS sensor with 30 mm source-detector separation was placed on the flexor digitorum profundus. Subjects performed two 5-min bouts of rhythmic handgrip exercise (2 s contraction/1 s relaxation) at 15% and 30% of maximal voluntary contraction. Venous blood was sampled before each bout, during the last minute of exercise, and after 5 minutes of recovery. Venous oxygen saturation (SvO2) was measured with a I-stat CG-4+ cartridge. Spectra were collected between 700-900 nm. A modified Beer's Law formula was used to calculate the absolute concentration of oxyhemoglobin (HbO2), deoxyhemoglobin (Hb) and water, as well as effective pathlength for each spectrum. Muscle oxygen saturation (SmO2) was calculated from the HbO2 and Hb results. The correlation between SvO2 and SmO2 was determined. Optical pathlength and water varied significantly during each exercise bout, with pathlength increasing approximately 20% and water increasing about 2%. R2 between blood and muscle SO2 was found to be 0.74, the figure shows the relationship over SvO2 values between 22% and 82%. The NIRS measurement was, on average, 6% lower than the blood measurement. It was concluded that pathlength changes during exercise because muscle contraction causes variation in optical scattering. Water concentration also changes, but only slightly. A new NIRS algorithm which accounts for exercise-induced variation in water and pathlength provided an accurate assessment of muscle oxygen saturation before, during and after exercise.
NOAA platforms provide the longest period of terrestrial observation since the 1980s. The progress in calibration, atmospheric corrections and physically based land retrieval offers the opportunity to reprocess these data for extending terrestrial product time series. Within the Quality Assurance for Essential Climate Variables (QA4ECV) project, the black-sky Joint Research Centre (JRC)-fraction of absorbed photosynthetically active radiation (FAPAR) algorithm was developed for the AVHRR sensors on-board NOAA-07 to -16 using the Land Surface Reflectance Climate Data Record. The retrieval algorithm was based on the radiative transfer theory, and uncertainties were included in the products. We proposed a time and spatial composite for providing both 10-day and monthly products at 0.05º × 0.05º. Quality control and validation were achieved through benchmarking against third-party products, including Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) datasets produced with the same retrieval algorithm. Past ground-based measurements, providing a proxy of FAPAR, showed good agreement of seasonality values over short homogeneous canopies and mixed vegetation. The average difference between SeaWiFS and QA4ECV monthly products over 2002–2005 is about 0.075 with a standard deviation of 0.091. We proposed a monthly linear bias correction that reduced these statistics to 0.02 and 0.001. The complete harmonized long-term time series was then used to address its fitness for the purpose of analysis of global terrestrial change.
The use of vectorial parameterization to create geometrical representations in computational models has a large number of applications. One particular application is the calculation of the 3D rotational motion of rigid bodies, that could be used for the spatial location estimation from objects. Provided the algebraic nature of this problem, it could benefit from Quantum Computing, in particular several vectors could be superposed to be transformed with a single operation, providing a quantum processing advantage. In this article, we propose an implementation of a Quantum Computing algorithm to compute Euler-Rodrigues Parameters to model rigid body rotations to transform arbitrary functions, rotating multiple vectors in superposition. We developed this algorithm using Qiskit, taking into account the limitations imposed by the current Noisy Intermediate Scale Quantum (NISQ) devices, such as the reduced number of qubits available and the limited coherence time.
The application of recently developed function-space Davidon-type techniques to the shuttle ascent trajectory optimization problem is discussed along with an investigation of the recently developed PRAXIS algorithm for parameter optimization. At the outset of this analysis, the major deficiency of the function-space algorithms was their potential storage problems. Since most previous analyses of the methods were with relatively low-dimension problems, no storage problems were encountered. However, in shuttle trajectory optimization, storage is a problem, and this problem was handled efficiently. Topics discussed include: the shuttle ascent model and the development of the particular optimization equations; the function-space algorithms; the operation of the algorithm and typical simulations; variable final-time problem considerations; and a modification of Powell's algorithm.
Deadbeat control theory gives a feedback solution to the time optimal control of discrete time systems. Experience has shown the results to be impractical because they ignore bounds on the actuator strength. This paper develops two algorithms for generating time optimal control in feedback form for discrete systems with bounded controls. The results are also applicable for generating recovery regions and the set of reachable states. For multiple control problems a method of generating sublayers is developed which decreases off-line and on-line computational effort. Two algorithms are presented with somewhat different computational and storage requirements. The algorithms are practical within certain dimension constraints, and are natural for implementation with parallel processing.
Gynecologic cancers are often asymmetric, yet current Ir-192 brachytherapy techniques provide only limited radial modulation of the dose. The shielded solutions investigated here solve this by providing the ability to modulate between highly asymmetric and radially symmetric dose distributions at a given location. To find applicator designs that can modulate between full dose and less than 50% dose, at the dimensions of the urethra, a 2D calculation algorithm was developed to narrow down the search space. Two shielding design types were then further investigated using Monte Carlo and Boltzmann-solver dose calculation algorithms. 3D printing techniques using ISO 10993 certified biocompatible plastics and 3D printable tungsten-loaded plastics were tested. It was also found that shadowing effects set by the shape of the shielding cannot be easily modulated out, hence careful design is required. The shielded applicator designs investigated here, allow for reduction of the dose by over 50% at 5 mm from the applicator surface in desired regions, while also allowing radially symmetric dose with isodose line deviations less than 0.5 mm from circular. The shielding designs were also chosen with treatment delivery time in mind. Treatment times for these shielded designs were found to be less than 1.4 times longer than a 6-channel unshielded cylinder for the equivalent fully symmetric dose distribution. The 2D calculation methods developed here provide a simple way to rapidly evaluate shielding designs, while the 3D printing techniques also allow for devices with novel shapes to be rapidly prototyped. Both TOPAS Monte Carlo and Acuros BV calculations show that significant dose shaping and organ at risk sparing can be achieved without significantly compromising the plan in regions that require the full dose.
Low density, high strength, and high creep and oxidation resistance properties of ceramic matrix composites (CMCs) make them an ideal choice for use in extreme environments in space and military applications. This paper presents a detailed characterization study of structural and manufacturing flaws in Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs at different length-scales. Energy-dispersive spectroscopy (EDS) is used for the chemical characterization of the material’s elemental constituents. High-resolution multiscale graphs obtained from scanning electron microscope (SEM) and confocal laser scanning microscope (LSM) are used to characterize the distribution and morphology of defects at different length scales. This is followed by the classification and quantification of the common manufacturing defects. An image processing algorithm based on the image segmentation process is developed to quantify the variability of various scale-dependent architectural parameters. Finally, a three-dimensional stochastic representative volume element (SRVE) generation algorithm is developed to provide precise representations of material textures at multiple length scales. The developed algorithm accurately accounts for material features and flaws based on a range of multiscale structural and defects characterization results.