Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational principles”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Statistical Analysis of Factors Riving Surface Ozone Variability over Continental South Africa

Statistical relationships between surface ozone (O3) concentration, precursor species and meteorological conditions in continental South Africa were examined from data obtained from measurement stations in north-eastern South Africa. Three multivariate statistical methods were applied in the investigation, i.e. multiple linear regression (MLR), principal component analysis (PCA) and –regression (PCR), and generalised additive model (GAM) analysis. The daily maximum 8-h moving average O3 concentrations were considered in these statistical models (dependent variable). MLR models indicated that meteorology and precursor species concentrations are able to explain ~50% of the variability in daily maximum O3 levels. MLR analysis revealed that atmospheric carbon monoxide (CO), temperature and relative humidity were the strongest factors affecting the daily O3 variability. In summer, daily O3 variances were mostly associated with relative humidity, while winter O3 levels were mostly linked to temperature and CO. PCA indicated that CO, temperature and relative humidity were not strongly collinear. GAM also identified CO, temperature and relative humidity as the strongest factors affecting the daily variation of O3. Partial residual plots found that temperature, radiation and nitrogen oxides most likely have a non-linear relationship with O3,while the relationship with relative humidity and CO is probably linear. An inter-comparison between O3 levels modelled with the three statistical models compared to measured O3 concentrations showed that the GAM model offered a slight improvement over the MLR model. These findings emphasise the critical role of regional-scale O3 precursors coupled with meteorological conditions in daily variances of O3 levels in continental South Africa.

multiple linear regression (MLR)↗
NASA NTRSDate not supplied

Synergy of Satellite Radiation, Precipitation, and Other Meteorological Variable Observations for Global Mean Sea Surface Turbulent Heat Flux Estimation

Sea surface turbulent heat flux is one of the key components in global freshwater and energy balances and plays an important role in atmospheric dynamics, thermodynamics and general circulation. This turbulent heat flux is dominantly decided by sea surface latent heat release with some contribution from sensible heat exchange. The flux and its anomaly could significantly affect ocean heat storage and ocean circulation. They are part of and have great impacts on climate variability. Though the turbulent heat flux is extremely important for the climate and weather systems, only very limited ship and buoy observations of the flux are available over open oceans. There are no global, operational, direct measurements of this crucial meteorological variable. Global estimates are, basically, indirectly calculated from bulk turbulent flux parameterization with a combination of satellite column water vapor, sea surface water temperature, and wind speed observations. Critical parameters such as sea surface air temperature and humidity are estimated from empirical relations of theses variables with the water vapor and sea surface water temperature, respectively. Lacking accurate knowledge on surface air temperature and humidity, uncertainties in the parameterization and potential changes in the non-linear turbulent processes with long-term climate variations could cause large errors in estimated long-term turbulent fluxes from the indirect method as shown in current global sea surface turbulent flux datasets. This study uses synergized data of satellite global radiation, precipitation, and other meteorological variable observations to estimate sea surface turbulent heat fluxes. The radiation observations are made by the satellite Clouds and the Earth’s Radiant Energy System (CERES) sensors, while the global precipitation data is from NASA’s satellite Global Precipitation Climatology Project (GPCP). Other data includes satellite sea surface water temperature and wind speed observations. These datasets are obtained from a wide range of space sensors from passive to active instruments and from visible and near infrared to thermal infrared and microwave spectral sounders. One significant feature of these datasets is that they have multi-decades long climate records. Top-of-atmosphere (TOA) radiation and its anomaly represents the net heat energy input to the climate system, and the oceanic precipitation and ocean-land moisture transport can be used to quantify sea surface latent heat release. Based on the synergized datasets and the principle of global water and energy balances, global mean turbulent heat fluxes are estimated. The turbulent heat anomalies for the first two decades of the 21st century are, then, obtained mainly from global CERES radiation and GPCP precipitation anomalies, along with CERES derived ocean-land heat transports. Analysis indicates that the uncertainties in the estimated turbulent flux anomalies may be reduced considerably. The results suggest a strong needs in synchronized and synergized observations of atmospheric radiation, precipitation, and oceanic meteorological variables for long-term climate studies.

Bing Lin↗
NASA NTRSDate not supplied

Next Generation Exercise Device (NGED): Advancing Exercise Capabilities for Future Space Missions Through Biomechanical Modeling

BACKGROUND As space exploration extends to long-duration missions on the Moon and Mars, maintaining astronaut health and fitness becomes increasingly critical. The Next Generation Exercise Device (NGED), developed and tested by the HumanWorks Lab in NASA Johnson Space Center's (JSC) Software, Robotics, and Simulation Division, aims to address this challenge through innovative approaches. This study presents the development and evaluation of an NGED system, focusing on its adaptability to various mission scenarios, including prospective use in a Lunar Pressurized Rover (LPR). Central to this project is the application of biomechanical modeling to optimize exercise efficacy and safety in microgravity and partial gravity environments. The project is a collaborative effort with the Human Health and Performance group at Johnson Space Center, ensuring a comprehensive approach to astronaut well-being that integrates biomechanical principles with practical exercise solutions. The NGED represents the next generation of exercise capabilities for missions in space, on the Moon and Mars, with a specific focus on applications such as the LPR. METHODS AND RESULTS Data collection for NGED development was conducted with two motor-driven Beyond Power Voltra I [1] systems and a custom test structure to allow placement of the cable-based devices on the ground, at shoulder height, and overhead. The collection was performed in JSC’s Prototype Immersive Technology (PIT) Lab, utilizing an OptiTrack motion capture system and AMTI force platform, to enable detailed biomechanical analysis via OpenSim [2,3]. Motion capture data were collected for three subjects representing different body types and statures. The marker set used was an enhanced version of the full-body Plug-in Gait marker set [4], with additional markers strategically placed for the primary objective of informing exercise volume requirements. Subjects performed a series of 17 exercises, carefully selected to engage various muscle groups, including novel spaceflight exercises such as skiing (ergometer style), lateral pulldowns, wood chops, triceps extensions, and flies, with load variations ranging from 10 to 90 pounds to maintain kinematic form. This comprehensive approach allowed for a thorough evaluation of the NGED's performance across a wide range of motions and loads. The biomechanical modeling and analysis were conducted using a modified OpenSim Full Body Rajagopal Model [4,5] and also scaled to the maximum and minimum anthropometry provided in NASA-STD-3001 [6]. Volumetric convex hulls were generated based on model marker trajectories and aggregated into geometric assemblies. These can be placed in models of vehicle designs to assess fit to protect for exercise as well as to adapt NGED exercise to fit available space. Preliminary findings from the collection indicate that the NGED prototype demonstrates significant adaptability across varying user anthropometrics and exercise types. The device showed consistent performance in load-bearing exercises, with subjects able to perform exercises effectively while maintaining proper biomechanical form. CONCLUSION NGED represents a forward-looking advancement in exercise capabilities for future space missions. In the future, this system can be used to capture valuable metrics (e.g., isometric mid-thigh pull for force output measurements, assessments of postural muscle strength, overall isometric strength). Its versatility in accommodating various exercises and user physiques, coupled with the ability to provide targeted biomechanical loading, makes it a promising approach for maintaining astronaut health during long-duration missions to the Moon and Mars. Future work will focus on refining the NGED based on initial biomechanical findings, leveraging the detailed insights provided by motion capture and analysis techniques. Particular emphasis will be placed on optimizing its use within the confined spaces of a LPR and other space habitats. This work contributes significantly to NASA's goals of supporting human health and performance in deep space exploration, paving the way for sustainable long-term presence beyond Low Earth Orbit through advanced, biomechanically-informed exercise solutions.

C Wang↗
NASA NTRS2011

Receive Mode Analysis and Design of Microstrip Reflectarrays

Traditionally microstrip or printed reflectarrays are designed using the transmit mode technique. In this method, the size of each printed element is chosen so as to provide the required value of the reflection phase such that a collimated beam results along a given direction. The reflection phase of each printed element is approximated using an infinite array model. The infinite array model is an excellent engineering approximation for a large microstrip array since the size or orientation of elements exhibits a slow spatial variation. In this model, the reflection phase from a given printed element is approximated by that of an infinite array of elements of the same size and orientation when illuminated by a local plane wave. Thus the reflection phase is a function of the size (or orientation) of the element, the elevation and azimuth angles of incidence of a local plane wave, and polarization. Typically, one computes the reflection phase of the infinite array as a function of several parameters such as size/orientation, elevation and azimuth angles of incidence, and in some cases for vertical and horizontal polarization. The design requires the selection of the size/orientation of the printed element to realize the required phase by interpolating or curve fitting all the computed data. This is a substantially complicated problem, especially in applications requiring a computationally intensive commercial code to determine the reflection phase. In dual polarization applications requiring rectangular patches, one needs to determine the reflection phase as a function of five parameters (dimensions of the rectangular patch, elevation and azimuth angles of incidence, and polarization). This is an extremely complex problem. The new method employs the reciprocity principle and reaction concept, two well-known concepts in electromagnetics to derive the receive mode analysis and design techniques. In the "receive mode design" technique, the reflection phase is computed for a plane wave incident on the reflectarray from the direction of the beam peak. In antenna applications with a single collimated beam, this method is extremely simple since all printed elements see the same angles of incidence. Thus the number of parameters is reduced by two when compared to the transmit mode design. The reflection phase computation as a function of five parameters in the rectangular patch array discussed previously is reduced to a computational problem with three parameters in the receive mode. Furthermore, if the beam peak is in the broadside direction, the receive mode design is polarization independent and the reflection phase computation is a function of two parameters only. For a square patch array, it is a function of the size, one parameter only, thus making it extremely simple.

Rengarajan, Sembiam↗
NASA NTRS1994

A mechanism for magnetospheric substorms

Energy-principle analysis performed on two-dimensional, self-consistent solutions for magnetospheric convection indicates that the magnetosphere is unstable to isobaric (yet still frozen-in) fluctuations of plasma-sheet flux tubes. Normally, pdV work associated with compression maintains stability of the inward/outward oscillating normal mode. However, if Earth's ionosphere can provide sufficient mass flux, isobaric expansion of flux tubes can occur. The growth of a field-aligned potential drop in the near-Earth, midnight portion of the plasma sheet, associated with upward field-aligned currents responsible for the Harang discontinuity, redistributes plasma along field lines in a manner that destabilizes the normal mode. The growth of this unstable mode results in an out-of-equilibrium situation near the inner edge. When this occurs over a downtail extent comparable to the half-thickness of the plasma sheet, collapse ensues and forces thinning of the plasma sheet whereby conditions favorable to reconnection occur. This scenario for substorm onset is consistent with observed upward fluxes of ions, parallel potential drops, and observations of substorm onset. These observations include near Earth onset, pseudobreakups, the substorm current wedge, and local variations of plasma-sheet thickness.

Erickson, G. M.↗
NASA NTRS2000

A Parallel Genetic Algorithm for Automated Electronic Circuit Design

Parallelized versions of genetic algorithms (GAs) are popular primarily for three reasons: the GA is an inherently parallel algorithm, typical GA applications are very compute intensive, and powerful computing platforms, especially Beowulf-style computing clusters, are becoming more affordable and easier to implement. In addition, the low communication bandwidth required allows the use of inexpensive networking hardware such as standard office ethernet. In this paper we describe a parallel GA and its use in automated high-level circuit design. Genetic algorithms are a type of trial-and-error search technique that are guided by principles of Darwinian evolution. Just as the genetic material of two living organisms can intermix to produce offspring that are better adapted to their environment, GAs expose genetic material, frequently strings of 1s and Os, to the forces of artificial evolution: selection, mutation, recombination, etc. GAs start with a pool of randomly-generated candidate solutions which are then tested and scored with respect to their utility. Solutions are then bred by probabilistically selecting high quality parents and recombining their genetic representations to produce offspring solutions. Offspring are typically subjected to a small amount of random mutation. After a pool of offspring is produced, this process iterates until a satisfactory solution is found or an iteration limit is reached. Genetic algorithms have been applied to a wide variety of problems in many fields, including chemistry, biology, and many engineering disciplines. There are many styles of parallelism used in implementing parallel GAs. One such method is called the master-slave or processor farm approach. In this technique, slave nodes are used solely to compute fitness evaluations (the most time consuming part). The master processor collects fitness scores from the nodes and performs the genetic operators (selection, reproduction, variation, etc.). Because of dependency issues in the GA, it is possible to have idle processors. However, as long as the load at each processing node is similar, the processors are kept busy nearly all of the time. In applying GAs to circuit design, a suitable genetic representation 'is that of a circuit-construction program. We discuss one such circuit-construction programming language and show how evolution can generate useful analog circuit designs. This language has the desirable property that virtually all sets of combinations of primitives result in valid circuit graphs. Our system allows circuit size (number of devices), circuit topology, and device values to be evolved. Using a parallel genetic algorithm and circuit simulation software, we present experimental results as applied to three analog filter and two amplifier design tasks. For example, a figure shows an 85 dB amplifier design evolved by our system, and another figure shows the performance of that circuit (gain and frequency response). In all tasks, our system is able to generate circuits that achieve the target specifications.

Long, Jason D.↗
NASA NTRS2015

TTEthernet for Integrated Spacecraft Networks

Aerospace projects have traditionally employed federated avionics architectures, in which each computer system is designed to perform one specific function (e.g. navigation). There are obvious downsides to this approach, including excessive weight (from so much computing hardware), and inefficient processor utilization (since modern processors are capable of performing multiple tasks). There has therefore been a push for integrated modular avionics (IMA), in which common computing platforms can be leveraged for different purposes. This consolidation of multiple vehicle functions to shared computing platforms can significantly reduce spacecraft cost, weight, and design complexity. However, the application of IMA principles introduces significant challenges, as the data network must accommodate traffic of mixed criticality and performance levels - potentially all related to the same shared computer hardware. Because individual network technologies are rarely so competent, the development of truly integrated network architectures often proves unreasonable. Several different types of networks are utilized - each suited to support a specific vehicle function. Critical functions are typically driven by precise timing loops, requiring networks with strict guarantees regarding message latency (i.e. determinism) and fault-tolerance. Alternatively, non-critical systems generally employ data networks prioritizing flexibility and high performance over reliable operation. Switched Ethernet has seen widespread success filling this role in terrestrial applications. Its high speed, flexibility, and the availability of inexpensive commercial off-the-shelf (COTS) components make it desirable for inclusion in spacecraft platforms. Basic Ethernet configurations have been incorporated into several preexisting aerospace projects, including both the Space Shuttle and International Space Station (ISS). However, classical switched Ethernet cannot provide the high level of network determinism required by real-time spacecraft applications. Even with modern advancements, the uncoordinated (i.e. event-driven) nature of Ethernet communication unavoidably leads to message contention within network switches. The arbitration process used to resolve such conflicts introduces variation in the time it takes for messages to be forwarded. TTEthernet1 introduces decentralized clock synchronization to switched Ethernet, enabling message transmission according to a time-triggered (TT) paradigm. A network planning tool is used to allocate each device a finite amount of time in which it may transmit a frame. Each time slot is repeated sequentially to form a periodic communication schedule that is then loaded onto each TTEthernet device (e.g. switches and end systems). Each network participant references the synchronized time in order to dispatch messages at predetermined instances. This schedule guarantees that no contention exists between time-triggered Ethernet frames in the network switches, therefore eliminating the need for arbitration (and the timing variation it causes). Besides time-triggered messaging, TTEthernet networks may provide two additional traffic classes to support communication of different criticality levels. In the rate-constrained (RC) traffic class, the frame payload size and rate of transmission along each communication channel are limited to predetermined maximums. The network switches can therefore be configured to accommodate the known worst-case traffic pattern, and buffer overflows can be eliminated. The best-effort (BE) traffic class behaves akin to classical Ethernet. No guarantees are provided regarding transmission latency or successful message delivery. TTEthernet coordinates transmission of all three traffic classes over the same physical connections, therefore accommodating the full spectrum of traffic criticality levels required in IMA architectures. Common computing platforms (e.g. LRUs) can share networking resources in such a way that failures in non-critical systems (using BE or RC communication modes) cannot impact flight-critical functions (using TT communication). Furthermore, TTEthernet hardware (e.g. switches, cabling) can be shared by both TTEthernet and classical Ethernet traffic.

Loveless, Andrew↗
NASA NTRS1991

System and component design and test of a 10 hp, 18,000 rpm AC dynamometer utilizing a high frequency AC voltage link, part 1

Hard and soft switching test results conducted with one of the samples of first generation MOS-controlled thyristor (MCTs) and similar test results with several different samples of second generation MCT's are reported. A simple chopper circuit is used to investigate the basic switching characteristics of MCT under hard switching and various types of resonant circuits are used to determine soft switching characteristics of MCT under both zero voltage and zero current switching. Next, operation principles of a pulse density modulated converter (PDMC) for three phase (3F) to 3F two-step power conversion via parallel resonant high frequency (HF) AC link are reviewed. The details for the selection of power switches and other power components required for the construction of the power circuit for the second generation 3F to 3F converter system are discussed. The problems encountered in the first generation system are considered. Design and performance of the first generation 3F to 3F power converter system and field oriented induction moter drive based upon a 3 kVA, 20 kHz parallel resonant HF AC link are described. Low harmonic current at the input and output, unity power factor operation of input, and bidirectional flow capability of the system are shown via both computer and experimental results. The work completed on the construction and testing of the second generation converter and field oriented induction motor drive based upon specifications for a 10 hp squirrel cage dynamometer and a 20 kHz parallel resonant HF AC link is discussed. The induction machine is designed to deliver 10 hp or 7.46 kW when operated as an AC-dynamo with power fed back to the source through the converter. Results presented reveal that the proposed power level requires additional energy storage elements to overcome difficulties with a peak link voltage variation problem that limits reaching to the desired power level. The power level test of the second generation converter after the addition of extra energy storage elements to the HF link are described. The importance of the source voltage level to achieve a better current regulation for the source side PDMC is also briefly discussed. The power levels achieved in the motoring mode of operation show that the proposed power levels achieved in the generating mode of operation can also be easily achieved provided that no mechanical speed limitation were present to drive the induction machine at the proposed power level.

Lipo, Thomas A.↗
NASA NTRS2006

High-Performance Scanning Acousto-Ultrasonic System

A high-performance scanning acousto-ultrasonic system, now undergoing development, is designed to afford enhanced capabilities for imaging microstructural features, including flaws, inside plate specimens of materials. The system is expected to be especially helpful in analyzing defects that contribute to failures in polymer- and ceramic-matrix composite materials, which are difficult to characterize by conventional scanning ultrasonic techniques and other conventional nondestructive testing techniques. Selected aspects of the acousto-ultrasonic method have been described in several NASA Tech Briefs articles in recent years. Summarizing briefly: The acousto-ultrasonic method involves the use of an apparatus like the one depicted in the figure (or an apparatus of similar functionality). Pulses are excited at one location on a surface of a plate specimen by use of a broadband transmitting ultrasonic transducer. The stress waves associated with these pulses propagate along the specimen to a receiving transducer at a different location on the same surface. Along the way, the stress waves interact with the microstructure and flaws present between the transducers. The received signal is analyzed to evaluate the microstructure and flaws. The specific variant of the acousto-ultrasonic method implemented in the present developmental system goes beyond the basic principle described above to include the following major additional features: Computer-controlled motorized translation stages are used to automatically position the transducers at specified locations. Scanning is performed in the sense that the measurement, data-acquisition, and data-analysis processes are repeated at different specified transducer locations in an array that spans the specimen surface (or a specified portion of the surface). A pneumatic actuator with a load cell is used to apply a controlled contact force. In analyzing the measurement data for each pair of transducer locations in the scan, the total (multimode) acousto-ultrasonic response of the specimen is utilized. The analysis is performed by custom software that extracts parameters of signals in the time and frequency domains. The computer hardware and software provide both real-time and postscan processing and display options. For example, oscilloscope displays of waveforms and power spectral densities are available in real time. Images can be computed while scanning continues. Signals can be digitally preprocessed and/or post-processed by filtering, windowing, time-segmenting, and running-waveform-averaging algorithms. In addition, the software affords options for off-line simulation of the waveform-data-acquisition and scanning processes. In tests, the system has been shown to be capable of characterizing microstructural changes and defects in SiC/SiC and C/SiC ceramic-matrix composites. Delaminations, variations in density, microstructural changes attributable to infiltration by silicon, and crack-space indications (defined in the next sentence) have been revealed in images formed from several time- and frequency-domain parameters of scanning acousto-ultrasonic signals. The crack-space indications were image features that were not revealed by other nondestructive testing methods and are so named because they turned out to mark locations where cracking eventually occurred.

Roth, Don↗
NASA NTRSDate not supplied

Preliminary Results Cycling GEOS-JEDI with GSI-based Background Errors

The first phase of transitioning the NASA GMAO GEOS atmospheric data assimilation capabilities to JEDI involves the replacement of the Grid-point Statistical Interpolation (GSI) with a corresponding JEDI analysis. This includes taking JEDI's Unified Observation Operator (UFO), its underlying dependencies, and the JEDI solver that enables a hybrid 4DEnVar strategy similar to what is used in the current GEOS-GSI system. Variational analysis involves at least two main components associated with the observation and background cost function terms. The first is directly related to the UFO, which is being carefully validated in a joint collaboration between GMAO and NCEP to demonstrate consistency with corresponding observations usage in GSI. The second component is the background term, which in a hybrid system involves the ability to set up both a climatologically-based term and an ensemble-based term. JEDI provides the means to implement both terms through its BUMP component. Use of BUMP would require a complete re-tune of both climatological and ensemble, which is a non-trivial exercise we would prefer to avoid. As an alternative, the work here studies the results of interfacing the GSI-background error capability (GSIBEC) into JEDI through SABER. With this, the exact same background error covariance formulation used in GSI can be employed in JEDI without need for re-tuning. This brief summary covers the work done to interface GSIBEC into JEDI and shows preliminary results where the background error covariances of the control (GEOS-GSI) and experiment (GEOS-JEDI) are identical in corresponding cycling experiments. The cycling exercise is obviously preliminary and so much can be expected from GEOS-JEDI when compared to GEOS-GSI. There is still a number of features that need closer attention and although in some cases in principle ready to cycle have been intentionally either turned off or not fully exercised (e.g., VarBC is applied but not cycled). Other features are still pending implementation, one such example is the implementation of the Tangent Linear Normal Mode Constraint. Still, results are quite encouraging as hopefully the discussion here illustrates.

Ricardo Todling↗
NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗
NASA NTRS2013

Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot

The optical extinction of a laser due to scattering of particles is a well-known phenomenon. In a laboratory environment, this physical principle is known as the Beer-Lambert law, and is often used to measure the concentration of scattering particles in a fluid or gas. This method has been experimentally shown to be a usable means to measure the dust density from a rocket plume interaction with the lunar surface. Using the same principles and experimental arrangement, this technique can be applied to hydrometeor size distributions, and for launch-pad operations, specifically as a passive hail detection and measurement system. Calibration of a hail monitoring system is a difficult process. In the past, it has required comparison to another means of measuring hydrometeor size and density. Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through hail (or dust in the rocket case) yields an estimate of the second moment of the particle cloud, and hydrometeor size distribution in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rain and hail fall make indirect measurements of the drop-size distribution. Instruments that scatter microwaves off of hydrometeors, such as the WSR-88D (Weather Surveillance Radar 88 Doppler), vertical wind profilers, and microwave disdrometers, measure the sixth moment of the drop size distribution (DSD). By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain and hail yield a measurement of the DSD's second moment by way of the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required. Depending on the intensity of the hail fall rate for moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. For hail fall only, the laser path may be shorter because of greater scattering due to the properties of hailstones versus raindrops. A photodetector may replace the camera in automated installations. Laser-based rain and hail measurement systems are available, but they are based on measuring the interruption of a thin laser beam, thus counting individual hydrometeors. These systems are true disdrometers since they also measure size and velocity. The method reported here is a simple method, requiring far less processing, but it is not a disdrometer.

Lane, John↗
NASA NTRS1994

Proposed law of nature linking impacts, plume volcanism, and Milankovitch cycles to terrestrial vertebrate mass extinctions via greenhouse-embryo death coupling

A greenhouse-physiological coupling killing mechanism active among mammals, birds, and reptiles has been identified. Operating via environmental thermal effects upon the maternal core-skin blood flow critical to the survival and development of embryos, it reduces the flow of blood to the uterine tract. Today, during hot summers, this phenomena kills embryos on a vast, global scale. Because of sensitivity of many mammals to modern heat, a major modern greenhouse could reduce population numbers on a global scale, and potentially trigger population collapses in the more vulnerable parts of the world. In the geological past, the killing mechanism has likely been triggered into action by greenhouse warming via impact events, plume volcanism, and Earth orbital variations (Milankovitch cycles). Earth's biosphere is maintained and molded by the flow of energy from the solar energy source to Earth and on to the space energy sink (SES). This SES energy flow maintains Earth's biosphere and its living components, as open, intermediate, dissipative, nonequilibrium systems whose states are dependent upon the rate of energy flowing through them. Greenhouse gases such as CO2 in the atmosphere influence the SES energy flow rate. Steady-state flow is necessary for global ecological stability (autopoiesis). Natural fluctuations of the C cycle such as rapid releases of CO2 from the mantle, or oceans, disrupt steady-state SES flow. These fluctuations constantly challenge the biosphere; slowdown of SES energy flow drives it toward thermodynamical equilibrium and stagnation. Fluctuations induced by impact event, mantle plume volcanism, and Milankovitch cycles can grow into structure-breaking waves triggering major perturbations of Earth's C cycle and mass extinctions. A major C cycle perturbation involves readjustment of the outer physiochemical spheres of the Earth: the atmosphere, hydrosphere, and lithosphere; and by necessity, the biosphere. A greenhouse, one manifestation of a major C cycle perturbation, is the most dangerous natural phenomenon that life on Earth can experience. Greenhouse conditions existed during the KT mass extinctions of 65 m.y. ago, and the Pleistocene-Holocene (P-H) mammalian extinctions of 10,000-12,000 yr ago. Coupling climatology to reproductive physiology via effects of ambient air temperature upon uterine blood flow to developing embryos accounts for the extinctions via established physiological principles.

Mclean, D. M.↗
NASA NTRS2009

Spatial Modulation Improves Performance in CTIS

Suitably formulated spatial modulation of a scene imaged by a computed-tomography imaging spectrometer (CTIS) has been found to be useful as a means of improving the imaging performance of the CTIS. As used here, "spatial modulation" signifies the imposition of additional, artificial structure on a scene from within the CTIS optics. The basic principles of a CTIS were described in "Improvements in Computed- Tomography Imaging Spectrometry" (NPO-20561) NASA Tech Briefs, Vol. 24, No. 12 (December 2000), page 38 and "All-Reflective Computed-Tomography Imaging Spectrometers" (NPO-20836), NASA Tech Briefs, Vol. 26, No. 11 (November 2002), page 7a. To recapitulate: A CTIS offers capabilities for imaging a scene with spatial, spectral, and temporal resolution. The spectral disperser in a CTIS is a two-dimensional diffraction grating. It is positioned between two relay lenses (or on one of two relay mirrors) in a video imaging system. If the disperser were removed, the system would produce ordinary images of the scene in its field of view. In the presence of the grating, the image on the focal plane of the system contains both spectral and spatial information because the multiple diffraction orders of the grating give rise to multiple, spectrally dispersed images of the scene. By use of algorithms adapted from computed tomography, the image on the focal plane can be processed into an image cube a three-dimensional collection of data on the image intensity as a function of the two spatial dimensions (x and y) in the scene and of wavelength (lambda). Thus, both spectrally and spatially resolved information on the scene at a given instant of time can be obtained, without scanning, from a single snapshot; this is what makes the CTIS such a potentially powerful tool for spatially, spectrally, and temporally resolved imaging. A CTIS performs poorly in imaging some types of scenes in particular, scenes that contain little spatial or spectral variation. The computed spectra of such scenes tend to approximate correct values to within acceptably small errors near the edges of the field of view but to be poor approximations away from the edges. The additional structure imposed on a scene according to the present method enables the CTIS algorithms to reconstruct acceptable approximations of the spectral data throughout the scene.

Bearman, Gregory H.↗
NASA NTRS2012

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew↗
NASA NTRS2005

Fiber-Optic Sensor Would Monitor Growth of Polymer Film

A proposed optoelectronic sensor system would measure the increase in thickness of a film of parylene (a thermoplastic polymer made from para-xylene) during growth of the film in a vapor deposition process. By enabling real-time monitoring of film thickness, the system would make it possible to identify process conditions favorable for growth and to tailor the final thickness of the film with greater precision than is now possible. The heart of the sensor would be a pair of fiber-optic Fabry-Perot interferometers, depicted schematically in the figure. (In principle, a single such interferometer would suffice. The proposal calls for the use of two interferometers for protective redundancy and increased accuracy.) Each interferometer would include a light source, a fiber-optic coupler, and photodetectors in a control box outside the deposition chamber. A single-mode optical fiber for each interferometer would run from inside the control box to a fused-silica faceplate in a sensor head. The sensory tips of the optical fibers would be polished flush with the free surface of the faceplate. In preparation for use, the sensor head would be mounted with a hermetic seal in a feed-through port in the deposition chamber, such that free face of the faceplate and the sensory tips of the optical fibers would be exposed to the deposition environment. During operation, light would travel along each optical fiber from the control box to the sensor head. A small portion of the light would be reflected toward the control box from the end face of each fiber. Once growth of the parylene film started, a small portion of the light would also be reflected toward the control box from the outer surface of the film. In the control box, the two reflected portions of the light beam would interfere in one of the photodetectors. The difference between the phases of the interfering reflected portions of the light beam would vary in proportion to the increasing thickness of the film and the known index of refraction of the film, causing the photodetector reading to vary in proportion to a known sinusoidal function of film thickness. Electronic means of monitoring this variation and the corresponding variation in phase and thickness are well established in the art of interferometry. Hence, by tracking the cumulative change in phase difference from the beginning of deposition, one could track the growing thickness of the film to within a small fraction of a wavelength of light.

Beamesderfer, Michael↗
NASA NTRSDate not supplied

Organic-rich fluids in Pluto’s interior—A source of cryovolcanic activity

The highly variegated surface of Pluto revealed by the New Horizons spacecraft is composed of both volatile ices (N2, CH4, CO) and H2O ice, all of which are intrinsically colorless. The colors on Pluto reveal a non-ice component generally acknowledged to be a refractory complex organic material (tholins) produced by photolysis in the planet’s atmosphere [1], and by photolysis and radiolysis of the surface ices [2-4]. These processes have produced a range of colors from pale yellow to red to brown, and result in a variation of a factor of 10in albedo across the surface. Tholins are disordered polymer-like materials made of repeating chains of linked subunits and complex combinations of functional groups containing carbon. Nitrogen, oxygen, and other elements can be substituted in both the aliphatic and aromatic subunits. The detection of the spectral absorption of an ammoniated compound associated with surface exposures of red-colored H2O ice (herein, RAW) in the LEISA spectral mapping data from New Horizons suggests that a third source of a colored material is ejected from Pluto’s interior as a water-rich fluid(cryomagma) that rapidly freezes at the low temperature (~40K) of the surface [5,6]. The ammoniated compound is thought to be an ammonia hydrate or an ammoniated salt; the spectral data cover a limited wavelength range, resulting in ambiguity in the identification of the form of the ammonia. The distribution of RAW in the vicinity of tectonically stressed and deformed geological structures, particularly graben (fossae), is seen as the result of both fluid effusion and fountaining ejection of cryoclastic materials in Pluto’s Virgil Fossae region [5,6]. RAW is seen elsewhere on Pluto, both in association with graben complexes and on broad expanses of old surface, but primarily in zones of tectonic stress. The nature of the putative subsurface fluid is unknown, but it must have existed at relatively shallow depths such that the graben faulting could reach it. Geochemical modeling of fluids in the interiors of small planetary bodies in the outer Solar System[7] shows that ammonia-bearing H2O in contact with rock and metal components produce a rich chemical soup, including gases that can, in principle, help propel the fluid from shallow reservoirs through crustal fractures to reach the surface. An organic component of the fluid cryomagma is predicted to originate from thermal processing of the components [ref] and may be an example of one-pot synthesis of such complexes as amino acid precursors [ref].

Dale P Cruikshank↗
NASA NTRS2011

Frame-Transfer Gating Raman Spectroscopy for Time-Resolved Multiscalar Combustion Diagnostics

Accurate experimental measurement of spatially and temporally resolved variations in chemical composition (species concentrations) and temperature in turbulent flames is vital for characterizing the complex phenomena occurring in most practical combustion systems. These diagnostic measurements are called multiscalar because they are capable of acquiring multiple scalar quantities simultaneously. Multiscalar diagnostics also play a critical role in the area of computational code validation. In order to improve the design of combustion devices, computational codes for modeling turbulent combustion are often used to speed up and optimize the development process. The experimental validation of these codes is a critical step in accepting their predictions for engine performance in the absence of cost-prohibitive testing. One of the most critical aspects of setting up a time-resolved stimulated Raman scattering (SRS) diagnostic system is the temporal optical gating scheme. A short optical gate is necessary in order for weak SRS signals to be detected with a good signal- to-noise ratio (SNR) in the presence of strong background optical emissions. This time-synchronized optical gating is a classical problem even to other spectroscopic techniques such as laser-induced fluorescence (LIF) or laser-induced breakdown spectroscopy (LIBS). Traditionally, experimenters have had basically two options for gating: (1) an electronic means of gating using an image intensifier before the charge-coupled-device (CCD), or (2) a mechanical optical shutter (a rotary chopper/mechanical shutter combination). A new diagnostic technology has been developed at the NASA Glenn Research Center that utilizes a frame-transfer CCD sensor, in conjunction with a pulsed laser and multiplex optical fiber collection, to realize time-resolved Raman spectroscopy of turbulent flames that is free from optical background noise (interference). The technology permits not only shorter temporal optical gating (down to <1 s, in principle), but also higher optical throughput, thus resulting in a substantial increase in measurement SNR.

Nguyen, Quang-Viet↗