Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

LUVOIR Primary Mirror Segment Alignment Control With Joint Laser Metrology and Segment Edge Sensing

An approach is developed for the alignment and stability maintenance of the LUVOIR segmented primary mirror using a segment state estimation and wavefront control method based on a hybrid segment motion sensing architecture of laser truss metrology and segment edge sensors. Our current computer model was generated for LUVOIR Architecture Option A with a 15-meter aperture, 120segment primary mirror. The methodology and simulation results will be presented and analyzed. JPL has a long history of technology development in laser metrology and edge sensors, including work in SIM [7], Keck and TMT [8], CCAT [3] and LUVOIR [1]. We will discuss our current efforts of LUVOIR laser metrology and edge-sensor models development, showing sensitivities of sensor measurements to various mirror eigenmodes, removing global modes and strengthening weak modes by performing joint (hybrid) laser-metrology and edge sensing. We will define and derive an important performance metric called wavefront error multiplier (WEM), and show that WEM provides a simple link between sensor errors and the closed-loop (controlled) system wavefront error. We will show WEM values for several hybrid sensor configuration options studied. We will discuss an algorithm for mirror shape control and maintenance through segment state and wavefront estimations using joint edge-metrology sensing. We will compare simulated performance of mirror state estimation, wavefront estimation and wavefront control based on joint edge-metrology sensing among several sensor configurations, and show the impact of sensor error distributions on the segmented mirror alignment performance. Mirror shape control performance will be also be evaluated in terms of imaging contrast between inner working angles (IWA) and outer working angles (OWA) of a LUVOIR coronagraph. The research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. LUVOIR is a mission concept under study for the next astronomy Decadal Survey—the information presented in this paper is pre-decisional and is provided for planning and discussion purposes only.

Shelton, Chris↗
NASA NTRS2014

Performance Verification of the Gravity and Extreme Magnetism Small Explorer GEMS X-Ray Polarimeter

olarimetry is a powerful tool for astrophysical observations that has yet to be exploited in the X-ray band. For satellite-borne and sounding rocket experiments, we have developed a photoelectric gas polarimeter to measure X-ray polarization in the 2-10 keV range utilizing a time projection chamber (TPC) and advanced micro-pattern gas electron multiplier (GEM) techniques. We carried out performance verification of a flight equivalent unit (1/4 model) which was planned to be launched on the NASA Gravity and Extreme Magnetism Small Explorer (GEMS) satellite. The test was performed at Brookhaven National Laboratory, National Synchrotron Light Source (NSLS) facility in April 2013. The polarimeter was irradiated with linearly-polarized monochromatic X-rays between 2.3 and 10.0 keV and scanned with a collimated beam at 5 different detector positions. After a systematic investigation of the detector response, a modulation factor greater than or equal to 35% above 4 keV was obtained with the expected polarization angle. At energies below 4 keV where the photoelectron track becomes short, diffusion in the region between the GEM and readout strips leaves an asymmetric photoelectron image. A correction method retrieves an expected modulation angle, and the expected modulation factor, approximately 20% at 2.7 keV. Folding the measured values of modulation through an instrument model gives sensitivity, parameterized by minimum detectable polarization (MDP), nearly identical to that assumed at the preliminary design review (PDR).

Time Projection↗
NASA NTRS2014

AztechSat-1, a First Collaborative CubeSat Between NASA and Mexico

An international collaboration program for capacity building in space technology is being established between NASA and the Mexican Space Agency (AEM). Its objective is to develop a series of CubeSats collaboratively in a way that benefits both agencies. For AEM, this experience will support the development of human capital required for its space program. This project provides an opportunity for NASA future missions to benefit from commercially available technologies demonstrated in space. AEM was recently created by the Mexican Government to use and develop space technology in Mexico. Current AEM plans call for having the capability of developing its own communication satellites by 2030, so the need for having a Mexican cadre of space experts is mandatory. The program will start with the development and launching of AztechSat-1, the first AEM CubeSat project to be deployed from the International Space Station (ISS). Several public Mexican universities will participate in this project to demonstrate GlobalStar, a satellite communications provider using a constellation of LEO satellites, as an option for small satellite missions. AztechSat-1 mission objectives are:(1) Develop a flight ready CubeSat for deployment from the ISS; (2) Demonstrate GlobalStar´s Network as a method to advance communications in CubeSat missions. AztechSat-1 is envisioned as a 1U CubeSat and will be designed and constructed by Mexican graduate students from two major universities in Mexico: UNAM and IPN. The students will be mentored and supervised by researchers from NASA Ames Research Center (NARC). Flight acceptance tests are going to be done both at Mexican university labs as well as at NARC. Funding for the project will be provided concurrently by AEM, the Mexican participating universities, and other Mexican entities. The development of AztechSat-1 will emphasize the transfer of knowledge in space mission analysis and design, flight acceptance testing and other areas of space technology which are not well developed in Mexico yet. Mexican students will benefit from their mentoring from NASA researchers and will provide a multiplier effect in the development of the next generation of Mexican aerospace engineers. AztechSat-1 will be developed in a period of twelve months and will be launched by late 2015. AztechSat´s success will define the future of collaboration between NASA and AEM in capacity building. A plan for developing a series of AztechSat’s in a one per year basis, with increasing capabilities and complexity is currently being negotiated between NASA an AEM.

Martinez, Andres↗
NASA NTRS2004

Video Image Tracking Engine

A method and system for processing an image including capturing an image and storing the image as image pixel data. Each image pixel datum is stored in a respective memory location having a corresponding address. Threshold pixel data is selected from the image pixel data and linear spot segments are identified from the threshold pixel data selected.. Ihe positions of only a first pixel and a last pixel for each linear segment are saved. Movement of one or more objects are tracked by comparing the positions of fust and last pixels of a linear segment present in the captured image with respective first and last pixel positions in subsequent captured images. Alternatively, additional data for each linear data segment is saved such as sum of pixels and the weighted sum of pixels i.e., each threshold pixel value is multiplied by that pixel's x-location).

Howard, Richard T.↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗
NASA NTRS1998

A Synthesis of Fluid Dynamics and Quantum Chemistry for the Design of Nanoelectronics

In 1959, during a famous lecture entitled "There's Plenty of Room at the Bottom", Richard Feynman focused on the startling technical possibilities that would exist at the limit of miniaturization, that being atomically precise devices with dimensions in the nanometer range. A nanometer is both a convenient unit of length for medium to large sized molecules, and the root of the name of the new interdisciplinary field of "nanotechnology". Essentially, "nanoelectronics" denotes the goal of shrinking electronic devices, such as diodes and transistors, as well as integrated circuits of such devices that can perform logical operations, down to dimensions in the range of 100 nanometers. The thirty-year hiatus in the development of nanotechnology can figuratively be seen as a period of waiting for the bottom-up and atomically precise construction skills of synthetic chemistry to meet the top-down reductionist aspirations of device physics. The sub-nanometer domain of nineteenth-century classical chemistry has steadily grown, and state-of-the-art supramolecular chemistry can achieve atomic precision in non-repeating molecular assemblies of the size desired for nanotechnology. For nanoelectronics in particular, a basic understanding of the electron transport properties of molecules must also be developed. Quantum chemistry provides powerful computational methods that can accurately predict the properties of small to medium sized molecules on a desktop workstation, and those of large molecules if one has access to a supercomputer. Of the many properties of a molecule that quantum chemistry routinely predicts, the ability to carry a current is one that had not even been considered until recently. "Currently", there is a controversy over just how to define this key property. Reminiscent of the situation in high-Tc superconductivity, much of the difficulty arises from the different models that are used to simplify the complex electronic structure of real materials. A model-independent approach has been proposed, that sacrifices the plentiful molecular orbitals and Bloch functions of conventional approaches, for a single three-dimensional observable quantity, the electron momentum density Pi(sub rho). This quantity is simply the probability of any electron having momentum rho, multiplied by the total number of electrons in the system (the position of the electron is uncertain). We have explored the utility of this new approach in providing a fundamental understanding of the electron transport properties of molecules that have provi been nominated as candidates for components in the design of nanoelectronics; phenylene-ethynylene oligomers. Some of the molecular systems that have been studied are sketched below.

MacDougall, Preston J.↗
NASA NTRS2002

Submillimeter Laboratory Investigations: Spectroscopy and Collisions

Currently, millimeter-wave and submillimeter-wave spectroscopy is conducted in our laboratory on several different types of spectrometers. Our standard spectrometer utilizes the output of a phase-locked klystron operating in the 40-60 GHz region, which is sent into a crossed-waveguide harmonic generator, or "multiplier". The high frequency millimeter-and submillimeter-wave radiation is transmitted via quasi-optical techniques through an absorption cell and then onto a detector, which is either an InSb hot electron bolometer cooled to 1.4 K or a Si bolometer cooled to 0.3 K. The detector response is sent to a computer for measurement and analysis. The frequency range produced and detected in this manner goes from 80 GHz to upwards of 1 THz. Spectra are normally taken with source modulation, with line frequencies typically measured to an accuracy of 50-100 kHz. Higher accuracy is available when needed. Recently, we developed a new, broad-band spectrometer in our laboratory based on a free-running backward wave oscillator (BWO) of Russian manufacture as the primary source of radiation. The so-called FASSST (fast-scan submillimeter spectroscopic technique) system uses fast-scan and optical calibration methods rather than the traditional locking techniques. The output power from the BWO is split such that 90% goes into the absorption cell while 10% is coupled to a 40-meter Fabry-Perot cavity, which yields fringe? for frequency measurement. Results from this spectrometer on the spectrum of nitric acid (HNO3) show that 100 GHz of spectral data can be obtained in 5 seconds with a measurement accuracy of 50 kHz. Currently, the frequency range of the FASSST system in our laboratory is roughly 100-700 GHz.

Herbst, Eric↗
NASA NTRS1993

Techniques for estimating Space Station aerodynamic characteristics

A method was devised and calculations were performed to determine the effects of reflected molecules on the aerodynamic force and moment coefficients for a body in free molecule flow. A procedure was developed for determining the velocity and temperature distributions of molecules reflected from a surface of arbitrary momentum and energy accommodation. A system of equations, based on momentum and energy balances for the surface, incident, and reflected molecules, was solved by a numerical optimization technique. The minimization of a 'cost' function, developed from the set of equations, resulted in the determination of the defining properties of the flow reflected from the arbitrary surface. The properties used to define both the incident and reflected flows were: average temperature of the molecules in the flow, angle of the flow with respect to a vector normal to the surface, and the molecular speed ratio. The properties of the reflected flow were used to calculate the contribution of multiply reflected molecules to the force and moments on a test body in the flow. The test configuration consisted of two flat plates joined along one edge at a right angle to each other. When force and moment coefficients of this 90 deg concave wedge were compared to results that did not include multiple reflections, it was found that multiple reflections could nearly double lift and drag coefficients, with nearly a 50 percent increase in pitching moment for cases with specular or nearly specular accommodation. The cases of diffuse or nearly diffuse accommodation often had minor reductions in axial and normal forces when multiple reflections were included. There were several cases of intermediate accommodation where the addition of multiple reflection effects more than tripled the lift coefficient over the convex technique.

Thomas, Richard E.↗
NASA NTRS1992

Failure detection and fault management techniques for flush airdata sensing systems

Methods based on chi-squared analysis are presented for detecting system and individual-port failures in the high-angle-of-attack flush airdata sensing system on the NASA F-18 High Alpha Research Vehicle. The HI-FADS hardware is introduced, and the aerodynamic model describes measured pressure in terms of dynamic pressure, angle of attack, angle of sideslip, and static pressure. Chi-squared analysis is described in the presentation of the concept for failure detection and fault management which includes nominal, iteration, and fault-management modes. A matrix of pressure orifices arranged in concentric circles on the nose of the aircraft indicate the parameters which are applied to the regression algorithms. The sensing techniques are applied to the F-18 flight data, and two examples are given of the computed angle-of-attack time histories. The failure-detection and fault-management techniques permit the matrix to be multiply redundant, and the chi-squared analysis is shown to be useful in the detection of failures.

Whitmore, Stephen A.↗
NASA NTRS2016

Precision Oxygen Isotope Measurements of Two C-Rich Hydrated Interplanetary Dust Particles

Introduction: Chondritic-smooth IDPs (Interplanetary Dust Particles) are low porosity objects whose mineralogy is dominated by aqueous alteration products such as Mg-rich phyllosilicates (smectite and serpentine group) and Mg-Fe carbonate minerals. Their hydrated mineralogy combined with low atmospheric entry velocities have been used to infer an origin largely from asteroidal sources. Spectroscopic studies show that the types and abundance of organic matter in CS IDPs is similar to that in CP IDPs. Although CS IDPs show broad similarities to primitive carbonaceous chondrites, only a few particles have been directly linked to specific meteorite groups such as CM and CI chondrites based on the presence of diagnostic minerals. Many CS IDPs however, have carbon contents that greatly exceed that of known meteorite groups suggesting that they either may derive from comets or represent samples of more primitive parent bodies than do meteorites. It is now recognized that many large, dark primitive asteroids in the outer main belt, as well as some trans-Neptunian objects, show spectroscopic evidence for aqueous alteration products on their surfaces. Some CS IDPs exhibit large bulk D enrichments similar to those observed in the cometary CP IDPs. While hydrated minerals in comets have not been unambiguously identified to date, the presence of the smectite group mineral nontronite has been inferred from infrared spectra obtained from the ejecta from comet 9P/Tempel 1 during the Deep Impact mission. Recent observations of low temperature sulfide minerals in Stardust mission samples suggest that limited aqueous activity occurred on comet Wild-2. All of these observations, taken together, suggest that the high-carbon hydrated IDPs are abundant and important samples of primitive solar system objects not represented in meteorite collections. Oxygen isotopic compositions of chondrites reflect mixing between a 16O-rich reservoir and a 17O,18O-rich reservoir produced via mass-independent fractionation. The composition of the 16O-rich reservoir is well constrained but material representing the 17O,18O-rich end-member is rare. Self-shielding models predict that cometary water, presumed to represent this reservoir, should be enriched in 17O and 18O by greater than 200 per mille. The high-carbon hydrated IDPs may be among the best materials available to search for preserved "cometary" H2O signatures. In order to better understand the origin and evolution of these particles, we have obtained 10 hydrated interplanetary dust particles for coordinated mineralogical, isotopic and organic analyses. We have previously reported the results of mineralogical and O isotopic measurements of two hydrated IDPs; here we present results of O isotopic measurements of three additional IDPs. Samples and Methods: Three interplanetary dust particles (L2079C35, L2083D46 and L2083E46) were embedded in S and partially ultramicrotomed into approximately 70 nanometer sections for analysis via transmission electron microscopy (TEM). The remainders of the unsliced particles were removed from S and pressed into high purity Au foil that was cleaned with HF acid and annealed at 800 degrees Centigrade. The pressed IDPs were analyzed via electron microprobe analysis (EPMA) for quantitative bulk chemical analysis. After EPMA analysis, the IDPs were subjected to precision O isotope analysis with the UCLA Cameca IMS-1270 ion probe. A 20 kiloelectronvolt, 0.5 nanoangstrom Cs+ primary beam of approximately 15 micrometers diameter was used for each measurement. Small particles of San Carlos olivine and Burma spinel were pressed into the Au foil for use as standards to correct for instrumental mass fractionation. The detection system was configured for multicollection, with 16O measured on a Faraday cup, and 17O and 18O measured on electron multipliers (EMs). Individual analyses consisted of 15 cycles of 10 seconds per cycle. Additionally, two microtome thin sections were measured for H isotopic compositions with the JSC NanoSIMS 50L ion probe. An 8 picoangstrom, 16 kiloelectronvolt Cs plus primary beam was used. Measurements consisted of H, C, 12C, 16O, and 18O collected with EMs in multicollection. Terrestrial biotite and kerogen were used for isotopic standards. A significant challenge in O isotope measurement of hydrated minerals is the interference from 16OH at mass 17O. We ensured that the 17O and 16OH peaks were fully resolved by using a mass resolution of greater than 7000 and by careful analyses of San Carlos olivine, Burma spinel and chlorite hydrated mineral standards. The hydride was further suppressed with a cold finger attached to an LN2 dewar to trap volatiles in the sample chamber. All sputtered ions were counted (i.e. presputtering was not used); after applying background, yield and deadtime corrections, we performed a change-point analysis on our data via R in order to determine when the sample reached sputtering equilibrium; data points collected prior to the change point were excluded. Change-point analysis was also used to determine whether the IDP had completely sputtered. Results: Mineralogy. IDPs C35 and E46 exhibited hydrated mineralogies, Fe-Ni sulfide grains, nanoglobules and occasional enstatite grains distributed throughout a fine-grained Mg-Fe saponite matrix. C35 also contained breunnerite (Mg,Fe)CO3; solar flare tracks were observed in enstatite, indicating minimal atmospheric entry heating. The mineralogy of D46 is dominated by a large FeS grain with a minor component of adhering silicate material. D46 was strongly heated during atmorpheric entry as evidenced by a well-developed magnetite rim. EPMA analyses show that both C35 and E46 have high carbon contents of 20 weight percentage (approximately 6X CI). D46 contains approximately 6 wt.weight percentage C. A significant amount of the carbon is present as carbon nanoglobules. Results: Hydrogen isotopes: Although the bulk delta D values of both sections of L2079C35 were within error of SMOW (minus 33 plus or minus 19 per mille, 1 plus or minus 14 per mille 1 sigma), several delta D-rich hotspots were also identified, reaching 2000 per mille. As shown in Fig. 1, these hotspots are clearly associated with discrete carbonaceous inclusions that are akin to nanoglobules found in many meteorites and other IDPs. Oxygen Isotopes. Results of the oxygen isotope measurements are shown in Figure 1. The oxygen isotope composition for L2079C35 was delta 18O equals plus11.6 plus or minus 1.9per mille, delta 17O equals plus 7.9 plus or minus 1.9per mille (2 standard errors). The oxygen isotope composition for L2083D46 was delta 18O equals minus 8.1 plus or minus 1.9 per mille, delta 17O equals minus 6.4 plus or minus 3.1 per mille (2 standard errors). The oxygen isotope composition for L2083E46 was delta 18O equals plus 12.0 plus or minus 1.9 per mille, delta 17O equals plus 9.2 plus or minus 2.0 per mille (2 standard errors). Discussion: Despite mineralogical similarities to highly aqueously altered carbonaceous chondrites, the hydrated IDPs we analyzed have oxygen isotopic compositions that are distinct from matrix materials in the CI, CM, and CR chondrites. The IDPs plot along the Young-Russell line, with delta 17O values for C35 and E46 suggestive of interaction with a 16O-poor reservoir. However, we have thus far not observed evidence of extreme 16O-poor reservoirs expected from self-shielding models and observed in Acfer 094 simplectite. The high carbon contents of the IDPs also set them apart from known meteoritic samples. The lack of atmospheric entry heating effects are consistent with low encounter velocities and suggest either an asteroidal source, or a low inclination, low eccentricity cometary origin. Conclusions: The unusual oxygen isotopic compositions, high carbon contents, and the abundance of Drich nanoglobules, together, suggest that the high-carbon, hydrated IDPs are derived from a primitive source that is not yet represented in meteorite collections.

Snead, C. J.↗
NASA NTRS2003

System for Testing Thermal Insulation of Pipes

An apparatus and method have been developed for measuring the rates of leakage of heat into pipes carrying liquids, the purpose of the measurements being to quantify the thermal performance of the insulation system. The apparatus is designed primarily for testing pipes used to carry cryogenic liquids, but can also be used for measuring the thermal performance of other insulated pipes or piping systems. The basic measurement principle is straightforward: The outer surface of the pipe insulation is maintained at a fixed warmer temperature. The interior of the pipe is maintained in a narrow fixed lower-temperature range by means of a regular liquid (e.g., water) that is pumped through the pipe at a known flow rate or a cryogenic liquid (e.g., nitrogen) that is saturated at atmospheric pressure and replenished until steady-state conditions are achieved. In the case of water or another liquid pumped through, the inlet and outlet temperatures are measured and heat-leak power is calculated as the mass flow rate of the liquid multiplied by the specific heat of the liquid multiplied by the inlet-to-outlet temperature rise of the liquid. In the case of liquid nitrogen or another low-temperature boiling liquid, the heat-leak power is calculated as the rate of boil-off multiplied by the latent heat of vaporization of the liquid. Then the thermal-insulation performance of the pipe system can be calculated as a function of the measured heat-leak power, the inner and outer boundary temperatures, and the dimensions of the pipe. The apparatus can test as many as three pipes simultaneously. The pipes can have inner diameters up to .15 cm and outer diameters up to .20 cm. The lengths of the pipes may vary; typical lengths are of the order of 18 m. Two thermal guard boxes . one for each end of the pipe(s) under test . are used to make the inlet and outlet fluid connections to the pipe(s) (see figure). The connections include bellows that accommodate thermal expansion and contraction of the pipes. The guard boxes and pipe(s) are positioned so that the pipe(s) slope upward from the upstream to the downstream end at an angle of at least 2 . The upward slope allows vapor bubbles to accumulate at the downstream end. The thermal guard boxes keep the ends of the pipes at the lower interior temperature to prevent spurious lengthwise leakage of heat into the pipes. It is important to prevent this spurious heat leakage because, if it were allowed to occur, it could contribute a large error in the measured heat-leak power. The upstream thermal guard box includes a heat exchanger through which liquid flowing into the pipe(s) is subcooled to the saturation temperature corresponding to the ambient pressure. Conversely, this heat exchanger can also be used to warm the flowing liquid to a desired fixed temperature. The apparatus includes a temperature control device that is placed around each pipe under test. Each device is operated under thermostatic control to maintain the outer surface of the pipe insulation at the specified test temperature. All measurements are recorded on a portable data-acquisition system.

Fesmire, James E.↗
NASA NTRS1992

Archiving data from ground-based telescopes

The scientific throughput of a particular observing facility has been demonstrated to be multiplied with the operation of a data archive and its corresponding retrieval system. A requisite to achieve such an exploitation is a well structured observations catalog, i.e. a catalog that includes all information necessary to reduce and analyze the data even many years after its acquisition. At the same time, an information system is required that allow users to browse through the catalog at different levels of detail, adapting the amount of information presented to the actual needs of the user. Archiving data acquired with ground-based telescopes is particularly difficult because of the relative short life-time of instruments and detectors in comparison to the expected life-time of the archive. This feature differentiates ground-based originated archives radically from its spaceborne counterparts. The organization of the observations catalog becomes highly dependent on the capability of the archive to deal with new instrumental configurations. We introduce in this paper, the concept of a catalog database as opposed to the static catalog design currently in use in many archiving facilities, as a method to deal with this problem. We also present a brief review of activities currently in progress in this area.

Albrecht, M. A.↗
NASA NTRS2006

The X-ray Halo of GX5-1

Using Chandra observations we have measured the energy-resolved dust-scattered X-ray halo around the low-mass X-ray binary GX5-1, detecting for the first time multiply scattered X-rays from interstellar dust. % e compared the observed X-ray halo at various energies to predictions from a range of dust models. These fits used both smoothly-distributed dust as well as dust in clumped clouds, with CO and 21 cm observations helping to determine the position of the clouds along the line of sight. We found that the BARE-GR-B model of Zubko, Dwek & Arendt (2004) generally led to the best results, although inadequacies in both the overall model and the data limit our conclusions. We did find that the composite dust models of Zubko, Dwek & Arendt (2004), especially the "no carbon" models, gave uniformly poor results. Although models using cloud positions and densities derived naively from CO and 21 cm data gave generally poor results, plausible adjustments to the distance of the largest cloud and the mass of a cloud in the expanding 3 kpc Arm lead to significantly improved fits. We suggest that combining X-ray halo, CO, and 21 cm observations will be a fruitful method to improve our understanding of both the gas and dust phases of the interstellar medium.

Smith, Randall K.↗
NASA NTRS2012

Modeling Programs Increase Aircraft Design Safety

Flutter may sound like a benign word when associated with a flag in a breeze, a butterfly, or seaweed in an ocean current. When used in the context of aerodynamics, however, it describes a highly dangerous, potentially deadly condition. Consider the case of the Lockheed L-188 Electra Turboprop, an airliner that first took to the skies in 1957. Two years later, an Electra plummeted to the ground en route from Houston to Dallas. Within another year, a second Electra crashed. In both cases, all crew and passengers died. Lockheed engineers were at a loss as to why the planes wings were tearing off in midair. For an answer, the company turned to NASA s Transonic Dynamics Tunnel (TDT) at Langley Research Center. At the time, the newly renovated wind tunnel offered engineers the capability of testing aeroelastic qualities in aircraft flying at transonic speeds near or just below the speed of sound. (Aeroelasticity is the interaction between aerodynamic forces and the structural dynamics of an aircraft or other structure.) Through round-the-clock testing in the TDT, NASA and industry researchers discovered the cause: flutter. Flutter occurs when aerodynamic forces acting on a wing cause it to vibrate. As the aircraft moves faster, certain conditions can cause that vibration to multiply and feed off itself, building to greater amplitudes until the flutter causes severe damage or even the destruction of the aircraft. Flutter can impact other structures as well. Famous film footage of the Tacoma Narrows Bridge in Washington in 1940 shows the main span of the bridge collapsing after strong winds generated powerful flutter forces. In the Electra s case, faulty engine mounts allowed a type of flutter known as whirl flutter, generated by the spinning propellers, to transfer to the wings, causing them to vibrate violently enough to tear off. Thanks to the NASA testing, Lockheed was able to correct the Electra s design flaws that led to the flutter conditions and return the aircraft to safe flight. Today, all aircraft must have a flutter boundary 15 percent beyond the aircraft s expected maximum speed to ensure that flutter conditions are not encountered in flight. NASA continues to support research in new aircraft designs to improve knowledge of aeroelasticity and flutter. Through platforms such as Dryden Flight Research Center s Active Aeroelastic Wing (AAW) research aircraft, the Agency researches methods for in-flight validation of predictions and for controlling and taking advantage of aeroelastic conditions to enhance aircraft performance.

Source record↗
NASA NTRS1992

Method and apparatus for frequency spectrum analysis

A method for frequency spectrum analysis of an unknown signal in real-time is discussed. The method is based upon integration of 1-bit samples of signal voltage amplitude corresponding to sine or cosine phases of a controlled center frequency clock which is changed after each integration interval to sweep the frequency range of interest in steps. Integration of samples during each interval is carried out over a number of cycles of the center frequency clock spanning a number of cycles of an input signal to be analyzed. The invention may be used to detect the frequency of at least two signals simultaneously. By using a reference signal of known frequency and voltage amplitude (added to the two signals for parallel processing in the same way, but in a different channel with a sampling at the known frequency and phases of the reference signal), the absolute voltage amplitude of the other two signals may be determined by squaring the sine and cosine integrals of each channel and summing the squares to obtain relative power measurements in all three channels and, from the known voltage amplitude of the reference signal, obtaining an absolute voltage measurement for the other two signals by multiplying the known voltage of the reference signal with the ratio of the relative power of each of the other two signals to the relative power of the reference signal.

Cole, Steven W.↗
NASA NTRS2017

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗
NASA NTRS2004

Finite Element Modeling Used to Study Stress Distribution on the Foot

A method to study the stress distribution inside the forefoot during walking was developed at the Cleveland Clinic Foundation by a researcher from the NASA Glenn Research Center. In this method, a semiautomated process was outlined to create a three-dimensional, patient-specific, finite element model (FEM) of the forefoot using magnetic resonance images (MRI). The images were processed in Matlab using the k-nearest neighbor (k-NN) classification algorithm and Sobel edge detection to separate the different tissue types: bone, skin, fat, and muscle. This information was used to create curves and surfaces that were exported to an FEM preprocessor known as Truegrid. In Truegrid, eight-noded or brick elements were created by using surface mapping. The FEM was processed and postprocessed in Abaqus. Material properties of the models were obtained from past experiments such as fat pad confined compression, skin axial and biaxial tests, muscle in vivo compressive tests, and reference literature (bone properties). Nonlinear (hyperelastic) material models were used for the skin (epidermis and dermis), fat, and muscles; and a linear elastic model was used for the bones. Muscle activation during walking yielded uncertainties in the muscle material model since contracted muscles are stiffer than relaxed muscles. These uncertainties were resolved by performing a sensitivity analysis of the muscle material properties. The original properties were multiplied by arbitrary factors of 2, 3, 0.5, and 0.33. The strain and stress distributions, as well as the locations of peak values, were similar in all cases. The peak contact pressure P obtained for each case varied with respect to the applied factor f as follows:

Morales, Nelson↗
NASA NTRS2001

Portfolios in Stochastic Local Search: Efficiently Computing Most Probable Explanations in Bayesian Networks

Portfolio methods support the combination of different algorithms and heuristics, including stochastic local search (SLS) heuristics, and have been identified as a promising approach to solve computationally hard problems. While successful in experiments, theoretical foundations and analytical results for portfolio-based SLS heuristics are less developed. This article aims to improve the understanding of the role of portfolios of heuristics in SLS. We emphasize the problem of computing most probable explanations (MPEs) in Bayesian networks (BNs). Algorithmically, we discuss a portfolio-based SLS algorithm for MPE computation, Stochastic Greedy Search (SGS). SGS supports the integration of different initialization operators (or initialization heuristics) and different search operators (greedy and noisy heuristics), thereby enabling new analytical and experimental results. Analytically, we introduce a novel Markov chain model tailored to portfolio-based SLS algorithms including SGS, thereby enabling us to analytically form expected hitting time results that explain empirical run time results. For a specific BN, we show the benefit of using a homogenous initialization portfolio. To further illustrate the portfolio approach, we consider novel additive search heuristics for handling determinism in the form of zero entries in conditional probability tables in BNs. Our additive approach adds rather than multiplies probabilities when computing the utility of an explanation. We motivate the additive measure by studying the dramatic impact of zero entries in conditional probability tables on the number of zero-probability explanations, which again complicates the search process. We consider the relationship between MAXSAT and MPE, and show that additive utility (or gain) is a generalization, to the probabilistic setting, of MAXSAT utility (or gain) used in the celebrated GSAT and WalkSAT algorithms and their descendants. Utilizing our Markov chain framework, we show that expected hitting time is a rational function - i.e. a ratio of two polynomials - of the probability of applying an additive search operator. Experimentally, we report on synthetically generated BNs as well as BNs from applications, and compare SGSs performance to that of Hugin, which performs BN inference by compilation to and propagation in clique trees. On synthetic networks, SGS speeds up computation by approximately two orders of magnitude compared to Hugin. In application networks, our approach is highly competitive in Bayesian networks with a high degree of determinism. In addition to showing that stochastic local search can be competitive with clique tree clustering, our empirical results provide an improved understanding of the circumstances under which portfolio-based SLS outperforms clique tree clustering and vice versa.

Mengshoel, Ole J.↗