Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 3”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Conductance distribution near the Anderson transition

Using a modification of the Shapiro approach, we introduce the two-parameter family of conductance distributions W(g), defined by simple differential equations, which are in the one-to-one correspondence with conductance distributions for quasi-one-dimensional systems of size L{sup d–1} × L{sub z}, characterizing by parameters L/ξ and L{sub z}/L (ξ is the correlation length, d is the dimension of space). This family contains the Gaussian and log-normal distributions, typical for the metallic and localized phases. For a certain choice of parameters, we reproduce the results for the cumulants of conductance in the space dimension d = 2 + ϵ obtained in the framework of the σ-model approach. The universal property of distributions is existence of two asymptotic regimes, log-normal for small g and exponential for large g. In the metallic phase they refer to remote tails, in the critical region they determine practically all distribution, in the localized phase the former asymptotics forces out the latter. A singularity at g = 1, discovered in numerical experiments, is admissible in the framework of their calculational scheme, but related with a deficient definition of conductance. Apart of this singularity, the critical distribution for d = 3 is well described by the present theory. One-parameter scaling for the whole distribution takes place under condition, that two independent parameters characterizing this distribution are functions of the ratio L/ξ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Simultaneous mitigation of density and energy errors in approximate DFT for transition metal chemistry (Final Technical Report)

There were three major goals and objectives of this project: 1) Develop tools to understand density-driven errors in transition metal complexes, 2) Evaluate and minimize energy delocalization error and static correlation error through judicious functional choice, and 3) applying this workflow to machine-learning accelerated screening of redox couples. Over the reporting period, we developed a framework for eliminating flat plane errors. We introduced fully non-empirical coefficients. We demonstrated the approach on both molecules and solids. We investigated and eliminated density driven errors and demonstrated their impact on potential energy surfaces. We trained machine learning models both in a method-dependent fashion and to predict errors in method accuracy. We built large data sets of small molecule energetics and multi-reference character.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

A review of some problems in global-local stress analysis

The various types of local-global finite-element problems point out the need to develop a new generation of software. First, this new software needs to have a complete analysis capability, encompassing linear and nonlinear analysis of 1-, 2-, and 3-dimensional finite-element models, as well as mixed dimensional models. The software must be capable of treating static and dynamic (vibration and transient response) problems, including the stability effects of initial stress, and the software should be able to treat both elastic and elasto-plastic materials. The software should carry a set of optional diagnostics to assist the program user during model generation in order to help avoid obvious structural modeling errors. In addition, the program software should be well documented so the user has a complete technical reference for each type of element contained in the program library, including information on such topics as the type of numerical integration, use of underintegration, and inclusion of incompatible modes, etc. Some packaged information should also be available to assist the user in building mixed-dimensional models. An important advancement in finite-element software should be in the development of program modularity, so that the user can select from a menu various basic operations in matrix structural analysis.

Nelson, Richard B.↗
NASA NTRS1992

A brief overview of computational structures technology related activities at NASA Lewis Research Center

The presentation gives a partial overview of research and development underway in the Structures Division of LeRC, which collectively is referred to as the Computational Structures Technology Program. The activities in the program are diverse and encompass four major categories: (1) composite materials and structures; (2) probabilistic analysis and reliability; (3) design optimization and expert systems; and (4) computational methods and simulation. The approach of the program is comprehensive and entails exploration of fundamental theories of structural mechanics to accurately represent the complex physics governing engine structural performance, formulation, and implementation of computational techniques and integrated simulation strategies to provide accurate and efficient solutions of the governing theoretical models by exploiting the emerging advances in computer technology, and validation and verification through numerical and experimental tests to establish confidence and define the qualities and limitations of the resulting theoretical models and computational solutions. The program comprises both in-house and sponsored research activities. The remainder of the presentation provides a sample of activities to illustrate the breadth and depth of the program and to demonstrate the accomplishments and benefits that have resulted.

Hopkins, Dale A.↗
DOE OSTI2021

Prospects for Surrogate Neutron Capture Measurements with Radioactive Ion Beams and GODDESS

Neutron capture reactions are responsible for the synthesis of almost all of the elements heavier than iron through the slow s-process, that proceeds close to the line of stability, and the rapid r-process, with very neutron-rich waiting points. Uncertainties in (n,γ) rates in neutron rich nuclei, especially near closed neu- tron shells, can have significant impact [1] on the predictions of final abundances for different astrophysical scenarios for the r process. Understanding (n,γ) rates on neutron-rich fission fragments is also important for nuclear forensics and stockpile stewardship science. Ratkiewicz et al. [2 and references therein] has recently demonstrated that the (d,pγ) reaction is a valid surrogate for (n,γ), where the formation of the compound nucleus from the breakup of the deuteron has been calculated in a reaction model and the subsequent measured gamma-decay probabilities are reproduced with standard level density and strength functions in a Bayesian approach. In parallel to the surrogate validation efforts, we have demonstrated that the (d,pγ) reaction can be measured in inverse kinematics with Gammasphere ORRUBA: Dual Detectors for Experimental Structure Studies (GODDESS) [3] where the Gammasphere array of Compton-suppressed HPGe detectors is coupled to the Oak Ridge Rutgers University Barrel Array of position-sensitive silicon strip detectors. During the commissioning campaign we measured the (d,pγ) reaction with 134 Xe and 95 Mo beams, the latter to demonstrate the surrogate method in inverse kinematics. The present talk will present preliminary results from this campaign including γ-decay probabilities and prospects for surrogate (n,γ) measurements with 143 Ba fission-fragment beams.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗
DOE OSTI2020

Development of new baseline models for U.S. medium office buildings based on commercial buildings energy consumption survey data

Building energy estimation for the building sector under various scenarios are needed for building energy regulation and policy making. This often starts with representative baselines (either empirical baseline or modeled baseline). Commercial Buildings Energy Consumption Survey (CBECS) data is a widely used empirical baseline for U.S. commercial buildings, but none of the existing baseline model are developed to represent the CBECS data. This paper aims to develop new baseline models for the U.S. medium office buildings, which can produce modeled baselines consistent with the CBECS data. Here, we introduced the methodology to create baseline models and the criteria to evaluate the performance of baseline models. The methodology consists of three phases: (1) identification of model inputs, (2) model calibration, and (3) model validation with uncertainty analysis. The evaluation index is the coefficient of variation of the root-mean-square deviation (CV(RMSD)) of site energy use intensities (EUIs) between the modeled baseline and empirical baseline. Then 30 new baseline models for two vintages (pre- and post-1980) and 15 climate zones were created. The evaluation shows that the CV(RMSD) is lower than 0.05 for the modeled baselines produced by the new baseline models. As a comparison, the CV(RMSD) is higher than 0.1 for the existing modeled baselines generated by DOE Commercial Reference Building Models. Further analysis shows that the new baseline models are able to capture the uncertainties of the representative features of existing buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

SU-E-T-05: A 2D EPID Transit Dosimetry Model Based On An Empirical Quadratic Formalism

Purpose: To describe a 2D electronic portal imaging device (EPID) transit dosimetry model, based on an empirical quadratic formalism, that can predict either EPID or in-phantom dose distribution for comparisons with EPID captured image or treatment planning system (TPS) dose respectively. Methods: A quadratic equation can be used to relate the reduction in intensity of an exit beam to the equivalent path length of the attenuator. The calibration involved deriving coefficients from a set of dose planes measured for homogeneous phantoms with known thicknesses under reference conditions. In this study, calibration dose planes were measured with EPID and ionisation chamber (IC) in water for the same reference beam (6MV, 100mu, 20×20cm{sup 2}) and set of thicknesses (0–30cm). Since the same calibration conditions were used, the EPID and IC measurements can be related through the quadratic equation. Consequently, EPID transit dose can be predicted from TPS exported dose planes and in-phantom dose can be predicted using EPID distribution captured during treatment as an input. The model was tested with 4 open fields, 6 wedge fields, and 7 IMRT fields on homogeneous and heterogeneous phantoms. Comparisons were done using 2D absolute gamma (3%/3mm) and results were validated against measurements with a commercial 2D array device. Results: The gamma pass rates for comparisons between EPID measured and predicted ranged from 93.6% to 100.0% for all fields and phantoms tested. Results from this study agreed with 2D array measurements to within 3.1%. Meanwhile, comparisons in-phantom between TPS computed and predicted ranged from 91.6% to 100.0%. Validation with 2D array device was not possible for inphantom comparisons. Conclusion: A 2D EPID transit dosimetry model for treatment verification was described and proven to be accurate. The model has the advantage of being generic and allows comparisons at the EPID plane as well as multiple planes in-phantom.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1985

Circulation of the atmosphere from the surface to 100 km

A reference model of the atmospheric circulation on Venus based on available observations is presented. The reference atmosphere has the following main features: (1) the entire atmosphere below 85 km moves predominantly from east to west in the planet's reference system (in the same direction as the rotation of the solid planet itself) with the possible exception of the lowest 10 km where velocities are low, (2) a jet is present near the cloud-top level at 45 deg latitude in both hemispheres, with a magnitude of approximately 100 m/s, (3) a weak meridional (north-south component) flow directed towards either pole is superimposed on the zonal (east-west motion at cloud-top level (about 68 km) altitudes, and, (4) eddies or wave motions are present in the atmosphere, with amplitudes of less than 16/ms in the upper atmosphere.

Kerzhanovich, V. V.↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Zero Voltage Switching AC-DC Converter Based on Zero State Modulation (ZSM)

In this paper, a half-bridge building block that can achieve soft switching across the complete operating range while operating under fixed frequency is proposed. In addition, dead-beat model predictive controls to control the same are also proposed. The proposed converter, referred to as zero-state modulated (ZSM) converter, uses a zero state, enabled by a small auxiliary switch, to achieve soft-switching. The approach allows lower switching loss, controlled dv/dt, small filter inductor, fixed frequency operation allowing interleaving to reduce filter capacitor. A single half-bridge building block can used to realize a ZVS DC-DC converter or multiple of these building blocks can be used to achieve soft switching 1-phase or 3-phase AC-DC converters.

Kandula, Prasad↗
NASA NTRS1992

Collected Software Engineering Papers, Volume 10

This document is a collection of selected technical papers produced by participants in the Software Engineering Laboratory (SEL) from Oct. 1991 - Nov. 1992. The purpose of the document is to make available, in one reference, some results of SEL research that originally appeared in a number of different forums. Although these papers cover several topics related to software engineering, they do not encompass the entire scope of SEL activities and interests. Additional information about the SEL and its research efforts may be obtained from the sources listed in the bibliography at the end of this document. For the convenience of this presentation, the 11 papers contained here are grouped into 5 major sections: (1) the Software Engineering Laboratory; (2) software tools studies; (3) software models studies; (4) software measurement studies; and (5) Ada technology studies.

Source record↗
NASA NTRS1981

Numerical Prediction Methods (Reynolds-Averaged Navier-Stokes Simulations of Transonic Separated Flows)

During the past five years, numerous pioneering archival publications have appeared that have presented computer solutions of the mass-weighted, time-averaged Navier-Stokes equations for transonic problems pertinent to the aircraft industry. These solutions have been pathfinders of developments that could evolve into a major new technological capability, namely the computational Navier-Stokes technology, for the aircraft industry. So far these simulations have demonstrated that computational techniques, and computer capabilities have advanced to the point where it is possible to solve forms of the Navier-Stokes equations for transonic research problems. At present there are two major shortcomings of the technology: limited computer speed and memory, and difficulties in turbulence modelling and in computation of complex three-dimensional geometries. These limitations and difficulties are the pacing items of the continuing developments, although the one item that will most likely turn out to be the most crucial to the progress of this technology is turbulence modelling. The objective of this presentation is to discuss the state of the art of this technology and suggest possible future areas of research. We now discuss some of the flow conditions for which the Navier-Stokes equations appear to be required. On an airfoil there are four different types of interaction of a shock wave with a boundary layer: (1) shock-boundary-layer interaction with no separation, (2) shock-induced turbulent separation with immediate reattachment (we refer to this as a shock-induced separation bubble), (3) shock-induced turbulent separation without reattachment, and (4) shock-induced separation bubble with trailing edge separation.

Mehta, Unmeel↗
NASA NTRS2013

The Effect of Bypass Nozzle Exit Area on Fan Aerodynamic Performance and Noise in a Model Turbofan Simulator

The design of effective new technologies to reduce aircraft propulsion noise is dependent on identifying and understanding the noise sources and noise generation mechanisms in the modern turbofan engine, as well as determining their contribution to the overall aircraft noise signature. Therefore, a comprehensive aeroacoustic wind tunnel test program was conducted called the Fan Broadband Source Diagnostic Test as part of the NASA Quiet Aircraft Technology program. The test was performed in the anechoic NASA Glenn 9- by 15-Foot Low Speed Wind Tunnel using a 1/5 scale model turbofan simulator which represented a current generation, medium pressure ratio, high bypass turbofan aircraft engine. The investigation focused on simulating in model scale only the bypass section of the turbofan engine. The test objectives were to: identify the noise sources within the model and determine their noise level; investigate several component design technologies by determining their impact on the aerodynamic and acoustic performance of the fan stage; and conduct detailed flow diagnostics within the fan flow field to characterize the physics of the noise generation mechanisms in a turbofan model. This report discusses results obtained for one aspect of the Source Diagnostic Test that investigated the effect of the bypass or fan nozzle exit area on the bypass stage aerodynamic performance, specifically the fan and outlet guide vanes or stators, as well as the farfield acoustic noise level. The aerodynamic performance, farfield acoustics, and Laser Doppler Velocimeter flow diagnostic results are presented for the fan and four different fixed-area bypass nozzle configurations. The nozzles simulated fixed engine operating lines and encompassed the fan stage operating envelope from near stall to cruise. One nozzle was selected as a baseline reference, representing the nozzle area which would achieve the design point operating conditions and fan stage performance. The total area change from the smallest to the largest nozzle was 12.9 percent of the baseline nozzle area. The results will show that there are significant changes in aerodynamic performance and farfield acoustics as the fan nozzle area is increased. The weight flow through the fan model increased between 7 and 9 percent, the fan and stage pressure dropped between 8 and 10 percent, and the adiabatic efficiency increased between 2 and 3 percent--the magnitude of the change dependent on the fan speed. Results from force balance measurements of fan and outlet guide vane thrust will show that as the nozzle exit area is increased the combined thrust of the fan and outlet guide vanes together also increases, between 2 and 3.5 percent, mainly due to the increase in lift from the outlet guide vanes. In terms of farfield acoustics, the overall sound power level produced by the fan stage dropped nearly linearly between 1 dB at takeoff condition and 3.5 dB at approach condition, mainly due to a decrease in the broadband noise levels. Finally, fan swirl angle survey and Laser Doppler Velocimeter mean velocity and turbulence data obtained in the fan wake will show that the swirl angles and turbulence levels within the wake decrease as the fan nozzle area increases, which helps to explain the drop in the fan broadband noise at all fan speeds.

Hughes, Christopher E.↗
DOE OSTI2014

Retrospective Analysis of NIST Standard Reference Material 1450, Fibrous Glass Board, for Thermal Insulation Measurements

Thermal conductivity data acquired previously for the establishment of Standard Reference Material (SRM) 1450, Fibrous Glass Board, as well as subsequent renewals 1450a, 1450b, 1450c, and 1450d, are re-analyzed collectively and as individual data sets. Additional data sets for proto-1450 material lots are also included in the analysis. The data cover 36 years of activity by the National Institute of Standards and Technology (NIST) in developing and providing thermal insulation SRMs, specifically high-density molded fibrous-glass board, to the public. Collectively, the data sets cover two nominal thicknesses of 13 mm and 25 mm, bulk densities from 60 kg∙m -3 to 180 kg∙m -3 , and mean temperatures from 100 K to 340 K. The analysis repetitively fits six models to the individual data sets. The most general form of the nested set of multilinear models used is given in the following equation: λ(ρ, T) = α 0 + α 1 ρ +α 2 T + α 3 T 3 +α 4 e - ($\frac{T-a_{5}}{a_{6}}$) 2 = where λ(ρ,T) is the predicted thermal conductivity (W∙m -1 ∙K -1 ), ρ is the bulk density (kg∙m -3 ), T is the mean temperature (K) and $a_i$ (for i = 1, 2, … 6) are the regression coefficients. The least squares fit results for each model across all data sets are analyzed using both graphical and analytic techniques. The prevailing generic model for the majority of data sets is the bilinear model in ρ and T. λ(ρ, T) = α 0 + α 1 ρ +α 2 T One data set supports the inclusion of a cubic temperature term and two data sets with low-temperature data support the inclusion of an exponential term in T to improve the model predictions. Physical interpretations of the model function terms are described. Recommendations for future renewals of SRM 1450 are provided. An Addendum provides historical background on the origin of this SRM and the influence of the SRM on external measurement programs.

regression analysis↗
DOE OSTI2016

Nitride Fuel Potential for Power Uprates

This work presents the methodology and current results for design analysis of grid supported pressurized water reactor cores using uranium nitride (UN) fuel in square arrays. Uranium nitride has a higher theoretical density and higher thermal conductivity than uranium oxide, allowing for lower fuel temperatures, longer cycle length and higher power densities with similar enrichment. In order to avoid the UN and high temperature water reaction, silicon carbide cladding has been proposed in place of Zircaloy cladding. The analysis aims at finding the maximum thermal power over a range of bundle geometries (rod diameter and pitch). The same methodology is applied to the design of UO{sub 2} oxide fueled cores to provide a fair comparison of the achievable power rating between the two fuel types. Steady-state limits are considered, including fuel bundle pressure drop, minimum departure from nucleate boiling ratio (MDNBR), fuel temperature, fuel rod vibrations and pin fretting and sliding wear. Transients such as loss of flow and overpower accident will be considered. A steady state analysis was undertaken to determine the maximum achievable power that can be sustained by square array fuel assemblies loaded with nitride and oxide fuels with same or better safety performances as the reference design. Steady-state design limits were applied including MDNBR, fuel bundle pressure drop, fuel temperature, cumulative sliding and fretting wear. The major conclusions of this study are as follows: 1) It is possible to achieve significant power increase for UN and UO{sub 2} compared to the reference core power output of 3400 MWth with smaller but more numerous rods and higher pitch over diameter ratios. A power uprate of 25% can be reached according to the VIPRE-MATLAB interface for full core analysis with the same number of spacer grids, and 32% according to the new MATLAB script for single hot channel analysis, with five more spacer grids. For this latter core, a higher enrichment would be required. 2) The nitride fuel temperature constraint from dissociation doesn't influence steady state performances as other design constraints such as MDNBR and pressure loss are more restrictive. The oxide fuel temperature constraint from fission gas release is reached in part of the study domain, but not in the region that allows for possible uprates. 3) As a result, the thermal hydraulic performance of UN and UO{sub 2} are similar. Differences in possible power uprates will arise from fretting and sliding wear. Wear depends of the time spent by the fuel in the core, which depends on burnup and heavy metal inventory. UN being denser, it can allow for power uprates or increased cycle length. There is a trade off between power uprates and cycle length with regards to fuel pin wear. Reducing fuel pin wear by reducing cycle length can allow for higher thermal power. This affects the plants' economics as well. This work will be completed by a transient analysis using CASMO4E, SIMULATE3 and SIMULATE3k to model the neutronics of a nitride fueled core. The transient performances will be used to determine the maximum power uprates while performing as well or better during transients as the reference design. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2000

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗
DOE OSTI2015

TEM Video Compressive Sensing

One of the main limitations of imaging at high spatial and temporal resolution during in-situ TEM experiments is the frame rate of the camera being used to image the dynamic process. While the recent development of direct detectors has provided the hardware to achieve frame rates approaching 0.1ms, the cameras are expensive and must replace existing detectors. In this paper, we examine the use of coded aperture compressive sensing methods [1, 2, 3, 4] to increase the framerate of any camera with simple, low-cost hardware modifications. The coded aperture approach allows multiple sub-frames to be coded and integrated into a single camera frame during the acquisition process, and then extracted upon readout using statistical compressive sensing inversion. Our simulations show that it should be possible to increase the speed of any camera by at least an order of magnitude. Compressive Sensing (CS) combines sensing and compression in one operation, and thus provides an approach that could further improve the temporal resolution while correspondingly reducing the electron dose rate. Because the signal is measured in a compressive manner, fewer total measurements are required. When applied to TEM video capture, compressive imaging couled improve acquisition speed and reduce the electron dose rate. CS is a recent concept, and has come to the forefront due the seminal work of Candès [5]. Since the publication of Candès, there has been enormous growth in the application of CS and development of CS variants. For electron microscopy applications, the concept of CS has also been recently applied to electron tomography [6], and reduction of electron dose in scanning transmission electron microscopy (STEM) imaging [7]. To demonstrate the applicability of coded aperture CS video reconstruction for atomic level imaging, we simulate compressive sensing on observations of Pd nanoparticles and Ag nanoparticles during exposure to high temperatures and other environmental conditions. Figure 1 highlights the results from the Pd nanoparticle experiment. On the left, 10 frames are reconstructed from a single coded frame—the original frames are shown for comparison. On the right a selection of three frames are shown from reconstructions at compression levels 10,20,30. The reconstructions, which are not post-processed, are true to the original and degrade in a straightforward manner. The final choice of compression level will obviously depend on both the temporal and spatial resolution required for a specific imaging task, but the results indicate that an increase in speed of better than an order of magnitude should be possible for all experiments. References: [1] P Llull, X Liao, X Yuan et al. Optics express 21(9), (2013), p. 10526. [2] J Yang, X Yuan, X Liao et al. Image Processing, IEEE Trans 23(11), (2014), p. 4863. [3] X Yuan, J Yang, P Llull et al. In ICIP 2013 (IEEE), p. 14. [4] X Yuan, P Llull, X Liao et al. In CVPR 2014. p. 3318. [5] EJ Candès, J Romberg and T Tao. Information Theory, IEEE Trans 52(2), (2006), p. 489. [6] P Binev, W Dahmen, R DeVore et al. In Modeling Nanoscale Imaging in Electron Microscopy, eds. T Vogt, W Dahmen and P Binev (Springer US), Nanostructure Science and Technology (2012). p. 73. [7] A Stevens, H Yang, L Carin et al. Microscopy 63(1), (2014), pp. 41.

Stevens, Andrew↗