Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DME-Propane Autoignition Measurements inside a Shock Tube

Dimethyl ether (DME) is a biofuel that has the potential to replace diesel in heavy-duty engines. A blend of propane (C3H8) and DME could reduce emissions when compared to diesel at heavy-duty compression engine-relevant conditions. Testing a potential new mixture in a compression engine requires a high-fidelity chemical kinetics model that accurately predicts autoignition delay times, which is especially necessary for a compression engine. In this work, autoignition data has been gathered at an equivalence ratio of 2.0 for pressures of 60 and 80 bar. DME and propane were combusted in synthetic air while excited hydroxyl (OH*) chemiluminescence was used to gather ignition delay times. Data was compared to recent chemical kinetic mechanisms for the two different pressures. Using the mechanism with the best fit of autoignition data, a sensitivity analysis was conducted to analyze the chemical kinetics of the DME/C3H8 combustion at elevated pressures.

Mohammed, Zuhayr Pasha↗
NASA NTRS1971

A study of the flexible body dynamics of the RAE-A satellite

Spacecraft flight data was compared with simulated data computed by the RAE and IMP dynamics simulator computer programs. There was not sufficient data available to provide definitive data on antenna boom position to the Radio Astronomy experimenter. There was, however, enough data to indicate the trend of antenna boom stiffness and temperature gradient. It was possible to obtain a good fit of simulated data to measured spacecraft attitude data when the antenna booms were extended to 450 feet. The value of the RAE-A antenna boom stiffness appears to be significantly less than the preflight ground measured value. It was confirmed that an improved simulation can be obtained by using values of the root angles of the antenna booms and by using an angle of the damper plane of motion, relative to the plane of the antennas, which differed from preflight ground measurements.

Source record↗
NASA NTRS1974

A global magnetic anomaly map

A subset of POGO satellite magnetometer data has been formed that is suitable for analysis of crustal magnetic anomalies. Using a thirteenth order field model, fit to these data, magnetic residuals have been calculated over the world to latitude limits of plus 50 deg. These residuals averaged over one degree latitude-longitude blocks represent a detailed global magnetic anomaly map derived solely from satellite data. Preliminary analysis of the map indicates that the anomalies are real and of geological origin.

Regan, R. D.↗
NASA NTRS2003

Data Access Tools And Services At The Goddard Distributed Active Archive Center (GDAAC)

As one of the largest providers of Earth Science data from the Earth Observing System, GDAAC provides the latest data from the Moderate Resolution Imaging Spectroradiometer (MODIS), Atmospheric Infrared Sounder (AIRS), Solar Radiation and Climate Experiment (SORCE) data products via GDAAC's data pool (50TB of disk cache). In order to make this huge volume of data more accessible to the public and science communities, the GDAAC offers multiple data access tools and services: Open Source Project for Network Data Access Protocol (OPeNDAP), Grid Analysis and Display System (GrADS/DODS) (GDS), Live Access Server (LAS), OpenGlS Web Map Server (WMS) and Near Archive Data Mining (NADM). The objective is to assist users in retrieving electronically a smaller, usable portion of data for further analysis. The OPeNDAP server, formerly known as the Distributed Oceanographic Data System (DODS), allows the user to retrieve data without worrying about the data format. OPeNDAP is capable of server-side subsetting of HDF, HDF-EOS, netCDF, JGOFS, ASCII, DSP, FITS and binary data formats. The GrADS/DODS server is capable of serving the same data formats as OPeNDAP. GDS has an additional feature of server-side analysis. Users can analyze the data on the server there by decreasing the computational load on their client's system. The LAS is a flexible server that allows user to graphically visualize data on the fly, to request different file formats and to compare variables from distributed locations. Users of LAS have options to use other available graphics viewers such as IDL, Matlab or GrADS. WMS is based on the OPeNDAP for serving geospatial information. WMS supports OpenGlS protocol to provide data in GIs-friendly formats for analysis and visualization. NADM is another access to the GDAAC's data pool. NADM gives users the capability to use a browser to upload their C, FORTRAN or IDL algorithms, test the algorithms, and mine data in the data pool. With NADM, the GDAAC provides an environment physically close to the data source. NADM will benefit users with mining or offer data reduction algorithms by reducing large volumes of data before transmission over the network to the user.

Pham, Long↗
NASA NTRS2010

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗
NASA NTRS1992

An improved Venus gravity field from Doppler tracking of the Pioneer Venus Orbiter and Magellan spacecraft

A total of 365,000 Pioneer Venus Orbiter (PVO) S-band Doppler data points were fit in arcs ranging from one to four revolutions in length by estimating the spacecraft position and velocity, atmospheric density, and a solar pressure parameter at each arc epoch. The resultant converged PVO orbits were used to produce 194 information arrays for the Venus spherical harmonic coefficients complete to degree and order 21. A similar procedure was used to produce a second set of 72 information arrays for a 21 x 21 gravitational field from a limited, but global, sampling of Magellan S- and X-band data from cycles 1 and 2. The PVO and Magellan information arrays were combined and solved to obtain the 21 x 21 spherical harmonic expansion designated JPL-Venus Gravity Model 6A (VGM6A). The VGM6A harmonic field field produced significant improvements in the Doppler residual statistics obtained from field validation fits to both PVO and Magellan tracking data in comparison with alternative representations of Venus gravity.

Mcnamee, J. B.↗
DOE OSTI2024

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Robust Interpretation of Electrochemical Impedance Spectra Using Numerical Complex Analysis

Electrochemical Impedance Spectroscopy (EIS) is a noninvasive technique widely used for understanding charge transfer and charge transport processes in electrochemical systems and devices. Standard approaches for the interpretation of EIS data involve starting with a hypothetical circuit model for the physical processes in the device based on experience/intuition and then fitting the EIS data to this circuit model. This work explores a mathematical approach for extracting key characteristic features from EIS data by relying on fundamental principles of complex analysis. These characteristic features can suggest the presence of inductors and constant phase elements (nonideal capacitors) from impedance data and enable us to answer questions about the identifiability and nonuniqueness of equivalent circuit models. In certain scenarios such as models with only resistors and capacitors, we are able to enumerate all possible families of circuit models. Finally, we apply the mathematical framework presented here to real-world electrochemical systems and highlight results using impedance measurements from a lithium-ion battery coin cell.

charge transfer↗
DOE OSTI2024

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Time-dependent mapping of the magnetic field at the core-mantle boundary

We consider the problem of constructing a time-dependent map of the magnetic field at the core-mantle boundary. We use almost all the available data from the last 300 years to produce two maps, one for the period 1690-1840 and the other for 1840-1990. We represent the spatial dependency of the field using spherical harmonics, the time dependency using a cubic B-spline basis, and seek the smoothest solutions compatible with the observations. We argue that, for observations from permanent magnetic observatories, the most efficient strategy is to use the first differences of annual means; for satellite data, the most efficient strategy is simply to limit the number of data used so as to minimize any tendency to map the crustal field into the core field. The resulting model fits the observatory data better than any previous model. The resulting time-dependent field map exhibits much of the same structure in the field and its secular variation identified in earlier studies.

Bloxham, Jeremy↗
NASA NTRS2005

Influence of Nucleation Mechanisms on the Radiative Properties of Deep Convective Clouds and Subvisible Cirrus in CRYSTAL/FACE

During the past few years we have conducted work on several different topics, as reflected by our publications. As one of the Co-Project scientists for The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL FACE) we worked to help design the mission and then conduct it in the field. Another major activity during the past two years has been to pull together various groups to formulate plans for follow on missions to CRYSTAL FACE. We organized a workshop at the University of Colorado during the summer of 2003 to assess the best locations for future missions. Working with a group of about 10 scientists from around the country we prepared a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. In conjunction with the group at NASA Ames we have helped build a new numerical model for deep convection and have applied that model to simulate the CRYSTAL data. Our goal in particular has been to better understand how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results have just been submitted to Geophys. Res. Lett (Smith et al., 2005). Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity and isotope abundances with this model. We have also been heavily involved in trying to improve our understanding of nitric acid condensation on ice. Gao et al (2004) have shown that water supersaturations above ice occur when the atmosphere is supersaturated with respect to nitric acid trihydrate. As one of the co-authors of that work, we suggested the mechanism that may explain why this is occurring. Essentially, ice does not like to grow near unit supersaturation, but does so because the water molecules can find sites on the ice surface to attach themselves to before they fly off the ice surface. This phenomena was well known in the 1960s when it was a source of debate about whether condensation and evaporation coefficients for ice would be the same. Evaporation does not require any molecular orientation, while condensation does, so it was possible that the coefficients would differ. They don't differ because the water molecules rapidly move across the surface and find places to attach. Nitric acid may be occupying these preferred sites and therefore the water molecules can't find a desirable place to attach. We anticipate that this research will be the subject of laboratory work during the coming few years. Another possibility that has been suggested is that cubic ice is forming in clouds. We have measured the vapor pressure of cubic ice, and plan to publish that result in the next few months. We have also been working on additional aspects of the condensation of nitric acid on ice. With Y. Kondo we studied the condensation of NOy on ice using the SOLVE data. Gamblin et al. have continued this work. The CRYSTAL NOy and HNO, groups have shown th their data can be fit using standard Langmuir isotherms as suggested in some, but not all, laboratory studies. We have found in the SOLVE data set that this is not the case. Moreover some laboratory studies show there are important kinetic effects that may be occurring in the atmosphere limiting the transfer of nitric acid to the ice. The SOLVE data seem consistent with these studies. We are currently re-analyzing the CRYSTAL data to look for these kinetic effects. There are a number of implications of these studies. One of the more interesting is that the nitric acid coating on ice can be used as a cloud clock to determine how long the cloud parcel has been in existence. We have also been involved with several laboratory studies. We have worked to improve the database on ice optical constants, which are critical for remote sensing. We have also studied the ways in which ice nucleates on clays. We suspect now that the standard theories used for depositional ice nucleation are completely incorrect. Further work will be needed to develop a new theory.

Toon, Owen B.↗
NASA NTRS2019

Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

space-borne imagery↗
DOE OSTI2024

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1987

Microwave dielectric behavior of vegetation material

The microwave dielectric behavior of vegetation was examined through the development of theoretical models involving dielectric dispersion by both bound and free water and supported by extensive dielectric measurements conducted over a wide range of conditions. The experimental data were acquired using an open-ended coaxial probe that was developed for sensing the dielectric constant of thin layers of materials, such as leaves, from measurements of the complex reflection coefficient using a network analyzer. The probe system was successfully used to record the spectral variation of the dielectric constant over a wide frequency range extending from 0.5 to 20.4 GHz at numerous temperatures between -40 to +40 C. The vegetation samples were measured over a wide range of moisture conditions. To model the dielectric spectrum of the bound water component of the water included in vegetation, dielectric measurements were made for several sucrose-water solutions as analogs for the situation in vegetation. The results were used in conjunction with the experimental data for leaves to determine some of the constant coefficients in the theoretical models. Two models, both of which provide good fit to the data, are proposed.

Elrayes, Mohamed A.↗
NASA NTRS2006

Rotating Rake Mode Measurements Over Passive Treatment in a Ducted Fan

The NASA Glenn Research Center s Rotating Rake mode measurement system has been successful in measuring the modal content propagating in hardwall ducts. This paper proposes an extension of the Rotating Rake measurement and analysis technique to treated sections by developing basis functions based on wall impedance boundary conditions for flow conditions (i.e., constant duct area and Mach number) where the closed form analytical solution exists. Analytical equations developed to estimate mode power are incorporated. This method is verified by decomposing and analyzing radial pressure profiles generated numerically by the Eversman propagation code. Several modes, frequencies and impedances are evaluated. Data from a low-speed ducted fan with several different impedance conditions was acquired and reduced to determine the best fit to the data. Using the impedance boundary conditions result in better mode measurement solutions.

Sutliff, Daniel L.↗
NASA NTRS2007

On Problems Associated with Modeling Wing-Tail Configurations from Wind Tunnel Data

This paper considers factors that contribute to poor identification of unsteady aerodynamics from wind tunnel data for an airliner configuration. One approach to modeling a wing-tail configuration is considered and applied to both steady and large-amplitude forced pitch oscillation wind tunnel data taken over a wide range of angles of attack but a limited range of amplitude and frequencies. The identified models fit the measured data well but in some cases with inaccurate parameters. Only limited conclusions can be drawn from analysis of the current data set until further experiments can be performed to resolve the identification issues. The analysis of measured and simulated data provides some insights and guidance on how an effective experiment may be designed for wing-tail configurations with nonlinear unsteady aerodynamics.

Murphy, Patrick C.↗