Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface shape fitting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

CRUQS: A Miniature Fine Sun Sensor for Nanosatellites

A new miniature fine Sun sensor has been developed that uses a quadrant photodiode and housing to determine the Sun vector. Its size, mass, and power make it especially suited to small satellite applications, especially nanosatellites. Its accuracy is on the order of one arcminute, and it will enable new science in the area of nanosatellites. The motivation for this innovation was the need for high-performance Sun sensors in the nanosatellite category. The design idea comes out of the LISS (Lockheed Intermediate Sun Sensor) used by the sounding rocket program on their solar pointing ACS (Attitude Control System). This system uses photodiodes and a wall between them. The shadow cast by the Sun is used to determine the Sun angle. The new sensor takes this concept and miniaturizes it. A cruciform shaped housing and a surface-mount quadrant photodiode package allow for a two-axis fine Sun sensor to be packaged into a space approx.1.25xl x0.25 in. (approx.3.2x2.5x0.6 cm). The circuitry to read the photodiodes is a simple trans-impedance operational amplifier. This is much less complex than current small Sun sensors for nanosatellites that rely on photo-arrays and processing of images to determine the Sun center. The simplicity of the circuit allows for a low power draw as well. The sensor consists of housing with a cruciform machined in it. The cruciform walls are 0.5-mm thick and the center of the cruciform is situated over the center of the quadrant photodiode sensor. This allows for shadows to be cast on each of the four photodiodes based on the angle of the Sun. A simple operational amplifier circuit is used to read the output of the photodiodes as a voltage. The voltage output of each photodiode is summed based on rows and columns, and then the values of both rows or both columns are differenced and divided by the sum of the voltages for all four photodiodes. The value of both difference over sums for the rows and columns is compared to a table or a polynomial fit (depending on processor power and accuracy requirements) to determine the angle of the Sun in the sensor frame.

Heatwole, Scott↗
NASA NTRS1993

Interferometric detection of freeze-thaw displacements of Alaskan permafrost using ERS-1 data

The possibility of making large scale (50 km) measurements of motions of the earth's surface with high resolution (10 m) and very high accuracy (1 cm) from multipass SAR interferometry was established in 1989. Other experiments have confirmed the viability and usefulness of the method. Work is underway in various groups to measure displacements from volcanic activity, seismic events, glacier motion, and in the present study, freeze-thaw cycles in Alaskan permafrost. The ground is known to move significantly in these cycles, and provided that freezing does not cause image decorrelation, it should be possible to measure both ground swelling and subsidence. The authors have obtained data from multiple passes of ERS-1 over the Toolik Lake region of northern Alaska of suitable quality for interferometry. The data are processed into images, and single interferograms are formed in the usual manner. Phase unwrapping is performed, and the multipass baselines are estimated from the images using both orbit ephemerides and scene tie points. The phases are scaled by the baseline ratio, and a double-difference interferogram (DDI) is formed. It is found that there is a residual 'saddle-shape' phase error across the image, which is postulated to be caused by a small divergence (10(exp -2) deg.) in the orbits. A simulation of a DDI from divergent orbits confirms the shape and magnitude of the error. A two-dimensional least squares fit to the error is performed, which is used to correct the DDI. The final, corrected DDI shows significant phase (altitude) changes over the period of the observation.

Werner, Charles L.↗
DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗
NASA NTRS2002

Compositionally-Graded Shape Memory Film for Self-Deployment of Membrane Reflectors and Optics

The next generation of space systems will require large apertures in order to image faint targets or cover large areas of Earth. These large apertures must be able to fit inside a launch vehicle fairing, be light enough for launch into orbit, and deploy on orbit with repeatability and reliability. The current state-of-the-art in flight optics is represented by the 4 meter LAMP telescope, with an areal density of 10 km sq m. Development of a Beryllium mirror demonstration article for NGST (Next Generation Space Telescope) at the University of Arizona indicate areal densities of 0.5 kg sq m with flight hardware in the 12 meter range. With progressive improvements in existing deployment, packaging, and structural technologies, the size of optics and reflectors will continue to increase, while mass is reduced. However, without a breakthrough in materials, packaging and/or deployment technologies, the goal for Gossamer structures of 0.1 kg sq m is unachievable for the near and mid-term NASA missions. Membrane technology provides the best hope of achieving such low areal densities. In combination with advances in membrane materials and structures, development of revolutionary techniques for deployment systems can provide significant improvements in large aperture technology. In this paper, the results of a six-month Phase I research effort to demonstrate the application of thin film NiTi to aerospace-qualified membrane and mesh materials are presented. Deposition of shape memory thin film was achieved Astromesh (trademark) metal mesh and CP-1, and optical-quality polymer membrane. Not only was full-coating deposition demonstrated, but also small segment deposition which holds potential for local surface control. Deployment of these materials was also demonstrated, setting the stage for the development of a larger test article.

Hill, Lisa↗
NASA NTRS1994

Excitation of solar p-modes

We investigate the rates at which energy is supplied to individual p-modes as a function of their frequencies nu and angular degrees l. The observationally determined rates are compared with those calculated on the hypothesis that the modes are stochastically excited by turbulent convection. The observationally determined excitation rate is assumed to be equal to the product of the mode's energy E and its (radian) line width Gamma. We obtain E from the mode's mean square surface velocity with the aid of its velocity eigenfuction. We assume that Gamma measures the mode's energy decay rate, even though quasi-elastic scattering may dominate true absorption. At fixed l, E(Gamma) arises as nu(exp 7) at low nu, reaches a peak at nu approximately equal 3.5 mHz, and then declines as nu(exp 4.4) at higher nu . At fixed nu, E(Gamma) exhibits a slow decline with increasing l. To calculate energy input rates, P(sub alpha), we rely on the mixing-length model of turbulent convection. We find entropy fluctuations to be about an order of magnitude more effective than the Reynolds stress in exciting p-modes . The calculated P(sub alpha) mimic the nu(exp 7) dependence of E(Gamma) at low nu and the nu(exp -4.4) dependence at high nu. The break of 11.4 powers in the nu-dependence of E(Gamma) across its peak is attributed to a combination of (1) the reflection of high-frequency acoustic waves just below the photosphere where the scale height drops precipitously and (2) the absence of energy-bearing eddies with short enough correlation times to excite high-frequency modes. Two parameters associated with the eddy correlation time are required to match the location and shape of the break. The appropriate values of these parameters, while not unnatural, are poorly constrained by theory. The calculated P(sub alpha) can also be made to fit the magnitude of E(Gamma) with a reasonable value for the eddy aspect ratio. Our resutls suggest a possible explanation for the decline of mode energy with increasing l at fixed nu. Entropy fluctuations couple to changes in volume associated with the oscillation mode. These decrease with decreasing n at fixed nu, becoming almost zero for the f-mode.

Goldreich, Peter↗
NASA NTRS2012

Photogrammetric and Global Positioning System Measurements of Active Pahoehoe Lava Lobe Emplacement on Kilauea, Hawaii

Basalt is the most common rock type on the surface of terrestrial bodies throughout the solar system and -- by total volume and areal coverage -- pahoehoe flows are the most abundant form of basaltic lava in subaerial and submarine environments on Earth. A detailed understanding of pahoehoe emplacement processes is necessary for developing accurate models of flow field development, assessing hazards associated with active lava flows, and interpreting the significance of lava flow morphology on Earth and other planetary bodies. Here, we examine the active emplacement of pahoehoe lobes along the margins of the Hook Flow from Pu'u 'O'o on Kilauea, Hawaii. Topographic data were acquired between 21 and 23 February 2006 using stereo-imaging and differential global positing system (DGPS) measurements. During this time, the average discharge rate for the Hook Flow was 0.01-0.05 cubic m/s. Using stereogrammetric point clouds and interpolated digital terrain models (DTMs), active flow fronts were digitized at 1 minute intervals. These areal spreading maps show that the lava lobe grew by a series of breakouts tha t broadly fit into two categories: narrow (0.2-0.6 m-wide) toes that grew preferentially down-slope, and broad (1.4-3.5 m-wide) breakouts that formed along the sides of the lobe, nearly perpendicular to the down-flow axis. These lobes inflated to half of their final thickness within approx 5 minutes, with a rate of inflation that generally deceased with time. Through a combination of down-slope and cross-slope breakouts, lobes developed a parabolic cross-sectional shape within tens of minutes. We also observed that while the average local discharge rate for the lobe was generally constant at 0.0064 +/- 0.0019 cubic m/s, there was a 2 to 6 fold increase in the areal coverage rate every 4.1 +/- 0.6 minutes. We attribute this periodicity to the time required for the dynamic pressurization of the liquid core of the lava lobe to exceed the cooling-induced strength of the lobe margins. Using DGPS-derived DTMs of the topography before and after pahoehoe lobe emplacement, we observed that the lava typically concentrated within existing topographic lows, with the lobe reaching a maximum thickness of approx 1.2 m above the lowest points of the initial topography and above reverse-facing slopes. Lobe margins were typically controlled by high-standing topography, with the zone directly adjacent to the final flow margin having average relief that is approx 4 cm higher than the lava-inundated region. This suggests that irregularities approx 25% of the height of the smallest breakout elements (i.e., toes) can exert a strong control on the paths of low-discharge pahoehoe lobes, with stagnated toes forming confining margins that allow interior portions of flow to topographically invert the landscape by inflation.

Hamilton, Christopher W.↗
NASA NTRS2002

The Evolution of Disks and Winds in Dwarf Nova Outbursts - FUSE

This project was a project to study the FUV spectra of two proto-typical dwarf novae, U Gem and SS Cygni, through an outburst cycle. The luminosity of the boundary layer in the two systems, as evidenced by earlier EUVE observations, is different in the two systems. Our intensive study of the two systems was intended to (1) probe the ionization and kinematic structure of the wind as a function of system brightness, (2) isolate the contributions of the disk to the FUV spectra, and (3) examine physical conditions and abundances of material just being accreted onto the disk from the secondary. The U Gem and SS Cyg observations took place in March and October 2000, respectively. The data obtained with FUSE was of excellent quality. Analysis of the both observations is now essentially complete, although some modeling of the SS Cyg spectra is ongoing, as we complete an ApJ manuscript on this object. Our main results for U Gem are as follows: The plateau spectra have continuum shapes and fluxes that are approximated by steady state accretion disk model spectra with an accretion rate 7x10(exp 9) Msolar/yr. The spectra also show numerous absorption lines of H I, He II, and 2-5 times ionized transitions of C, N, O, P, S, and Si. There are no emission features in the spectra, with the possible exception of a weak feature on the red wing of the 0 VI doublet. The absorption lines are narrow (FWHM approx. 50 km/s), too narrow to arise from the disk photosphere, and at low velocities (less than or equal to 700 km/s). The S VI and O VI doublets are optically thick. The absorption lines in the plateau spectra show orbital variability: in spectra obtained at orbital phases between 0.53 and 0.79, low-ionization absorption lines appear and the central depths of the preexisting lines increase. The increase in line absorption occurs at the same orbital phases as previously observed EUV and X-ray light-curve dips. If the absorbing material is in (near-) Keplerian rotation around the disk, it must be located at large disk radii. The final observation occurred when U Gem was about 2 mag from optical quiescence. The spectra are dominated by emission from an approx. 43,000 K, metal-enriched white dwarf (WD). The inferred radius of the WD is 4.95x10(exp 8) cm, close to that observed in quiescence. Allowing for a hot heated region on the surface of the WD improves the fit to the spectrum at short (less than 960 A) wavelengths. Our main results for SS Cyg are as follows: The first two of four observations of SS Cyg show disk dominated spectra with accretion rates of order 10(exp -8) Msolar/yr. Except for narrow interstellar features (atomic and molecular H), the lines are all broad consistent with a disk or wind origin. The O VI line in the spectra is mostly of wind origin as detailed modeling with our Monte Carlo code (developed in part using funds from this project) show. The continua from spectra in observations 3 and 4, observed during the decline phase, are not well fit with steady-state disks, and show considerable resemblance to quiescent spectra obtained with HUT. The most probable interpretation for the emission features seen in the spectrum in the last two observations is that they arise from a photo-illuminated choronosphere above the disk, rather than a wind.

Long, Knox↗
NASA NTRS2012

Metrology of Large Parts

As discussed in the first chapter of this book, there are many different methods to measure a part using optical technology. Chapter 2 discussed the use of machine vision to measure macroscopic features such as length and position, which was extended to the use of interferometry as a linear measurement tool in chapter 3, and laser or other trackers to find the relation of key points on large parts in chapter 4. This chapter looks at measuring large parts to optical tolerances in the sub-micron range using interferometry, ranging, and optical tools discussed in the previous chapters. The purpose of this chapter is not to discuss specific metrology tools (such as interferometers or gauges), but to describe a systems engineering approach to testing large parts. Issues such as material warpage and temperature drifts that may be insignificant when measuring a part to micron levels under a microscope, as will be discussed in later chapters, can prove to be very important when making the same measurement over a larger part. In this chapter, we will define a set of guiding principles for successfully overcoming these challenges and illustrate the application of these principles with real world examples. While these examples are drawn from specific large optical testing applications, they inform the problems associated with testing any large part to optical tolerances. Manufacturing today relies on micrometer level part performance. Fields such as energy and transportation are demanding higher tolerances to provide increased efficiencies and fuel savings. By looking at how the optics industry approaches sub-micrometer metrology, one can gain a better understanding of the metrology challenges for any larger part specified to micrometer tolerances. Testing large parts, whether optical components or precision structures, to optical tolerances is just like testing small parts, only harder. Identical with what one does for small parts, a metrologist tests large parts and optics in particular to quantify their mechanical properties (such as dimensions, mass, etc); their optical prescription or design (i.e. radius of curvature, conic constant, vertex location, size); and their full part shape. And, just as with small parts, a metrologist accomplishes these tests using distance measuring instruments such as tape measures, inside micrometers, coordinate measuring machines, distance measuring interferometers; angle measuring instruments such as theodolites, autocollimators; and surface measuring instruments including interferometers, stylus profilers, interference microscopes, photogrammetric cameras, or other tools. However, while the methodology may be similar, it is more difficult to test a large object for the simple reason that most metrologists do not have the necessary intuition. The skills used to test small parts or optics in a laboratory do not extrapolate to testing large parts in an industrial setting any more than a backyard gardener might successfully operate a farm. But first, what is a large part? A simple definition might be the part's size or diameter. For optics and diffuse surface parts alike, the driving constraint is ability to illuminate the part's surface. For reflective convex mirrors, large is typically anything greater than 1 meter. But, for refractive optics, flats or convex mirrors, large is typically greater than 0.5 meter. While a size definition is simple, it may be less than universal. A more nuanced definition might be that a large part is any component which cannot be easily tested in a standard laboratory environment, on a standard vibration isolated table using standard laboratory infrastructure. A micro-switch or a precision lens might be easily measured to nanometer levels under a microscope in a lab, but a power turbine spline or a larger telescope mirror will not fit under that microscope and may not even fit on the table.

Stahl, H. Philip↗
NASA NTRS2007

Enhanced Elliptic Grid Generation

An enhanced method of elliptic grid generation has been invented. Whereas prior methods require user input of certain grid parameters, this method provides for these parameters to be determined automatically. "Elliptic grid generation" signifies generation of generalized curvilinear coordinate grids through solution of elliptic partial differential equations (PDEs). Usually, such grids are fitted to bounding bodies and used in numerical solution of other PDEs like those of fluid flow, heat flow, and electromagnetics. Such a grid is smooth and has continuous first and second derivatives (and possibly also continuous higher-order derivatives), grid lines are appropriately stretched or clustered, and grid lines are orthogonal or nearly so over most of the grid domain. The source terms in the grid-generating PDEs (hereafter called "defining" PDEs) make it possible for the grid to satisfy requirements for clustering and orthogonality properties in the vicinity of specific surfaces in three dimensions or in the vicinity of specific lines in two dimensions. The grid parameters in question are decay parameters that appear in the source terms of the inhomogeneous defining PDEs. The decay parameters are characteristic lengths in exponential- decay factors that express how the influences of the boundaries decrease with distance from the boundaries. These terms govern the rates at which distance between adjacent grid lines change with distance from nearby boundaries. Heretofore, users have arbitrarily specified decay parameters. However, the characteristic lengths are coupled with the strengths of the source terms, such that arbitrary specification could lead to conflicts among parameter values. Moreover, the manual insertion of decay parameters is cumbersome for static grids and infeasible for dynamically changing grids. In the present method, manual insertion and user specification of decay parameters are neither required nor allowed. Instead, the decay parameters are determined automatically as part of the solution of the defining PDEs. Depending on the shape of the boundary segments and the physical nature of the problem to be solved on the grid, the solution of the defining PDEs may provide for rates of decay to vary along and among the boundary segments and may lend itself to interpretation in terms of one or more physical quantities associated with the problem.

Kaul, Upender K.↗
NASA NTRS2016

Unsteady Aerodynamic Force Sensing from Measured Strain

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm. A cantilevered rectangular wing built and tested at the NASA Langley Research Center (Hampton, Virginia, USA) in 1959 is used to validate the simple approach. Unsteady aerodynamic forces as well as wing deflections, velocities, accelerations, and strains are computed using the CFL3D computational fluid dynamics (CFD) code and an MSC/NASTRAN code (MSC Software Corporation, Newport Beach, California, USA), and these CFL3D-based results are assumed as measured quantities. Based on the measured strains, wing deflections, velocities, accelerations, and aerodynamic forces are computed using the proposed approach. These computed deflections, velocities, accelerations, and unsteady aerodynamic forces are compared with the CFL3D/NASTRAN-based results. In general, computed aerodynamic forces based on the lifting surface theory in subsonic speeds are in good agreement with the target aerodynamic forces generated using CFL3D code with the Euler equation. Excellent aeroelastic responses are obtained even with unsteady strain data under the signal to noise ratio of -9.8dB. The deflections, velocities, and accelerations at each sensor location are independent of structural and aerodynamic models. Therefore, the distributed strain data together with the current proposed approaches can be used as distributed deflection, velocity, and acceleration sensors. This research demonstrates the feasibility of obtaining induced drag and lift forces through the use of distributed sensor technology with measured strain data. An active induced drag control system thus can be designed using the two computed aerodynamic forces, induced drag and lift, to improve the fuel efficiency of an aircraft. Interpolation elements between structural finite element grids and the CFD grids and centroids are successfully incorporated with the unsteady aeroelastic computation scheme. The most critical technology for the success of the proposed approach is the robust on-line parameter estimator, since the least-squares curve fitting method depends heavily on aeroelastic system frequencies and damping factors.

shape sensing↗
NASA NTRS2010

Fast-Response-Time Shape-Memory-Effect Foam Actuators

Bulk shape memory alloys, such as Nitinol or CuAlZn, display strong recovery forces undergoing a phase transformation after being strained in their martensitic state. These recovery forces are used for actuation. As the phase transformation is thermally driven, the response time of the actuation can be slow, as the heat must be passively inserted or removed from the alloy. Shape memory alloy TiNi torque tubes have been investigated for at least 20 years and have demonstrated high actuation forces [3,000 in.-lb (approximately equal to 340 N-m) torques] and are very lightweight. However, they are not easy to attach to existing structures. Adhesives will fail in shear at low-torque loads and the TiNi is not weldable, so that mechanical crimp fits have been generally used. These are not reliable, especially in vibratory environments. The TiNi is also slow to heat up, as it can only be heated indirectly using heater and cooling must be done passively. This has restricted their use to on-off actuators where cycle times of approximately one minute is acceptable. Self-propagating high-temperature synthesis (SHS) has been used in the past to make porous TiNi metal foams. Shape Change Technologies has been able to train SHS derived TiNi to exhibit the shape memory effect. As it is an open-celled material, fast response times were observed when the material was heated using hot and cold fluids. A methodology was developed to make the open-celled porous TiNi foams as a tube with integrated hexagonal ends, which then becomes a torsional actuator with fast response times. Under processing developed independently, researchers were able to verify torques of 84 in.-lb (approximately equal to 9.5 Nm) using an actuator weighing 1.3 oz (approximately equal to 37 g) with very fast (less than 1/16th of a second) initial response times when hot and cold fluids were used to facilitate heat transfer. Integrated structural connections were added as part of the net shape process, eliminating the need for welding, adhesives, or mechanical crimping. Inexpensive net-shape processing was used, which reduces the cost of the actuator by over a factor of 10 over nonporous TiNi made by hot drawing of tube or electrical discharge machining. By forming the alloy as an open-celled foam, the surface area for heat transfer is dramatically increased, allowing for much faster response times. The technology also allows for netshape fabrication of the actuator, which allows for structural connections to be integrated into the actuator material, making these actuators significantly less expensive. Commercial applications include actuators for concepts such as the variable area chevron and nozzle in jet aircraft. Lightweight tube or rod components can be supplied to interested parties.

Jardine, Peter↗
NASA NTRS1995

A new model for Mars atmospheric dust based upon analysis of ultraviolet through infrared observations from Mariner 9, Viking, and Phobos

We propose key modifications to the Toon et al. (1977) model of the particle size distribution and composition of Mars atmospheric dust, based on a variety of spacecraft and wavelength observations of the dust. A much broader (r(sub eff)variance-0.8 micron), smaller particle size (r(sub mode)-0.02 microns) distribution coupled with a "palagonite-like" composition is argued to fit the complete ultraviolet-to-30-micron absorption properties of the dust better than the montmorillonite-basalt r(sub eff)variance= 0.4 micron, r(sub mode)= 0.40 micron dust model of Toon et al. Mariner 9 (infrared interferometer spectrometer) IRIS spectra of high atmospheric dust opacities during the 1971 - 1972 Mars global dust storm are analyzed in terms of the Toon et al. dust model, and a Hawaiian palagonite sample with two different size distribution models incorporating smaller dust particle sizes. Viking Infrared Thermal Mapper (IRTM) emission-phase-function (EPF) observations at 9 microns are analyzed to retrieve 9-micron dust opacities coincident with solar band dust opacities obtained from the same EPF sequences. These EPF dust opacities provide an independent measurement of the visible/9-microns extinction opacity ratio (> or equal to 2) for Mars atmospheric dust, which is consistent with a previous measurement by Martin (1986). Model values for the visible/9-microns opacity ratio and the ultraviolet and visible single-scattering albedos are calculated for the palagonite model with the smaller particle size distributions and compared to the same properties for the Toon et al. model of dust. The montmorillonite model of the dust is found to fit the detailed shape of the dust 9-micron absorption well. However, it predicts structured, deep absorptions at 20 microns which are not observed and requires a separate ultraviolet-visible absorbing component to match the observed behavior of the dust in this wavelength region. The modeled palagonite does not match the 8- to 9-micron absorption presented by the dust in the IRIS spectra, probably due to its low SiO2 content (31%). However, it does provide consistent levels of ultraviolet/visible absorption, 9- to 12-micron absorption, and a lack of structured absorption at 20 microns. The ratios of dust extinction opacities at visible, 9 microns, and 30 microns are strongly affected by the dust particle size distribution. The Toon et al. dust size distribution (r(sub mode)= 0.40, r(sub eff)variance= 0.4 microns, r(sub cw mu)= 2.7 microns) predicts the correct ratio of the 9- to 30-micron opacity, but underpredicts the visible/9-micron opacity ratio considerably (1 versus > or equal to 2). A similar particle distribution width with smaller particle sizes (r(sub mode)= 0.17, r(sub eff)variance= 0.4 microns, r(sub cw mu)=1.2 microns) will fit the observed visible/9-micron opacity ratio, but overpredicts the observed 9-micron/30-micron opacity ratio. A smaller and much broader particle size distribution (r(sub mode)= 0.02, r(sub eff)variance= 0.8 microns, r(sub cw mu)= 1.8 microns) can fit both dust opacity ratios. Overall, the nanocrystalline structure of palagonite coupled with a smaller, broader distribution of dust particle sizes provides a more consistent fit than the Toon et al. model of the dust to the IRIS spectra, the observed visible/9-micron dust opacity ratio, the Phobos occultation measurements of dust particle sizes, and the weakness of surface near IR absorptions expected for clay minerals.

Clancy, R. T.↗
NASA NTRS2004

Evidence for exclusively inorganic formation of magnetite in Martian meteorite ALH84001

Magnetite crystals produced by terrestrial magnetotactic bacterium MV-1 are elongated on a [111] crystallographic axis, in a so-called truncated hexa-Octahedral shape. This morphology has been proposed to constitute a biomarker (i.e., formed only in biogenic processes). A subpopulation of magnetite crystals associated with carbonate globules in Martian meteorite ALH84001 is reported to have this morphology, and the observation has been taken as evidence for biological activity on Mars. In this study, we present evidence for the exclusively inorganic origin of [111]-elongated magnetite crystals in ALH84001. We report three-dimensional(3-D) morphologies for approx.1000 magnetite crystals extracted from: (1) thermal decomposition products of Fe-rich carbonate produced by inorganic hydrothermal precipitation in laboratory experiments; (2) carbonate globules in Martian meteoriteeALH84001; and (3) cells of magnetotactic bacterial strain MV-1. The 3-D morphologies were derived by fitting 3-D shape models to two-dimensional bright-field transmission-electron microscope (TEAM) images obtained at a series of viewing angles. The view down the {110} axes closest to the [111] elongation axis of magnetite crystals ([111]x{110) not equal to 0) provides a 2-D projection that uniquely discriminates among the three [111]-elongated magnetite morphologies found in these samples: [111]-elongated truncated hexaoctahedron ([111]-THO), [111]-elongated cubo-octahedron ([111]-ECO), and [111]-elongated simple octahedron ([111]-ESO). All [111] -elongated morphologies are present in the three types of sample, but in different proportions. In the ALH84001 Martian meteorite and in our inorganic laboratory products, the most common [111]-elongated magnetite crystal morphology is [111]-ECO. In contrast, the most common morphology for magnetotactic bacterial strain MV-1 is [111]-THO. These results show that: (1) the morphology of [111]-elongated magnetite crystals associated with the carbonate globules in Martian meteorite ALH84001 is replicated by an inorganic process; and (2) the most common crystal morphology for biogenic (MV-1) magnetite is distinctly different from that in both ALH84001 and our inorganic laboratory products. Therefore, [111]-elongated magnetite crystals in ALH84001 do not constitute, as previously claimed, a robust biosignature and, in fact, an exclusively inorganic origin for the magnetite is fully consistent with our results. Furthermore, the inorganic synthesis method, i.e., the thermal decomposition of hydrothermally precipitated Fe-rich carbonate, is a process analogue for formation of the magnetite on Mars. Namely, precipitation of carbonate globules from carbonate-rich hydrothermal solutions followed at some later time by a thermal pulse, perhaps in association with meteoritic impact or volcanic processes on the Martian surface.

Golden, D. C.↗
NASA NTRS1995

A new model for Mars atmospheric dust based upon analysis of ultraviolet through infrared observations from Mariner 9, Viking, and Phobos

We propose key modifications to the Toon et al. (1977) model of the particle size distribution and composition of Mars atmospheric dust, based on a variety of spacecraft and wavelength observations of the dust. A much broader (r(sub eff) variance approximately 0.8 micrometers), smaller particle size (r(sub mode) approximately 0.02 micrometers) distribution coupled with a 'palagonite-like' composition is argued to fit the complete ultraviolet-to-30-micrometer absorption properties of the dust better than the montmorillonite-basalt, r(sub eff) variance = 0.4 micrometers, r(sub mode) = 0.40 dust model of Toon et al. Mariner 9 (infrared interferometer spectrometer) IRIS spectra of high atmospheric dust opacities during the 1971-1972 Mars global dust storm are analyzed in terms of the Toon et al. dust model, and a Hawaiian palagonite sample (Rousch et al., 1991) with two different size distribution models incorporating smaller dust particle sizes. Viking Infrared Thermal Mapper (IRTM) emmission-phase-function (EPF) observations at 9 micrometers are analyzed to retrieve 9-micrometer dust opacities coincident with solar band dust opacities obtained from the same EPF sequences (Clancy and Lee, 1991). These EPF dust opacities provide an independent measurement of the visible/9-micrometer extinction opacity ratio (greater than or = 2) for Mars atmospheric dust, which is consistent with a previous measurement by Martin (1986). Model values for the visible/9-micrometer opacity ratio and the ultraviolet and visible single-scattering albedos are calculated for the palagonite model with the smaller particle size distributions compared to the same properties for the Toon et al. model of dust. The montmorillonite model of the dust is found to fit the detailed shape of the dust 9-micrometer absorption well. However, it predicts structured, deep aborptions at 20 micrometers which are not observed and requires a separate ultraviolet-visible absorbing component to match the observed behavior of the dust in this wavelength region. The modeled palagonite does not match the 8-to 9-micrometer absorption presented by the dust in the IRIS spectra, probably due to its low SiO2 content (31%). However, it does provide consistent levels of ultraviolet/visible absorption, 9-to 12-micrometer absorption, and a lack of structured absorption at 20 micrometers. The ratios of dust extinction opacities at visible, 9 micrometers, and 30 micrometers are strongly affected by the dust particle size distribution. The Toon et al. dust size distribution (r(sub mode) = 0.40,r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 2.7 micrometers) predicts the correct ratio of the 9- to 30-micrometer opacity, but underpredicts the visible/9-micrometer opacity ratio considerably (1 versus greater than or = 2). A similar particle distribution width with smaller particle sizes (r(sub mode) = 0.17, r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 1.2 micrometers) will fit the observed visible/9-micrometer opacity ratio, but overpredicts the observed 9-micrometer/30-micrometer opacity ratio. A smaller and much broader particle size distribution (r(sub mode) = 0.002, r(sub eff) variance = 0.8 micrometers, r(sub cwmu) = 1.8 micrometers) can fit both dust opacity ratios. Overall, the nanocrystalline structure of palagonite coupled with a smaller, broader distribution of dust particle sizes provides a more consistent fit than the Toon et al. model of the dust to the IRIS spectra, the observed visible/9-micrometer dust opacity ratio, the Phobos occulation measurements of the dust particle sizes (Chassefiere et al., 1992), and the weakness of surface near IR absorptions expected for clay minerals (Clark, 1992; Bell and Crisp, 1993).

Clancy, R. T.↗
NASA NTRSDate not supplied

Dynamic Stability Analysis Using Free-Flight CFD

Blunt body probes are the primary design for atmospheric entry. The blunted capsule shape, while desirable for the decreased forebody heating, typically experiences dynamic instability at low supersonic and transonic flight. At moderate to low supersonic conditions, the unsteady wake begins to adversely interact with the vehicle in free-flight resulting in attitude oscillations which can grow in time and result in tumbling. Understanding vehicle dynamics is thus important when considering vehicle entry design and guidance navigation and control (GN&C) considerations. Historically, dynamic stability has been carried out exclusively using experimental techniques. The wide range of experimental approaches provide reasonable coverage in Mach and dynamic pressure space for a given trajectory, however achieving flight similitude is often difficult or impossible. Advancements in Computational Fluid Dynamics (CFD) techniques allow for dynamic motion of the vehicle surface during the simulation with either forced/prescribed motion or free motion in response to fluid forces. Several examples of CFD software with the ability to simulate fluid dynamics and vehicle motion are available in literature. The Free-Flight CFD (FF-CFD) capability developed by NASA Ames within the US3D flow solver has shown extensive validation and verification (V&V) efforts from moderate to low supersonic flow across a range of vehicle architectures. The solver is able to simulate full rotational and translation motion resulting in full 6 degrees-of-freedom (DoF) motion. The validation and verification efforts the solver's ability to predict the motion of a ballistic range model as compared to experimental results. Similar to other computational capabilities in literature, FF-CFD relied on heritage data reduction methodologies and processes for reducing dynamic data into aerodynamic coefficients intended for use in design and trajectory codes. These heritage approaches, initially developed with the assumptions (and restrictions) of ground test facilities, resulted in the significant culling of the full 6-DoF dynamic data set to discrete values which are then fit using an aerodynamic model. Recent development of simulation and data reduction methodologies has resulted in an improved approach to FF-CFD simulations of vehicle dynamics which leverage the control of dynamic constraints within a simulation set-up. The presentation of this work will cover the constrained dynamic simulation approach for deriving dynamic coefficients for a Genesis capsule. The methodology will be outlined and a discussion of implications in motion restriction will be discussed. Finally, reconstructed vehicle trajectories using FF-CFD derived aerodynamic coefficients will be presented.

Joseph M Brock↗
NASA NTRSDate not supplied

The Instrumented Walking and Turning Test to Evaluate Suited Gait Dynamics and Performance in Extravehicular Activity Training Environments

Background and aims: Walking will be required for many exploration tasks on the Moon during the Artemis program. Walking in a straight line on the confined floorspace of a testing area, and repetitive treadmill walking that requires no change in direction may not adequately reflect the balance and coordination required during ambulation. Also, performance of turning maneuvers may be affected differently in different extravehicular activity (EVA) training facilities that simulate partial gravity. For example, the Neutral Buoyancy Lab (NBL) simulates lunar gravity by adding weight to underwater subjects to alter buoyancy and achieve the equivalent ground reaction force of 1/6 of Earth’s gravity (1/6G), whereas the Active Response Gravity Offload System (ARGOS) uses a computer controlled overhead suspension system programmed to continuously offload a percentage of a subject’s weight to simulate 1/6G. The degree to which dynamic movements such as turning are comparable across these EVA training facilities has not yet been evaluated. The instrumented gait test helps NASA scientists and engineers evaluate gait dynamics and performance in suited conditions, and this test demonstrates the unique characteristics and limitations of EVA training facilities. We developed an instrumented walking and turning test using inertial measurement units (IMUs) and conducted the test at NASA’s EVA training facilities. Results were used to compare suited walking and turning characteristics in the ARGOS and the NBL. Methods: Subjects donned the Mark III space suit during offloading with the ARGOS spreader bar gimbal and donned the Z2.5 space suit while underwater in the NBL with weights and floatation added to achieve realistic suit center of gravity. The test team securely attached three Opal (APDM, OR, USA) wireless IMUs on the space suit for each test run: one on the middle of the hard upper torso, and one on the left and on the right ankle bearings. During the NBL tests, the IMUs were encased in a waterproof housing (GoPro) with foam added to create a tighter fit. At both testing facilities, 6.3 m x 1.0 m (LxW) walking lines were marked, and a cone for turning or walking around was located at the end of the walking path with another line on the other side of the cone to indicate the stopping point after walking around the cone. Under simulated 1/6G, subjects began by standing at the marked line with their arms folded across the chest, they then walked at a preferred speed along the straight walking path until they reached the end, turned 180 degrees around the cone, and finally stopped at the marked stopping point. All IMU data recorded during testing were automatically saved to the internal memory. Then, raw IMU signals were processed using custom MATLAB (Mathworks, MA, USA) code to compare gait parameters during both the walking and the turning components of the task. These parameters included time (s), speed (m/s for walking and rad/s for turning), step number (n) and walk:turn time ratio (% time spent straight walking versus turning). Results: Less time, faster gait, fewer steps, and higher walk:turn ratio during both walking and turning components were exhibited during tests performed at the ARGOS versus those performed at the NBL. During the NBL tests, the slower walking speed continued at the same rate throughout a U-shape turn. During the ARGOS tests, the subjects performed shorter and tighter turns at 4 times the speed of the NBL turns because they walked 30% faster and the vertical offloading system gave them more support. Conclusion: Our data show that the differences in walking and turning parameters during the NBL tests may be due to the high viscosity in the water environment where the motion of the lower limbs was slow and did not reach full flexion and extension. These tests improve the current knowledge of testing environments in preparation for EVAs on the lunar surface.

Kyoung Jae Kim↗
NASA NTRS2002

A Proposed Model for Protein Crystal Nucleation and Growth

How does one take a molecule, strongly asymmetric in both shape and charge distribution, and assemble it into a crystal? We propose a model for the nucleation and crystal growth process for tetragonal lysozyme, based upon fluorescence, light, neutron, and X-ray scattering data, size exclusion chromatography experiments, dialysis kinetics, AFM, and modeling of growth rate data, from this and other laboratories. The first species formed is postulated to be a 'head to side' dimer. Through repeating associations involving the same intermolecular interactions this grows to a 4(sub 3) helix structure, that in turn serves as the basic unit for nucleation and subsequent crystal growth. High salt attenuates surface charges while promoting hydrophobic interactions. Symmetry facilitates subsequent helix-helix self-association. Assembly stability is enhanced when a four helix structure is obtained, with each bound to two neighbors. Only two unique interactions are required. The first are those for helix formation, where the dominant interaction is the intermolecular bridging anion. The second is the anti-parallel side-by-side helix-helix interaction, guided by alternating pairs of symmetry related salt bridges along each side. At this stage all eight unique positions of the P4(sub3)2(sub 1),2(sub 1) unit cell are filled. The process is one of a) attenuating the most strongly interacting groups, such that b) the molecules begin to self-associate in defined patterns, so that c) symmetry is obtained, which d) propagates as a growing crystal. Simple and conceptually obvious in hindsight, this tells much about what we are empirically doing when we crystallize macromolecules. By adjusting the growth parameters we are empirically balancing the intermolecular interactions, preferentially attenuating the dominant strong (for lysozyme the charged groups) while strengthening the lesser strong (hydrophobic) interactions. In the general case for proteins the lack of a singularly defined association pathway may lead to formation of multiple species, i.e., amorphous precipitation. Weak interactions, such as hydrogen bonds, are promiscuous, serving to strengthen rather than define specific interactions. Participation in an interaction sequesters that surface from subsequent interactions, and we expect the strongest bonds to form first. This model, its basis, how it fits into the currently understood osmotic second virial coefficient approach to crystallization, and what it suggests will be discussed.

Pusey, Marc↗
NASA NTRS2011

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗