Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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 1,009 records · Page 56

Determination of 3D coronal structures from UVCS/SOHO synoptic observations

Daily coronal synoptic observations with the ultraviolet coronagraph spectrometer (UVCS) on SOHO are used to create Carrington maps of line of sight intensities and spectral line widths for H I Lyman alpha and the O 6 1032/1037 doublet. These 2D maps can then be used as inputs to a tomographic inversion routine to produce 3D models of the corona out to the 2.5 solar radius in the polar regions (3.0 solar radius in equatorial regions). The initial results for the determination of O(5+) outflow velocities as a function of heliographic latitude is presented.

Strachan, L.↗

Ocean-Atmosphere Interaction in Climate Changes

The diagram, which attests the El Nino teleconnection observed by the NASA Scatterometer (NSCAT) in 1997, is an example of the results of our research in air-sea interaction - the core component of our three-part contribution to the Climate Variability Program. We have established an interplay among scientific research, which turns spacebased data into knowledge, a push in instrument technology, which improves observations of climate variability, and an information system, which produces and disseminates new data to support our scientific research. Timothy Liu led the proposal for advanced technology, in response to the NASA Post-2002 Request for Information. The sensor was identified as a possible mission for continuous ocean surface wind measurement at higher spatial resolution, and with the unique capability to measure ocean surface salinity. He is participating in the Instrument Incubator Program to improve the antenna technology, and is initiating a study to integrate the concept on Japanese missions. He and his collaborators have set up a system to produce and disseminate high level (gridded) ocean surface wind/stress data from NSCAT and European missions. The data system is being expanded to produce real-time gridded ocean surface winds from Quikscat, and precipitation and evaporation from the Tropical Rain Measuring Mission. It will form the basis for a spacebased data analysis system which will include momentum, heat and water fluxes. The study on 1997 El Nino teleconnection illustrates our interdisciplinary and multisensor approach to study climate variability. The diagram shows that the collapse of trade wind and the westerly wind anomalies in the central equatorial Pacific led to the equatorial ocean warming. The equatorial wind anomalies are connected to the anomalous cyclonic wind pattern in the northeast Pacific. The anomalous warming along the west coast of the United States is the result of the movement of the pre-existing warm sea surface temperature anomalies with the cyclonic wind anomalies toward the coast. The results led to a new study which identifies decadal ocean variations in the Northeast Pacific. Three studies of oceanic responses to wind forcing caused by the seasonal change of monsoons, the passage of a typhoon, and the 1997 El Nino, were successfully conducted. Besides wind forcing, we continue to examine new techniques for estimating thermal and hydrologic fluxes, through the inverse ocean mixed-layer model, through divergence of atmospheric water transport, and by direct retrieval from radiances observed by microwave radiometers. Greenhouse warming has been linked to water vapor measured by two spaceborne sensors in two studies. In the first study, strong baroclinicity and deep convection were found to transport water vapor to the upper atmosphere and increase greenhouse trapping over the storm tracks of the North Pacific and Atlantic. In another study, the annual cycle of greenhouse warming were related to sea surface temperature (SST) and integrated water vapor, and the latitudinal dependence of the magnitudes and phases of the annual cycles were compared.

Liu, W. Timothy↗

Global Distribution and Sources of Volatile and Nonvolatile Aerosol In the Remote Troposphere

Airborne measurements of aerosol (Condensation Nuclei, CN) and selected trace gases made in the areas of the North Atlantic Ocean during SONEX (October/November 1997), and the south tropical Pacific Ocean during PEM-Tropics A (September/October 1996) and PEM-Tropics B (March/April 1999) have been analyzed. Emphasis is on the interpretations of variations in the number densities of the fine (>17 nm) and ultrafine (>8 nm) CN in the upper troposphere (8-12 km). These data suggest that large number densities of highly volatile CN (10(exp 4)-10(exp 5)/cu cm) are present in the clean upper troposphere with highest values over the tropical1subtropical region. Through marine convection and long-distance horizontal transport, volatile CN originating from the tropical/subtropical regions can frequently impact the abundance of aerosol in the middle and upper troposphere at mid to high latitudes. Nonvolatile aerosol particles behave in a manner similar to tracers of combustion (CO) and photochemical pollution (PAN), implying a source from continental pollution of industrial or biomass burning origin. In the upper troposphere, we find that volatile and nonvolatile partials number densities are inversely correlated. An aerosol microphysical model is used to suggest that coagulation of fine volatile particles with fewer larger nonvolatile particles provides one possible mechanism for this relationship. It appears that nonvolatile particles, of principally anthropogenic origin,provide a highly efficient removal process for the fine volatile aerosol.

Singh, Hanwant B.↗

In-Flight Suppression of an Unstable F/A-18 Structural Mode Using the Space Launch System Adaptive Augmenting Control System

NASA's Space Launch System (SLS) Flight Control System (FCS) includes an Adaptive Augmenting Control (AAC) component which employs a multiplicative gain update law to enhance the performance and robustness of the baseline control system for extreme off-nominal scenarios. The SLS FCS algorithm including AAC has been flight tested utilizing a specially outfitted F/A-18 fighter jet in which the pitch axis control of the aircraft was performed by a Non-linear Dynamic Inversion (NDI) controller, SLS reference models, and the SLS flight software prototype. This paper describes test cases from the research flight campaign in which the fundamental F/A-18 airframe structural mode was identified using post-flight frequency-domain reconstruction, amplified to result in closed loop instability, and suppressed in-flight by the SLS adaptive control system.

VanZwieten, Tannen S.↗

In-Flight Suppression of a Destabilized F/A-18 Structural Mode Using the Space Launch System Adaptive Augmenting Control System

NASA's Space Launch System (SLS) Flight Control System (FCS) includes an Adaptive Augmenting Control (AAC) component which employs a multiplicative gain update law to enhance the performance and robustness of the baseline control system for extreme off nominal scenarios. The SLS FCS algorithm including AAC has been flight tested utilizing a specially outfitted F/A-18 fighter jet in which the pitch axis control of the aircraft was performed by a Non-linear Dynamic Inversion (NDI) controller, SLS reference models, and the SLS flight software prototype. This paper describes test cases from the research flight campaign in which the fundamental F/A-18 airframe structural mode was identified using frequency-domain reconstruction of flight data, amplified to result in closed loop instability, and suppressed in-flight by the SLS adaptive control system.

Wall, John H.↗

Aquarius Retrieval of Sea Ice Thickness: Initial Results

Aquarius brightness temperature data are used to calculate sea ice thickness in the Arctic region. The method is based on the inversion of a radiative transfer model for icecovered sea. Using this technique, the initial sea ice thickness values retrieved from Aquarius data are compared to the SMOSIce Data as well as to estimates from NASA's Operation IceBridge. The results show similar trends between the SMOS- and Aquarius-derived sea ice thickness, however the Aquarius estimates tend to be higher and noisier than the corresponding SMOS values. The accuracy of retrieved Aquarius ice thickness is possibly influenced by uncertainties in the ancillary input parameters and by the coarser resolutions of Aquarius.

Sea Ice Thickness↗

EFT-1 Heatshield Aerothermal Environment Reconstruction

The EFT-1 Avcoat heatshield was instrumented with 34 plugs containing multiple thermocouples in-depth with an objective being to measure the flight aerothermal environment. This paper presents a discussion of the instrumentation and the techniques used to reconstruct the heating environment from the measured in-depth temperatures. The inverse heat transfer problem algorithms, models and assumptions will be outlined, and available results will be presented.

Oliver, A. Brandon↗

The Role of Magnetic Helicity in Structuring the Solar Corona

Two of the most widely observed and striking features of the Suns magnetic field are coronal loops, which are smooth and laminar, and prominences or filaments, which are strongly sheared. Loops are puzzling because they show little evidence of tangling or braiding, at least on the quiet Sun, despite the chaotic nature of the solar surface convection. Prominences are mysterious because the origin of their underlying magnetic structure filament channels is poorly understood at best. These two types of features would seem to be quite unrelated and wholly distinct. We argue that, on the contrary, they are inextricably linked and result from a single process: the injection of magnetic helicity into the corona by photospheric motions and the subsequent evolution of this helicity by coronal reconnection. In this paper, we present numerical simulations of the response of a Parker (1972) corona to photospheric driving motions that have varying degrees of helicity preference. We obtain four main conclusions: (1) in agreement with the helicity condensation model of Antiochos (2013), the inverse cascade of helicity by magnetic reconnection in the corona results in the formation of filament channels localized about polarity inversion lines; (2) this same process removes most complex fine structure from the rest of the corona, resulting in smooth and laminar coronal loops; (3) the amount of remnant tangling in coronal loops is inversely dependent on the net helicity injected by the driving motions; and (4) the structure of the solar corona depends only on the helicity preference of the driving motions and not on their detailed time dependence. We discuss the implications of our results for high-resolution observations of the corona.

Knizhnik, K. J.↗

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (−30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (−19%) and ecosystem (−7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Félicien Meunier↗

Designing Learning Experiences With A Low-Cost Robotic Arm

Robots have gained immense popularity in Hollywood and growth globally in both industrial manufacturing and non-manufacturing environments and applications such as healthcare, service sector, and space exploration. To date, there are several examples of simple and low-cost educational robotic platforms and commercially available platforms (e.g., Lego Mindstorms, VEX Robotics) for incorporation into existing curriculum. However, most low-cost examples require access to rapid-prototyping tools, such as 3D printers to manufacture the structure of the robot, while commercially available platforms are relatively expensive (> $1,000). Although the low-cost, open-source examples provide increased access, these examples require design and manufacturing tasks that would be considered outside the learning objectives of an upper-level robotics course and are better suited for other introductory courses. Thus, we asked, how can a robotic platform be incorporated into existing robotics curriculum to enhance students' learning experiences? To explore this research question, we introduce a low-cost (<$200), untethered, and transportable robotic platform that is easy to assemble using off-the-shelf components. This kit can be powered from a laptop computer and does not rely on access to rapid-prototyping tools such as 3D printers or laser cutters, making this a more accessible option in undergraduate engineering courses. Specifically, we aimed to investigate the design of experiential learning experiences for the mathematical modeling of the forward and inverse kinematics of a serial robotic arm that complements existing robotics curriculum. The experiential learning experience focuses on traditional written answer, simulation in MATLAB, and finally implementation on a robotic platform. Few studies examined this accessible option when evaluating experiential learning experiences that complement existing robotics curricula. To assess the impact of this robotic arm kit in an undergraduate course, we implemented an educational intervention that allowed us insight into student perceptions, takeaways on the course and activities involving the robotic arm, and the impact of the course on their career outlook when comparing activities that involved use of the robotic arm and those that did not. Details on the design, development and implementation of the learning activities is provided. Both quantitative and qualitative data were collected and analyzed. The results presented in this paper discuss students finding the learning activities on the robotic arms more helpful than those without, and that students found high value in the hands-on experiences and real-world scenarios offered by the activities using the robotic arm. Challenges to implementation of the robotic arms are discussed, including students’ prior knowledge of using robotic arms.

Eric J. Markvicka↗

Integrating deep neural networks with full-waveform inversion: Reparameterization, regularization, and uncertainty quantification

Full-waveform inversion (FWI) is an accurate imaging approach for modeling the velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong nonlinearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We have adopted a neural-network-based full-waveform inversion (NNFWI) method that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. Furthermore, the velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model determine that NNFWI can mitigate local minima, especially for imaging high-contrast features such as salt bodies, and it significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.

58 GEOSCIENCES↗

Examination of Soil Moisture Retrieval Using SIR-C Radar Data and a Distributed Hydrological Model

A major objective of soil moisture-related hydrological-research during NASA's SIR-C/X-SAR mission was to determine and compare soil moisture patterns within humid watersheds using SAR data, ground-based measurements, and hydrologic modeling. Currently available soil moisture-inversion methods using active microwave data are only accurate when applied to bare and slightly vegetated surfaces. Moreover, as the surface dries down, the number of pixels that can provide estimated soil moisture by these radar inversion methods decreases, leading to less accuracy and, confidence in the retrieved soil moisture fields at the watershed scale. The impact of these errors in microwave- derived soil moisture on hydrological modeling of vegetated watersheds has yet to be addressed. In this study a coupled water and energy balance model operating within a topographic framework is used to predict surface soil moisture for both bare and vegetated areas. In the first model run, the hydrological model is initialized using a standard baseflow approach, while in the second model run, soil moisture values derived from SIR-C radar data are used for initialization. The results, which compare favorably with ground measurements, demonstrate the utility of combining radar-derived surface soil moisture information with basin-scale hydrological modeling.

Hsu, A. Y.↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume numerical techniques with finite thickness models for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the one-dimensional semi -infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody were investigated. An array of streamwise orientated heating striations were generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients due to the striation patterns two-dimensional heat transfer techniques were necessary to obtain accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗

TPSAS-NF1676L-33433-DND

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic scattering of geometric bodies treated with acoustic liners. This work considers using either an impedance boundary condition or admittance boundary condition to allow acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces; admittance is defined to be the inverse of impedance. Three different acoustic liner models will be discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and a Broadband Model. In the Extended Helmholtz Resonator Model and Three-Parameter Impedance Model, the liner impedance is specified at a single frequency; the Broadband Model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for each model will be derived and coupled with a time-domain boundary integral equation. The scattering solution will be obtained iteratively using spatial and temporal basis functions. Stability will be demonstrated through eigenvalue analysis. Moreover, fast algorithms and high performance computing are considered to both assess and reduce the computational cost of the simulation.

Michelle E Rodio↗

Is There an Optimal CO2 Partial Column for Flux Inversions?

The fidelity of flux estimates from an atmospheric inversion depends on the ability of atmospheric transport models to simulate the measured quantity. For species such as CO2, with surface fluxes and a large network of surface measurements, this means correctly simulating the dynamics of the planetary boundary layer (PBL), which is one of the most uncertain aspects of atmospheric transport modeling. In contrast, the simulated total column average mole fraction of CO2 (XCO2) is largely insensitive to simulated PBL dynamics. Therefore, measurements of XCO2 provided by current and future short wave infrared (SWIR) greenhouse gas (GHG) satellites such as GOSAT and the OCO family would seem to be more appropriate to flux inversions, as far as minimizing transport model errors (the "noise") is concerned. Unfortunately, the flux-induced variation of CO2 (the "signal") is the largest within the PBL and smallest in XCO2. Therefore, assimilating XCO2 as opposed to PBL CO2 need not give us the strongest "signal to noise" in flux inversions. Recent work on GOSAT and OCO2 retrievals suggest that SWIR satellite spectra may be used to estimate a lower partial column CO2, which could be assimilated in a flux inversion, instead of XCO2.Here we report on a study to assess whether there is an optimal partial column average CO2, intermediate between PBL CO2 and XCO2, whose assimilation might yield the best signal to noise in flux inversions, where (as before) "signal" is the flux-induced variation and "noise" is the error in transport modeling. We simulate atmospheric CO2 with five different global transport models and a common surface CO2 flux over ten years. We consider the spread across the five models to be a proxy for transport model error (the "noise"), and the common variation of CO2 in all five models to be a proxy for the "signal". We compare these signals and noises at different spatiotemporal scales for different partial column specifications to investigate whether there exists an optimal partial column that has large surface flux-driven variations and yet is relatively insensitive to errors in transport models. Finally, we comment on the feasibility of estimating such a partial column from current and future SWIR GHG satellites in the light of recent work on vertically resolved CO2 from current SWIR GHG satellites.

Basu, Sourish↗

Bayesian operator inference for data-driven reduced-order modeling

This work proposes a Bayesian inference method for the reduced-order modeling of time-dependent systems. Informed by the structure of the governing equations, the task of learning a reduced-order model from data is posed as a Bayesian inverse problem with Gaussian prior and likelihood. The resulting posterior distribution characterizes the operators defining the reduced-order model, hence the predictions subsequently issued by the reduced-order model are endowed with uncertainty. The statistical moments of these predictions are estimated via a Monte Carlo sampling of the posterior distribution. Since the reduced models are fast to solve, this sampling is computationally efficient. Furthermore, the proposed Bayesian framework provides a statistical interpretation of the regularization term that is present in the deterministic operator inference problem, and the empirical Bayes approach of maximum marginal likelihood suggests a selection algorithm for the regularization hyperparameters. The proposed method is demonstrated on two examples: the compressible Euler equations with noise-corrupted observations, and a single-injector combustion process.

97 MATHEMATICS AND COMPUTING↗

Reactor antineutrino flux and anomaly

Reactor antineutrinos have played a significant role in establishing the standard model of particle physics and the theory of neutrino oscillations. In this article, we review the reactor antineutrino flux and in particular the reactor antineutrino anomaly (RAA) coined over a decade ago. RAA refers to a deficit of the measured antineutrino inverse beta decay rates at very short-baseline reactor experiments compared to the theoretically improved predictions (i.e. the Huber–Mueller model). Since the resolution of several previous experimental anomalies have led to the discovery of non-zero neutrino mass and mixing, many efforts have been invested to study the origin of RAA both experimentally and theoretically. The progress includes the observation of discrepancies in antineutrino energy spectrum between data and the Huber–Mueller model, the re-evaluation of the Huber–Mueller model uncertainties, the potential isotope-dependent rate deficits, and the better agreement between data and new model predictions using the improved summation method. Importantly, these developments disfavor the hypothesis of a light sterile neutrino as the explanation of RAA and supports the deficiencies of Huber–Mueller model as the origin. Looking forward, more effort from both the theoretical and experimental sides is needed to fully understand the root of RAA and to make accurate predictions of reactor antineutrino flux and energy spectrum for future discoveries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗