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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

On Characterizing Particle Shape

It is well known that particle shape affects flow characteristics of granular materials, as well as a variety of other solids processing issues such as compaction, rheology, filtration and other two-phase flow problems. The impact of shape crosses many diverse and commercially important applications, including pharmaceuticals, civil engineering, metallurgy, health, and food processing. Two applications studied here include the dry solids flow of lunar simulants (e.g. JSC-1, NU-LHT-2M, OB-1), and the flow properties of wet concrete, including final compressive strength. A multi-dimensional generalized, engineering method to quantitatively characterize particle shapes has been developed, applicable to both single particle orientation and multi-particle assemblies. The two-dimension, three dimension inversion problem is also treated, and the application of these methods to DEM model particles will be discussed. In the case of lunar simulants, flow properties of six lunar simulants have been measured, and the impact of particle shape on flowability - as characterized by the shape method developed here -- is discussed, especially in the context of three simulants of similar size range. In the context of concrete processing, concrete construction is a major contributor to greenhouse gas production, of which the major contributor is cement binding loading. Any optimization in concrete rheology and packing that can reduce cement loading and improve strength loading can also reduce currently required construction safety factors. The characterization approach here is also demonstrated for the impact of rock aggregate shape on concrete slump rheology and dry compressive strength.

Ennis, Bryan J.↗
NASA NTRS2018

An Active Subspace Method for Accelerating Convergence in Delaunay-Based Optimization via Dimension Reduction

Delaunay-based derivative-free optimization, ∆DOGS, is an efficient and provably-convergent global optimization method for the problems which has computationally expensive objection function and the analytical expression for the objective function is not available. ∆-DOGS is a novel optimization scheme in the family of response surface methods (RSMs); however, it suffers from the curse of dimensionality since the computational cost increases dramatically as the number of design parameters increases. As a result, the number of design parameters in ∆-DOGS algorithm is relatively low (n.10). To avoid such problems, this paper proposes a combination of derivative-free optimization, seeking the global minimizer of an expensive and nonconvex objective function f(x) and active subspace method, detecting the directions of the most variability using evaluations of the gradient. The contribution of other directions to the objective function is bounded by a sufficiently small constant. This new algorithm iteratively applied Delaunay-based derivative-free optimization to seek the minimizer on the d-dimensional active subspace that has most function variation. Inverse mapping is needed to project data from active subspace to full-model for evaluating function values. This task is overcome by solving an inequality constrained problem that curves the response surface of the objective function. The test results show that this strategy is effective on a handful of optimization problems.

Bewley, Thomas R.↗
NASA NTRSDate not supplied

LOFTID Surface Heating Reconstruction

On November 10, 2022, the Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) reentry vehicle launched to low-Earth orbit aboard a United Launch Alliance Atlas V rocket out of Vandenberg Air Force Base. The aeroshell, the largest Hypersonic Inflatable Aerodynamic Decelerator (HIAD) ever flown, was inflated to its full 6-meter diameter before successful re-entry into the atmosphere. The aeroshell was heavily instrumented in order to understand its behavior during entry. There were 82 thermocouples (TCs) distributed across the aeroshell, with 22 integrated into the flexible thermal protection system (FTPS) on the rigid nose, 36 in the FTPS on the deployable structure, and 24 on the inflatable structure. TCs were placed at different depths throughout the FTPS. Those nearest to the surface were located just beneath the two SiC outer fabric layers. The near-surface TCs on the rigid nose were Type R with flame spray alumina insulation, while those on the flank were Type N with mica/ceramic insulation. Additionally, a radiometer was placed at the center of the nose surrounded by four total heat flux gauges in a cruciform configuration at a radius of 0.41 m. The nose instrumentation is shown in Fig. 1 and a cross-section of the aeroshell with all TC locations is shown in Fig. 2. The objective of this work was to use the temperatures measured by the TCs during flight to estimate the surface heat rate across the aeroshell throughout the period of re-entry by inverse analysis methodology. The results were used to evaluate the fidelity of measurements from the total heat flux gauges on the nose, determine the surface heat flux at aeroshell locations where gauges were not present, and compare to pre-flight CFD-based heating predictions. Inversely estimated surface heat flux continues to be used to correlate FTPS thermal models to reconstruct in-flight thermal response.

H S Alpert↗
NASA NTRSDate not supplied

LOFTID Surface Heating Reconstruction

On November 10, 2022, the Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) reentry vehicle launched to low-Earth orbit aboard a United Launch Alliance Atlas V rocket out of Vandenberg Air Force Base. The aeroshell, the largest Hypersonic Inflatable Aerodynamic Decelerator (HIAD) ever flown, was inflated to its full 6-meter diameter before successful re-entry into the atmosphere. The aeroshell was heavily instrumented in order to understand its behavior during entry. There were 82 thermocouples (TCs) distributed across the aeroshell, with 22 integrated into the flexible thermal protection system (FTPS) on the rigid nose, 36 in the FTPS on the deployable structure, and 24 on the inflatable structure. TCs were placed at different depths throughout the FTPS. Those nearest to the surface were located just beneath the two SiC outer fabric layers. The near-surface TCs on the rigid nose were Type R with flame spray alumina insulation, while those on the flank were Type N with mica/ceramic insulation. Additionally, a radiometer was placed at the center of the nose surrounded by four total heat flux gauges in a cruciform configuration at a radius of 0.41 m. The nose instrumentation is shown in Fig. 1 and a cross-section of the aeroshell with all TC locations is shown in Fig. 2. The objective of this work was to use the temperatures measured by the TCs during flight to estimate the surface heat rate across the aeroshell throughout the period of re-entry by inverse analysis methodology. The results were used to evaluate the fidelity of measurements from the total heat flux gauges on the nose, determine the surface heat flux at aeroshell locations where gauges were not present, and compare to pre-flight CFD-based heating predictions. Inversely estimated surface heat flux continues to be used to correlate FTPS thermal models to reconstruct in-flight thermal response.

LOFTID↗
NASA NTRS2016

Satellite-Based Spatiotemporal Trends in PM2.5 Concentrations: China 2004-2013

Three decades of rapid economic development is causing severe and widespread PM2.5(particulate matter (is) less than 2.5 ) pollution in China. However, research on the health impacts of PM2.5 exposure has been hindered by limited historical PM2.5 concentration data. We estimated ambient PM2.5 concentrations from 2004 to 2013 in China at 0.1 deg resolution using the most recent satellite data and evaluated model performance with available ground observations. We developed a two-stage spatial statistical model using the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 aerosol optical depth (AOD) and assimilated meteorology, land use data, and PM2.5 concentrations from China's recently established ground monitoring network. An inverse variance weighting (IVW) approach was developed to combine MODIS Dark Target and Deep Blue AOD to optimize data coverage. We evaluated model predicted PM2.5 concentrations from 2004 to early 2014 using ground observations. The overall model cross-validation R(sup 2) and relative prediction error were 0.79 and 35.6%, respectively. Validation beyond the model year (2013) indicated that it accurately predicted PM(sub 2.5) concentrations with little bias at the monthly (R(sup 2) = 0.73), regression slope = 0.91) and seasonal (R(sup 2) = 0.79), regression slope = 0.92) levels. Seasonal variations revealed that winter was the most polluted season and that summer was the cleanest season. Analysis of predicted PM2.5 levels showed a mean annual increase of 1.97 micro-g/cu cm between 2004 and 2007 and a decrease of 0.46 micro-g/cu cm between 2008 and 2013. Our satellite-driven model can provide reliable historical PM2.5 estimates in China at a resolution comparable to those used in epidemiologic studies on the health effects of long-term PM2.5 exposure in North America. This data source can potentially advance research on PM2.5 health effects in China.

satellite remote sensing↗
NASA NTRS2004

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗
NASA NTRS1976

Global atomic hydrogen density derived from OGO-6 Lyman-alpha measurements

The paper analyzes a one-year set of Lyman-alpha airglow data measured in the local zenith at altitudes from 400 to 1100 km by a UV photometer aboard OGO-6. The zenith-intensity data are fitted to theoretical airglow calculations in four spherically symmetric models of the hydrogen geocorona to determine both the Ly-alpha solar flux at line center and the average atomic hydrogen column density. After correcting for a loss of instrument sensitivity, the Ly-alpha flux is found to be linearly correlated with daily Zurich sunspot number. It is also found that the hydrogen density is inversely correlated with Jacchia exospheric temperature, but the dependence is not that predicted by steady-state models with Jeans evaporative escape as the only loss mechanism. It is suggested that charge-exchange production of fast hydrogen atoms from 'hot' ionospheric protons might provide the additional loss this result requires.

Thomas, G. E.↗
NASA NTRS1994

Multi-application controls: Robust nonlinear multivariable aerospace controls applications

This viewgraph presentation describes the general methodology used to apply Honywell's Multi-Application Control (MACH) and the specific application to the F-18 High Angle-of-Attack Research Vehicle (HARV) including piloted simulation handling qualities evaluation. The general steps include insertion of modeling data for geometry and mass properties, aerodynamics, propulsion data and assumptions, requirements and specifications, e.g. definition of control variables, handling qualities, stability margins and statements for bandwidth, control power, priorities, position and rate limits. The specific steps include choice of independent variables for least squares fits to aerodynamic and propulsion data, modifications to the management of the controls with regard to integrator windup and actuation limiting and priorities, e.g. pitch priority over roll, and command limiting to prevent departures and/or undesirable inertial coupling or inability to recover to a stable trim condition. The HARV control problem is characterized by significant nonlinearities and multivariable interactions in the low speed, high angle-of-attack, high angular rate flight regime. Systematic approaches to the control of vehicle motions modeled with coupled nonlinear equations of motion have been developed. This paper will discuss the dynamic inversion approach which explicity accounts for nonlinearities in the control design. Multiple control effectors (including aerodynamic control surfaces and thrust vectoring control) and sensors are used to control the motions of the vehicles in several degrees-of-freedom. Several maneuvers will be used to illustrate performance of MACH in the high angle-of-attack flight regime. Analytical methods for assessing the robust performance of the multivariable control system in the presence of math modeling uncertainty, disturbances, and commands have reached a high level of maturity. The structured singular value (mu) frequency response methodology is presented as a method for analyzing robust performance and the mu-synthesis method will be presented as a method for synthesizing a robust control system. The paper concludes with the author's expectations regarding future applications of robust nonlinear multivariable controls.

Enns, Dale F.↗
NASA NTRS1983

Electrodynamic coupling in magnetically confined X-ray plasmas of astrophysical origin

Explicitly magnetic field-independent relationships that quantitatively connect the thermodynamic state of the confined beta-less-than-one X-ray plasma to the confining magnetic loop characteristics and the stressing velocity at beta greater than about one are presented. The derived relations fall into two categories: (1) diffuse heating resulting from efficient transport of energy throughout the magnetic loop volume, and (2) isolated heating within the electrodynamic dissipation shell. For each of these cases the Joule dissipation, shear viscous dissipation, and compressional dissipation are considered. Each of these models is characterized by three thermodynamic scalings that are parametrized by the magnetic loop properties and by the stressing velocity. Eliminating the magnetic field intermediary from these relations and normalizing the inverse aspect ratio and the magnetic expansion factor with respect to solar conditions results in a single, unique velocity scaling law for each model.

Ionson, J. A.↗
NASA NTRS2022

Idaho Wildfires II: Assessing the Relationship Between Drought Indicators and Wildfire Risk to Enhance Hazard Modeling and Inform Mitigation Planning

The western United States has experienced twenty years of increased and prolonged drought which have exacerbated wildfire hazards. These jeopardize population centers through increased risks to ecosystem services, local economies, and livelihoods. The Idaho Office of Emergency Management, Water Resources, and Department of Lands are seeking methods to dynamically monitor these conditions and update models that inform hazard mitigation planning and resource allocation. Towards this, these agencies partnered with NASA DEVELOP to produce drought-enhanced wildfire hazard models. Part of a two-term project, the two teams revised the state’s static wildfire hazard model with refined data layers and remotely-sensed data to reflect dynamic ecosystem responses to drought conditions and wildfire potential. Our team distinguished between rangeland and forestland ecosystems, and investigated relationships between drought metrics and vegetation condition using TerrSet Earth Trends Modeler. This analysis determined that total precipitation at a 5-month lag interval (r 2 = 0.72) along with the Evaporative Stress Index (r 2 = 0.69); and precipitation at a 5-month interval (r 2 = 0.42) were important drivers in rangeland and forestland, respectively. These driver variables were incorporated into a temporally dynamic wildfire hazard map. Our team used linear regression to correlate hazard ratings with wildfire frequency. For the year 2020, neither the enhanced hazard model (p < 0.10, r 2 = 0.01) nor the state’s static model (p < 0.05, r 2 = 0.03) were strongly correlated with actual wildfire frequency as they expressed an inverse relationship between wildfire hazard and frequency. This suggests wildfire occurrence is complex and not necessarily driven by the variables used.

Wildfire↗
NASA NTRS1989

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z=(z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y (sup 0) = (y (sub 1) (sup 0), ..., y (sub D (sup 0)), using full or partial knowledge of the statistical distribution of the random errors in y (sup 0). The data space Y containing y(sup 0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x, Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. Confidence set inference is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface.

Backus, George E.↗
NASA NTRSDate not supplied

Comparing Acoustic Prediction Methods for Additively Manufactured Porous Strutures

While macroscale methods for predicting the acoustic properties of porous structures have been popular in the past, they often require time-consuming manufacturing and testing workflows. Meanwhile, microscale approaches allow the prediction of transport parameters based exclusively on a periodic structure’s unit cell geometry. Here, we compare these methods to predict the characteristic impedance of additively manufactured porous structures. We use the microscale approach to estimate the geometry’s transport parameters, then predict the characteristic properties using the Johnson-Champoux-Allard (JCA) model. We measure the acoustic properties of the printed structures using a normal incidence impedance tube and estimate the transport parameters using an inverse characterization approach. We use the two-thickness method as a macroscale approach to predict the characteristic properties from the measured surface impedances of two sample thicknesses. Finally, we compare these characteristic prediction methods. Our results show that the inverse characterization and two-thickness methods offer the closest match to the measured values at low frequencies.

impedance↗
NASA NTRS1977

Three-dimensional vector modeling and restoration of flat finite wave tank radiometric measurements

In this paper, a three-dimensional Fourier transform inversion method describing the interaction between water surface emitted radiation from a flat finite wave tank and antenna radiation characteristics is reported. The transform technique represents the scanning of the antenna mathematically as a correlation. Computation time is reduced by using the efficient and economical fast Fourier transform algorithm. To verify the inversion method, computations have been made and compared with known data and other available results. The technique has been used to restore data of the finite wave tank system and other available antenna temperature measurements made at the Cape Cod Canal. The restored brightness temperatures serve as better representations of the emitted radiation than the measured antenna temperatures.

Truman, W. M.↗
NASA NTRS2010

A Parameter Identification Method for Helicopter Noise Source Identification and Physics-Based Semi-Empirical Modeling

A new physics-based parameter identification method for rotor harmonic noise sources is developed using an acoustic inverse simulation technique. This new method allows for the identification of individual rotor harmonic noise sources and allows them to be characterized in terms of their individual non-dimensional governing parameters. This new method is applied to both wind tunnel measurements and ground noise measurements of two-bladed rotors. The method is shown to match the parametric trends of main rotor Blade-Vortex Interaction (BVI) noise, allowing accurate estimates of BVI noise to be made for operating conditions based on a small number of measurements taken at different operating conditions.

Greenwood, Eric, II↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗
NASA NTRS2020

Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters: Experimental Validation and Model Investigations

The recently developed collisional radiative (CR) models for neutral xenon (Xe I) and singly charged xenon ion (Xe II) described by Chaplin et al. in the paper "Insights from Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters" also submitted to the 2021 AIAA Propulsion and Energy Forum are more thoroughly investigated. Comparisons between the existing and recently developed JPL Xe I CR models are done and reasons for discrepancies between the two models are explored. An experiment to validate the Xe II CR model is carried out in the anode region of the UCLA Plasma interactions facility. The Xe II CR model is applied to measured xenon ion emission line intensities to determine electron temperature and density. Spatial resolution of these optical measurements is obtained through an Abel inversion for direct comparison to Langmuir probe results. Although the measured electron densities were on the order of $10^{17}$ m$^{-3}$ and the Xe II CR model is most suitable for use with densities above $10^{18}$ m$^{-3}$ due to the strong dependence of the line ratios on electron temperature in that regime, initial results from the validation experiment show reasonable agreements between the optical and Langmuir probe results. The 441 nm to 605 nm line ratio appears to have the most promise for non-intrusive electron temperature measurements with the Xe II CR model. Shortcomings of the Xe II CR model revealed in the validation experiment are to be explored further.

Lobbia, Robert B.↗