Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing technique”

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.

808 records · Page 45

Silicon Micromachined Microlens Array for THz Antennas

5 5 silicon microlens array was developed using a silicon micromachining technique for a silicon-based THz antenna array. The feature of the silicon micromachining technique enables one to microfabricate an unlimited number of microlens arrays at one time with good uniformity on a silicon wafer. This technique will resolve one of the key issues in building a THz camera, which is to integrate antennas in a detector array. The conventional approach of building single-pixel receivers and stacking them to form a multi-pixel receiver is not suited at THz because a single-pixel receiver already has difficulty fitting into mass, volume, and power budgets, especially in space applications. In this proposed technique, one has controllability on both diameter and curvature of a silicon microlens. First of all, the diameter of microlens depends on how thick photoresist one could coat and pattern. So far, the diameter of a 6- mm photoresist microlens with 400 m in height has been successfully microfabricated. Based on current researchers experiences, a diameter larger than 1-cm photoresist microlens array would be feasible. In order to control the curvature of the microlens, the following process variables could be used: 1. Amount of photoresist: It determines the curvature of the photoresist microlens. Since the photoresist lens is transferred onto the silicon substrate, it will directly control the curvature of the silicon microlens. 2. Etching selectivity between photoresist and silicon: The photoresist microlens is formed by thermal reflow. In order to transfer the exact photoresist curvature onto silicon, there needs to be etching selectivity of 1:1 between silicon and photoresist. However, by varying the etching selectivity, one could control the curvature of the silicon microlens. The figure shows the microfabricated silicon microlens 5 x5 array. The diameter of the microlens located in the center is about 2.5 mm. The measured 3-D profile of the microlens surface has a smooth curvature. The measured height of the silicon microlens is about 280 microns. In this case, the original height of the photoresist was 210 microns. The change was due to the etching selectivity of 1.33 between photoresist and silicon. The measured surface roughness of the silicon microlens shows the peak-to-peak surface roughness of less than 0.5 microns, which is adequate in THz frequency. For example, the surface roughness should be less than 7 microns at 600 GHz range. The SEM (scanning electron microscope) image of the microlens confirms the smooth surface. The beam pattern at 550 GHz shows good directivity.

Lee, Choonsup↗

Solar Cycle #24 and the Solar Dynamo

We focus on two solar aspects related to flight dynamics. These are the solar dynamo and long-term solar activity predictions. The nature of the solar dynamo is central to solar activity predictions, and these predictions are important for orbital planning of satellites in low earth orbit (LEO). The reason is that the solar ultraviolet (UV) and extreme ultraviolet (EUV) spectral irradiances inflate the upper atmospheric layers of the Earth, forming the thermosphere and exosphere through which these satellites orbit. Concerning the dynamo, we discuss some recent novel approaches towards its understanding. For solar predictions we concentrate on a solar precursor method, in which the Sun's polar field plays a major role in forecasting the next cycle s activity based upon the Babcock-Leighton dynamo. With a current low value for the Sun s polar field, this method predicts that solar cycle #24 will be one of the lowest in recent times, with smoothed F10.7 radio flux values peaking near 130 plus or minus 30 (2 sigma), in the 2013 timeframe. One may have to consider solar activity as far back as the early 20th century to find a cycle of comparable magnitude. Concomitant effects of low solar activity upon satellites in LEO will need to be considered, such as enhancements in orbital debris. Support for our prediction of a low solar cycle #24 is borne out by the lack of new cycle sunspots at least through the first half of 2007. Usually at the present epoch in the solar cycle (approx. 7+ years after the last solar maximum), for a normal size following cycle, new cycle sunspots would be seen. The lack of their appearance at this time is only consistent with a low cycle #24. Polar field observations of a weak magnitude are consistent with unusual structures seen in the Sun s corona. Polar coronal holes are the hallmarks of the Sun's open field structures. At present, it appears that the polar coronal holes are relatively weak, and there have been many equatorial coronal holes. This appears consistent with a weakening polar field, but coronal hole data must be scrutinized carefully as observing techniques have changed. We also discuss new solar dynamo ideas, and the SODA (SOlar Dynamo Amplitude) index, which provides the user with the ability to track the Sun's hidden, dynamo magnetic fields throughout the various stages of the Sun's cycle. Our solar dynamo ideas are a modernization and rejuvenation of the Babcock-Leighton original idea of a shallow solar dynamo, using modern observations that appear to support their shallow dynamo viewpoint. We are in awe of being able to see an object the size of the Sun undergoing as dramatic a change as our model provides in a few short years. The Sun, however, has undergone changes as rapid as this before! The weather on the Sun is at least as fickle as the weather on the Earth.

Schatten, Kenneth↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

The three-dimensional structure of cumulus clouds over the ocean. 1: Structural analysis

Thermal channel (channel 6, 10.4-12.5 micrometers) images of five Landsat thematic mapper cumulus scenes over the ocean are examined. These images are thresholded using the standard International Satellite Cloud Climatology Project (ISCCP) thermal threshold algorithm. The individual clouds in the cloud fields are segmented to obtain their structural statistics which include size distribution, orientation angle, horizontal aspect ratio, and perimeter-to-area (PtA) relationship. The cloud size distributions exhibit a double power law with the smaller clouds having a smaller absolute exponent. The cloud orientation angles, horizontal aspect ratios, and PtA exponents are found in good agreement with earlier studies. A technique also is developed to recognize individual cells within a cloud so that statistics of cloud cellular structure can be obtained. Cell structural statistics are computed for each cloud. Unicellular clouds are generally smaller (less than or equal to 1 km) and have smaller PtA exponents, while multicellular clouds are larger (greater than or equal to 1 km) and have larger PtA exponents. Cell structural statistics are similar to those of the smaller clouds. When each cell is approximated as a quadric surface using a linear least squares fit, most cells have the shape of a hyperboloid of one sheet, but about 15% of the cells are best modeled by a hyperboloid of two sheets. Less than 1% of the clouds are ellipsoidal. The number of cells in a cloud increases slightly faster than linearly with increasing cloud size. The mean nearest neighbor distance between cells in a cloud, however, appears to increase linearly with increasing cloud size and to reach a maximum when the cloud effective diameter is about 10 km; then it decreases with increasing cloud size. Sensitivity studies of threshold and lapse rate show that neither has a significant impact upon the results. A goodness-of-fit ratio is used to provide a quantitative measure of the individual cloud results. Significantly improved results are obtained after applying a smoothing operator, suggesting the eliminating subresolution scale variations with higher spatial resolution may yield even better shape analyses.

Kuo, Kwo-Sen↗

Carbon-Nanotube-Carpet Heat-Transfer Pads

Microscopic thermal-contact pads that include carpet-like arrays of carbon nanotubes have been invented for dissipating heat generated in integrated circuits and similarly sized single electronic components. The need for these or other innovative thermal-contact pads arises because the requisite high thermal conductances cannot be realized by scaling conventional macroscopic thermal-contact pads down to microscopic sizes. Overcoming limitations of conventional thermal-contact materials and components, the carbon-nanotube thermal-contact pads offer the high thermal conductivities needed to accommodate the high local thermal power densities of modern electronic circuits, without need for large clamping pressures, extreme smoothness of surfaces in contact, or gap-filling materials (e.g., thermally conductive greases) to ensure adequate thermal contact. Moreover, unlike some conventional thermal-contact components, these pads are reusable. The figure depicts a typical pad according to the invention, in contact with a rough surface on an electronic component that is to be cooled. Through reversible bending and buckling of carbon nanotubes at asperities on the rough surface, the pad yields sufficiently, under relatively low contact pressure, that thermal contact is distributed to many locations on the surface to be cooled, including valleys where contact would not ordinarily occur in conventional clamping of rigid surfaces. Hence, the effective thermal-contact area is greater than that achievable through scaling down of a macroscopic thermal-contact pad. The extremely high longitudinal thermal conductivities of the carbon nanotubes are utilized to conduct heat away from potential hot spots on the surface to be cooled. The fibers protrude from a layer of a filler material (Cu, Ag, Au, or metal-particle- filled gels), which provides both mechanical support to maintain the carbon nanotubes in alignment and thermal conductivity to enhance the diffusion of concentrated heat from the nanotubes into the larger adjacent volume of a heat sink. The array of carbon nanotubes, the filler material, and the heat sink are parts of a unitary composite structure that is fabricated as follows: 1. Using techniques that have been reported previously, the array of substantially perpendicularly oriented carbon nanotubes is grown on a metal, silicon, or other suitable thermally conductive substrate that is intended to become the heat sink. 2. By means of chemical vapor deposition, physical vapor deposition, plasma deposition, ion sputtering, electrochemical deposition, or casting from a liquid phase, some or all of the interstitial volume between carbon nanotubes is filled with the aforementioned layer of mechanically supporting, thermally conductive material. 3. To cause the carbon nanotubes to protrude the desired length from the filler material, an outer layer of filler is removed by mechanical polishing, chemical mechanical polishing, wet chemical etching, electrochemical etching, or dry plasma etching.

Li, Jun↗

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗

Calibration of the SNPP and NOAA 20 VIIRS Sensors for Continuity of the MODIS Climate Data Records

Accurate long-term sensor calibration and periodic re-processing to ensure consistency and continuity of atmospheric, land and ocean geophysical retrievals from space within the mission period and across different missions is a major requirement of climate data records. In this work, we applied the Multi-Angle Implementation of Atmospheric Correction (MAIAC)-based vicarious calibration technique over Libya-4 desert site to perform calibration analysis of Visible Infrared Imaging Radiometer Suite (VIIRS) on Suomi National Polar-orbiting Partnership (SNPP) and NOAA-20 satellites. For both VIIRS sensors we characterized residual linear calibration trends and cross-calibrated both sensors to MODerate resolution Imaging Spectroradiometer (MODIS) Aqua regarded as a calibration standard. The relative spectral response (RSR) differences were accounted for using the German Aerospace Center (DLR) Earth Sensing Imaging Spectrometer (DESIS) hyperspectral surface reflectance data. Our results agree with independent vicarious calibration results of both the MODIS/VIIRS Characterization Support Team as well as the CERES Imager and Geostationary Calibration Group within estimated uncertainty of 1–2%. Analysis of MAIAC geophysical products with the new calibration shows a high level of agreement of MAIAC aerosol, surface reflectance and NDVI records between MODIS and VIIRS. Excluding high aerosol optical depth (AOD), all three sensors agree in AOD with mean difference (MD) less than 0.01 and residual mean squared difference rmsd ∼ 0.04. Spectral geometrically normalized surface reflectance agrees within rmsd of 0.003–0.005 in the visible and 0.01–0.012 at longer wavelengths. The residual surface reflectance differences are fully explained by differences in spectral filter functions. Finally, difference in NDVI is characterized by rmsd ∼ 0.02 and MD less than 0.003 for NDVI based on VIIRS imagery bands I1/I2 and less than 0.01 for NDVI based on VIIRS radiometric bands M5/M7. In practical sense, these numbers indicate consistency and continuity in MAIAC records ensuring the smooth transition from MODIS to VIIRS.

MAIAC↗

Introduction of the M-85 high-speed rotorcraft concept

As a result of studying possible requirements for high-speed rotorcraft and studying many high-speed concepts, a new high-speed rotorcraft concept, designated as M-85, was derived. The M-85 is a helicopter that is reconfigured to a fixed-wing aircraft for high-speed cruise. The concept was derived as an approach to enable smooth, stable conversion between fixed-wing and rotary-wing while retaining hover and low-speed flight characteristics of a low disk loading helicopter. The name, M-85, reflects the high-speed goals of 0.85 Mach number at high altitude. For a high-speed rotorcraft, it is expected that a viable concept must be a cruise-efficient, fixed-wing aircraft so it may be attractive for a multiplicity of missions. It is also expected that a viable high-speed rotorcraft concept must be cruise efficient first and secondly, efficient in hover. What makes the M-85 unique is the large circular hub fairing that is large enough to support the aircraft during conversion between rotary-wind and fixed-wing modes. With the aircraft supported by this hub fairing, the rotor blades can be unloaded during the 100 percent change in rotor rpm. With the blades unloaded, the potential for vibratory loads would be lessened. In cruise, the large circular hub fairing would be part of the lifting system with additional lifting panels deployed for better cruise efficiency. In hover, the circular hub fairing would slightly reduce lift potential and/or decrease hover efficiency of the rotor system. The M-85 concept is described and estimated forward flight performance characteristics are presented in terms of thrust requirements and L/D with airspeed. The forward flight performance characteristics reflect recent completed wind tunnel tests of the wing concept. Also presented is a control system technique that is critical to achieving low oscillatory loads in rotary-wing mode. Hover characteristics, C(sub p) versus C(sub T) from test data, is discussed. Other techniques pertinent to the M-85 concept such as passively controlling inplane vibration during starting and stopping of the rotor system, aircraft control system, and rotor drive technologies are discussed.

Stroub, Robert H.↗

Adaptive Morphological Feature-Based Object Classifier for a Color Imaging System

Utilizing a Compact Color Microscope Imaging System (CCMIS), a unique algorithm has been developed that combines human intelligence along with machine vision techniques to produce an autonomous microscope tool for biomedical, industrial, and space applications. This technique is based on an adaptive, morphological, feature-based mapping function comprising 24 mutually inclusive feature metrics that are used to determine the metrics for complex cell/objects derived from color image analysis. Some of the features include: Area (total numbers of non-background pixels inside and including the perimeter), Bounding Box (smallest rectangle that bounds and object), centerX (x-coordinate of intensity-weighted, center-of-mass of an entire object or multi-object blob), centerY (y-coordinate of intensity-weighted, center-of-mass, of an entire object or multi-object blob), Circumference (a measure of circumference that takes into account whether neighboring pixels are diagonal, which is a longer distance than horizontally or vertically joined pixels), . Elongation (measure of particle elongation given as a number between 0 and 1. If equal to 1, the particle bounding box is square. As the elongation decreases from 1, the particle becomes more elongated), . Ext_vector (extremal vector), . Major Axis (the length of a major axis of a smallest ellipse encompassing an object), . Minor Axis (the length of a minor axis of a smallest ellipse encompassing an object), . Partial (indicates if the particle extends beyond the field of view), . Perimeter Points (points that make up a particle perimeter), . Roundness [(4(pi) x area)/perimeter(squared)) the result is a measure of object roundness, or compactness, given as a value between 0 and 1. The greater the ratio, the rounder the object.], . Thin in center (determines if an object becomes thin in the center, (figure-eight-shaped), . Theta (orientation of the major axis), . Smoothness and color metrics for each component (red, green, blue) the minimum, maximum, average, and standard deviation within the particle are tracked. These metrics can be used for autonomous analysis of color images from a microscope, video camera, or digital, still image. It can also automatically identify tumor morphology of stained images and has been used to detect stained cell phenomena (see figure).

McDowell, Mark↗

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗

Biomaterials Out of Thin Air: In Situ, On-Demand Printing of Advanced Biocomposites: A New Materials Design and Production Technique Using 3D-Printed Arrays of Bioengineered Cells

We have completed the proof of concept described in our Phase I proposal, a two-material array of nonstructural proteins. We created an implementation of each step in our technology concept and demonstrated its critical functionality. The biological chassis and printing hardware we created as part of this work can be re-used for future work by inserting a material coding region upstream of the fluorescent tag. Overall, we showed that our technology concept is sound. The mission benefit analyses, as described in our Phase I proposal, are complete and contained in this report. These calculations show that our technology can save hundreds of kilograms of upmass for a potential planetary human habit construction mission: the mass per habitat module can be reduced by approximately one third if the biomaterials are manufactured on Earth and included in the mission upmass, and the full 240 kg per module can be saved if the materials are derived entirely from in situ resources. Mass savings between these two extremes is expected for an actual mission, depending on the level of in situ resource extraction technology. We have shown that continued advancement of this technology concept for use in a space mission environment is justified. Our survey of future development pathways proved extremely informative in light of the lessons learned from our proof of concept work and mission scenario analyses. For example, we were able for the first time to distinguish between the levels of functionality provided by production of structural proteins, other polymers such as polysaccharides, and true organic-inorganic composites such as bone and mineralized shell. This new information represents a significant advance in formulating specific applications, and key enabling technologies, for our proposed concept. We surveyed potential collaborations with other projects and synergies with enabling technologies that are developing. We have received requests for collaboration from other institutions, including labs at Stanford University and Drexel University. We have also received visits from industry, including Organovo, a tissue engineering company, and Autodesk, a major 3D and materials design software company. Finally, we have been in touch with the team behind the 2013 NIAC Phase ll 'Super Ball Bot-Structures for Planetary Landing and Exploration' and are planning to develop our biomaterial printing technology with the goal of enabling tensegrity-based rovers such as theirs to use lighter, more robust materials. A smooth transition from TRL 2 to TRL 3 assumes that the implementations of the technology concept which demonstrate critical functionality are also pathways for future development; while this is the case for most hardware or software projects, the multidisciplinary nature of our project, particularly the biological aspect of it, means that this is not always true. For example, as part of this work we showed that although there are large number of known genetic parts that correspond to non-structural materials, this is not true for sequences for structural organic proteins, let alone biominerals. These realizations allowed us to further subdivide our concept into more detailed development areas, some of which are clearly established at TRL 3, others of which were newly identified sub-technologies moved from TRL 1 to TRL 2. Similarly, although a single feasibility /benefit analysis is sufficient for advancement from TRL 2 to TRL 3, not all potential benefits to a technology concept as broad in scope as ours are apparent at TRL 2. Both our future pathways survey and our proof of concept work highlighted that the true mass savings potential of our technology concept cannot be quantified without modification of existing materials modelling tools to take into account the possibility of positional materials properties customization. Therefore, we have simultaneously both advanced one potential set of applications of our technology concept from TRL 2 to TRL 3 and also identified a previously unknown set of applications and advanced it from TRL 1 to TRL 2. Overall, we have moved the original formulation of our concept forward from TRL 2 to TRL 3, and the expanded formulation of it presented in this document has been advanced from a combination of TRL 1 and early 1RL 2 to an overall late TRL 2. We have also identified the key areas necessary for both short-term and long-term advancement, and made recommendations for specific future work in the most promising directions. With future work on a 1-2 year timeframe to continue advancement to overall TRL 3, we will be well positioned to begin work on a specific space mission technology insertion path.

Biology↗

Locomotion training of legged robots using hybrid machine learning techniques

In this study artificial neural networks and fuzzy logic are used to control the jumping behavior of a three-link uniped robot. The biped locomotion control problem is an increment of the uniped locomotion control. Study of legged locomotion dynamics indicates that a hierarchical controller is required to control the behavior of a legged robot. A structured control strategy is suggested which includes navigator, motion planner, biped coordinator and uniped controllers. A three-link uniped robot simulation is developed to be used as the plant. Neurocontrollers were trained both online and offline. In the case of on-line training, a reinforcement learning technique was used to train the neurocontroller to make the robot jump to a specified height. After several hundred iterations of training, the plant output achieved an accuracy of 7.4%. However, when jump distance and body angular momentum were also included in the control objectives, training time became impractically long. In the case of off-line training, a three-layered backpropagation (BP) network was first used with three inputs, three outputs and 15 to 40 hidden nodes. Pre-generated data were presented to the network with a learning rate as low as 0.003 in order to reach convergence. The low learning rate required for convergence resulted in a very slow training process which took weeks to learn 460 examples. After training, performance of the neurocontroller was rather poor. Consequently, the BP network was replaced by a Cerebeller Model Articulation Controller (CMAC) network. Subsequent experiments described in this document show that the CMAC network is more suitable to the solution of uniped locomotion control problems in terms of both learning efficiency and performance. A new approach is introduced in this report, viz., a self-organizing multiagent cerebeller model for fuzzy-neural control of uniped locomotion is suggested to improve training efficiency. This is currently being evaluated for a possible patent by NASA, Johnson Space Center. An alternative modular approach is also developed which uses separate controllers for each stage of the running stride. A self-organizing fuzzy-neural controller controls the height, distance and angular momentum of the stride. A CMAC-based controller controls the movement of the leg from the time the foot leaves the ground to the time of landing. Because the leg joints are controlled at each time step during flight, movement is smooth and obstacles can be avoided. Initial results indicate that this approach can yield fast, accurate results.

Simon, William E.↗

Comparative measurement of visual stability in Earth and cosmic space (L-4)

The frequency and the intensity of space motion sickness was reviewed and investigated. Unusual induced-gravity situations, such as rotation, linear acceleration, parallel swinging etc. were investigated as well. Visually-induced motion sensation or distorted perception were evaluated with respect to visual stability by many investigators. These studies are all concerned with the effect of gravity and the induction of motion sickness through human visual perception. Direct investigation under the microgravity situation, such as parabolic flight, Skylab, and Spacelab, was carried out. These studies focused on how human visual stability is established through various sensory afferents in specific gravity conditions. Results from these investigations indicate that sensory mismatch probably plays an important role in space motion sickness. The interaction of visual, vestibular, and somatosensory perception is smoothly coordinated under normal gravitational conditions on the Earth in our daily life. When the cooperation is destroyed or a mismatch occurs among them, motion sickness may develop not only on the Earth but also in microgravity. The latter case may be the cause of space motion sickness or space adaption syndrome. How human beings obtain visual stability even with posture changes on the ground was investigated. Visual stability can be categorized as static or dynamic. Static visual stability is concerned with orientation and dynamic stability is concerned with object motion perception. The perception of visual stability is modified by many other sensations, such as somatosensory, vestibular, and muscle tension. We will mainly focus on modifications by vestibular inputs to visual perception produced by eye movements in microgravity. The Vestibular-Oscular Reflex (VOR) is a well-known characteristic which results from the relationship between eye mobility and vestibular afferent inputs. Eye movements also modify dynamic visual perception, such as perceived object motion velocity. The VOR is constantly simulated under 1-g conditions on Earth. In fact, human beings have been habituated and 'programmed' for orientation (visual stability) in their everyday, 1-g environment. When humans are exposed to a different gravity situation, this programmed behavior must change; that is, it is reprogrammed. This is called habituation or familiarization. We hope to examine how object motion perception is perturbed and subsequently adapted in the microgravity environment. This experiment is focused on the cooperation of visual, vestibular, and somatosensory perception coordination and how it is changed or reduced in space compared to 1-g environment. We will obtain information on the coordination between eye movement and neck muscle activity by using EOG and EMG. We will also collect data from Payload Specialists using a self-diagnostic questionnaire concerned with perceptual abnormality. When each sensory input function and its integration in the higher nervous system are well-characterized, then more effective techniques to control SAS may be developed.

Koga, Kazugo↗

Fabric geometry distortion during composites processing

Waviness and tow misalignment are often cited as possible causes of data scatter and lower compression stiffness and strength in textile composites. Strength differences of as much as 40 percent have been seen in composites that appear to have the same basic material and structural properties -- i.e., yarn orientation, yarn size, interlacing geometry. Fabric geometry distortion has been suggested as a possible reason for this discrepancy, but little quantitative data or substantial evidence exists. The focus of this research is to contribute to the present understanding of the causes and effects of geometric distortion in textile composites. The initial part of the study was an attempt to gather qualitative information on a variety of textile structures. Existing and new samples confirmed that structures with a significant direction presence would be more susceptible to distortion due to the compaction process. Thus, uniweaves (fiber vol frac: 54-72 percent) biaxial braids (vf: 34-58 percent) demonstrated very little fabric geometry distortion. In stitched panels, only slight buckling of z-direction stitches was observed, primarily near the surface. In contrast, for structures with high compaction ratios -- e.g., large cylindrical yarns (2.5:1) orpowder towpreg (4:1) -- there were visible distortions where previously smooth and periodic undulations were transformed to abrupt changes in direction. A controlled study of the effect of forming pressure on distortion was conducted on type 162 glass plain weave fabrics. Panels (6 x 6 in) were produced via a resin infusion type setup, but with an EPON 815 epoxy resin. Pressures ranging from hand layup to 200 psi were used (vf: 34-54 percent). Photomicrographs indicated that at pressures up to 50 psi, large changes in thickness were due primarily to resin squeeze out. At higher pressures, when intimate contact was made between the layers, there was some tow flattening and in-plane shifting to optimize nesting. However, even at 200 psi the period and amplitude of the tow undulation remained constant, suggesting that for this relatively fine fabric, distortions from compaction were not a problem. Because of the interest in using larger tows (to reduce cost) and more complex structures, tests were also run on 2D triaxial glass braid (113 yd/lb at 0, 225 yd/lb at +/- 45). Forming pressures of 20, 50, 200, and 500 psi were used, and short block compression tests were run. The 500 psi specimen had a 10 percent decrease in modulus and an almost 50 percent decrease in strength (vs. 20 psi). Because the total fiber wgt/panel was kept constant, the thickness varied from 0.32 to 0.22 in (49-70 percent vf). Yet, the strength value is clearly below what would be expected, even with the decrease in thickness. Photomicrographs of these samples will be taken to determine if more fabric distortion exists in the 500 psi specimens. Finally, because the ultimate goal is to be able to predict and control distortion in a variety of textile structures, a model compaction test was developed to directly measure the deformation of the tows during compaction. Layers of dry glass fabric were placed in a mold with a clear plexiglass window. The yarn amplitude and period was then calculated using image analysis of the videotaped deformation. Preliminary tests demonstrated the feasibility of this technique for simple fabrics with large tows.

Chen, Julie↗

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗