Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Comparisons of the Solutions for the Category 3--Problem 2: Cascade-Gust Interaction

Six different solutions were submitted for this benchmark problem. These were obtained using a variety of methods that can be conveniently categorized in two main groups, a nonlinear time-domain group and a linearized frequency-domain group. The first includes solutions submitted by (1) Hixon, (2) Nallasamy et. al, (3) Shieh et. al, and (4) Wang et. al, and the second includes solutions submitted by (5) Coupland and (6) Serrano et. al. Methods (1) and (2) use sixth order compact differencing schemes and the rest are essentially second order in space. With the exception of the solution submitted by Shieh et. al, all are individually discussed in great detail in the workshop proceedings. Comparisons of the submitted solutions with the benchmark solution are presented below. Due to differences in the level of solution detail provided to the author by the participants, the comparisons do not always include results from all submissions. It should be noted at the outset that, since the benchmark solution itself was numerically computed, the comparisons are somewhat subjective. In order to provide maximum latitude for the participants of the workshop, no restrictions were placed on the type of method that could be used to solve the problem. Neither was there were any stipulations to use a particular grid topology or grid density. Therefore, without a detailed study of the critical features of the computed solutions, it is not possible to make concrete statements about the relative merits of one method over another. Such a study is beyond the scope of the current exercise, especially since complete flowfield details were not provided to the author by all participants. Instead a package, containing the information about the benchmark solution (both the steady and unsteady parts of it), is included on the proceedings CD should the authors who submitted solutions for this problem wish to examine in detail the benchmark solution and compare their results to it.

Envia, Edmane↗
NASA NTRS2021

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗
NASA NTRS1990

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 4: Maintenance document (version 3.0)

The Maintenance Document Version 3.0 is a guide to the PAN AIR software system, a system which computes the subsonic or supersonic linear potential flow about a body of nearly arbitrary shape, using a higher order panel method. The document describes the overall system and each program module of the system. Sufficient detail is given for program maintenance, updating, and modification. It is assumed that the reader is familiar with programming and CRAY computer systems. The PAN AIR system was written in FORTRAN 4 language except for a few CAL language subroutines which exist in the PAN AIR library. Structured programming techniques were used to provide code documentation and maintainability. The operating systems accommodated are COS 1.11, COS 1.12, COS 1.13, and COS 1.14 on the CRAY 1S, 1M, and X-MP computing systems. The system is comprised of a data base management system, a program library, an execution control module, and nine separate FORTRAN technical modules. Each module calculates part of the posed PAN AIR problem. The data base manager is used to communicate between modules and within modules. The technical modules must be run in a prescribed fashion for each PAN AIR problem. In order to ease the problem of supplying the many JCL cards required to execute the modules, a set of CRAY procedures (PAPROCS) was created to automatically supply most of the JCL cards. Most of this document has not changed for Version 3.0. It now, however, strictly applies only to PAN AIR version 3.0. The major changes are: (1) additional sections covering the new FDP module (which calculates streamlines and offbody points); (2) a complete rewrite of the section on the MAG module; and (3) strict applicability to CRAY computing systems.

Purdon, David J.↗
NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗
NASA NTRS2015

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗
NASA NTRS1995

3-D Reconstruction of Macular Type II Cell Innervation Patterns in Space-Flight and Control Rats

A semiautomated method for reconstructing objects from serial thin sections has been developed in the Biocomputation Center. The method is being used to completely, for the first time, type II hair cells and their innervations. The purposes are to learn more about the fundamental circuitry of the macula on Earth and to determine whether changes in connectivities occur under space flight conditions. Data captured directly from a transmission electron microscope via a video camera are sent to a graphics workstation. There, the digitized micrographs are mosaicked into sections and contours are traced, registered and displayed by semiautomated methods. Current reconstructions are of type II cells from the medial part of rat maculas collected in-flight on the Space Life Sciences-2 mission, 4.5 hrs post-flight, and from a ground control. Results show that typical type II cells receive processes from tip to six nearby calyces or afferents. Nearly all processes are elongated and have bouton-like enlargements; some have numerous vesicles. Multiple (2 to 4) processes from a single calyx to a type II cell are common, and approximately 1/3 of the processes innervale 2 or 3 type II cells or a neighboring cluster. From 2% to 6% of the cells resemble type I cells morphologically but have demi-calyces. Thus far, increments in synaptic number in type II cells of flight rats are prominent along processes that supply two hair cells. It is clear that reconstruction methods provide insights into details of macular circuitry not obtainable by other techniques. The results demonstrate a morphological basis for interactions between adjacent receptive fields through feed back-feed forward connections, and for dynamic alterations in receptive field range and activity during preprocessing of linear acceleratory information by the maculas. The reconstruction method we have developed will find further applications in the study of the details of neuronal architecture of more complex systems, to seek out shared organizational properties or neuronal networks and to understand better localization of synaptic changes in altered environments.

Ross, Muriel Dorothy↗
NASA NTRS2012

Analyses of Cometary Silicate Crystals: DDA Spectral Modeling of Forsterite

Comets are the Solar System's deep freezers of gases, ices, and particulates that were present in the outer protoplanetary disk. Where comet nuclei accreted was so cold that CO ice (approximately 50K) and other supervolatile ices like ethane (C2H2) were preserved. However, comets also accreted high temperature minerals: silicate crystals that either condensed (greater than or equal to 1400 K) or that were annealed from amorphous (glassy) silicates (greater than 850-1000 K). By their rarity in the interstellar medium, cometary crystalline silicates are thought to be grains that formed in the inner disk and were then radially transported out to the cold and ice-rich regimes near Neptune. The questions that comets can potentially address are: How fast, how far, and over what duration were crystals that formed in the inner disk transported out to the comet-forming region(s)? In comets, the mass fractions of silicates that are crystalline, f_cryst, translate to benchmarks for protoplanetary disk radial transport models. The infamous comet Hale-Bopp has crystalline fractions of over 55%. The values for cometary crystalline mass fractions, however, are derived assuming that the mineralogy assessed for the submicron to micron-sized portion of the size distribution represents the compositional makeup of all larger grains in the coma. Models for fitting cometary SEDs make this assumption because models can only fit the observed features with submicron to micron-sized discrete crystals. On the other hand, larger (0.1-100 micrometer radii) porous grains composed of amorphous silicates and amorphous carbon can be easily computed with mixed medium theory wherein vacuum mixed into a spherical particle mimics a porous aggregate. If crystalline silicates are mixed in, the models completely fail to match the observations. Moreover, models for a size distribution of discrete crystalline forsterite grains commonly employs the CDE computational method for ellipsoidal platelets (c:a:b=8.14x8.14xl in shape with geometrical factors of x:y:z=1:1:10, Fabian et al. 2001; Harker et al. 2007). Alternatively, models for forsterite employ statistical methods like the Distribution of Hollow Spheres (Min et al. 2008; Oliveira et al. 2011) or Gaussian Random Spheres (GRS) or RGF (Gielen et al. 200S). Pancakes, hollow spheres, or GRS shapes similar to wheat sheaf crystal habit (e.g., Volten et al. 2001; Veihelmann et al. 2006), however, do not have the sharp edges, flat faces, and vertices seen in images of cometary crystals in interplanetary dust particles (IDPs) or in Stardust samples. Cometary forsterite crystals often have equant or tabular crystal habit (J. Bradley). To simulate cometary crystals, we have computed absorption efficiencies of forsterite using the Discrete Dipole Approximation (DDA) DDSCAT code on NAS supercomputers. We compute thermal models that employ a size distribution of discrete irregularly shaped forsterite crystals (nonspherical shapes with faces and vertices) to explore how crystal shape affects the shape and wavelength positions of the forsterite spectral features and to explore whether cometary crystal shapes support either condensation or annealing scenarios (Lindsay et al. 2012a, b). We find forsterite crystal shapes that best-fit comet Hale-Bopp are tetrahedron, bricks or brick platelets, essentially equant or tabular (Lindsay et al. 2012a,b), commensurate with high temperature condensation experiments (Kobatake et al. 2008). We also have computed porous aggregates with crystal monomers and find that the crystal resonances are amplified. i.e., the crystalline fraction is lower in the aggregate than is derived by fitting a linear mix of spectral features from discrete subcomponents, and the crystal resonances 'appear' to be from larger crystals (Wooden et al. 2012). These results may indicate that the crystalline mass fraction in comets with comae dominated by aggregates may be lower than deduced by popular methods that only emoy ensembles of discrete crystals.

Wooden, Diane↗
NASA NTRS1991

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗
NASA NTRS2009

Highly Scalable Matching Pursuit Signal Decomposition Algorithm

Matching Pursuit Decomposition (MPD) is a powerful iterative algorithm for signal decomposition and feature extraction. MPD decomposes any signal into linear combinations of its dictionary elements or atoms . A best fit atom from an arbitrarily defined dictionary is determined through cross-correlation. The selected atom is subtracted from the signal and this procedure is repeated on the residual in the subsequent iterations until a stopping criterion is met. The reconstructed signal reveals the waveform structure of the original signal. However, a sufficiently large dictionary is required for an accurate reconstruction; this in return increases the computational burden of the algorithm, thus limiting its applicability and level of adoption. The purpose of this research is to improve the scalability and performance of the classical MPD algorithm. Correlation thresholds were defined to prune insignificant atoms from the dictionary. The Coarse-Fine Grids and Multiple Atom Extraction techniques were proposed to decrease the computational burden of the algorithm. The Coarse-Fine Grids method enabled the approximation and refinement of the parameters for the best fit atom. The ability to extract multiple atoms within a single iteration enhanced the effectiveness and efficiency of each iteration. These improvements were implemented to produce an improved Matching Pursuit Decomposition algorithm entitled MPD++. Disparate signal decomposition applications may require a particular emphasis of accuracy or computational efficiency. The prominence of the key signal features required for the proper signal classification dictates the level of accuracy necessary in the decomposition. The MPD++ algorithm may be easily adapted to accommodate the imposed requirements. Certain feature extraction applications may require rapid signal decomposition. The full potential of MPD++ may be utilized to produce incredible performance gains while extracting only slightly less energy than the standard algorithm. When the utmost accuracy must be achieved, the modified algorithm extracts atoms more conservatively but still exhibits computational gains over classical MPD. The MPD++ algorithm was demonstrated using an over-complete dictionary on real life data. Computational times were reduced by factors of 1.9 and 44 for the emphases of accuracy and performance, respectively. The modified algorithm extracted similar amounts of energy compared to classical MPD. The degree of the improvement in computational time depends on the complexity of the data, the initialization parameters, and the breadth of the dictionary. The results of the research confirm that the three modifications successfully improved the scalability and computational efficiency of the MPD algorithm. Correlation Thresholding decreased the time complexity by reducing the dictionary size. Multiple Atom Extraction also reduced the time complexity by decreasing the number of iterations required for a stopping criterion to be reached. The Course-Fine Grids technique enabled complicated atoms with numerous variable parameters to be effectively represented in the dictionary. Due to the nature of the three proposed modifications, they are capable of being stacked and have cumulative effects on the reduction of the time complexity.

Christensen, Daniel↗
NASA NTRS1999

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗
NASA NTRS2013

Microgravity Drill and Anchor System

This work is a method to drill into a rock surface regardless of the gravitational field or orientation. The required weight-on-bit (WOB) is supplied by a self-contained anchoring mechanism. The system includes a rotary percussive coring drill, forming a complete sampling instrument usable by robot or human. This method of in situ sample acquisition using micro - spine anchoring technology enables several NASA mission concepts not currently possible with existing technology, including sampling from consolidated rock on asteroids, providing a bolt network for astronauts visiting a near-Earth asteroid, and sampling from the ceilings or vertical walls of lava tubes and cliff faces on Mars. One of the most fundamental parameters of drilling is the WOB; essentially, the load applied to the bit that allows it to cut, creating a reaction force normal to the surface. In every drilling application, there is a minimum WOB that must be maintained for the system to function properly. In microgravity (asteroids and comets), even a small WOB could not be supported conventionally by the weight of the robot or astronaut. An anchoring mechanism would be needed to resist the reactions, or the robot or astronaut would push themselves off the surface and into space. The ability of the system to anchor itself to a surface creates potential applications that reach beyond use in low gravity. The use of these anchoring mechanisms as end effectors on climbing robots has the potential of vastly expanding the scope of what is considered accessible terrain. Further, because the drill is supported by its own anchor rather than by a robotic arm, the workspace is not constrained by the reach of such an arm. Yet, if the drill is on a robotic arm, it has the benefit of not reflecting the forces of drilling back to the arm s joints. Combining the drill with the anchoring feet will create a highly mobile, highly stable, and highly reliable system. The drilling system s anchor uses hundreds of microspine toes that independently find holes and ledges on a rock to create an anchor. Once the system is anchored, a linear translation mechanism moves the drill axially into the surface while maintaining the proper WOB. The linear translation mechanism is composed of a ball screw and stepper motor that can translate a carriage with high precision and applied load. The carriage slides along rails using self-aligning linear bearings that correct any axial misalignment caused by bending and torsion. The carriage then compresses a series of springs that simultaneously transmit the load to the drill along the bit axis and act as a suspension that compensates for the vibration caused by percussive drilling. The drill is a compacted, modified version of an off-the-shelf rotary percussive drill, which uses a custom carbide-tipped coring bit. By using rotary percussive drilling, the drill time is greatly reduced. The percussive action fractures the rock debris, which is removed during rotation. The final result is a 0.75-in. (.1.9- cm) diameter hole and a preserved 0.5- in. (.1.3-cm) diameter rock core. This work extends microspine technology, making it applicable to astronaut missions to asteroids and a host of robotic sampling concepts. At the time of this reporting, it is the first instrument to be demonstrated using microspine anchors, and is the first self-contained drill/anchor system to be demonstrated that is capable of drilling in inverted configurations and would be capable of drilling in microgravity.

Parness, Aaron↗
NASA NTRS2010

RF Device for Acquiring Images of the Human Body

A safe, non-invasive method for forming images through clothing of large groups of people, in order to search for concealed weapons either made of metal or not, has been developed. A millimeter wavelength scanner designed in a unique, ring-shaped configuration can obtain a full 360 image of the body with a resolution of less than a millimeter in only a few seconds. Millimeter waves readily penetrate normal clothing, but are highly reflected by the human body and concealed objects. Millimeter wave signals are nonionizing and are harmless to human tissues when used at low power levels. The imager (see figure) consists of a thin base that supports a small-diameter vertical post about 7 ft (=2.13 m) tall. Attached to the post is a square-shaped ring 2 in. (=5 cm) wide and 3 ft (=91 cm) on a side. The ring is oriented horizontally, and is supported halfway along one side by a connection to a linear bearing on the vertical post. A planar RF circuit board is mounted to the inside of each side of the ring. Each circuit board contains an array of 30 receivers, one transmitter, and digitization electronics. Each array element has a printed-circuit patch antenna coupled to a pair of mixers by a 90 coupler. The mixers receive a reference local oscillator signal to a subharmonic of the transmitter frequency. A single local oscillator line feeds all 30 receivers on the board. The resulting MHz IF signals are amplified and carried to the edge of the board where they are demodulated and digitized. The transmitted signal is derived from the local oscillator at a frequency offset determined by a crystal oscillator. One antenna centrally located on each side of the square ring provides the source illumination power. The total transmitted power is less than 100 mW, resulting in an exposure level that is completely safe to humans. The output signals from all four circuit boards are fed via serial connection to a data processing computer. The computer processes the approximately 1-MB data set into a three-dimensional image in a matter of seconds. The innovation is to configure the receiver array in a ring topology surrounding the scanned object. The ring is then scanned vertically to cover the necessary two-dimensional surface. This fabrication of the ring is made possible by using planar antenna and circuit technology. A planar circuit board serves as a medium for both antennas and signal processing components. Using this technique, parts counts are kept low, and the cost per element is a small fraction of a waveguide-based system.

Gaier, Todd C.↗
NASA NTRS1995

Supporting technology for the development of Controlled Ecological Life Support Systems (CELSS)

To support the development of Controlled Ecological Life Support Systems (CELSS) in the space program, a metabolic simulator has been selected for use in a closed chamber to test functions of the CELSS. This metabolic simulator is a catalytic reactor which oxidizes the methyl acetate to produce carbon dioxide and water vapor. In this project, kinetic studies of catalytic oxidation of methyl acetate were conducted using monolithic and pellet catalysts with 0.5% (by weight) platinum (Pt) on aluminum oxide (Al2O3). The reaction was studied at a pressure of one atmosphere and at temperatures varying from 160 C to 420 C. By-products were identified at the exit of the preheater and reactor. For the kinetic study with the monolithic catalyst, a linear regression method was used to correlate the kinetic data with zero-order, first-order and Langmuir-Hinshelwood models. Results indicate that the first-order model represents the data adequately at low concentrations of methyl acetate. For higher concentrations of methyl acetate, the Langmuir-Hinshelwood model best represents the kinetic data. Both rate constant and adsorption equilibrium constants were estimated from the regression. A Taguchi orthogonal array (L(sub 9)) was used to investigate the effects of temperature, flow rate, and concentration on the catalytic oxidation of methyl acetate. For the monolithic catalyst, temperature exerts the most significant effect, followed by concentration of methyl acetate. For the pellet catalyst, reaction temperature is the most significant factor, followed by gas flow rate and methyl acetate concentration. Concentrations of either carbon dioxide or oxygen were seen to have insignificant effect on the methyl acetate conversion process. Experimental results indicate that the preheater with glass beads can accomplish thermal cracking and catalytic reaction of methyl acetate to produce acetic acid, methanol, methyl formate, and 1-propanol. The concentration of all by-products was measured in ppmv (parts per million by volume). At higher temperatures, greater amounts of these products are produced, as expected. In all cases, methanol was the predominant concentration detected, followed by methyl formate. At temperatures lower than 320 C for the P-type monolithic catalyst, methanol, acetic acid, and acetone were detected, whereas, for the E-type monolithic catalyst, only methanol was detected at 160 C. Both P and E types of the monolithic catalyst were specified with the same substrates (ceramic), washcoat (Al2O3), and promoter (Pt). However, the manufacturing and treatment procedures were quite different. It was therefore concluded that the performance of the E-type monolithic catalyst is superior to that of the P-type for oxidation of methyl acetate. At higher reaction temperatures, e.g., above 420 C, all reactants and byproducts were completely oxidized using these two types of monolithic catalyst to produce carbon dioxide and water vapor. A complex heterogenous catalytic reaction mechanism was proposed to explain the formation of the byproducts (methanol, acetic acid, and methyl formate) as the methyl acetate traveled through the preheater packed with glass beads. The by-product, 1-propanol, may be formed only through a homogeneous reaction, since it is difficult to develop a reasonable sequence of heterogeneous reaction steps to explain its formation. The homogeneous thermal decomposition of methyl acetate to form free radicals was proposed to explain the formation of 1-propanol, and also methanol, in the preheater. A dual-site catalytic reaction mechanism was proposed for the oxidation of methyl acetate over Pt/Al2O3 monolithic catalyst. The dual-site mechanism describes the chemisorption of oxygen molecules as well as a physical adsorption of methyl acetate on the active sites. On the active sites, methyl acetate is oxidized rapidly to form carbon dioxide and water vapor. A rate equation derived from this mechanism gives the Langmuir-Hinshelwood rate formula which has been observed from the experimental data obtained in this project for high methyl acetate concentration (greater than 1000 ppmv) over a monolithic catalyst. If the oxygen concentration is very high and methyl acetate concentration is very low, the reaction rate equation is then reduced to a first-order with respect to methyl acetate concentration. The first-order model has also been observed from the experimental data obtained in this project for low methyl acetate concentration (less than 1000 ppmv).

Li, Ku-Yen↗
NASA NTRS2018

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗
NASA NTRS2013

Robust Platinum Resistor Thermometer (PRT) Sensors and Reliable Bonding for Space Missions

Platinum resistance thermometers (PRTs) provide accurate temperature measurements over a wide temperature range and are used extensively on space missions due to their simplicity and linearity. A standard on spacecraft, PRTs are used to provide precision temperature control and vehicle health assessment. This paper reviews the extensive reliability testing of platinum resistor thermometer sensors (PRTs) and bonding methods used on the Mars Science Laboratory (MSL) mission and for the upcoming Soil Moisture Active Passive (SMAP) mission. During the Mars Exploration Rover (MER) mission, several key, JPL-packaged PRTs failed on those rovers prior to and within 1-Sol of landing due to thermally induced stresses. Similar failures can be traced back to other JPL missions dating back thirty years. As a result, MSL sought out a PRT more forgiving to the packaging configurations used at JPL, and extensively tested the Honeywell HRTS-5760-B-U-0-12 sensor to successfully demonstrate suitable robustness to thermal cycling. Specifically, this PRT was cycled 2,000 times, simulating three Martian winters and summers. The PRTs were bonded to six substrate materials (Aluminum 7050, treated Magnesium AZ231-B, Stainless Steel 304, Albemet, Titanium 6AL4V, and G-10), using four different aerospace adhesives--two epoxies and two silicones--that conformed to MSL's low out-gassing requirements. An additional epoxy was tested in a shorter environmental cycling test, when the need for a different temperature range adhesive was necessary for mobility and actuator hardware late in the fabrication process. All of this testing, along with electrostatic discharge (ESD) and destructive part analyses, demonstrate that this PRT is highly robust, and not subject to the failure of PRTs on previous missions. While there were two PRTs that failed during fabrication, to date there have been no in-flight PRT failures on MSL, including those on the Curiosity rover. Since MSL, the sensor has gone through a change in construction such that the manufacturer significantly restricts the minimum temperature. However, significant subsequent testing was performed with this new version of the part to show that it indeed is still robust to at least Mars minimum temperatures of -135 degrees Centigrade. The additional completed testing will be described. This work has resulted in a successful sensor package qualification and a reliable bonding method suitable for use over large temperature extremes

Platinum Resistor Thermometer (PRT)↗
NASA NTRS2013

Robust Platinum Resistor Thermometer (PRT) Sensors and Reliable Bonding for Space Missions

Platinum resistance thermometers (PRTs) provide accurate temperature measurements over a wide temperature range and are used extensively on space missions due to their simplicity and linearity. A standard on spacecraft, PRTs are used to provide precision temperature control and vehicle health assessment. This paper reviews the extensive reliability testing of platinum resistor thermometer sensors (PRTs) and bonding methods used on the Mars Science Laboratory (MSL) mission and for the upcoming Soil Moisture Active Passive (SMAP) mission. During the Mars Exploration Rover (MER) mission, several key, JPL-packaged PRTs failed on those rovers prior to and within 1-Sol of landing due to thermally induced stresses. Similar failures can be traced back to other JPL missions dating back thirty years. As a result, MSL sought out a PRT more forgiving to the packaging configurations used at JPL, and extensively tested the Honeywell HRTS-5760-B-U-0-12 sensor to successfully demonstrate suitable robustness to thermal cycling. Specifically, this PRT was cycled 2,000 times, simulating three Martian winters and summers. The PRTs were bonded to six substrate materials (Aluminum 7050, treated Magnesium AZ231-B, Stainless Steel 304, Albemet, Titanium 6AL4V, and G-10), using four different aerospace adhesives--two epoxies and two silicones--that conformed to MSL's low out-gassing requirements. An additional epoxy was tested in a shorter environmental cycling test, when the need for a different temperature range adhesive was necessary for mobility and actuator hardware late in the fabrication process. All of this testing, along with electrostatic discharge (ESD) and destructive part analyses, demonstrate that this PRT is highly robust, and not subject to the failure of PRTs on previous missions. While there were two PRTs that failed during fabrication, to date there have been no in-flight PRT failures on MSL, including those on the Curiosity rover. Since MSL, the sensor has gone through a change in construction such that the manufacturer significantly restricts the minimum temperature. However, significant subsequent testing was performed with this new version of the part to show that it indeed is still robust to at least Mars minimum temperatures of -135(sup o)C. The additional completed testing will be described. This work has resulted in a successful sensor package qualification and a reliable bonding method suitable for use over large temperature extremes.

dectector↗
NASA NTRS2004

Strain-Gauge Measurement of Weight of Fluid in a Tank

A method of determining the amount of fluid in a tank is based on measurement of strains induced in tank supports by the weight of the fluid. Unlike most prior methods, this method is nonintrusive: there is no need to insert instrumentation in the tank and, hence, no need to run wires, cables, or tubes through the tank wall. Also unlike most prior methods, this method is applicable even if the fluid in the tank is at supercritical pressure and temperature, because it does not depend on the presence of a liquid/gas interface (as in liquid-level-measuring methods). The strain gauges used in this method may be of two types: foil and fiber-optic. Four foil gauges (full bridge) are mounted on each of the tank-supporting legs. As the tank is filled or emptied, the deformation in each leg increases or decreases, respectively. Measured deformations of all legs are added to obtain a composite deformation indicative of the change in weight of the tank plus fluid. An initial calibration is performed by recording data at two points (usually, empty and full) for which the mass or weight of fluid is known. It is assumed that the deformations are elastic, so that the line passing through the two points can be used as a calibration curve of mass (or weight) of fluid versus deformation. One or more fiber-optic gauges may be used instead of the foil gauges. The resolution of the fiber-optic and foil gauges is approximately the same, but the fiber-optic gauges are immune to EMI (electromagnetic interference), are linear with respect to temperature over their entire dynamic range (as defined by the behavior of the sample), and measure thermally induced deformations as predictable signals. Conversely, long term testing has demonstrated that the foil gauges exhibit an erratic behavior whenever subjected to direct sun radiation (even if protected with a rubberized cover). Henceforth, for deployment in outdoor conditions, fiber-optic gauges are the only option if one is to rely on the system for an extended period of time when a recalibration procedure may not be acceptable. A set of foil gauges had been tested on the supports of a 500-gallon (1,900-liter) tank. The gauges were found to be capable of measuring the deformations (up to 22 micro-strain) that occurred during filling and emptying of the tank. The fluid masses calculated from the gauge readings were found to be accurate within 4.5 percent. However, the reliability of the foil gauges over a few hours was not acceptable. Therefore, the foil sensor system is acceptable for use only in controlled environments (complete shade, or indoors).

Figueroa, Jorge↗
NASA NTRS1994

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗